1
|
Eberhardt J, Forli S. WaterKit: Thermodynamic Profiling of Protein Hydration Sites. J Chem Theory Comput 2023; 19:2535-2556. [PMID: 37094087 PMCID: PMC10732097 DOI: 10.1021/acs.jctc.2c01087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023]
Abstract
Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.
Collapse
Affiliation(s)
- Jerome Eberhardt
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| |
Collapse
|
2
|
Banerjee M, Hatial I, Keegan BM, Blagg BSJ. Assay design and development strategies for finding Hsp90 inhibitors and their role in human diseases. Pharmacol Ther 2021; 221:107747. [PMID: 33245994 PMCID: PMC8744950 DOI: 10.1016/j.pharmthera.2020.107747] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Revised: 11/11/2020] [Accepted: 11/16/2020] [Indexed: 12/30/2022]
Abstract
Heat shock protein 90 (Hsp90) is a molecular chaperone that facilitates the maturation of its client proteins including protein kinases, transcription factors, and steroid hormone receptors which are structurally and functionally diverse. These client proteins are involved in various cellular signaling pathways, and Hsp90 is implicated in various human diseases including cancer, inflammation, and diseases associated with protein misfolding; thus making Hsp90 a promising target for drug discovery. Some of its client proteins are well-known cancer targets. Instead of targeting these client proteins individually, however, targeting Hsp90 is more practical for cancer drug development. Efforts have been invested in recognizing potential drugs for clinical use that inhibit Hsp90 activity and result in the prevention of Hsp90 client maturation and dampening of subsequent signaling cascades. Here, we discuss current assays and technologies used to find and characterize Hsp90 inhibitors that include biophysical, biochemical, cell-based assays and computational modeling. This review highlights recent discoveries that N-terminal isoform-selective compounds and inhibitors that target the Hsp90 C-terminus that may offer the potential to overcome some of the detriments observed with pan Hsp90 inhibitors. The tools and assays summarized in this review should be used to develop Hsp90-targeting drugs with high specificity, potency, and drug-like properties that may prove immensely useful in the clinic.
Collapse
Affiliation(s)
- Monimoy Banerjee
- Department of Chemistry & Biochemistry, Warren Family Research Center for Drug Discovery and Development, 305 McCourtney Hall, University of Norte Dame, Norte Dame, IN 46656, USA
| | - Ishita Hatial
- Department of Chemistry & Biochemistry, Warren Family Research Center for Drug Discovery and Development, 305 McCourtney Hall, University of Norte Dame, Norte Dame, IN 46656, USA
| | - Bradley M Keegan
- Department of Chemistry & Biochemistry, Warren Family Research Center for Drug Discovery and Development, 305 McCourtney Hall, University of Norte Dame, Norte Dame, IN 46656, USA
| | - Brian S J Blagg
- Department of Chemistry & Biochemistry, Warren Family Research Center for Drug Discovery and Development, 305 McCourtney Hall, University of Norte Dame, Norte Dame, IN 46656, USA.
| |
Collapse
|
3
|
Shadrack DM, Swai HS, Hassanali A. A computational study on the role of water and conformational fluctuations in Hsp90 in response to inhibitors. J Mol Graph Model 2019; 96:107510. [PMID: 31877402 DOI: 10.1016/j.jmgm.2019.107510] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2019] [Revised: 12/04/2019] [Accepted: 12/11/2019] [Indexed: 10/25/2022]
Abstract
Molecular chaperone Heat Shock Protein 90 (Hsp90) represents an interesting chemotherapeutic target for cancer treatments as it plays a role in cancer proliferation. Thus, continued effort to identify novel inhibitors of this target is an important task. Drug design using computational approach has gained significant attention in recent years. This work aims to propose docking protocols to re-purpose FDA-approved drugs targeting Hsp90. Sensitivity of results to different docking protocols such apo, holo and receptor ensembles (relaxed complex) structures, the role of water and conformational changes of Hsp90, are described. We show that the protein conformation and water have effects on drug binding. Holo relaxed complex receptors ensembles improves the binding energy of ligands to the protein. We also compare and contrast structural stability of three drugs namely: ezetimibe, pitavastatin and vilazodon in the Hsp90 protein. The results obtained serves as a possible basis towards developing Hsp90 inhibitors.
Collapse
Affiliation(s)
- Daniel M Shadrack
- Department of Health and Biomedical Sciences, School of Life Science and Bioengieering, The Nelson Mandela African Institution of Science and Technology, P.O.Box 447, Arusha, Tanzania; International Centre for Theoretical Physics, Strada Costiera, 11, 34151, Trieste, Italy; Department of Chemistry, FaNAS, St John's University of Tanzania, P.O.Box 47, Dodoma, Tanzania.
| | - Hulda S Swai
- Department of Health and Biomedical Sciences, School of Life Science and Bioengieering, The Nelson Mandela African Institution of Science and Technology, P.O.Box 447, Arusha, Tanzania
| | - Ali Hassanali
- International Centre for Theoretical Physics, Strada Costiera, 11, 34151, Trieste, Italy.
| |
Collapse
|
4
|
Improved pose and affinity predictions using different protocols tailored on the basis of data availability. J Comput Aided Mol Des 2016; 30:817-828. [DOI: 10.1007/s10822-016-9982-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Accepted: 09/28/2016] [Indexed: 12/22/2022]
|
5
|
Santos-Martins D. Interaction with specific HSP90 residues as a scoring function: validation in the D3R Grand Challenge 2015. J Comput Aided Mol Des 2016; 30:731-742. [PMID: 27549813 DOI: 10.1007/s10822-016-9943-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2016] [Accepted: 08/17/2016] [Indexed: 01/19/2023]
Abstract
Here is reported the development of a novel scoring function that performs remarkably well at identifying the native binding pose of a subset of HSP90 inhibitors containing aminopyrimidine or resorcinol based scaffolds. This scoring function is called PocketScore, and consists of the interaction energy between a ligand and three residues in the binding pocket: Asp93, Thr184 and a water molecule. We integrated PocketScore into a molecular docking workflow, and used it to participate in the Drug Design Data Resource (D3R) Grand Challenge 2015 (GC2015). PocketScore was able to rank 180 molecules of the GC2015 according to their binding affinity with satisfactory performance. These results indicate that the specific residues considered by PocketScore are determinant to properly model the interaction between HSP90 and its subset of inhibitors containing aminopyrimidine or resorcinol based scaffolds. Moreover, the development of PocketScore aimed at improving docking power while neglecting the prediction of binding affinities, suggesting that accurate identification of native binding poses is a determinant factor for the performance of virtual screens.
Collapse
Affiliation(s)
- Diogo Santos-Martins
- UCIBIO, REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto, 4169-007, Porto, Portugal.
| |
Collapse
|
6
|
Gupta SD, Bommaka MK, Mazaira GI, Galigniana MD, Subrahmanyam CVS, Gowrishankar NL, Raghavendra NM. Molecular docking study, synthesis and biological evaluation of Mannich bases as Hsp90 inhibitors. Int J Biol Macromol 2015; 80:253-9. [DOI: 10.1016/j.ijbiomac.2015.06.039] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2015] [Revised: 06/18/2015] [Accepted: 06/20/2015] [Indexed: 12/17/2022]
|
7
|
Ioannidis D, Papadopoulos GE, Anastassopoulos G, Kortsaris A, Anagnostopoulos K. Structural properties and interaction energies affecting drug design. An approach combining molecular simulations, statistics, interaction energies and neural networks. Comput Biol Chem 2015; 56:7-12. [DOI: 10.1016/j.compbiolchem.2015.02.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2014] [Revised: 02/17/2015] [Accepted: 02/22/2015] [Indexed: 11/26/2022]
|
8
|
Anighoro A, Stumpfe D, Heikamp K, Beebe K, Neckers LM, Bajorath J, Rastelli G. Computational polypharmacology analysis of the heat shock protein 90 interactome. J Chem Inf Model 2015; 55:676-86. [PMID: 25686391 DOI: 10.1021/ci5006959] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The design of a single drug molecule that is able to simultaneously and specifically interact with multiple biological targets is gaining major consideration in drug discovery. However, the rational design of drugs with a desired polypharmacology profile is still a challenging task, especially when these targets are distantly related or unrelated. In this work, we present a computational approach aimed at the identification of suitable target combinations for multitarget drug design within an ensemble of biologically relevant proteins. The target selection relies on the analysis of activity annotations present in molecular databases and on ligand-based virtual screening. A few target combinations were also inspected with structure-based methods to demonstrate that the identified dual-activity compounds are able to bind target combinations characterized by remote binding site similarities. Our approach was applied to the heat shock protein 90 (Hsp90) interactome, which contains several targets of key importance in cancer. Promising target combinations were identified, providing a basis for the computational design of compounds with dual activity. The approach may be used on any ensemble of proteins of interest for which known inhibitors are available.
Collapse
Affiliation(s)
- Andrew Anighoro
- †Life Sciences Department, University of Modena and Reggio Emilia, Via Campi 183, 41125 Modena, Italy
| | - Dagmar Stumpfe
- ‡Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstr. 2, D-53113 Bonn, Germany
| | - Kathrin Heikamp
- ‡Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstr. 2, D-53113 Bonn, Germany
| | - Kristin Beebe
- §Urological Oncology Branch, Center for Cancer Research, National Cancer Institute, Bethesda, Maryland 20892, United States
| | - Leonard M Neckers
- §Urological Oncology Branch, Center for Cancer Research, National Cancer Institute, Bethesda, Maryland 20892, United States
| | - Jürgen Bajorath
- ‡Department of Life Science Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstr. 2, D-53113 Bonn, Germany
| | - Giulio Rastelli
- †Life Sciences Department, University of Modena and Reggio Emilia, Via Campi 183, 41125 Modena, Italy
| |
Collapse
|
9
|
Sun H, Zhao L, Peng S, Huang N. Incorporating replacement free energy of binding-site waters in molecular docking. Proteins 2014; 82:1765-76. [DOI: 10.1002/prot.24530] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 01/17/2014] [Accepted: 01/28/2014] [Indexed: 12/24/2022]
Affiliation(s)
- Hanzi Sun
- College of Life Sciences; Beijing Normal University; Beijing 100875 China
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Lifeng Zhao
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Shiming Peng
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| | - Niu Huang
- College of Life Sciences; Beijing Normal University; Beijing 100875 China
- National Institute of Biological Sciences, Beijing, Zhongguancun Life Science Park; Beijing 102206 China
| |
Collapse
|
10
|
Dutta Gupta S, Snigdha D, Mazaira GI, Galigniana MD, Subrahmanyam CVS, Gowrishankar NL, Raghavendra NM. Molecular docking study, synthesis and biological evaluation of Schiff bases as Hsp90 inhibitors. Biomed Pharmacother 2014; 68:369-76. [PMID: 24486109 DOI: 10.1016/j.biopha.2014.01.003] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2013] [Accepted: 01/01/2014] [Indexed: 11/25/2022] Open
Abstract
Heat shock protein 90 (Hsp90) is an emerging attractive target for the discovery of novel cancer therapeutic agents. Docking methods are powerful in silico tools for lead generation and optimization. In our mission to rationally develop novel effective small molecules against Hsp90, we predicted the potency of our designed compounds by Sybyl surflex Geom X docking method. The results of the above studies revealed that Schiff bases derived from 2,4-dihydroxy benzaldehyde/5-chloro-2,4-dihydroxy benzaldehyde demonstrated effective binding with the protein. Subsequently, a few of them were synthesized (1-10) and characterized by IR, (1)HNMR and mass spectral analysis. The synthesized molecules were evaluated for their potential to suppress Hsp90 ATPase activity by Malachite green assay. The anticancer studies were performed by 3-(4,5-dimethythiazol- 2-yl)-2,5-diphenyl tetrazolium bromide (MTT) assay method. The software generated results was in satisfactory agreement with the evaluated biological activity.
Collapse
Affiliation(s)
- Sayan Dutta Gupta
- Department of Pharmaceutical Chemistry, Gokaraju Rangaraju College of Pharmacy, Osmania University, Bachupally, Hyderabad, India; R&D centre, Department of Pharmaceutical Sciences, Jawaharlal Nehru Technological University, Hyderabad, India.
| | - D Snigdha
- Department of Pharmaceutical Chemistry, Gokaraju Rangaraju College of Pharmacy, Osmania University, Bachupally, Hyderabad, India
| | - Gisela I Mazaira
- Department of Biological Chemistry, Faculty of Natural Sciences, University of Buenos Aires, Buenos Aires, Argentina
| | - Mario D Galigniana
- Department of Biological Chemistry, Faculty of Natural Sciences, University of Buenos Aires, Buenos Aires, Argentina; Institute of Experimental Biology and Medicine-CONICET, Buenos Aires, Argentina
| | - C V S Subrahmanyam
- Department of Pharmaceutical Chemistry, Gokaraju Rangaraju College of Pharmacy, Osmania University, Bachupally, Hyderabad, India
| | - N L Gowrishankar
- Swami Vivekananda Institute of Pharmaceutical Sciences, Nalgonda, Andhrapradesh, India
| | - N M Raghavendra
- Department of Pharmaceutical Chemistry, Gokaraju Rangaraju College of Pharmacy, Osmania University, Bachupally, Hyderabad, India
| |
Collapse
|
11
|
Haider K, Huggins DJ. Combining solvent thermodynamic profiles with functionality maps of the Hsp90 binding site to predict the displacement of water molecules. J Chem Inf Model 2013; 53:2571-86. [PMID: 24070451 PMCID: PMC3840717 DOI: 10.1021/ci4003409] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Intermolecular interactions in the aqueous phase must compete with the interactions between the two binding partners and their solvating water molecules. In biological systems, water molecules in protein binding sites cluster at well-defined hydration sites and can form strong hydrogen-bonding interactions with backbone and side-chain atoms. Displacement of such water molecules is only favorable when the ligand can form strong compensating hydrogen bonds. Conversely, water molecules in hydrophobic regions of protein binding sites make only weak interactions, and the requirements for favorable displacement are less stringent. The propensity of water molecules for displacement can be identified using inhomogeneous fluid solvation theory (IFST), a statistical mechanical method that decomposes the solvation free energy of a solute into the contributions from different spatial regions and identifies potential binding hotspots. In this study, we employed IFST to study the displacement of water molecules from the ATP binding site of Hsp90, using a test set of 103 ligands. The predicted contribution of a hydration site to the hydration free energy was found to correlate well with the observed displacement. Additionally, we investigated if this correlation could be improved by using the energetic scores of favorable probe groups binding at the location of hydration sites, derived from a multiple copy simultaneous search (MCSS) method. The probe binding scores were not highly predictive of the observed displacement and did not improve the predictivity when used in combination with IFST-based hydration free energies. The results show that IFST alone can be used to reliably predict the observed displacement of water molecules in Hsp90. However, MCSS can augment IFST calculations by suggesting which functional groups should be used to replace highly displaceable water molecules. Such an approach could be very useful in improving the hit-to-lead process for new drug targets.
Collapse
Affiliation(s)
- Kamran Haider
- Department of Biology, Syed Babar Ali School of Science and Engineering, Lahore University of Management Sciences , Lahore, 54792, Pakistan
| | | |
Collapse
|
12
|
Huggins DJ, Tidor B. Systematic placement of structural water molecules for improved scoring of protein-ligand interactions. Protein Eng Des Sel 2011; 24:777-89. [PMID: 21771870 PMCID: PMC3170077 DOI: 10.1093/protein/gzr036] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2011] [Revised: 06/03/2011] [Accepted: 06/15/2011] [Indexed: 11/13/2022] Open
Abstract
Structural water molecules are found in many protein-ligand complexes. They are known to be vital in mediating hydrogen-bonding interactions and, in some cases, key for facilitating tight binding. It is thus very important to consider water molecules when attempting to model protein-ligand interactions for cognate ligand identification, virtual screening and drug design. While the rigid treatment of water molecules present in structures is feasible, the more relevant task of treating all possible positions and orientations of water molecules with each possible ligand pose is computationally daunting. Current methods in molecular docking provide partial treatment for such water molecules, with modest success. Here we describe a new method employing dead-end elimination to place water molecules within a binding site, bridging interactions between protein and ligand. Dead-end elimination permits a thorough, though still incomplete, treatment of water placement. The results show that this method is able to place water molecules correctly within known complexes and to create physically reasonable hydrogen bonds. The approach has also been incorporated within an inverse molecular design approach, to model a variety of compounds in the process of de novo ligand design. The inclusion of structural water molecules, combined with ranking based on the electrostatic contribution to binding affinity, improves a number of otherwise poor energetic predictions.
Collapse
Affiliation(s)
- David J. Huggins
- Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA 02139–4307, USA
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139–4307, USA
| | - Bruce Tidor
- Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA 02139–4307, USA
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139–4307, USA
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA 02139–4307, USA
| |
Collapse
|
13
|
Lu SY, Jiang YJ, Lv J, Zou JW, Wu TX. Role of bridging water molecules in GSK3β-inhibitor complexes: insights from QM/MM, MD, and molecular docking studies. J Comput Chem 2011; 32:1907-18. [PMID: 21469159 DOI: 10.1002/jcc.21775] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2010] [Revised: 01/02/2011] [Accepted: 01/21/2011] [Indexed: 02/01/2023]
Abstract
The role of water molecules is increasingly gaining interest in drug design, and several studies have highlighted their paramount contributions to the specificity and the affinity of ligand binding. In this study, we employ the two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations, molecular dynamics (MD) simulations, and molecular docking studies to investigate the effect of bridging water molecules at the GSK3β-inhibitors interfaces. The results obtained from the ONIOM geometry optimization and AIM analysis corroborated the presence of bridging water molecules that form hydrogen bonds with protein side chain of Thr138 and/or backbone of Gln185, and mediate interactions with inhibitors in the 10 selected GSK3β-inhibitor complexes. Subsequently, MD simulations carried out on a representative system of 1R0E demonstrated that the bridging water molecule is stable at the GSK3β-inhibitor interface and appears to contribute to the stability of the protein-inhibitor interactions. Furthermore, molecular docking studies of GSK3β-inhibitor complexes indicated that the inhibitors can increase binding affinities and the better docked conformation of inhibitors can be obtained by inclusion of the bridging water molecules, especially for the flexible inhibitors, in docking experiments into individual protein conformations. Our results elucidate the importance of bridging water molecules at the GSK3β-inhibitor interfaces and suggest that they might prove useful in rational drug design.
Collapse
Affiliation(s)
- Shao-Yong Lu
- Key Laboratory for Molecular Design and Nutrition Engineering, Ningbo Institute of Technology, Zhejiang University, Ningbo 315104, China
| | | | | | | | | |
Collapse
|
14
|
Florová P, Sklenovský P, Banáš P, Otyepka M. Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact. J Chem Theory Comput 2010; 6:3569-79. [PMID: 26617103 DOI: 10.1021/ct1003687] [Citation(s) in RCA: 92] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Conventional molecular dynamics simulations on 50 ns to 1 μs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a β-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified. Chignolin was highly solvated in TIP5P, whereas reduced specific solvation was found in the TIP4P, SPC/E, and TIP3P models for kinetic, thermodynamic, and both kinetic and thermodynamic reasons, respectively. The differences in specific solvation did not influence the dynamics of structured parts of the folded peptide. However, substantial differences between TIP5P and the other models were observed in the dynamics of unfolded chignolin, stability of salt bridges, and specific solvation of the backbone carbonyls of EK-peptide. Thus, we conclude that the choice of water model may affect the dynamics of flexible parts of proteins that are solvent-exposed. On the other hand, all water models should perform similarly for well-structured folded protein regions. The merits of the TIP3P model include its high and overestimated mobility, which accelerates simulation processes and thus effectively increases sampling.
Collapse
Affiliation(s)
- Petra Florová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Sklenovský
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, tr. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
15
|
Huggins DJ, McKenzie GJ, Robinson DD, Narváez AJ, Hardwick B, Roberts-Thomson M, Venkitaraman AR, Grant GH, Payne MC. Computational analysis of phosphopeptide binding to the polo-box domain of the mitotic kinase PLK1 using molecular dynamics simulation. PLoS Comput Biol 2010; 6:e1000880. [PMID: 20711360 PMCID: PMC2920843 DOI: 10.1371/journal.pcbi.1000880] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2010] [Accepted: 07/13/2010] [Indexed: 11/27/2022] Open
Abstract
The Polo-Like Kinase 1 (PLK1) acts as a central regulator of mitosis and is over-expressed in a wide range of human tumours where high levels of expression correlate with a poor prognosis. PLK1 comprises two structural elements, a kinase domain and a polo-box domain (PBD). The PBD binds phosphorylated substrates to control substrate phosphorylation by the kinase domain. Although the PBD preferentially binds to phosphopeptides, it has a relatively broad sequence specificity in comparison with other phosphopeptide binding domains. We analysed the molecular determinants of recognition by performing molecular dynamics simulations of the PBD with one of its natural substrates, CDC25c. Predicted binding free energies were calculated using a molecular mechanics, Poisson-Boltzmann surface area approach. We calculated the per-residue contributions to the binding free energy change, showing that the phosphothreonine residue and the mainchain account for the vast majority of the interaction energy. This explains the very broad sequence specificity with respect to other sidechain residues. Finally, we considered the key role of bridging water molecules at the binding interface. We employed inhomogeneous fluid solvation theory to consider the free energy of water molecules on the protein surface with respect to bulk water molecules. Such an analysis highlights binding hotspots created by elimination of water molecules from hydrophobic surfaces. It also predicts that a number of water molecules are stabilized by the presence of the charged phosphate group, and that this will have a significant effect on the binding affinity. Our findings suggest a molecular rationale for the promiscuous binding of the PBD and highlight a role for bridging water molecules at the interface. We expect that this method of analysis will be very useful for probing other protein surfaces to identify binding hotspots for natural binding partners and small molecule inhibitors.
Collapse
Affiliation(s)
- David J Huggins
- TCM Group, Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom.
| | | | | | | | | | | | | | | | | |
Collapse
|
16
|
Abstract
'Everything that living things do can be understood in terms of the jigglings and wigglings of atoms' as Richard Feynman provocatively stated nearly 50 years ago. But how can we 'see' this wiggling and jiggling and understand how it drives biology? Increasingly, computer simulations of biological macromolecules are helping to meet this challenge.
Collapse
Affiliation(s)
- Adrian J Mulholland
- School of Chemistry, Centre for Computational Chemistry, University of Bristol, Bristol, UK.
| |
Collapse
|