1
|
Stampolaki M, Stylianakis I, Zgurskaya HI, Kolocouris A. Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations. J Comput Aided Mol Des 2023; 37:245-264. [PMID: 37129848 DOI: 10.1007/s10822-023-00504-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Accepted: 04/03/2023] [Indexed: 05/03/2023]
Abstract
N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol-1 compared to the experimental values.
Collapse
Affiliation(s)
- Marianna Stampolaki
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece
- Department of NMR-Based Structural Biology, Max Planck Institute for Multidisciplinary Sciences, Am Faßberg 11, 37077, Göttingen, Germany
| | - Ioannis Stylianakis
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece
| | - Helen I Zgurskaya
- Department of Chemistry and Biochemistry, University of Oklahoma, Stephenson Life Sciences Research Center, 101 Stephenson Parkway, Norman, OK, 73019-5251, USA
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece.
| |
Collapse
|
2
|
Boittier ED, Burns JM, Gandhi NS, Ferro V. GlycoTorch Vina: Docking Designed and Tested for Glycosaminoglycans. J Chem Inf Model 2020; 60:6328-6343. [PMID: 33152249 DOI: 10.1021/acs.jcim.0c00373] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Glycosaminoglycans (GAGs) are a family of anionic carbohydrates that play an essential role in the physiology and pathology of all eukaryotic life forms. Experimental determination of GAG-protein complexes is challenging due to their difficult isolation from biological sources, natural heterogeneity, and conformational flexibility-including possible ring puckering of sulfated iduronic acid from 1C4 to 2SO conformation. To overcome these challenges, we present GlycoTorch Vina (GTV), a molecular docking tool based on the carbohydrate docking program VinaCarb (VC). Our program is unique in that it contains parameters to model 2SO sugars while also supporting glycosidic linkages specific to GAGs. We discuss how crystallographic models of carbohydrates can be biased by the choice of refinement software and structural dictionaries. To overcome these variations, we carefully curated 12 of the best available GAG and GAG-like crystal structures (ranging from tetra- to octasaccharides or longer) obtained from the PDB-REDO server and refined using the same protocol. Both GTV and VC produced pose predictions with a mean root-mean-square deviation (RMSD) of 3.1 Å from the native crystal structure-a statistically significant improvement when compared to AutoDock Vina (4.5 Å) and the commercial software Glide (5.9 Å). Examples of how real-space correlation coefficients can be used to better assess the accuracy of docking pose predictions are given. Comparisons between statistical distributions of empirical "salt bridge" interactions, relevant to GAGs, were compared to density functional theory (DFT) studies of model salt bridges, and water-mediated salt bridges; however, there was generally a poor agreement between these data. Water bridges appear to play an important, yet poorly understood, role in the structures of GAG-protein complexes. To aid in the rapid prototyping of future pose scoring functions, we include a module that allows users to include their own torsional and nonbonded parameters.
Collapse
Affiliation(s)
- Eric D Boittier
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Jed M Burns
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Neha S Gandhi
- Chemistry and Physics, Centre for Genomics and Personalised Health, Faculty of Science and Engineering, Queensland University of Technology, Brisbane, Queensland 4000, Australia
| | - Vito Ferro
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland 4072, Australia.,Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland 4072, Australia
| |
Collapse
|
3
|
Borbulevych O, Martin RI, Westerhoff LM. High-throughput quantum-mechanics/molecular-mechanics (ONIOM) macromolecular crystallographic refinement with PHENIX/DivCon: the impact of mixed Hamiltonian methods on ligand and protein structure. Acta Crystallogr D Struct Biol 2018; 74:1063-1077. [PMID: 30387765 PMCID: PMC6213575 DOI: 10.1107/s2059798318012913] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2018] [Accepted: 09/12/2018] [Indexed: 12/28/2022] Open
Abstract
Conventional macromolecular crystallographic refinement relies on often dubious stereochemical restraints, the preparation of which often requires human validation for unusual species, and on rudimentary energy functionals that are devoid of nonbonding effects owing to electrostatics, polarization, charge transfer or even hydrogen bonding. While this approach has served the crystallographic community for decades, as structure-based drug design/discovery (SBDD) has grown in prominence it has become clear that these conventional methods are less rigorous than they need to be in order to produce properly predictive protein-ligand models, and that the human intervention that is required to successfully treat ligands and other unusual chemistries found in SBDD often precludes high-throughput, automated refinement. Recently, plugins to the Python-based Hierarchical ENvironment for Integrated Xtallography (PHENIX) crystallographic platform have been developed to augment conventional methods with the in situ use of quantum mechanics (QM) applied to ligand(s) along with the surrounding active site(s) at each step of refinement [Borbulevych et al. (2014), Acta Cryst D70, 1233-1247]. This method (Region-QM) significantly increases the accuracy of the X-ray refinement process, and this approach is now used, coupled with experimental density, to accurately determine protonation states, binding modes, ring-flip states, water positions and so on. In the present work, this approach is expanded to include a more rigorous treatment of the entire structure, including the ligand(s), the associated active site(s) and the entire protein, using a fully automated, mixed quantum-mechanics/molecular-mechanics (QM/MM) Hamiltonian recently implemented in the DivCon package. This approach was validated through the automatic treatment of a population of 80 protein-ligand structures chosen from the Astex Diverse Set. Across the entire population, this method results in an average 3.5-fold reduction in ligand strain and a 4.5-fold improvement in MolProbity clashscore, as well as improvements in Ramachandran and rotamer outlier analyses. Overall, these results demonstrate that the use of a structure-wide QM/MM Hamiltonian exhibits improvements in the local structural chemistry of the ligand similar to Region-QM refinement but with significant improvements in the overall structure beyond the active site.
Collapse
Affiliation(s)
- Oleg Borbulevych
- QuantumBio Inc., 2790 West College Avenue, State College, PA 16801, USA
| | - Roger I. Martin
- QuantumBio Inc., 2790 West College Avenue, State College, PA 16801, USA
| | | |
Collapse
|
4
|
Malde AK, Stroet M, Caron B, Visscher KM, Mark AE. Predicting the Prevalence of Alternative Warfarin Tautomers in Solution. J Chem Theory Comput 2018; 14:4405-4415. [PMID: 29999318 DOI: 10.1021/acs.jctc.8b00453] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined. It is shown that varying the method used to estimate the heat of formation in vacuum (direct or via homodesmic reactions), whether entropic corrections were included, and the method used to estimate the free enthalpy of solvation (i.e., PCM, COSMO, or SMD implicit models or explicit solvent) lead to large differences in the predicted rank and relative populations of the tautomers. In this case, only a combination of the enthalpy of formation using homodesmic reactions and explicit solvent to estimate the free enthalpy of solvation yielded results compatible with the available experimental data. The work also suggests that a small but significant subset of the possible Warfarin tautomers are likely to be physiologically relevant.
Collapse
Affiliation(s)
- Alpeshkumar K Malde
- School of Chemistry and Molecular Biosciences , The University of Queensland , Brisbane , QLD 4072 , Australia
| | - Martin Stroet
- School of Chemistry and Molecular Biosciences , The University of Queensland , Brisbane , QLD 4072 , Australia
| | - Bertrand Caron
- School of Chemistry and Molecular Biosciences , The University of Queensland , Brisbane , QLD 4072 , Australia
| | - Koen M Visscher
- Division of Molecular Toxicology , VU University , Amsterdam , The Netherlands
| | - Alan E Mark
- School of Chemistry and Molecular Biosciences , The University of Queensland , Brisbane , QLD 4072 , Australia.,Institute for Molecular Bioscience , The University of Queensland , Brisbane , QLD 4072 , Australia
| |
Collapse
|
5
|
Ariz-Extreme I, Hub JS. Assigning crystallographic electron densities with free energy calculations-The case of the fluoride channel Fluc. PLoS One 2018; 13:e0196751. [PMID: 29771936 PMCID: PMC5957342 DOI: 10.1371/journal.pone.0196751] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Accepted: 04/18/2018] [Indexed: 11/25/2022] Open
Abstract
Approximately 90% of the structures in the Protein Data Bank (PDB) were obtained by X-ray crystallography or electron microscopy. Whereas the overall quality of structure is considered high, thanks to a wide range of tools for structure validation, uncertainties may arise from density maps of small molecules, such as organic ligands, ions or water, which are non-covalently bound to the biomolecules. Even with some experience and chemical intuition, the assignment of such disconnected electron densities is often far from obvious. In this study, we suggest the use of molecular dynamics (MD) simulations and free energy calculations, which are well-established computational methods, to aid in the assignment of ambiguous disconnected electron densities. Specifically, estimates of (i) relative binding affinities, for instance between an ion and water, (ii) absolute binding free energies, i.e., free energies for transferring a solute from bulk solvent to a binding site, and (iii) stability assessments during equilibrium simulations may reveal the most plausible assignments. We illustrate this strategy using the crystal structure of the fluoride specific channel (Fluc), which contains five disconnected electron densities previously interpreted as four fluoride and one sodium ion. The simulations support the assignment of the sodium ion. In contrast, calculations of relative and absolute binding free energies as well as stability assessments during free MD simulations suggest that four of the densities represent water molecules instead of fluoride. The assignment of water is compatible with the loss of these densities in the non-conductive F82I/F85I mutant of Fluc. We critically discuss the role of the ion force fields for the calculations presented here. Overall, these findings indicate that MD simulations and free energy calculations are helpful tools for modeling water and ions into crystallographic density maps.
Collapse
Affiliation(s)
- Igor Ariz-Extreme
- Institute for Microbiology and Genetics, University of Goettingen, Göttingen, Germany
| | - Jochen S. Hub
- Institute for Microbiology and Genetics, University of Goettingen, Göttingen, Germany
- * E-mail:
| |
Collapse
|
6
|
Peach ML, Cachau RE, Nicklaus MC. Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding. J Mol Recognit 2017; 30:10.1002/jmr.2618. [PMID: 28233410 PMCID: PMC5553890 DOI: 10.1002/jmr.2618] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2016] [Revised: 01/31/2017] [Accepted: 01/31/2017] [Indexed: 12/25/2022]
Abstract
In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.
Collapse
Affiliation(s)
- Megan L Peach
- Basic Science Program, Chemical Biology Laboratory, Leidos Biomedical Research Inc., Frederick National Laboratory for Cancer Research, Frederick, MD, USA
| | - Raul E Cachau
- Data Science and Information Technology Program, Advanced Biomedical Computing Center, Leidos Biomedical Research Inc., Frederick National Laboratory for Cancer Research, Frederick, MD, USA
| | - Marc C Nicklaus
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, MD, USA
| |
Collapse
|
7
|
Kim RR, Malde AK, Nematollahi A, Scott KF, Church WB. Molecular dynamics simulations reveal structural insights into inhibitor binding modes and functionality in human Group IIA phospholipase A 2. Proteins 2017; 85:827-842. [PMID: 28056488 DOI: 10.1002/prot.25235] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2016] [Revised: 12/14/2016] [Accepted: 12/21/2016] [Indexed: 11/09/2022]
Abstract
Human Group IIA phospholipase A2 (hGIIA) promotes inflammation in immune-mediated pathologies by regulating the arachidonic acid pathway through both catalysis-dependent and -independent mechanisms. The hGIIA crystal structure, both alone and inhibitor-bound, together with structures of closely related snake-venom-derived secreted phospholipase enzymes has been well described. However, differentiation of biological and nonbiological contacts and the relevance of structures determined from snake venom enzymes to human enzymes are not clear. We employed molecular dynamics (MD) and docking approaches to understand the binding of inhibitors that selectively or nonselectively block the catalysis-independent mechanism of hGIIA. Our results indicate that hGIIA behaves as a monomer in the solution environment rather than a dimer arrangement that is in the asymmetric unit of some crystal structures. The binding mode of a nonselective inhibitor, KH064, was validated by a combination of the experimental electron density and MD simulations. The binding mode of the selective pentapeptide inhibitor FLSYK to hGIIA was stipulated to be different to that of the snake venom phospholipases A2 of Daboia russelli pulchella (svPLA2 ). Our data suggest that the application of MD approaches to crystal structure data is beneficial in evaluating the robustness of conclusions drawn based on crystal structure data alone. Proteins 2017; 85:827-842. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Ryung Rae Kim
- Group in Biomolecular Structure and Informatics, Faculty of Pharmacy, The University of Sydney, Sydney, NSW, 2006, Australia
| | - Alpeshkumar K Malde
- School of Chemistry and Molecular Biosciences, University of Queensland, Brisbane, QLD, Australia, 4072
| | - Alireza Nematollahi
- Group in Biomolecular Structure and Informatics, Faculty of Pharmacy, The University of Sydney, Sydney, NSW, 2006, Australia
| | - Kieran F Scott
- School of Medicine, Western Sydney University, The Ingham Institute for Applied Medical Research, Liverpool, NSW, Australia, 2170.,Centre for Oncology Education and Research Translation (CONCERT), The Ingham Institute for Applied Medical Research, Liverpool, NSW, Australia, 2170
| | - W Bret Church
- Group in Biomolecular Structure and Informatics, Faculty of Pharmacy, The University of Sydney, Sydney, NSW, 2006, Australia
| |
Collapse
|
8
|
Adams PD, Aertgeerts K, Bauer C, Bell JA, Berman HM, Bhat TN, Blaney JM, Bolton E, Bricogne G, Brown D, Burley SK, Case DA, Clark KL, Darden T, Emsley P, Feher VA, Feng Z, Groom CR, Harris SF, Hendle J, Holder T, Joachimiak A, Kleywegt GJ, Krojer T, Marcotrigiano J, Mark AE, Markley JL, Miller M, Minor W, Montelione GT, Murshudov G, Nakagawa A, Nakamura H, Nicholls A, Nicklaus M, Nolte RT, Padyana AK, Peishoff CE, Pieniazek S, Read RJ, Shao C, Sheriff S, Smart O, Soisson S, Spurlino J, Stouch T, Svobodova R, Tempel W, Terwilliger TC, Tronrud D, Velankar S, Ward SC, Warren GL, Westbrook JD, Williams P, Yang H, Young J. Outcome of the First wwPDB/CCDC/D3R Ligand Validation Workshop. Structure 2016; 24:502-508. [PMID: 27050687 DOI: 10.1016/j.str.2016.02.017] [Citation(s) in RCA: 53] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2015] [Revised: 02/24/2016] [Accepted: 02/25/2016] [Indexed: 10/22/2022]
Abstract
Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent an important source of information concerning drug-target interactions, providing atomic level insights into the physical chemistry of complex formation between macromolecules and ligands. Of the more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ∼75% include at least one non-polymeric ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30-31, 2015. Experts in protein crystallography from academe and industry came together with non-profit and for-profit software providers for crystallography and with experts in computational chemistry and data archiving to discuss and make recommendations on best practices, as framed by a series of questions central to structural studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived in the PDB? How should the ligands be best represented? How should structural models of macromolecule-ligand complexes be validated? What supplementary information should accompany publications of structural studies of biological macromolecules? Consensus recommendations on best practices developed in response to each of these questions are provided, together with some details regarding implementation. Important issues addressed but not resolved at the workshop are also enumerated.
Collapse
Affiliation(s)
- Paul D Adams
- Molecular Biophysics & Integrated Bioimaging Division, Lawrence Berkeley Laboratory, Department of Bioengineering, UC Berkeley, Berkeley, CA 94720-8235, USA
| | | | - Cary Bauer
- Bruker AXS, Inc., Madison, WI 53711, USA
| | | | - Helen M Berman
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Talapady N Bhat
- Biosystems and Biomaterials Division, NIST, Gaithersburg, MD 20899, USA
| | | | - Evan Bolton
- National Center for Biotechnology Information, U.S. National Library of Medicine, Bethesda, MD 20894, USA
| | | | - David Brown
- School of Biosciences, University of Kent, Canterbury CT2 7NH, UK; Charles River Ltd., Structural Biology and Biophysics, Cambridge CB10 1XL, UK
| | - Stephen K Burley
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Skaggs School of Pharmacy and Pharmaceutical Sciences and San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA 92093, USA.
| | - David A Case
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Kirk L Clark
- Novartis Institutes for BioMedical Research, Cambridge, MA 02139, USA
| | - Tom Darden
- OpenEye Scientific, Cambridge, MA 02142, USA
| | - Paul Emsley
- MRC Laboratory of Molecular Biology, Cambridge CB2 0QH, UK
| | - Victoria A Feher
- Drug Design Data Resource and Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, CA 92093, USA.
| | - Zukang Feng
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Colin R Groom
- Cambridge Crystallographic Data Centre, Cambridge CB2 1EZ, UK.
| | | | - Jorg Hendle
- Structural Biology, Lilly Biotechnology Center, San Diego, CA 92121, USA
| | | | - Andrzej Joachimiak
- Structural Biology Center, Biosciences, Argonne National Laboratory, Argonne, IL 60439, USA
| | - Gerard J Kleywegt
- Protein Data Bank in Europe, European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Tobias Krojer
- Structural Genomics Consortium, University of Oxford, Oxford OX3 7DQ, UK
| | - Joseph Marcotrigiano
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Center for Advanced Biotechnology and Medicine, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Alan E Mark
- School of Chemistry & Molecular Biosciences, University of Queensland, St Lucia, QLD 4072, Australia
| | - John L Markley
- BioMagResBank, Department of Biochemistry, University of Wisconsin-Madison, Madison, WI 53706-1544, USA
| | - Matthew Miller
- Center for Advanced Biotechnology and Medicine, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Wladek Minor
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA 22908, USA
| | - Gaetano T Montelione
- Center for Advanced Biotechnology and Medicine, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Molecular Biology and Biochemistry, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | | | - Atsushi Nakagawa
- Protein Data Bank Japan, Institute for Protein Research, Osaka University, Osaka 565-0871, Japan
| | - Haruki Nakamura
- Protein Data Bank Japan, Institute for Protein Research, Osaka University, Osaka 565-0871, Japan
| | | | - Marc Nicklaus
- Computer-Aided Drug Design Group, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Frederick, MD 21702, USA
| | | | | | | | - Susan Pieniazek
- Bristol-Myers Squibb Research and Development, Pennington, NJ 08534, USA
| | - Randy J Read
- Department of Haematology, Cambridge Institute for Medical Research, University of Cambridge, Cambridge CB2 0XY, UK
| | - Chenghua Shao
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Steven Sheriff
- Bristol-Myers Squibb Research and Development, Princeton, NJ 08543, USA
| | - Oliver Smart
- Protein Data Bank in Europe, European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | | | - John Spurlino
- Janssen Pharmaceuticals, Inc., Spring House, PA 19002, USA
| | - Terry Stouch
- Science For Solutions, LLC, West Windsor, NJ 08550, USA
| | - Radka Svobodova
- CEITEC-Central European Institute of Technology and National Centre for Biomolecular Research, Masaryk University Brno, 625 00 Brno, Czech Republic
| | - Wolfram Tempel
- Structural Genomics Consortium, University of Toronto, Toronto, ON M5G 1L7, Canada
| | | | - Dale Tronrud
- Department of Biochemistry and Biophysics, Oregon State University, Corvallis, OR 97331, USA
| | - Sameer Velankar
- Protein Data Bank in Europe, European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Suzanna C Ward
- Cambridge Crystallographic Data Centre, Cambridge CB2 1EZ, UK
| | | | - John D Westbrook
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | | | - Huanwang Yang
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| | - Jasmine Young
- Research Collaboratory for Structural Bioinformatics Protein Data Bank, Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA; Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, USA
| |
Collapse
|
9
|
Borbulevych O, Martin RI, Tickle IJ, Westerhoff LM. XModeScore: a novel method for accurate protonation/tautomer-state determination using quantum-mechanically driven macromolecular X-ray crystallographic refinement. Acta Crystallogr D Struct Biol 2016; 72:586-98. [PMID: 27050137 PMCID: PMC4822566 DOI: 10.1107/s2059798316002837] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2015] [Accepted: 02/17/2016] [Indexed: 11/16/2022] Open
Abstract
Gaining an understanding of the protein-ligand complex structure along with the proper protonation and explicit solvent effects can be important in obtaining meaningful results in structure-guided drug discovery and structure-based drug discovery. Unfortunately, protonation and tautomerism are difficult to establish with conventional methods because of difficulties in the experimental detection of H atoms owing to the well known limitations of X-ray crystallography. In the present work, it is demonstrated that semiempirical, quantum-mechanics-based macromolecular crystallographic refinement is sensitive to the choice of a protonation-state/tautomer form of ligands and residues, and can therefore be used to explore potential states. A novel scoring method, called XModeScore, is described which enumerates the possible protomeric/tautomeric modes, refines each mode against X-ray diffraction data with the semiempirical quantum-mechanics (PM6) Hamiltonian and scores each mode using a combination of energetic strain (or ligand strain) and rigorous statistical analysis of the difference electron-density distribution. It is shown that using XModeScore it is possible to consistently distinguish the correct bound protomeric/tautomeric modes based on routine X-ray data, even at lower resolutions of around 3 Å. These X-ray results are compared with the results obtained from much more expensive and laborious neutron diffraction studies for three different examples: tautomerism in the acetazolamide ligand of human carbonic anhydrase II (PDB entries 3hs4 and 4k0s), tautomerism in the 8HX ligand of urate oxidase (PDB entries 4n9s and 4n9m) and the protonation states of the catalytic aspartic acid found within the active site of an aspartic protease (PDB entry 2jjj). In each case, XModeScore applied to the X-ray diffraction data is able to determine the correct protonation state as defined by the neutron diffraction data. The impact of QM-based refinement versus conventional refinement on XModeScore is also discussed.
Collapse
Affiliation(s)
- Oleg Borbulevych
- QuantumBio Inc., 2790 West College Avenue, State College, PA 16801, USA
| | - Roger I. Martin
- QuantumBio Inc., 2790 West College Avenue, State College, PA 16801, USA
| | - Ian J. Tickle
- Astex Pharmaceuticals, 436 Science Park, Milton Road, Cambridge CB4 0QA, England
| | | |
Collapse
|
10
|
|
11
|
Smart OS, Bricogne G. Achieving High Quality Ligand Chemistry in Protein-Ligand Crystal Structures for Drug Design. MULTIFACETED ROLES OF CRYSTALLOGRAPHY IN MODERN DRUG DISCOVERY 2015. [DOI: 10.1007/978-94-017-9719-1_13] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
12
|
Abstract
INTRODUCTION X-ray crystallography plays an important role in structure-based drug design (SBDD), and accurate analysis of crystal structures of target macromolecules and macromolecule-ligand complexes is critical at all stages. However, whereas there has been significant progress in improving methods of structural biology, particularly in X-ray crystallography, corresponding progress in the development of computational methods (such as in silico high-throughput screening) is still on the horizon. Crystal structures can be overinterpreted and thus bias hypotheses and follow-up experiments. As in any experimental science, the models of macromolecular structures derived from X-ray diffraction data have their limitations, which need to be critically evaluated and well understood for structure-based drug discovery. AREAS COVERED This review describes how the validity, accuracy and precision of a protein or nucleic acid structure determined by X-ray crystallography can be evaluated from three different perspectives: i) the nature of the diffraction experiment; ii) the interpretation of an electron density map; and iii) the interpretation of the structural model in terms of function and mechanism. The strategies to optimally exploit a macromolecular structure are also discussed in the context of 'Big Data' analysis, biochemical experimental design and structure-based drug discovery. EXPERT OPINION Although X-ray crystallography is one of the most detailed 'microscopes' available today for examining macromolecular structures, the authors would like to re-emphasize that such structures are only simplified models of the target macromolecules. The authors also wish to reinforce the idea that a structure should not be thought of as a set of precise coordinates but rather as a framework for generating hypotheses to be explored. Numerous biochemical and biophysical experiments, including new diffraction experiments, can and should be performed to verify or falsify these hypotheses. X-ray crystallography will find its future application in drug discovery by the development of specific tools that would allow realistic interpretation of the outcome coordinates and/or support testing of these hypotheses.
Collapse
Affiliation(s)
- Heping Zheng
- University of Virginia, Department of Molecular Physiology and Biological Physics, 1340 Jefferson Park Avenue, Charlottesville, VA 22908, USA
- Center for Structural Genomics of Infectious Diseases (CSGID)
- Midwest Center for Structural Genomics (MCSG), USA
- New York Structural Genomics Research Consortium (NYSGRC), USA
- Specializes in Protein Crystallography, Data Analytics and Data Mining, Research Scientist
| | - Jing Hou
- University of Virginia, Department of Molecular Physiology and Biological Physics, 1340 Jefferson Park Avenue, Charlottesville, VA 22908, USA
- Center for Structural Genomics of Infectious Diseases (CSGID)
- Enzyme Structure Initiative (EFI), USA
- New York Structural Genomics Research Consortium (NYSGRC), USA
- Specializes in Protein Crystallography, Research Associate
| | - Matthew D Zimmerman
- University of Virginia, Department of Molecular Physiology and Biological Physics, 1340 Jefferson Park Avenue, Charlottesville, VA 22908, USA
- Center for Structural Genomics of Infectious Diseases (CSGID)
- Enzyme Structure Initiative (EFI), USA
- Midwest Center for Structural Genomics (MCSG), USA
- New York Structural Genomics Research Consortium (NYSGRC), USA
- Specializes in Protein Crystallography, Data Mining and Management, Instructor of Research
| | - Alexander Wlodawer
- National Cancer Institute, Center for Cancer Research, Frederick, MD 21702, USA
- Specializes in Macromolecular Structure and Function, Chief of the Macromolecular Crystallography Laboratory
| | - Wladek Minor
- University of Virginia, Department of Molecular Physiology and Biological Physics, 1340 Jefferson Park Avenue, Charlottesville, VA 22908, USA
- Center for Structural Genomics of Infectious Diseases (CSGID)
- Enzyme Structure Initiative (EFI), USA
- Midwest Center for Structural Genomics (MCSG), USA
- New York Structural Genomics Research Consortium (NYSGRC), USA
- Specializes in Structural Biology, Data Mining and Management, Professor
| |
Collapse
|
13
|
Domagalski MJ, Zheng H, Zimmerman MD, Dauter Z, Wlodawer A, Minor W. The quality and validation of structures from structural genomics. Methods Mol Biol 2014; 1091:297-314. [PMID: 24203341 DOI: 10.1007/978-1-62703-691-7_21] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Quality control of three-dimensional structures of macromolecules is a critical step to ensure the integrity of structural biology data, especially those produced by structural genomics centers. Whereas the Protein Data Bank (PDB) has proven to be a remarkable success overall, the inconsistent quality of structures reveals a lack of universal standards for structure/deposit validation. Here, we review the state-of-the-art methods used in macromolecular structure validation, focusing on validation of structures determined by X-ray crystallography. We describe some general protocols used in the rebuilding and re-refinement of problematic structural models. We also briefly discuss some frontier areas of structure validation, including refinement of protein-ligand complexes, automation of structure redetermination, and the use of NMR structures and computational models to solve X-ray crystal structures by molecular replacement.
Collapse
Affiliation(s)
- Marcin J Domagalski
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA
| | | | | | | | | | | |
Collapse
|
14
|
Gioiello A, Venturoni F, Tamimi S, Custodi C, Pellicciari R, Macchiarulo A. Conformational properties of cholic acid, a lead compound at the crossroads of bile acid inspired drug discovery. MEDCHEMCOMM 2014. [DOI: 10.1039/c4md00024b] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
DFT and NMR spectroscopy studies unveil three major minima conformations of cholic acid that may affect its biological properties.
Collapse
Affiliation(s)
- Antimo Gioiello
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
| | - Francesco Venturoni
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
| | - Sara Tamimi
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
| | - Chiara Custodi
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
| | - Roberto Pellicciari
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
- TES Pharma S.r.l
- Perugia, Italy
| | - Antonio Macchiarulo
- Dipartimento di Chimica e Tecnologia del Farmaco
- Università degli Studi di Perugia
- 06123 Perugia, Italy
| |
Collapse
|
15
|
Salani B, Marini C, Rio AD, Ravera S, Massollo M, Orengo AM, Amaro A, Passalacqua M, Maffioli S, Pfeffer U, Cordera R, Maggi D, Sambuceti G. Metformin impairs glucose consumption and survival in Calu-1 cells by direct inhibition of hexokinase-II. Sci Rep 2013; 3:2070. [PMID: 23797762 PMCID: PMC3691576 DOI: 10.1038/srep02070] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2013] [Accepted: 06/03/2013] [Indexed: 12/27/2022] Open
Abstract
The anti-hyperglycaemic drug metformin has important anticancer properties as shown by the direct inhibition of cancer cells proliferation. Tumor cells avidly use glucose as a source for energy production and cell building blocks. Critical to this phenotype is the production of glucose-6-phosphate (G6P), catalysed by hexokinases (HK) I and II, whose role in glucose retention and metabolism is highly advantageous for cell survival and proliferation. Here we show that metformin impairs the enzymatic function of HKI and II in Calu-1 cells. This inhibition virtually abolishes cell glucose uptake and phosphorylation as documented by the reduced entrapment of 18F-fluorodeoxyglucose. In-silico models indicate that this action is due to metformin capability to mimic G6P features by steadily binding its pocket in HKII. The impairment of this energy source results in mitochondrial depolarization and subsequent cell death. These results could represent a starting point to open effective strategies in cancer prevention and treatment.
Collapse
Affiliation(s)
- Barbara Salani
- Department of Internal Medicine-DIMI, University of Genova, 16132 Genova, Italy
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
16
|
Yuriev E, Ramsland PA. Latest developments in molecular docking: 2010-2011 in review. J Mol Recognit 2013; 26:215-39. [PMID: 23526775 DOI: 10.1002/jmr.2266] [Citation(s) in RCA: 193] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2012] [Revised: 01/16/2013] [Accepted: 01/19/2013] [Indexed: 12/28/2022]
Affiliation(s)
- Elizabeth Yuriev
- Medicinal Chemistry, Monash Institute of Pharmaceutical Sciences; Monash University; Parkville; VIC; 3052; Australia
| | | |
Collapse
|
17
|
Davis BJ, Erlanson DA. Learning from our mistakes: the 'unknown knowns' in fragment screening. Bioorg Med Chem Lett 2013; 23:2844-52. [PMID: 23562240 DOI: 10.1016/j.bmcl.2013.03.028] [Citation(s) in RCA: 97] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2013] [Revised: 03/01/2013] [Accepted: 03/08/2013] [Indexed: 12/27/2022]
Abstract
In the past 15 years, fragment-based lead discovery (FBLD) has been adopted widely throughout academia and industry. The approach entails discovering very small molecular fragments and growing, merging, or linking them to produce drug leads. Because the affinities of the initial fragments are often low, detection methods are pushed to their limits, leading to a variety of artifacts, false positives, and false negatives that too often go unrecognized. This Digest discusses some of these problems and offers suggestions to avoid them. Although the primary focus is on FBLD, many of the lessons also apply to more established approaches such as high-throughput screening.
Collapse
Affiliation(s)
- Ben J Davis
- Vernalis Ltd., Granta Park, Great Abington, Cambridge CB21 6GB, UK.
| | | |
Collapse
|
18
|
Jones PM, George AM. Mechanism of the ABC transporter ATPase domains: catalytic models and the biochemical and biophysical record. Crit Rev Biochem Mol Biol 2012; 48:39-50. [PMID: 23131203 DOI: 10.3109/10409238.2012.735644] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
ABC transporters comprise a large, diverse, and ubiquitous superfamily of membrane active transporters. Their core architecture is a dimer of dimers, comprising two transmembrane domains that bind substrate and form the channel, and two ATP-binding cassettes, which bind and hydrolyze ATP to energize the translocase function. The prevailing paradigm for the ABC transport mechanism is the Switch Model, in which the nucleotide binding domains are proposed to dimerise upon binding of two ATP molecules, and thence dissociate upon sequential hydrolysis of the ATP. This idea appears consistent with crystal structures of both isolated subunits and whole transporters, as well as with a significant body of biochemical data. Nonetheless, an alternative Constant Contact Model has been proposed, in which the nucleotide binding domains do not fully dissociate, and ATP hydrolysis occurs alternately at each of the two active sites. Here, we review the biochemical and biophysical data relating to the ABC catalytic mechanism, to show how they may be construed as consistent with a Constant Contact Model, and to assess to what extent they support the Switch Model.
Collapse
Affiliation(s)
- Peter M Jones
- School of Medical and Molecular Biosciences, University of Technology Sydney, Broadway, NSW, Australia
| | | |
Collapse
|
19
|
Abstract
Fragment-based drug discovery (FBDD) concerns the screening of low-molecular weight compounds against macromolecular targets of clinical relevance. These compounds act as starting points for the development of drugs. FBDD has evolved and grown in popularity over the past 15 years. In this paper, the rationale and technology behind the use of X-ray crystallography in fragment based screening (FBS) will be described, including fragment library design and use of synchrotron radiation and robotics for high-throughput X-ray data collection. Some recent uses of crystallography in FBS will be described in detail, including interrogation of the drug targets β-secretase, phenylethanolamine N-methyltransferase, phosphodiesterase 4A and Hsp90. These examples provide illustrations of projects where crystallography is straightforward or difficult, and where other screening methods can help overcome the limitations of crystallography necessitated by diffraction quality.
Collapse
Affiliation(s)
- Zorik Chilingaryan
- School of Chemistry, University of Wollongong, Northfields Ave, Wollongong 2522, NSW, Australia.
| | | | | |
Collapse
|
20
|
Ullmann RT, Andrade SLA, Ullmann GM. Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations. J Phys Chem B 2012; 116:9690-703. [PMID: 22804733 DOI: 10.1021/jp305440f] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Amt-1 from Archaeoglobus fulgidus (AfAmt-1) belongs to the Amt/Rh family of ammonium/ammonia transporting membrane proteins. The transport mode and the precise microscopic permeation mechanism utilized by these proteins are intensely debated. Open questions concern the identity of the transported substrate (ammonia and/or ammonium) and whether the transport is passive or active. To address these questions, we studied the overall thermodynamics of the different transport modes as a function of the environmental conditions. Then, we investigated the thermodynamics of the underlying microscopic transport mechanisms with free energy calculations within a continuum electrostatics model. The formalism developed for this purpose is of general utility in the calculation of binding free energies for ligands with multiple protonation forms or other binding forms. The results of our calculations are compared to the available experimental and theoretical data on Amt/Rh proteins and discussed in light of the current knowledge on the physiological conditions experienced by microorganisms and plants. We found that microscopic models of electroneutral and electrogenic transport modes are in principle thermodynamically viable. However, only the electrogenic variants have a net thermodynamic driving force under the physiological conditions experienced by microorganisms and plants. Thus, the transport mechanism of AfAmt-1 is most likely electrogenic.
Collapse
Affiliation(s)
- R Thomas Ullmann
- Structural Biology/Bioinformatics, University of Bayreuth, Universitätsstrasse 30, BGI, 95447 Bayreuth, Germany.
| | | | | |
Collapse
|
21
|
Gleeson D, Tehan B, Gleeson MP, Limtrakul J. Evaluating the enthalpic contribution to ligand binding using QM calculations: effect of methodology on geometries and interaction energies. Org Biomol Chem 2012; 10:7053-61. [PMID: 22858758 DOI: 10.1039/c2ob25657f] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
As a result of research on ligand efficiency in the pharmaceutical industry, there is greater focus on optimizing the strength of polar interactions within receptors, so that the contribution of overall size and lipophilicity to binding can be decreased. A number of quantum mechanical (QM) methods involving simple probes are available to assess the H-bonding potential of different heterocycles or functional groups. However, in most receptors, multiple features are present, and these have distinct directionality, meaning very minimalist models may not be so ideal to describe the interactions. We describe how the use of gas phase QM models of kinase protein-ligand complex, which can more closely mimic the polar features of the active site region, can prove useful in assessing alterations to a core template, or different substituents. We investigate some practical issues surrounding the use of QM cluster models in structure based design (SBD). These include the choice of the method; semi-empirical, density functional theory or ab-initio, the choice of the basis set, whether to include implicit or explicit solvation, whether BSSE should be included, etc. We find a combination of the M06-2X method and the 6-31G* basis set is sufficiently rapid, and accurate, for the computation of structural and energetic parameters for this system.
Collapse
Affiliation(s)
- Duangkamol Gleeson
- Department of Chemistry, Faculty of Science, King Mongkut's Institute of Technology Ladkrabang, Bangkok 10520, Thailand
| | | | | | | |
Collapse
|
22
|
Nair PC, Malde AK, Drinkwater N, Mark AE. Missing fragments: detecting cooperative binding in fragment-based drug design. ACS Med Chem Lett 2012; 3:322-6. [PMID: 24900472 DOI: 10.1021/ml300015u] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2012] [Accepted: 02/14/2012] [Indexed: 11/28/2022] Open
Abstract
The aim of fragment-based drug design (FBDD) is to identify molecular fragments that bind to alternate subsites within a given binding pocket leading to cooperative binding when linked. In this study, the binding of fragments to human phenylethanolamine N-methyltransferase is used to illustrate how (a) current protocols may fail to detect fragments that bind cooperatively, (b) theoretical approaches can be used to validate potential hits, and (c) apparent false positives obtained when screening against cocktails of fragments may in fact indicate promising leads.
Collapse
Affiliation(s)
| | | | - Nyssa Drinkwater
- Randall Division
of Cell and Molecular Biophysics, King's College London, London, SE1 1UL, United Kingdom
| | | |
Collapse
|
23
|
Sitzmann M, Weidlich IE, Filippov IV, Liao C, Peach ML, Ihlenfeldt WD, Karki RG, Borodina YV, Cachau RE, Nicklaus MC. PDB ligand conformational energies calculated quantum-mechanically. J Chem Inf Model 2012; 52:739-56. [PMID: 22303903 PMCID: PMC7491702 DOI: 10.1021/ci200595n] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present here a greatly updated version of an earlier study on the conformational energies of protein-ligand complexes in the Protein Data Bank (PDB) [Nicklaus et al. Bioorg. Med. Chem. 1995, 3, 411-428], with the goal of improving on all possible aspects such as number and selection of ligand instances, energy calculations performed, and additional analyses conducted. Starting from about 357,000 ligand instances deposited in the 2008 version of the Ligand Expo database of the experimental 3D coordinates of all small-molecule instances in the PDB, we created a "high-quality" subset of ligand instances by various filtering steps including application of crystallographic quality criteria and structural unambiguousness. Submission of 640 Gaussian 03 jobs yielded a set of about 415 successfully concluded runs. We used a stepwise optimization of internal degrees of freedom at the DFT level of theory with the B3LYP/6-31G(d) basis set and a single-point energy calculation at B3LYP/6-311++G(3df,2p) after each round of (partial) optimization to separate energy changes due to bond length stretches vs bond angle changes vs torsion changes. Even for the most "conservative" choice of all the possible conformational energies-the energy difference between the conformation in which all internal degrees of freedom except torsions have been optimized and the fully optimized conformer-significant energy values were found. The range of 0 to ~25 kcal/mol was populated quite evenly and independently of the crystallographic resolution. A smaller number of "outliers" of yet higher energies were seen only at resolutions above 1.3 Å. The energies showed some correlation with molecular size and flexibility but not with crystallographic quality metrics such as the Cruickshank diffraction-component precision index (DPI) and R(free)-R, or with the ligand instance-specific metrics such as occupancy-weighted B-factor (OWAB), real-space R factor (RSR), and real-space correlation coefficient (RSCC). We repeated these calculations with the solvent model IEFPCM, which yielded energy differences that were generally somewhat lower than the corresponding vacuum results but did not produce a qualitatively different picture. Torsional sampling around the crystal conformation at the molecular mechanics level using the MMFF94s force field typically led to an increase in energy.
Collapse
Affiliation(s)
- Markus Sitzmann
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, National Institutes of Health, DHHS
| | - Iwona E. Weidlich
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, National Institutes of Health, DHHS
| | - Igor V. Filippov
- Basic Science Program, SAIC-Frederick, Inc., NCI-Frederick, 376 Boyles Street, Frederick, Maryland 21702, United States
| | - Chenzhong Liao
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, National Institutes of Health, DHHS
| | - Megan L. Peach
- Basic Science Program, SAIC-Frederick, Inc., NCI-Frederick, 376 Boyles Street, Frederick, Maryland 21702, United States
| | | | - Rajeshri G. Karki
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, National Institutes of Health, DHHS
| | - Yulia V. Borodina
- National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health, 8600 Rockville Pike, Bethesda Maryland 20894, United States
| | - Raul E. Cachau
- Advanced Structure Analysis Collaboratory, Information Systems Program, SAIC-Frederick, Inc., National Cancer Institute at Frederick, Frederick, Maryland 21702, United States
| | - Marc C. Nicklaus
- Chemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, National Institutes of Health, DHHS
| |
Collapse
|
24
|
Malde AK, Zuo L, Breeze M, Stroet M, Poger D, Nair PC, Oostenbrink C, Mark AE. An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. J Chem Theory Comput 2011; 7:4026-37. [DOI: 10.1021/ct200196m] [Citation(s) in RCA: 1122] [Impact Index Per Article: 86.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Alpeshkumar K. Malde
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - Le Zuo
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - Matthew Breeze
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - Martin Stroet
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - David Poger
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - Pramod C. Nair
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
| | - Chris Oostenbrink
- Leiden/Amsterdam Center for Drug Research, Division of Molecular Toxicology, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Alan E. Mark
- School of Chemistry and Molecular Biosciences, University of Queensland, St. Lucia, Australia
- Institute for Molecular Bioscience, University of Queensland, St. Lucia, QLD 4072, Australia
| |
Collapse
|
25
|
Cooper DR, Porebski PJ, Chruszcz M, Minor W. X-ray crystallography: Assessment and validation of protein-small molecule complexes for drug discovery. Expert Opin Drug Discov 2011; 6:771-782. [PMID: 21779303 PMCID: PMC3138648 DOI: 10.1517/17460441.2011.585154] [Citation(s) in RCA: 49] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
INTRODUCTION: Crystallography is the key initial component for structure-based and fragment-based drug design and can often generate leads that can be developed into high potency drugs. Therefore, huge sums of money are committed based on the outcome of crystallography experiments and their interpretation. AREAS COVERED: This review discusses how to evaluate the correctness of an X-ray structure, focusing on the validation of small molecule-protein complexes. Various types of inaccuracies found within the PDB are identified and the ramifications of these errors are discussed. The reader will gain an understanding of the key parameters that need to be inspected before a structure can be used in drug discovery efforts, as well as an appreciation of the difficulties of correctly interpreting electron density for small molecules. The reader will also be introduced to methods for validating small molecules within the context of a macromolecular structure. EXPERT OPINION: One of the reasons that ligand identification and positioning, within a macromolecular crystal structure, is so difficult is that the quality of small molecules widely varies in the PDB. For this reason, the PDB can not always be considered a reliable repository of structural information pertaining to small molecules, and this makes the derivation of general principles that govern small molecule-protein interactions more difficult.
Collapse
Affiliation(s)
- David R Cooper
- Department of Molecular Physiology and Biological Physics, University of Virginia, 1340 Jefferson Park Avenue, Charlottesville, VA 22908, USA
| | | | | | | |
Collapse
|
26
|
Nair PC, Malde AK, Mark AE. Using Theory to Reconcile Experiment: The Structural and Thermodynamic Basis of Ligand Recognition by Phenylethanolamine N-Methyltransferase (PNMT). J Chem Theory Comput 2011; 7:1458-68. [DOI: 10.1021/ct1007229] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Pramod C. Nair
- School of Chemistry and Molecular Biosciences (SCMB) and ‡Institute for Molecular Bioscience (IMB), The University of Queensland (UQ), St. Lucia Campus, Brisbane, QLD 4072 Australia
| | - Alpeshkumar K. Malde
- School of Chemistry and Molecular Biosciences (SCMB) and ‡Institute for Molecular Bioscience (IMB), The University of Queensland (UQ), St. Lucia Campus, Brisbane, QLD 4072 Australia
| | - Alan E. Mark
- School of Chemistry and Molecular Biosciences (SCMB) and ‡Institute for Molecular Bioscience (IMB), The University of Queensland (UQ), St. Lucia Campus, Brisbane, QLD 4072 Australia
| |
Collapse
|