Browse
Submit Data
Databases
API
Help

Dataset Information

35 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

GEM*: A Molecular Electronic Density-Based Force Field for Molecular Dynamics Simulations.

ABSTRACT: GEM*, a force field that combines Coulomb and Exchange terms calculated with Hermite Gaussians with the polarization, bonded, and modified van der Waals terms from AMOEBA is presented. GEM* is tested on an initial water model fitted at the same level as AMOEBA. The integrals required for the evaluation of the intermolecular Coulomb interactions are efficiently evaluated by means of reciprocal space methods. The GEM* water model is tested by comparing energies and forces for a series of water oligomers and MD simulations. Timings for GEM* compared to AMOEBA are presented and discussed.

SUBMITTER: Duke RE

PROVIDER: S-EPMC5207213 | biostudies-literature | 2014 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

GEM*: A Molecular Electronic Density-Based Force Field for Molecular Dynamics Simulations.

Duke Robert E RE Starovoytov Oleg N ON Piquemal Jean-Philip JP Cisneros G Andrés GA

Journal of chemical theory and computation 20140303 4

GEM*, a force field that combines Coulomb and Exchange terms calculated with Hermite Gaussians with the polarization, bonded, and modified van der Waals terms from AMOEBA is presented. GEM* is tested on an initial water model fitted at the same level as AMOEBA. The integrals required for the evaluation of the intermolecular Coulomb interactions are efficiently evaluated by means of reciprocal space methods. The GEM* water model is tested by comparing energies and forces for a series of water oli ...[more]

PMID: 26580355

Similar Datasets

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids.

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

Nanocomposite of Fullerenes and Natural Rubbers: MARTINI Force Field Molecular Dynamics Simulations.

Project description:The mechanical properties of natural rubber (NR) composites depend on many factors, including the filler loading, filler size, filler dispersion, and filler-rubber interfacial interactions. Thus, NR composites with nano-sized fillers have attracted a great deal of attention for improving properties such as stiffness, chemical resistance, and high wear resistance. Here, a coarse-grained (CG) model based on the MARTINI force field version 2.1 has been developed and deployed for simulations of cis-1,4-polyisoprene (cis-PI). The model shows qualitative and quantitative agreement with the experiments and atomistic simulations. Interestingly, only a 0.5% difference with respect to the experimental result of the glass transition temperature (Tg) of the cis-PI in the melts was observed. In addition, the mechanical and thermodynamical properties of the cis-PI-fullerene(C60) composites were investigated. Coarse-grained molecular dynamics (MD) simulations of cis-PI-C60 composites with varying fullerene concentrations (0-32 parts per hundred of rubber; phr) were performed over 200 microseconds. The structural, mechanical, and thermal properties of the composites were determined. The density, bulk modulus, thermal expansion, heat capacity, and Tg of the NR composites were found to increase with increasing C60 concentration. The presence of C60 resulted in a slight increasing of the end-to-end distance and radius of the gyration of the cis-PI chains. The contribution of C60 and cis-PI interfacial interactions led to an enhancement of the bulk moduli of the composites. This model should be helpful in the investigations and design of effective fillers of NR-C60 composites for improving their properties.

| S-EPMC8626026 | biostudies-literature

Molecular dynamics simulations of Guanine quadruplex loops: advances and force field limitations.

Project description:A computational analysis of d(GGGGTTTTGGGG)(2) guanine quadruplexes containing either lateral or diagonal four-thymidine loops was carried out using molecular dynamics (MD) simulations in explicit solvent, locally enhanced sampling (LES) simulations, systematic conformational search, and free energy molecular-mechanics, Poisson Boltzmann, surface area (MM-PBSA) calculations with explicit inclusion of structural monovalent cations. The study provides, within the approximations of the applied all-atom additive force field, a qualitatively complete analysis of the available loop conformational space. The results are independent of the starting structures. Major conformational transitions not seen in conventional MD simulations are observed when LES is applied. The favored LES structures consistently provide lower free energies (as estimated by molecular-mechanics, Poisson Boltzmann, surface area) than other structures. Unfortunately, the predicted optimal structure for the diagonal loop arrangement differs substantially from the atomic resolution experiments. This result is attributed to force field deficiencies, such as the potential misbalance between solute-cation and solvent-cation terms. The MD simulations are unable to maintain the stable coordination of the monovalent cations inside the diagonal loops as reported in a recent x-ray study. The optimal diagonal and lateral loop arrangements appear to be close in energy although a proper inclusion of the loop monovalent cations could stabilize the diagonal architecture.

| S-EPMC1304345 | biostudies-literature

Molecular Dynamics Simulations of Phosphorylated Intrinsically Disordered Proteins: A Force Field Comparison.

Project description:Phosphorylation is a common post-translational modification among intrinsically disordered proteins and regions, which helps regulate function by changing the protein conformations, dynamics, and interactions with binding partners. To fully comprehend the effects of phosphorylation, computer simulations are a helpful tool, although they are dependent on the accuracy of the force field used. Here, we compared the conformational ensembles produced by Amber ff99SB-ILDN+TIP4P-D and CHARMM36m, for four phosphorylated disordered peptides ranging in length from 14-43 residues. CHARMM36m consistently produced more compact conformations with a higher content of bends, mainly due to more stable salt bridges. Based on comparisons with experimental size estimates for the shortest and longest peptide, CHARMM36m appeared to overestimate the compactness. The difference between the force fields was largest for the peptide showing the greatest separation between positively charged and phosphorylated residues, in line with the importance of charge distribution. For this peptide, the conformational ensemble did not change significantly upon increasing the ionic strength from 0 mM to 150 mM, despite a reduction of the salt-bridging probability in the CHARMM36m simulations, implying that salt concentration has negligible effects in this study.

| S-EPMC8470740 | biostudies-literature

Parameterization and optimization of the menthol force field for molecular dynamics simulations.

Project description:Menthol's various biological properties render it a useful component for medical and cosmetological applications, while its three centers of asymmetry mean that it can be used in a range of organic reactions. Menthol-substituted ionic liquids (ILs) have been found to exhibit promising antimicrobial and antielectrostatic properties, as well as being useful in organic catalysis and biochemical studies. However, so far, a force field designed and validated specifically for the menthol molecule has not been constructed. In the present work, the validation and optimization of force field parameters with regard to the ability to reproduce the macroscopic properties of menthol is presented. The set of optimized potentials for liquid simulations all atom (OPLS-AA) compatible parameters was tested and carefully tuned. The refinement of parameters included fitting of partial atomic charges, optimization of Lennard-Jones parameters, and recalculation of the dihedral angle parameters needed to reproduce quantum energy profiles. To validate the force field, a variety of physicochemical properties were calculated for liquid menthol. Both thermodynamic and kinetic properties were taken into account, including density, surface tension, enthalpy of vaporization, and shear viscosity. The obtained force field was proven to accurately reproduce the properties of the investigated compound while being fully compatible with the OPLS-AA force field.

| S-EPMC5014899 | biostudies-other

ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations.

Project description:Combined quantum mechanical/molecular mechanical (QM/MM) models using semiempirical and ab initio methods have been extensively reported on over the past few decades. These methods have been shown to be capable of providing unique insights into a range of problems, but they are still limited to relatively short time scales, especially QM/MM models using ab initio methods. An intermediate approach between a QM based model and classical mechanics could help fill this time-scale gap and facilitate the study of a range of interesting problems. Reactive force fields represent the intermediate approach explored in this paper. A widely used reactive model is ReaxFF, which has largely been applied to materials science problems and is generally used as a stand-alone (i.e., the full system is modeled using ReaxFF). We report a hybrid ReaxFF/AMBER molecular dynamics (MD) tool, which introduces ReaxFF capabilities to capture bond breaking and formation within the AMBER MD software package. This tool enables us to study local reactive events in large systems at a fraction of the computational costs of QM/MM models. We describe the implementation of ReaxFF/AMBER, validate this implementation using a benzene molecule solvated in water, and compare its performance against a range of similar approaches. To illustrate the predictive capabilities of ReaxFF/AMBER, we carried out a Claisen rearrangement study in aqueous solution. In a first for ReaxFF, we were able to use AMBER's potential of mean force (PMF) capabilities to perform a PMF study on this organic reaction. The ability to capture local reaction events in large systems using combined ReaxFF/AMBER opens up a range of problems that can be tackled using this model to address both chemical and biological processes.

| S-EPMC8145783 | biostudies-literature

High-performance scalable molecular dynamics simulations of a polarizable force field based on classical Drude oscillators in NAMD.

Project description:Incorporating the influence of induced polarization in large-scale atomistic molecular dynamics (MD) simulations is a critical challenge in the progress toward computations of increased accuracy. One computationally efficient treatment is based on the classical Drude oscillator in which an auxiliary charged particle is attached by a spring to each nucleus. Here, we report the first implementation of this model in the program NAMD. An extended Lagrangian dynamics with a dual-Langevin thermostat scheme applied to the Drude-nucleus pairs is employed to efficiently generate classical dynamic propagation near the self-consistent field limit. Large-scale MD simulations based on the Drude polarizable force field scale very well on massively distributed supercomputing platforms, the computational demand being only about 50-100% higher than for nonpolarizable models. As an illustration, a large-scale 150 mM NaCl aqueous salt solution is simulated, and the calculated ionic conductivity is shown to be in excellent agreement with experiment.

| S-EPMC3092300 | biostudies-literature

A polarizable force field of dipalmitoylphosphatidylcholine based on the classical Drude model for molecular dynamics simulations of lipids.

Project description:A polarizable force field of saturated phosphatidylcholine-containing lipids based on the classical Drude oscillator model is optimized and used in molecular dynamics simulations of bilayer and monolayer membranes. The hierarchical parametrization strategy involves the optimization of parameters for small molecules representative of lipid functional groups, followed by their application in larger model compounds and full lipids. The polar headgroup is based on molecular ions tetramethyl ammonium and dimethyl phosphate, the esterified glycerol backbone is based on methyl acetate, and the aliphatic lipid hydrocarbon tails are based on linear alkanes. Parameters, optimized to best represent a collection of gas and liquid properties for these compounds, are assembled into a complete model of dipalmitoylphosphatidylcholine (DPPC) lipids that is tested against the experimental properties of bilayer and monolayer membranes. The polarizable model yields average structural properties that are in broad accord with experimental data. The area per lipid of the model is 60 Å(2), slightly smaller than the experimental value of 63 Å(2). The order parameters from nuclear magnetic resonance deuterium quadrupolar splitting measures, the electron density profile, and the monolayer dipole potential are in reasonable agreement with experimental data, and with the nonpolarizable CHARMM C36 lipid force field.

| S-EPMC3799809 | biostudies-literature

General Force-Field Parametrization Scheme for Molecular Dynamics Simulations of Conjugated Materials in Solution.

Project description:We describe a general scheme to obtain force-field parameters for classical molecular dynamics simulations of conjugated polymers. We identify a computationally inexpensive methodology for calculation of accurate intermonomer dihedral potentials and partial charges. Our findings indicate that the use of a two-step methodology of geometry optimization and single-point energy calculations using DFT methods produces potentials which compare favorably to high level theory calculation. We also report the effects of varying the conjugated backbone length and alkyl side-chain lengths on the dihedral profiles and partial charge distributions and determine the existence of converged lengths above which convergence is achieved in the force-field parameter sets. We thus determine which calculations are required for accurate parametrization and the scope of a given parameter set for variations to a given molecule. We perform simulations of long oligomers of dioctylfluorene and hexylthiophene in explicit solvent and find peristence lengths and end-length distributions consistent with experimental values.

| S-EPMC4980687 | biostudies-literature

Mapping the Protein Fold Universe Using the CamTube Force Field in Molecular Dynamics Simulations.

Project description:It has been recently shown that the coarse-graining of the structures of polypeptide chains as self-avoiding tubes can provide an effective representation of the conformational space of proteins. In order to fully exploit the opportunities offered by such a 'tube model' approach, we present here a strategy to combine it with molecular dynamics simulations. This strategy is based on the incorporation of the 'CamTube' force field into the Gromacs molecular dynamics package. By considering the case of a 60-residue polyvaline chain, we show that CamTube molecular dynamics simulations can comprehensively explore the conformational space of proteins. We obtain this result by a 20 μs metadynamics simulation of the polyvaline chain that recapitulates the currently known protein fold universe. We further show that, if residue-specific interaction potentials are added to the CamTube force field, it is possible to fold a protein into a topology close to that of its native state. These results illustrate how the CamTube force field can be used to explore efficiently the universe of protein folds with good accuracy and very limited computational cost.

| S-EPMC4624779 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data