Reconciling NMR Structures of the HIV-1 Nucleocapsid Protein NCp7 Using Extensive Polarizable Force Field Free-Energy Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: El Khoury L
PROVIDER: S-EPMC7375347 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
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:A balance of van der Waals, electrostatic, and hydrophobic forces drive the folding and packing of protein side chains. Although such interactions between residues are often approximated as being pairwise additive, in reality, higher-order many-body contributions that depend on environment drive hydrophobic collapse and cooperative electrostatics. Beginning from dead-end elimination, we derive the first algorithm, to our knowledge, capable of deterministic global repacking of side chains compatible with many-body energy functions. The approach is applied to seven PCNA x-ray crystallographic data sets with resolutions 2.5-3.8 Å (mean 3.0 Å) using an open-source software. While PDB_REDO models average an Rfree value of 29.5% and MOLPROBITY score of 2.71 Å (77th percentile), dead-end elimination with the polarizable AMOEBA force field lowered Rfree by 2.8-26.7% and improved mean MOLPROBITY score to atomic resolution at 1.25 Å (100th percentile). For structural biology applications that depend on side-chain repacking, including x-ray refinement, homology modeling, and protein design, the accuracy limitations of pairwise additivity can now be eliminated via polarizable or quantum mechanical potentials.
| S-EPMC4547145 | biostudies-literature
Project description:In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.
| S-EPMC6017617 | biostudies-literature
Project description:The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling.
| S-EPMC3948463 | biostudies-literature
Project description:Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.
| S-EPMC5595662 | biostudies-other
Project description:We have recently introduced a quantum mechanical polarizable force field (QMPFF) fitted solely to high-level quantum mechanical data for simulations of biomolecular systems. Here, we present an improved form of the force field, QMPFF2, and apply it to simulations of liquid water. The results of the simulations show excellent agreement with a variety of experimental thermodynamic and structural data, as good or better than that provided by specialized water potentials. In particular, QMPFF2 is the only ab initio force field to accurately reproduce the anomalous temperature dependence of water density to our knowledge. The ability of the same force field to successfully simulate the properties of both organic molecules and water suggests it will be useful for simulations of proteins and protein-ligand interactions in the aqueous environment.
| S-EPMC1466545 | 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 previously introduced polarizable simulations with second-order interaction model (POSSIM) force field has been extended to include parameters for small molecules serving as models for peptide and protein side-chains. Parameters have been fitted to permit reproducing many-body energies, gas-phase dimerization energies, and geometries and liquid-phase heats of vaporization and densities. Quantum mechanical and experimental data have been used as the target for the fitting. The POSSIM framework combines accuracy of a polarizable force field and computational efficiency of the second-order approximation of the full-scale induced point dipole polarization formalism. The resulting parameters can be used for simulations of the parameterized molecules themselves or their analogues. In addition to this, these force field parameters are currently being used in further development of the POSSIM fast polarizable force field for proteins.
| S-EPMC3633718 | biostudies-literature
Project description:Presented is the implementation of the Drude force field in the open-source OpenMM simulation package allowing for access to graphical processing unit (GPU) hardware. In the Drude model, electronic degrees of freedom are represented by negatively charged particles attached to their parent atoms via harmonic springs, such that extra computational overhead comes from these additional particles and virtual sites representing lone pairs on electronegative atoms, as well as the associated thermostat and integration algorithms. This leads to an approximately fourfold increase in computational demand over additive force fields. However, by making the Drude model accessible to consumer-grade desktop GPU hardware it will be possible to perform simulations of one microsecond or more in less than a month, indicating that the barrier to employ polarizable models has largely been removed such that polarizable simulations with the classical Drude model are readily accessible and practical.
| S-EPMC6031474 | biostudies-literature
Project description:Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models toward more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date, it has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with AMOEBA for demanding problems such as relative conformational energies of the alanine tetrapeptide and isomers of water sulfate complexes. AMOEBA is shown to be especially successful on protein-ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are critical.
| S-EPMC2918242 | biostudies-literature