scholarly journals Measurements of periodically perturbed dewetting force fields and their consequences on the symmetry of the resulting patterns

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Konstantinos Roumpos ◽  
Sarah Fontaine ◽  
Thomas Pfohl ◽  
Oswald Prucker ◽  
Jürgen Rühe ◽  
...  

AbstractWe studied the origin of breaking the symmetry for moving circular contact lines of dewetting polymer films suspended on a periodic array of pillars. There, dewetting force fields driving polymer flow were perturbed by elastic micro-pillars arranged in a regular square pattern. Elastic restoring forces of deformed pillars locally balance driving capillary forces and broke the circular symmetry of expanding dewetting holes. The observed envelope of the dewetting holes reflected the symmetry of the underlying pattern, even at sizes much larger than the characteristic period of the pillar array, demonstrating that periodic perturbations in a driving force field can establish a well-defined pattern of lower symmetry. For the presented system, we succeeded in squaring the circle.

2018 ◽  
Author(s):  
Maximiliano Riquelme ◽  
Alejandro Lara ◽  
David L. Mobley ◽  
Toon Vestraelen ◽  
Adelio R Matamala ◽  
...  

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


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

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


Author(s):  
Joshua Horton ◽  
Alice Allen ◽  
Leela Dodda ◽  
Daniel Cole

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


2019 ◽  
Author(s):  
Kateryna Goloviznina ◽  
José N. Canongia Lopes ◽  
Margarida Costa Gomes ◽  
Agilio Padua

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>


2017 ◽  
Vol 114 (31) ◽  
pp. 8265-8270 ◽  
Author(s):  
Simon Olsson ◽  
Hao Wu ◽  
Fabian Paul ◽  
Cecilia Clementi ◽  
Frank Noé

Accurate mechanistic description of structural changes in biomolecules is an increasingly important topic in structural and chemical biology. Markov models have emerged as a powerful way to approximate the molecular kinetics of large biomolecules while keeping full structural resolution in a divide-and-conquer fashion. However, the accuracy of these models is limited by that of the force fields used to generate the underlying molecular dynamics (MD) simulation data. Whereas the quality of classical MD force fields has improved significantly in recent years, remaining errors in the Boltzmann weights are still on the order of a few kT, which may lead to significant discrepancies when comparing to experimentally measured rates or state populations. Here we take the view that simulations using a sufficiently good force-field sample conformations that are valid but have inaccurate weights, yet these weights may be made accurate by incorporating experimental data a posteriori. To do so, we propose augmented Markov models (AMMs), an approach that combines concepts from probability theory and information theory to consistently treat systematic force-field error and statistical errors in simulation and experiment. Our results demonstrate that AMMs can reconcile conflicting results for protein mechanisms obtained by different force fields and correct for a wide range of stationary and dynamical observables even when only equilibrium measurements are incorporated into the estimation process. This approach constitutes a unique avenue to combine experiment and computation into integrative models of biomolecular structure and dynamics.


1999 ◽  
Vol 103 (33) ◽  
pp. 6998-7014 ◽  
Author(s):  
Carl S. Ewig ◽  
Thomas S. Thacher ◽  
Arnold T. Hagler

2018 ◽  
Vol 39 (28) ◽  
pp. 2360-2370 ◽  
Author(s):  
Adam K. Sieradzan ◽  
Artur Giełdoń ◽  
Yanping Yin ◽  
Yi He ◽  
Harold A. Scheraga ◽  
...  

2019 ◽  
Author(s):  
Siva Dasetty ◽  
John K. Barrows ◽  
Sapna Sarupria

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


Sign in / Sign up

Export Citation Format

Share Document