1
|
Bodnarchuk MS, Viner R, Michel J, Essex JW. Strategies to calculate water binding free energies in protein-ligand complexes. J Chem Inf Model 2014; 54:1623-33. [PMID: 24684745 DOI: 10.1021/ci400674k] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Water molecules are commonplace in protein binding pockets, where they can typically form a complex between the protein and a ligand or become displaced upon ligand binding. As a result, it is often of great interest to establish both the binding free energy and location of such molecules. Several approaches to predicting the location and affinity of water molecules to proteins have been proposed and utilized in the literature, although it is often unclear which method should be used under what circumstances. We report here a comparison between three such methodologies, Just Add Water Molecules (JAWS), Grand Canonical Monte Carlo (GCMC), and double-decoupling, in the hope of understanding the advantages and limitations of each method when applied to enclosed binding sites. As a result, we have adapted the JAWS scoring procedure, allowing the binding free energies of strongly bound water molecules to be calculated to a high degree of accuracy, requiring significantly less computational effort than more rigorous approaches. The combination of JAWS and GCMC offers a route to a rapid scheme capable of both locating and scoring water molecules for rational drug design.
Collapse
Affiliation(s)
- Michael S Bodnarchuk
- School of Chemistry, University of Southampton , Highfield, Southampton, SO17 1BJ, U.K
| | | | | | | |
Collapse
|
2
|
Lin YL, Aleksandrov A, Simonson T, Roux B. An Overview of Electrostatic Free Energy Computations for Solutions and Proteins. J Chem Theory Comput 2014; 10:2690-709. [PMID: 26586504 DOI: 10.1021/ct500195p] [Citation(s) in RCA: 104] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Free energy simulations for electrostatic and charging processes in complex molecular systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calculate the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liquid-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no physical boundary to the system. An important consequence is that the reference value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the volume fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.
Collapse
Affiliation(s)
- Yen-Lin Lin
- Department of Biochemistry and Molecular Biology, University of Chicago , Chicago, Illinois 60637, United States
| | - Alexey Aleksandrov
- Department of Biology, Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique , 91128 Palaiseau, France
| | - Thomas Simonson
- Department of Biology, Laboratoire de Biochimie (CNRS UMR7654), Ecole Polytechnique , 91128 Palaiseau, France
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago , Chicago, Illinois 60637, United States.,Biosciences Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| |
Collapse
|
3
|
Genheden S, Cabedo Martinez AI, Criddle MP, Essex JW. Extensive all-atom Monte Carlo sampling and QM/MM corrections in the SAMPL4 hydration free energy challenge. J Comput Aided Mol Des 2014; 28:187-200. [PMID: 24488307 DOI: 10.1007/s10822-014-9717-3] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2013] [Accepted: 01/21/2014] [Indexed: 12/14/2022]
Abstract
We present our predictions for the SAMPL4 hydration free energy challenge. Extensive all-atom Monte Carlo simulations were employed to sample the compounds in explicit solvent. While the focus of our study was to demonstrate well-converged and reproducible free energies, we attempted to address the deficiencies in the general Amber force field force field with a simple QM/MM correction. We show that by using multiple independent simulations, including different starting configurations, and enhanced sampling with parallel tempering, we can obtain well converged hydration free energies. Additional analysis using dihedral angle distributions, torsion-root mean square deviation plots and thermodynamic cycles support this assertion. We obtain a mean absolute deviation of 1.7 kcal mol(-1) and a Kendall's τ of 0.65 compared with experiment.
Collapse
Affiliation(s)
- Samuel Genheden
- School of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, UK,
| | | | | | | |
Collapse
|
4
|
Cisneros GA, Karttunen M, Ren P, Sagui C. Classical electrostatics for biomolecular simulations. Chem Rev 2014; 114:779-814. [PMID: 23981057 PMCID: PMC3947274 DOI: 10.1021/cr300461d] [Citation(s) in RCA: 192] [Impact Index Per Article: 19.2] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
5
|
Reif MM, Kräutler V, Kastenholz MA, Daura X, Hünenberger PH. Molecular dynamics simulations of a reversibly folding beta-heptapeptide in methanol: influence of the treatment of long-range electrostatic interactions. J Phys Chem B 2009; 113:3112-28. [PMID: 19228001 DOI: 10.1021/jp807421a] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Eight 100-ns molecular dynamics simulations of a beta-heptapeptide in methanol at 340 K (within cubic periodic computational boxes of about 6-nm edge) are reported and compared. These simulations were performed with three different charge-state combinations at the peptide termini, one of them with or without a neutralizing chloride counterion, and using either the lattice-sum (LS) or reaction-field (RF) scheme to handle electrostatic interactions. The choice of the electrostatic scheme has essentially no influence on the folding-unfolding equilibrium when the peptide termini are uncharged and only a small influence when the peptide is positively charged at its N-terminus (with or without inclusion of a neutralizing chloride counterion). However, when the peptide is zwitterionic, the LS scheme leads to preferential sampling of the high-dipole folded helical state, whereas the RF scheme leads to preferential sampling of a low-dipole unfolded salt-bridged state. A continuum electrostatics analysis based on the sampled configurations (zwitterionic case) suggests that the LS scheme stabilizes the helical state through artificial periodicity, but that the magnitude of this perturbation is essentially negligible (compared to the thermal energy) for the large box size and relatively polar solvent considered. The results thus provide clear evidence (continuum electrostatics analysis) for the absence of LS artifacts and some indications (still not definitive because of the limited sampling of the folding-unfolding transition) for the presence of RF artifacts in this specific system.
Collapse
Affiliation(s)
- Maria M Reif
- Laboratory of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | | | | | | | | |
Collapse
|
6
|
Brooks B, Brooks C, MacKerell A, Nilsson L, Petrella R, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner A, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor R, Post C, Pu J, Schaefer M, Tidor B, Venable RM, Woodcock HL, Wu X, Yang W, York D, Karplus M. CHARMM: the biomolecular simulation program. J Comput Chem 2009; 30:1545-614. [PMID: 19444816 PMCID: PMC2810661 DOI: 10.1002/jcc.21287] [Citation(s) in RCA: 6075] [Impact Index Per Article: 405.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.
Collapse
Affiliation(s)
- B.R. Brooks
- Laboratory of Computational Biology, National Heart, Lung, and
Blood Institute, National Institutes of Health, Bethesda, MD 20892
| | - C.L. Brooks
- Departments of Chemistry & Biophysics, University of
Michigan, Ann Arbor, MI 48109
| | - A.D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy,
University of Maryland, Baltimore, MD, 21201
| | - L. Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition,
SE-141 57, Huddinge, Sweden
| | - R.J. Petrella
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Department of Medicine, Harvard Medical School, Boston, MA
02115
| | - B. Roux
- Department of Biochemistry and Molecular Biology, University of
Chicago, Gordon Center for Integrative Science, Chicago, IL 60637
| | - Y. Won
- Department of Chemistry, Hanyang University, Seoul
133–792 Korea
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - M. Karplus
- Department of Chemistry and Chemical Biology, Harvard University,
Cambridge, MA 02138
- Laboratoire de Chimie Biophysique, ISIS, Université de
Strasbourg, 67000 Strasbourg France
| |
Collapse
|
7
|
Li H, Fajer M, Yang W. Simulated scaling method for localized enhanced sampling and simultaneous “alchemical” free energy simulations: A general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations. J Chem Phys 2007; 126:024106. [PMID: 17228942 DOI: 10.1063/1.2424700] [Citation(s) in RCA: 73] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A potential scaling version of simulated tempering is presented to efficiently sample configuration space in a localized region. The present "simulated scaling" method is developed with a Wang-Landau type of updating scheme in order to quickly flatten the distributions in the scaling parameter lambdam space. This proposal is meaningful for a broad range of biophysical problems, in which localized sampling is required. Besides its superior capability and robustness in localized conformational sampling, this simulated scaling method can also naturally lead to efficient "alchemical" free energy predictions when dual-topology alchemical hybrid potential is applied; thereby simultaneously, both of the chemically and conformationally distinct portions of two end point chemical states can be efficiently sampled. As demonstrated in this work, the present method is also feasible for the quantum mechanical and quantum mechanical/molecular mechanical simulations.
Collapse
Affiliation(s)
- Hongzhi Li
- School of Computational Science, Florida State University, Tallahassee, Florida 32306, USA
| | | | | |
Collapse
|
8
|
Ma A, Hu J, Karplus M, Dinner AR. Implications of Alternative Substrate Binding Modes for Catalysis by Uracil-DNA Glycosylase: An Apparent Discrepancy Resolved. Biochemistry 2006; 45:13687-96. [PMID: 17105188 DOI: 10.1021/bi061061y] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A theoretical study showed that the base excision repair enzyme uracil-DNA glycosylase (UDG) exploits electrostatic interactions with backbone phosphate groups in the substrate for catalysis. Although experiments performed to test the calculated results confirmed the predicted importance of the -2, -1, and +1 phosphate groups, there was an apparent disagreement with regard to the +2 phosphate group. The calculations indicated that it made an important contribution, while experimentally, the effect of its deletion or neutralization was small. The +2 phosphate group interacts directly with an active site histidine (H148 in humans) in the crystal structure of UDG in complex with double-stranded (ds) DNA. We previously calculated that H148 has a strong anticatalytic effect due to its protonation, and here we use alchemical free energy simulations to estimate its site-specific pKa. We find that it is positively charged over the entire experimental pH range (4-10), so its deprotonation cannot compensate for deletion or neutralization of the +2 phosphate group. The free energy simulations are facilitated by an efficient charge-scaling procedure that allows quantitative correction for the implicit treatment of solvent far from the active site; improvements are made to that method to account carefully for differences in the truncation of electrostatic interactions in the contributing molecular-mechanical and continuum-electrostatic approaches. Additional simulations are used to demonstrate that the +2 phosphate group is fully solvent exposed in complexes with single-stranded DNA substrates like those used in the experiments. In contrast, it is well-structured and buried in the dsDNA complex used in the original simulations. Differences in solvent shielding thus account for the apparent lack of an effect observed experimentally upon neutralization or deletion of this group.
Collapse
Affiliation(s)
- Ao Ma
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | | | | | | |
Collapse
|
9
|
Laino T, Mohamed F, Laio A, Parrinello M. An Efficient Linear-Scaling Electrostatic Coupling for Treating Periodic Boundary Conditions in QM/MM Simulations. J Chem Theory Comput 2006; 2:1370-8. [DOI: 10.1021/ct6001169] [Citation(s) in RCA: 133] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Teodoro Laino
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56125 Pisa, Italy, and Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Fawzi Mohamed
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56125 Pisa, Italy, and Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Alessandro Laio
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56125 Pisa, Italy, and Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56125 Pisa, Italy, and Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| |
Collapse
|
10
|
Kastenholz MA, Hünenberger PH. Computation of methodology-independent ionic solvation free energies from molecular simulations. II. The hydration free energy of the sodium cation. J Chem Phys 2006; 124:224501. [PMID: 16784292 DOI: 10.1063/1.2201698] [Citation(s) in RCA: 161] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The raw ionic solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions (finite or periodic system, system shape, and size) and treatment of electrostatic interactions (Coulombic, lattice sum, or cutoff based) used during these simulations. In the present article, it is shown that correction terms can be derived for the effect of (A) an incorrect solvent polarization around the ion due to the use of an approximate (not strictly Coulombic) electrostatic scheme; (B) the finite size or artificial periodicity of the simulated system; (C) an improper summation scheme to evaluate the potential at the ion site and the possible presence of a liquid-vacuum interface in the simulated system. Taking the hydration free energy of the sodium cation as a test case, it is shown that the raw solvation free energies obtained using seven different types of boundary conditions and electrostatic schemes commonly used in explicit-solvent simulations (for a total of 72 simulations differing in the corresponding simulation parameters) can be corrected so as to obtain a consistent value for this quantity.
Collapse
Affiliation(s)
- Mika A Kastenholz
- Laboratorium für Physikalische Chemie, ETH Zentrum, CH-8093 Zürich, Switzerland
| | | |
Collapse
|
11
|
Paliwal A, Asthagiri D, Pratt LR, Ashbaugh HS, Paulaitis ME. An analysis of molecular packing and chemical association in liquid water using quasichemical theory. J Chem Phys 2006; 124:224502. [PMID: 16784293 DOI: 10.1063/1.2202350] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We calculate the hydration free energy of liquid TIP3P water at 298 K and 1 bar using a quasi-chemical theory framework in which interactions between a distinguished water molecule and the surrounding water molecules are partitioned into chemical associations with proximal (inner-shell) waters and classical electrostatic-dispersion interactions with the remaining (outer-shell) waters. The calculated free energy is found to be independent of this partitioning, as expected, and in excellent agreement with values derived from the literature. An analysis of the spatial distribution of inner-shell water molecules as a function of the inner-shell volume reveals that water molecules are preferentially excluded from the interior of large volumes as the occupancy number decreases. The driving force for water exclusion is formulated in terms of a free energy for rearranging inner-shell water molecules under the influence of the field exerted by outer-shell waters in order to accommodate one water molecule at the center. The results indicate a balance between chemical association and molecular packing in liquid water that becomes increasingly important as the inner-shell volume grows in size.
Collapse
Affiliation(s)
- A Paliwal
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | | | | | | | | |
Collapse
|
12
|
Yang W, Bitetti-Putzer R, Karplus M. Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. J Chem Phys 2006; 120:2618-28. [PMID: 15268405 DOI: 10.1063/1.1638996] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method is proposed for improving the accuracy and efficiency of free energy simulations. The essential idea is that the convergence of the relevant measure (e.g., the free energy derivative in thermodynamic integration) is monitored in the reverse direction starting from the last frame of the trajectory, instead of the usual approach, which begins with the first frame and goes in the forward direction. This simple change in the use of the simulation data makes it straightforward to eliminate the contamination of the averages by contributions from the equilibrating region. A statistical criterion is introduced for distinguishing the equilibrated (production) region from the equilibrating region. The proposed method, called reverse cumulative averaging, is illustrated by its application to the well-studied case of the alchemical free energy simulation of ethane to methanol.
Collapse
Affiliation(s)
- Wei Yang
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | | | | |
Collapse
|
13
|
Schilling B, Brickmann J, Kast SM. Hybrid integral equation/simulation model for enhancing free energy computations. Phys Chem Chem Phys 2006; 8:1086-95. [PMID: 16633590 DOI: 10.1039/b514185k] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Integral equation theory is used for extrapolating free energy data from molecular simulations of a reference state with respect to a modification of the interaction potential. The methodology is applied to the correction of artefacts arising from potential shifting and truncation. Corrective contributions for the hydration free energy with respect to the full potential are analysed for the case that both the solute-solvent as well as the solvent-solvent potentials are truncated and modified by a shifted-force term, reaching beyond the range of the dielectric continuum approximation and simple long-range correction expressions. The model systems argon in water and pure water are used as examples for apolar and polar solutes, revealing significant correction contributions even for the short-ranged dispersive interactions and the magnitude of solute-solvent and solvent-solvent components. In comparison with simulation-based extrapolation techniques the integral equation method is shown to be capable of quantitatively predicting truncation artefacts at negligible computational overhead.
Collapse
Affiliation(s)
- Bernd Schilling
- Physikalische Chemie, Technische Universität Darmstadt, Petersenstrasse 20, 64287 Darmstadt, Germany
| | | | | |
Collapse
|
14
|
Heinz TN, Hünenberger PH. Combining the lattice-sum and reaction-field approaches for evaluating long-range electrostatic interactions in molecular simulations. J Chem Phys 2005; 123:34107. [PMID: 16080730 DOI: 10.1063/1.1955525] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A new scheme, the lattice-sum-emulated reaction-field (LSERF) method, is presented that combines the lattice-sum (LS) and reaction-field (RF) approaches for evaluating electrostatic interactions in molecular simulations. More precisely, the LSERF scheme emulates a RF calculation (based on an atomic cutoff) via the LS machinery. This is achieved by changing the form of the electrostatic interactions in a standard LS calculation (Coulombic) to the form corresponding to RF electrostatics (Coulombic plus quadratic reaction-field correction term, truncated at the cutoff distance). It is shown (both analytically and numerically) that in the limit of infinite reciprocal-space accuracy, (i) the LSERF scheme with a finite reaction-field cutoff and a given reaction-field permittivity is identical to the RF scheme with the same parameters (and an atomic cutoff), and (ii) the LSERF scheme is identical to the LS scheme in the limit of an infinite reaction-field cutoff, irrespective of the reaction-field permittivity. This new scheme offers two key advantages: (i) from a conceptual point of view, it shows that there is a continuity between the RF and LS schemes and unifies them into a common framework; (ii) from a practical point of view, it allows us to perform RF calculations with arbitrarily large reaction-field cutoff distances for the same computational costs as a corresponding LS calculation. The optimal choice for the cutoff will be the one that achieves the best compromise between artifacts arising from the dielectric heterogeneity of the system (short cutoff) and its artificial periodicity (long cutoff). The implementation of the LSERF method is extremely easy, requiring only very limited modifications of any standard LS code. For practical applications to biomolecular systems, the use of the LSERF scheme with large reaction-field cutoff distances is expected to represent a significant improvement over the current RF simulations involving comparatively much shorter cutoffs.
Collapse
Affiliation(s)
- Tim N Heinz
- Laboratorium für Physikalische Chemie, Eidgenössìsche Technische Hochschule Hönggerberg, Zürich, Switzerland
| | | |
Collapse
|
15
|
Dehez F, Martins-Costa MTC, Rinaldi D, Millot C. Long-range electrostatic interactions in hybrid quantum and molecular mechanical dynamics using a lattice summation approach. J Chem Phys 2005; 122:234503. [PMID: 16008458 DOI: 10.1063/1.1931667] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A method based on a lattice summation technique for treating long-range electrostatic interactions in hybrid quantum mechanics/molecular mechanics simulations is presented in this article. The quantum subsystem is studied at the semiempirical level, whereas the solvent is described by a two-body potential of molecular mechanics. Molecular dynamics simulations of a (quantum) chloride ion in (classical) water have been performed to test this technique. It is observed that the application of the lattice summations to solvent-solvent interactions as well as on solute-solvent ones has a significant effect on solvation energy and diffusion coefficient. Moreover, two schemes for the computation of the long-range contribution to the electrostatic interaction energy are investigated. The first one replaces the exact charge distribution of the quantum solute by a Mulliken charge distribution. The long-range electrostatic interactions are then calculated for this charge distribution that interacts with the solvent molecule charges. The second one is more accurate and involves a modified Fock operator containing long-range electron-charge interactions. It is shown here that both schemes lead to similar results, the method using Mulliken charges for the evaluation of long-range interactions being, however, much more computationally efficient.
Collapse
Affiliation(s)
- François Dehez
- Unité Mixte de Recherche, Centre National de la Recherche Scientifique-Université Henri Poincaré (CNRS-UHP) 7565, Institut Nancéien de Chimie Moléculaire, Université Henri Poincaré-Nancy I, BP 239, Vandoeuvre-lès-Nancy 54506, France
| | | | | | | |
Collapse
|
16
|
D'Abramo M, D'Alessandro M, Amadei A. On the use of the quasi-Gaussian entropy theory in the study of simulated dilute solutions. J Chem Phys 2004; 120:5226-34. [PMID: 15267394 DOI: 10.1063/1.1647530] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
In a recent paper [M. D'Alessandro, M. D'Abramo, G. Brancato, A. Di Nola, and A. Amadei, J. Phys. Chem. B 106, 11843 (2002)] we showed how to combine molecular dynamics simulations with the quasi-Gaussian entropy theory, in order to model the statistical mechanics and thermodynamics of ionic (water) solutions. In this paper we extend the method to treat nonspherical solutes, describe more thoroughly its theoretical basis and apply it to a set of more complex solute molecules in water (i.e., water, methane, ethane, methanol, and ethanol). Results show that this approach can really provide an excellent theoretical description of solute-solvent systems over a wide range of temperatures.
Collapse
Affiliation(s)
- M D'Abramo
- Dipartimento di Chimica, Universita di Roma La Sapienza, Italy
| | | | | |
Collapse
|
17
|
Zhang Z, Duan Z. Lithium chloride ionic association in dilute aqueous solution: a constrained molecular dynamics study. Chem Phys 2004. [DOI: 10.1016/j.chemphys.2003.10.030] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
18
|
Boresch S, Steinhauser O. The dielectric self-consistent field method. I. Highways, byways, and illustrative results. J Chem Phys 2001. [DOI: 10.1063/1.1421359] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
19
|
Boresch S, Steinhauser O. The dielectric self-consistent field method. II. Application to the study of finite range effects. J Chem Phys 2001. [DOI: 10.1063/1.1421360] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|