A critical assessment of force field accuracy against NMR data for cyclic peptides containing β-amino acids

2018 ◽  
Vol 20 (23) ◽  
pp. 15807-15816 ◽  
Author(s):  
C. Paissoni ◽  
F. Nardelli ◽  
S. Zanella ◽  
F. Curnis ◽  
L. Belvisi ◽  
...  
Keyword(s):  
Amino Acids ◽  
Force Field ◽  
Cyclic Peptides ◽  
Force Fields ◽  
Critical Assessment ◽  
Isoaspartic Acid ◽  
Nmr Parameters ◽  
Nmr Data ◽  
Conformational Properties

A critical assessment of the reproducibility of NMR parameters of β amino acids pinpoints the major weaknesses of eight widely used force fields in reproducing the equilibrium conformational properties of highly constrained cyclic peptides containing isoAspartic acid.

scholarly journals Adsorption of Amino Acids on Graphene: Assessment of Current Force Fields

2019 ◽  
Author(s):  
Siva Dasetty ◽  
John K. Barrows ◽  
Sapna Sarupria
Keyword(s):  
Amino Acids ◽  
Amino Acid ◽  
Force Field ◽  
Force Fields ◽  
Principal Component ◽  
Free Energies ◽  
Adsorbed State ◽  
Free Energy Of Adsorption ◽  
Density Based Clustering ◽  
Structural Preferences

<div> <div> <div> <p>We compare the free energies of adsorption (∆Aads) and the structural preferences of amino acids obtained using the force fields — Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid–graphene interactions are favorable irrespective of the force field. While the magnitudes of ∆Aads differ between the force fields, the trends in the free energy of adsorption with amino acids are similar across the studied force fields. ∆Aads positively correlates with amino acid–graphene and negatively correlates with graphene–water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that while the thermodynamics of adsorption of proteins and peptides would be described consistently across different force fields, the structural preferences of peptides and proteins on graphene will be force field dependent. </p> </div> </div> </div>

scholarly journals On the faithfulness of molecular mechanics representations of proteins towards quantum-mechanical energy surfaces

2020 ◽  
Vol 10 (6) ◽  
pp. 20190121
Author(s):  
Gerhard König ◽  
Sereina Riniker
Keyword(s):  
Amino Acids ◽  
Force Field ◽  
Molecular Mechanics ◽  
Force Fields ◽  
Mechanical Energy ◽  
Reorganization Energy ◽  
Quantum Mechanical ◽  
Energy Minimum ◽  
Structural Reorganization ◽  
Reorganization Energies

Force fields based on molecular mechanics (MM) are the main computational tool to study the relationship between protein structure and function at the molecular level. To validate the quality of such force fields, high-level quantum-mechanical (QM) data are employed to test their capability to reproduce the features of all major conformational substates of a series of blocked amino acids. The phase-space overlap between MM and QM is quantified in terms of the average structural reorganization energies over all energy minima. Here, the structural reorganization energy is the MM potential-energy difference between the structure of the respective QM energy minimum and the structure of the closest MM energy minimum. Thus, it serves as a measure for the relative probability of visiting the QM minimum during an MM simulation. We evaluate variants of the AMBER, CHARMM, GROMOS and OPLS biomolecular force fields. In addition, the two blocked amino acids alanine and serine are used to demonstrate the dependence of the measured agreement on the QM method, the phase, and the conformational preferences. Blocked serine serves as an example to discuss possible improvements of the force fields, such as including polarization with Drude particles, or using tailored force fields. The results show that none of the evaluated force fields satisfactorily reproduces all energy minima. By decomposing the average structural reorganization energies in terms of individual energy terms, we can further assess the individual weaknesses of the parametrization strategies of each force field. The dominant problem for most force fields appears to be the van der Waals parameters, followed to a lesser degree by dihedral and bonded terms. Our results show that performing a simple QM energy optimization from an MM-optimized structure can be a first test of the validity of a force field for a particular target molecule.

scholarly journals Adsorption of Amino Acids on Graphene: Assessment of Current Force Fields

2019 ◽  
Author(s):  
Siva Dasetty ◽  
John K. Barrows ◽  
Sapna Sarupria
Keyword(s):  
Amino Acids ◽  
Amino Acid ◽  
Force Field ◽  
Force Fields ◽  
Principal Component ◽  
Free Energies ◽  
Adsorbed State ◽  
Free Energy Of Adsorption ◽  
Density Based Clustering ◽  
Structural Preferences

<div> <div> <div> <p>We compare the free energies of adsorption (∆Aads) and the structural preferences of amino acids obtained using the force fields — Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid–graphene interactions are favorable irrespective of the force field. While the magnitudes of ∆Aads differ between the force fields, the trends in the free energy of adsorption with amino acids are similar across the studied force fields. ∆Aads positively correlates with amino acid–graphene and negatively correlates with graphene–water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that while the thermodynamics of adsorption of proteins and peptides would be described consistently across different force fields, the structural preferences of peptides and proteins on graphene will be force field dependent. </p> </div> </div> </div>

scholarly journals Adsorption of Amino Acids on Graphene: Assessment of Current Force Fields

2019 ◽  
Author(s):  
Siva Dasetty ◽  
John K. Barrows ◽  
Sapna Sarupria
Keyword(s):  
Amino Acids ◽  
Amino Acid ◽  
Force Field ◽  
Force Fields ◽  
Principal Component ◽  
Free Energies ◽  
Adsorbed State ◽  
Free Energy Of Adsorption ◽  
Density Based Clustering ◽  
Structural Preferences

<div> <div> <div> <p>We compare the free energies of adsorption (∆Aads) and the structural preferences of amino acids obtained using the force fields — Amberff99SB-ILDN/TIP3P, CHARMM36/modified-TIP3P, OPLS-AA/M/TIP3P, and Amber03w/TIP4P/2005. The amino acid–graphene interactions are favorable irrespective of the force field. While the magnitudes of ∆Aads differ between the force fields, the trends in the free energy of adsorption with amino acids are similar across the studied force fields. ∆Aads positively correlates with amino acid–graphene and negatively correlates with graphene–water interaction energies. Using a combination of principal component analysis and density-based clustering technique, we grouped the structures observed in the graphene adsorbed state. The resulting population of clusters, and the conformation in each cluster indicate that the structures of the amino acid in the graphene adsorbed state vary across force fields. The differences in the conformations of amino acids are more severe in the graphene adsorbed state compared to the bulk state for all the force fields. Our findings suggest that while the thermodynamics of adsorption of proteins and peptides would be described consistently across different force fields, the structural preferences of peptides and proteins on graphene will be force field dependent. </p> </div> </div> </div>

Hydration Free Energies of Organic Molecules in the FreeSolv Database Calculated with Polarized Atom In Molecules Atomic Charges and the GAFF Force Field.

2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  
Keyword(s):  
Force Field ◽  
Functional Groups ◽  
Electrostatic Interactions ◽  
Organic Molecules ◽  
Force Fields ◽  
Chemical Elements ◽  
Atomic Charges ◽  
Free Energies ◽  
Molecular Systems ◽  
Solvent Model

<div>Computer simulations of bio-molecular systems often use force fields, which are combinations of simple empirical atom-based functions to describe the molecular interactions. Even though polarizable force fields give a more detailed description of intermolecular interactions, nonpolarizable force fields, developed several decades ago, are often still preferred because of their reduced computation cost. Electrostatic interactions play a major role in bio-molecular systems and are therein described by atomic point charges.</div><div>In this work, we address the performance of different atomic charges to reproduce experimental hydration free energies in the FreeSolv database in combination with the GAFF force field. Atomic charges were calculated by two atoms-in-molecules approaches, Hirshfeld-I and Minimal Basis Iterative Stockholder (MBIS). To account for polarization effects, the charges were derived from the solute's electron density computed with an implicit solvent model and the energy required to polarize the solute was added to the free energy cycle. The calculated hydration free energies were analyzed with an error model, revealing systematic errors associated with specific functional groups or chemical elements. The best agreement with the experimental data is observed for the MBIS atomic charge method, including the solvent polarization, with a root mean square error of 2.0 kcal mol<sup>-1</sup> for the 613 organic molecules studied. The largest deviation was observed for phosphor-containing molecules and the molecules with amide, ester and amine functional groups.</div>

QUBEKit: Automating the Derivation of Force Field Parameters from Quantum Mechanics

2019 ◽  
Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole
Keyword(s):  
Quantum Mechanics ◽  
Force Field ◽  
Organic Molecules ◽  
Force Fields ◽  
Liquid Density ◽  
Small Organic Molecules ◽  
Automated Generation ◽  
Energy Of Hydration ◽  
Novel Method ◽  
Force Field Parameters

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>

QUBEKit: Automating the Derivation of Force Field Parameters from Quantum Mechanics

2019 ◽  
Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole
Keyword(s):  
Quantum Mechanics ◽  
Force Field ◽  
Organic Molecules ◽  
Force Fields ◽  
Liquid Density ◽  
Small Organic Molecules ◽  
Automated Generation ◽  
Energy Of Hydration ◽  
Novel Method ◽  
Force Field Parameters

<div><div><div><p>Modern molecular mechanics force fields are widely used for modelling the dynamics and interactions of small organic molecules using libraries of transferable force field parameters. For molecules outside the training set, parameters may be missing or inaccurate, and in these cases, it may be preferable to derive molecule-specific parameters. Here we present an intuitive parameter derivation toolkit, QUBEKit (QUantum mechanical BEspoke Kit), which enables the automated generation of system-specific small molecule force field parameters directly from quantum mechanics. QUBEKit is written in python and combines the latest QM parameter derivation methodologies with a novel method for deriving the positions and charges of off-center virtual sites. As a proof of concept, we have re-derived a complete set of parameters for 109 small organic molecules, and assessed the accuracy by comparing computed liquid properties with experiment. QUBEKit gives highly competitive results when compared to standard transferable force fields, with mean unsigned errors of 0.024 g/cm3, 0.79 kcal/mol and 1.17 kcal/mol for the liquid density, heat of vaporization and free energy of hydration respectively. This indicates that the derived parameters are suitable for molecular modelling applications, including computer-aided drug design.</p></div></div></div>

A Transferable, Polarisable Force Field for Ionic Liquids

2019 ◽  
Author(s):  
Kateryna Goloviznina ◽  
José N. Canongia Lopes ◽  
Margarida Costa Gomes ◽  
Agilio Padua
Keyword(s):  
Ionic Liquids ◽  
Force Field ◽  
Force Fields ◽  
Dipole Moments ◽  
Van Der Waals Interactions ◽  
First Principles Calculations ◽  
Ion Diffusion ◽  
Fixed Integer ◽  
Molecular Compounds ◽  
Induced Dipoles

A general, transferable polarisable force field for molecular simulation of ionic liquids and their mixtures with molecular compounds is developed. This polarisable model is derived from the widely used CL\&P fixed-charge force field that describes most families of ionic liquids, in a form compatible with OPLS-AA, one of the major force fields for organic compounds. Models for ionic liquids with fixed, integer ionic charges lead to pathologically slow dynamics, a problem that is corrected when polarisation effects are included explicitly. In the model proposed here, Drude induced dipoles are used with parameters determined from atomic polarisabilities. The CL\&P force field is modified upon inclusion of the Drude dipoles, to avoid double-counting of polarisation effects. This modification is based on first-principles calculations of the dispersion and induction contributions to the van der Waals interactions, using symmetry-adapted perturbation theory (SAPT) for a set of dimers composed of positive, negative and neutral fragments representative of a wide variety of ionic liquids. The fragment approach provides transferability, allowing the representation of a multitude of cation and anion families, including different functional groups, without need to re-parametrise. Because SAPT calculations are expensive an alternative predictive scheme was devised, requiring only molecular properties with a clear physical meaning, namely dipole moments and atomic polarisabilities. The new polarisable force field, CL\&Pol, describes a broad set set of ionic liquids and their mixtures with molecular compounds, and is validated by comparisons with experimental data on density, ion diffusion coefficients and viscosity. The approaches proposed here can also be applied to the conversion of other fixed-charged force fields into polarisable versions.<br>

Synthesis and characteristics of ethylenediamine-N,N'-di-S-α-isovalerate and isomers of its cobalt(III) complexes

1982 ◽  
Vol 47 (1) ◽  
pp. 210-216 ◽  
Author(s):  
Milan Strašák ◽  
František Bachratý ◽  
Jaroslav Majer
Keyword(s):  
Amino Acids ◽  
Absolute Configuration ◽  
Electronic Absorption ◽  
Electronic Absorption Spectra ◽  
Chemical Parameters ◽  
Methyl Group ◽  
Nmr Data ◽  
Physico Chemical ◽  
Natural Amino Acids ◽  
C Nmr

The synthesis and physico-chemical parameters are described of a new complexone based on natural amino acids, viz. ethylenediamine-N,N'-di-S-α-isovalerate (SS-EDDIV). 1H- and 13C-NMR data revealed that the methyl group in the substance are not equivalent. The isomers of the cobalt(III) complex with the asymmetric tetradentate SS-EDDIV ligand were prepared and separated; their characteristics are given. The absolute configuration of two of the five theoretically feasible isomers was determined based on their electronic absorption spectra and circular dichroism data.

Sign in / Sign up

Export Citation Format

Share Document