Roland Faller

Molecular Simulation Study of Phospholipid Bilayers and Insights of the Interactions with Disaccharides

Molecular simulations of hydrated dipalmitoylphosphatidylcholine lipid bilayers have been performed for temperatures in the range of 250-450 K. The area per headgroup increases with temperature from 58 to 77 A(2). Other properties such as hydration number, alkyl tail order parameter, diffusion coefficients, and radial distribution functions exhibit a clear dependence on temperature. Simulations of bilayers have also been performed in the presence of two disaccharides, namely trehalose and sucrose, at concentrations of up to 18 wt % (lipid-free basis). The simulated area per headgroup of the bilayer is not affected by the presence of the disaccharides, suggesting that the overall structure of the bilayer remains undisturbed. The results of simulations reveal that the interaction of disaccharide molecules with the bilayer occurs at the surface of the bilayer, and it is governed by the formation of multiple hydrogen bonds to specific groups of the lipid. Disaccharide molecules are observed to adopt specific conformations to fit onto the surface topology of the bilayer, often interacting with up to three different lipids simultaneously. At high disaccharide concentrations, the results of simulations indicate that disaccharides can serve as an effective replacement for water under anhydrous conditions, which helps explain their effectiveness as lyophilization agents for liposomes and cells.

Density-of-states Monte Carlo method for simulation of fluids

A Monte Carlo method based on a density-of-states sampling is proposed for study of arbitrary statistical mechanical ensembles in a continuum. A random walk in the two-dimensional space of particle number and energy is used to estimate the density of states of the system; this density of states is continuously updated as the random walk visits individual states. The validity and usefulness of the method are demonstrated by applying it to the simulation of a Lennard-Jones fluid. Results for its thermodynamic properties, including the vapor–liquid phase coexistence curve, are shown to be in good agreement with high-accuracy literature data.

Automatic parameterization of force fields for liquids by simplex optimization

In this study we demonstrate an automatic method of force field development for molecular simulations. Parameter tuning is taken as an optimization problem in many dimensions. The parameters are automatically adapted to reproduce known experimental data such as the density and the heat of vaporization. Our method is more systematic than guessing parameters and, at the same time, saves human labor in parameterization. It was applied successfully to several molecular liquids. As a test, force fields for 2‐methylpentane, tetrahydrofurane, cyclohexene, and cyclohexane were developed. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1009–1017, 1999

Coarse Graining of Nonbonded Inter-Particle Potentials Using Automatic Simplex Optimization to Fit Structural Properties

We implemented a coarse-graining procedure to construct mesoscopic models of complex molecules. The final aim is to obtain better results on properties depending on slow modes of the molecules. Therefore the number of particles considered in molecular dynamics simulations is reduced while conserving as many properties of the original substance as possible. We address the problem of finding nonbonded interaction parameters which reproduce structural properties from experiment or atomistic simulations. The approach consists of optimizing automatically nonbonded parameters using the simplex algorithm to fit structural properties like the radial distribution function as target functions. Moreover, any mix of structural and thermodynamic properties can be included in the target function. Different spherically symmetric inter-particle potentials are discussed. Besides demonstrating the method for Lennard-Jones liquids, it is applied to several more complex molecular liquids such as diphenyl carbonate, tetrahydrofurane, and monomers of poly(isoprene).

Automatic coarse graining of polymers

Several recently proposed semi-automatic and fully-automatic coarse-graining schemes for polymer simulations are discussed. All these techniques derive effective potentials for multi-atom units or super-atoms from atomistic simulations. These include techniques relying on single chain simulations in vacuum and self-consistent optimizations from the melt like the simplex method and the inverted Boltzmann method. The focus is on matching the polymer structure on different scales. Several ways to obtain a time-scale for dynamic mapping are discussed additionally. Finally, similarities to other simulation areas where automatic optimization are applied as well are pointed out.

Multicanonical parallel tempering

We present a novel implementation of the parallel tempering Monte Carlo method in a multicanonical ensemble. Multicanonical weights are derived by a self-consistent iterative process using a Boltzmann inversion of global energy histograms. This procedure gives rise to a much broader overlap of thermodynamic-property histograms; fewer replicas are necessary in parallel tempering simulations, and the acceptance of trial swap moves can be made arbitrarily high. We demonstrate the usefulness of the method in the context of a grand-multicanonical ensemble, where we use multicanonical simulations in energy space with the addition of an unmodified chemical potential term in particle-number space. Several possible implementations are discussed, and the best choice is presented in the context of the liquid–gas phase transition of the Lennard-Jones fluid. A substantial decrease in the necessary number of replicas can be achieved through the proposed method, thereby providing a higher efficiency and the possibility ...

Interactions of Lipid Bilayers with Supports: A Coarse-Grained Molecular Simulation Study

The study of lipid structure and phase behavior at the nanoscale is of utmost importance due to implications in understanding the role of the lipids in biochemical membrane processes. Supported lipid bilayers play a key role in understanding real biological systems, but they are vastly underrepresented in computational studies. In this paper, we discuss molecular dynamics simulations of supported lipid bilayers using a coarse-grained model. We first focus on the technical implications of modeling solid supports for biomembrane simulations. We then describe noticeable influences of the support on the systems. We are able to demonstrate that the bilayer system behavior changes when supported by a hydrophilic surface. We find that the thickness of the water layer between the support and the bilayer (the inner-water region in the latter part of this paper) adapts through water permeation on the microsecond time scale. Additionally, we discuss how different surface topologies affect the bilayer. Finally, we point out the differences between the two leaflets induced by the support.

Potential of mean force between a spherical particle suspended in a nematic liquid crystal and a substrate

The expanded ensemble density of states method (ExEDOS) is used to investigate the effective interaction of a spherical colloidal particle suspended in a confined liquid crystal (LC) with a substrate. The potential of mean force (PMF) is determined as a function of the normal distance between the particle and the substrates surface. The presence of the substrate induces a layered structure of the LC, which in turn greatly influences the PMF. We analyze the structure of the Saturn ring defect that accompanies the colloidal sphere, and find that the ring is displaced slightly towards the surface when the sphere is within the first LC surface layer. A transition occurs from an overall attraction of the colloid to the substrate to a global repulsion when the spheres radius is roughly twice the length of the LC molecules.

Coarse-grained modeling of lipids.

Molecular modeling of phospholipids on many scales has progressed significantly over the last years. Here we review several membrane models on intermediate to large length scales restricting ourselves to particle based coarse-grained models with implicit and explicit solvent. We explain similarities and differences as well as their connection to experiments and fine-grained models. We neglect any field descriptions on larger scales. We discuss then a few examples where we focus on studies of lipid phase behavior as well as supported lipid bilayers as these examples can only be meaningfully studied using large-scale models to date.

Examining the Contributions of Lipid Shape and Headgroup Charge on Bilayer Behavior

To better understand bilayer property dependency on lipid electrostatics and headgroup size, we use atomistic molecular dynamics simulations to study negatively charged and neutral lipid membranes. We compare the negatively charged phosphatidic acid (PA), which at physiological pH and salt concentration has a negative spontaneous curvature, with the negatively charged phosphatidylglycerol (PG) and neutrally charged phosphatidylcholine (PC), both of which have zero spontaneous curvature. The PA lipids are simulated using two different sets of partial charges for the headgroup and the varied charge distribution between the two PA systems results in significantly different locations for the Na(+) ions relative to the water/membrane interface. For one PA system, the Na(+) ions are localized around the phosphate group. In the second PA system, the Na(+) ions are located near the ester carbonyl atoms, which coincides with the preferred location site for the PG Na(+) ions. We find that the Na(+) ion location has a larger effect on bilayer fluidity properties than lipid headgroup size, where the A(lipid) and acyl chain order parameter values are more similar between the PA and PG bilayers that have Na(+) ions located near the ester groups than between the two PA bilayers.