Search

Adobe PDF - MCT-07-0149--Suppl_Data.pdf from Inhibition of the heat shock protein 90 molecular chaperone in vitro and in vivo by novel, synthetic, potent resorcinylic pyrazole/isoxazole amide analogues

Although the heat shock protein 90 (HSP90) inhibitor 17-allylamino-17-demethoxygeldanamycin (17-AAG) shows clinical promise, potential limitations encourage development of alternative chemotypes. We discovered the 3,4-diarylpyrazole resorcinol CCT018159 by high-throughput screening and used structure-based design to generate more potent pyrazole amide analogues, exemplified by VER-49009. Here, we describe the detailed biological properties of VER-49009 and the corresponding isoxazole VER-50589. X-ray crystallography showed a virtually identical HSP90 binding mode. However, the dissociation constant (Kd) of VER-50589 was 4.5 ± 2.2 nmol/L compared with 78.0 ± 10.4 nmol/L for VER-49009, attributable to higher enthalpy for VER-50589 binding. A competitive binding assay gave a lower IC50 of 21 ± 4 nmol/L for VER-50589 compared with 47 ± 9 nmol/L for VER-49009. Cellular uptake of VER-50589 was 4-fold greater than for VER-49009. Mean cellular antiproliferative GI50 values for VER-50589 and VER-49009 for a human cancer cell line panel were 78 ± 15 and 685 ± 119 nmol/L, respectively, showing a 9-fold potency gain for the isoxazole. Unlike 17-AAG, but as with CCT018159, cellular potency of these analogues was independent of NAD(P)H:quinone oxidoreductase 1/DT-diaphorase and P-glycoprotein expression. Consistent with HSP90 inhibition, VER-50589 and VER-49009 caused induction of HSP72 and HSP27 alongside depletion of client proteins, including C-RAF, B-RAF, and survivin, and the protein arginine methyltransferase PRMT5. Both caused cell cycle arrest and apoptosis. Extent and duration of pharmacodynamic changes in an orthotopic human ovarian carcinoma model confirmed the superiority of VER-50589 over VER-49009. VER-50589 accumulated in HCT116 human colon cancer xenografts at levels above the cellular GI50 for 24 h, resulting in 30% growth inhibition. The results indicate the therapeutic potential of the resorcinylic pyrazole/isoxazole amide analogues as HSP90 inhibitors. [Mol Cancer Ther 2007;6(4):1198–211]

The predictive power of simulation has become embedded in the infrastructure of modern economies. The scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations and how this might be achieved are discussed., The predictive power of simulation has become embedded in the infrastructure of modern economies. Computer-aided design is ubiquitous throughout industry. In aeronautical engineering, built infrastructure and materials manufacturing, simulations are routinely used to compute the performance of potential designs before construction. The ability to predict the behaviour of products is a driver of innovation by reducing the cost barrier to new designs, but also because radically novel ideas can be piloted with relatively little risk. Accurate weather forecasting is essential to guide domestic and military flight paths, and therefore the underpinning simulations are critical enough to have implications for national security. However, in the pharmaceutical and biotechnological industries, the application of computer simulations remains limited by the capabilities of the technology with respect to the complexity of molecular biology and human physiology. Over the last 30 years, molecular-modelling tools have gradually gained a degree of acceptance in the pharmaceutical industry. Drug discovery has begun to benefit from physics-based simulations. While such simulations have great potential for improved molecular design, much scepticism remains about their value. The motivations for such reservations in industry and areas where simulations show promise for efficiency gains in preclinical research are discussed. In this, the first of two complementary papers, the scientific and technical progress that needs to be made to improve the predictive power of biomolecular simulations, and how this might be achieved, is firstly discussed (Part 1). In Part 2, the status of computer simulations in pharma is contrasted with aerodynamics modelling and weather forecasting, and comments are made on the cultural changes needed for equivalent computational technologies to become integrated into life-science industries.

The serine/threonine kinase DYRK1A has been implicated in regulation of a variety of cellular processes associated with cancer progression, including cell cycle control, DNA damage repair, protection from apoptosis, cell differentiation, and metastasis. In addition, elevated-level DYRK1A activity has been associated with increased severity of symptoms in Down's syndrome. A selective inhibitor of DYRK1A could therefore be of therapeutic benefit. We have used fragment and structure-based discovery methods to identify a highly selective, well-tolerated, brain-penetrant DYRK1A inhibitor which showed in vivo activity in a tumor model. The inhibitor provides a useful tool compound for further exploration of the effect of DYRK1A inhibition in models of disease.

The intramolecular reorganization energy (ΔEReorg) of compounds upon binding to proteins is a component of the binding free energy, which has long received particular attention, for fundamental and practical reasons. Understanding ΔEReorg would benefit the science of molecular recognition and drug design. For instance, the tolerable strain energy of compounds upon binding has been elusive. Prior studies found some large ΔEReorg values (e.g. > 10 kcal/mol), received with skepticism since they imply excessive opposition to binding. Indeed, estimating ΔEReorg is technically difficult. Typically, ΔEReorg has been approached by taking two energy-minimized conformers representing the bound and unbound states, and subtracting their conformational energy. This is a drastic oversimplification, liable to conformational collapse of the unbound conformer. Instead, the present work applies extensive molecular dynamics (MD) and the modern OPLS3 force-field to simulate compounds bound and unbound states, in explicit solvent under physically relevant conditions. The thermalized unbound compounds populate multiple conformations, not reducible to one or a few energy-minimized conformers. The intramolecular energies in the bound and unbound states were averaged over pertinent conformational ensembles, and the reorganization enthalpy upon binding (ΔHReorg) deduced by subtraction. This was applied to 76 systems, including 43 approved drugs, carefully selected for i) the quality of the bioactive X-ray structures and ii) the diversity of the chemotypes, their properties and protein targets. It yielded comparatively low ΔHReorg values (median = 1.4 kcal/mol, mean = 3.0 kcal/mol). A new finding is the observation of negative ΔHReorg values. Indeed, reorganization energies do not have to oppose binding, e.g. when intramolecular interactions stabilize preferentially the bound state. Conversely, even with competing water molecules, intramolecular interactions can occur predominantly in the unbound compound, and be replaced by intermolecular counterparts upon protein binding. Such disruption of intramolecular interactions upon binding gives rise to occasional larger ΔHReorg values. Such counterintuitive larger ΔHReorg values may be rationalized as a redistribution of interactions upon binding, qualitatively compatible with binding.

The generation of 3D conformers of small molecules underpins most computational drug discovery. Thus, the conformer quality is critical and depends on their energetics. A key parameter is the empirical conformational energy window (ΔEw), since only conformers within ΔEw are retained. However, ΔEw values in use appear unrealistically large. We analyze the factors pertaining to the conformer energetics and ΔEw. We argue that more attention must be focused on the problem of collapsed low-energy conformers. That is due to artificial intramolecular stabilization and occurs even with continuum solvation. Consequently, the conformational energy of extended bioactive structures is artefactually increased, which inflates ΔEw. Thus, this Perspective highlights the issues arising from low-energy conformers and suggests improvements via empirical or physics-based strategies.

The kinase DYRK1A is an attractive target for drug discovery programs due to its implication in multiple diseases. Through a fragment screen, we identified a simple biaryl compound that is bound to the DYRK1A ATP site with very high efficiency, although with limited selectivity. Structure-guided optimization cycles enabled us to convert this fragment hit into potent and selective DYRK1A inhibitors. Exploiting the structural differences in DYRK1A and its close homologue DYRK2, we were able to fine-tune the selectivity of our inhibitors. Our best compounds potently inhibited DYRK1A in the cell culture and in vivo and demonstrated drug-like properties. The inhibition of DYRK1A in vivo translated into dose-dependent tumor growth inhibition in a model of ovarian carcinoma.

Most of eukaryotic cellular DNA is packed in nucleosome core particles (NCPs), in which the DNA (DNANCP) is wrapped around histones. The influence of this organization on the intrinsic local dynamics of DNA is largely unknown, in particular because capturing such information from experiments remains notoriously challenging. Given the importance of dynamical properties in DNA functions, we addressed this issue using CHARMM36 MD simulations of a nucleosome containing the NCP positioning 601 sequence and four related free dodecamers. Comparison between DNANCP and free DNA reveals a limited impact of the dense DNA-histone interface on correlated motions of dinucleotide constituents and on fluctuations of inter base pair parameters. A characteristic feature intimately associated with the DNANCP super-helical path is a set of structural periodicities that includes a marked alternation of regions enriched in backbone BI and BII conformers. This observation led to uncover a convincing correspondence between the sequence effect on BI/BII propensities in both DNANCP and free DNA, strengthening the idea that the histone preference for particular DNA sequences relies on those intrinsic structural properties. These results offer for the first time a detailed view of the DNA dynamical behavior within NCP. They show in particular that the DNANCP dynamics is substantial enough to preserve the ability to structurally adjust to external proteins, for instance remodelers. Also, fresh structural arguments highlight the relevance of relationships between DNA sequence and structural properties for NCP formation. Overall, our work offers a more rational framework to approach the functional, biological roles of NCP.

Drug discovery is increasingly tackling challenging protein binding sites regarding molecular recognition and druggability, including shallow and solvent-exposed protein-protein interaction interfaces. Macrocycles are emerging as promising chemotypes to modulate such sites. Despite their chemical complexity, macrocycles comprise important drugs and offer advantages compared to non-cyclic analogs, hence the recent impetus in the medicinal chemistry of macrocycles. Elaboration of macrocycles, or constituent fragments, can strongly benefit from knowledge of their binding mode to a target. When such information from X-ray crystallography is elusive, computational docking can provide working models. However, few studies have explored docking protocols for macrocycles, since conventional docking methods struggle with the conformational complexity of macrocycles, and also potentially with the shallower topology of their binding sites. Indeed, macrocycle binding mode prediction with the mainstream docking software GOLD has hardly been explored. Here, we present an in-depth study of macrocycle docking with GOLD and the ChemPLP scores. First, we summarize the thorough curation of a test set of 41 protein-macrocycle X-ray structures, raising the issue of lattice contacts with such systems. Rigid docking of the known bioactive conformers was successful (three top ranked poses) for 92.7% of the systems, in absence of crystallographic waters. Thus, without conformational search issues, scoring performed well. However, docking success dropped to 29.3% with the GOLD built-in conformational search. Yet, the success rate doubled to 58.5% when GOLD was supplied with extensive conformer ensembles docked rigidly. The reasons for failure, sampling or scoring, were analyzed, exemplified with particular cases. Overall, binding mode prediction of macrocycles remains challenging, but can be much improved with tailored protocols. The analysis of the interplay between conformational sampling and docking will be relevant to the prospective modelling of macrocycles in general.

International audience; Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (RDCs) involving C6/8-H6/8, C3-H3 and C4-H4 vectors are correlated to the 31 P chemical shifts (␦P), which reflect the populations of the BI and BII backbone states. The ␦Ps are also correlated to the internucleotide distances (D inter) involving H6/8, H2 and H2 protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied RDCs depend mostly on roll, while D inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how ␦P measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium , and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5 and 3 ends of the dode-camers also supplies new information on the fraying events, otherwise neglected.

We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI ↔ BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.

Cytochrome P450 enzymes (P450s) are important in drug metabolism and have been linked to adverse drug reactions. P450s display broad substrate reactivity, and prediction of metabolites is complex. QM/MM studies of P450 reactivity have provided insight into important details of the reaction mechanisms and have the potential to make predictions of metabolite formation. Here we present a comprehensive study of the oxidation of three widely used pharmaceutical compounds (S-ibuprofen, diclofenac, and S-warfarin) by one of the major drug-metabolizing P450 isoforms, CYP2C9. The reaction barriers to substrate oxidation by the iron-oxo species (Compound I) have been calculated at the B3LYP-D/CHARMM27 level for different possible metabolism sites for each drug, on multiple pathways. In the cases of ibuprofen and warfarin, the process with the lowest activation energy is consistent with the experimentally preferred metabolite. For diclofenac, the pathway leading to the experimentally observed metabolite is not the one with the lowest activation energy. This apparent inconsistency with experiment might be explained by the two very different binding modes involved in oxidation at the two competing positions. The carboxylate of diclofenac interacts strongly with the CYP2C9 Arg108 side chain in the transition state for formation of the observed metabolite-but not in that for the competing pathway. We compare reaction barriers calculated both in the presence and in the absence of the protein and observe a marked improvement in selectivity prediction ability upon inclusion of the protein for all of the substrates studied. The barriers calculated with the protein are generally higher than those calculated in the gas phase. This suggests that active-site residues surrounding the substrate play an important role in controlling selectivity in CYP2C9. The results show that inclusion of sampling (particularly) and dispersion effects is important in making accurate predictions of drug metabolism selectivity of P450s using QM/MM methods.

There has been an explosion of structural information for pharmaceutical compounds bound to biological targets, but the conformations and dynamics of compounds free in solution are poorly characterized, if at all. Yet, knowledge of the unbound state is essential to understand the fundamentals of molecular recognition, including the much debated conformational intramolecular reorganization energy of a compound upon binding (ΔEReorg). Also, dependable observation of the unbound compounds is important for ligand-based drug discovery, e.g. with pharmacophore modelling. Here, these questions are addressed with long (⩾0.5μs) state-of-the-art molecular dynamics (MD) simulations of 26 compounds (including 7 approved drugs) unbound in explicit solvent. These compounds were selected to be chemically diverse, with a range of flexibility, and good quality bioactive X-ray structures. The MD-simulated free compounds are compared to their bioactive structure and conformers generated with ad hoc sampling in vacuo or with implicit generalized Born (GB) aqueous solvation models. The GB conformational models clearly depart from those obtained in explicit solvent, and suffer from conformational collapse almost as severe as in vacuo. Thus, the global energy minima in vacuo or with GB are not suitable representations of the unbound state, which can instead be extensively sampled by MD simulations. Many, but not all, MD-simulated compounds displayed some structural similarity to their bioactive structure, supporting the notion of conformational pre-organization for binding. The ligand-protein complexes were also simulated in explicit solvent, to estimate ΔEReorg as an enthalpic difference ΔHReorg between the intramolecular energies in the bound and unbound states. This fresh approach yielded ΔHReorg values⩽6kcal/mol for 18 out of 26 compounds. For three particularly polar compounds 15⩽ΔHReorg⩽20kcal/mol, supporting the notion that ΔHReorg can be substantial. Those large ΔHReorg values correspond to a redistribution of electrostatic interactions upon binding. Overall, the study illustrates how MD simulations offer a promising avenue to characterize the unbound state of medicinal compounds.

Inhibitors of the Hsp90 molecular chaperone are showing promise as anti-cancer agents. Here we describe a series of 4-aryl-5-cyanopyrrolo[2,3-d]pyrimidine ATP competitive Hsp90 inhibitors that were identified following structure-driven optimization of purine hits revealed by NMR based screening of a proprietary fragment library. Ligand-Hsp90 X-ray structures combined with molecular modeling led to the rational displacement of a conserved water molecule leading to enhanced affinity for Hsp90 as measured by fluorescence polarization, isothermal titration calorimetry and surface plasmon resonance assays. This displacement was achieved with a nitrile group, presenting an example of efficient gain in binding affinity with minimal increase in molecular weight. Some compounds in this chemical series inhibit the proliferation of human cancer cell lines in vitro and cause depletion of oncogenic Hsp90 client proteins and concomitant elevation of the co-chaperone Hsp70. In addition, one compound was demonstrated to be orally bioavailable in the mouse. This work demonstrates the power of structure-based design for the rapid evolution of potent Hsp90 inhibitors and the importance of considering conserved water molecules in drug design.

The -Cys-X1-X2-Cys- active site motif is central to the function of enzymes of the thioredoxin superfamily, including glutaredoxins. Their chemistry depends on the lowered pK(a) of the N-terminal thiolate cysteine of the -Cys-X1-X2-Cys- sequence; therefore its structure, dynamics, and electrostatics matter. Much information about the glutaredoxin structures was obtained by nuclear magnetic resonance (NMR), yet these various NMR structures produced heterogeneous and discordant views of the -Cys-X1-X2-Cys- motifs. This study addresses these inconsistencies by a computational and experimental investigation of three diverse reduced -Cys-X1-X2-Cys- motifs, from human glutaredoxin 1 (hGrx1), Escherichia coli glutaredoxin 2 (EcGrx2), and T4 virus glutaredoxin (T4Grx). The NMR models do not account for the low pK(a) of the N-terminal cysteine. However, extensive investigations of the NMR conformers by simulations yielded consensus structures for the -Cys-X1-X2-Cys- motif, with well-defined orientations for the cysteines. pK(a) calculations indicated that the consensus structure stabilizes the thiolate by local hydrogen bonds. The consensus structures of EcGrx2 and T4Grx formed the basis for predicting low pK(a) values for their N-terminal cysteines. Subsequent experimental titrations showed that these pK(a) values are

Summary To survive hostile conditions, the bacterial pathogen Mycobacterium tuberculosis produces millimolar concentrations of mycothiol as a redox buffer against oxidative stress. The reductases that couple the reducing power of mycothiol to redox active proteins in the cell are not known. We report a novel mycothiol-dependent reductase (mycoredoxin-1) with a CGYC catalytic motif. With mycoredoxin-1 and mycothiol deletion strains of Mycobacterium smegmatis, we show that mycoredoxin-1 and mycothiol are involved in the protection against oxidative stress. Mycoredoxin-1 acts as an oxidoreductase exclusively linked to the mycothiol electron transfer pathway and it can reduce S-mycothiolated mixed disulphides. Moreover, we solved the solution structures of oxidized and reduced mycoredoxin-1, revealing a thioredoxin fold with a putative mycothiol-binding site. With HSQC snapshots during electron transport, we visualize the reduction of oxidized mycoredoxin-1 as a function of time and find that mycoredoxin-1 gets S-mycothiolated on its N-terminal nucleophilic cysteine. Mycoredoxin-1 has a redox potential of −218 mV and hydrogen bonding with neighbouring residues lowers the pKa of its N-terminal nucleophilic cysteine. Determination of the oxidized and reduced structures of mycoredoxin-1, better understanding of mycothiol-dependent reactions in general, will likely give new insights in how M. tuberculosis survives oxidative stress in human macrophages.

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.

With computational tools, it has become straightforward to obtain three-dimensional (3D) conformers for drug-like molecules, even for large numbers of compounds. This underpins many downstream molecular modeling activities. We summarize the principles behind the methods, explain why considering multiple conformers per compound is necessary, and discuss the tests used to assess the computed conformers. It leads us to suggest that conformational sampling per se is essentially a solved problem for small drug-like compounds, and that efforts should focus on improving the energy models. Drug Dev Res 72: 85–94, 2011. © 2010 Wiley-Liss, Inc.

Structure-based drug design could benefit greatly from computational methodologies that accurately predict the binding affinity of small compounds to target biomolecules. However, the current scoring functions used to rank compounds in virtual screens by docking are not sufficiently accurate to guide reliably the design of tight binding ligands. Thus, there is strong interest in methodologies based on molecular simulations of protein-ligand complexes which are perceived to be more accurate and, with advances in computing power, amenable to routine use. This report provides an overview of the technical details necessary to understand, execute and analyze binding free energy calculations, using free energy perturbation or thermodynamic integration methods. Examples of possible applications in structure-based drug design are discussed. Current methodological limitations are highlighted as well as a number of ongoing developments to improve the scope, reliability, and practicalities of free energy calculations. These efforts are paving the way for a more common use of free energy calculations in molecular design.

Virtual screening of a library of commercially available compounds versus the structure of Mycobacterium tuberculosis lumazine synthase identified 2-(2-oxo-1,2-dihydrobenzo[cd]indole-6-sulfonamido)acetic acid (9) as a possible lead compound. Compound 9 proved to be an effective inhibitor of M. tuberculosis lumazine synthase with a K(i) of 70microM. Lead optimization through replacement of the carboxymethylsulfonamide sidechain with sulfonamides substituted with alkyl phosphates led to a four-carbon phosphate 38 that displayed a moderate increase in enzyme inhibitory activity (K(i) 38microM). Molecular modeling based on known lumazine synthase/inhibitor crystal structures suggests that the main forces stabilizing the present benzindolone/enzyme complexes involve pi-pi stacking interactions with Trp27 and hydrogen bonding of the phosphates with Arg128, the backbone nitrogens of Gly85 and Gln86, and the side chain hydroxyl of Thr87.

On sustaining damage to their DNA, cells employ a sophisticated mechanism of detection and repair, termed the DNA damage response (DDR). As a critical component of the DDR and G2/M checkpoint, Chk1 kinase represents an attractive target for cancer therapy. We have utilized a structure-based drug design approach to identify and develop VER-250840, a novel, orally active inhibitor of the checkpoint kinase, Chk1. VER-250840 exhibited sub-nM potency against Chk1 kinase with exquisite selectivity over an extensive and diverse panel of kinases. In vitro, VER-250840 inhibited Chk1 autophosphorylation with an IC50 of 1.0 nM and increased the number of S-phase tumor cells staining positive for pan-nuclear γH2AX with an EC50 of 7 - 27 nM. Accumulated genomic DNA damage by Chk1 inhibition led to irreversible cell cycle arrest, inhibition of tumor cell proliferation, increased replication stress, and cell death in both 2D culture and multicellular tumor spheroids. In an in vivo A2058 tumor xenograft PD study, VER-250840 demonstrated rapid and sustained inhibition of Chk1 auto-phosphorylation within 30 minutes of oral administration. Doses of 10 mg/kg and higher PO resulted in greater than 90% inhibition of tumor pChk1 (S296) over 24 hours. In SKOV3 in vivo models, VER-250840 inhibited Chk1 auto-phosphorylation, modulated other biomarkers of replication stress and DNA damage, and exhibited moderate antitumor activity with minimal toxicity when administered orally on a 21-day once-daily schedule. Work is ongoing to further optimize in vivo efficacy. In conclusion, VER-250840 demonstrates potent and selective activity as a monotherapy both in vitro and in vivo. From these findings, further evaluation and optimization of this novel kinase inhibitor is justly merited. Citation Format: Joanne Wayne, Stephen Stokes, Nicolas Foloppe, Helen Browne, Teresa Brooks, Karen Benwell, Lisa Baker, Zoe Daniels, Andrea Fiumana, Christopher Graham, Alba Macias, Daniel Maddox, Sean McKenna, Christopher Northfield, Stuart Ray, Heather Simmonite, Emma Stefaniak, Paul Webb, Mike Wood, Andrew Massey. Identification and preclinical characterisation of VER-250840, a potent, selective Chk1 inhibitor with in vivo oral single-agent antitumor activity [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr B163.

B-DNA flexibility, crucial for DNA-protein recognition, is sequence dependent. Free DNA in solution would in principle be the best reference state to uncover the relation between base sequences and their intrinsic flexibility; however, this has long been hampered by a lack of suitable experimental data. We investigated this relationship by compiling and analyzing a large dataset of NMR (31)P chemical shifts in solution. These measurements reflect the BI--BII equilibrium in DNA, intimately correlated to helicoidal descriptors of the curvature, winding and groove dimensions. Comparing the ten complementary DNA dinucleotide steps indicates that some steps are much more flexible than others. This malleability is primarily controlled at the dinucleotide level, modulated by the tetranucleotide environment. Our analyses provide an experimental scale called TRX that quantifies the intrinsic flexibility of the ten dinucleotide steps in terms of Twist, Roll, and X-disp (base pair displacement). Applying the TRX scale to DNA sequences optimized for nucleosome formation reveals a 10 base-pair periodic alternation of stiff and flexible regions. Thus, DNA flexibility captured by the TRX scale is relevant to nucleosome formation, suggesting that this scale may be of general interest to better understand protein-DNA recognition.

Ligand-based virtual screening with a 3D pharmacophore led to the discovery of 30 novel, diverse and drug-like ligands of the human cannabinoid receptor 1 (hCB(1)). The pharmacophore was validated with a hit rate of 16%, binding selectivity versus hCB(2), and expected functional profiles. The discovered compounds provide new tools for exploring cannabinoid pharmacology.

// Andrew J. Massey 1 , Stephen Stokes 1 , Helen Browne 1 , Nicolas Foloppe 1 , Andrea Fiumana 1 , Simon Scrace 1,3 , Mandy Fallowfield 1 , Simon Bedford 1 , Paul Webb 1 , Lisa Baker 1 , Mark Christie 2 , Martin J. Drysdale 1,4 and Mike Wood 1 1 Vernalis Research, Granta Park, Cambridge, UK 2 Akranim Ltd, Maidstone, UK 3 Horizon Discovery, Cambridge Research Park, Waterbeach, Cambridge, UK 4 Cancer Research UK Beatson Institute, Garscube Estate, Bearsden, Glasgow, UK Correspondence to: Andrew J. Massey, email: // Keywords : Chk1, V158411, fragment, drug discovery, colon cancer Received : June 17, 2015 Accepted : September 14, 2015 Published : September 30, 2015 Abstract Chk1 kinase is a critical component of the DNA damage response checkpoint especially in cancer cells and targeting Chk1 is a potential therapeutic opportunity for potentiating the anti-tumor activity of DNA damaging chemotherapy drugs. Fragment elaboration by structure guided design was utilized to identify and develop a novel series of Chk1 inhibitors culminating in the identification of V158411, a potent ATP-competitive inhibitor of the Chk1 and Chk2 kinases. V158411 abrogated gemcitabine and camptothecin induced cell cycle checkpoints, resulting in the expected modulation of cell cycle proteins and increased cell death in cancer cells. V158411 potentiated the cytotoxicity of gemcitabine, cisplatin, SN38 and camptothecin in a variety of p53 deficient human tumor cell lines in vitro , p53 proficient cells were unaffected. In nude mice, V158411 showed minimal toxicity as a single agent and in combination with irinotecan. In tumor bearing animals, V158411 was detected at high levels in the tumor with a long elimination half-life; no pharmacologically significant in vivo drug-drug interactions with irinotecan were identified through analysis of the pharmacokinetic profiles. V158411 potentiated the anti-tumor activity of irinotecan in a variety of human colon tumor xenograft models without additional systemic toxicity. These results demonstrate the opportunity for combining V158411 with standard of care chemotherapeutic agents to potentiate the therapeutic efficacy of these agents without increasing their toxicity to normal cells. Thus, V158411 would warrant further clinical evaluation.

The accurate prediction of ligand-biopolymer binding affinities is of general interest to medicinal chemistry, as well as to the broader field of molecular recognition. The ability to predict computationally the thermodynamics of these molecular recognition processes has been relatively weak until recently, however, continued developments on several fronts are extending the scope of applicability of these methods. The rapid growth in the number of protein-ligand structures has initially led to the development of a range of empirical scoring functions based on relatively simple descriptions of intermolecular interactions. These methods have had some success in ranking binding affinities when tuned to particular protein systems or in rather qualitative estimates of molecular fit in fast docking calculations. However, they are too unreliable for more detailed, quantitative, assessment and comparison of binding affinities. Physics-based free energy calculations are in principle more general and have the potential to be significantly more accurate. These approaches have seen steady development over many years and rely on carefully calibrated molecular energy functions (force-fields), simulations of the systems with explicit solvent, and the coming-of-age of continuum solvation models. In addition to the initially developped Free Energy Perturbation (FEP) and Thermodynamic Integration (TI) methods, new approaches include the Molecular Mechanics-Poisson-Boltzmann Surface Area (MM-PBSA) and the Linear Interaction Energy (LIE) approaches. This review concentrates on MM-PBSA and LIE, and their variants. The routine application of these calculations is becoming possible because of enhanced computational hardware and the development of a range of computational chemistry tools. This review addresses: i) the basic principles behind free energy calculations ii) recent methodological advances iii) comparisons of predicted and experimentally determined affinities iv) the uncertainties and limitations of both the computational and experimental data v) areas where progress can be made vi) the practicality of applying the methods at the different stages of the drug discovery and optimization process.

Targeting RNA with small molecule drugs is an area of great potential for therapeutic treatment of infections and possibly genetic and autoimmune diseases. However, a mature set of precedents and established methodology is lacking. The physicochemical properties of RNA raise specific issues and obstacles to development, and contribute to explain the distinct characteristics of natural RNA ligands, including antibiotics. Yet, RNA-targeting strategies are being implemented to reinvigorate antibacterial discovery by using the ribosomal X-ray structures to modify known antibiotics. To exploit further these structures, we suggest the use of existing protein kinase-directed libraries of drug-like compounds to target the A-site of the bacterial ribosome, on the basis of a specific structural hypothesis.

The backbone states of B-DNA influence its helical parameters, groove dimensions, and overall curvature. Therefore, detection and fine characterization of these conformational states are desirable. Using routine NMR experiments on a nonlabeled B-DNA oligomer and analyzing high-resolution X-ray structures, we investigated the relationship between interproton distances and backbone conformational states. The three H2'i-H6/8i+1, H2' 'i-H6/8i+1, and H6/8i-H6/8i+1 sequential distances were found cross-correlated and linearly coupled to epsilon-zeta values in X-ray structures and 31P chemical shifts (deltaP) in NMR that reflect the interconversion between the backbone BI (epsilon-zeta0 degrees ) and BII (epsilon-zeta0 degrees) states. These relationships provide a detailed check of the NMR data consistency and the possibility to extend the set of restraints for structural refinement through various extrapolations. Furthermore, they allow translation of deltaP in terms of BI/BII ratios. Also, comparison of many published deltaP in solution to crystal data shows that the impact of sequence on the BI/BII propensities is similar in both environments and is therefore an intrinsic and general property of B-DNA. This quantification of the populations of BI and BII is of general interest because these sequence-dependent backbone states act on DNA overall structure, a key feature for DNA-protein-specific recognition.

Inhibition of the Chk1 kinase by small molecules binding to its active site is a strategy of great therapeutic interest for oncology. We report how computational modelling predicted the binding mode of ligands of special interest to the Chk1 ATP site, for representatives of an indazole series and debromohymenialdisine. These binding modes were subsequently confirmed by X-ray crystallography. The binding mode of a potent indazole derivative involves non-conventional C–H⋯O and N–H⋯π-aromatic interactions with the protein. These interactions are formed in a buried pocket at the periphery of the ATP-binding site, the importance of which has previously been overlooked for ligand design against Chk1. It is demonstrated that filling this pocket can confer ligands with dramatically enhanced affinity for Chk1. Structural arguments in conjunction with assay data explain why targeting this pocket is also advantageous for selective binding to Chk1. Structural overlays of known inhibitors complexed with Chk1 show that only the indazole series utilizes the pocket of interest. Therefore, the analysis presented here should prove helpful in guiding future structure-based ligand design efforts against Chk1.

The specific protonation state of amino acids is crucial for the physicochemical properties of proteins and their biological functions. These protonation states influence, for instance, properties related to hydrogen bonding, solubility, and folding. pKa calculations for proteins are, therefore, important and require, in principle, a specification of the most stable protonated and deprotonated forms of each titratable group. This is complicated by the existence of multiple tautomers, like the five neutral tautomers of the guanidine moiety in arginine. In this study, the compounds N-methyl-guanidine and N-ethyl-guanidine were used to model the charged and all neutral protonation states of the arginine side chain. The relative stabilities of all five neutral tautomers were investigated systematically for the first time, using quantum-mechanical calculations. These relative stabilities were obtained in vacuo, water and chloroform, by combining the quantum-mechanical calculations with a continuum solvation model. The water model was used to represent arginines exposed to an aqueous solution, whereas the chloroform model has a polarity representative of a protein core or a membrane. This allowed determining the relative pKa's associated with each neutral tautomer in these environments. A key result is that significant differences in stability are found between the neutral tautomers, in both water and chloroform. Some tautomers are consistently found to be the most stable. These findings will be helpful to refine pKa calculations in proteins.

We report the discovery, synthesis, and crystallographic binding mode of novel furanopyrimidine and pyrrolopyrimidine inhibitors of the Chk1 kinase, an oncology target. These inhibitors are synthetically tractable and inhibit Chk1 by competing for its ATP site. A chronological account allows an objective comparison of modeled compound docking modes to the subsequently obtained crystal structures. The comparison provides insights regarding the interpretation of modeling results, in relationship to the multiple reasonable docking modes which may be obtained in a kinase-ATP site. The crystal structures were used to guide medicinal chemistry efforts. This led to a thorough characterization of a pair of ligand-protein complexes which differ by a single hydrogen bond. An analysis indicates that this hydrogen bond is expected to contribute a fraction of the 10-fold change in binding affinity, adding a valuable observation to the debate about the energetic role of hydrogen bonding in molecular recognition.

The eight nucleoside constituents of nucleic acids were simulated for 50 ns in explicit water with molecular dynamics. This provides equilibrium populations of the torsional degrees of freedom, their kinetics of interconversion, their couplings, and how they are influenced by water. This is important, given that a full and quantitative characterization of the nucleosides in aqueous solution by experimental means has been elusive, despite immense efforts in that direction. It is with the anti/syn equilibrium that the simulations are most complementary to experiment, by accessing directly the influence of the sugar type, sugar pucker, and base on the anti/syn populations. The glycosidic torsion distributions in the anti conformation are strongly affected by water and depart from the corresponding X-ray modal values and the associated energy minima in vacuo. Water also preferentially stabilizes some sugar conformations, showing that potential energies in vacuo are not sufficient to understand the nucleosides. Deoxythymidine (but not other pyrimidines) significantly populates the syn orientation. Guanine favors the syn orientation more than adenine. The ribose favors the syn orientation significantly more than the deoxyribose. The NORTH pucker coexists with the syn conformers. A hydrogen bond is frequently formed between the 5'-OH group and the syn bases, despite competition by water. The rate of the anti/syn transitions with purines is on the nanosecond time scale, confirming a long held assumption underpinning the interpretation of ultrasonic relaxation studies. Therefore, our knowledge of the structure and dynamics of nucleosides in solvent is only limited by the accuracy of the potential used to simulate them, and it is shown that such simulations provide a distinct and unique test of nucleic acid force fields. This confirmed that the widely distributed CHARMM27 force field is, overall, well-balanced with a particularly good representation of the ribose. Specific improvements, however, are suggested for the deoxyribose and torsion gamma.

The need for novel antibiotics is widely recognized. A well validated target of antibiotics is the bacterial ribosome. Recent X-ray structures of the ribosome bound to antibiotics have shed new light on the binding sites of these antibiotics, providing fresh impetus for structure-based strategies aiming at identifying new ribosomal ligands. In that respect, the ribosomal decoding region of the aminoacyl-tRNA acceptor site (A-site) is of particular interest because oligonucleotide model systems of this site are available for crystallography, NMR and compound binding assays. This work presents how these different resources can be combined in a hierarchical screening strategy which has led to the identification of new A-site ligands. The approach exploits an X-ray structure of the A-site against which large and diverse libraries of compounds were computationally docked. The complementarity of the compounds to the A-site was assessed using a scoring function specifically calibrated for RNA targets. Starting from approximately 1 million compounds, the computational selection of candidate ligands allowed us to focus the experimental work on 129 compounds, 34 of which showed affinity for the A-site in a FRET-based binding assay. NMR experiments confirmed binding to the A-site for some compounds. For the most potent compound in the FRET assay, a tentative binding mode is suggested, which is compatible with the NMR data and the limited SAR in this series. Overall, the results validate the screening strategy.

We have implemented five drug-like filters, based on 1D and 2D molecular descriptors, and applied them to characterize the drug-like properties of commercially available chemical compounds. In addition to previously published filters (Lipinski and Veber), we implemented a filter for medicinal chemistry tractability based on lists of chemical features drawn up by a panel of medicinal chemists. A filter based on the modeling of aqueous solubility (1 microM) was derived in-house, as well as another based on the modeling of Caco-2 passive membrane permeability (10 nm/s). A library of 2.7 million compounds was collated from the 23 compound suppliers and analyzed with these filters, highlighting a tendency toward highly lipophilic compounds. The library contains 1.6 M unique structures, of which 37% (607,223) passed all five drug-like filters. None of the 23 suppliers provides all the members of the drug-like subset, emphasizing the benefit of considering compounds from various compound suppliers as a source of diversity for drug discovery.

The variety of cellular functions performed by proteins of the thioredoxin superfamily is made possible by the wide range of redox potential associated with their active site -Cys-X-X-Cys- motif. The determinants of these differences in redox potential are of considerable interest but are not well understood. E. coli Glutaredoxin 1 (Grx1) and 3 (Grx3) are important model systems with different redox properties, despite sharing the same -Cys-Pro-Tyr-Cys- motif, very similar overall structures, and 33% sequence identity. Very long molecular dynamics simulations (0.25 micros total) and electrostatic calculations provide a revised view of the reduced Grx1 active site, which now can be reconciled with biochemical and functional data. Comparison of this new model to Grx3 uncovers differences in the structure, dynamics, and electrostatics of these active sites. The influence of peripheral residues on the properties of the -Cys-X-X-Cys- motif is illustrated specifically with the effect of a Lys to Arg substitution.

Recent crystal structures of the bacterial ribosome have identified the complex molecular interactions involved in antibiotic-ribosome recognition. Insights into the binding of aminoglycosides, macrolides, tetracyclines and other antibiotics provide opportunities for computational, structure-based approaches to be used in the design of appropriate modifications to existing antibiotics as well as in the discovery of completely new drug classes.

The glycosyl torsion (chi) in nucleic acids has long been recognized to be a major determinant of their conformational properties. chi torsional energetics were systematically mapped in deoxyribonucleosides using high-level quantum mechanical methods, for north and south sugar puckers and with gamma in the g(+) and trans conformations. In all cases, the syn conformation is found higher in energy than the anti. When gamma is changed from g(+) to trans, the anti orientation of the base is strongly destabilized, and the energy difference and barrier between anti and syn are significantly decreased. The barrier between anti and syn in deoxyribonucleosides is found to be less than 10 kcal/mol and tends to be lower with purines than with pyrimidines. With gamma = g(+)/chi = anti, a south sugar yields a significantly broader energy well than a north sugar with no energy barrier between chi values typical of A or B DNA. Contrary to the prevailing view, the syn orientation is not more stable with south puckers than with north puckers. The syn conformation is significantly more energetically accessible with guanine than with adenine in 5-nucleotides but not in nucleosides. Analysis of nucleic acid crystal structures shows that gamma = trans/chi = anti is a minor but not negligible conformation. Overall, chi appears to be a very malleable structural parameter with the experimental chi distributions reflecting, to a large extent, the associated intrinsic torsional energetics.

Molecular dynamics simulations of E. coli glutaredoxin1 in water have been performed to relate the dynamical parameters and entropy obtained in NMR relaxation experiments, with results extracted from simulated trajectory data. NMR relaxation is the most widely used experimental method to obtain data on dynamics of proteins, but it is limited to relatively short timescales and to motions of backbone amides or in some cases (13)C-H vectors. By relating the experimental data to the all-atom picture obtained in molecular dynamics simulations, valuable insights on the interpretation of the experiment can be gained. We have estimated the internal dynamics and their timescales by calculating the generalized order parameters (O) for different time windows. We then calculate the quasiharmonic entropy (S) and compare it to the entropy calculated from the NMR-derived generalized order parameter of the amide vectors. Special emphasis is put on characterizing dynamics that are not expressed through the motions of the amide group. The NMR and MD methods suffer from complementary limitations, with NMR being restricted to local vectors and dynamics on a timescale determined by the rotational diffusion of the solute, while in simulations, it may be difficult to obtain sufficient sampling to ensure convergence of the results. We also evaluate the amount of sampling obtained with molecular dynamics simulations and how it is affected by the length of individual simulations, by clustering of the sampled conformations. We find that two structural turns act as hinges, allowing the α helix between them to undergo large, long timescale motions that cannot be detected in the time window of the NMR dipolar relaxation experiments. We also show that the entropy obtained from the amide vector does not account for correlated motions of adjacent residues. Finally, we show that the sampling in a total of 100 ns molecular dynamics simulation can be increased by around 50%, by dividing the trajectory into 10 replicas with different starting velocities.

The chemistry of active-site cysteine residues is central to the activity of thiol-disulfide oxidoreductases of the thioredoxin superfamily. In these reactions, a nucleophilic thiolate is required, but the associated pK(a) values differ vastly in the superfamily, from less than 4 in DsbA to greater than 7 in Trx. The factors that stabilize this thiolate are, however, not clearly established. The glutaredoxins (Grxs), which are members of this superfamily, contain a Cys-Pro-Tyr-Cys motif in their active site. In reduced Grxs, the pK(a) of the N-terminal active-site nucleophilic cysteine residue is lowered significantly, and the stabilization of the corresponding thiolate is expected to influence the redox potential of these enzymes. Here, we use a combination of long molecular dynamics (MD) simulations, pK(a) calculations, and experimental investigations to derive the structure and dynamics of the reduced active site from Escherichia coli Grx3, and investigate the factors that stabilize the thiolate. Several different MD simulations converged toward a consensus conformation for the active-site cysteine residues (Cys11 and Cys14), after a number of local conformational changes. Key features of the model were tested experimentally by measurement of NMR scalar coupling constants, and determination of pK(a) values of selected residues. The pK(a) values of the Grx3 active-site residues were calculated during the MD simulations, and support the underlying structural model. The structure of Grx3, in combination with the pK(a) calculations, indicate that the pK(a) of the N-terminal active-site cysteine residue in Grx3 is intermediate between that of its counterpart in DsbA and Trx. The pK(a) values in best agreement with experiment are obtained with a low (4) protein dielectric constant. The calculated pK(a) values fluctuate significantly in response to protein dynamics, which underscores the importance of the details of the underlying structures when calculating pK(a) values. The thiolate of Cys11 is stabilized primarily by direct hydrogen bonding with the amide protons of Tyr13 and Cys14 and the thiol proton of Cys14, rather than by long-range interactions from charged groups or from a helix macrodipole. From the comparison of reduced Grx3 with other members of the thioredoxin superfamily, a unifying theme for the structural basis of thiol pK(a) differences in this superfamily begins to emerge.

In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has provided new insights into the intrinsic conformational energetics of these compounds and its relationship to nucleic acids structure and dynamics. In this article we review recent developments and present new results. New data include comparison of various levels of theory on conformational properties of nucleic acid building blocks, calculations on the abasic sugar, known to occur in vivo in DNA, on the TA conformation of DNA observed in the complex with the TATA box binding protein, and on inosine. Tests of the Hartree–Fock (HF), second-order Moller–Plesset (MP2), and Density Functional Theory/Becke3, Lee, Yang and Par (DFT/B3LYP) levels of theory show the overall shape of backbone torsional energy profiles (for γ, e, and χ) to be similar for the different levels, though some systematic differences are identified between the MP2 and DFT/B3LYP profiles. The east pseudorotation energy barrier in deoxyribonucleosides is also sensitive to the level of theory, with the HF and DFT/B3LYP east barriers being significantly lower (∼2.5 kcal/mol) than the MP2 counterpart (∼4.0 kcal/mol). Additional calculations at various levels of theory suggest that the east barrier in deoxyribonucleosides is between 3.0 and 4.0 kcal/mol. In the abasic sugar, the west pseudorotation energy barrier is found to be slightly lower than the east barrier and the south pucker is favored more than in standard nucleosides. Results on the TA conformation suggest that, at the nucleoside level, this conformation is significantly destabilized relative to the global energy minimum, or relative to the A- and B-DNA conformations. Deoxyribocytosine would destabilize the TA conformation more than other bases relative to the A-DNA conformation, but not relative to the B-DNA conformation. © 2002 Wiley Periodicals, Inc. Biopoly (Nucleic Acid Sci) 61: 61–76, 2002; DOI 10.1002/bip.10047

Empirical force-field calculations on biological molecules represent an effective method to obtain atomic detail information on the relationship of their structure to their function. Results from those calculations depend on the quality of the force field. In this manuscript, optimization of the CHARMM27 all-atom empirical force field for nucleic acids is presented together with the resulting parameters. The optimization procedure is based on the reproduction of small molecule target data from both experimental and quantum mechanical studies and condensed phase structural properties of DNA and RNA. Via an iterative approach, the parameters were primarily optimized to reproduce macromolecular target data while maximizing agreement with small molecule target data. This approach is expected to ensure that the different contributions from the individual moieties in the nucleic acids are properly balanced to yield condensed phase properties of DNA and RNA, which are consistent with experiment. The quality of the presented force field in reproducing both crystal and solution properties are detailed in the present and an accompanying manuscript (MacKerell and Banavali, J Comput Chem, this issue). The resultant parameters represent the latest step in the continued development of the CHARMM all-atom biomolecular force field for proteins, lipids, and nucleic acids. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 86–104, 2000

The intrinsic conformational energetics of macromolecules are an important contribution to their structure and dynamics. This information, however, is lacking for most dihedrals about rotatable bonds in nucleic acids. To bridge this gap, high-level ab initio quantum mechanical calculations have been performed on a number of furanose-based model compounds designed to model the intrinsic torsional energetics of the dihedrals e, γ, β, and χ in deoxyribonucleic acid (DNA). Energy profiles have been obtained with furanose in both the C3‘ endo and C2‘ endo conformations, to model the torsions in the context of the two conformational ranges populated by the sugar in DNA. The resulting energy profiles are compared to the corresponding crystallographically determined population distributions, with which they can be reconciled to a large extent. This suggests that the models used to derive these energy surfaces are relevant to DNA. Torsion e is found to be intrinsically very flexible with the low-energy regions enc...

The present work analyzes the intrinsic conformational energetics associated with the puckering of the deoxyribose and ribose sugars in nucleic acids using high-level ab initio quantum mechanical calculations. A variety of model compounds have been designed to define the minimal structural unit suitable to model the sugar moiety in nucleic acids. Results suggest that all the structural features of a nucleoside are required to model the sugar moiety of nucleic acids. Stuctures calculated at the MP2 level of theory are in close agreement with experimental structural information. In deoxyribose, the south pucker (B form of double helices) is intrinsically favored over the north pucker (A form of double helices) by ∼1.0 kcal/mol. In contrast, for ribose, with torsion e in an RNA-like conformation, the north pucker is favored over the south pucker by ∼2.0 kcal/mol. For both the deoxyribose and ribose of nucleic acids, the lowest energy barrier between the north and south puckers is >4.0 kcal/mol. The present c...

Computational conformational sampling underpins much of molecular modeling and design in pharmaceutical work. The sampling of smaller drug-like compounds has been an active area of research. However, few studies have tested in details the sampling of larger more flexible compounds, which are also relevant to drug discovery, including therapeutic peptides, macrocycles, and inhibitors of protein–protein interactions. Here, we investigate extensively mainstream conformational sampling methods on three carefully curated compound sets, namely the ‘Drug-like’, larger ‘Flexible’, and ‘Macrocycle’ compounds. These test molecules are chemically diverse with reliable X-ray protein-bound bioactive structures. The compared sampling methods include Stochastic Search and the recent LowModeMD from MOE, all the low-mode based approaches from MacroModel, and MD/LLMOD recently developed for macrocycles. In addition to default settings, key parameters of the sampling protocols were explored. The performance of the computational protocols was assessed via (i) the reproduction of the X-ray bioactive structures, (ii) the size, coverage and diversity of the output conformational ensembles, (iii) the compactness/extendedness of the conformers, and (iv) the ability to locate the global energy minimum. The influence of the stochastic nature of the searches on the results was also examined. Much better results were obtained by adopting search parameters enhanced over the default settings, while maintaining computational tractability. In MOE, the recent LowModeMD emerged as the method of choice. Mixed torsional/low-mode from MacroModel performed as well as LowModeMD, and MD/LLMOD performed well for macrocycles. The low-mode based approaches yielded very encouraging results with the flexible and macrocycle sets. Thus, one can productively tackle the computational conformational search of larger flexible compounds for drug discovery, including macrocycles.

E. coli DsbA is an intensively studied enzyme of the thioredoxin superfamily of thiol-disulfide oxidoreductases. DsbA catalyzes the disulfide bond formation and folding of proteins in the bacterial periplasm. DsbA and its mutants have highlighted the strong and puzzling influence of the -C-X1-X2-C- active site variants, found across the thioredoxin superfamily, on the ionization and redox properties of this site. However, the interpretation of these observations remains wanting, largely due to a dearth of structural information. Here, molecular dynamics simulations are used to provide extensive information on the structure and dynamics of reduced -C30-X31-X32-C33- motifs in wild type DsbA and 13 of its mutants. These simulations are combined with calculations of the pK of H32 and of the very low pK of the catalytic cysteine C30. In wild type DsbA, the titrations of C30 and H32 are shown to be coupled; the protonation states and dynamics of H32 are examined. The thiolate of C30 is stabilized by hydrogen bonds with the protein. Modulation of these hydrogen bonds by alteration of residue X32 has the greatest impact on the pK of C30, which rationalizes its higher pK in thioredoxin and tryparedoxin. Because of structural constrains, residue X31 has only an indirect and weak influence on the pK of C30. The dynamics of C30 is clearly related to its stabilizing interactions and pK value. Although relatively small differences between pKs were not reproduced in the calculations, the major trends are explained, adding new insights to our understanding of enzymes in this family.