Human rhinovirus promotes STING trafficking to replication organelles to promote viral replication.

Human rhinovirus (HRV), like coronavirus (HCoV), are positive-strand RNA viruses that cause both upper and lower respiratory tract illness, with their replication facilitated by concentrating RNA-synthesizing machinery in intracellular compartments made of modified host membranes, referred to as replication organelles (ROs). Here we report a non-canonical, essential function for stimulator of interferon genes (STING) during HRV infections. While the canonical function of STING is to detect cytosolic DNA and activate inflammatory responses, HRV infection triggers the release of STIM1-bound STING in the ER by lowering Ca2+, thereby allowing STING to interact with phosphatidylinositol 4-phosphate (PI4P) and traffic to ROs to facilitates viral replication and transmission via autophagy. Our results thus hint a critical function of STING in HRV viral replication and transmission, with possible implications for other RO-mediated RNA viruses.

Automated high throughput pKa and distribution coefficient measurements of pharmaceutical compounds for the SAMPL8 blind prediction challenge.

The goal of the Statistical Assessment of the Modeling of Proteins and Ligands (SAMPL) challenge is to improve the accuracy of current computational models to estimate free energy of binding, deprotonation, distribution and other associated physical properties that are useful for the design of new pharmaceutical products. New experimental datasets of physicochemical properties provide opportunities for prospective evaluation of computational prediction methods. Here, aqueous pKa and a range of bi-phasic logD values for a variety of pharmaceutical compounds were determined through a streamlined automated process to be utilized in the SAMPL8 physical property challenge. The goal of this paper is to provide an in-depth review of the experimental methods utilized to create a comprehensive data set for the blind prediction challenge. The significance of this work involves the use of high throughput experimentation equipment and instrumentation to produce acid dissociation constants for twenty-three drug molecules, as well as distribution coefficients for eleven of those molecules.

Characterization of P. falciparum dipeptidyl aminopeptidase 3 specificity identifies differences in amino acid preferences between peptide-based substrates and covalent inhibitors.

Malarial dipeptidyl aminopeptidases (DPAPs) are cysteine proteases important for parasite development thus making them attractive drug targets. In order to develop inhibitors specific to the parasite enzymes, it is necessary to map the determinants of substrate specificity of the parasite enzymes and its mammalian homologue cathepsin C (CatC). Here, we screened peptide-based libraries of substrates and covalent inhibitors to characterize the differences in specificity between parasite DPAPs and CatC, and used this information to develop highly selective DPAP1 and DPAP3 inhibitors. Interestingly, while the primary amino acid specificity of a protease is often used to develop potent inhibitors, we show that equally potent and highly specific inhibitors can be developed based on the sequences of nonoptimal peptide substrates. Finally, our homology modelling and docking studies provide potential structural explanations of the differences in specificity between DPAP1, DPAP3, and CatC, and between substrates and inhibitors in the case of DPAP3. Overall, this study illustrates that focusing the development of protease inhibitors solely on substrate specificity might overlook important structural features that can be exploited to develop highly potent and selective compounds.

Author Correction: Design of amidobenzimidazole STING receptor agonists with systemic activity.

Change history: In this Letter, author Ana Puhl was inadvertently omitted; this error has been corrected online.An amendment to this paper has been published and can be accessed via a link at the top of the paper.

D3R Grand Challenge 3: blind prediction of protein-ligand poses and affinity rankings.

The Drug Design Data Resource aims to test and advance the state of the art in protein-ligand modeling by holding community-wide blinded, prediction challenges. Here, we report on our third major round, Grand Challenge 3 (GC3). Held 2017-2018, GC3 centered on the protein Cathepsin S and the kinases VEGFR2, JAK2, p38-α, TIE2, and ABL1, and included both pose-prediction and affinity-ranking components. GC3 was structured much like the prior challenges GC2015 and GC2. First, Stage 1 tested pose prediction and affinity ranking methods; then all available crystal structures were released, and Stage 2 tested only affinity rankings, now in the context of the available structures. Unique to GC3 was the addition of a Stage 1b self-docking subchallenge, in which the protein coordinates from all of the cocrystal structures used in the cross-docking challenge were released, and participants were asked to predict the pose of CatS ligands using these newly released structures. We provide an overview of the outcomes and discuss insights into trends and best-practices.

Design of amidobenzimidazole STING receptor agonists with systemic activity.

Stimulator of interferon genes (STING) is a receptor in the endoplasmic reticulum that propagates innate immune sensing of cytosolic pathogen-derived and self DNA1. The development of compounds that modulate STING has recently been the focus of intense research for the treatment of cancer and infectious diseases and as vaccine adjuvants2. To our knowledge, current efforts are focused on the development of modified cyclic dinucleotides that mimic the endogenous STING ligand cGAMP; these have progressed into clinical trials in patients with solid accessible tumours amenable to intratumoral delivery3. Here we report the discovery of a small molecule STING agonist that is not a cyclic dinucleotide and is systemically efficacious for treating tumours in mice. We developed a linking strategy to synergize the effect of two symmetry-related amidobenzimidazole (ABZI)-based compounds to create linked ABZIs (diABZIs) with enhanced binding to STING and cellular function. Intravenous administration of a diABZI STING agonist to immunocompetent mice with established syngeneic colon tumours elicited strong anti-tumour activity, with complete and lasting regression of tumours. Our findings represent a milestone in the rapidly growing field of immune-modifying cancer therapies.

Allosteric Wip1 phosphatase inhibition through flap-subdomain interaction.

Although therapeutic interventions of signal-transduction cascades with targeted kinase inhibitors are a well-established strategy, drug-discovery efforts to identify targeted phosphatase inhibitors have proven challenging. Herein we report a series of allosteric, small-molecule inhibitors of wild-type p53-induced phosphatase (Wip1), an oncogenic phosphatase common to multiple cancers. Compound binding to Wip1 is dependent on a 'flap' subdomain located near the Wip1 catalytic site that renders Wip1 structurally divergent from other members of the protein phosphatase 2C (PP2C) family and that thereby confers selectivity for Wip1 over other phosphatases. Treatment of tumor cells with the inhibitor GSK2830371 increases phosphorylation of Wip1 substrates and causes growth inhibition in both hematopoietic tumor cell lines and Wip1-amplified breast tumor cells harboring wild-type TP53. Oral administration of Wip1 inhibitors in mice results in expected pharmacodynamic effects and causes inhibition of lymphoma xenograft growth. To our knowledge, GSK2830371 is the first orally active, allosteric inhibitor of Wip1 phosphatase.

The amino-acid substituents of dipeptide substrates of cathepsin C can determine the rate-limiting steps of catalysis.

We examined the cathepsin C-catalyzed hydrolysis of dipeptide substrates of the form Yaa-Xaa-AMC, using steady-state and pre-steady-state kinetic methods. The substrates group into three kinetic profiles based upon the broad range observed for k(cat)/K(a) and k(cat) values, pre-steady-state time courses, and solvent kinetic isotope effects (sKIEs). The dipeptide substrate Gly-Arg-AMC displayed large values for k(cat)/K(a) (1.6 ± 0.09 µM(-1) s(-1)) and k(cat) (255 ± 6 s(-1)), an inverse sKIE on k(cat)/K(a) ((D)(k(cat)/K(a)) = 0.6 ± 0.15), a modest, normal sKIE on k(cat) ((D)k(cat) = 1.6 ± 0.2), and immeasurable pre-steady-state kinetics, indicating an extremely fast pre-steady-state rate (>400 s(-1)). (Errors on fitted values are omitted in the text for clarity but may be found in Table 2.) These results conformed to a kinetic model where the acylation (k(ac)) and deacylation (k(dac)) half-reactions are very fast and similar in value. The second substrate type, Gly-Tyr-AMC and Ser-Tyr-AMC, the latter the subject of a comprehensive kinetic study (Schneck et al. (2008) Biochemistry 47, 8697-8710), were found to be less active substrates compared to Gly-Arg-AMC, with respective k(cat)/K(a) values of 0.49 ± 0.07 µM(-1 )s(-1) and 5.3 ± 0.5 µM(-1 )s(-1), and k(cat) values of 28 ± 1 s(-1) and 25 ± 0.5 s(-1). Solvent kinetic isotope effects for Ser-Tyr-AMC were found to be inverse for k(cat)/K(a) ((D)(k(cat)/K(a)) = 0.74 ± 0.05) and normal for k(cat) ((D)k(cat) = 2.3 ± 0.1) but unlike Gly-Arg-AMC, pre-steady-state kinetics of Gly-Tyr-AMC and Ser-Tyr-AMC were measurable and characterized by a single-exponential burst, with fast transient rates (490 s(-1) and 390 s(-1), respectively), from which it was determined that k(ac) ≫ k(dac) ∼ k(cat). The third substrate type, Gly-Ile-AMC, gave very low values of k(cat)/K(a) (0.0015 ± 0.0001 µM(-1) s(-1)) and k(cat) (0.33 ± 0.02 s(-1)), no sKIEs, ((D)(k(cat)/K(a)) = 1.05 ± 0.5 and (D)k(cat) = 1.06 ± 0.4), and pre-steady-state kinetics exhibited a discernible, but negligible, transient phase. For this third class of substrate, kinetic modeling was consistent with a mechanism in which k(dac) > k(ac) ∼ k(cat), and for which an isotope-insensitive step in the acylation half-reaction is the slowest. The combined results of these studies suggested that the identity of the amino acid at the P(1) position of the substrate is the main determinant of catalysis. On the basis of these kinetic data, together with crystallographic studies of substrate analogues and molecular dynamics analysis with models of acyl-enzyme intermediates, we present a catalytic model derived from the relative rates of the acylation vs deacylation half-reactions of cathepsin C. The chemical steps of catalysis are proposed to be dependent upon the conformational freedom of the amino acid substituents for optimal alignment for thiolation (acylation) or hydrolysis (deacylation). These studies suggest ideas for inhibitor design for papain-family cysteine proteases and strategies to progress drug discovery for other classes of disease-relevant cysteine proteases.

Druggability Assessment of Allosteric Proteins by Dynamics Simulations in the Presence of Probe Molecules.

Druggability assessment of a target protein has emerged in recent years as an important concept in hit-to-lead optimization. A reliable and physically relevant measure of druggability would allow informed decisions on the risk of investing in a particular target. Here, we define "druggability" as a quantitative estimate of binding sites and affinities for a potential drug acting on a specific protein target. In the present study, we describe a new methodology that successfully predicts the druggability and maximal binding affinity for a series of challenging targets, including those that function through allosteric mechanisms. Two distinguishing features of the methodology are (i) simulation of the binding dynamics of a diversity of probe molecules selected on the basis of an analysis of approved drugs and (ii) identification of druggable sites and estimation of corresponding binding affinities on the basis of an evaluation of the geometry and energetics of bound probe clusters. The use of the methodology for a variety of targets such as murine double mutant-2, protein tyrosine phosphatase 1B (PTP1B), lymphocyte function-associated antigen 1, vertebrate kinesin-5 (Eg5), and p38 mitogen-activated protein kinase provides examples for which the method correctly captures the location and binding affinities of known drugs. It also provides insights into novel druggable sites and the target's structural changes that would accommodate, if not promote and stabilize, drug binding. Notably, the ability to identify high affinity spots even in challenging cases such as PTP1B or Eg5 shows promise as a rational tool for assessing the druggability of protein targets and identifying allosteric or novel sites for drug binding.

Azepanone-based inhibitors of human cathepsin S: optimization of selectivity via the P2 substituent.

A series of azepanone inhibitors of cathepsin S is described. Selectivity over both cathepsin K and cathepsin L was achieved by varying the P2 substituent. Ultimately, a balanced potency and selectivity profile was achieved in compound 39 possessing a 1-methylcyclohexyl alanine at P2 and nicotinamide as the P' substituent. The cellular potency of selected analogs is also described.

Discovery of novel cyanamide-based inhibitors of cathepsin C.

The discovery of potent and selective cyanamide-based inhibitors of the cysteine protease cathepsin C is detailed. Optimization of the template with regard to plasma stability led to the identification of compound 17, a potent cathepsin C inhibitor with excellent selectivity over other cathepsins and potent in vivo activity in a cigarette smoke mouse model.

Molecular shape and medicinal chemistry: a perspective.

The eight contributions here provide ample evidence that shape as a volume or as a surface is a vibrant and useful concept when applied to drug discovery. It provides a reliable scaffold for "decoration" with chemical intuition (or bias) for virtual screening and lead optimization but also has its unadorned uses, as in library design, ligand fitting, pose prediction, or active site description. Computing power has facilitated this evolution by allowing shape to be handled precisely without the need to reduce down to point descriptors or approximate metrics, and the diversity of resultant applications argues for this being an important step forward. Certainly, it is encouraging that as computation has enabled our intuition, molecular shape has consistently surprised us in its usefulness and adaptability. The first Aurelius question, "What is the essence of a thing?", seems well answered, however, the third, "What do molecules do?", only partly so. Are the topics covered here exhaustive, or is there more to come? To date, there has been little published on the use of the volumetric definition of shape described here as a QSAR variable, for instance, in the prediction or classification of activity, although other shape definitions have been successful applied, for instance, as embodied in the Compass program described above in "Shape from Surfaces". Crystal packing is a phenomenon much desired to be understood. Although powerful models have been applied to the problem, to what degree is this dominated purely by the shape of a molecule? The shape comparison described here is typically of a global nature, and yet some importance must surely be placed on partial shape matching, just as the substructure matching of chemical graphs has proved useful. The approach of using surfaces, as described here, offers some flavor of this, as does the use of metrics that penalize volume mismatch less than the Tanimoto, e.g., Tversky measures. As yet, there is little to go on as to how useful a paradigm this will be because there is less software and fewer concrete results.Finally, the distance between molecular shapes, or between any shapes defined as volumes or surfaces, is a metric property in the mathematical sense of the word. As yet, there has been little, if any, application of this observation. We cannot know what new application to the design and discovery of pharmaceuticals may yet arise from the simple concept of molecular shape, but it is fair to say that the progress so far is impressive.

Antitumor activity of an allosteric inhibitor of centromere-associated protein-E.

Centromere-associated protein-E (CENP-E) is a kinetochore-associated mitotic kinesin that is thought to function as the key receptor responsible for mitotic checkpoint signal transduction after interaction with spindle microtubules. We have identified GSK923295, an allosteric inhibitor of CENP-E kinesin motor ATPase activity, and mapped the inhibitor binding site to a region similar to that bound by loop-5 inhibitors of the kinesin KSP/Eg5. Unlike these KSP inhibitors, which block release of ADP and destabilize motor-microtubule interaction, GSK923295 inhibited release of inorganic phosphate and stabilized CENP-E motor domain interaction with microtubules. Inhibition of CENP-E motor activity in cultured cells and tumor xenografts caused failure of metaphase chromosome alignment and induced mitotic arrest, indicating that tight binding of CENP-E to microtubules is insufficient to satisfy the mitotic checkpoint. Consistent with genetic studies in mice suggesting that decreased CENP-E function can have a tumor-suppressive effect, inhibition of CENP-E induced tumor cell apoptosis and tumor regression.

ATP-competitive inhibitors of the mitotic kinesin KSP that function via an allosteric mechanism.

The mitotic kinesin KSP (kinesin spindle protein, or Eg5) has an essential role in centrosome separation and formation of the bipolar mitotic spindle. Its exclusive involvement in the mitotic spindle of proliferating cells presents an opportunity for developing new anticancer agents with reduced side effects relative to antimitotics that target tubulin. Ispinesib is an allosteric small-molecule KSP inhibitor in phase 2 clinical trials. Mutations that attenuate ispinesib binding to KSP have been identified, which highlights the need for inhibitors that target different binding sites. We describe a new class of selective KSP inhibitors that are active against ispinesib-resistant forms of KSP. These ATP-competitive KSP inhibitors do not bind in the nucleotide binding pocket. Cumulative data from generation of resistant cells, site-directed mutagenesis and photo-affinity labeling suggest that they compete with ATP binding via a novel allosteric mechanism.

A critical assessment of docking programs and scoring functions.

Docking is a computational technique that samples conformations of small molecules in protein binding sites; scoring functions are used to assess which of these conformations best complements the protein binding site. An evaluation of 10 docking programs and 37 scoring functions was conducted against eight proteins of seven protein types for three tasks: binding mode prediction, virtual screening for lead identification, and rank-ordering by affinity for lead optimization. All of the docking programs were able to generate ligand conformations similar to crystallographically determined protein/ligand complex structures for at least one of the targets. However, scoring functions were less successful at distinguishing the crystallographic conformation from the set of docked poses. Docking programs identified active compounds from a pharmaceutically relevant pool of decoy compounds; however, no single program performed well for all of the targets. For prediction of compound affinity, none of the docking programs or scoring functions made a useful prediction of ligand binding affinity.

Structure activity relationships of 5-, 6-, and 7-methyl-substituted azepan-3-one cathepsin K inhibitors.

The syntheses, in vitro characterizations, and rat and monkey in vivo pharmacokinetic profiles of a series of 5-, 6-, and 7-methyl-substituted azepanone-based cathepsin K inhibitors are described. Depending on the particular regiochemical substitution and stereochemical configuration, methyl-substituted azepanones were identified that had widely varied cathepsin K inhibitory potency as well as pharmacokinetic properties compared to the 4S-parent azepanone analogue, 1 (human cathepsin K, K(i,app) = 0.16 nM, rat oral bioavailability = 42%, rat in vivo clearance = 49.2 mL/min/kg). Of particular note, the 4S-7-cis-methylazepanone analogue, 10, had a K(i,app) = 0.041 nM vs human cathepsin K and 89% oral bioavailability and an in vivo clearance rate of 19.5 mL/min/kg in the rat. Hypotheses that rationalize some of the observed characteristics of these closely related analogues have been made using X-ray crystallography and conformational analysis. These examples demonstrate the potential for modulation of pharmacological properties of cathepsin inhibitors by substituting the azepanone core. The high potency for inhibition of cathepsin K coupled with the favorable rat and monkey pharmacokinetic characteristics of compound 10, also known as SB-462795 or relacatib, has made it the subject of considerable in vivo evaluation for safety and efficacy as an inhibitor of excessive bone resorption in rat, monkey, and human studies, which will be reported elsewhere.

Crossover point between dialkoxy disulfides (ROSSOR) and thionosulfites ((RO)2S=S): prediction, synthesis, and structure.

Isomeric preference between cyclic dialkoxy disulfides and thionosulfites is governed by the ring size of the heterocycle. Rings smaller than seven atoms prefer the thionosulfite connectivity, whereas larger rings or acyclic analogues favor the unbranched dialkoxy disulfide structure. Density functional calculations were employed to predict the crossover point at which both constitutional isomers are of comparable stability. Follow-up synthesis provides the previously unknown eight-membered ring dialkoxy disulfide 14 and seven-membered ring thionosulfite 15 from the same reaction. X-ray crystallography for all but one of the reaction products and complementary NMR analysis furnishes insights into both solid-state and solution conformations. A long-standing issue regarding the concerted vs catalyzed isomerization pathway between XSSX and X(2)S=S has been addressed for X = RO and shown to be acid dependent.

Calculated conformer energies for organic molecules with multiple polar functionalities are method dependent: Taxol (case study).

BACKGROUND: Molecular mechanics (MM) and quantum chemical (QM) calculations are widely applied and powerful tools for the stereochemical and conformational investigations of molecules. The same methods have been extensively used to probe the conformational profile of Taxol (Figure 1) both in solution and at the beta-tubulin protein binding site. RESULTS: In the present work, the relative energies of seven conformations of Taxol derived from NMR and X-ray analyses were compared with a set of widely used force fields and semiempirical MO methods coupled to a continuum solvent treatment. The procedures not only diverge significantly in their assessment of relative conformational energies, but none of them provide satisfactory agreement with experiment. CONCLUSIONS: For Taxol, molecular mechanics and semiempirical QM methods are unable to provide a consistent energetic ranking of side-chain conformations. For similar highly polar organic structures, "energy-free" conformational search methods are advised.