Probing microscopic material properties inside simulated membranes through spatially resolved three-dimensional local pressure fields and surface tensions

Influence of membrane lipid composition on the structure and activity of γ-secretase

γ-Secretase (GS) is a multi-subunit membrane-embedded aspartyl protease that cleaves more than 80 integral membrane proteins, including the amyloid precursor protein (APP) to produce the amyloid-β (Aβ) peptide. Oligomerization and aggregation of the 42-amino acid length Aβ isoform in the brain has been associated with the development and progression of Alzheimer's disease (AD). Based on recent experimental structural studies and using multiscale computational modeling approaches, the conformational states and protein-membrane interactions of the GS complex embedded in six homogeneous and six heterogeneous lipid bilayers were characterized. In order to identify potential lipid and cholesterol binding sites, GS regions with high lipid/cholesterol occupancy values were analyzed using atomistic and coarse-grained simulations. Long lipid residence times were observed to be correlated with a large number of hydrogen bonds between the charged headgroups and key GS amino acids. This observation provides a plausible explanation for the inhibition of GS by charged lipids observed in previous experimental studies. Computed lateral pressure profiles suggest that higher transmembrane pressures favor active state conformations of the catalytic subunit. A probable mechanism for the regulation of the local stress response in cholesterol-rich multicomponent lipid bilayers is identified. Finally, it is demonstrated that interactions between the nicastrin extracellular domain and lipid headgroups leads to a compact structural conformation of the GS complex. Overall, this study provides valuable insight into the effect of bilayer lipid composition on the GS structural ensemble and its function.

Membrane stress profiles from self-consistent field theory

Using self-consistent field theory (SCFT), we develop an accurate, local expression for the stress profiles in membranes and soft matter interfaces, in general. The bond stresses are expressed in terms of pre-computed chain propagators, which are used to describe the statistical weight of the molecules and therefore require minimal additional calculations. In addition, we overcome the resolution limit of the molecular bond length by including the Irving and Kirkwood bond assignment and recover a constant normal stress profile across an interface. Using this theory, we find that the membrane lateral stress profile contains repulsive (positive) stresses in the regions of the head and tail groups, and attractive (negative) stresses near the hydrophobic/hydrophilic interface. We also verify that the zeroth and first moments of the stress profile correspond to the thermodynamic tension and product of the bending modulus and the spontaneous curvature, respectively.

Maximally asymmetric transbilayer distribution of anionic lipids alters the structure and interaction with lipids of an amyloidogenic protein dimer bound to the membrane surface

We used molecular dynamics simulations to explore the effects of asymmetric transbilayer distribution of anionic phosphatidylserine (PS) lipids on the structure of a protein on the membrane surface and subsequent protein-lipid interactions. Our simulation systems consisted of an amyloidogenic, beta-sheet rich dimeric protein (D42) absorbed to the phosphatidylcholine (PC) leaflet, or protein-contact PC leaflet, of two membrane systems: a single-component PC bilayer and double PC/PS bilayers. The latter comprised of a stable but asymmetric transbilayer distribution of PS in the presence of counterions, with a 1-component PC leaflet coupled to a 1-component PS leaflet in each bilayer. The maximally asymmetric PC/PS bilayer had a non-zero transmembrane potential (TMP) difference and higher lipid order packing, whereas the symmetric PC bilayer had a zero TMP difference and lower lipid order packing under physiologically relevant conditions. Analysis of the adsorbed protein structures revealed weaker protein binding, more folding in the N-terminal domain, more aggregation of the N- and C-terminal domains and larger tilt angle of D42 on the PC leaflet surface of the PC/PS bilayer versus the PC bilayer. Also, analysis of protein-induced membrane structural disruption revealed more localized bilayer thinning in the PC/PS versus PC bilayer. Although the electric field profile in the non-protein-contact PS leaflet of the PC/PS bilayer differed significantly from that in the non-protein-contact PC leaflet of the PC bilayer, no significant difference in the electric field profile in the protein-contact PC leaflet of either bilayer was evident. We speculate that lipid packing has a larger effect on the surface adsorbed protein structure than the electric field for a maximally asymmetric PC/PS bilayer. Our results support the mechanism that the higher lipid packing in a lipid leaflet promotes stronger protein-protein but weaker protein-lipid interactions for a dimeric protein on membrane surfaces.

Coupled diffusion in lipid bilayers upon close approach

Biomembrane interfaces create regions of slowed water dynamics in their vicinity. When two lipid bilayers come together, this effect is further accentuated, and the associated slowdown can affect the dynamics of larger-scale processes such as membrane fusion. We have used molecular dynamics simulations to examine how lipid and water dynamics are affected as two lipid bilayers approach each other. These two interacting fluid systems, lipid and water, both slow and become coupled when the lipid membranes are separated by a thin water layer. We show in particular that the water dynamics become glassy, and diffusion of lipids in the apposed leaflets becomes coupled across the water layer, while the "outer" leaflets remain unaffected. This dynamic coupling between bilayers appears mediated by lipid-water-lipid hydrogen bonding, as it occurs at bilayer separations where water-lipid hydrogen bonds become more common than water-water hydrogen bonds. We further show that such coupling occurs in simulations of vesicle-vesicle fusion prior to the fusion event itself. Such altered dynamics at membrane-membrane interfaces may both stabilize the interfacial contact and slow fusion stalk formation within the interface region.

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

Lipid tail protrusion in simulations predicts fusogenic activity of influenza fusion peptide mutants and conformational models

Fusion peptides from influenza hemagglutinin act on membranes to promote membrane fusion, but the mechanism by which they do so remains unknown. Recent theoretical work has suggested that contact of protruding lipid tails may be an important feature of the transition state for membrane fusion. If this is so, then influenza fusion peptides would be expected to promote tail protrusion in proportion to the ability of the corresponding full-length hemagglutinin to drive lipid mixing in fusion assays. We have performed molecular dynamics simulations of influenza fusion peptides in lipid bilayers, comparing the X-31 influenza strain against a series of N-terminal mutants. As hypothesized, the probability of lipid tail protrusion correlates well with the lipid mixing rate induced by each mutant. This supports the conclusion that tail protrusion is important to the transition state for fusion. Furthermore, it suggests that tail protrusion can be used to examine how fusion peptides might interact with membranes to promote fusion. Previous models for native influenza fusion peptide structure in membranes include a kinked helix, a straight helix, and a helical hairpin. Our simulations visit each of these conformations. Thus, the free energy differences between each are likely low enough that specifics of the membrane environment and peptide construct may be sufficient to modulate the equilibrium between them. However, the kinked helix promotes lipid tail protrusion in our simulations much more strongly than the other two structures. We therefore predict that the kinked helix is the most fusogenic of these three conformations.

Membrane fusion is a common process critical to both cellular function and infection by enveloped viruses. Influenza is a particularly useful model system for studying fusion because the fusion reaction is accomplished by a single protein, hemagglutinin. Furthermore, mutations to the membrane-inserted portion of hemagglutinin have been identified that do not detectably alter the rest of the protein but can either arrest fusion halfway or block it entirely. For influenza at least, it seems that the membrane-inserted hemagglutinin peptide plays a critical role in promoting fusion, perhaps by increasing the local disorder of lipid bilayers. However, we lack a mechanistic understanding sufficient to predict the activity of fusion peptide mutants from their sequence. Here, we have used lipid tail protrusion as a way to measure how much fusion peptides disorder their surrounding bilayer; we see a strong relationship between lipid tail protrusion and the ability of fusion peptide mutants to promote lipid mixing between membranes. Our simulations also predict that this lipid tail protrusion is much more common when the peptides adopt a kinked helix structure than when they are straight or hairpin-like. We therefore hypothesize that, while all three types of structure likely undergo conformational exchange, the kinked helix structure is most active in promoting fusion.