Hydrogen Bonding and π-Stacking: How Reliable are Force Fields? A Critical Evaluation of Force Field Descriptions of Nonbonded Interactions

2009 ◽  
Vol 49 (4) ◽  
pp. 944-955 ◽  
Author(s):  
Robert S. Paton ◽  
Jonathan M. Goodman
Keyword(s):  
Hydrogen Bonding ◽  
Force Field ◽  
Force Fields ◽  
Critical Evaluation ◽  
Nonbonded Interactions ◽  
Π Stacking

scholarly journals Adsorption of CO2, N2, CH4, and their mixtures on silicalite: A critical evaluation of force fields

2020 ◽  
Vol 26 (3) ◽  
pp. 295-308
Author(s):  
Sarah Arvelos ◽  
Thalles Diógenes ◽  
Eponina Hori ◽  
Romanielo Lobato
Keyword(s):  
Force Field ◽  
Molecular Simulation ◽  
Force Fields ◽  
Critical Evaluation ◽  
Computational Power ◽  
Potential Accuracy ◽  
Simulation Results ◽  
Adsorption Of Co2 ◽  
Reliable Software ◽  
Potential Use

The use of molecular simulation has been growing in the field of engineering, fueled not just by the advances in computational power but also on the availability of reliable software. One potential use of molecular simulation is related to the screening of materials for a specific application. The reliability of molecular simulation results depends on the trustworthiness of the force field used, which for engineering purposes should be as simple as possible. This work provides an evaluation of the potential accuracy cost of using simple generic force fields to predict the adsorption of CO2, CH4, N2 and their mixtures on MFI. We employed the GCMC technique for this investigation. Different models and force fields to describe adsorbates and adsorbent were tested. The force fields performances were estimated through comparison with available adsorption experimental data. Transferability was evaluated on the prediction of pure and mixtures adsorption on CHA, LTA and FER. The results showed that a simple force field presented similar performance when compared to a more sophisticated one.

Hierarchy of π-Stacking Determines the Conformational Preference of Bis-Squaraines

2020 ◽  
Author(s):  
Abhishek Singh ◽  
Reman K. Singh ◽  
G Naresh Patwari
Keyword(s):  
Hydrogen Bonding ◽  
Rational Design ◽  
Biological Molecules ◽  
Conformational Space ◽  
X Ray Crystallography ◽  
Simple Strategy ◽  
Induced Folding ◽  
Π Stacking ◽  
Varying Length ◽  
Dynamics Simulations

The rational design of conformationally controlled foldable modules can lead to a deeper insight into the conformational space of complex biological molecules where non-covalent interactions such as hydrogen bonding and π-stacking are known to play a pivotal role. Squaramides are known to have excellent hydrogen bonding capabilities and hence, are ideal molecules for designing foldable modules that can mimic the secondary structures of bio-molecules. The π-stacking induced folding of bis-squaraines tethered using aliphatic primary and secondary-diamine linkers of varying length is explored with a simple strategy of invoking small perturbations involving the length linkers and degree of substitution. Solution phase NMR investigations in combination with molecular dynamics simulations suggest that bis-squaraines predominantly exist as extended conformations. Structures elucidated by X-ray crystallography confirmed a variety of folded and extended secondary conformations including hairpin turns and 𝛽-sheets which are determined by the hierarchy of π-stacking relative to N–H···O hydrogen bonds.

Hierarchy of π-Stacking Determines the Conformational Preference of Bis-Squaraines

2020 ◽  
Author(s):  
Abhishek Singh ◽  
Reman K. Singh ◽  
G Naresh Patwari
Keyword(s):  
Hydrogen Bonding ◽  
Rational Design ◽  
Biological Molecules ◽  
Conformational Space ◽  
X Ray Crystallography ◽  
Simple Strategy ◽  
Induced Folding ◽  
Π Stacking ◽  
Varying Length ◽  
Dynamics Simulations

The rational design of conformationally controlled foldable modules can lead to a deeper insight into the conformational space of complex biological molecules where non-covalent interactions such as hydrogen bonding and π-stacking are known to play a pivotal role. Squaramides are known to have excellent hydrogen bonding capabilities and hence, are ideal molecules for designing foldable modules that can mimic the secondary structures of bio-molecules. The π-stacking induced folding of bis-squaraines tethered using aliphatic primary and secondary-diamine linkers of varying length is explored with a simple strategy of invoking small perturbations involving the length linkers and degree of substitution. Solution phase NMR investigations in combination with molecular dynamics simulations suggest that bis-squaraines predominantly exist as extended conformations. Structures elucidated by X-ray crystallography confirmed a variety of folded and extended secondary conformations including hairpin turns and 𝛽-sheets which are determined by the hierarchy of π-stacking relative to N–H···O hydrogen bonds.

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>

The 1-Naphthol Dimer and Its Surprising Preference for π–π Stacking over Hydrogen Bonding

2019 ◽  
Vol 10 (11) ◽  
pp. 2836-2841 ◽  
Author(s):  
Nathan A. Seifert ◽  
Arsh S. Hazrah ◽  
Wolfgang Jäger
Keyword(s):  
Hydrogen Bonding ◽  
Π Stacking
Sign in / Sign up

Export Citation Format

Share Document