Effect of High Pressure on Fully Hydrated DPPC and POPC Bilayers

Molecular dynamics insights on temperature and pressure effects on electroporation

Electroporation is a cell-level phenomenon caused by an ionic imbalance in the membrane, being of great relevance in various fields of knowledge. A dependence of the pore formation kinetics on the environmental conditions (temperature and pressure) of the cell membrane has already been reported, but further clarification regarding how these variables affect the pore formation/resealing dynamics and the transport of molecules through the membrane is still lacking. The objective of the present study was to investigate the temperature (288-348 K) and pressure (1-5000 atm) effects on the electroporation kinetics using coarse-grained molecular dynamics simulations. Results shown that the time for pore formation and resealing increased with pressure and decreased with temperature, whereas the maximum pore radius increased with temperature and decreased with pressure. This behavior influenced the ion migration through the bilayer, and the higher ionic mobility was obtained in the 288 K/1000 atm simulations, i.e., a combination of low temperature and (not excessively) high pressure. These results were used to discuss some experimental observations regarding the extraction of intracellular compounds applying this technique. This study contributes to a better understanding of electroporation under different thermodynamic conditions and to an optimal selection of processing parameters in practical applications which exploit this phenomenon.

Can Urea and Trimethylamine-N-oxide Prevent the Pressure-Induced Phase Transition of Lipid Membrane?

Organisms dwelling in ocean trenches are exposed to the high hydrostatic pressure of ocean water. Increasing pressure can alter the membrane packing density and fluidity and trigger the fluid-to-gel phase transition. To combat environmental stress, the organisms synthesize small polar solutes, which are known as osmolytes. Urea and trimethylamine-N-oxide (TMAO) are two such solutes found in deep-sea creatures. While TMAO stabilizes protein, urea induces protein denaturation. These solutes strongly influence the packing density and membrane fluidity of the lipid bilayer at different conditions. But can these solutes affect the pressure-induced phase transition of the lipid membrane? In the present work, we have studied the effect of these two solutes on pressure-induced fluid-to-gel phase transition based on the all-atom molecular dynamics (MD) simulation approach. A high-pressure-stimulated fluid-to-gel phase transition of the membrane is seen at 800 bar, which is consistent with previous experiments. We have also observed that in the low-pressure region (1-400 bar), urea slightly increases the membrane fluidity where TMAO decreases the same. However, the phase transition pressure remains almost unchanged on the addition of urea while TMAO shifts the phase transition toward a lower pressure. We have found that the hydrogen (H)-bond interaction between lipid and urea plays an important role in preserving the fluidity of the membrane in the low-pressure zone. However, at a higher pressure, both water and urea are excluded from the membrane surface. TMAO is also excluded from the interfacial region of the membrane at all pressures. Exclusion from the membrane surface further triggers the phase transition of the lipid membrane from the fluid to gel phase at a high pressure.

Ripple-like instability in the simulated gel phase of finite size phosphocholine bilayers

Atomistic molecular dynamics simulations have reached a degree of maturity that makes it possible to investigate the lipid polymorphism of model bilayers over a wide range of temperatures. However if both the fluid L_α and tilted gel [Formula: see text] states are routinely obtained, the [Formula: see text] ripple phase of phosphatidylcholine lipid bilayers is still unsatifactorily described. Performing simulations of lipid bilayers made of different numbers of DPPC (1,2-dipalmitoylphosphatidylcholine) molecules ranging from 32 to 512, we demonstrate that the tilted gel phase [Formula: see text] expected below the pretransition cannot be obtained for large systems (equal or larger than 94 DPPC molecules) through common simulations settings or temperature treatments. Large systems are instead found in a disordered gel phase which display configurations, topography and energies reminiscent from the ripple phase [Formula: see text] observed between the pretransition and the main melting transition. We show how the state of the bilayers below the melting transition can be controlled and depends on thermal history and conditions of preparations. A mechanism for the observed topographic instability is suggested.

Structure of lateral heterogeneities in a coarse-grained model for multicomponent membranes

We study the lateral domain structure in a coarse-grained molecular model for multicomponent lipid bilayers by semi-grandcanonical Monte Carlo simulations. The membranes are filled with liquid ordered (lo) domains surrounded by a liquid disordered (ld) matrix. Depending on the membrane composition and temperature, we identify different morphological regimes: one regime (I) where the lo domains are small and relatively compact, and two regimes (II, II') where they are larger and often interconnected. In the latter two regimes, the ld matrix forms a network of disordered trenches separating the lo domains, with a relatively high content of interdigitated line defects. Since such defects are also a structural element of the modulated ripple phase in one component membranes, we argue that the regimes II, II' may be amorphous equivalents of the ripple phase in multicomponent membranes. We also analyze the local structure and provide evidence that the domains in regime I are stabilized by a monolayer curvature mechanism postulated in earlier work [S. Meinhardt et al., PNAS, 2013, 110, 4476].

Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane

We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.

The effect of H3O+ on the membrane morphology and hydrogen bonding of a phospholipid bilayer

At the 2017 meeting of the Australian Society for Biophysics, we presented the combined results from two recent studies showing how hydronium ions (H₃O⁺) modulate the structure and ion permeability of phospholipid bilayers. In the first study, the impact of H₃O⁺ on lipid packing had been identified using tethered bilayer lipid membranes in conjunction with electrical impedance spectroscopy and neutron reflectometry. The increased presence of H₃O⁺ (i.e. lower pH) led to a significant reduction in membrane conductivity and increased membrane thickness. A first-order explanation for the effect was assigned to alterations in the steric packing of the membrane lipids. Changes in packing were described by a critical packing parameter (CPP) related to the interfacial area and volume and shape of the membrane lipids. We proposed that increasing the concentraton of H₃O⁺ resulted in stronger hydrogen bonding between the phosphate oxygens at the water-lipid interface leading to a reduced area per lipid and slightly increased membrane thickness. At the meeting, a molecular model for these pH effects based on the result of our second study was presented. Multiple μs-long, unrestrained molecular dynamic (MD) simulations of a phosphatidylcholine lipid bilayer were carried out and showed a concentration dependent reduction in the area per lipid and an increase in bilayer thickness, in agreement with experimental data. Further, H₃O⁺ preferentially accumulated at the water-lipid interface, suggesting the localised pH at the membrane surface is much lower than the bulk bathing solution. Another significant finding was that the hydrogen bonds formed by H₃O⁺ ions with lipid headgroup oxygens are, on average, shorter in length and longer-lived than the ones formed in bulk water. In addition, the H₃O⁺ ions resided for longer periods in association with the carbonyl oxygens than with either phosphate oxygen in lipids. In summary, the MD simulations support a model where the hydrogen bonding capacity of H₃O⁺ for carbonyl and phosphate oxygens is the origin of the pH-induced changes in lipid packing in phospholipid membranes. These molecular-level studies are an important step towards a better understanding of the effect of pH on biological membranes.

The effect of hydronium ions on the structure of phospholipid membranes

This work studies the mechanisms by which hydronium ions modulate the structure of phospholipid bilayers.

Effects of High Pressure on Phospholipid Bilayers

The response of lipid membranes to changes in external pressure is important for many biological processes, and it can also be exploited for technological applications. In this work, we employ all-atom molecular dynamics simulations to characterize the changes in the physical properties of phospholipid bilayers brought about by high pressure (1000 bar). In particular, we study how the response differs, in relation to different chain unsaturation levels, by comparing monounsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and biunsaturated dioleoyl-phosphatidylcholine (DOPC) bilayers. Various structural, mechanical, and dynamical features are found to be altered by the pressure increase in both bilayers. Notably, for most properties, including bilayer area and thickness, lipid order parameters, lateral pressure profile, and curvature frustration energy, we observe significantly more pronounced effects for monounsaturated POPC than biunsaturated DOPC. Possible biological implications of the results obtained are discussed, especially in relation to how different lipids can control the structure and function of membrane proteins.

Effect of benzocaine and propranolol on phospholipid-based bilayers

Drug/bilayer interactions are fundamental in determining the action mechanism of active ingredients. Neutron techniques represent unique tools for having a clear comprehension of such interactions.

Some Like It Hot: The Effect of Sterols and Hopanoids on Lipid Ordering at High Temperature

Sterols and hopanoids have been suggested to reinforce membranes and protect against unfavorable environmental conditions. In particular, hopanoids are found in high concentrations in membranes of thermotolerant and thermophilic bacteria. However, the mechanism whereby sterols and hopanoids stabilize membranes at elevated temperatures is poorly understood. Here, the effect of temperature on the ordering of lipids in bilayers containing cholesterol or the hopanoids bacteriohopanetetrol and diplopterol was explored using molecular dynamics simulations. It is shown that cholesterol induces a high level of ordering over a wide range of temperatures. Bacteriohopanetetrol promotes order within the lipid tails but enhances fluid-like properties of the head groups at high temperatures. In contrast, diplopterol partitions in the midplane of the bilayer. This suggests that individual hopanoids fulfill distinct functions in membranes, with the ordering properties of bacteriohopanetetrol being particularly well suited to maintain the integrity of membranes at temperatures preferred by thermotolerant and thermophilic bacteria.

Long-timescale motions in glycerol-monopalmitate lipid bilayers investigated using molecular dynamics simulation

The occurrence of long-timescale motions in glycerol-1-monopalmitate (GMP) lipid bilayers is investigated based on previously reported 600 ns molecular dynamics simulations of a 2×8×8 GMP bilayer patch in the temperature range 302-338 K, performed at three different hydration levels, or in the presence of the cosolutes methanol or trehalose at three different concentrations. The types of long-timescale motions considered are: (i) the possible phase transitions; (ii) the precession of the relative collective tilt-angle of the two leaflets in the gel phase; (iii) the trans-gauche isomerization of the dihedral angles within the lipid aliphatic tails; and (iv) the flipping of single lipids across the two leaflets. The results provide a picture of GMP bilayers involving a rich spectrum of events occurring on a wide range of timescales, from the 100-ps range isomerization of single dihedral angles, via the 100-ns range of tilt precession motions, to the multi-μs range of phase transitions and lipid-flipping events.

Quantifying the diffusion of membrane proteins and peptides in black lipid membranes with 2-focus fluorescence correlation spectroscopy

Protein diffusion in lipid membranes is a key aspect of many cellular signaling processes. To quantitatively describe protein diffusion in membranes, several competing theoretical models have been proposed. Among these, the Saffman-Delbrück model is the most famous. This model predicts a logarithmic dependence of a protein's diffusion coefficient on its inverse hydrodynamic radius (D ∝ ln 1/R) for small radius values. For large radius values, it converges toward a D ∝ 1/R scaling. Recently, however, experimental data indicate a Stokes-Einstein-like behavior (D ∝ 1/R) of membrane protein diffusion at small protein radii. In this study, we investigate protein diffusion in black lipid membranes using dual-focus fluorescence correlation spectroscopy. This technique yields highly accurate diffusion coefficients for lipid and protein diffusion in membranes. We find that despite its simplicity, the Saffman-Delbrück model is able to describe protein diffusion extremely well and a Stokes-Einstein-like behavior can be ruled out.

Response to hydrostatic pressure of bicellar dispersions containing an anionic lipid: pressure-induced interdigitation

Bicellar dispersions of chain perdeuterated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC-d54), 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG), and 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC), with molar ratios of 3:1:1, were studied using variable-pressure (2)H NMR spectroscopy at hydrostatic pressures up to 125 MPa. Upon warming of the dispersions, spectra at ambient pressure indicated a progressive coalescence from small bilayered disks undergoing isotropic reorientation to more extended micellar structures in which spectra indicated anisotropic reorientation and, under some conditions, magnetic orientation and finally to randomly oriented lamellae or multilamellar vesicles. Temperatures for the onsets of anisotropic reorientation and random lamellar orientation increased with pressure at rates of 0.22 and 0.15 °C/MPa, respectively. In the 3.5-T magnetic field used for this work, magnetic orientation within the intermediate phase was not observed at 83 MPa or higher pressures. Comparison of spectra obtained at fixed pressure showed significant asymmetry between behaviors upon warming and cooling. For samples of DMPC-d54/DMPG/DHPC (3:1:1), but not DMPC-d54/DHPC (4:1), a persistent interdigitated phase was formed after repeated cooling from high temperature at 83 MPa. This is likely a metastable phase and might reflect kinetic trapping of the short-chain lipid component, DHPC, in a nonequilibrium spatial distribution as temperature is lowered at high pressure. Bicellar dispersions typically behave differently upon warming and cooling, and these observations could provide some insight into the observed behaviors in such systems. This work also suggests the possibility of trapping bicellar dispersions in persistent nonequilibrium morphologies.

The human transporter associated with antigen processing: molecular models to describe peptide binding competent states

The human transporter associated with antigen processing (TAP) is a member of the ATP binding cassette (ABC) transporter superfamily. TAP plays an essential role in the antigen presentation pathway by translocating cytosolic peptides derived from proteasomal degradation into the endoplasmic reticulum lumen. Here, the peptides are loaded into major histocompatibility class I molecules to be in turn exposed at the cell surface for recognition by T-cells. TAP is a heterodimer formed by the association of two half-transporters, TAP1 and TAP2, with a typical ABC transporter core that consists of two nucleotide binding domains and two transmembrane domains. Despite the availability of biological data, a full understanding of the mechanism of action of TAP is limited by the absence of experimental structures of the full-length transporter. Here, we present homology models of TAP built on the crystal structures of P-glycoprotein, ABCB10, and Sav1866. The models represent the transporter in inward- and outward-facing conformations that could represent initial and final states of the transport cycle, respectively. We described conserved regions in the endoplasmic reticulum-facing loops with a role in the opening and closing of the cavity. We also identified conserved π-stacking interactions in the cytosolic part of the transmembrane domains that could explain the experimental data available for TAP1-Phe-265. Electrostatic potential calculations gave structural insights into the role of residues involved in peptide binding, such as TAP1-Val-288, TAP2-Cys-213, TAP2-Met-218. Moreover, these calculations identified additional residues potentially involved in peptide binding, in turn verified with replica exchange simulations performed on a peptide bound to the inward-facing models.