Molecular dynamics simulations reliably identify vibrational modes in far-IR spectra of phospholipids.

The properties of self-assembled phospholipid membranes are of essential importance in biochemistry and physical chemistry, providing a platform for many cellular life functions. Far-infrared (far-IR) vibrational spectroscopy, on the other hand, is a highly information-rich method to characterize intermolecular interactions and collective behaviour of lipids that can help explain, e.g., chain packing, thermodynamic phase behaviour, and sequestration. However, reliable interpretation of the far-IR spectra is still lacking. Here we present a molecular dynamics (MD) based approach to simulate vibrational modes of individual lipids and in an ensemble. The results are a good match to synchrotron far-IR measurements and enable identification of the molecular motions corresponding to each vibrational mode, thus allowing the correct interpretation of membrane spectra with high accuracy and resolving the longstanding ambiguities in the literature in this regard. Our results demonstrate the feasibility of using MD simulations for interpreting far-IR spectra broadly, opening new avenues for practical use of this powerful method.

Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes.

Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.

Modulation of lipid vesicle-membrane interactions by cholesterol.

Nanoscale lipid vesicles are attractive vehicles for drug delivery. Although often considered as soft nanoparticles in terms of mechanical deformability, the fluidic nature of the lipid membrane makes their interactions with another lipid membrane much more complex. Cholesterol is a key molecule that not only effectively stiffens lipid bilayer membranes but also induces membrane fusion. As such, how cholesterol modulates lipid vesicle-membrane interactions during endocytosis remains elusive. Through systematic molecular dynamics simulations, we find that membrane stiffening upon incorporating cholesterol reduces vesicle wrapping by a planar membrane, hindering endocytosis. Membrane fusion is also accelerated when either the vesicle or the planar membrane is cholesterol-rich, but fusion becomes minimal when both the vesicle and planar membrane are cholesterol-rich. This study provides insights into vesicle-membrane interactions in the presence of cholesterol and enlightens how cholesterol may be used to direct the cellular uptake pathways of nanoliposomes.

Ferroptosis induces membrane blebbing in placental trophoblasts.

Ferroptosis is a regulated, non-apoptotic form of cell death, characterized by hydroxy-peroxidation of discrete phospholipid hydroperoxides, particularly hydroperoxyl (Hp) forms of arachidonoyl- and adrenoyl-phosphatidylethanolamine, with a downstream cascade of oxidative damage to membrane lipids, proteins and DNA, culminating in cell death. We recently showed that human trophoblasts are particularly sensitive to ferroptosis caused by depletion or inhibition of glutathione peroxidase 4 (GPX4) or the lipase PLA2G6. Here, we show that trophoblastic ferroptosis is accompanied by a dramatic change in the trophoblast plasma membrane, with macro-blebbing and vesiculation. Immunofluorescence revealed that ferroptotic cell-derived blebs stained positive for F-actin, but negative for cytoplasmic organelle markers. Transfer of conditioned medium that contained detached macrovesicles or co-culture of wild-type target cells with blebbing cells did not stimulate ferroptosis in target cells. Molecular modeling showed that the presence of Hp-phosphatidylethanolamine in the cell membrane promoted its cell ability to be stretched. Together, our data establish that membrane macro-blebbing is characteristic of trophoblast ferroptosis and can serve as a useful marker of this process. Whether or not these blebs are physiologically functional remains to be established.

Chloroplast membrane lipid remodeling protects against dehydration by limiting membrane fusion and distortion.

Dehydration damages the structural integrity of the chloroplast membrane and, consequently, the normal photosynthetic function of this organelle. Remodeling of galactolipids by converting monogalactosyl-diacylglycerol (MGDG) to digalactosyl-diacylglycerol (DGDG) and oligo-galactolipids is an effective adaptation strategy for protecting against dehydration damage to the chloroplast membrane. However, detailed molecular mechanisms are missing. In this study, by performing molecular-level simulations of bi-lamellar membranes under various dehydration conditions, we find that MGDG-to-DGDG remodeling protects the chloroplast membrane in a unique manner by simultaneously dictating both the extent and the pattern of fusion stalks formed with the apposed membrane. Specifically, MGDG-rich membranes form elongated stalks at a moderate dehydration level, whereas DGDG-rich membranes form smaller, rounded stalks. Simulations of wild-type and mutant Arabidopsis (Arabidopsis thaliana) outer chloroplast membranes further confirm that the mutant membrane without galactolipid remodeling is more susceptible to membrane fusion due to its higher MGDG content. Our work reveals the underlying physical mechanisms that govern the pattern and extent of membrane fusion structures, paving the way for rational genetic engineering of crops with improved dehydration tolerance.

Potentiation of plant defense by bacterial outer membrane vesicles is mediated by membrane nanodomains.

Outer membrane vesicles (OMVs) are released from the outer membranes of Gram-negative bacteria during infection and modulate host immunity during host-pathogen interactions. The mechanisms by which OMVs are perceived by plants and affect host immunity are unclear. Here, we used the pathogen Xanthomonas campestris pv. campestris to demonstrate that OMV-plant interactions at the Arabidopsis thaliana plasma membrane (PM) modulate various host processes, including endocytosis, innate immune responses, and suppression of pathogenesis by phytobacteria. The lipid phase of OMVs is highly ordered and OMVs directly insert into the Arabidopsis PM, thereby enhancing the plant PM's lipid order; this also resulted in strengthened plant defenses. Strikingly, the integration of OMVs into the plant PM is host nanodomain- and remorin-dependent. Using coarse-grained simulations of molecular dynamics, we demonstrated that OMV integration into the plant PM depends on the membrane lipid order. Our computational simulations further showed that the saturation level of the OMV lipids could fine-tune the enhancement of host lipid order. Our work unraveled the mechanisms underlying the ability of OMVs produced by a plant pathogen to insert into the host PM, alter host membrane properties, and modulate plant immune responses.

Role of Membrane Stretch in Adsorption of Antiviral Peptides onto Lipid Membranes and Membrane Pore Formation.

Many medically important viruses are enveloped viruses, which are surrounded by a structurally conserved, host-derived lipid membrane coating. Agents that target and disrupt this membrane coating could potentially function as broad-spectrum antiviral drugs. The amphipathic α-helical (AH) peptide derived from the N-terminus of the hepatitis C virus NS5A protein is one such candidate and has been demonstrated to be able to selectively rupture lipid vesicles in the size range of viruses (<160 nm diameter). However, the mechanism underlying this membrane curvature selectivity remains elusive. In this study, we have performed molecular dynamics simulations to study the binding of the AH peptide to model membranes that are stretched to resemble the looser lipid headgroup packing present on highly curved outer membranes of nanoscale vesicles. We found that the AH peptide binds more favorably to membranes that are stretched. In addition, a tetrameric placement of peptides across the membrane induced stable pore formation in the stretched membrane. Thus, our results suggest that the AH peptide senses the high curvature of nanoscale vesicles via the enhanced exposure of lipid packing defects induced by membrane area strain.

Site-Specific Peroxidation Modulates Lipid Bilayer Mechanics.

Peroxidation of plasma membranes, characterized by oxidative attack of lipidic carbon-carbon double-bonds in unsaturated fatty acids, has been identified as an important biochemical event in multiple pathological conditions, including neurodegenerative diseases, atherosclerosis, diabetes, preeclampsia, aging, cancer, etc. Changes to the lipid bilayer structure as a result of lipid peroxidation may lead to lipid membrane malfunction, and consequently initiate further downstream biochemical cascades. However, how lipid peroxidation modulates the mechanical properties of lipid membranes remains largely controversial. In this study, we investigate the peroxidation of lipids with polyunsaturated fatty acid tails using molecular dynamics simulations. By systematically varying the oxidation site, we find that lipid peroxidation alters the biophysical properties of bilayer membrane in a peroxidation site-specific manner. Specifically, our results suggest that peroxidation at sites in the bilayer interior disturbs and softens the membrane, whereas peroxidation at sites near the membrane-water interface results in a more ordered and stiffer membrane. Such a peroxidation site-specific modulation of lipid membrane mechanics provides an explanation for the contradictory results obtained in previous experiments. Our study paves the way for an improved understanding of the initiation of the downstream cellular dysfunction caused by lipid peroxidation.

Curvature-Regulated Lipid Membrane Softening of Nano-Vesicles.

The physico-mechanical properties of nanoscale lipid vesicles (e.g., natural nano-vesicles and artificial nano-liposomes) dictate their interaction with biological systems. Understanding the interplay between vesicle size and stiffness is critical to both the understanding of the biological functions of natural nano-vesicles and the optimization of nano-vesicle-based diagnostics and therapeutics. It has been predicted that, when vesicle size is comparable to its membrane thickness, the effective bending stiffness of the vesicle increases dramatically due to both the entropic effect as a result of reduced thermal undulation and the nonlinear curvature elasticity effect. Through systematic molecular dynamics simulations, we show that the vesicle membrane thins and softens with the decrease in vesicle size, which effectively counteracts the stiffening effects as already mentioned. Our simulations indicate that the softening of nano-vesicles results from a change in the bilayer's interior structure - a decrease in lipid packing order - as the membrane curvature increases. Our work thus leads to a more complete physical framework to understand the physico-mechanical properties of nanoscale lipid vesicles, paving the way to further advances in the biophysics of nano-vesicles and their biomedical applications.

The binding interface of kindlin-2 and ILK involves Asp344/Asp352/Thr356 in kindlin-2 and Arg243/Arg334 in ILK.

Focal adhesion (FA) proteins, kindlin-2 and integrin-linked kinase (ILK), regulate cell adhesion and migration. ILK interacts with and promotes kindlin-2 targeting to FAs. Leu353 and Leu357 in kindlin-2 have been reported to be important for the interaction between kindlin-2 and ILK. However, the binding interface between kindlin-2 and ILK remains unclear. Using molecular modeling and molecular dynamics simulations, we show that Asp344, Asp352, and Thr356 in kindlin-2 and Arg243 and Arg334 in ILK kinase domain (KD) are important in kindlin-2/ILK complex formation. Mutations that disrupt these interactions abrogate kindlin-2 and ILK colocalization in HeLa cells. The interactions are direct based on data from pull-down assays using purified recombinant kindlin-2 F2-pleckstrin homology and ILK KDs. These data provide additional insights into the binding interface between kindlin-2 and ILK.

Lipid-associated aggregate formation of superoxide dismutase-1 is initiated by membrane-targeting loops.

Copper-Zinc superoxide dismutase 1 (SOD1) is a homodimeric enzyme that protects cells from oxidative damage. Hereditary and sporadic amyotrophic lateral sclerosis may be linked to SOD1 when the enzyme is destabilized through mutation or environmental stress. The cytotoxicity of demetallated or apo-SOD1 aggregates may be due to their ability to cause defects within cell membranes by co-aggregating with phospholipids. SOD1 monomers may associate with the inner cell membrane to receive copper ions from membrane-bound copper chaperones. But how apo-SOD1 interacts with lipids is unclear. We have used atomistic molecular dynamics simulations to reveal that flexible electrostatic and zinc-binding loops in apo-SOD1 dimers play a critical role in the binding of 1-octanol clusters and phospholipid bilayer, without any significant unfolding of the protein. The apo-SOD1 monomer also associates with phospholipid bilayer via its zinc-binding loop rather than its exposed hydrophobic dimerization interface. Our observed orientation of the monomer on the bilayer would facilitate its association with a membrane-bound copper chaperone. The orientation also suggests how membrane-bound monomers could act as seeds for membrane-associated SOD1 aggregation.

Extending neurites sense the depth of the underlying topography during neuronal differentiation and contact guidance.

The topography of the extracellular microenvironment influences cell morphology, provides conduct guidance and directs cell differentiation. Aspect ratio and dimension of topography have been shown to affect cell behaviours, but the ability and mechanism of depth-sensing is not clearly understood. We showed that murine neural progenitor cells (mNPCs) can sense the depth of the micro-gratings. Neurite elongation, alignment and neuronal differentiation were observed to increase with grating depth. We proposed a mechanism for depth-sensing by growing neurites: filopodial adhesion in the growth cones favour elongation but the bending rigidity of the neurite cytoskeleton resists it. Thus, perpendicular extension on deeper grooves is unfavourable as neurites need to bend over a larger angle. A quantitative model was developed and its prediction of neurite growth on gratings fit well with the experimental data. The results indicated that mNPC fate can be directed by appropriately designed patterned surfaces.

Leukocyte integrin αLß2 transmembrane association dynamics revealed by coarse-grained molecular dynamics simulations.

Integrins are transmembrane (TM) proteins that mediate bidirectional mechanical signaling between the extracellular matrix and the cellular cytoskeletal network. Each integrin molecule consists of non-covalently associated α- and ß-subunits, with each subunit consisting of a large ectodomain, a single-pass TM helix, and a short cytoplasmic tail. Previously we found evidence for a polar interaction (hydrogen bond) in the outer membrane clasp (OMC) of the leukocyte integrin αLß2 TMs that is absent in the platelet integrin αIIß3 OMC. Here, we compare the self-assembly dynamics of αLß2 and αIIß3 TM helices in a model membrane using coarse-grained molecular dynamics simulations. We found that although αIIß3 TM helices associate more easily, packing is suboptimal. In contrast, αLß2 TM helices achieve close-to-optimal packing. This suggests that αLß2 TM packing is more specific, possibly due to the interhelix hydrogen bond. Theoretical association free energy profiles show a deeper minimum at a smaller helix-helix separation for αLß2 compared with αIIß3. The αIIß3 profile is also more rugged with energetic barriers whereas that of αLß2 is almost without barriers. Disruption of the interhelix hydrogen bond in αLß2 via the ß2T686G mutation results in poorer association and a similar profile as αIIß3. The OMC polar interaction in αLß2 thus plays a significant role in the packing of the TM helices.

A transmembrane polar interaction is involved in the functional regulation of integrin alpha L beta 2.

Integrins are heterodimeric transmembrane (TM) receptors formed by noncovalent associations of alpha and beta subunits. Each subunit contains a single alpha-helical TM domain. Inside-out activation of an integrin involves the separation of its cytoplasmic tails, leading to disruption of alphabeta TM packing. The leukocyte integrin alpha L beta 2 is required for leukocyte adhesion, migration, proliferation, cytotoxic function, and antigen presentation. In this study, we show by mutagenesis experiments that the packing of alpha L beta 2 TMs is consistent with that of the integrin alpha IIb beta 3 TMs. However, molecular dynamics simulations of alpha L beta 2 TMs in lipids predicted a polar interaction involving the side chains of alpha L Ser1071 and beta2 Thr686 in the outer-membrane association clasp (OMC). This is supported by carbonyl vibrational shifts observed in isotope-labeled alpha L beta 2 TM peptides that were incorporated into lipid bilayers. Molecular dynamics studies simulating the separation of alpha L beta 2 tails showed the presence of polar interaction during the initial perturbation of the inner-membrane association clasp. When the TMs underwent further separation, the polar interaction was disrupted. OMC polar interaction is important in regulating the functions of beta2 integrins because mutations that disrupt the OMC polar interaction generated constitutively activated alpha L beta 2, alpha M beta 2, and alpha X beta 2 in 293T transfectants. We also show that the expression of mutant beta2 Thr686Gly in beta2-deficient T cells rescued cell adhesion to intercellular adhesion molecule 1, but the cells showed overt elongated morphologies in response to chemokine stromal-cell-derived factor 1 alpha treatment as compared to wild-type beta2-expressing cells. These two TM polar residues are totally conserved in other members of the beta2 integrins in humans and across different species. Our results provide an example of the stabilizing effect of polar interactions within the low dielectric environment of the membrane interior and demonstrate its importance in the regulation of alpha L beta 2 function.

Mechanical unfolding of bacterial flagellar filament protein by molecular dynamics simulation.

Bacterial flagellum is a nano-scale motility device constructed by self-assembly. During construction of the cell-exterior filament (the 'propeller'), subunit proteins (called flagellin) are thought to be exported through the hollow flagellum to the growing filament tip in an unfolded state. To gain insight into the unfolded state preceding any force-spectroscopy experiments on flagellin, we employed force-probe molecular dynamics simulations. Two schemes to attain an unfolded state suitable for efficient transport were examined: (i) stretching flagellin along its length; (ii) unzipping flagellin from its adjacently placed termini. Atomic-level unfolding pathways and the mechanical efforts involved under each scheme were obtained for the four-domain flagellin from S. typhimurium. Flagellin appeared stiffer and required larger unfolding forces when stretched as the relative sliding of beta-strands require the breaking of multiple hydrogen bonds at once. In contrast, unzipping requires lower unfolding forces as it mainly involves unraveling beta-sheets by breaking hydrogen bonds one by one.

Thermal unfolding simulations of bacterial flagellin: insight into its refolding before assembly.

Flagellin is the subunit of the bacterial filament, the micrometer-long propeller of a bacterial flagellum. The protein is believed to undergo unfolding for transport through the channel of the filament and to refold in a chamber at the end of the channel before being assembled into the growing filament. We report a thermal unfolding simulation study of S. typhimurium flagellin in aqueous solution as an attempt to gain atomic-level insight into the refolding process. Each molecule comprises two filament-core domains {D0, D1} and two hypervariable-region domains {D2, D3}. D2 can be separated into subdomains D2a and D2b. We observed a similar unfolding order of the domains as reported in experimental thermal denaturation. D2a and D3 exhibited high thermal stability and contained persistent three-stranded beta-sheets in the denatured state which could serve as folding cores to guide refolding. A recent mutagenesis study on flagellin stability seems to suggest the importance of the folding cores. Using crude size estimates, our data suggests that the chamber might be large enough for either denatured hypervariable-region domains or filament-core domains, but not whole flagellin; this implicates a two-staged refolding process.

Coarse-grained models reveal functional dynamics--I. Elastic network models--theories, comparisons and perspectives.

In this review, we summarize the progress on coarse-grained elastic network models (CG-ENMs) in the past decade. Theories were formulated to allow study of conformational dynamics in time/space frames of biological interest. Several highlighted models and their underlined hypotheses are introduced in physical depth. Important ENM offshoots, motivated to reproduce experimental data as well as to address the slow-mode-encoded configurational transitions, are also introduced. With the theoretical developments, computational cost is significantly reduced due to simplified potentials and coarse-grained schemes. Accumulating wealth of data suggest that ENMs agree equally well with experiment in describing equilibrium dynamics despite their distinct potentials and levels of coarse-graining. They however do differ in the slowest motional components that are essential to address large conformational changes of functional significance. The difference stems from the dissimilar curvatures of the harmonic energy wells described for each model. We also provide our views on the predictability of 'open to close' (open-->close) transitions of biomolecules on the basis of conformational selection theory. Lastly, we address the limitations of the ENM formalism which are partially alleviated by the complementary CG-MD approach, to be introduced in the second paper of this two-part series.

Coarse-grained models reveal functional dynamics--II. Molecular dynamics simulation at the coarse-grained level--theories and biological applications.

Molecular dynamics (MD) simulation has remained the most indispensable tool in studying equilibrium/non-equilibrium conformational dynamics since its advent 30 years ago. With advances in spectroscopy accompanying solved biocomplexes in growing sizes, sampling their dynamics that occur at biologically interesting spatial/temporal scales becomes computationally intractable; this motivated the use of coarse-grained (CG) approaches. CG-MD models are used to study folding and conformational transitions in reduced resolution and can employ enlarged time steps due to the absence of some of the fastest motions in the system. The Boltzmann-Inversion technique, heavily used in parameterizing these models, provides a smoothed-out effective potential on which molecular conformation evolves at a faster pace thus stretching simulations into tens of microseconds. As a result, a complete catalytic cycle of HIV-1 protease or the assembly of lipid-protein mixtures could be investigated by CG-MD to gain biological insights. In this review, we survey the theories developed in recent years, which are categorized into Folding-based and Molecular-Mechanics-based. In addition, physical bases in the selection of CG beads/time-step, the choice of effective potentials, representation of solvent, and restoration of molecular representations back to their atomic details are systematically discussed.

Hydration of a hydrophobic cavity and its functional role: a simulation study of human interleukin-1beta.

Molecular dynamics simulations reveal that the hydrophobic cavity in human cytokine Interleukin-1beta is hydrated and can dynamically accommodate between one and four water molecules. These waters have residence times >> 500 ps and can give rise to detectable NOEs, in agreement with NMR observations of Ernst et al. (Science 1995; 267:1813-1817). The waters also display high positional disorder within the cavity, which explains why they have not been resolved crystallographically. The average distribution of water molecules over time within the cavity matches well the low resolution electron density extracted by Yu et al. (Proc Natl Acad Sci 1999; 96:103-108). The water molecules hydrate the hydrophobic cavity preferentially as complex clusters. These clusters result from a combination of hydrogen bonds between the waters and stabilizing interactions between the waters and aromatic rings forming the cavity. Free energy estimates suggest that it takes 4-waters to hydrate the cavity in a thermodynamically stable manner leading to a gain in free energy of transfer from bulk of approximately approximately 3.6 kcal/mol. This arises from the existence of the water clusters in multiple hydrogen bonded states. In addition, the waters are found to migrate either individually or as clusters out of the cavity through several pathways. The upper limit for one-dimensional diffusion of the waters within the protein matrix is 4 A/ps (relative to 6 A/ps for bulk). Simulations reveal pathways in addition to those identified crystallographically, with motions controlled by the rotations of sidechains. We find that only when the hydrophobic cavity is hydrated, do correlated motions couple distant sites with the sites that make contact with the receptor and this data partly offers an explanation of experimental mutagenesis data. Simulations, together with recent observations based on mutagenesis by Heidary et al. (J Mol Biol 2005; 353:1187-1198) that hydrogen bond networks couple motions across long distances in interleukin-1beta, lead us to hypothesize that the hydration of the cavity (conserved across mammals) can thermodynamically enhance hydrogen bond networks to enable coupling across long distances by acting as a plug and this in turn enables a kinetic control of the rate of transmission of signals.