Browse
Submit Data
Databases
API
Help

Dataset Information

30 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Optimal Solution to the Torsional Coefficient Fitting Problem in Force Field Parametrization.

ABSTRACT: 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.

SUBMITTER: Kania A

PROVIDER: S-EPMC8041298 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

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

Parametrization of macrolide antibiotics using the force field toolkit.

Project description:Macrolides are an important class of antibiotics that target the bacterial ribosome. Computer simulations of macrolides are limited as specific force field parameters have not been previously developed for them. Here, we determine CHARMM-compatible force field parameters for erythromycin, azithromycin, and telithromycin, using the force field toolkit (ffTK) plugin in VMD. Because of their large size, novel approaches for parametrizing them had to be developed. Two methods for determining partial atomic charges, from interactions with TIP3P water and from the electrostatic potential, as well as several approaches for fitting the dihedral parameters were tested. The performance of the different parameter sets was evaluated by molecular dynamics simulations of the macrolides in ribosome, with a distinct improvement in maintenance of key interactions observed after refinement of the initial parameters. Based on the results of the macrolide tests, recommended procedures for parametrizing very large molecules using ffTK are given.

| S-EPMC4575642 | biostudies-literature

Force Field Parametrization of Metal Ions from Statistical Learning Techniques.

Project description:A novel statistical procedure has been developed to optimize the parameters of nonbonded force fields of metal ions in soft matter. The criterion for the optimization is the minimization of the deviations from ab initio forces and energies calculated for model systems. The method exploits the combination of the linear ridge regression and the cross-validation techniques with the differential evolution algorithm. Wide freedom in the choice of the functional form of the force fields is allowed since both linear and nonlinear parameters can be optimized. In order to maximize the information content of the data employed in the fitting procedure, the composition of the training set is entrusted to a combinatorial optimization algorithm which maximizes the dissimilarity of the included instances. The methodology has been validated using the force field parametrization of five metal ions (Zn2+, Ni2+, Mg2+, Ca2+, and Na+) in water as test cases.

| S-EPMC5763284 | biostudies-literature

Force-field parametrization based on radial and energy distribution functions.

Project description:We propose a novel force-field-parametrization procedure that fits the parameters of potential functions in a manner that the pair distribution function (DF) of molecules derived from candidate parameters can reproduce the given target DF. Conventionally, approaches to minimize the difference between the candidate and target DFs employ radial DFs (RDF). RDF itself has been reported to be insufficient for uniquely identifying the parameters of a molecule. To overcome the weakness, we introduce energy DF (EDF) as a target DF, which describes the distribution of the pairwise energy of molecules. We found that the EDF responds more sensitively to a small perturbation in the pairwise potential parameters and provides better fitting accuracy compared to that of RDF. These findings provide valuable insights into a wide range of coarse graining methods, which determine parameters using information obtained from a higher-level calculation than that of the developed force field. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

| S-EPMC6790663 | biostudies-literature

Automated conformational energy fitting for force-field development.

Project description:We present a general conformational-energy fitting procedure based on Monte Carlo simulated annealing (MCSA) for application in the development of molecular mechanics force fields. Starting with a target potential energy surface and an unparametrized molecular mechanics potential energy surface, an optimized set of either dihedral or grid-based correction map (CMAP) parameters is produced that minimizes the root mean squared error RMSE between the parametrized and targeted energies. The fitting is done using an MCSA search in parameter space and consistently converges to the same RMSE irrespective of the randomized parameters used to seed the search. Any number of dihedral parameters can be simultaneously parametrized, allowing for fitting to multi-dimensional potential energy scans. Fitting options for dihedral parameters include non-uniform weighting of the target data, constraining multiple optimized parameters to the same value, constraining parameters to be no greater than a user-specified maximum value, including all or only a subset of multiplicities defining the dihedral Fourier series, and optimization of phase angles in addition to force constants. The dihedral parameter fitting algorithm's performance is characterized through multi-dimensional fitting of cyclohexane, tetrahydropyran, and hexopyranose monosaccharide energetics, with the latter case having a 30-dimensional parameter space. The CMAP fitting is applied in the context of polypeptides, and is used to develop a parametrization that simultaneously captures the phi,psi energetics of the alanine dipeptide and the alanine tetrapeptide. Because the dihedral energy term is common to many force fields, we have implemented the dihedral-fitting algorithm in the portable Python scripting language and have made it freely available as "fit_dihedral.py" for download at http://mackerell.umaryland.edu.

| S-EPMC2864003 | biostudies-literature

Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field.

Project description:The development and validation of new peptide dihedral parameters are reported for the OPLS-AA force field. High accuracy quantum chemical methods were used to scan φ, ψ, χ1, and χ2 potential energy surfaces for blocked dipeptides. New Fourier coefficients for the dihedral angle terms of the OPLS-AA force field were fit to these surfaces, utilizing a Boltzmann-weighted error function and systematically examining the effects of weighting temperature. To prevent overfitting to the available data, a minimal number of new residue-specific and peptide-specific torsion terms were developed. Extensive experimental solution-phase and quantum chemical gas-phase benchmarks were used to assess the quality of the new parameters, named OPLS-AA/M, demonstrating significant improvement over previous OPLS-AA force fields. A Boltzmann weighting temperature of 2000 K was determined to be optimal for fitting the new Fourier coefficients for dihedral angle parameters. Conclusions are drawn from the results for best practices for developing new torsion parameters for protein force fields.

| S-EPMC4504185 | biostudies-literature

Infinite Dilution Activity Coefficients as Constraints for Force Field Parametrization and Method Development.

Project description:Molecular simulations begin with an underlying energy model or force field and from this can predict diverse physical properties. However, force fields were often developed with relatively limited data sets, yet accuracy for diverse properties across a broad chemical space is desirable; therefore, tests of such accuracy are particularly important. Here, to this end, we calculated 237 infinite dilution activity coefficients (IDACs), comparing with experimental values from NIST's ThermoML database. We found that calculated IDAC values correlate strongly with experiment (Pearson R of 0.92 ± 0.01) and allow us to identify specific functional groups that appear to present challenges to the force field employed. One potentially valuable aspect of IDACs, as compared to solvation free energies, which have been frequently employed as force field tests, is that the same molecules serve both as solutes and solvents in different cases, allowing us to ensure that force fields are not overly tuned to one particular environment or solvent.

| S-EPMC6663494 | biostudies-literature

Parametrization of an Orbital-Based Linear-Scaling Quantum Force Field for Noncovalent Interactions.

Project description:We parametrize a linear-scaling quantum mechanical force field called mDC for the accurate reproduction of nonbonded interactions. We provide a new benchmark database of accurate ab initio interactions between sulfur-containing molecules. A variety of nonbond databases are used to compare the new mDC method with other semiempirical, molecular mechanical, ab initio, and combined semiempirical quantum mechanical/molecular mechanical methods. It is shown that the molecular mechanical force field significantly and consistently reproduces the benchmark results with greater accuracy than the semiempirical models and our mDC model produces errors twice as small as the molecular mechanical force field. The comparisons between the methods are extended to the docking of drug candidates to the Cyclin-Dependent Kinase 2 protein receptor. We correlate the protein-ligand binding energies to their experimental inhibition constants and find that the mDC produces the best correlation. Condensed phase simulation of mDC water is performed and shown to produce O-O radial distribution functions similar to TIP4P-EW.

| S-EPMC3985928 | biostudies-literature

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model.

Project description:In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

| S-EPMC7147463 | biostudies-literature

FFParam: Standalone package for CHARMM additive and Drude polarizable force field parametrization of small molecules.

Project description:Accurate force-field (FF) parameters are key to reliable prediction of properties obtained from molecular modeling (MM) and molecular dynamics (MD) simulations. With ever-widening applicability of MD simulations, robust parameters need to be generated for a wider range of chemical species. The CHARMM General Force Field program (CGenFF, https://cgenff.umaryland.edu/) is a tool for obtaining initial parameters for a given small molecule based on analogy with the available CGenFF parameters. However, improvement of these parameters is often required and performing their optimization remains tedious and time consuming. In addition, tools for optimization of small molecule parameters in the context of the Drude polarizable FF are not yet available. To overcome these issues, the FFParam package has been designed to facilitate the parametrization process. The package includes a graphical user interface (GUI) created using Qt libraries. FFParam supports Gaussian and Psi4 for performing quantum mechanical calculations and CHARMM and OpenMM for MM calculations. A Monte Carlo simulated annealing (MCSA) algorithm has been implemented for automated fitting of partial atomic charge, atomic polarizabilities and Thole scale parameters. The LSFITPAR program is called for automated fitting of bonded parameters. Accordingly, FFParam provides all the features required for generation and analysis of CHARMM and Drude FF parameters for small molecules. FFParam-GUI includes a text editor, graph plotter, molecular visualization, and text to table converter to meet various requirements of the parametrization process. It is anticipated that FFParam will facilitate wider use of CGenFF as well as promote future use of the Drude polarizable FF.

| S-EPMC7323454 | 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