W. F. Drew Bennett

Improved Parameters for the Martini Coarse-Grained Protein Force Field

The Martini coarse-grained force field has been successfully used for simulating a wide range of (bio)molecular systems. Recent progress in our ability to test the model against fully atomistic force fields, however, has revealed some shortcomings. Most notable, phenylalanine and proline were too hydrophobic, and dimers formed by polar residues in apolar solvents did not bind strongly enough. Here, we reparametrize these residues either through reassignment of particle types or by introducing embedded charges. The new parameters are tested with respect to partitioning across a lipid bilayer, membrane binding of Wimley-White peptides, and dimerization free energy in solvents of different polarity. In addition, we improve some of the bonded terms in the Martini protein force field that lead to a more realistic length of α-helices and to improved numerical stability for polyalanine and glycine repeats. The new parameter set is denoted Martini version 2.2.

Distribution of Amino Acids in a Lipid Bilayer from Computer Simulations

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.

Molecular View of Cholesterol Flip-Flop and Chemical Potential in Different Membrane Environments

The relative stability of cholesterol in cellular membranes and the thermodynamics of fluctuations from equilibrium have important consequences for sterol trafficking and lateral domain formation. We used molecular dynamics computer simulations to investigate the partitioning of cholesterol in a systematic set of lipid bilayers. In addition to atomistic simulations, we undertook a large set of coarse grained simulations, which allowed longer time and length scales to be sampled. Our results agree with recent experiments (Steck, T. L.; et al. Biophys. J. 2002, 83, 2118-2125) that the rate of cholesterol flip-flop can be fast on physiological time scales, while extending our understanding of this process to a range of lipids. We predicted that the rate of flip-flop is strongly dependent on the composition of the bilayer. In polyunsaturated bilayers, cholesterol undergoes flip-flop on a submicrosecond time scale, while flip-flop occurs in the second range in saturated bilayers with high cholesterol content. We also calculated the free energy of cholesterol desorption, which can be equated to the excess chemical potential of cholesterol in the bilayer compared to water. The free energy of cholesterol desorption from a DPPC bilayer is 80 kJ/mol, compared to 67 kJ/mol for a DAPC bilayer. In general, cholesterol prefers more ordered and rigid bilayers and has the lowest affinity for bilayers with two polyunsaturated chains. Overall, the simulations provide a detailed molecular level thermodynamic description of cholesterol interactions with lipid bilayers, of fundamental importance to eukaryotic life.

Statistical Convergence of Equilibrium Properties in Simulations of Molecular Solutes Embedded in Lipid Bilayers

In recent years, atomistic molecular simulations have become a method of choice for studying the interaction of small molecules, peptides, and proteins with biological membranes. Here, we critically examine the statistical convergence of equilibrium properties in molecular simulations of two amino acid side-chain analogs, leucine and arginine, in the presence of a hydrated phospholipid bilayer. To this end, the convergence of the standard binding free energy for the reversible insertion of the solutes in the bilayer is systematically assessed by evaluating dozens of separate sets of umbrella sampling calculations for a total simulation time exceeding 400 μs. We identify rare and abrupt transitions in bilayer structure as a function of solute insertion depth. These transitions correspond to the slow reorganization of ionic interactions involving zwitterionic phospholipid headgroups when the solutes penetrate the lipid-water interface and when arginine is forced through the bilayer center. These rare events are shown to constitute hidden sampling barriers that limit the rate of convergence of equilibrium properties and result in systematic sampling errors. Our analysis demonstrates that the difficulty of attaining convergence for lipid bilayer-embedded solutes has, in general, been drastically underestimated. This information will assist future studies in improving accuracy by selecting a more appropriate reaction coordinate or by focusing computational resources on those regions of the reaction coordinate that exhibit slow convergence of equilibrium properties.

Thermodynamic analysis of the effect of cholesterol on dipalmitoylphosphatidylcholine lipid membranes.

Cholesterol is an important component of eukaryotic cellular membranes. Despite extensive literature on the physiochemical effects of cholesterol on membranes, much remains unknown about the precise role of cholesterol and its molecular interactions in membranes. Regular thermal fluctuations of lipids normal to the plane of the membrane are biologically relevant for many processes, such as interactions with enzymes, elastic properties, and hydrophobic matching, while larger fluctuations are involved in vesicle budding and fusion, passive lipid flip-flop, and pore formation. Here we used molecular dynamics simulations to investigate the thermodynamic effect of the cholesterol concentration on dipalmitoylphosphatidylcholine (DPPC) bilayers. We calculated the potentials of mean force for DPPC partitioning in DPPC bilayers containing 20 and 40 mol % cholesterol. Increasing the cholesterol content increases the free energy barrier for transferring the headgroup of DPPC to the center of the bilayer and slows the rate of DPPC flip-flop by orders of magnitude. Cholesterol increases the order, thickness, and rigidity of the bilayers, which restricts bilayer deformations and prevents pore formation. While DPPC flip-flop is pore-mediated in a pure bilayer, we do not observe pores in the 20 and 40 mol % bilayers. Increasing the cholesterol concentration causes a decrease in the free energy to transfer DPPC from its equilibrium position into bulk waterindicating that DPPC prefers to be in cholesterol-free bilayers. We also observe a reduction in small fluctuations of DPPC normal to the bilayer as the cholesterol concentration is increased.

Computer simulations of lipid membrane domains

There is great diversity in the composition and structure of biological lipid membranes. We are beginning to appreciate the crucial role of lipids in many cellular processes, and characterize some of the lateral structures within membranes that could play a role in the activity of lipids. Simulations probe molecular level interactions between single molecules, which provide complementary information to experiments. Lipid membrane simulations have reached an exciting point, where the time and length scales of our simulations are approaching experimental resolutions and can be used to interpret experiments on increasingly complex model membranes. The focus of this review is on recent molecular simulations of domain formation in lipid bilayers. We highlight a number of recent examples where simulations are used in collaboration with experiments. We review recent simulation studies on lipid-lipid interactions related to domain formation and on lipid-protein interactions relevant for lipid raft function.

Structural arrangement of the transmission interface in the antigen ABC transport complex TAP

The transporter associated with antigen processing (TAP) represents a focal point in the immune recognition of virally or malignantly transformed cells by translocating proteasomal degradation products into the endoplasmic reticulum–lumen for loading of MHC class I molecules. Based on a number of experimental data and the homology to the bacterial ABC exporter Sav1866, we constructed a 3D structural model of the core TAP complex and used it to examine the interface between the transmembrane and nucleotide-binding domains (NBD) by cysteine-scanning and cross-linking approaches. Herein, we demonstrate the functional importance of the newly identified X-loop in the NBD in coupling substrate binding to downstream events in the transport cycle. We further verified domain swapping in a heterodimeric ABC half-transporter complex by cysteine cross-linking. Strikingly, either substrate binding or translocation can be blocked by cross-linking the X-loop to coupling helix 2 or 1, respectively. These results resolve the structural arrangement of the transmission interface and point to different functions of the cytosolic loops and coupling helices in substrate binding, signaling, and transport.

Transfer of Arginine into Lipid Bilayers Is Nonadditive

Computer simulations suggest that the translocation of arginine through the hydrocarbon core of a lipid membrane proceeds by the formation of a water-filled defect that keeps the arginine molecule hydrated even at the center of the bilayer. We show here that adding additional arginine molecules into one of these water defects causes only a small change in free energy. The barrier for transferring multiple arginines through the membrane is approximately the same as for a single arginine and may even be lower depending on the exact geometry of the system. We discuss these results in the context of arginine-rich peptides such as antimicrobial and cell-penetrating peptides.

Atomistic Simulations of Pore Formation and Closure in Lipid Bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

Conical Lipids in Flat Bilayers Induce Packing Defects Similar to that Induced by Positive Curvature

In biological membranes, changes in lipid composition or mechanical deformations produce defects in the geometrical arrangement of lipids, thus allowing the adsorption of certain peripheral proteins. Here, we perform molecular dynamics simulations on bilayers containing a cylindrical lipid (PC) and a conical lipid (DOG). Profiles of atomic density and lateral pressure across the bilayer show differences in the acyl chain region due to deeper partitioning of DOG compared to PC. However, such analyses are less informative for the interfacial region where peripheral proteins adsorb. To circumvent this limitation, we develop, to our knowledge, a new method of membrane surface analysis. This method allows the identification of chemical defects, where hydrocarbon chains are accessible to the solvent, and geometrical defects, i.e., voids deeper than the glycerol backbone. The size and number of both types of defects increase with the number of monounsaturated acyl chains in PC and with the introduction of DOG, although the defects do not colocalize with the conical lipid. Interestingly, the size and probability of the defects promoted by DOG resemble those induced by positive curvature, thus explaining why conical lipids and positive curvature can both drive the adsorption of peripheral proteins that use hydrophobic residues as membrane anchors.