Enhanced stability of the model mini-protein in amino acid ionic liquids and their aqueous solutions.

Using molecular dynamics simulations, the structure of model mini-protein was thoroughly characterized in the imidazolium-based amino acid ionic liquids and their aqueous solutions. Complete substitution of water by organic cations and anions further results in hindered conformational flexibility of the mini-protein. This observation suggests that amino acid-based ionic liquids are able to defend proteins from thermally induced denaturation. We show by means of radial distributions that the mini-protein is efficiently solvated by both solvents due to a good mutual miscibility. Amino acid-based anions prevail in the first coordination sphere of positively charged sites of the mini-protein whereas water molecules prevail in the first coordination sphere of negatively charged sites of the mini-protein.

Model-free simulation approach to molecular diffusion tensors.

In the present work, we propose a simple model-free approach for the computation of molecular diffusion tensors from molecular dynamics trajectories. The method uses a rigid body trajectory of the molecule under consideration, which is constructed a posteriori by an accumulation of quaternion-based superposition fits of consecutive conformations. From the rigid body trajectory, we compute the translational and angular velocities of the molecule and by integration of the latter also the corresponding angular trajectory. All quantities can be referred to the laboratory frame and a molecule-fixed frame. The 6 × 6 diffusion tensor is computed from the asymptotic slope of the tensorial mean square displacement and, for comparison, also from the Kubo integral of the velocity correlation tensor. The method is illustrated for two simple model systems - a water molecule and a lysozyme molecule in bulk water. We give estimations of the statistical accuracy of the calculations.

Impact of anisotropic atomic motions in proteins on powder-averaged incoherent neutron scattering intensities.

This paper addresses the question to which extent anisotropic atomic motions in proteins impact angular-averaged incoherent neutron scattering intensities, which are typically recorded for powder samples. For this purpose, the relevant correlation functions are represented as multipole series in which each term corresponds to a different degree of intrinsic motional anisotropy. The approach is illustrated by a simple analytical model and by a simulation-based example for lysozyme, considering in both cases the elastic incoherent structure factor. The second example shows that the motional anisotropy of the protein atoms is considerable and contributes significantly to the scattering intensity.

Molecular dynamics and kinetic study of carbon coagulation in the release wave of detonation products.

We present a combined molecular dynamics and kinetic study of a carbon cluster aggregation process in thermodynamic conditions relevant for the detonation products of oxygen deficient explosives. Molecular dynamics simulations with the LCBOPII potential under gigapascal pressure and high temperatures indicate that (i) the cluster motion in the detonation gas is compatible with Brownian diffusion and (ii) the coalescence probability is 100% for two clusters entering the interaction cutoff distance. We used these results for a subsequent kinetic study with the Smoluchowski model, with realistic models applied for the physical parameters such as viscosity and cluster size. We found that purely aggregational kinetics yield too fast clustering, with moderate influence of the model parameters. In agreement with previous studies, the introduction of surface reactivity through a simple kinetic model is necessary to approach the clustering time scales suggested by experiments (1000 atoms after 100 ns, 10 000 atoms after 1 µs). However, these models fail to reach all experimental criteria simultaneously and more complex modelling of the surface process seems desirable to go beyond these current limitations.

Least constraint approach to the extraction of internal motions from molecular dynamics trajectories of flexible macromolecules.

We propose a rigorous method for removing rigid-body motions from a given molecular dynamics trajectory of a flexible macromolecule. The method becomes exact in the limit of an infinitesimally small sampling step for the input trajectory. In a recent paper [G. Kneller, J. Chem. Phys. 128, 194101 (2008)], one of us showed that virtual internal atomic displacements for small time increments can be derived from Gauss' principle of least constraint, which leads to a rotational superposition problem for the atomic coordinates in two consecutive time frames of the input trajectory. Here, we demonstrate that the accumulation of these displacements in a molecular-fixed frame, which evolves in time according to the virtual rigid-body motions, leads to the desired trajectory for internal motions. The atomic coordinates in the input and output trajectory are related by a roto-translation, which guarantees that the internal energy of the molecule is left invariant. We present a convenient implementation of our method, in which the accumulation of the internal displacements is performed implicitly. Two numerical examples illustrate the difference to the classical approach for removing macromolecular rigid-body motions, which consists of aligning its configurations in the input trajectory with a fixed reference structure.

The role of caveolin-1 in lipid droplets and their biogenesis.

We address unresolved questions of the energetics and mechanism of lipid droplet (LD) biogenesis, and of the role of caveolins in the endoplasmic reticulum (ER) and in mature LDs. LDs are eukaryotic repositories of neutral lipids, which are believed to be synthesised in the ER. We investigate the effects of a curvature-inducing protein, caveolin-1, on the formation and structure of a spontaneously aggregated triolein (TO) lipid lens in a flat lipid bilayer using molecular dynamics (MD) simulations. A truncated form of caveolin-1 (Cav1) localises on the interface between the spontaneously formed TO aggregate and the bulk bilayer, and thins the bilayer at the edge of the aggregate, which may contribute to lowering the energy barrier for pinching off the aggregate from the host bilayer. Simulations of fully mature LDs do not conclusively establish the optimal localisation of Cav1 in LDs, but when Cav1 is in the LD core, the distribution of both neutral lipids in the LD core, and of phospholipids on the engulfing monolayer are altered significantly. Our simulations provide an unprecedented molecular description of the distribution and dynamics of various lipid species in both mature LDs and in the nascent LD inside the bilayer.

Protein remains stable at unusually high temperatures when solvated in aqueous mixtures of amino acid based ionic liquids.

Using molecular dynamics simulations, we investigated the thermal stability and real-time denaturation of a model mini-protein in four solvents: (1) water, (2) 1-ethyl-3-methylimidazolium alaninate [EMIM][ALA] (5 mol% in water), (3) methioninate [EMIM][MET] (5 mol% in water), and (4) tryptophanate [EMIM][TRP] (5 mol% in water). Upon analyzing the radius of gyration, the solvent-accessible surface area, root-mean-squared deviations, and inter- and intramolecular hydrogen bonds, we found that the mini-protein remains stable at 30-40 K higher temperatures in aqueous amino acid based ionic liquids (AAILs) than in water. This thermal stability was correlated with the thermodynamics and shear viscosity of the AAIL-containing mixtures. These results suggest that AAILs are generally favorable for protein conservation. Graphical Abstract Conformation of the [TRP]-cage mini-protein in an aqueous amino acid based ionic liquid (AAIL).