51
|
Ryazantsev MN, Nikolaev DM, Struts AV, Brown MF. Quantum Mechanical and Molecular Mechanics Modeling of Membrane-Embedded Rhodopsins. J Membr Biol 2019; 252:425-449. [PMID: 31570961 DOI: 10.1007/s00232-019-00095-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Accepted: 09/10/2019] [Indexed: 12/20/2022]
Abstract
Computational chemistry provides versatile methods for studying the properties and functioning of biological systems at different levels of precision and at different time scales. The aim of this article is to review the computational methodologies that are applicable to rhodopsins as archetypes for photoactive membrane proteins that are of great importance both in nature and in modern technologies. For each class of computational techniques, from methods that use quantum mechanics for simulating rhodopsin photophysics to less-accurate coarse-grained methodologies used for long-scale protein dynamics, we consider possible applications and the main directions for improvement.
Collapse
Affiliation(s)
- Mikhail N Ryazantsev
- Institute of Chemistry, Saint Petersburg State University, 26 Universitetskii pr, Saint Petersburg, Russia, 198504
| | - Dmitrii M Nikolaev
- Saint-Petersburg Academic University - Nanotechnology Research and Education Centre RAS, Saint Petersburg, Russia, 194021
| | - Andrey V Struts
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ, 85721, USA.,Laboratory of Biomolecular NMR, Saint Petersburg State University, Saint Petersburg, Russia, 199034
| | - Michael F Brown
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, AZ, 85721, USA. .,Department of Physics, University of Arizona, Tucson, AZ, 85721, USA.
| |
Collapse
|
52
|
Pal RK, Gadhiya S, Ramsey S, Cordone P, Wickstrom L, Harding WW, Kurtzman T, Gallicchio E. Inclusion of enclosed hydration effects in the binding free energy estimation of dopamine D3 receptor complexes. PLoS One 2019; 14:e0222902. [PMID: 31568493 PMCID: PMC6768453 DOI: 10.1371/journal.pone.0222902] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Accepted: 08/30/2019] [Indexed: 01/04/2023] Open
Abstract
Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.
Collapse
Affiliation(s)
- Rajat Kumar Pal
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, NY 11210, United States of America
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
| | - Satishkumar Gadhiya
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Steven Ramsey
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Lehman College, 250 Bedford Park Blvd. West, Bronx, NY 10468, United States of America
| | - Pierpaolo Cordone
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Lauren Wickstrom
- Department of Science, Borough of Manhattan Community College, 199 Chambers Street, New York, NY 10007, United States of America
| | - Wayne W. Harding
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Tom Kurtzman
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Lehman College, 250 Bedford Park Blvd. West, Bronx, NY 10468, United States of America
| | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, NY 11210, United States of America
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- * E-mail:
| |
Collapse
|
53
|
Jones L, Tynes M, Smith P. Prediction of models for ordered solvent in macromolecular structures by a classifier based upon resolution-independent projections of local feature data. Acta Crystallogr D Struct Biol 2019; 75:696-717. [PMID: 31373570 PMCID: PMC6677017 DOI: 10.1107/s2059798319008933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Accepted: 06/22/2019] [Indexed: 11/11/2022] Open
Abstract
Current software tools for the automated building of models for macromolecular X-ray crystal structures are capable of assembling high-quality models for ordered macromolecule and small-molecule scattering components with minimal or no user supervision. Many of these tools also incorporate robust functionality for modelling the ordered water molecules that are found in nearly all macromolecular crystal structures. However, no current tools focus on differentiating these ubiquitous water molecules from other frequently occurring multi-atom solvent species, such as sulfate, or the automated building of models for such species. PeakProbe has been developed specifically to address the need for such a tool. PeakProbe predicts likely solvent models for a given point (termed a `peak') in a structure based on analysis (`probing') of its local electron density and chemical environment. PeakProbe maps a total of 19 resolution-dependent features associated with electron density and two associated with the local chemical environment to a two-dimensional score space that is independent of resolution. Peaks are classified based on the relative frequencies with which four different classes of solvent (including water) are observed within a given region of this score space as determined by large-scale sampling of solvent models in the Protein Data Bank. Designed to classify peaks generated from difference density maxima, PeakProbe also incorporates functionality for identifying peaks associated with model errors or clusters of peaks likely to correspond to multi-atom solvent, and for the validation of existing solvent models using solvent-omit electron-density maps. When tasked with classifying peaks into one of four distinct solvent classes, PeakProbe achieves greater than 99% accuracy for both peaks derived directly from the atomic coordinates of existing solvent models and those based on difference density maxima. While the program is still under development, a fully functional version is publicly available. PeakProbe makes extensive use of cctbx libraries, and requires a PHENIX licence and an up-to-date phenix.python environment for execution.
Collapse
Affiliation(s)
- Laurel Jones
- Department of Chemistry, Fordham University, Bronx, NY 10458, USA
| | - Michael Tynes
- Department of Computer and Information Science, Fordham University, Bronx, NY 10458, USA
| | - Paul Smith
- Department of Chemistry, Fordham University, Bronx, NY 10458, USA
| |
Collapse
|
54
|
Labbé CM, Pencheva T, Jereva D, Desvillechabrol D, Becot J, Villoutreix BO, Pajeva I, Miteva MA. AMMOS2: a web server for protein-ligand-water complexes refinement via molecular mechanics. Nucleic Acids Res 2019; 45:W350-W355. [PMID: 28486703 PMCID: PMC5570140 DOI: 10.1093/nar/gkx397] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Accepted: 04/28/2017] [Indexed: 12/20/2022] Open
Abstract
AMMOS2 is an interactive web server for efficient computational refinement of protein-small organic molecule complexes. The AMMOS2 protocol employs atomic-level energy minimization of a large number of experimental or modeled protein-ligand complexes. The web server is based on the previously developed standalone software AMMOS (Automatic Molecular Mechanics Optimization for in silico Screening). AMMOS utilizes the physics-based force field AMMP sp4 and performs optimization of protein-ligand interactions at five levels of flexibility of the protein receptor. The new version 2 of AMMOS implemented in the AMMOS2 web server allows the users to include explicit water molecules and individual metal ions in the protein-ligand complexes during minimization. The web server provides comprehensive analysis of computed energies and interactive visualization of refined protein-ligand complexes. The ligands are ranked by the minimized binding energies allowing the users to perform additional analysis for drug discovery or chemical biology projects. The web server has been extensively tested on 21 diverse protein-ligand complexes. AMMOS2 minimization shows consistent improvement over the initial complex structures in terms of minimized protein-ligand binding energies and water positions optimization. The AMMOS2 web server is freely available without any registration requirement at the URL: http://drugmod.rpbs.univ-paris-diderot.fr/ammosHome.php.
Collapse
Affiliation(s)
- Céline M Labbé
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico, INSERM UMR-S 973, Paris, France.,INSERM, U973 Paris, France
| | - Tania Pencheva
- Department of QSAR and Molecular Modelling, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, 105 Acad G. Bonchev Str., 1113 Sofia, Bulgaria
| | - Dessislava Jereva
- Department of QSAR and Molecular Modelling, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, 105 Acad G. Bonchev Str., 1113 Sofia, Bulgaria
| | - Dimitri Desvillechabrol
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico, INSERM UMR-S 973, Paris, France.,INSERM, U973 Paris, France
| | - Jérôme Becot
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico, INSERM UMR-S 973, Paris, France.,INSERM, U973 Paris, France
| | - Bruno O Villoutreix
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico, INSERM UMR-S 973, Paris, France.,INSERM, U973 Paris, France
| | - Ilza Pajeva
- Department of QSAR and Molecular Modelling, Institute of Biophysics and Biomedical Engineering, Bulgarian Academy of Sciences, 105 Acad G. Bonchev Str., 1113 Sofia, Bulgaria
| | - Maria A Miteva
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico, INSERM UMR-S 973, Paris, France.,INSERM, U973 Paris, France
| |
Collapse
|
55
|
Pradhan MR, Nguyen MN, Kannan S, Fox SJ, Kwoh CK, Lane DP, Verma CS. Characterization of Hydration Properties in Structural Ensembles of Biomolecules. J Chem Inf Model 2019; 59:3316-3329. [DOI: 10.1021/acs.jcim.8b00453] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Mohan R. Pradhan
- Bioinformatics Institute, A*STAR (Agency for Science, Technology and Research), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
- School of Computer Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798
| | - Minh N. Nguyen
- Bioinformatics Institute, A*STAR (Agency for Science, Technology and Research), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Srinivasaraghavan Kannan
- Bioinformatics Institute, A*STAR (Agency for Science, Technology and Research), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Stephen J. Fox
- Bioinformatics Institute, A*STAR (Agency for Science, Technology and Research), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
| | - Chee Keong Kwoh
- School of Computer Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798
| | - David P. Lane
- p53 Laboratory, A*STAR (Agency for Science, Technology and Research), 8A Biomedical Grove, #06-04/05, Neuros/Immunos, Singapore 138648
| | - Chandra S. Verma
- Bioinformatics Institute, A*STAR (Agency for Science, Technology and Research), 30 Biopolis Street, #07-01 Matrix, Singapore 138671
- School of Biological Sciences, Nanyang Technological University, 50 Nanyang Drive, Singapore 637551
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543
| |
Collapse
|
56
|
Wei W, Luo J, Waldispühl J, Moitessier N. Predicting Positions of Bridging Water Molecules in Nucleic Acid-Ligand Complexes. J Chem Inf Model 2019; 59:2941-2951. [PMID: 30998377 DOI: 10.1021/acs.jcim.9b00163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Over the past two decades, interests in DNA and RNA as drug targets have been growing rapidly. Following the trends observed with protein drug targets, computational approaches for drug design have been developed for this new class of molecules. Our efforts toward the development of a universal docking program, Fitted, led us to focus on nucleic acids. Throughout the development of this docking program, efforts were directed toward displaceable water molecules which must be accurately located for optimal docking-based drug discovery. However, although there is a plethora of methods to place water molecules in and around protein structures, there is, to the best of our knowledge, no such fully automated method for nucleic acids, which are significantly more polar and solvated than proteins. We report herein a new method, Splash'Em (Solvation Potential Laid around Statistical Hydration on Entire Macromolecules) developed to place water molecules within the binding cavity of nucleic acids. This fast method was shown to have high agreement with water positions in crystal structures and will therefore provide essential information to medicinal chemists.
Collapse
|
57
|
Slater O, Kontoyianni M. The compromise of virtual screening and its impact on drug discovery. Expert Opin Drug Discov 2019; 14:619-637. [PMID: 31025886 DOI: 10.1080/17460441.2019.1604677] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Introduction: Docking and structure-based virtual screening (VS) have been standard approaches in structure-based design for over two decades. However, our understanding of the limitations, potential, and strength of these techniques has enhanced, raising expectations. Areas covered: Based on a survey of reports in the past five years, we assess whether VS: (1) predicts binding poses in agreement with crystallographic data (when available); (2) is a superior screening tool, as often claimed; (3) is successful in identifying chemical scaffolds that can be starting points for subsequent lead optimization cycles. Data shows that knowledge of the target and its chemotypes in postprocessing lead to viable hits in early drug discovery endeavors. Expert opinion: VS is capable of accurate placements in the pocket for the most part, but does not consistently score screening collections accurately. What matters is capitalization on available resources to get closer to a viable lead or optimizable series. Integration of approaches, subjective hit selection guided by knowledge of the receptor or endogenous ligand, libraries driven by experimental guides, validation studies to identify the best docking/scoring that reproduces experimental findings, constraints regarding receptor-ligand interactions, thoroughly designed methodologies, and predefined cutoff scoring criteria strengthen VS's position in pharmaceutical research.
Collapse
Affiliation(s)
- Olivia Slater
- a Department of Pharmaceutical Sciences , Southern Illinois University Edwardsville , Edwardsville , IL , USA
| | - Maria Kontoyianni
- a Department of Pharmaceutical Sciences , Southern Illinois University Edwardsville , Edwardsville , IL , USA
| |
Collapse
|
58
|
Li J, Fu A, Zhang L. An Overview of Scoring Functions Used for Protein-Ligand Interactions in Molecular Docking. Interdiscip Sci 2019; 11:320-328. [PMID: 30877639 DOI: 10.1007/s12539-019-00327-w] [Citation(s) in RCA: 166] [Impact Index Per Article: 33.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2018] [Revised: 02/06/2019] [Accepted: 03/06/2019] [Indexed: 12/17/2022]
Abstract
Currently, molecular docking is becoming a key tool in drug discovery and molecular modeling applications. The reliability of molecular docking depends on the accuracy of the adopted scoring function, which can guide and determine the ligand poses when thousands of possible poses of ligand are generated. The scoring function can be used to determine the binding mode and site of a ligand, predict binding affinity and identify the potential drug leads for a given protein target. Despite intensive research over the years, accurate and rapid prediction of protein-ligand interactions is still a challenge in molecular docking. For this reason, this study reviews four basic types of scoring functions, physics-based, empirical, knowledge-based, and machine learning-based scoring functions, based on an up-to-date classification scheme. We not only discuss the foundations of the four types scoring functions, suitable application areas and shortcomings, but also discuss challenges and potential future study directions.
Collapse
Affiliation(s)
- Jin Li
- College of Computer and Information Science, Southwest University, Chongqing, 400715, China.,School of Medical Information and Engineering, Southwest Medical University, Luzhou, 646000, China
| | - Ailing Fu
- College of Pharmaceutical Sciences, Southwest University, Chongqing, 400715, China
| | - Le Zhang
- College of Computer and Information Science, Southwest University, Chongqing, 400715, China. .,College of Computer Science, Sichuan University, Chengdu, 610065, China. .,Medical Big Data Center, Sichuan University, Chengdu, 610065, China. .,Zdmedical, Information Polytron Technologies Inc Chongqing, Chongqing, 401320, China.
| |
Collapse
|
59
|
Molecular insights of newly identified potential peptide inhibitors of hypoxia inducible factor 1α causing breast cancer. J Mol Struct 2019. [DOI: 10.1016/j.molstruc.2018.09.072] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
60
|
Bruce Macdonald HE, Cave-Ayland C, Ross GA, Essex JW. Ligand Binding Free Energies with Adaptive Water Networks: Two-Dimensional Grand Canonical Alchemical Perturbations. J Chem Theory Comput 2018; 14:6586-6597. [PMID: 30451501 PMCID: PMC6293443 DOI: 10.1021/acs.jctc.8b00614] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
![]()
Computational methods
to calculate ligand binding affinities are
increasing in popularity, due to improvements in simulation algorithms,
computational resources, and easy-to-use software. However, issues
can arise in relative ligand binding free energy simulations if the
ligands considered have different active site water networks, as simulations
are typically performed with a predetermined number of water molecules
(fixed N ensembles) in preassigned locations. If an alchemical perturbation
is attempted where the change should result in a different active
site water network, the water molecules may not be able to adapt appropriately
within the time scales of the simulations—particularly if the
active site is occluded. By combining the grand canonical ensemble
(μVT) to sample active site water molecules, with conventional
alchemical free energy methods, the water network is able to dynamically
adapt to the changing ligand. We refer to this approach as grand canonical
alchemical perturbation (GCAP). In this work we demonstrate GCAP for
two systems; Scytalone Dehydratase (SD) and Adenosine A2A receptor. For both systems, GCAP is shown to perform
well at reproducing experimental binding affinities. Calculating the
relative binding affinities with a naïve, conventional
attempt to solvate the active site illustrates how poor results can
be if proper consideration of water molecules in occluded pockets
is neglected. GCAP results are shown to be consistent with time-consuming
double decoupling simulations. In addition, by obtaining the free
energy surface for ligand perturbations, as a function of both the
free energy coupling parameter and water chemical potential, it is
possible to directly deconvolute the binding energetics in terms of
protein–ligand direct interactions and protein binding site
hydration.
Collapse
Affiliation(s)
| | | | - Gregory A Ross
- Computational and Systems Biology Program, Sloan Kettering Institute , Memorial Sloan Kettering Cancer Center , New York, New York 10065 , United States
| | - Jonathan W Essex
- School of Chemistry , University of Southampton , Southampton , SO17 1BJ , U.K
| |
Collapse
|
61
|
Brysbaert G, Blossey R, Lensink MF. The Inclusion of Water Molecules in Residue Interaction Networks Identifies Additional Central Residues. Front Mol Biosci 2018; 5:88. [PMID: 30364190 PMCID: PMC6193073 DOI: 10.3389/fmolb.2018.00088] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2018] [Accepted: 09/18/2018] [Indexed: 12/13/2022] Open
Abstract
The relevance of water molecules for the recognition and the interaction of biomolecules is widely appreciated. In this paper we address the role that water molecules associated to protein complexes play for the functional relevance of residues by considering their residue interaction networks (RINs). These are commonly defined on the basis of the amino acid composition of the proteins themselves, disregarding the solvation state of the protein. We determine properties of the RINs of two protein complexes, colicin E2/Im2 and barnase/barstar, with and without associated water molecules, using a previously developed methodology and its associated application RINspector. We find that the inclusion of water molecules in RINs leads to an increase in the number of central residues which adds a novel mechanism to the relevance of water molecules for protein function.
Collapse
Affiliation(s)
- Guillaume Brysbaert
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| | - Ralf Blossey
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| | - Marc F Lensink
- CNRS UMR8576 UGSF, Institute for Structural and Functional Glycobiology, University of Lille, Lille, France
| |
Collapse
|
62
|
Farahvash A, Stuchebrukhov A. Investigating the Many Roles of Internal Water in Cytochrome c Oxidase. J Phys Chem B 2018; 122:7625-7635. [PMID: 30011995 DOI: 10.1021/acs.jpcb.7b11920] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.
Collapse
Affiliation(s)
- Ardavan Farahvash
- Department of Chemistry , University of California-Davis , One Shields Avenue , Davis , California 95616 , United States
| | - Alexei Stuchebrukhov
- Department of Chemistry , University of California-Davis , One Shields Avenue , Davis , California 95616 , United States
| |
Collapse
|
63
|
Flavonoids as Putative Epi-Modulators: Insight into Their Binding Mode with BRD4 Bromodomains Using Molecular Docking and Dynamics. Biomolecules 2018; 8:biom8030061. [PMID: 30041464 PMCID: PMC6164663 DOI: 10.3390/biom8030061] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2018] [Revised: 07/14/2018] [Accepted: 07/18/2018] [Indexed: 02/07/2023] Open
Abstract
Flavonoids are widely recognized as natural polydrugs, given their anti-inflammatory, antioxidant, sedative, and antineoplastic activities. Recently, different studies showed that flavonoids have the potential to inhibit bromodomain and extraterminal (BET) bromodomains. Previous reports suggested that flavonoids bind between the Z and A loops of the bromodomain (ZA channel) due to their orientation and interactions with P86, V87, L92, L94, and N140. Herein, a comprehensive characterization of the binding modes of fisetin and the biflavonoid, amentoflavone, is discussed. To this end, both compounds were docked with BET bromodomain 4 (BRD4) using four docking programs. The results were post-processed with protein–ligand interaction fingerprints. To gain further insight into the binding mode of the two natural products, the docking results were further analyzed with molecular dynamics simulations. The results showed that amentoflavone makes numerous contacts in the ZA channel, as previously described for flavonoids and kinase inhibitors. It was also found that amentoflavone can potentially make contacts with non-canonical residues for BET inhibition. Most of these contacts were not observed with fisetin. Based on these results, amentoflavone was experimentally tested for BRD4 inhibition, showing activity in the micromolar range. This work may serve as the basis for scaffold optimization and the further characterization of flavonoids as BET inhibitors.
Collapse
|
64
|
Nittinger E, Flachsenberg F, Bietz S, Lange G, Klein R, Rarey M. Placement of Water Molecules in Protein Structures: From Large-Scale Evaluations to Single-Case Examples. J Chem Inf Model 2018; 58:1625-1637. [DOI: 10.1021/acs.jcim.8b00271] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Eva Nittinger
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Florian Flachsenberg
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Stefan Bietz
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Gudrun Lange
- Bayer CropScience AG, Industriepark Hoechst G836, 65926 Frankfurt am Main, Germany
| | - Robert Klein
- Bayer CropScience AG, Industriepark Hoechst G836, 65926 Frankfurt am Main, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH − Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| |
Collapse
|
65
|
Takemura K, Matubayasi N, Kitao A. Binding free energy analysis of protein-protein docking model structures by evERdock. J Chem Phys 2018; 148:105101. [PMID: 29544320 DOI: 10.1063/1.5019864] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
To aid the evaluation of protein-protein complex model structures generated by protein docking prediction (decoys), we previously developed a method to calculate the binding free energies for complexes. The method combines a short (2 ns) all-atom molecular dynamics simulation with explicit solvent and solution theory in the energy representation (ER). We showed that this method successfully selected structures similar to the native complex structure (near-native decoys) as the lowest binding free energy structures. In our current work, we applied this method (evERdock) to 100 or 300 model structures of four protein-protein complexes. The crystal structures and the near-native decoys showed the lowest binding free energy of all the examined structures, indicating that evERdock can successfully evaluate decoys. Several decoys that show low interface root-mean-square distance but relatively high binding free energy were also identified. Analysis of the fraction of native contacts, hydrogen bonds, and salt bridges at the protein-protein interface indicated that these decoys were insufficiently optimized at the interface. After optimizing the interactions around the interface by including interfacial water molecules, the binding free energies of these decoys were improved. We also investigated the effect of solute entropy on binding free energy and found that consideration of the entropy term does not necessarily improve the evaluations of decoys using the normal model analysis for entropy calculation.
Collapse
Affiliation(s)
- Kazuhiro Takemura
- Institute of Molecular and Cellular Biosciences, University of Tokyo, Bunkyo, Tokyo 113-0032, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Akio Kitao
- Institute of Molecular and Cellular Biosciences, University of Tokyo, Bunkyo, Tokyo 113-0032, Japan
| |
Collapse
|
66
|
Lee AG. A Database of Predicted Binding Sites for Cholesterol on Membrane Proteins, Deep in the Membrane. Biophys J 2018; 115:522-532. [PMID: 30007584 DOI: 10.1016/j.bpj.2018.06.022] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2018] [Revised: 06/18/2018] [Accepted: 06/19/2018] [Indexed: 01/13/2023] Open
Abstract
The outer membranes of animal cells contain high concentrations of cholesterol, of which a small proportion is located deep within the hydrophobic core of the membrane. An automated docking procedure is described that allows the characterization of binding sites for these deep cholesterol molecules on the membrane-spanning surfaces of membrane proteins and in protein cavities or pores, driven by hydrogen bond formation. A database of this class of predicted binding site is described, covering 397 high-resolution structures. The database includes sites on the transmembrane surfaces of many G-protein coupled receptors; within the fenestrations of two-pore K+ channels and ATP-gated P2X3 channels; in the central cavities of a number of transporters, including Glut1, Glut5, and P-glycoprotein; and in deep clefts in mitochondrial complexes III and IV.
Collapse
Affiliation(s)
- Anthony G Lee
- Centre for Biological Sciences, University of Southampton, Southampton, United Kingdom.
| |
Collapse
|
67
|
Jeszenői N, Schilli G, Bálint M, Horváth I, Hetényi C. Analysis of the influence of simulation parameters on biomolecule-linked water networks. J Mol Graph Model 2018; 82:117-128. [DOI: 10.1016/j.jmgm.2018.04.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2017] [Revised: 03/27/2018] [Accepted: 04/21/2018] [Indexed: 12/11/2022]
|
68
|
Maffucci I, Hu X, Fumagalli V, Contini A. An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings. Front Chem 2018; 6:43. [PMID: 29556494 PMCID: PMC5844977 DOI: 10.3389/fchem.2018.00043] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Accepted: 02/19/2018] [Indexed: 01/05/2023] Open
Abstract
Nwat-MMGBSA is a variant of MM-PB/GBSA based on the inclusion of a number of explicit water molecules that are the closest to the ligand in each frame of a molecular dynamics trajectory. This method demonstrated improved correlations between calculated and experimental binding energies in both protein-protein interactions and ligand-receptor complexes, in comparison to the standard MM-GBSA. A protocol optimization, aimed to maximize efficacy and efficiency, is discussed here considering penicillopepsin, HIV1-protease, and BCL-XL as test cases. Calculations were performed in triplicates on both classic HPC environments and on standard workstations equipped by a GPU card, evidencing no statistical differences in the results. No relevant differences in correlation to experiments were also observed when performing Nwat-MMGBSA calculations on 4 or 1 ns long trajectories. A fully automatic workflow for structure-based virtual screening, performing from library set-up to docking and Nwat-MMGBSA rescoring, has then been developed. The protocol has been tested against no rescoring or standard MM-GBSA rescoring within a retrospective virtual screening of inhibitors of AmpC β-lactamase and of the Rac1-Tiam1 protein-protein interaction. In both cases, Nwat-MMGBSA rescoring provided a statistically significant increase in the ROC AUCs of between 20 and 30%, compared to docking scoring or to standard MM-GBSA rescoring.
Collapse
Affiliation(s)
- Irene Maffucci
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica "Alessandro Marchesini," Università degli Studi di Milano, Milan, Italy
| | - Xiao Hu
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica "Alessandro Marchesini," Università degli Studi di Milano, Milan, Italy
| | - Valentina Fumagalli
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica "Alessandro Marchesini," Università degli Studi di Milano, Milan, Italy
| | - Alessandro Contini
- Dipartimento di Scienze Farmaceutiche, Sezione di Chimica Generale e Organica "Alessandro Marchesini," Università degli Studi di Milano, Milan, Italy
| |
Collapse
|
69
|
Graham SE, Smith RD, Carlson HA. Predicting Displaceable Water Sites Using Mixed-Solvent Molecular Dynamics. J Chem Inf Model 2018; 58:305-314. [PMID: 29286658 PMCID: PMC6190669 DOI: 10.1021/acs.jcim.7b00268] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Water molecules are an important factor in protein-ligand binding. Upon binding of a ligand with a protein's surface, waters can either be displaced by the ligand or may be conserved and possibly bridge interactions between the protein and ligand. Depending on the specific interactions made by the ligand, displacing waters can yield a gain in binding affinity. The extent to which binding affinity may increase is difficult to predict, as the favorable displacement of a water molecule is dependent on the site-specific interactions made by the water and the potential ligand. Several methods have been developed to predict the location of water sites on a protein's surface, but the majority of methods are not able to take into account both protein dynamics and the interactions made by specific functional groups. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that explicitly accounts for the interaction of both water and small molecule probes with a protein's surface, allowing for their direct competition. This method has previously been shown to identify both active and allosteric sites on a protein's surface. Using a test set of eight systems, we have developed a method using MixMD to identify conserved and displaceable water sites. Conserved sites can be determined by an occupancy-based metric to identify sites which are consistently occupied by water even in the presence of probe molecules. Conversely, displaceable water sites can be found by considering the sites which preferentially bind probe molecules. Furthermore, the inclusion of six probe types allows the MixMD method to predict which functional groups are capable of displacing which water sites. The MixMD method consistently identifies sites which are likely to be nondisplaceable and predicts the favorable displacement of water sites that are known to be displaced upon ligand binding.
Collapse
Affiliation(s)
- Sarah E. Graham
- Department of Biophysics, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, Michigan, 48109-1065
| | - Richard D. Smith
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, Michigan, 48109-1065
| | - Heather A. Carlson
- Department of Biophysics, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, Michigan, 48109-1065
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan, 428 Church St., Ann Arbor, Michigan, 48109-1065
| |
Collapse
|
70
|
Mota F, Fotinou C, Rana RR, Chan AWE, Yelland T, Arooz MT, O'Leary AP, Hutton J, Frankel P, Zachary I, Selwood D, Djordjevic S. Architecture and hydration of the arginine-binding site of neuropilin-1. FEBS J 2018; 285:1290-1304. [PMID: 29430837 PMCID: PMC5947257 DOI: 10.1111/febs.14405] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2017] [Revised: 01/26/2018] [Accepted: 02/06/2018] [Indexed: 12/15/2022]
Abstract
Neuropilin‐1 (NRP1) is a transmembrane co‐receptor involved in binding interactions with variety of ligands and receptors, including receptor tyrosine kinases. Expression of NRP1 in several cancers correlates with cancer stages and poor prognosis. Thus, NRP1 has been considered a therapeutic target and is the focus of multiple drug discovery initiatives. Vascular endothelial growth factor (VEGF) binds to the b1 domain of NRP1 through interactions between the C‐terminal arginine of VEGF and residues in the NRP1‐binding site including Tyr297, Tyr353, Asp320, Ser346 and Thr349. We obtained several complexes of the synthetic ligands and the NRP1‐b1 domain and used X‐ray crystallography and computational methods to analyse atomic details and hydration profile of this binding site. We observed side chain flexibility for Tyr297 and Asp320 in the six new high‐resolution crystal structures of arginine analogues bound to NRP1. In addition, we identified conserved water molecules in binding site regions which can be targeted for drug design. The computational prediction of the VEGF ligand‐binding site hydration map of NRP1 was in agreement with the experimentally derived, conserved hydration structure. Displacement of certain conserved water molecules by a ligand's functional groups may contribute to binding affinity, whilst other water molecules perform as protein–ligand bridges. Our report provides a comprehensive description of the binding site for the peptidic ligands’ C‐terminal arginines in the b1 domain of NRP1, highlights the importance of conserved structural waters in drug design and validates the utility of the computational hydration map prediction method in the context of neuropilin. Database The structures were deposited to the PDB with accession numbers PDB ID: 5IJR, 5IYY, 5JHK, 5J1X, 5JGQ, 5JGI.
Collapse
Affiliation(s)
- Filipa Mota
- Magnus Life, Magnus Life Science, London, UK
| | | | | | - A W Edith Chan
- Wolfson Institute for Biomedical Research, University College London, UK
| | | | - Mohamed T Arooz
- The Institute of Structural and Molecular Biology, University College London, UK
| | | | | | - Paul Frankel
- Magnus Life, Magnus Life Science, London, UK.,Centre for Cardiovascular Biology & Medicine, BHF Laboratories at University College London, UK
| | - Ian Zachary
- Centre for Cardiovascular Biology & Medicine, BHF Laboratories at University College London, UK
| | - David Selwood
- Wolfson Institute for Biomedical Research, University College London, UK
| | - Snezana Djordjevic
- The Institute of Structural and Molecular Biology, University College London, UK
| |
Collapse
|
71
|
Rudling A, Orro A, Carlsson J. Prediction of Ordered Water Molecules in Protein Binding Sites from Molecular Dynamics Simulations: The Impact of Ligand Binding on Hydration Networks. J Chem Inf Model 2018; 58:350-361. [PMID: 29308882 PMCID: PMC6716772 DOI: 10.1021/acs.jcim.7b00520] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Water
plays a major role in ligand binding and is attracting increasing
attention in structure-based drug design. Water molecules can make
large contributions to binding affinity by bridging protein–ligand
interactions or by being displaced upon complex formation, but these
phenomena are challenging to model at the molecular level. Herein,
networks of ordered water molecules in protein binding sites were
analyzed by clustering of molecular dynamics (MD) simulation trajectories.
Locations of ordered waters (hydration sites) were first identified
from simulations of high resolution crystal structures of 13 protein–ligand
complexes. The MD-derived hydration sites reproduced 73% of the binding
site water molecules observed in the crystal structures. If the simulations
were repeated without the cocrystallized ligands, a majority (58%)
of the crystal waters in the binding sites were still predicted. In
addition, comparison of the hydration sites obtained from simulations
carried out in the absence of ligands to those identified for the
complexes revealed that the networks of ordered water molecules were
preserved to a large extent, suggesting that the locations of waters
in a protein–ligand interface are mainly dictated by the protein.
Analysis of >1000 crystal structures showed that hydration sites
bridged
protein–ligand interactions in complexes with different ligands,
and those with high MD-derived occupancies were more likely to correspond
to experimentally observed ordered water molecules. The results demonstrate
that ordered water molecules relevant for modeling of protein–ligand
complexes can be identified from MD simulations. Our findings could
contribute to development of improved methods for structure-based
virtual screening and lead optimization.
Collapse
Affiliation(s)
- Axel Rudling
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Adolfo Orro
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC , Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
72
|
Cuzzolin A, Deganutti G, Salmaso V, Sturlese M, Moro S. AquaMMapS: An Alternative Tool to Monitor the Role of Water Molecules During Protein-Ligand Association. ChemMedChem 2018; 13:522-531. [DOI: 10.1002/cmdc.201700564] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2017] [Revised: 11/21/2017] [Indexed: 12/11/2022]
Affiliation(s)
- Alberto Cuzzolin
- Molecular Modeling Section, MMS, Department of Pharmaceutical and Pharmacological Sciences; University of Padova; via Marzolo 5 35131 Padova Italy
| | - Giuseppe Deganutti
- Molecular Modeling Section, MMS, Department of Pharmaceutical and Pharmacological Sciences; University of Padova; via Marzolo 5 35131 Padova Italy
| | - Veronica Salmaso
- Molecular Modeling Section, MMS, Department of Pharmaceutical and Pharmacological Sciences; University of Padova; via Marzolo 5 35131 Padova Italy
| | - Mattia Sturlese
- Molecular Modeling Section, MMS, Department of Pharmaceutical and Pharmacological Sciences; University of Padova; via Marzolo 5 35131 Padova Italy
| | - Stefano Moro
- Molecular Modeling Section, MMS, Department of Pharmaceutical and Pharmacological Sciences; University of Padova; via Marzolo 5 35131 Padova Italy
| |
Collapse
|
73
|
Heifetz A, Southey M, Morao I, Townsend-Nicholson A, Bodkin MJ. Computational Methods Used in Hit-to-Lead and Lead Optimization Stages of Structure-Based Drug Discovery. Methods Mol Biol 2018; 1705:375-394. [PMID: 29188574 DOI: 10.1007/978-1-4939-7465-8_19] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
GPCR modeling approaches are widely used in the hit-to-lead (H2L) and lead optimization (LO) stages of drug discovery. The aims of these modeling approaches are to predict the 3D structures of the receptor-ligand complexes, to explore the key interactions between the receptor and the ligand and to utilize these insights in the design of new molecules with improved binding, selectivity or other pharmacological properties. In this book chapter, we present a brief survey of key computational approaches integrated with hierarchical GPCR modeling protocol (HGMP) used in hit-to-lead (H2L) and in lead optimization (LO) stages of structure-based drug discovery (SBDD). We outline the differences in modeling strategies used in H2L and LO of SBDD and illustrate how these tools have been applied in three drug discovery projects.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK. .,Division of Biosciences, Research Department of Structural and Molecular Biology, Institute of Structural and Molecular Biology, University College London, London, WC1E 6BT, UK.
| | - Michelle Southey
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Andrea Townsend-Nicholson
- Division of Biosciences, Research Department of Structural and Molecular Biology, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Mike J Bodkin
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| |
Collapse
|
74
|
Abstract
Advances in the structural biology of G-protein Coupled Receptors have resulted in a significant step forward in our understanding of how this important class of drug targets function at the molecular level. However, it has also become apparent that they are very dynamic molecules, and moreover, that the underlying dynamics is crucial in shaping the response to different ligands. Molecular dynamics simulations can provide unique insight into the dynamic properties of GPCRs in a way that is complementary to many experimental approaches. In this chapter, we describe progress in three distinct areas that are particularly difficult to study with other techniques: atomic level investigation of the conformational changes that occur when moving between the various states that GPCRs can exist in, the pathways that ligands adopt during binding/unbinding events and finally, the influence of lipids on the conformational dynamics of GPCRs.
Collapse
Affiliation(s)
- Naushad Velgy
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | - George Hedger
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | - Philip C Biggin
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK.
| |
Collapse
|
75
|
Chen D, Li Y, Zhao M, Tan W, Li X, Savidge T, Guo W, Fan X. Effective lead optimization targeting the displacement of bridging receptor–ligand water molecules. Phys Chem Chem Phys 2018; 20:24399-24407. [DOI: 10.1039/c8cp04118k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Enhancing the binding affinities of ligands by means of lead modifications that displace bridging water molecules at protein–ligand interfaces is an important and widely studied lead optimization strategy.
Collapse
Affiliation(s)
- Deliang Chen
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Yibao Li
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Mingming Zhao
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Wen Tan
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Xun Li
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Tor Savidge
- Department of Pathology & Immunology
- Baylor College of Medicine
- Houston
- USA
- Texas Children's Microbiome Center
| | - Wei Guo
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Xiaolin Fan
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| |
Collapse
|
76
|
Cournia Z, Allen B, Sherman W. Relative Binding Free Energy Calculations in Drug Discovery: Recent Advances and Practical Considerations. J Chem Inf Model 2017; 57:2911-2937. [PMID: 29243483 DOI: 10.1021/acs.jcim.7b00564] [Citation(s) in RCA: 395] [Impact Index Per Article: 56.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Accurate in silico prediction of protein-ligand binding affinities has been a primary objective of structure-based drug design for decades due to the putative value it would bring to the drug discovery process. However, computational methods have historically failed to deliver value in real-world drug discovery applications due to a variety of scientific, technical, and practical challenges. Recently, a family of approaches commonly referred to as relative binding free energy (RBFE) calculations, which rely on physics-based molecular simulations and statistical mechanics, have shown promise in reliably generating accurate predictions in the context of drug discovery projects. This advance arises from accumulating developments in the underlying scientific methods (decades of research on force fields and sampling algorithms) coupled with vast increases in computational resources (graphics processing units and cloud infrastructures). Mounting evidence from retrospective validation studies, blind challenge predictions, and prospective applications suggests that RBFE simulations can now predict the affinity differences for congeneric ligands with sufficient accuracy and throughput to deliver considerable value in hit-to-lead and lead optimization efforts. Here, we present an overview of current RBFE implementations, highlighting recent advances and remaining challenges, along with examples that emphasize practical considerations for obtaining reliable RBFE results. We focus specifically on relative binding free energies because the calculations are less computationally intensive than absolute binding free energy (ABFE) calculations and map directly onto the hit-to-lead and lead optimization processes, where the prediction of relative binding energies between a reference molecule and new ideas (virtual molecules) can be used to prioritize molecules for synthesis. We describe the critical aspects of running RBFE calculations, from both theoretical and applied perspectives, using a combination of retrospective literature examples and prospective studies from drug discovery projects. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative binding free energy simulations, with a focus on real-world drug discovery applications. We offer guidelines for improving the accuracy of RBFE simulations, especially for challenging cases, and emphasize unresolved issues that could be improved by further research in the field.
Collapse
Affiliation(s)
- Zoe Cournia
- Biomedical Research Foundation, Academy of Athens , 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Bryce Allen
- Silicon Therapeutics , 300 A Street, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Silicon Therapeutics , 300 A Street, Boston, Massachusetts 02210, United States
| |
Collapse
|
77
|
Ross GA, Bruce Macdonald HE, Cave-Ayland C, Cabedo Martinez AI, Essex JW. Replica-Exchange and Standard State Binding Free Energies with Grand Canonical Monte Carlo. J Chem Theory Comput 2017; 13:6373-6381. [PMID: 29091438 PMCID: PMC5729546 DOI: 10.1021/acs.jctc.7b00738] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
![]()
The
ability of grand canonical Monte Carlo (GCMC) to create and
annihilate molecules in a given region greatly aids the identification
of water sites and water binding free energies in protein cavities.
However, acceptance rates without the application of biased moves
can be low, resulting in large variations in the observed water occupancies.
Here, we show that replica-exchange of the chemical potential significantly
reduces the variance of the GCMC data. This improvement comes at a
negligible increase in computational expense when simulations comprise
of runs at different chemical potentials. Replica-exchange GCMC is
also found to substantially increase the precision of water binding
free energies as calculated with grand canonical integration, which
has allowed us to address a missing standard state correction.
Collapse
Affiliation(s)
- Gregory A Ross
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center , New York, New York 10065, United States
| | | | | | - Ana I Cabedo Martinez
- Department of Chemistry, University of Southampton , Southampton, SO17 1BJ, United Kingdom
| | - Jonathan W Essex
- Department of Chemistry, University of Southampton , Southampton, SO17 1BJ, United Kingdom
| |
Collapse
|
78
|
Bartl J, Palazzesi F, Parrinello M, Hommers L, Riederer P, Walitza S, Grünblatt E. The impact of methylphenidate and its enantiomers on dopamine synthesis and metabolism in vitro. Prog Neuropsychopharmacol Biol Psychiatry 2017; 79:281-288. [PMID: 28690202 DOI: 10.1016/j.pnpbp.2017.07.002] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/01/2017] [Revised: 06/02/2017] [Accepted: 07/03/2017] [Indexed: 12/19/2022]
Abstract
Methylphenidate (MPH), a psychostimulant, is an effective first-line treatment for the symptoms associated with Attention-Deficit/Hyperactivity Disorder (ADHD). Although most MPH formulations are composed of the racemic 1:1 mixture of the two enantiomers (d- and l-threo), converging lines of evidence indicate that d-threo MPH seems to be superior to the l-isomer. We aimed to investigate whether MPH racemic mixture or pure enantiomers influence the enzyme activity of tyrosine hydroxylase (TH), monoamine oxidase B (MAO-B), catechol-O-methyltransferase (COMT), and aldehyde dehydrogenase (ALDH) in vitro in homogenates of rat PC12 cells incubated with racemic, d- and l-threo MPH (1nM up to 100μM), or a vehicle for control. We could observe dose dependent enhancement of TH activity with d-threo MPH, probably due to its higher affinity to the enzyme, which we could confirm for d-threo versus l-threo MPH via docking and molecular dynamic simulations analysis. MAO-B enzyme activity was found to be enhanced when incubated with both d- and l-isomers but not with the racemic mixture. This conflicting result was hypothesized to be due to possible aggregation of the two enantiomers or other molecular conformations. Such a possible interaction was observed indirectly, when TH was incubated with constant d-threo MPH while increasing l-isomer (increasing total MPH concentrations). Hence, TH activity was slightly decreased with increased l-isomer. In conclusion, the current in vitro investigation points to the stereoselectivity of the investigated enzymes and pharmacological effects of MPH enantiomers.
Collapse
Affiliation(s)
- Jasmin Bartl
- Department of Child and Adolescent Psychiatry and Psychotherapy, University Hospital of Psychiatry Zurich, University of Zurich, Switzerland
| | - Ferruccio Palazzesi
- Department of Chemistry and Applied Biosciences, Eidgenössische Technische Hochschule (ETH) Zurich, 8093 Zurich, Switzerland; Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, 6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, Eidgenössische Technische Hochschule (ETH) Zurich, 8093 Zurich, Switzerland; Facoltà di Informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, 6900 Lugano, Switzerland
| | - Leif Hommers
- Center of Mental Health, Department of Psychiatry, Psychosomatics and Psychotherapy, University Hospital of Wuerzburg, Wuerzburg, Germany
| | - Peter Riederer
- Center of Mental Health, Department of Psychiatry, Psychosomatics and Psychotherapy, University Hospital of Wuerzburg, Wuerzburg, Germany; Unit of Psychiatry, University of Southern Denmark, Odense, Denmark
| | - Susanne Walitza
- Department of Child and Adolescent Psychiatry and Psychotherapy, University Hospital of Psychiatry Zurich, University of Zurich, Switzerland; Center for Integrative Human Physiology, University of Zurich, Zurich, Switzerland; Neuroscience Center Zurich, University of Zurich and ETH Zurich, Switzerland
| | - Edna Grünblatt
- Department of Child and Adolescent Psychiatry and Psychotherapy, University Hospital of Psychiatry Zurich, University of Zurich, Switzerland; Center for Integrative Human Physiology, University of Zurich, Zurich, Switzerland; Neuroscience Center Zurich, University of Zurich and ETH Zurich, Switzerland.
| |
Collapse
|
79
|
Verma R, Mitchell-Koch K. In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function. Catalysts 2017; 7:212. [PMID: 30464857 PMCID: PMC6241538 DOI: 10.3390/catal7070212] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.
Collapse
Affiliation(s)
- Rajni Verma
- Department of Chemistry, McKinley Hall, Wichita State University, 1845 Fairmount, Wichita, KS 67260-0051, USA
| | - Katie Mitchell-Koch
- Department of Chemistry, McKinley Hall, Wichita State University, 1845 Fairmount, Wichita, KS 67260-0051, USA
| |
Collapse
|
80
|
Xiao W, He Z, Sun M, Li S, Li H. Statistical Analysis, Investigation, and Prediction of the Water Positions in the Binding Sites of Proteins. J Chem Inf Model 2017; 57:1517-1528. [DOI: 10.1021/acs.jcim.6b00620] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Wei Xiao
- School
of Information Science and Engineering, East China University of Science and Technology, Shanghai 200237, China
- Shanghai
Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Zenghui He
- Shanghai
Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Meijian Sun
- Shanghai
Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Shiliang Li
- Shanghai
Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| | - Honglin Li
- School
of Information Science and Engineering, East China University of Science and Technology, Shanghai 200237, China
- Shanghai
Key Laboratory of New Drug Design, School of Pharmacy, East China University of Science and Technology, 130 Meilong Road, Shanghai 200237, China
| |
Collapse
|
81
|
Spyrakis F, Ahmed MH, Bayden AS, Cozzini P, Mozzarelli A, Kellogg GE. The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery. J Med Chem 2017; 60:6781-6827. [PMID: 28475332 DOI: 10.1021/acs.jmedchem.7b00057] [Citation(s) in RCA: 98] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.
Collapse
Affiliation(s)
- Francesca Spyrakis
- Dipartimento di Scienza e Tecnologia del Farmaco, Università degli Studi di Torino , Via Pietro Giuria 9, 10125 Torino, Italy
| | - Mostafa H Ahmed
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| | - Alexander S Bayden
- CMD Bioscience , 5 Science Park, New Haven, Connecticut 06511, United States
| | - Pietro Cozzini
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Modellistica Molecolare, Università degli Studi di Parma , Parco Area delle Scienze 59/A, 43121 Parma, Italy
| | - Andrea Mozzarelli
- Dipartimento di Scienze degli Alimenti e del Farmaco, Laboratorio di Biochimica, Università degli Studi di Parma , Parco Area delle Scienze 23/A, 43121 Parma, Italy.,Istituto di Biofisica, Consiglio Nazionale delle Ricerche , Via Moruzzi 1, 56124 Pisa, Italy
| | - Glen E Kellogg
- Department of Medicinal Chemistry & Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University , Richmond, Virginia 23298-0540, United States
| |
Collapse
|
82
|
Estimation of kinetic and thermodynamic ligand-binding parameters using computational strategies. Future Med Chem 2017; 9:507-523. [DOI: 10.4155/fmc-2016-0224] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Kinetic and thermodynamic ligand–protein binding parameters are gaining growing importance as key information to consider in drug discovery. The determination of the molecular structures, using particularly x-ray and NMR techniques, is crucial for understanding how a ligand recognizes its target in the final binding complex. However, for a better understanding of the recognition processes, experimental studies of ligand–protein interactions are needed. Even though several techniques can be used to investigate both thermodynamic and kinetic profiles for a ligand–protein complex, these procedures are very often laborious, time consuming and expensive. In the last 10 years, computational approaches have enormous potential in providing insights into each of the above effects and in parsing their contributions to the changes in both kinetic and thermodynamic binding parameters. The main purpose of this review is to summarize the state of the art of computational strategies for estimating the kinetic and thermodynamic parameters of a ligand–protein binding.
Collapse
|
83
|
Sridhar A, Ross GA, Biggin PC. Waterdock 2.0: Water placement prediction for Holo-structures with a pymol plugin. PLoS One 2017; 12:e0172743. [PMID: 28235019 PMCID: PMC5325533 DOI: 10.1371/journal.pone.0172743] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2016] [Accepted: 02/08/2017] [Indexed: 12/30/2022] Open
Abstract
Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.
Collapse
Affiliation(s)
- Akshay Sridhar
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - Gregory A. Ross
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| | - Philip C. Biggin
- Department of Biochemistry, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
84
|
Nikolaev DM, Emelyanov A, Boitsov VM, Panov MS, Ryazantsev MN. A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling. F1000Res 2017; 6:33. [PMID: 28435665 PMCID: PMC5381632 DOI: 10.12688/f1000research.10541.3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Accepted: 11/23/2017] [Indexed: 12/25/2022] Open
Abstract
It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.
Collapse
Affiliation(s)
- Dmitrii M Nikolaev
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Anton Emelyanov
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Vitaly M Boitsov
- Saint-Petersburg National Research Academic University of the Russian Academy of Science, St. Petersburg, Russian Federation
| | - Maxim S Panov
- Saint-Petersburg State University, St. Petersburg, Russian Federation
| | - Mikhail N Ryazantsev
- Saint-Petersburg State University, St. Petersburg, Russian Federation.,Saint-Petersburg Scientific Center of the Russian Academy of Sciences, St. Petersburg, Russian Federation
| |
Collapse
|
85
|
Pal RK, Haider K, Kaur D, Flynn W, Xia J, Levy RM, Taran T, Wickstrom L, Kurtzman T, Gallicchio E. A combined treatment of hydration and dynamical effects for the modeling of host-guest binding thermodynamics: the SAMPL5 blinded challenge. J Comput Aided Mol Des 2017; 31:29-44. [PMID: 27696239 PMCID: PMC5477994 DOI: 10.1007/s10822-016-9956-6] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2016] [Accepted: 08/25/2016] [Indexed: 01/02/2023]
Abstract
As part of the SAMPL5 blinded experiment, we computed the absolute binding free energies of 22 host-guest complexes employing a novel approach based on the BEDAM single-decoupling alchemical free energy protocol with parallel replica exchange conformational sampling and the AGBNP2 implicit solvation model specifically customized to treat the effect of water displacement as modeled by the Hydration Site Analysis method with explicit solvation. Initial predictions were affected by the lack of treatment of ionic charge screening, which is very significant for these highly charged hosts, and resulted in poor relative ranking of negatively versus positively charged guests. Binding free energies obtained with Debye-Hückel treatment of salt effects were in good agreement with experimental measurements. Water displacement effects contributed favorably and very significantly to the observed binding affinities; without it, the modeling predictions would have grossly underestimated binding. The work validates the implicit/explicit solvation approach employed here and it shows that comprehensive physical models can be effective at predicting binding affinities of molecular complexes requiring accurate treatment of conformational dynamics and hydration.
Collapse
Affiliation(s)
- Rajat Kumar Pal
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, New York, 11210, USA
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - Kamran Haider
- Department of Chemistry, Lehman College, The City University of New York, 250 Bedford Park Blvd. West, Bronx, New York, NY, 10468, USA
| | - Divya Kaur
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - William Flynn
- Center for Biophysics and Computational Biology, Institute of Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, PA, USA
- Department of Physics and Astronomy, Rutgers University, Piscataway, NJ, 08854, USA
| | - Junchao Xia
- Center for Biophysics and Computational Biology, Institute of Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Ronald M Levy
- Center for Biophysics and Computational Biology, Institute of Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Tetiana Taran
- Borough of Manhattan Community College, Department of Science, The City University of New York, 199 Chambers Street, New York, NY, 10007, USA
| | - Lauren Wickstrom
- Borough of Manhattan Community College, Department of Science, The City University of New York, 199 Chambers Street, New York, NY, 10007, USA
| | - Tom Kurtzman
- Department of Chemistry, Lehman College, The City University of New York, 250 Bedford Park Blvd. West, Bronx, New York, NY, 10468, USA
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA
| | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, New York, 11210, USA.
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY, 10016, USA.
| |
Collapse
|
86
|
McRobb FM, Negri A, Beuming T, Sherman W. Molecular dynamics techniques for modeling G protein-coupled receptors. Curr Opin Pharmacol 2016; 30:69-75. [DOI: 10.1016/j.coph.2016.07.001] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 06/28/2016] [Accepted: 07/03/2016] [Indexed: 11/17/2022]
|
87
|
Heifetz A, James T, Morao I, Bodkin MJ, Biggin PC. Guiding lead optimization with GPCR structure modeling and molecular dynamics. Curr Opin Pharmacol 2016; 30:14-21. [DOI: 10.1016/j.coph.2016.06.004] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Revised: 06/16/2016] [Accepted: 06/17/2016] [Indexed: 01/04/2023]
|
88
|
Morozenko A, Stuchebrukhov AA. Dowser++, a new method of hydrating protein structures. Proteins 2016; 84:1347-57. [PMID: 27273373 DOI: 10.1002/prot.25081] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2016] [Revised: 05/19/2016] [Accepted: 05/30/2016] [Indexed: 12/30/2022]
Abstract
A new method of hydrating protein structures, which we call Dowser++, is presented. The method is based on a semi-empirical modification of a popular program for protein hydration Dowser, and the usage of protocols AutoDock Vina, and WaterDock. The positions of water molecules predicted by Dowser++ were compared with experimental data for a set of 14 high-resolution crystal structures of oligopeptide-binding protein (OppA) containing a large number of resolved internal water molecules, as well as for the D- and K-channels of cytochrome c oxidase, and the recent data on PSII. Comparison is also made with the predictions of the original Dowser, and its improved version, Dowser+, described in our previous publication. We also present a model for quantitative estimation of the quality of water molecules placement made by a program, which includes an assumption of possible false negative data from the crystallographic analysis. The comparison of predictions made by Dowser++, Dowser and Dowser+ demonstrates significant improvement of predictive power of the new method. Proteins 2016; 84:1347-1357. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- A Morozenko
- Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California, 95616.
| | - A A Stuchebrukhov
- Department of Chemistry, University of California Davis, One Shields Avenue, Davis, California, 95616.
| |
Collapse
|
89
|
Bodnarchuk MS. Water, water, everywhere… It's time to stop and think. Drug Discov Today 2016; 21:1139-46. [DOI: 10.1016/j.drudis.2016.05.009] [Citation(s) in RCA: 62] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Revised: 04/15/2016] [Accepted: 05/13/2016] [Indexed: 12/11/2022]
|
90
|
Haider K, Wickstrom L, Ramsey S, Gilson MK, Kurtzman T. Enthalpic Breakdown of Water Structure on Protein Active-Site Surfaces. J Phys Chem B 2016; 120:8743-56. [PMID: 27169482 DOI: 10.1021/acs.jpcb.6b01094] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The principles underlying water reorganization around simple nonpolar solutes are well understood and provide the framework for the classical hydrophobic effect, whereby water molecules structure themselves around solutes so that they maintain favorable energetic contacts with both the solute and the other water molecules. However, for certain solute surface topographies, water molecules, due to their geometry and size, are unable to simultaneously maintain favorable energetic contacts with both the surface and neighboring water molecules. In this study, we analyze the solvation of ligand-binding sites for six structurally diverse proteins using hydration site analysis and measures of local water structure, in order to identify surfaces at which water molecules are unable to structure themselves in a way that maintains favorable enthalpy relative to bulk water. These surfaces are characterized by a high degree of enclosure, weak solute-water interactions, and surface constraints that induce unfavorable pair interactions between neighboring water molecules. Additionally, we find that the solvation of charged side chains in an active site generally results in favorable enthalpy but can also lead to pair interactions between neighboring water molecules that are significantly unfavorable relative to bulk water. We find that frustrated local structure can occur not only in apolar and weakly polar pockets, where overall enthalpy tends to be unfavorable, but also in charged pockets, where overall water enthalpy tends to be favorable. The characterization of local water structure in these terms may prove useful for evaluating the displacement of water from diverse protein active-site environments.
Collapse
Affiliation(s)
- Kamran Haider
- Department of Chemistry, Lehman College, The City University of New York , 250 Bedford Park Boulevard West, Bronx, New York 10468, United States
| | - Lauren Wickstrom
- Borough of Manhattan Community College, Department of Science, The City University of New York , 199 Chambers Street, New York, New York 10007, United States
| | - Steven Ramsey
- Department of Chemistry, Lehman College, The City University of New York , 250 Bedford Park Boulevard West, Bronx, New York 10468, United States.,Ph.D. Program in Biochemistry, The Graduate Center of The City University of New York , 365 Fifth Avenue, New York, New York 10016, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego , 9500 Gilman Drive, La Jolla, California 92093-0736, United States
| | - Tom Kurtzman
- Department of Chemistry, Lehman College, The City University of New York , 250 Bedford Park Boulevard West, Bronx, New York 10468, United States.,Ph.D. Program in Biochemistry, The Graduate Center of The City University of New York , 365 Fifth Avenue, New York, New York 10016, United States.,Ph.D. Program in Chemistry, The Graduate Center of The City University of New York , 365 Fifth Avenue, New York, New York 10016, United States
| |
Collapse
|
91
|
Spinozzi F, Mariani P, Ortore MG. Proteins in binary solvents. Biophys Rev 2016; 8:87-106. [PMID: 28510051 PMCID: PMC5425779 DOI: 10.1007/s12551-016-0193-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2015] [Accepted: 02/01/2016] [Indexed: 01/09/2023] Open
Abstract
Proteins in living organisms exist in complex aqueous solutions or embedded in membranes. In solution, proteins are surrounded by a tightly bound hydration layer, which is more ordered and less mobile than bulk water. As a consequence, water plays a major role in controlling protein structure stability, conformational flexibility, dynamics, and functionality, but it also appears that protein surface regulates the structuring of the surrounding water. The presence of cosolvents can modify the hydration layer characteristics and then the whole protein structural and dynamical properties. Because cytoplasm or biological liquids are complex solutions, the knowledge of the solvation shell characteristics in mixed solvents should be considered as a crucial step in describing biological processes at molecular level. This review reports on recent studies on the structural and thermodynamic properties of model proteins dissolved in binary solvent mixtures by small-angle neutron scattering (SANS) and differential scanning microcalorimetry (DSC) techniques. We will show that contrast variation SANS experiments allow to acquire a direct knowledge of both protein structure and protein solvation shell (in terms of low-resolution shape and solvent/cosolvent composition), while DSC experiments provide information on all the relevant thermodynamic properties. We will focus on two main points. First, an extended description of the thermodynamic model used to define the equilibria between water and cosolvent molecules in the protein solvation shell will be presented. Second, the determination of the peculiar characteristics of the protein solvation layer, which will be illustrated by considering different systems. As a conclusion, we will show that the investigation of structure and thermodynamics of proteins in binary aqueous mixtures is an important way to understand the role of hydration in protein stability and activity.
Collapse
Affiliation(s)
- Francesco Spinozzi
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy.
| | - Paolo Mariani
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy
| | - Maria Grazia Ortore
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy
| |
Collapse
|
92
|
Hydration of proteins and nucleic acids: Advances in experiment and theory. A review. Biochim Biophys Acta Gen Subj 2016; 1860:1821-35. [PMID: 27241846 DOI: 10.1016/j.bbagen.2016.05.036] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2015] [Revised: 05/20/2016] [Accepted: 05/26/2016] [Indexed: 11/21/2022]
Abstract
BACKGROUND Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics. SCOPE OF REVIEW This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties. MAJOR CONCLUSIONS In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration. GENERAL SIGNIFICANCE An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.
Collapse
|
93
|
Fierro A, Edmondson DE, Celis-Barros C, Rebolledo-Fuentes M, Zapata-Torres G. Why p-OMe- and p-Cl-β-Methylphenethylamines Display Distinct Activities upon MAO-B Binding. PLoS One 2016; 11:e0154989. [PMID: 27152414 PMCID: PMC4859490 DOI: 10.1371/journal.pone.0154989] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2016] [Accepted: 04/23/2016] [Indexed: 12/16/2022] Open
Abstract
Despite their structural and chemical commonalities, p-chloro-β-methylphenethylamine and p-methoxy-β-methylphenethylamine display distinct inhibitory and substrate activities upon MAO-B binding. Density Functional Theory (DFT) quantum chemical calculations reveal that β-methylation and para-substitution underpin the observed activities sustained by calculated transition state energy barriers, attained conformations and key differences in their interactions in the enzyme’s substrate binding site. Although both compounds meet substrate requirements, it is clear that β-methylation along with the physicochemical features of the para-substituents on the aromatic ring determine the activity of these compounds upon binding to the MAO B-isoform. While data for a larger set of compounds might lend generality to our conclusions, our experimental and theoretical results strongly suggest that the contrasting activities displayed depend on the conformations adopted by these compounds when they bind to the enzyme.
Collapse
Affiliation(s)
- Angélica Fierro
- Facultad de Química, Departamento de Química Orgánica, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Dale E. Edmondson
- Departments of Chemistry and Biochemistry, Emory University, Atlanta, GA, United States of America
| | - Cristian Celis-Barros
- Universidad Andres Bello, Facultad de Ciencias Exactas, Departamento de Ciencias Quimicas, Santiago, Chile
- * E-mail:
| | | | - Gerald Zapata-Torres
- Faculty of Chemical and Pharmaceutical Sciences, University of Chile, Santiago, Chile
- * E-mail:
| |
Collapse
|
94
|
Robinson D, Bertrand T, Carry JC, Halley F, Karlsson A, Mathieu M, Minoux H, Perrin MA, Robert B, Schio L, Sherman W. Differential Water Thermodynamics Determine PI3K-Beta/Delta Selectivity for Solvent-Exposed Ligand Modifications. J Chem Inf Model 2016; 56:886-94. [DOI: 10.1021/acs.jcim.5b00641] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniel Robinson
- Schrodinger, 120 W 45th St, New York, New York 10036, United States
| | | | | | | | | | | | | | | | | | | | - Woody Sherman
- Schrodinger, 120 W 45th St, New York, New York 10036, United States
| |
Collapse
|
95
|
Chen D, Oezguen N, Urvil P, Ferguson C, Dann SM, Savidge TC. Regulation of protein-ligand binding affinity by hydrogen bond pairing. SCIENCE ADVANCES 2016; 2:e1501240. [PMID: 27051863 PMCID: PMC4820369 DOI: 10.1126/sciadv.1501240] [Citation(s) in RCA: 360] [Impact Index Per Article: 45.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2015] [Accepted: 02/05/2016] [Indexed: 05/20/2023]
Abstract
Hydrogen (H)-bonds potentiate diverse cellular functions by facilitating molecular interactions. The mechanism and the extent to which H-bonds regulate molecular interactions are a largely unresolved problem in biology because the H-bonding process continuously competes with bulk water. This interference may significantly alter our understanding of molecular function, for example, in the elucidation of the origin of enzymatic catalytic power. We advance this concept by showing that H-bonds regulate molecular interactions via a hitherto unappreciated donor-acceptor pairing mechanism that minimizes competition with water. On the basis of theoretical and experimental correlations between H-bond pairings and their effects on ligand binding affinity, we demonstrate that H-bonds enhance receptor-ligand interactions when both the donor and acceptor have either significantly stronger or significantly weaker H-bonding capabilities than the hydrogen and oxygen atoms in water. By contrast, mixed strong-weak H-bond pairings decrease ligand binding affinity due to interference with bulk water, offering mechanistic insight into why indiscriminate strengthening of receptor-ligand H-bonds correlates poorly with experimental binding affinity. Further support for the H-bond pairing principle is provided by the discovery and optimization of lead compounds targeting dietary melamine and Clostridium difficile toxins, which are not realized by traditional drug design methods. Synergistic H-bond pairings have therefore evolved in the natural design of high-affinity binding and provide a new conceptual framework to evaluate the H-bonding process in biological systems. Our findings may also guide wider applications of competing H-bond pairings in lead compound design and in determining the origin of enzymatic catalytic power.
Collapse
Affiliation(s)
- Deliang Chen
- Key Laboratory of Organo-Pharmaceutical Chemistry of Jiangxi Province, Chemistry and Chemical Engineering College, Gannan Normal University, Ganzhou, Jiangxi 341000, P. R. China
| | - Numan Oezguen
- Department of Pathology and Immunology, Baylor College of Medicine, Houston, TX 77030, USA
- Texas Children’s Microbiome Center, Texas Children’s Hospital, Houston, TX 77030, USA
| | - Petri Urvil
- Department of Pathology and Immunology, Baylor College of Medicine, Houston, TX 77030, USA
- Texas Children’s Microbiome Center, Texas Children’s Hospital, Houston, TX 77030, USA
| | | | - Sara M. Dann
- Department of Internal Medicine, University of Texas Medical Branch, Galveston, TX 77555, USA
| | - Tor C. Savidge
- Department of Pathology and Immunology, Baylor College of Medicine, Houston, TX 77030, USA
- Texas Children’s Microbiome Center, Texas Children’s Hospital, Houston, TX 77030, USA
| |
Collapse
|
96
|
Sridhar A, Johnston AJ, Varathan L, McLain SE, Biggin PC. The solvation structure of alprazolam. Phys Chem Chem Phys 2016; 18:22416-25. [DOI: 10.1039/c6cp02645a] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Alprazolam is a benzodiazepine that is commonly prescribed for the treatment of anxiety and other related disorders.
Collapse
|
97
|
Heifetz A, Chudyk EI, Gleave L, Aldeghi M, Cherezov V, Fedorov DG, Biggin PC, Bodkin MJ. The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions. J Chem Inf Model 2015; 56:159-72. [PMID: 26642258 DOI: 10.1021/acs.jcim.5b00644] [Citation(s) in RCA: 79] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Ewa I. Chudyk
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Laura Gleave
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Matteo Aldeghi
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Vadim Cherezov
- Department
of Chemistry, Bridge Institute, University of Southern California, Los Angeles, California 90089, United States
- Laboratory
for Structural Biology of GPCRs, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Dmitri G. Fedorov
- NMRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Philip C. Biggin
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Mike J. Bodkin
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| |
Collapse
|
98
|
Ross GA, Bodnarchuk MS, Essex JW. Water Sites, Networks, And Free Energies with Grand Canonical Monte Carlo. J Am Chem Soc 2015; 137:14930-43. [PMID: 26509924 DOI: 10.1021/jacs.5b07940] [Citation(s) in RCA: 85] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Water molecules play integral roles in the formation of many protein-ligand complexes, and recent computational efforts have been focused on predicting the thermodynamic properties of individual waters and how they may be exploited in rational drug design. However, when water molecules form highly coupled hydrogen-bonding networks, there is, as yet, no method that can rigorously calculate the free energy to bind the entire network or assess the degree of cooperativity between waters. In this work, we report theoretical and methodological developments to the grand canonical Monte Carlo simulation technique. Central to our results is a rigorous equation that can be used to calculate efficiently the binding free energies of water networks of arbitrary size and complexity. Using a single set of simulations, our methods can locate waters, estimate their binding affinities, capture the cooperativity of the water network, and evaluate the hydration free energy of entire protein binding sites. Our techniques have been applied to multiple test systems and compare favorably to thermodynamic integration simulations and experimental data. The implications of these methods in drug design are discussed.
Collapse
Affiliation(s)
- Gregory A Ross
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| | - Michael S Bodnarchuk
- School of Mechanical Engineering, Imperial College London , Exhibition Road, London, SW1 2AZ, United Kingdom
| | - Jonathan W Essex
- School of Chemistry, University of Southampton , Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
99
|
Setny P. Prediction of Water Binding to Protein Hydration Sites with a Discrete, Semiexplicit Solvent Model. J Chem Theory Comput 2015; 11:5961-72. [PMID: 26642995 DOI: 10.1021/acs.jctc.5b00839] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Buried water molecules are ubiquitous in protein structures and are found at the interface of most protein-ligand complexes. Determining their distribution and thermodynamic effect is a challenging yet important task, of great of practical value for the modeling of biomolecular structures and their interactions. In this study, we present a novel method aimed at the prediction of buried water molecules in protein structures and estimation of their binding free energies. It is based on a semiexplicit, discrete solvation model, which we previously introduced in the context of small molecule hydration. The method is applicable to all macromolecular structures described by a standard all-atom force field, and predicts complete solvent distribution within a single run with modest computational cost. We demonstrate that it indicates positions of buried hydration sites, including those filled by more than one water molecule, and accurately differentiates them from sterically accessible to water but void regions. The obtained estimates of water binding free energies are in fair agreement with reference results determined with the double decoupling method.
Collapse
Affiliation(s)
- Piotr Setny
- Centre of New Technologies, University of Warsaw , Banacha 2c, 02-097 Warsaw, Poland
| |
Collapse
|
100
|
Bradley AR, Wall ID, von Delft F, Green DVS, Deane CM, Marsden BD. WONKA: objective novel complex analysis for ensembles of protein-ligand structures. J Comput Aided Mol Des 2015; 29:963-73. [PMID: 26387008 PMCID: PMC4621702 DOI: 10.1007/s10822-015-9866-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2015] [Accepted: 09/04/2015] [Indexed: 01/16/2023]
Abstract
WONKA is a tool for the systematic analysis of an ensemble of protein-ligand structures. It makes the identification of conserved and unusual features within such an ensemble straightforward. WONKA uses an intuitive workflow to process structural co-ordinates. Ligand and protein features are summarised and then presented within an interactive web application. WONKA's power in consolidating and summarising large amounts of data is described through the analysis of three bromodomain datasets. Furthermore, and in contrast to many current methods, WONKA relates analysis to individual ligands, from which we find unusual and erroneous binding modes. Finally the use of WONKA as an annotation tool to share observations about structures is demonstrated. WONKA is freely available to download and install locally or can be used online at http://wonka.sgc.ox.ac.uk.
Collapse
Affiliation(s)
- A R Bradley
- SGC, Nuffield Department of Medicine, University of Oxford, Old Road Campus Research Building, Roosevelt Drive, Headington, Oxford, OX3 7DQ, UK
- Oxford Protein Informatics Group, Department of Statistics, University of Oxford, 1 South Parks Road, Oxford, OX1 TG, UK
| | - I D Wall
- Computational & Structural Chemistry, GlaxoSmithKline, Medicines Research Centre, Gunnels Wood Road, Stevenage, Hertfordshire, SG1 2NY, UK
| | - F von Delft
- SGC, Nuffield Department of Medicine, University of Oxford, Old Road Campus Research Building, Roosevelt Drive, Headington, Oxford, OX3 7DQ, UK
- Diamond Light Source Ltd, Harwell Science and Innovation Campus, Didcot, OX11 0QX, UK
- Department of Biochemistry, University of Johannesburg, Aukland Park, 2006, South Africa
| | - D V S Green
- Computational & Structural Chemistry, GlaxoSmithKline, Medicines Research Centre, Gunnels Wood Road, Stevenage, Hertfordshire, SG1 2NY, UK
| | - C M Deane
- Oxford Protein Informatics Group, Department of Statistics, University of Oxford, 1 South Parks Road, Oxford, OX1 TG, UK
| | - B D Marsden
- SGC, Nuffield Department of Medicine, University of Oxford, Old Road Campus Research Building, Roosevelt Drive, Headington, Oxford, OX3 7DQ, UK.
- Kennedy Institute of Rheumatology, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences, University of Oxford, Roosevelt Drive, Headington, Oxford, OX3 7FY, UK.
| |
Collapse
|