Optimizing potentials for a liquid mixture: a new force field for a tert-butanol and water solution.
Ontology highlight
ABSTRACT:
SUBMITTER: Di Pierro M
PROVIDER: S-EPMC4306501 | biostudies-other | 2015 Jan
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Molecular simulations can be used to study disordered polypeptide systems and to generate hypotheses on the underlying structural and thermodynamic mechanisms that govern their function. As the number of disordered protein systems investigated with simulations increase, it is important to understand how particular force fields affect the structural properties of disordered polypeptides in solution. To this end, we performed a comparative structural analysis of Gly(3) and Gly(10) in aqueous solution from all atom, microsecond molecular dynamics (MD) simulations using the CHARMM 27 (C27), CHARMM 36 (C36), and Amber ff12SB force fields. For each force field, Gly(3) and Gly(10) were simulated for at least 300 ns and 1 ?s, respectively. Simulating oligoglycines of two different lengths allows us to evaluate how force field effects depend on polypeptide length. Using a variety of structural metrics (e.g., end-to-end distance, radius of gyration, dihedral angle distributions), we characterize the distribution of oligoglycine conformers for each force field and show that each sample conformation space differently, yielding considerably different structural tendencies of the same oligoglycine model in solution. Notably, we find that C36 samples more extended oligoglycine structures than both C27 and ff12SB.
| S-EPMC4450816 | biostudies-literature
Project description:We employ the popular all-atom optimized potential for liquid simulations, OPLSAA, force-field to model 17 different alcohols in the liquid state. Using the standard simulation protocol for few hundred nanosecond time periods, we find that 1-octanol, 1-nonanol, and 1-decanol undergo spontaneous transition to a crystalline state at temperatures which are 35-55 K higher than the experimental melting temperatures. Nevertheless, the crystal structures obtained from the simulations are very similar to those determined by X-ray powder diffraction data for several n-alcohols. Although some degree of deviations from the experimental freezing points are to be expected, for 1-nonanol and 1-decanol, the elevation of the freezing temperature warrants special attention because at room temperature, these alcohols are liquids; however, if simulated by the OPLSAA force-field, they will crystallize. This behavior is likely a consequence of exaggerated attractive interactions between the alkane chains of the alcohols. To circumvent this problem, we combined the OPLSAA model with the L-OPLS force-field. We adopted the L-OPLS parameters to model the hydrocarbon tail of the alcohols, whereas the hydroxyl head group remained as in the original OPLSAA force-field. The resulting alcohols stayed in the liquid state at temperatures above their experimental melting points, thus, resolving the enhanced freezing observed with the OPLSAA force-field. In fact, the mixed-model alcohols did not exhibit any spontaneous freezing even at temperatures much lower than the experimental values. However, a series of simulations in which these mixed-OPLSAA alcohols started from a coexistence configuration of the liquid and solid phases resulted in freezing transitions at temperatures 14-25 K lower than the experimental values, confirming the validity of the proposed model. For all of the other alcohols, the mixed model yields results very similar to the OPLSAA force-field and is in good agreement with the experimental data. Thus, for simulating alcohols in the liquid phase, the mixed-OPLSAA model is necessary for large (7 carbons and above) hydrocarbon chains.
| S-EPMC6643393 | biostudies-literature
Project description:Recent advances in hardware and software have enabled increasingly long molecular dynamics (MD) simulations of biomolecules, exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, we further advance simulation accuracy by improving the amino acid side-chain torsion potentials of the Amber ff99SB force field. First, we used simulations of model alpha-helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, we optimized the side-chain torsion potentials of these residues to match new, high-level quantum-mechanical calculations. Finally, we used microsecond-timescale MD simulations in explicit solvent to validate the resulting force field against a large set of experimental NMR measurements that directly probe side-chain conformations. The new force field, which we have termed Amber ff99SB-ILDN, exhibits considerably better agreement with the NMR data.
| S-EPMC2970904 | biostudies-literature
Project description:We developed a novel neural network-based force field for water based on training with high-level ab initio theory. The force field was built based on an electrostatically embedded many-body expansion method truncated at binary interactions. The many-body expansion method is a common strategy to partition the total Hamiltonian of large systems into a hierarchy of few-body terms. Neural networks were trained to represent electrostatically embedded one-body and two-body interactions, which require as input only one and two water molecule calculations at the level of ab initio electronic structure method CCSD/aug-cc-pVDZ embedded in the molecular mechanics water environment, making it efficient as a general force field construction approach. Structural and dynamic properties of liquid water calculated with our force field show good agreement with experimental results. We constructed two sets of neural network based force fields: nonpolarizable and polarizable force fields. Simulation results show that the nonpolarizable force field using fixed TIP3P charges has already behaved well, since polarization effects and many-body effects are implicitly included due to the electrostatic embedding scheme. Our results demonstrate that the electrostatically embedded many-body expansion combined with neural network provides a promising and systematic way to build next-generation force fields at high accuracy and low computational costs, especially for large systems.
| S-EPMC6241215 | biostudies-literature
Project description:In this paper, we present a new molecular model that can accurately predict thermodynamic liquid state and phase-change properties for organosilicon molecules including several functional groups (alkylsilane, alkoxysilane, siloxane, and silanol). These molecules are of great importance in geological processes, biological systems, and material science, yet no force field currently exists that is widely applicable to organosilicates. The model is parametrized according to the recent Polarization-Consistent Approach (PolCA), which allows for polarization effects to be incorporated into a nonpolarizable model through post facto correction terms and is therefore consistent with previous parametrizations of the PolCA force field. Alkyl groups are described by the United-Atom approach, bond and angle parameters were taken from previous literature studies, dihedral parameters were fitted to new quantum chemical energy profiles, point charges were calculated from quantum chemical optimizations in a continuum solvent, and Lennard-Jones dispersion/repulsion parameters were fitted to match the density and enthalpy of vaporization of a small number of selected compounds. Extensive validation efforts were carried out, after careful collection and curation of experimental data for organosilicates. Overall, the model performed quite well for the density, enthalpy of vaporization, dielectric constant, and self-diffusion coefficient, but it slightly overestimated the magnitude of self-solvation free energies. The modular and transferable nature of the PolCA force field allows for further extensions to other types of silicon-containing compounds.
| S-EPMC8679648 | biostudies-literature
Project description:The basal ganglia are known to be involved in the planning, execution and control of gripping force and movement vigour. Here we aim to define the nature of the basal ganglia control signal for force and to decode gripping force based on local field potential (LFP) activities recorded from the subthalamic nucleus (STN) in patients with deep brain stimulation (DBS) electrodes. We found that STN LFP activities in the gamma (55-90 Hz) and beta (13-30m Hz) bands were most informative about gripping force, and that a first order dynamic linear model with these STN LFP features as inputs can be used to decode the temporal profile of gripping force. Our results enhance the understanding of how the basal ganglia control gripping force, and also suggest that deep brain LFPs could potentially be used to decode movement parameters related to force and movement vigour for the development of advanced human-machine interfaces.
| S-EPMC5148608 | biostudies-literature
Project description:Current Brain Computer Interface (BCI) systems are limited by relying on neuronal spikes and decoding limited to kinematics only. For a BCI system to be practically useful, it should be able to decode brain information on a continuous basis with low latency. This study investigates if force can be decoded from local field potentials (LFP) recorded with deep brain electrodes located at the Subthalamic nucleus (STN) using data from 5 patients with Parkinson's disease, on a continuous basis with low latency. A Wiener-Cascade (WC) model based decoder was proposed using both time-domain and frequency-domain features. The results suggest that high gamma band (300-500Hz) activity, in addition to the beta (13-30Hz) and gamma band (55-90Hz) activity is the most informative for force prediction but combining all features led to better decoding performance. Furthermore, LFP signals preceding the force output by up to 1256 milliseconds were found to be predictive of the force output.
| S-EPMC5557341 | biostudies-other
Project description:Molecular modeling is an excellent tool for studying biological systems on the atomic scale. Depending on objects, which may be proteins, nucleic acids, or lipids, different force fields are recommended. The phospholipid bilayers constitute an example, in which behavior is extensively studied using molecular dynamics simulations due to limitations of experimental methods. The reliability of the results is strongly dependent on an appropriate description of these compounds. There are some deficiencies in the parametrization of intra- and intermolecular interactions that result in incorrect reproduction of phospholipid bilayer properties known from experimental studies, such as temperatures of phase transitions. Refinement of the force field parameters of nonbonded interactions present in the studied system is required to close these discrepancies. Such parameters as partial charges and torsional potential coefficients are crucial in this issue and not obtainable from experimental studies. This work presents a new fitting procedure for torsional coefficients that employs linear algebra theory and compares it with the Monte Carlo method. The proposed algebraic approach can be applied to any considered molecular system. In the manuscript, it is presented on the example of dimethyl phosphoric acid molecule. The advantages of our method encompass finding an optimal solution, the lack of additional parameters required by the algorithm, and significantly shorter computational time. Additionally, we indicate the importance of proper assignment of the partial charges.
| S-EPMC8041298 | 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:To investigate whether water relay plays an important role in phosphoryl transfer reactions, we have used several theoretical approaches to compare key properties of uridine 3'-m-nitrobenzyl phosphate (UNP) in aqueous and tert-butanol solutions. Previous kinetic experiments found that the isomerization reaction of UNP is abolished in tert-butanol, which was interpreted as the direct evidence that supports the role of water relay in phosphoryl transfer. We have analyzed solute flexibility and solvent structure near the solute using equilibrium molecular dynamics simulations and a combined quantum mechanical/molecular mechanism (QM/MM) potential function for the solute. Snapshots from the simulations are then used in minimum energy path calculations to compare the energetics of direct nucleophilic attack and water-mediated nucleophilic attack pathways. QM/MM simulations are also used to compare the pseudorotation barriers for the pentavalent intermediate formed following the nucleophilic attack, another key step for the isomerization reaction. Combined results from these calculations suggest that water relay does not offer any significant energetic advantage over the direct nucleophilic attack. Unfortunately, the lack of isomerization in tert-butanol solution cannot be straightforwardly explained based on the results we have obtained here and therefore requires additional analysis. This study, nevertheless, has provided new insights into several most commonly discussed possibilities.
| S-EPMC2705241 | biostudies-literature