Fritz Jähnig

Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature

Molecular dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 molecules of the lipid dipalmitoylphosphatidylcholine and 23 water molecules per lipid at an isotropic pressure of 1 atm and 50 degrees C. Special attention was devoted to reproduce the correct density of the lipid, because this quantity is known experimentally with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane molecules and varying the Lennard-Jones parameters until the experimental density and heat of vaporization were obtained. With these parameters the lipid density resulted in perfect agreement with the experimental density. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the experimental values, which, because of its correlation with the area per lipid, makes it possible to give a proper estimate of the area per lipid of 0.61 +/- 0.01 nm2.

Models for the structure of outer-membrane proteins of Escherichia coli derived from Raman spectroscopy and prediction methods

The secondary structure of porin, maltoporin and OmpA protein reconstituted in lipid membranes is determined by Raman spectroscopy. The three proteins have similar structures consisting of 50 to 60% beta-strand, about 20% beta-turn, and less than 15% alpha-helix. Employing a method for structural prediction that accounts for amphipathic beta-strands, folding models are developed for porin and for the segment of OmpA protein incorporated into the membrane. In the model, the OmpA fragment consists of eight amphipathic membrane-spanning beta-strands that form a beta-barrel. Similarly, porin is folded into ten amphipathic membrane-spanning beta-strands that are located at the surface of the trimer towards the lipids and eight predominantly hydrophilic strands in the interior.

The structure of melittin in membranes.

The conformation of the polypeptide melittin in lipid membranes as determined by Raman spectroscopy is a bent alpha-helix formed by the mainly hydrophobic residues 1-21, and a nonhelical COOH-terminal segment of the hydrophilic residues 22-26. Fluorescence quenching experiments on residue Trp19 reveal that all COOH-termini are located on that side of a vesicular membrane to which melittin was added. By means of fluorescence energy transfer between unmodified and modified Trp19 residues, melittin is shown to aggregate in membranes predominantly in the form of tetramers. These and previous results on the location and orientation of melittin permit the development of a model for the structure of melittin tetramers in membranes. The hydrophilic sides of four bilayer-spanning helices face each other to form a hydrophilic pore through the membrane.

Lipid vesicle adsorption versus formation of planar bilayers on solid surfaces

The absorption and spreading behavior of lipid vesicles composed of either palmitoyloleoylphosphatidylcholine (POPC) or Escherichia coli lipid upon contact with a glass surface was examined by fluorescence measurements. Fluorescently labeled lipids were used to determine 1) the amount of lipid adsorbed at the surface, 2) the extent of fusion of the vesicles upon contact with the surface, 3) the ability of the adsorbed lipids to undergo lateral diffusion, and 4) the accessibility of the adsorbed lipids by external water soluble molecules. The results of these measurements indicate that POPC vesicles spread on the surface and form a supported planar bilayer, whereas E. coli lipid vesicles adsorb to the surface and form a supported vesicle layer. Supported planar bilayers were found to be permeable for small molecules, whereas supported vesicles were impermeable and thus represented immobilized, topologically separate compartments.

Adhesion forces of lipids in a phospholipid membrane studied by molecular dynamics simulations

Lipid adhesion forces can be measured using several experimental techniques, but none of these techniques provide insight on the atomic level. Therefore, we performed extensive nonequilibrium molecular dynamics simulations of a phospholipid membrane in the liquid-crystalline phase out of which individual lipid molecules were pulled. In our method, as an idealization of the experimental setups, we have simply attached a harmonic spring to one of the lipid headgroup atoms. Upon retraction of the spring, the force needed to drag the lipid out of the membrane is recorded. By simulating different retraction rates, we were able to investigate the high pull rate part of the dynamical spectrum of lipid adhesion forces. We find that the adhesion force increases along the unbinding path, until the point of rupture is reached. The maximum value of the adhesion force, the rupture force, decreases as the pull rate becomes slower, and eventually enters a friction-dominated regime. The computed bond lengths depend on the rate of rupture, and show some scatter due to the nonequilibrium nature of the experiment. On average, the bond length increases from approximately 1.7 nm to 2.3 nm as the rates go down. Conformational analyses elucidate the detailed mechanism of lipid-membrane bond rupture. We present results of over 15 ns of membrane simulations. Implications for the interpretation and understanding of experimental rupture data are discussed.

Structure predictions of membrane proteins are not that bad

A simple generalization of the well-known hydrophobicity analysis is sufficient to predict membrane-spanning amphiphilic alpha-helices and beta-strands. The use of this method is illustrated for bacteriorhodopsin, OmpA protein and the photosynthetic reaction centre.

Proton transport across transient single-file water pores in a lipid membrane studied by molecular dynamics simulations

To test the hypothesis that water pores in a lipid membrane mediate the proton transport, molecular dynamic simulations of a phospholipid membrane, in which the formation of a water pore is induced, are reported. The probability density of such a pore in the membrane was obtained from the free energy of formation of the pore, which was computed from the average force needed to constrain the pore in the membrane. It was found that the free energy of a single file of water molecules spanning the bilayer is 108(+/-10) kJ/mol. From unconstrained molecular dynamic simulations it was further deduced that the nature of the pore is very transient, with a mean lifetime of a few picoseconds. The orientations of water molecules within the pore were also studied, and the spontaneous translocation of a turning defect was observed. The combined data allowed a permeability coefficient for proton permeation across the membrane to be computed, assuming that a suitable orientation of the water molecules in the pore allows protons to permeate the membrane relatively fast by means of a wirelike conductance mechanism. The computed value fits the experimental data only if it is assumed that the entry of the proton into the pore is not rate limiting.

Electrostatic free energy and shift of the phase transition for charged lipid membranes

For a charged membrane in an electrolyte solution the electrostatic free energy is derived treating the system as a diffuse double layer. The dependence of the energy on external parameters like surface charge density and temperature is obtained and the physical basis discussed. As an application the charges are shown to exert an electrostatic surface pressure on the lipid chain packing which leads to a shift in the phase transition of membranes. The results confirm the interpretation of experimental data as given by Träuble et al. in the accompanying paper.

The structure of the lactose permease derived from Raman spectroscopy and prediction methods

The secondary structure of the lactose permease of Escherichia coli reconstituted in lipid membranes was determined by Raman spectroscopy. The alpha‐helix content is approximately 70%, the beta‐strand content below 10% and beta‐turns contribute 15%. About 1/3 of the residues in alpha‐helices and most other residues are exposed to water. Employing a method for structural prediction which accounts for amphipathic helices, 10 membrane‐spanning helices are predicted which are either hydrophobic or amphipathic. They are expected to form an outer ring of helices in the membrane. The interior of the ring would be made of residues which are predominantly hydrophilic and, evoking the analogy to sugar‐binding proteins, suited to provide the sugar binding site.

The orientation of melittin in lipid membranes. A polarized infrared spectroscopy study

Polarized infrared spectra of melittin incorporated into macroscopically oriented lipid membranes are reported. From the linear dichroism of the amide I and amide II vibrational bands, the spatial orientation of the melittin helices was determined as being preferentially parallel to the membrane normal, under our experimental condition of low water content and an ordered lipid phase. Considering the various models for the orientation of melittin in lipid membranes proposed in the literature, we conclude that our data are in accord with an arrangement whereby the hydrophobic part of the polypeptide either spans the bilayer in the form of two bent helix segments, or is folded back within one monolayer in the form of a wedge.