Bottom-up fabrication of a proteasome-nanopore that unravels and processes single proteins

The precise assembly and engineering of molecular machines capable of handling biomolecules play crucial roles in most single-molecule methods. In this work we use components from all three domains of life to fabricate an integrated multiprotein complex that controls the unfolding and threading of individual proteins across a nanopore. This 900 kDa multicomponent device was made in two steps. First, we designed a stable and low-noise β-barrel nanopore sensor by linking the transmembrane region of bacterial protective antigen to a mammalian proteasome activator. An archaeal 20S proteasome was then built into the artificial nanopore to control the unfolding and linearized transport of proteins across the nanopore. This multicomponent molecular machine opens the door to two approaches in single-molecule protein analysis, in which selected substrate proteins are unfolded, fed to into the proteasomal chamber and then addressed either as fragmented peptides or intact polypeptides.

Thermostable D-amino acid decarboxylases derived from Thermotoga maritima diaminopimelate decarboxylase

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

Adaption to glucose limitation is modulated by the pleotropic regulator CcpA, independent of selection pressure strength

Background

A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized.

Results

We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments.

Conclusion

This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1331-x) contains supplementary material, which is available to authorized users.

cgHeliParm: analysis of dsDNA helical parameters for coarse-grained MARTINI molecular dynamics simulations

Summary

We introduce cgHeliParm, a python program that provides the conformational analysis of Martini-based coarse-grained double strand DNA molecules. The software calculates the helical parameters such as base, base pair and base pair step parameters. cgHeliParm can be used for the analysis of coarse grain Martini molecular dynamics trajectories without transformation into atomistic models.

Availability and implementation

This package works with Python 2.7 on MacOS and Linux. The program is freely available for download from https://github.com/ifaust83/cgheliparm. Together with the main script, the base reference files CG_X_std.lib, a number of examples and R scripts are also available from the same website. A tutorial on the use and application is also available at http://cgmartini.nl/index.php/tutorials-general-introduction/tutorial-martini-dna.

Contact

i.faustino@rug.nl.

Supplementary information

Supplementary data are available at Bioinformatics online.

Insight into the complete substrate-binding pocket of ThiT by chemical and genetic mutations

Energy-coupling factor (ECF) transporters are involved in the uptake of micronutrients in bacteria. The transporters capture the substrate by high-affinity binding proteins, the so-called S-components. Here, we present the analysis of two regions of the substrate-binding pocket of the thiamine-specific S-component in Lactococcus lactis, ThiT. First, interaction of the thiazolium ring of thiamine with residues Trp34, His125 and Glu84 by π-π-stacking and cation-π is studied, and second, the part of the binding pocket that extends from the hydroxyl group. We mutated either the transported ligand (chemically) or the protein (genetically). Surprisingly, modifications in the thiazolium ring by introducing substituents with opposite electronic effects had similar effects on the binding affinity. We hypothesize that the electronic effects are superseeded by steric effects of the added substituents, which renders the study of isolated interactions difficult. Amino acid substitutions in ThiT indicate that the electrostatic interaction facilitated by residue Glu84 of ThiT and thiamine is necessary for picomolar affinity. Deazathiamine derivatives that explore the subpocket of the binding site extending from the hydroxyl group of thiamine bind with high affinity to ThiT and may be developed into selective inhibitors of thiamine transport by ECF transporters. Molecular-dynamics simulations suggest that two of these derivatives may not only bind to ThiT, but could also be transported.

Parameters for Martini sterols and hopanoids based on a virtual-site description

Sterols play an essential role in modulating bilayer structure and dynamics. Coarse-grained molecular dynamics parameters for cholesterol and related molecules are available for the Martini force field and have been successfully used in multiple lipid bilayer studies. In this work, we focus on the use of virtual sites as a means of increasing the stability of cholesterol and cholesterol-like structures. We improve and extend the Martini parameterization of sterols in four different ways: 1-the cholesterol parameters were adapted to make use of virtual interaction sites, which markedly improves numerical stability; 2-cholesterol parameters were also modified to address reported shortcomings in reproducing correct lipid phase behavior in mixed membranes; 3-parameters for ergosterol were created and adapted from cholesterols; and 4-parameters for the hopanoid class of bacterial polycyclic molecules were created, namely, for hopane, diploptene, bacteriohopanetetrol, and for their polycyclic base structure.

Molecular mechanism of cardiolipin-mediated assembly of respiratory chain supercomplexes

Mitochondria produce most of the ATP consumed by cells through the respiratory chain in their inner membrane. This process involves protein complexes assembled into larger structures, the respiratory supercomplexes (SCs). Cardiolipin (CL), the mitochondrial signature phospholipid, is crucial for the structural and functional integrity of these SCs, but it is as yet unclear by what mechanism it operates. Our data disclose the mechanism for bulk CL in gluing SCs, steering their formation, and suggest how it may stabilize specific interfaces. We describe self-assembly molecular dynamics simulations of 9 cytochrome bc1 (CIII) dimers and 27 cytochrome c oxidase (CIV) monomers from bovine heart mitochondria embedded in a CL-containing model lipid bilayer, aimed at mimicking the crowdedness and complexity of mitochondrial membranes. The simulations reveal a large diversity of interfaces, including those of existing experimental CIII/CIV SC models and an alternative interface with CIV rotated by 180°. SC interfaces enclose 4 to 12 CLs, a ∼10 fold enrichment from the bulk. Half of these CLs glue complexes together using CL binding sites at the surface of both complexes. Free energy calculations demonstrate a larger CL binding strength, compared to other mitochondrial lipids, that is exclusive to these binding sites and results from non-additive electrostatic and van der Waals forces. This study provides a key example of the ability of lipids to selectively mediate protein-protein interactions by altering all ranges of forces, lubricate protein interfaces and act as traffic control agents steering proteins together.

Altered secondary structure of Dynorphin A associates with loss of opioid signalling and NMDA-mediated excitotoxicity in SCA23

Spinocerebellar ataxia type 23 (SCA23) is caused by missense mutations in prodynorphin, encoding the precursor protein for the opioid neuropeptides α-neoendorphin, Dynorphin (Dyn) A and Dyn B, leading to neurotoxic elevated mutant Dyn A levels. Dyn A acts on opioid receptors to reduce pain in the spinal cord, but its cerebellar function remains largely unknown. Increased concentration of or prolonged exposure to Dyn A is neurotoxic and these deleterious effects are very likely caused by an N-methyl-d-aspartate-mediated non-opioid mechanism as Dyn A peptides were shown to bind NMDA receptors and potentiate their glutamate-evoked currents. In the present study, we investigated the cellular mechanisms underlying SCA23-mutant Dyn A neurotoxicity. We show that SCA23 mutations in the Dyn A-coding region disrupted peptide secondary structure leading to a loss of the N-terminal α-helix associated with decreased κ-opioid receptor affinity. Additionally, the altered secondary structure led to increased peptide stability of R6W and R9C Dyn A, as these peptides showed marked degradation resistance, which coincided with decreased peptide solubility. Notably, L5S Dyn A displayed increased degradation and no aggregation. R6W and wt Dyn A peptides were most toxic to primary cerebellar neurons. For R6W Dyn A, this is likely because of a switch from opioid to NMDA- receptor signalling, while for wt Dyn A, this switch was not observed. We propose that the pathology of SCA23 results from converging mechanisms of loss of opioid-mediated neuroprotection and NMDA-mediated excitotoxicity.

Hsc70-4 Deforms Membranes to Promote Synaptic Protein Turnover by Endosomal Microautophagy

Synapses are often far from their cell bodies and must largely independently cope with dysfunctional proteins resulting from synaptic activity and stress. To identify membrane-associated machines that can engulf synaptic targets destined for degradation, we performed a large-scale in vitro liposome-based screen followed by functional studies. We identified a presynaptically enriched chaperone Hsc70-4 that bends membranes based on its ability to oligomerize. This activity promotes endosomal microautophagy and the turnover of specific synaptic proteins. Loss of microautophagy slows down neurotransmission while gain of microautophagy increases neurotransmission. Interestingly, Sgt, a cochaperone of Hsc70-4, is able to switch the activity of Hsc70-4 from synaptic endosomal microautophagy toward chaperone activity. Hence, Hsc70-4 controls rejuvenation of the synaptic protein pool in a dual way: either by refolding proteins together with Sgt, or by targeting them for degradation by facilitating endosomal microautophagy based on its membrane deforming activity.

Dual action of BPC194: a membrane active peptide killing bacterial cells

Membrane active peptides can perturb the lipid bilayer in several ways, such as poration and fusion of the target cell membrane, and thereby efficiently kill bacterial cells. We probe here the mechanistic basis of membrane poration and fusion caused by membrane-active, antimicrobial peptides. We show that the cyclic antimicrobial peptide, BPC194, inhibits growth of Gram-negative bacteria and ruptures the outer and inner membrane at the onset of killing, suggesting that not just poration is taking place at the cell envelope. To simplify the system and to better understand the mechanism of action, we performed Förster resonance energy transfer and cryogenic transmission electron microscopy studies in model membranes and show that the BPC194 causes fusion of vesicles. The fusogenic action is accompanied by leakage as probed by dual-color fluorescence burst analysis at a single liposome level. Atomistic molecular dynamics simulations reveal how the peptides are able to simultaneously perturb the membrane towards porated and fused states. We show that the cyclic antimicrobial peptides trigger both fusion and pore formation and that such large membrane perturbations have a similar mechanistic basis.

Efficient Algorithms for Langevin and DPD Dynamics

In this article, we present several algorithms for stochastic dynamics, including Langevin dynamics and different variants of Dissipative Particle Dynamics (DPD), applicable to systems with or without constraints. The algorithms are based on the impulsive application of friction and noise, thus avoiding the computational complexity of algorithms that apply continuous friction and noise. Simulation results on thermostat strength and diffusion properties for ideal gas, coarse-grained (MARTINI) water, and constrained atomic (SPC/E) water systems are discussed. We show that the measured thermal relaxation rates agree well with theoretical predictions. The influence of various parameters on the diffusion coefficient is discussed.

Rhodopsin forms a dimer with cytoplasmic helix 8 contacts in native membranes

G protein-coupled receptors form dimers and higher-order oligomers in membranes, but the precise mode of receptor-receptor interaction remains unknown. To probe the intradimeric proximity of helix 8 (H8), we conducted chemical cross-linking of endogenous cysteines in rhodopsin in disk membranes. We identified a Cys316-Cys316 cross-link using partial proteolysis and liquid chromatography with mass spectrometry. These results show that a symmetric dimer interface mediated by H1 and H8 contacts is present in native membranes.

Sphere-to-rod transitions of nonionic surfactant micelles in aqueous solution modeled by molecular dynamics simulations

Control of the size and agglomeration of micellar systems is important for pharmaceutical applications such as drug delivery. Although shape-related transitions in surfactant solutions are studied experimentally, their molecular mechanisms are still not well understood. In this study, we use coarse-grained molecular dynamics simulations to describe micellar assemblies of pentaethylene glycol monododecyl ether (C(12)E(5)) in aqueous solution at different concentrations. The obtained size and aggregation numbers of the aggregates formed are in very good agreement with the available experimental data. Importantly, increase of the concentration leads to a second critical micelle concentration where a transition to rod-like aggregates is observed. This transition is quantified in terms of shape anisotropy, together with a detailed structural analysis of the micelles as a function of aggregation number.

A tool for the morphological analysis of mixtures of lipids and water in computer simulations

When analyzing computer simulations of mixtures of lipids and water, the questions to be answered are often of a morphological nature. They can deal with global properties, like the kind of phase that is adopted or the presence or absence of certain key features like a pore or stalk, or with local properties, like the local curvature present at a particular part of the lipid/water interface. While in principle all of the information relating to the global and local morphological properties of a system can be obtained from the set of atomic coordinates generated by a computer simulation, the extraction of this information is a tedious task that usually involves using a visualization program and performing the analysis by eye. Here we present a tool that employs the technique of morphological image analysis (MIA) to automatically extract the global morphology—as given by Minkowski functionals—from a set of atomic coordinates, and creates an image of the system onto which the local curvatures are mapped as a color code.

Role of lipids in spheroidal high density lipoproteins

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.

Areas of monounsaturated diacylphosphatidylcholines

We have studied the structural properties of monounsaturated diacylphosphatidylcholine lipid bilayers (i.e., diCn:1PC, where n = 14, 16, 18, 20, 22, and 24 is the number of acyl chain carbons). High-resolution x-ray scattering data were analyzed in conjunction with contrast-varied neutron scattering data using a technique we recently developed. Analyses of the data show that the manner by which bilayer thickness increases with increasing n in monounsaturated diacylphosphatidylcholines is dependent on the double bond's position. For commonly available monounsaturated diacylphosphatidylcholines, this results in the nonlinear behavior of both bilayer thickness and lipid area, whereas for diC18:1PC bilayers, lipid area assumes a maximum value. It is worthwhile to note that compared to previous data, our results indicate that lipid areas are smaller by approximately 10%. This observation highlights the need to revisit lipid areas, as they are often used in comparisons with molecular dynamics simulations. Moreover, simulators are encouraged to compare their results not only to x-ray scattering data, but to neutron data as well.

A coarse-grained model for polyethylene oxide and polyethylene glycol: conformation and hydrodynamics

A coarse-grained (CG) model for polyethylene oxide (PEO) and polyethylene glycol (PEG) developed within the framework of the MARTINI CG force field (FF) using the distributions of bonds, angles, and dihedrals from the CHARMM all-atom FF is presented. Densities of neat low molecular weight PEO agree with experiment, and the radius of gyration R(g) = 19.1 A +/- 0.7 for 76-mers of PEO (M(w) approximately 3400), in excellent agreement with neutron scattering results for an equal sized PEG. Simulations of 9, 18, 27, 36, 44, 67, 76, 90, 112, 135, and 158-mers of the CG PEO (442 < M(w) < 6998) at low concentration in water show the experimentally observed transition from ideal chain to real chain behavior at 1600 < M(w) < 2000, in excellent agreement with the dependence of experimentally observed hydrodynamic radii of PEG. Hydrodynamic radii of PEO calculated from diffusion coefficients of the higher M(w) PEO also agree well with experiment. R(g) calculated from both all-atom and CG simulations of PEO76 at 21 and 148 mg/cm(3) are found to be nearly equal. This lack of concentration dependence implies that apparent R(g) from scattering experiments at high concentration should not be taken to be the chain dimension. Simulations of PEO grafted to a nonadsorbing surface yield a mushroom to brush transition that is well described by the Alexander-de Gennes formalism.

Combining an Elastic Network With a Coarse-Grained Molecular Force Field: Structure, Dynamics, and Intermolecular Recognition

Structure-based and physics-based coarse-grained molecular force fields have become attractive approaches to gain mechanistic insight into the function of large biomolecular assemblies. Here, we study how both approaches can be combined into a single representation, that we term ELNEDIN. In this representation an elastic network is used as a structural scaffold to describe and maintain the overall shape of a protein and a physics-based coarse-grained model (MARTINI-2.1) is used to describe both inter- and intramolecular interactions in the system. The results show that when used in molecular dynamics simulations ELNEDIN models can be built so that the resulting structural and dynamical properties of a protein, including its collective motions, are comparable to those obtained using atomistic protein models. We then evaluate the behavior of such models in (1) long, microsecond time-scale, simulations, (2) the modeling of very large macromolecular assemblies, a viral capsid, and (3) the study of a protein-protein association process, the reassembly of the ROP homodimer. The results for this series of tests indicate that ELNEDIN models allow microsecond time-scale molecular dynamics simulations to be carried out readily, that large biological entities such as the viral capsid of the cowpea mosaic virus can be stably modeled as assemblies of independent ELNEDIN models, and that ELNEDIN models show significant promise for modeling protein-protein association processes.

Alternative mechanisms for the interaction of the cell-penetrating peptides penetratin and the TAT peptide with lipid bilayers

Cell-penetrating peptides (CPPs) have recently attracted much interest due to their apparent ability to penetrate cell membranes in an energy-independent manner. Here molecular-dynamics simulation techniques were used to study the interaction of two CPPs: penetratin and the TAT peptide with 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospolipid bilayers shed light on alternative mechanisms by which these peptides might cross biological membranes. In contrast to previous simulation studies of charged peptides interacting with lipid bilayers, no spontaneous formation of transmembrane pores was observed. Instead, the simulations suggest that the peptides may enter the cell by micropinocytosis, whereby the peptides induce curvature in the membrane, ultimately leading to the formation of small vesicles within the cell that encapsulate the peptides. Specifically, multiple peptides were observed to induce large deformations in the lipid bilayer that persisted throughout the timescale of the simulations (hundreds of nanoseconds). Pore formation could be induced in simulations in which an external potential was used to pull a single penetratin or TAT peptide into the membrane. With the use of umbrella-sampling techniques, the free energy of inserting a single penetratin peptide into a DPPC bilayer was estimated to be approximately 75 kJmol(-1), which suggests that the spontaneous penetration of single peptides would require a timescale of at least seconds to minutes. This work also illustrates the extent to which the results of such simulations can depend on the initial conditions, the extent of equilibration, the size of the system, and the conditions under which the simulations are performed. The implications of this with respect to the current systems and to simulations of membrane-peptide interactions in general are discussed.

G Protein-Coupled Receptors Self-Assemble in Dynamics Simulations of Model Bilayers

Many integral membrane proteins assemble to form oligomeric structures in biological membranes. In particular, seven-transmembrane helical G protein-coupled receptors (GPCRs) appear to self-assemble constitutively in membranes, but the mechanism and physiological role of this assembly are unknown. We developed and employed coarse-grain molecular dynamics (CGMD) models to investigate the molecular basis of how the physicochemical properties of the phospholipid bilayer membrane affect self-assembly of visual rhodopsin, a prototypical GPCR. The CGMD method is a mesoscopic simulation technique in which groups of atoms are mapped to particles on the basis of a four-to-one rule. This systematic reduction of the degrees of freedom allows for computationally efficient calculation of the structure and dynamics of molecular assemblies for larger time and length scales than accessible to atomistic models, providing here an unprecedented view of spontaneous protein assembly in biomembranes. Systems with up to 16 rhodopsin molecules at a protein-to-lipid ratio of 1:100 were simulated for time scales of up to 8 micros. The results obtained for four different phospholipid environments showed that localized adaptation of the membrane bilayer to the presence of receptors is reproducibly most pronounced near transmembrane helices 2, 4, and 7. This local membrane deformation appears to be a key factor defining the rate, extent, and orientational preference of protein-protein association. The implications of our findings are discussed within a framework of a generalized mechanism of membrane protein self-assembly.

Molecular dynamics simulations of lipid vesicle fusion in atomic detail

The fusion of a membrane-bounded vesicle with a target membrane is a key step in intracellular trafficking, exocytosis, and drug delivery. Molecular dynamics simulations have been used to study the fusion of small unilamellar vesicles composed of a dipalmitoyl-phosphatidylcholine (DPPC)/palmitic acid 1:2 mixture in atomic detail. The simulations were performed at 350-370 K and mimicked the temperature- and pH-induced fusion of DPPC/palmitic acid vesicles from experiments by others. To make the calculations computationally feasible, a vesicle simulated at periodic boundary conditions was fused with its periodic image. Starting from a preformed stalk between the outer leaflets of the vesicle and its periodic image, a hemifused state formed within 2 ns. In one out of six simulations, a transient pore formed close to the stalk, resulting in the mixing of DPPC lipids between the outer and the inner leaflet. The hemifused state was (meta)stable on a timescale of up to 11 ns. Forcing a single lipid into the interior of the hemifusion diaphragm induced the formation and expansion of a fusion pore on a nanosecond timescale. This work opens the perspective to study a wide variety of mesoscopic biological processes in atomic detail.

Ion transport across transmembrane pores

To study the pore-mediated transport of ionic species across a lipid membrane, a series of molecular dynamics simulations have been performed of a dipalmitoyl-phosphatidyl-choline bilayer containing a preformed water pore in the presence of sodium and chloride ions. It is found that the stability of the transient water pores is greatly reduced in the presence of the ions. Specifically, the binding of sodium cations at the lipid/water interface increases the pore line tension, resulting in a destabilization of the pore. However, the application of mechanical stress opposes this effect. The flux of ions through these mechanically stabilized pores has been analyzed. Simulations indicate that the transport of the ions through the pores depends strongly on the size of the water channel. In the presence of small pores (radius <1.5 nm) permeation is slow, with both sodium and chloride permeating at similar rates. In the case in which the pores are larger (radius >1.5 nm), a crossover is observed to a regime where the anion flux is greatly enhanced. Based on these observations, a mechanism for the basal membrane permeability of ions is discussed.

Lipids Out of Equilibrium: Energetics of Desorption and Pore Mediated Flip-Flop

The potential of mean force (PMF) of a phospholipid in a bilayer is a key thermodynamic property that describes the energetic cost of localized lipid defects. We have calculated the PMF by umbrella sampling using molecular dynamics simulations. The profile has a deep minimum at the equilibrium position in the bilayer and steeply rises for displacements both deeper into the bilayer and moving away from the bilayer. As the lipid loses contact with the bilayer, the profile abruptly flattens without a significant barrier. The calculated free energy difference of 80 kJ/mol between the minimum of the PMF and the value in water agrees well with the free energy difference calculated from the experimentally measured critical micelle concentration. Significant water/lipid defects form when a lipid is forced into the bilayer interior, in the form of a small water pore that spans the membrane. The energy required to form such a water pore is also found to be 80 kJ/mol. On the basis of this energy, we estimate the lipid flip-flop rate and permeability rate of sodium ions. The resulting rates are in good agreement with experimental measurements, suggesting lipid flip-flop and basal permeability of ions are pore mediated.

Phase behavior of a phospholipid/fatty acid/water mixture studied in atomic detail

Molecular dynamics simulations have been used to study the phase behavior of a dipalmitoylphosphatidylcholine (DPPC)/palmitic acid (PA)/water 1:2:20 mixture in atomic detail. Starting from a random solution of DPPC and PA in water, the system adopts either a gel phase at temperatures below approximately 330 K or an inverted hexagonal phase above approximately 330 K in good agreement with experiment. It has also been possible to observe the direct transformation from a gel to an inverted hexagonal phase at elevated temperature (approximately 390 K). During this transformation, a metastable fluid lamellar intermediate is observed. Interlamellar connections or stalks form spontaneously on a nanosecond time scale and subsequently elongate, leading to the formation of an inverted hexagonal phase. This work opens the possibility of studying in detail how the formation of nonlamellar phases is affected by lipid composition and (fusion) peptides and, thus, is an important step toward understanding related biological processes, such as membrane fusion.

Simulation studies of pore and domain formation in a phospholipid monolayer

Despite extensive study the phase behavior of phospholipid monolayers at an air-water interface is still not fully understood. In particular recent vibrational sum-frequency generation (VSFG) spectra of DPPC monolayers as a function of area density show a sharp transition in the order of the lipid chains at 1.10 nm2/molecule. This is in a region where the lateral pressure as a function of area is effectively constant. We have investigated the nature of this transition by studying the phase behavior of DPPC monolayers as a function of area density using molecular-dynamics simulations. The changes in order within the monolayer as a function of area density correlate well with the experimental signal. At 0.58 nm2/molecule we observe the onset of lateral separation of highly ordered and disordered lipids, indicating the coexistence of a gel-like liquid condensed and a fluidlike liquid expanded phase. At 0.97 nm2/molecule the monolayer ruptures, marking the onset of the liquid-gas (G) coexistence region. This is much earlier than suggested by fluorescence microscopy results and implies that at the point of rupture, the initial pores have an equilibrium size smaller than approximately 500 nm in diameter. The rupture of the monolayer leads to a sharp increase in the overall lipid order that explains the sharp transition observed in the VSFG measurements. VSFG measurements thus may represent a sensitive means to determine the onset of the liquid-gas (G) coexistence region for such systems.

Molecular view of hexagonal phase formation in phospholipid membranes

Important biological processes, such as vesicle fusion or budding, require the cell matrix to undergo a transition from a lamellar to a nonlamellar state. Although equilibrium properties of membranes are amenable to detailed theoretical studies, collective rearrangements involved in phase transitions have thus far only been modeled on a qualitative level. Here, for the first time, the complete transition pathway from a multilamellar to an inverted hexagonal phase is elucidated at near-atomic detail using a recently developed coarse-grained molecular dynamics simulation model. Insight is provided into experimentally inaccessible data such as the molecular structure of the intermediates and the kinetics involved. Starting from multilamellar configurations, the spontaneous formation of stalks between the bilayers is observed on a nanosecond timescale at elevated temperatures or reduced hydration levels. The stalks subsequently elongate in a cooperative manner leading to the formation of an inverted hexagonal phase. The rate of stalk elongation is approximately 0.1 nm ns(-1). Within a narrow hydration/temperature/composition range the stalks appear stable and rearrange into the rhombohedral phase.

Simulation of domain formation in DLPC-DSPC mixed bilayers

Binary mixtures of two phosphatidylcholines of different chain lengths are simulated in the bilayer state. We find a phase transition between a liquid state and a gel state at all concentrations. This phase transition is characterized by the area per lipid headgroup, the order parameter, and a change in dynamics. At concentrations with a majority of the longer lipid, we find phase separation into a gel and a liquid state in a small temperature window. This leads to a strong dynamic heterogeneity. Experimental phase transition temperatures are reproduced semiquantitatively. We see a clear shift in the phase transition to higher temperatures with increasing concentration of the longer lipid.