To test and validate the Automated force field Topology Builder and Repository (ATB; http://compbio.biosci.uq.edu.au/atb/ ) the hydration free enthalpies for a set of 214 drug-like molecules, including 47 molecules that form part of the SAMPL4 challenge have been estimated using thermodynamic integration and compared to experiment. The calculations were performed using a fully automated protocol that incorporated a dynamic analysis of the convergence and integration error in the selection of intermediate points. The system has been designed and implemented such that hydration free enthalpies can be obtained without manual intervention following the submission of a molecule to the ATB. The overall average unsigned error (AUE) using ATB 2.0 topologies for the complete set of 214 molecules was 6.7 kJ/mol and for molecules within the SAMPL4 7.5 kJ/mol. The root mean square error (RMSE) was 9.5 and 10.0 kJ/mol respectively. However, for molecules containing functional groups that form part of the main GROMOS force field the AUE was 3.4 kJ/mol and the RMSE was 4.0 kJ/mol. This suggests it will be possible to further refine the parameters provided by the ATB based on hydration free enthalpies.

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. This in turn enables for a range of target data from multiple molecules to be combined in a robust manner and for the optimal region of parameter space to be trivially identified. The effectiveness of the approach is illustrated by using the method to refine the chlorine 6-12 Lennard-Jones parameters against experimental solvation free enthalpies in water and hexane as well as the density and heat of vaporization of the liquid at atmospheric pressure for a set of 10 aromatic-chloro compounds simultaneously. Single-step perturbation is used to efficiently calculate solvation free enthalpies for a wide range of parameter combinations. The capacity of this approach to parametrize accurate and transferrable force fields is discussed.

Biominerals are crucial materials that play a vital role in many forms of life. Understanding the various steps through which ions in aqueous environment associate to form increasingly structured particles that eventually transform into the final crystalline or amorphous poly(a)morph in the presence of biologically active molecules is therefore of great significance. In this context, computer modeling is now able to provide an accurate atomistic picture of the dynamics and thermodynamics of possible association events in solution, as well as to make predictions as to particle stability and possible alternative nucleation pathways, as a complement to experiment. This review provides a general overview of the most significant computational methods and of their achievements in this field, with a focus on calcium carbonate as the most abundant biomineral.

Since Pasteur first successfully separated right-handed and left-handed tartrate crystals in 1848, the understanding of how homochirality is achieved from enantiomeric mixtures has long been incomplete. Here, we report on a chirality dominance effect where organized, three-dimensional homochiral suprastructures of the biomineral calcium carbonate (vaterite) can be induced from a mixed nonracemic amino acid system. Right-handed (counterclockwise) homochiral vaterite helicoids are induced when the amino acid l -Asp is in the majority, whereas left-handed (clockwise) homochiral morphology is induced when d -Asp is in the majority. Unexpectedly, the Asp that incorporates into the homochiral vaterite helicoids maintains the same enantiomer ratio as that of the initial growth solution, thus showing chirality transfer without chirality amplification. Changes in the degree of chirality of the vaterite helicoids are postulated to result from the extent of majority enantiomer assembly on the mineral surface. These mechanistic insights potentially have major implications for high-level advanced materials synthesis.

The interaction of organic molecules with mineral systems is relevant to a wide variety of scientific problems both in the environment and minerals processing.