On the Validation of Protein Force Fields Based on Structural Criteria.

Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, -coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small.

Facilitating the structural characterisation of non-canonical amino acids in biomolecular NMR.

Peptides and proteins containing non-canonical amino acids (ncAAs) are a large and important class of biopolymers. They include non-ribosomally synthesised peptides, post-translationally modified proteins, expressed or synthesised proteins containing unnatural amino acids, and peptides and proteins that are chemically modified. Here, we describe a general procedure for generating atomic descriptions required to incorporate ncAAs within popular NMR structure determination software such as CYANA, CNS, Xplor-NIH and ARIA.

Engineering Transferable Atomic Force Fields: Empirical Optimization of Hydrocarbon Lennard-Jones Interactions by Direct Mapping of Parameter Space.

The utility of atomistic simulations depends on the accuracy of the force field used to represent the potential energy landscape, the consistency with which interaction parameters can be assigned, and the extent to which parameters can be transferred between chemical entities. Here, parameter space mapping, a simple and robust procedure for atom typing (parameter assignment) and parameter optimization, is used to identify a minimal set of parameters capable of simultaneously reproducing the density, heat of vaporization, and solvation free energies for a targeted set of simple hydrocarbons. Using an atom-centered fixed charge model and a 6-12 Lennard-Jones potential, the experimental densities and the heats of vaporization for 22 hydrocarbons (linear, cyclic, and aromatic) could be predicted with high precision: average unsigned error (AUE) of 6.

OFraMP: a fragment-based tool to facilitate the parametrization of large molecules.

An Online tool for Fragment-based Molecule Parametrization (OFraMP) is described. OFraMP is a web application for assigning atomic interaction parameters to large molecules by matching sub-fragments within the target molecule to equivalent sub-fragments within the Automated Topology Builder (ATB, atb.uq.

Self-cyclisation as a general and efficient platform for peptide and protein macrocyclisation.

Macrocyclisation of proteins and peptides results in a remarkable increase in structural stability, making cyclic peptides and proteins of great interest in drug discovery-either directly as drug leads or as in the case of cyclised nanodiscs (cNDs), as tools for studies of trans-membrane receptors and membrane-active peptides. Various biological methods have been developed that are capable of yielding head-to-tail macrocyclised products. Recent advances in enzyme-catalysed macrocyclisation include discovery of new enzymes or design of new engineered enzymes.

PyThinFilm: Automated Molecular Dynamics Simulation Protocols for the Generation of Thin Film Morphologies.

The performance of organic optoelectronic devices, such as organic light-emitting diodes (OLEDs) and organic solar cells (OSCs), is intrinsically related to the molecular-scale morphology of the thin films from which they are composed. However, the experimental characterization of morphology at the molecular level is challenging due to the often amorphous or at best semicrystalline nature of these films. Classical molecular modeling techniques, such as molecular dynamics (MD) simulation, are increasingly used to understand the relationship between morphology and the properties of thin-film devices.

Understanding the Effect of pH on the Solubility and Aggregation Extent of Humic Acid in Solution by Combining Simulation and the Experiment.

Molecular dynamics (MD) simulations were performed to investigate the dynamics of humic acid (HA) in an aqueous solution and the influence of pH, temperature, and HA concentration. The HA model employed in MD simulations was chosen and validated using experimental chemical composition data and Fourier transform infrared (FTIR) spectra. The simulations showed that the HA molecule has a strong propensity to adopt a compact conformation in water independent of pH, while the aggregation of HA was found to be pH-dependent.

Curved or linear? Predicting the 3-dimensional structure of α-helical antimicrobial peptides in an amphipathic environment.

α-Helical membrane-active antimicrobial peptides (AMPs) are known to act via a range of mechanisms, including the formation of barrel-stave and toroidal pores and the micellisation of the membrane (carpet mechanism). Different mechanisms imply that the peptides adopt different 3D structures when bound at the water-membrane interface, a highly amphipathic environment. Here, an evolutionary algorithm is used to predict the 3D structure of a range of α-helical membrane-active AMPs at the water-membrane interface by optimising amphipathicity.

Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane.

The ability of atomic interaction parameters generated using the Automated Topology Builder and Repository version 3.0 (ATB3.0) to predict experimental hydration free enthalpies (Δ G) and solvation free enthalpies in the apolar solvent hexane (Δ G) is presented.

Predicting the Prevalence of Alternative Warfarin Tautomers in Solution.

Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined.

Optimization of Empirical Force Fields by Parameter Space Mapping: A Single-Step Perturbation Approach.

A general method for parametrizing atomic interaction functions is presented. The method is based on an analysis of surfaces corresponding to the difference between calculated and target data as a function of alternative combinations of parameters (parameter space mapping). The consideration of surfaces in parameter space as opposed to local values or gradients leads to a better understanding of the relationships between the parameters being optimized and a given set of target data.

The Molecular Origin of Anisotropic Emission in an Organic Light-Emitting Diode.

Atomistic nonequilibrium molecular dynamics simulations have been used to model the induction of molecular orientation anisotropy within the emission layer of an organic light-emitting diode (OLED) formed by vapor deposition. Two emitter species were compared: racemic fac-tris(2-phenylpyridine)iridium(III) (Ir(ppy)) and trans-bis(2-phenylpyridine)(acetylacetonate)iridium(III) (Ir(ppy)(acac)). The simulations show that the molecular symmetry axes of both emitters preferentially align perpendicular to the surface during deposition.

Real Cost of Speed: The Effect of a Time-Saving Multiple-Time-Stepping Algorithm on the Accuracy of Molecular Dynamics Simulations.

To enhance efficiency in molecular dynamics simulations, forces that vary slowly are often evaluated less often than those that vary rapidly. We show that the multiple-time-step algorithm implemented in recent versions of GROMACS results in significant differences in the collective properties of a system under conditions where the system was previously stable. The implications of changing the simulation algorithm without assessment of potential artifacts on the parametrization and transferability of effective force fields are discussed.

Elucidating the Spatial Arrangement of Emitter Molecules in Organic Light-Emitting Diode Films.

The effect of varying the emitter concentration on the structural properties of an archetypal phosphorescent blend consisting of 4,4'-bis(N-carbazolyl)biphenyl and tris(2-phenylpyridyl)iridium(III) has been investigated using non-equilibrium molecular dynamics (MD) simulations that mimic the process of vacuum deposition. By comparison with reflectometry measurements, we show that the simulations provide an accurate model of the average density of such films. The emitter molecules were found not to be evenly distributed throughout film, but rather they can form networks that provide charge and/or energy migration pathways, even at emitter concentrations as low as ≈5 weight percent.

Testing and validation of the Automated Topology Builder (ATB) version 2.0: prediction of hydration free enthalpies.

To test and validate the Automated force field Topology Builder and Repository (ATB; http://compbio.biosci.uq.

An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0.

The Automated force field Topology Builder (ATB, http://compbio.biosci.uq.