Brad Sherborne

Selective small-molecule inhibition of an RNA structural element.

Riboswitches are non-coding RNA structures located in messenger RNAs that bind endogenous ligands, such as a specific metabolite or ion, to regulate gene expression. As such, riboswitches serve as a novel, yet largely unexploited, class of emerging drug targets. Demonstrating this potential, however, has proven difficult and is restricted to structurally similar antimetabolites and semi-synthetic analogues of their cognate ligand, thus greatly restricting the chemical space and selectivity sought for such inhibitors. Here we report the discovery and characterization of ribocil, a highly selective chemical modulator of bacterial riboflavin riboswitches, which was identified in a phenotypic screen and acts as a structurally distinct synthetic mimic of the natural ligand, flavin mononucleotide, to repress riboswitch-mediated ribB gene expression and inhibit bacterial cell growth. Our findings indicate that non-coding RNA structural elements may be more broadly targeted by synthetic small molecules than previously expected.

De Novo Design of Protein Kinase Inhibitors by in Silico Identification of Hinge Region-Binding Fragments

Protein kinases constitute an attractive family of enzyme targets with high relevance to cell and disease biology. Small molecule inhibitors are powerful tools to dissect and elucidate the function of kinases in chemical biology research and to serve as potential starting points for drug discovery. However, the discovery and development of novel inhibitors remains challenging. Here, we describe a structure-based de novo design approach that generates novel, hinge-binding fragments that are synthetically feasible and can be elaborated to small molecule libraries. Starting from commercially available compounds, core fragments were extracted, filtered for pharmacophoric properties compatible with hinge-region binding, and docked into a panel of protein kinases. Fragments with a high consensus score were subsequently short-listed for synthesis. Application of this strategy led to a number of core fragments with no previously reported activity against kinases. Small libraries around the core fragments were synthesized, and representative compounds were tested against a large panel of protein kinases and subjected to co-crystallization experiments. Each of the tested compounds was active against at least one kinase, but not all kinases in the panel were inhibited. A number of compounds showed high ligand efficiencies for therapeutically relevant kinases; among them were MAPKAP-K3, SRPK1, SGK1, TAK1, and GCK for which only few inhibitors are reported in the literature.

Fragment-based discovery of 6-substituted isoquinolin-1-amine based ROCK-I inhibitors.

Fragment-based NMR screening of a small literature focused library led to identification of a historical thrombin/FactorXa building block, 17A, that was found to be a ROCK-I inhibitor. In the absence of an X-ray structure, fragment growth afforded 6-substituted isoquinolin-1-amine derivatives which were profiled in the primary ROCK-I IMAP assay. Compounds 23A and 23E were selected as fragment optimized hits for further profiling. Compound 23A has similar ROCK-1 affinity, potency and cell based efficacy to the first generation ROCK inhibitors, however, it has a superior PK profile in C57 mouse. Compound 23E demonstrates the feasibility of improving ROCK-1 affinity, potency and cell based efficacy for the series, however, it has a poor PK profile relative to 23A.

Optimisation of 6-substituted isoquinolin-1-amine based ROCK-I inhibitors

Rho kinase is an important target implicated in a variety of cardiovascular diseases. Herein, we report the optimisation of the fragment derived ATP-competitive ROCK inhibitors 1 and 2 into lead compound 14A. The initial goal of improving ROCK-I potency relative to 1, whilst maintaining a good PK profile, was achieved through removal of the aminoisoquinoline basic centre. Lead 14A was equipotent against both ROCK-I and ROCK-II, showed good in vivo efficacy in the spontaneous hypertensive rat model, and was further optimised to demonstrate the scope for improving selectivity over PKA versus hydroxy Fasudil 3.

The importance of protonation and tautomerization in relative binding affinity prediction: a comparison of AMBER TI and Schrödinger FEP.

AbstractIn drug discovery, protonation states and tautomerization are easily overlooked. Through a Merck–Rutgers collaboration, this paper re-examined the initial settings and preparations for the Thermodynamic Integration (TI) calculation in AMBER Free-Energy Workflows, demonstrating the value of careful consideration of ligand protonation and tautomer state. Finally, promising results comparing AMBER TI and Schrödinger FEP+ are shown that should encourage others to explore the value of TI in routine Structure-based Drug Design.

Discovery and optimisation of a selective non-steroidal glucocorticoid receptor antagonist.

High-throughput screening of 3.87 million compounds delivered a novel series of non-steroidal GR antagonists. Subsequent rounds of optimisation allowed progression from a non-selective ligand with a poor ADMET profile to an orally bioavailable, selective, stable, glucocorticoid receptor antagonist.

Workflows and performances in the ranking prediction of 2016 D3R Grand Challenge 2: lessons learned from a collaborative effort

The 2016 D3R Grand Challenge 2 includes both pose and affinity or ranking predictions. This article is focused exclusively on affinity predictions submitted to the D3R challenge from a collaborative effort of the modeling and informatics group. Our submissions include ranking of 102 ligands covering 4 different chemotypes against the FXR ligand binding domain structure, and the relative binding affinity predictions of the two designated free energy subsets of 15 and 18 compounds. Using all the complex structures prepared in the same way allowed us to cover many types of workflows and compare their performances effectively. We evaluated typical workflows used in our daily structure-based design modeling support, which include docking scores, force field-based scores, QM/MM, MMGBSA, MD-MMGBSA, and MacroModel interaction energy estimations. The best performing methods for the two free energy subsets are discussed. Our results suggest that affinity ranking still remains very challenging; that the knowledge of more structural information does not necessarily yield more accurate predictions; and that visual inspection and human intervention are considerably important for ranking. Knowledge of the mode of action and protein flexibility along with visualization tools that depict polar and hydrophobic maps are very useful for visual inspection. QM/MM-based workflows were found to be powerful in affinity ranking and are encouraged to be applied more often. The standardized input and output enable systematic analysis and support methodology development and improvement for high level blinded predictions.

Macrocyclic α helical peptide therapeutic modality: A perspective of learnings and challenges

Macrocyclic α-helical peptides have emerged as a compelling new therapeutic modality to tackle targets confined to the intracellular compartment. Within the scope of hydrocarbon-stapling there has been significant progress to date, including the first stapled α-helical peptide to enter into clinical trials. The principal design concept of stapled α-helical peptides is to mimic a cognate (protein) ligand relative to binding its target via an α-helical interface. However, it was the proclivity of such stapled α-helical peptides to exhibit cell permeability and proteolytic stability that underscored their promise as unique macrocyclic peptide drugs for intracellular targets. This perspective highlights key learnings as well as challenges in basic research with respect to structure-based design, innovative chemistry, cell permeability and proteolytic stability that are essential to fulfill the promise of stapled α-helical peptide drug development.

Use of molecular modeling aided design to dial out hERG liability in adenosine A2A receptor antagonists

Molecular modeling was performed on a triazolo quinazoline lead compound to help develop a series of adenosine A2A receptor antagonists with improved hERG profile. Superposition of the lead compound onto MK-499, a benchmark hERG inhibitor, combined with pKa calculations and measurement, identified terminal fluorobenzene to be responsible for hERG activity. Docking of the lead compound into an A2A crystal structure suggested that this group is located at a flexible, spacious, and solvent-exposed opening of the binding pocket, making it possible to tolerate various functional groups. Transformation analysis (MMP, matched molecular pair) of in-house available experimental data on hERG provided suggestions for modifications in order to mitigate this liability. This led to the synthesis of a series of compounds with significantly reduced hERG activity. The strategy used in the modeling work can be applied to other medicinal chemistry programs to help improve hERG profile.

Inhibitor of intramembrane protease RseP blocks the σE response causing lethal accumulation of unfolded outer membrane proteins

Significance The σE stress response monitors outer membrane protein (OMP) assembly. Uninduced, σE is sequestered to the plasma membrane by its anti-sigma factor, RseA. Mutations perturbing OMP biogenesis induce degradation of RseA by the proteases DegS and RseP, liberating σE so it can activate gene expression. σE activity is essential in wild-type Escherichia coli, although why it is essential has remained unclear. We report that batimastat is an inhibitor of RseP, preventing it from cleaving RseA and thereby causing a lethal decrease in σE activity. Surprisingly, lethality is caused by the accumulation of unfolded OMPs, despite a wild-type OMP biogenesis pathway. Hence, σE is essential because it must fine-tune OMP synthesis, assembly, and degradation to prevent the appearance of toxic unfolded OMPs. The outer membrane (OM) of Gram-negative bacteria forms a robust permeability barrier that blocks entry of toxins and antibiotics. Most OM proteins (OMPs) assume a β-barrel fold, and some form aqueous channels for nutrient uptake and efflux of intracellular toxins. The Bam machine catalyzes rapid folding and assembly of OMPs. Fidelity of OMP biogenesis is monitored by the σE stress response. When OMP folding defects arise, the proteases DegS and RseP act sequentially to liberate σE into the cytosol, enabling it to activate transcription of the stress regulon. Here, we identify batimastat as a selective inhibitor of RseP that causes a lethal decrease in σE activity in Escherichia coli, and we further identify RseP mutants that are insensitive to inhibition and confer resistance. Remarkably, batimastat treatment allows the capture of elusive intermediates in the OMP biogenesis pathway and offers opportunities to better understand the underlying basis for σE essentiality.