Structural Insights into Self-Assembled Aerosol-OT Aggregates in Aqueous Media Using Atomistic Molecular Dynamics.

In water, the surfactant dioctyl sulfosuccinate (Aerosol-OT or AOT) exhibits diverse aggregate structures, ranging from micelles to lamella. An atomic-level understanding, however, of the formation and structure of these aggregates is lacking. Herein, using atomistic molecular dynamics (MD) with microsecond-long simulations, self-assembly of AOT in water is studied for concentrations of 1, 7.2, and 20 wt % at 293 K and for 7.2 wt % at 353 K. Assembly proceeds through stepwise association and dissociation of single AOT molecules, and the fusion and fission of AOT clusters. At 293 K, AOT self-assembles into either (i) spherical micelles (1 wt %), (ii) biphasic systems consisting of rod-like and prolate spheroidal micelles (7.2 wt %), or (iii) bilayers (20 wt %). We hypothesize that the observed rod-like structure is a precursor to lamellar microdomains found experimentally in biphasic dispersions. Increasing temperature to 353 K at 7.2 wt % results in a system consisting of prolate micelles but no rod-like micelles. Simulated phase behavior agrees with previously published experimental observations. Individual aggregates formed during self-assembly are identified using graph theory. Structural metrics of these aggregates like the radius of gyration, shape anisotropy, and prolateness are presented. Trends in structural metrics quantitatively reflect how shapes and sizes of AOT aggregates vary with surfactant concentration and temperature. These simulations provide deeper insight into open questions in the scientific community and demonstrate a method to generate physics-based micelle structures that can be used to rationalize experimental observations.

Interpretation of Interfacial Protein Spectra with Enhanced Molecular Simulation Ensembles.

An atomistically detailed picture of protein folding at interfaces can effectively be obtained by comparing interface-sensitive spectroscopic techniques to molecular simulations. Here, we present an extensive evaluation of the capability of contemporary force fields to model protein folding at air-water interfaces with a general scheme for sampling and reweighting theoretical conformational ensembles of interfacial peptides. Force field combinations of CHARMM22*/TIP3P and AMBER99SB*-ILDN/SPC/E were found to reproduce experimental observations best.

The Structure of the Diatom Silaffin Peptide R5 within Freestanding Two-Dimensional Biosilica Sheets.

The silaffin peptide R5 is instrumental to the mineralization of silica cell walls of diatom organisms. The peptide is also widely employed in biotechnology, for example, in the encapsulation of enzymes and for fusion proteins in tissue regeneration. Despite its scientific and technological importance, the interfacial structure of R5 during silica precipitation remains poorly understood. We herein elucidate the conformation of the peptide in its active form within silica sheets by interface-specific vibrational spectroscopy in combination with molecular dynamics simulations. Contrary to previous solution-state NMR studies, our data confirm that R5 maintains a defined structure when interacting with extended silica sheets. We show that the entire amino acid sequence of R5 interacts with silica during silica formation, leading to the intercalation of silica into the assembled peptide film.

Acetylation dictates the morphology of nanophase biosilica precipitated by a 14-amino acid leucine-lysine peptide.

N-terminal acetylation is a commonly used modification technique for synthetic peptides, mostly applied for reasons of enhanced stability, and in many cases regarded as inconsequential. In engineered biosilification - the controlled deposition of silica for nanotechnology applications by designed peptides - charged groups often play a deciding role. Here we report that changing the charge by acetylation of a 14-amino acid leucine-lysine (LK) peptide dramatically changes the morphology of precipitated biosilica; acetylated LK peptides produce nano-spheres, whereas nano-wires are precipitated by the same peptide in a non-acetylated form. By using interface-specific vibrational spectroscopy and coarse-grained molecular simulations, we show that this change in morphology is not the result of modified peptide-silica interactions, but rather caused by the stabilization of the hydrophobic core of peptide aggregates created by the removal of a peptide charge upon acetylation. These results should raise awareness of the potential impact of N-terminal modifications in peptide applications. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.

Destabilization of Human Serum Albumin by Ionic Liquids Studied Using Enhanced Molecular Dynamics Simulations.

Ionic liquid (IL) containing solvents can change the structure, dynamics, function, and stability of proteins. In order to investigate the mechanisms by which ILs induce structural changes in a large multidomain protein, we study the interactions of human serum albumin (HSA) with two different ILs, 1-butyl-3-methylimidazolium tetrafluoroborate and choline dihydrogen phosphate. Root mean square deviation and fluctuation calculations indicate that high concentrations of ILs in mixtures with water lead to protein structures that remain close to their crystallographic structures on time scales of hundreds of nanoseconds. To overcome potential time scale limitations due to the high viscosity of the solvent, we employed enhanced sampling techniques to estimate the free energy of an experimentally determined important transition within the protein structure. Metadynamics simulations show that the free energy landscape of the unfolding of loop 1 of domain I is different in the presence of ILs than it is in water, consistent with previously published experimental evidence. We then apply essential dynamics coarse graining to systematically predict differences in the dynamics of proteins solvated in IL-water mixtures versus pure water systems. We also demonstrate that the presence of ionic liquids changes the distribution of intermolecular distances among several ligands, indicating that the protein structure swells in the presence of certain ILs, consistent with experimental evidence.

Structure Characterization and Properties of Metal-Surfactant Complexes Dispersed in Organic Solvents.

This work describes the synthesis and characterization of metal-surfactant complexes. Dioctyl sulfosuccinate and dodecylbenzenesulfonate are associated with multivalent aluminum, iron, and vanadium ions using an ion exchange reaction. The metal complexes are dispersible in various organic solvents. In solvents with low polarity, the complexes form "inverse" macromolecular structures with multiple metal ions. In contrast, in alcohols, the complex size is reduced, showing a more disperse conformation. The metal and surfactant ions are still strongly bonded to each other in all the solvents probed. Small-angle X-ray and neutron scattering (SAXS and SANS) are used to characterize the structures. Simultaneous fitting of neutron and X-ray scattering spectra is performed in order to obtain an accurate description of the system. Scattering results are also validated by performing molecular dynamics (MD) simulations. The conductive and electrochemical properties of the complexes in solution are also evaluated. The dispersion of metal-organic complexes significantly increases electric conductivity, and some metal ions in the core of the complexes are shown to be electrochemically active in apolar solvents.

The general AMBER force field (GAFF) can accurately predict thermodynamic and transport properties of many ionic liquids.

We have applied molecular dynamics to calculate thermodynamic and transport properties of a set of 19 room-temperature ionic liquids. Since accurately simulating the thermophysical properties of solvents strongly depends upon the force field of choice, we tested the accuracy of the general AMBER force field, without refinement, for the case of ionic liquids. Electrostatic point charges were developed using ab initio calculations and a charge scaling factor of 0.8 to more accurately predict dynamic properties. The density, heat capacity, molar enthalpy of vaporization, self-diffusivity, and shear viscosity of the ionic liquids were computed and compared to experimentally available data, and good agreement across a wide range of cation and anion types was observed. Results show that, for a wide range of ionic liquids, the general AMBER force field, with no tuning of parameters, can reproduce a variety of thermodynamic and transport properties with similar accuracy to that of other published, often IL-specific, force fields.

Comparison of three ionic liquid-tolerant cellulases by molecular dynamics.

We have employed molecular dynamics to investigate the differences in ionic liquid tolerance among three distinct family 5 cellulases from Trichoderma viride, Thermogata maritima, and Pyrococcus horikoshii. Simulations of the three cellulases were conducted at a range of temperatures in various binary mixtures of the ionic liquid 1-ethyl-3-methyl-imidazolium acetate with water. Our analysis demonstrates that the effects of ionic liquids on the enzymes vary in each individual case from local structural disturbances to loss of much of one of the enzyme's secondary structure. Enzymes with more negatively charged surfaces tend to resist destabilization by ionic liquids. Specific and unique structural changes in the enzymes are induced by the presence of ionic liquids. Disruption of the secondary structure, changes in dynamical motion, and local changes in the binding pocket are observed in less tolerant enzymes. Ionic-liquid-induced denaturation of one of the enzymes is indicated over the 500 ns timescale. In contrast, the most tolerant cellulase behaves similarly in water and in ionic-liquid-containing mixtures. Unlike the heuristic approaches that attempt to predict enzyme stability using macroscopic properties, molecular dynamics allows us to predict specific atomic-level structural and dynamical changes in an enzyme's behavior induced by ionic liquids and other mixed solvents. Using these insights, we propose specific experimentally testable hypotheses regarding the origin of activity loss for each of the systems investigated in this study.

Diatom mimics: directing the formation of biosilica nanoparticles by controlled folding of lysine-leucine peptides.

Silaffins, long chain polyamines, and other biomolecules found in diatoms are involved in the assembly of a large number of silica nanostructures under mild, ambient conditions. Nanofabrication researchers have sought to mimic the diatom's biosilica production capabilities by engineering proteins to resemble aspects of naturally occurring biomolecules. Such mimics can produce monodisperse biosilica nanospheres, but in vitro production of the variety of intricate biosilica nanostructures that compose the diatom frustule is not yet possible. In this study we demonstrate how LK peptides, composed solely of lysine (K) and leucine (L) amino acids arranged with varying hydrophobic periodicities, initiate the formation of different biosilica nanostructures in vitro. When L and K residues are arranged with a periodicity of 3.5 the α-helical form of the LK peptide produces monodisperse biosilica nanospheres. However, when the LK periodicity is changed to 3.0, corresponding to a 310 helix, the morphology of the nanoparticles changes to elongated rod-like structures. ß-strand LK peptides with a periodicity of 2.0 induce wire-like silica morphologies. This study illustrates how the morphology of biosilica can be changed simply by varying the periodicity of polar and nonpolar amino acids.

Structure, dynamics, and activity of xylanase solvated in binary mixtures of ionic liquid and water.

We have discovered that a family 11 xylanase from Trichoderma longibrachiatum maintains significant activity in low concentrations of the ionic liquids (IL) 1-ethyl-3-methyl-imidazolium acetate ([EMIM][OAc]) or 1-ethyl-3-methyl-imidazolium ethyl sulfate ([EMIM][EtSO4]) in water. In order to understand the mechanisms by which the ionic liquids affect the activity of xylanase, we conducted molecular dynamics simulations of the enzyme in various concentrations of the cosolvent. The simulations show that higher concentrations of ionic liquid correlate with less deviation from the starting crystallographic structure. Dynamic motion of the protein is severely dampened by even the lowest tested concentrations of ionic liquid as measured by root-mean-square fluctuation. Principal component analysis shows that the characteristics of the main modes of enzyme motion are greatly affected by the choice of solvent. Cations become kinetically trapped in the binding pocket, allowing them to act as a competitive inhibitor to the natural substrate. Dynamic light scattering and kinetic studies evaluated the stability of the enzyme in the new solvents. These studies indicate that likely factors in the loss of enzyme activity for this xylanase are the dampening of dynamic motion and kinetic trapping of cations in the binding pocket as opposed to the denaturing of the protein.