1
|
Lee TS, Allen BK, Giese TJ, Guo Z, Li P, Lin C, McGee TD, Pearlman DA, Radak BK, Tao Y, Tsai HC, Xu H, Sherman W, York DM. Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery. J Chem Inf Model 2020; 60:5595-5623. [PMID: 32936637 PMCID: PMC7686026 DOI: 10.1021/acs.jcim.0c00613] [Citation(s) in RCA: 175] [Impact Index Per Article: 43.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Predicting protein-ligand binding affinities and the associated thermodynamics of biomolecular recognition is a primary objective of structure-based drug design. Alchemical free energy simulations offer a highly accurate and computationally efficient route to achieving this goal. While the AMBER molecular dynamics package has successfully been used for alchemical free energy simulations in academic research groups for decades, widespread impact in industrial drug discovery settings has been minimal because of the previous limitations within the AMBER alchemical code, coupled with challenges in system setup and postprocessing workflows. Through a close academia-industry collaboration we have addressed many of the previous limitations with an aim to improve accuracy, efficiency, and robustness of alchemical binding free energy simulations in industrial drug discovery applications. Here, we highlight some of the recent advances in AMBER20 with a focus on alchemical binding free energy (BFE) calculations, which are less computationally intensive than alternative binding free energy methods where full binding/unbinding paths are explored. In addition to scientific and technical advances in AMBER20, we also describe the essential practical aspects associated with running relative alchemical BFE calculations, along with recommendations for best practices, highlighting the importance not only of the alchemical simulation code but also the auxiliary functionalities and expertise required to obtain accurate and reliable results. This work is intended to provide a contemporary overview of the scientific, technical, and practical issues associated with running relative BFE simulations in AMBER20, with a focus on real-world drug discovery applications.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Bryce K. Allen
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Timothy J. Giese
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Zhenyu Guo
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Pengfei Li
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Charles Lin
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - T. Dwight McGee
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - David A. Pearlman
- QSimulate Incorporated, Cambridge, Massachusetts 02139, United States
| | - Brian K. Radak
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Yujun Tao
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Hsu-Chun Tsai
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| | - Huafeng Xu
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Silicon Therapeutics, Boston, Massachusetts 02210, United States
| | - Darrin M. York
- Rutgers, the State University of New Jersey, Laboratory for Biomolecular Simulation Research, and Department of Chemistry and Chemical Biology, United States
| |
Collapse
|
2
|
Kraml J, Hofer F, Kamenik AS, Waibl F, Kahler U, Schauperl M, Liedl KR. Solvation Thermodynamics in Different Solvents: Water-Chloroform Partition Coefficients from Grid Inhomogeneous Solvation Theory. J Chem Inf Model 2020; 60:3843-3853. [PMID: 32639731 PMCID: PMC7460078 DOI: 10.1021/acs.jcim.0c00289] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Indexed: 11/28/2022]
Abstract
Reliable information on partition coefficients plays a key role in drug development, as solubility decisively affects bioavailability. In a physicochemical context, the partition coefficient of a solute between two different solvents can be described as a function of solvation free energies. Hence, substantial scientific efforts have been made toward accurate predictions of solvation free energies in various solvents. The grid inhomogeneous solvation theory (GIST) facilitates the calculation of solvation free energies. In this study, we introduce an extended version of the GIST algorithm, which enables the calculation for chloroform in addition to water. Furthermore, GIST allows localization of enthalpic and entropic contributions. We test our approach by calculating partition coefficients between water and chloroform for a set of eight small molecules. We report a Pearson correlation coefficient of 0.96 between experimentally determined and calculated partition coefficients. The capability to reliably predict partition coefficients between water and chloroform and the possibility to localize their contributions allow the optimization of a compound's partition coefficient. Therefore, we presume that this methodology will be of great benefit for the efficient development of pharmaceuticals.
Collapse
Affiliation(s)
- Johannes Kraml
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Florian Hofer
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Anna S. Kamenik
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Franz Waibl
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Ursula Kahler
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Michael Schauperl
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| | - Klaus R. Liedl
- Institute of General, Inorganic and
Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck
(CMBI), University of Innsbruck, Innrain 80−82, A-6020 Innsbruck, Austria
| |
Collapse
|
3
|
Dasari S, Mallik BS. Biosolvation Nature of Ionic Liquids: Molecular Dynamics Simulation of Methylated Nucleobases in Hydrated 1-Ethyl-3-methylimidazolium Acetate. ACS OMEGA 2018; 3:8344-8354. [PMID: 31458966 PMCID: PMC6644902 DOI: 10.1021/acsomega.8b01231] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/04/2018] [Accepted: 07/18/2018] [Indexed: 06/10/2023]
Abstract
Solvation free energies of methylated nucleobases were calculated in pure and hydrated 1-ethyl-3-methylimidazolium acetate, [Emim][Ac], ionic liquid, and pure water using classical molecular dynamics simulations using multistate Bennett's acceptance ratio method. The calculated solvation free energies in pure water were compared with the previous experimental and theoretical findings and found to be in agreement. We observe that the solvation free energy of methylated nucleobases is more in the pure ionic liquid compared to that in the pure water and on changing the mole fraction of water in the ionic liquid, the solvation free energy decreases gradually. Comparing the Coulombic and van der Waals contribution to the solvation free energy, electrostatic contribution is more compared to that of the latter for all nucleobases. To obtain the atomistic details and explain the solvation mechanism, we calculated radial distribution functions (RDFs), spatial distribution functions (SDFs), and stacking angle distribution of cations to the nucleobases. From RDFs and SDFs, we find that the acetate anions of the ionic liquid are forming strong hydrogen bonds with the amine hydrogen atoms of the nucleobases. These hydrogen bonds contribute to the major part of the Coulombic contribution to the solvation free energy. Stacking of cations to the nucleobases is primarily due to the van der Waals contribution to the solvation free energy.
Collapse
|
4
|
Lara A, Riquelme M, Vöhringer-Martinez E. Partition coefficients of methylated DNA bases obtained from free energy calculations with molecular electron density derived atomic charges. J Comput Chem 2018; 39:1728-1737. [PMID: 29752734 DOI: 10.1002/jcc.25347] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Revised: 03/15/2018] [Accepted: 04/08/2018] [Indexed: 01/25/2023]
Abstract
Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have also been used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study, atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in a vacuum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases, the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform-water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER-99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account, and when the energy needed to polarize the electron density of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogs. Comparison of the two partitioning methods, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS), revealed some deficiencies in the Hirshfeld-I method related to the unstable isolated anionic nitrogen pro-atom used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- A Lara
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Chile
| | - M Riquelme
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Chile
| | - E Vöhringer-Martinez
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Chile
| |
Collapse
|
5
|
Zhang C, Bell D, Harger M, Ren P. Polarizable Multipole-Based Force Field for Aromatic Molecules and Nucleobases. J Chem Theory Comput 2017; 13:666-678. [PMID: 28030769 PMCID: PMC5312700 DOI: 10.1021/acs.jctc.6b00918] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Aromatic molecules with π electrons
are commonly involved
in chemical and biological recognitions. For example, nucleobases
play central roles in DNA/RNA structure and their interactions with
proteins. The delocalization of the π electrons is responsible
for the high polarizability of aromatic molecules. In this work, the
AMOEBA force field has been developed and applied to 5 regular nucleobases
and 12 aromatic molecules. The permanent electrostatic energy is expressed
as atomic multipole interactions between atom pairs, and many-body
polarization is accounted for by mutually induced atomic dipoles.
We have systematically investigated aromatic ring stacking and aromatic-water
interactions for nucleobases and aromatic molecules, as well as base–base
hydrogen-bonding pair interactions, all at various distances and orientations.
van der Waals parameters were determined by comparison to the quantum
mechanical interaction energy of these dimers and fine-tuned using
condensed phase simulation. By comparing to quantum mechanical calculations,
we show that the resulting classical potential is able to accurately
describe molecular polarizability, molecular vibrational frequency,
and dimer interaction energy of these aromatic systems. Condensed
phase properties, including hydration free energy, liquid density,
and heat of vaporization, are also in good overall agreement with
experimental values. The structures of benzene liquid phase and benzene-water
solution were also investigated by simulation and compared with experimental
and PDB structure derived statistical results.
Collapse
Affiliation(s)
- Changsheng Zhang
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - David Bell
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Matthew Harger
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
6
|
Wolf MG, Groenhof G. Evaluating nonpolarizable nucleic acid force fields: A systematic comparison of the nucleobases hydration free energies and chloroform-to-water partition coefficients. J Comput Chem 2012; 33:2225-32. [DOI: 10.1002/jcc.23055] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2012] [Revised: 05/29/2012] [Accepted: 06/04/2012] [Indexed: 02/01/2023]
|
7
|
Gao J. Methods and Applications of Combined Quantum Mechanical and Molecular Mechanical Potentials. REVIEWS IN COMPUTATIONAL CHEMISTRY 2007. [DOI: 10.1002/9780470125847.ch3] [Citation(s) in RCA: 115] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
8
|
Huang N, Kalyanaraman C, Bernacki K, Jacobson MP. Molecular mechanics methods for predicting protein–ligand binding. Phys Chem Chem Phys 2006; 8:5166-77. [PMID: 17203140 DOI: 10.1039/b608269f] [Citation(s) in RCA: 150] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ligand binding affinity prediction is one of the most important applications of computational chemistry. However, accurately ranking compounds with respect to their estimated binding affinities to a biomolecular target remains highly challenging. We provide an overview of recent work using molecular mechanics energy functions to address this challenge. We briefly review methods that use molecular dynamics and Monte Carlo simulations to predict absolute and relative ligand binding free energies, as well as our own work in which we have developed a physics-based scoring method that can be applied to hundreds of thousands of compounds by invoking a number of simplifying approximations. In our previous studies, we have demonstrated that our scoring method is a promising approach for improving the discrimination between ligands that are known to bind and those that are presumed not to, in virtual screening of large compound databases. In new results presented here, we explore several improvements to our computational method including modifying the dielectric constant used for the protein and ligand interiors, and empirically scaling energy terms to compensate for deficiencies in the energy model. Future directions for further improving our physics-based scoring method are also discussed.
Collapse
Affiliation(s)
- Niu Huang
- Department of Pharmaceutical Chemistry, University of California San Francisco, UCSF MC 2240, Genentech Hall, Room N472C, 600 16th St., San Francisco, CA 94158-2517, USA
| | | | | | | |
Collapse
|
9
|
Burda JV, Shukla MK, Leszczynski J. Theoretical model of the aqua-copper [Cu(H2O)5]+cation interactions with guanine. J Mol Model 2005; 11:362-9. [PMID: 15928918 DOI: 10.1007/s00894-005-0269-2] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2004] [Accepted: 01/31/2005] [Indexed: 10/25/2022]
Abstract
Pentaaqua complexes of Cu(I) with guanine were optimized at the DFT B3PW91/6-31G(d) level. For the most stable structures, vibration frequencies and NBO charges were computed followed by energy analyses. The order of individual conformers was very sensitive to the method and basis sets used for the calculation. Several conformers are practically degenerated in energy. The inclusion of an entropy term changes the order of the conformers' stability. Water molecules associated at the N9 position of guanine are favored by the inclusion of the entropy correction. Bonding energies of Cu-O(aqua) interactions were estimated to be about 60 kcal mol(-1) and for Cu-N7 bonding in the range of 75-83 kcal mol(-1). The broad range in Cu-N interaction energies demonstrates the role of induction effects caused by water molecules associated at the various sites of guanine. The charge distribution of the guanine molecule is changed remarkably by the coordination of a Cu(I) cation, which can also change the base-pairing pattern of the guanine.
Collapse
Affiliation(s)
- Jaroslav V Burda
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16, Prague 2, Czech Republic.
| | | | | |
Collapse
|
10
|
Biot C, Buisine E, Rooman M. Free-Energy Calculations of Protein−Ligand Cation−π and Amino−π Interactions: From Vacuum to Proteinlike Environments. J Am Chem Soc 2003; 125:13988-94. [PMID: 14611235 DOI: 10.1021/ja035223e] [Citation(s) in RCA: 72] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
To probe the role of cation-pi and amino-pi interactions in the context of protein-ligand interactions, the stability of 55 X-ray cation/amino-pi motifs involving the Ade moieties of cofactor molecules and Arg, Lys, Asn, or Gln side chains of their host protein was evaluated using quantum chemistry calculations. The conjunction of vacuum interaction energies, vibrational entropy, and solvation contributions led to identify Arg-Ade as the most favorable cation/amino-pi complex in the solvents considered, followed by Asn/Gln-Ade and Lys-Ade: their minimum interaction free energies are approximately equal to -7, -4, and -2 kcal/mol, respectively, in the solvents of dielectric constant similar to that estimated for proteins (i.e., acetone, THF, and CCl(4)). Remarkably, these free-energy values of cation/amino-pi interactions correlate well with their frequency of occurrences in protein-ligand structures, which corroborates our approach in the absence of experimental data.
Collapse
Affiliation(s)
- Christophe Biot
- Ingénierie Biomoléculaire, Université Libre de Bruxelles, CP 165/64, 50 Avenue Roosevelt, B-1050 Brussels, Belgium.
| | | | | |
Collapse
|
11
|
Barsky D, Colvin ME. Guanine−Cytosine Base Pairs in Parallel-Stranded DNA: An ab Initio Study of the Keto−Amino Wobble Pair versus the Enol−Imino Minor Tautomer Pair. J Phys Chem A 2000. [DOI: 10.1021/jp001420d] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Daniel Barsky
- Biology and Biotechnology Research Program, Lawrence Livermore National Lab, L-448, Livermore, California 94550
| | - Michael E. Colvin
- Biology and Biotechnology Research Program, Lawrence Livermore National Lab, L-448, Livermore, California 94550
| |
Collapse
|
12
|
Zhu T, Li J, Hawkins GD, Cramer CJ, Truhlar DG. Density functional solvation model based on CM2 atomic charges. J Chem Phys 1998. [DOI: 10.1063/1.476521] [Citation(s) in RCA: 98] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
Tazaki K, Shimizu K. Molecular Dynamics Simulations in Aqueous Solution: Application to Free Energy Calculation of Oligopeptides. J Phys Chem B 1998. [DOI: 10.1021/jp9803788] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Koichi Tazaki
- Department of Biotechnology, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
| | - Kentaro Shimizu
- Department of Biotechnology, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
| |
Collapse
|
14
|
Ikeguchi M, Shimizu S, Tazaki K, Nakamura S, Shimizu K. Calculation of temperature dependence of free energy caused by potential function changes. Chem Phys Lett 1998. [DOI: 10.1016/s0009-2614(98)00282-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
15
|
Hawkins GD, Cramer CJ, Truhlar DG. Universal Quantum Mechanical Model for Solvation Free Energies Based on Gas-Phase Geometries. J Phys Chem B 1998. [DOI: 10.1021/jp973306+] [Citation(s) in RCA: 149] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
16
|
Eksterowicz JE, Miller JL, Kollman PA. Calculation of Chloroform/Water Partition Coefficients for the N-Methylated Nucleic Acid Bases. J Phys Chem B 1997. [DOI: 10.1021/jp972440r] [Citation(s) in RCA: 36] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John E. Eksterowicz
- Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94143
| | - Jennifer L. Miller
- Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94143
| | - Peter A. Kollman
- Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94143
| |
Collapse
|
17
|
Hawkins GD, Cramer CJ, Truhlar DG. Parametrized Model for Aqueous Free Energies of Solvation Using Geometry-Dependent Atomic Surface Tensions with Implicit Electrostatics. J Phys Chem B 1997. [DOI: 10.1021/jp971035x] [Citation(s) in RCA: 69] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
18
|
Giesen DJ, Chambers CC, Cramer CJ, Truhlar DG. What Controls Partitioning of the Nucleic Acid Bases between Chloroform and Water? J Phys Chem B 1997. [DOI: 10.1021/jp971229c] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David J. Giesen
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| | - Candee C. Chambers
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| | - Christopher J. Cramer
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
19
|
Schön JC. A thermodynamic distance criterion of optimality for the calculation of free energy changes from computer simulations. J Chem Phys 1996. [DOI: 10.1063/1.472836] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
20
|
Affiliation(s)
- Jennifer L. Miller
- Department of Pharmaceutical Chemistry, University of California−San Francisco, San Francisco, California 94143
| | - Peter A. Kollman
- Department of Pharmaceutical Chemistry, University of California−San Francisco, San Francisco, California 94143
| |
Collapse
|
21
|
Poltev VI, Malenkov GG, Gonzalez EJ, Teplukhin AV, Rein R, Shibata M, Miller JH. Modeling DNA hydration: comparison of calculated and experimental hydration properties of nuclic acid bases. J Biomol Struct Dyn 1996; 13:717-26. [PMID: 8906892 DOI: 10.1080/07391102.1996.10508884] [Citation(s) in RCA: 24] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (iii) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9-methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methylguanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.
Collapse
Affiliation(s)
- V I Poltev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Sciences, Pushchino, Moscow Region, Russia
| | | | | | | | | | | | | |
Collapse
|
22
|
|
23
|
Abstract
A combined Monte Carlo quantum mechanical and molecular mechanical (QM/MM) simulation method is used to determine the free energy of hydration and the solvent polarization effect for the nucleotide bases. In the present AM1/TIP3P model, the solute molecule is characterized by valence electrons and effective nucleus cores with Hartree-Fock molecular orbital theory incorporating a solute-solvent interaction Hamiltonian. It is found that polarization energy contributes up to 37%-61% of the total solute-solvent interaction for the systems considered. The computed free energies of hydration are compared with previous theoretical results.
Collapse
Affiliation(s)
- J Gao
- Department of Chemistry, State University of New York at Buffalo 14214
| |
Collapse
|
24
|
Bradley EK, Ng SC, Simon RJ, Spellmeyer DC. Synthesis, molecular modelling, and NMR structure determination of four cyclic peptide antagonists of endothelin. Bioorg Med Chem 1994; 2:279-96. [PMID: 7922139 DOI: 10.1016/s0968-0896(00)82171-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
A combined distance geometry and molecular mechanics/dynamics (MM/MD) protocol was unable to predict the active conformation of the cyclic pentapeptide inhibitor of endothelin-1 receptor, BQ-123, and two analogues. However, the MM/MD method alone is sufficient to predict the solution conformation of a third analogue. In that one case, the combination of proline at residue 3 and an N alpha-methyl substitution at residue 5 provides enough internal constraints to eliminate conformational flexibility seen in the other three analogues. For this constrained analogue, the 50 lowest energy conformations (out of a set of 500 DGEOM-generated, MM/MD minimized conformations) differ by no more than 3.9 kcal/mol. Thirty three of these 50 conformations have backbone atom RMSDs of less than 0.33 A, relative to the lowest energy conformation. The accuracy of this MM/MD model is verified by determining the solution structure of each of the four analogues with 2D NMR techniques. Each of the cyclic pentapeptides has a well defined solution conformation where a proline residue is clearly in a gamma-turn, leaving the remaining residues in a loose beta-turn. All four experimental NMR conformations agree closely with the MM/MD model.
Collapse
|
25
|
Miyamoto S, Kollman PA. Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches. Proteins 1993; 16:226-45. [PMID: 8346190 DOI: 10.1002/prot.340160303] [Citation(s) in RCA: 198] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
We present calculations of the absolute and relative binding free energies of complexation of streptavidin with biotin and its analogs by means of a thermodynamic free energy perturbation method implemented with molecular dynamics. Using the recently solved crystal structure of the streptavidin-biotin complex, biotin was mutated into a dummy molecule as well as thiobiotin and iminobiotin both in the protein and in solution. The calculated absolute binding free energy was dependent on the simulation model used. Encouragingly, the "best models" provided a reasonable semiquantitative reproduction (-20 to -22 kcal/mol) of the experimental free energy (-18.3 kcal/mol). Furthermore, the calculated results give clear insights into the binding nature of the protein-ligand complex, showing that the van der Waals energy dominates the electrostatic and hydrogen bonding energies in the binding of biotin by streptavidin. Specifically, the mutation of biotin into a dummy molecule in solution has a delta G (van der Waals) approximately -4 kcal/mol, due to the cancellation of dispersion and repulsion "cavity" effects. On the other hand, in the protein, a very small free energy price must be paid to create a cavity since one already exists and the mutation of biotin into a dummy molecule has a delta G (van der Waals) approximately 15 kcal/mol. These results are also consistent with the interpretation that the entropy increase to be expected from hydrophobic interactions from desolvation of biotin is counterbalanced by a decrease in entropy accompanying the formation of buried hydrogen bonds, which have been derived from the apparently conflicting experimental data. They provide an alternative interpretation of the reason for the extremely high affinity of the biotin-streptavidin interaction than that recently proposed by Weber et al. (J. Am. Chem. Soc. 114:3197, 1992). In the case of the relative binding free energies, the calculated values of 3.8 +/- 0.6 and 7.2 +/- 0.6 kcal/mol compare well with the experimental values of 3.6 and 6.2 kcal/mol for the perturbation of biotin to thiobiotin and iminobiotin, respectively in the related protein avidin. The calculations indicate that desolvation of the ligand is important in understanding the relative affinity of the ligands with the protein. The above successful simulations suggest that the molecular dynamics/free energy perturbation method is useful for understanding the energetic features affecting the binding between proteins and ligands, since it is generally difficult to determine these factors unambiguously by experiment.(ABSTRACT TRUNCATED AT 400 WORDS)
Collapse
Affiliation(s)
- S Miyamoto
- Department of Pharmaceutical Chemistry, University of California, San Francisco 94143
| | | |
Collapse
|