Water properties from first principles: simulations by a general-purpose quantum mechanical polarizable force field.
Ontology highlight
ABSTRACT:
SUBMITTER: Donchev AG
PROVIDER: S-EPMC1466545 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We introduce a quantum mechanical polarizable force field (QMPFF) fitted solely to QM data at the MP2/aTZ(-hp) level. Atomic charge density is modeled by point-charge nuclei and floating exponentially shaped electron clouds. The functional form of interaction energy parallels quantum mechanics by including electrostatic, exchange, induction, and dispersion terms. Separate fitting of each term to the counterpart calculated from high-quality QM data ensures high transferability of QMPFF parameters to different molecular environments, as well as accurate fit to a broad range of experimental data in both gas and liquid phases. QMPFF, which is much more efficient than ab initio QM, is optimized for the accurate simulation of biomolecular systems and the design of drugs.
| S-EPMC1142387 | biostudies-literature
Project description:Conspectus There is need in the molecular simulation community to develop new quantum mechanical (QM) methods that can be routinely applied to the simulation of large molecular systems in complex, heterogeneous condensed phase environments. Although conventional methods, such as the hybrid quantum mechanical/molecular mechanical (QM/MM) method, are adequate for many problems, there remain other applications that demand a fully quantum mechanical approach. QM methods are generally required in applications that involve changes in electronic structure, such as when chemical bond formation or cleavage occurs, when molecules respond to one another through polarization or charge transfer, or when matter interacts with electromagnetic fields. A full QM treatment, rather than QM/MM, is necessary when these features present themselves over a wide spatial range that, in some cases, may span the entire system. Specific examples include the study of catalytic events that involve delocalized changes in chemical bonds, charge transfer, or extensive polarization of the macromolecular environment; drug discovery applications, where the wide range of nonstandard residues and protonation states are challenging to model with purely empirical MM force fields; and the interpretation of spectroscopic observables. Unfortunately, the enormous computational cost of conventional QM methods limit their practical application to small systems. Linear-scaling electronic structure methods (LSQMs) make possible the calculation of large systems but are still too computationally intensive to be applied with the degree of configurational sampling often required to make meaningful comparison with experiment. In this work, we present advances in the development of a quantum mechanical force field (QMFF) suitable for application to biological macromolecules and condensed phase simulations. QMFFs leverage the benefits provided by the LSQM and QM/MM approaches to produce a fully QM method that is able to simultaneously achieve very high accuracy and efficiency. The efficiency of the QMFF is made possible by partitioning the system into fragments and self-consistently solving for the fragment-localized molecular orbitals in the presence of the other fragment's electron densities. Unlike a LSQM, the QMFF introduces empirical parameters that are tuned to obtain very accurate intermolecular forces. The speed and accuracy of our QMFF is demonstrated through a series of examples ranging from small molecule clusters to condensed phase simulation, and applications to drug docking and protein-protein interactions. In these examples, comparisons are made to conventional molecular mechanical models, semiempirical methods, ab initio Hamiltonians, and a hybrid QM/MM method. The comparisons demonstrate the superior accuracy of our QMFF relative to the other models; nonetheless, we stress that the overarching role of QMFFs is not to supplant these established computational methods for problems where their use is appropriate. The role of QMFFs within the toolbox of multiscale modeling methods is to extend the range of applications to include problems that demand a fully quantum mechanical treatment of a large system with extensive configurational sampling.
| S-EPMC4165466 | biostudies-other
Project description:In this work, we report two polarizable molecular mechanics (polMM) force field models for estimating the polarization energy in hybrid quantum mechanical molecular mechanical (QM/MM) calculations. These two models, named the potential of atomic charges (PAC) and potential of atomic dipoles (PAD), are formulated from the ab initio quantum mechanical (QM) response kernels for the prediction of the QM density response to an external molecular mechanical (MM) environment (as described by external point charges). The PAC model is similar to fluctuating charge (FQ) models because the energy depends on external electrostatic potential values at QM atomic sites; the PAD energy depends on external electrostatic field values at QM atomic sites, resembling induced dipole (ID) models. To demonstrate their uses, we apply the PAC and PAD models to 12 small molecules, which are solvated by TIP3P water. The PAC model reproduces the QM/MM polarization energy with a R2 value of 0.71 for aniline (in 10,000 TIP3P water configurations) and 0.87 or higher for other 11 solute molecules, while the PAD model has a much better performance with R2 values of 0.98 or higher. The PAC model reproduces reference QM/MM hydration free energies for 12 solute molecules with a RMSD of 0.59 kcal/mol. The PAD model is even more accurate, with a much smaller RMSD of 0.12 kcal/mol, with respect to the reference. This suggests that polarization effects, including both local charge distortion and intramolecular charge transfer, can be well captured by induced dipole type models with proper parametrization.
| S-EPMC5315702 | biostudies-literature
Project description:We present an ab initio polarizable representation of classical molecular mechanics (MM) atoms by employing an angular momentum-based expansion scheme of the point charges into partial wave orbitals. The charge density represented by these orbitals can be fully polarized, and for hybrid quantum-mechanical-molecular-mechanical (QM/MM) calculations, mutual polarization within the QM/MM Hamiltonian can be obtained. We present the mathematical formulation and the analytical expressions for the energy and forces pertaining to the method. We further develop a variational scheme to appropriately determine the expansion coefficients and then validate the method by considering polarizations of ions by the QM system employing the hybrid GROMACS-CPMD QM/MM program. Finally, we present a simpler prescription for adding isotropic polarizability to MM atoms in a QM/MM simulation. Employing this simpler scheme, we present QM/MM energy minimization results for the classic case of a water dimer and a hydrogen sulfide dimer. Also, we present single-point QM/MM results with and without the polarization to study the change in the ionization potential of tetrahydrobiopterin (BH(4)) in water and the change in the interaction energy of solvated BH(4) (described by MM) with the P(450) heme described by QM. The model can be employed for the development of an extensive classical polarizable force-field.
| S-EPMC2736615 | biostudies-literature
Project description:Although a great number of computational models of water are available today, the majority of current biological simulations are done with simple models, such as TIP3P and SPC, developed almost thirty years ago and only slightly modified since then. The reason is that the non-polarizable force fields that are mostly used to describe proteins and other biological molecules are incompatible with more sophisticated modern polarizable models of water. The issue is electronic polarizability: in liquid state, in protein, and in vacuum the water molecule is polarized differently, and therefore has different properties; thus the only way to describe all these different media with the same model is to use a polarizable water model. However, to be compatible with the force field of the rest of the system, e.g. a protein, the latter should be polarizable as well. Here we describe a novel model of water that is in effect polarizable, and yet compatible with the standard non-polarizable force fields such as AMBER, CHARMM, GROMOS, OPLS, etc. Thus the model resolves the outstanding problem of incompatibility.
| S-EPMC4285689 | biostudies-literature
Project description:A quantum mechanical force field (QMFF) for water is described. Unlike traditional approaches that use quantum mechanical results and experimental data to parameterize empirical potential energy functions, the present QMFF uses a quantum mechanical framework to represent intramolecular and intermolecular interactions in an entire condensed-phase system. In particular, the internal energy terms used in molecular mechanics are replaced by a quantum mechanical formalism that naturally includes electronic polarization due to intermolecular interactions and its effects on the force constants of the intramolecular force field. As a quantum mechanical force field, both intermolecular interactions and the Hamiltonian describing the individual molecular fragments can be parameterized to strive for accuracy and computational efficiency. In this work, we introduce a polarizable molecular orbital model Hamiltonian for water and for oxygen- and hydrogen-containing compounds, whereas the electrostatic potential responsible for intermolecular interactions in the liquid and in solution is modeled by a three-point charge representation that realistically reproduces the total molecular dipole moment and the local hybridization contributions. The present QMFF for water, which is called the XP3P (explicit polarization with three-point-charge potential) model, is suitable for modeling both gas-phase clusters and liquid water. The paper demonstrates the performance of the XP3P model for water and proton clusters and the properties of the pure liquid from about 900 × 10(6) self-consistent-field calculations on a periodic system consisting of 267 water molecules. The unusual dipole derivative behavior of water, which is incorrectly modeled in molecular mechanics, is naturally reproduced as a result of an electronic structural treatment of chemical bonding by XP3P. We anticipate that the XP3P model will be useful for studying proton transport in solution and solid phases as well as across biological ion channels through membranes.
| S-EPMC3747793 | biostudies-other
Project description:Accurately predicting protein behavior across diverse pH environments remains a significant challenge in biomolecular simulations. Existing constant-pH molecular dynamics (CpHMD) algorithms are limited to fixed-charge force fields, hindering their application to biomolecular systems described by permanent atomic multipoles or induced dipoles. This work overcomes these limitations by introducing the first polarizable CpHMD algorithm in the context of the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. Additionally, our implementation in the open-source Force Field X (FFX) software has the unique ability to handle titration state changes for crystalline systems including flexible support for all 230 space groups. The evaluation of constant-pH molecular dynamics (CpHMD) with the AMOEBA force field was performed on 11 crystalline peptide systems that span the titrating amino acids (Asp, Glu, His, Lys, and Cys). Titration states were correctly predicted for 15 out of the 16 amino acids present in the 11 systems, including for the coordination of Zn2+ by cysteines. The lone exception was for a HIS-ALA peptide where CpHMD predicted both neutral histidine tautomers to be equally populated, whereas the experimental model did not consider multiple conformers and diffraction data are unavailable for rerefinement. This work demonstrates the promise polarizable CpHMD simulations for pKa predictions, the study of biochemical mechanisms such as the catalytic triad of proteases, and for improved protein-ligand binding affinity accuracy in the context of pharmaceutical lead optimization.
| S-EPMC11008096 | biostudies-literature
Project description:Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models.
| S-EPMC2883601 | biostudies-literature
Project description:We present here the first application of the quantum chemical topology force field FFLUX to condensed matter simulations. FFLUX offers many-body potential energy surfaces learnt exclusively from ab initio data using Gaussian process regression. FFLUX also includes high-rank, polarizable multipole moments (up to quadrupole moments in this work) that are learnt from the same ab initio calculations as the potential energy surfaces. Many-body effects (where a body is an atom) and polarization are captured by the machine learning models. The choice to use machine learning in this way allows the force field's representation of reality to be improved (e.g., by including higher order many-body effects) with little detriment to the computational scaling of the code. In this manner, FFLUX is inherently future-proof. The "plug and play" nature of the machine learning models also ensures that FFLUX can be applied to any system of interest, not just liquid water. In this work we study liquid water across a range of temperatures and compare the predicted bulk properties to experiment as well as other state-of-the-art force fields AMOEBA(+CF), HIPPO, MB-Pol and SIBFA21. We find that FFLUX finds a place amongst these.
| S-EPMC9476653 | biostudies-literature
Project description:Recently synthetized iron complexes have achieved long-lived excited states and stabilities which are comparable, or even superior, to their ruthenium analogues, thus representing an eco-friendly and cheaper alternative to those materials based on rare metals. Most of computational tools which could help unravel the origin of this large efficiency rely on ab-initio methods which are not able, however, to capture the nanosecond time scale underlying these photophysical processes and the influence of their realistic environment. Therefore, it exists an urgent need of developing new low-cost, but still accurate enough, computational methodologies capable to deal with the steady-state and transient spectroscopy of transition metal complexes in solution. Following this idea, here we focus on the comparison between general-purpose transferable force-fields (FFs), directly available from existing databases, and specific quantum mechanical derived FFs (QMD-FFs), obtained in this work through the Joyce procedure. We have chosen a recently reported FeIII complex with nanosecond excited-state lifetime as a representative case. Our molecular dynamics (MD) simulations demonstrated that the QMD-FF nicely reproduces the structure and the dynamics of the complex and its chemical environment within the same precision as higher cost QM methods, whereas general-purpose FFs failed in this purpose. Although in this particular case the chemical environment plays a minor role on the photo physics of this system, these results highlight the potential of QMD-FFs to rationalize photophysical phenomena provided an accurate QM method to derive its parameters is chosen.
| S-EPMC7411876 | biostudies-literature