Simulating POPC and POPC/POPG Bilayers: Conserved Packing and Altered Surface Reactivity

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Intrinsically disordered membrane anchors of Rheb, RhoA and DiRas3 small GTPases: Molecular dynamics, membrane organization, and interactions

Protein structure has been well established to play a key role in determining function; however, intrinsically disordered proteins and regions (IDPs and IDRs) defy this paradigm. IDPs and IDRs exist as an ensemble of structures rather than a stable 3D structure yet play essential roles in many cell signaling processes. Nearly all Ras Superfamily GTPases are tethered to membranes by a lipid tail at the end of a flexible IDR. The sequence of these IDRs are key determinants of membrane localization, and interactions between the IDR and the membrane have been shown to affect signaling in RAS proteins through modulation of dynamic membrane organization. Here we utilized atomistic molecular dynamics simulations to study the membrane interactions, conformational dynamics, and lipid sorting of three IDRs from small GTPases Rheb, RhoA and DiRas3 in model membranes representing their physiological target membranes. We found that complementarity between lipidated IDR sequence and target membrane lipid composition is a determinant of conformational plasticity. We also show that electrostatic interactions between anionic lipids and basic residues on IDRs generate semi-stable conformational sub-states, and a lack of these residues leads to greater conformational diversity. Finally, we show that small GTPase IDRs with a polybasic domain alter local lipid composition by segregating anionic membrane lipids, and, in some cases, excluding other lipids from their immediate proximity in favor of anionic lipids.

Development of the Integrated Computer Simulation Model of the Intracellular, Transmembrane, and Extracellular Domain of Platelet Integrin α IIb β 3 (Platelet Membrane Glycoprotein: GPIIb-IIIa)

Background  The structure and functions of the extracellular domain of platelet integrin α IIb β 3 (platelet membrane glycoprotein: GPIIb-IIIa) change substantially upon platelet activation. However, the stability of the integrated model of extracellular/transmembrane/intracellular domains of integrin α IIb β 3 with the inactive state of the extracellular domain has not been clarified. Methods  The integrated model of integrin α IIb β 3 was developed by combining the extracellular domain adopted from the crystal structure and the transmembrane and intracellular domain obtained by Nuclear Magnetic Resonace (NMR). The transmembrane domain was settled into the phosphatidylcholine (2-oleoyl-1-palmitoyl-sn-glycerol-3-phosphocholine (POPC)) lipid bilayer model. The position coordinates and velocity vectors of all atoms and water molecules around them were calculated by molecular dynamic (MD) simulation with the use of Chemistry at Harvard Macromolecular Mechanics force field in every 2 × 10 -15 seconds. Results  The root-mean-square deviations (RMSDs) of atoms constructing the integrated α IIb β 3 model apparently stabilized at approximately 23 Å after 200 ns of calculation. However, minor fluctuation persisted during the entire calculation period of 650 ns. The RMSDs of both α IIb and β 3 showed similar trends before 200 ns. The RMSD of β 3 apparently stabilized approximately at 15 Å at 400 ns with persisting minor fluctuation afterward, while the structural fluctuation in α IIb persisted throughout the 650 ns calculation period. Conclusion  In conclusion, the integrated model of the intracellular, transmembrane, and extracellular domain of integrin α IIb β 3 suggested persisting fluctuation even after convergence of MD calculation.

Effect of phospholipid liposomes on prion fragment (106-128) amyloid formation

Misfolding and aggregation of cellular prion protein (PrPc) is a major molecular process involved in the pathogenesis of prion diseases. Here, we studied the aggregation properties of a prion fragment peptide PrP(106-128). The results show that the peptide aggregates in a concentration-dependent manner in an aqueous solution and that the aggregation is sensitive to pH and the preformed amyloid seeds. Furthermore, we show that the zwitterionic POPC liposomes moderately inhibit the aggregation of PrP(106-128), whereas POPC/cholesterol (8:2) vesicles facilitate peptide aggregation likely due to the increase of the lipid packing order and membrane rigidity in the presence of cholesterol. In addition, anionic lipid vesicles of POPG and POPG/cholesterol above a certain concentration accelerate the aggregation of the peptide remarkably. The strong electrostatic interactions between the N-terminal region of the peptide and POPG may constrain the conformational plasticity of the peptide, preventing insertion of the peptide into the inner side of the membrane and thus promoting fibrillation on the membrane surface. The results suggest that the charge properties of the membrane, the composition of the liposomes, and the rigidity of lipid packing are critical in determining peptide adsorption on the membrane surface and the efficiency of the membrane in catalyzing peptide oligomeric nucleation and amyloid formation. The peptide could be used as an improved model molecule to investigate the mechanistic role of the crucial regions of PrP in aggregation in a membrane-rich environment and to screen effective inhibitors to block key interactions between these regions and membranes for preventing PrP aggregation.

Nucleoside Analog Reverse-Transcriptase Inhibitors in Membrane Environment: Molecular Dynamics Simulations

The behavior of four drugs from the family of nucleoside analog reverse-transcriptase inhibitors (zalcitabine, stavudine, didanosine, and apricitabine) in a membrane environment was traced using molecular dynamics simulations. The simulation models included bilayers and monolayers composed of POPC and POPG phospholipids. It was demonstrated that the drugs have a higher affinity towards POPG membranes than POPC membranes due to attractive long-range electrostatic interactions. The results obtained for monolayers were consistent with those obtained for bilayers. The drugs accumulated in the phospholipid polar headgroup region. Two adsorption modes were distinguished. They differed in the degree of penetration of the hydrophilic headgroup region. Hydrogen bonds between drug molecules and phospholipid heads were responsible for adsorption. It was shown that apricitabine penetrated the hydrophilic part of the POPC and POPG membranes more effectively than the other drugs. Van der Waals interactions between S atoms and lipids were responsible for this.

Lipid acyl chain protrusion induced by the influenza virus hemagglutinin fusion peptide detected by NMR paramagnetic relaxation enhancement

The glycoprotein spikes of membrane-enveloped viruses include a subunit that catalyzes fusion (joining) of the viral and target cell membranes. For influenza virus, this is subunit 2 of hemagglutinin which has a ∼ 20-residue N-terminal fusion peptide (Fp) region that binds target membrane. An outstanding question is whether there are associated membrane changes important for fusion. Several computational studies have found increased "protrusion" of lipid acyl chains near Fp, i.e. one or more chain carbons are closer to the aqueous region than the headgroup phosphorus. Protrusion may accelerate initial joining of outer leaflets of the two membranes into a stalk intermediate. In this study, higher protrusion probability in membrane with vs. without Fp is convincingly detected by larger Mn2+-associated increases in chain 13C NMR transverse relaxation rates (Γ2's). Data analysis provides a ratio Γ2,neighbor/Γ2,distant for lipids neighboring vs. more distant from the Fp. The calculated ratio depends on the number of Fp-neighboring lipids and the experimentally-derived range of 4 to 24 matches the range of increased protrusion probabilities from different simulations. For samples either with or without Fp, the Γ2 values are well-fitted by an exponential decay as the 13C site moves closer to the chain terminus. The decays correlate with free-energy of protrusion proportional to the number of protruded -CH2 groups, with free energy per -CH2 of ∼0.25 kBT. The NMR data support one major fusion role of the Fp to be much greater protrusion of lipid chains, with highest protrusion probability for chain regions closest to the headgroups.

Protein Corona of Anionic Fluid-Phase Liposomes Compromises Their Integrity Rather than Uptake by Cells

Despite the undisputable role of the protein corona in the biointeractions of liposome drug carriers, the field suffers from a lack of knowledge regarding the patterns of protein deposition on lipid surfaces with different compositions. Here, we investigated the protein coronas formed on liposomes of basic compositions containing combinations of egg phosphatidylcholine (PC), palmitoyloleoyl phosphatidylglycerol (POPG), and cholesterol. Liposome-protein complexes isolated by size-exclusion chromatography were delipidated and analyzed using label-free LC-MS/MS. The addition of the anionic lipid and cholesterol both affected the relative protein abundances (and not the total bound proteins) in the coronas. Highly anionic liposomes, namely those containing 40% POPG, carried corona enriched with cationic proteins (apolipoprotein C1, beta-2-glycoprotein 1, and cathelicidins) and were the least stable in the calcein release assay. Cholesterol improved the liposome stability in the plasma. However, the differences in the corona compositions had little effect on the liposome uptake by endothelial (EA.hy926) and phagocytic cells in the culture (U937) or ex vivo (blood-derived monocytes and neutrophils). The findings emphasize that the effect of protein corona on the performance of the liposomes as drug carriers occurs through compromising particle stability rather than interfering with cellular uptake.

Antimicrobial and Cell-Penetrating Peptides: Understanding Penetration for the Design of Novel Conjugate Antibiotics

Antimicrobial peptides (AMPs) are short oligopeptides that can penetrate the bacterial inner and outer membranes. Together with cell-penetrating peptides (CPPs), they are called membrane active peptides; peptides which can translocate across biological membranes. Over the last fifty years, attempts have been made to understand the molecular features that drive the interactions of membranes with membrane active peptides. This review examines the features of a membrane these peptides exploit for translocation, as well as the physicochemical characteristics of membrane active peptides which are important for translocation. Moreover, it presents examples of how these features have been used in recent years to create conjugates consisting of a membrane active peptide, called a "vector", attached to either a current or novel antibiotic, called a "cargo" or "payload". In addition, the review discusses what properties may contribute to an ideal peptide vector able to deliver cargoes across the bacterial outer membrane as the rising issue of antimicrobial resistance demands new strategies to be employed to combat this global public health threat.

Molecular Dynamics Simulation of Transport Mechanism of Graphene Quantum Dots Through Different Cell Membranes

Exploring the mechanisms underlying the permeation of graphene quantum dots (GQDs) through different cell membranes is key for the practical application of GQDs in medicine. Here, the permeation process of GQDs through different lipid membranes was evaluated using molecular dynamics (MD) simulations. Our results showed that GQDs can easily permeate into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE) lipid membranes with low phospholipid molecule densities but cannot permeate into 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE) lipid membranes with high phospholipid densities. Free energy calculation showed that a high-energy barrier exists on the surface of the POPE lipid membrane, which prevents GQDs from entering the cell membrane interior. Further analysis of the POPE membrane structure showed that sparsely arranged phospholipid molecules of the low-density lipid membrane facilitated the entry of GQDs into the interior of the membrane, compared to compactly arranged molecules in the high-density lipid membrane. Our simulation study provides new insights into the transmembrane transport of GQDs.

Salt effects on mixed composition membranes containing an antioxidant lipophilic edaravone derivative: a computational-experimental study

The protection of lipid membranes against oxidation avoids diseases associated with oxidative stress. As a strategy to contrast it, functionalized lipids with antioxidant activity are used to become part of membranes thus protecting them. For this purpose, a lipophilic edaravone derivative has been synthesized, adding a C18 saturated chain to the original structure. The antioxidant activity of C18-Edv has been demonstrated in our previous work. In this study, molecular dynamics simulations have been performed to define the effects of NaCl, MgCl₂, KCl, and CaCl₂ salts on a palmitoyl-oleoyl-sn-glycero-phosphocholine (POPC) lipid bilayer encapsulating C18-Edv. The results showed how different salts influence POPC lateral diffusion, and the movements of C18-Edv heads, which are antioxidant moieties, were correlated to the ability of C18-Edv molecules to protect membranes. MgCl₂ showed a negative impact leading to C18-Edv clusterization and membrane stretching, while KCl and NaCl showed a moderate influence on the mixed lipid membrane structure. CaCl₂ increased the exposure of the C18-Edv heads to the lipid-water interface, resulting in the salt with a higher propensity to guarantee protection against radicals in the aqueous phase. Finally, C18-Edv-POPC liposomes have been prepared following the simulation conditions, and then an experimental Oxygen Radical Absorbance Capacity (ORAC) assay has been performed to confirm the in silico predicted results.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

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.

Small molecules enhance the potency of natural antimicrobial peptides

The skin-associated microbiome plays an important role in general well-being and in a variety of treatable skin conditions. In this regard, endogenous antimicrobial peptides have both a direct and indirect role in determining the composition of the microbiota. We demonstrate here that certain small molecular species can amplify the antimicrobial potency of naturally occurring antimicrobial peptides. In this study, we have used niacinamide, a form of vitamin B3 naturally found in foods and widely used in cosmetic skincare products, and two of its structural analogs, to investigate their cooperativity with the human antimicrobial peptide LL37 on the bacterium Staphylococcus aureus. We observed a clear synergistic effect of niacinamide and, to some extent, N-methylnicotinamide, whereas isonicotinamide showed no significant cooperativity with LL37. Adaptively biased molecular dynamics simulations using simplified model membrane substrates and single peptides revealed that these molecules partition into the headgroup region of an anionic bilayer used to mimic the bacterial membrane. The simulated effects on the physical properties of the simulated model membrane are well correlated with experimental activity observed in real biological assays despite the simplicity of the model. In contrast, these molecules have little effect on zwitterionic bilayers that mimic a mammalian membrane. We conclude that niacinamide and N-methylnicotinamide can therefore potentiate the activity of host peptides by modulating the physical properties of the bacterial membrane, and to a lesser extent through direct interactions with the peptide. The level of cooperativity is strongly dependent on the detailed chemistry of the additive, suggesting an opportunity to fine-tune the behavior of host peptides.

Impact of Lipid Ratio on the Permeability of Mixed Phosphatidylcholine/Phosphatidylglycerol Membranes in the Presence of 1-Dodecyl-3-methylimidazolium Bromide Ionic Liquid

We have studied the impact of the lipid ratio on the membrane permeability of mixed phosphatidylcholine (POPC)/phosphatidylglycerol (POPG) membranes induced by 1-dodecyl-3-methylimidazolium bromide ([C₁₂MIM]⁺Br^-) ionic liquid by evaluating the role of affinity and architecture of the phospholipid bilayer. Nine different model membranes composed of negatively charged POPG and zwitterionic POPC lipids mixed in molar ratios of 9:1, 8:2, 7:3, 6:4, 5:5, 4:6, 3:7, 2:8, and 1:9 have been studied. The membrane permeability of each composition has been evaluated using fluorescence-based dye leakage assays. Despite having the highest membrane affinity, POPG-rich membranes doped with 10 and 20 mol % POPC are found to be the least permeable. ³¹P- and ²H-based solid-state NMR investigations reveal that the minor POPC component is homogeneously dispersed in the PG/PC (8:2) membrane. In contrast, the lipids seem to be segregated into POPG- and POPC-rich domains in the complementary PG/PC (2:8) composition. Although [C₁₂MIM]⁺ cations have a stronger interaction with the POPG component in the mixed membranes, their insertion has a limited impact on the overall structure and dynamics of the PG/PC (8:2) composition.

Clustering of Stearic Acids in Model Phospholipid Membranes Revealed by Double Electron-Electron Resonance

Free fatty acids play various important roles in biological membranes. Double electron-electron resonance spectroscopy (DEER, also known as PELDOR) of spin-labeled biomolecules is capable of studying magnetic dipole-dipole (d-d) interactions between spin labels at the nanoscale range of distances. Here, DEER is applied to study intermolecular d-d interactions between doxyl-spin-labeled stearic acids (DSA) in gel-phase phospholipid bilayers composed either of an equimolecular mixture of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1,2-dioleoyl-sn-glycero-3-phosphocholine or of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine. DEER data obtained for different DSA concentrations showed that DSA molecules at their concentration in the bilayer χ larger than 0.5 mol % are assembled into lateral lipid-mediated clusters, with a characteristic intermolecular distance of 2 nm. Some evidences were obtained indicating that clusters may consist of "subclusters", alternatively appearing in two opposite leaflets. Conventional electron paramagnetic resonance (EPR) spectra for the gel-phase bilayers showed that for χ larger than 2 mol % the molecules in the clusters stick together, forming oligomers. Room-temperature EPR spectra for the liquid-crystalline phase were found to change noticeably for χ larger than 0.5 mol %, which may indicate the clustering in a liquid-crystalline phase similar to that observed by DEER in the gel phase.

Membrane electrostatics sensed by tryptophan anchors in hydrophobic model peptides depends on non-aromatic interfacial amino acids: implications in hydrophobic mismatch

WALPs are synthetic α-helical membrane-spanning peptides that constitute a well-studied system for exploring hydrophobic mismatch. These peptides represent a simplified consensus motif for transmembrane domains of intrinsic membrane proteins due to their hydrophobic core of alternating leucine and alanine flanked by membrane-anchoring aromatic tryptophan residues. Although the modulation of mismatch responses in WALPs by tryptophan anchors has been reported earlier, there have been limited attempts to utilize the intrinsic tryptophan fluorescence of this class of peptides in mismatch sensors. We have previously shown, utilizing the red edge excitation shift (REES) approach, that interfacial WALP tryptophan residues in fluid phase bilayers experience a dynamically constrained membrane microenvironment. Interestingly, emerging reports suggest the involvement of non-aromatic interfacially localized residues in modulating local structure and dynamics in WALP analogs. In this backdrop, we have explored the effect of interfacial amino acids, such as lysine (in KWALPs) and glycine (in GWALPs), on the tryptophan microenvironment of WALP analogs in zwitterionic and negatively charged membranes. We show that interfacial tryptophans in KWALP and GWALP experience a more restricted microenvironment, as reflected in the substantial increase in magnitude of REES and apparent rotational correlation time, relative to those in WALP in zwitterionic membranes. Interestingly, in contrast to WALP, the tryptophan anchors in KWALP and GWALP appear insensitive to the presence of negatively charged lipids in the membrane. These results reveal a subtle interplay between non-aromatic flanking residues in transmembrane helices and negatively charged lipids at the membrane interface, which could modulate the membrane microenvironment experienced by interfacially localized tryptophan residues. Since interfacial tryptophans are known to influence mismatch responses in WALPs, our results highlight the possibility of utilizing the fluorescence signatures of tryptophans in membrane proteins or model peptides such as WALP as markers for assessing protein responses to hydrophobic mismatch. More importantly, these results constitute one of the first reports on the influence of lipid headgroup charge in fine-tuning hydrophobic mismatch in membrane bilayers, thereby enriching the existing framework of hydrophobic mismatch.

Hydrodynamic coupling for particle-based solvent-free membrane models

The great challenge with biological membrane systems is the wide range of scales involved, from nanometers and picoseconds for individual lipids to the micrometers and beyond millisecond for cellular signaling processes. While solvent-free coarse-grained membrane models are convenient for large-scale simulations and promising to provide insight into slow processes involving membranes, these models usually have unrealistic kinetics. One major obstacle is the lack of an equally convenient way of introducing hydrodynamic coupling without significantly increasing the computational cost of the model. To address this, we introduce a framework based on anisotropic Langevin dynamics, for which major in-plane and out-of-plane hydrodynamic effects are modeled via friction and diffusion tensors from analytical or semi-analytical solutions to Stokes hydrodynamic equations. Using this framework, in conjunction with our recently developed membrane model, we obtain accurate dispersion relations for planar membrane patches, both free-standing and in the vicinity of a wall. We briefly discuss how non-equilibrium dynamics is affected by hydrodynamic interactions. We also measure the surface viscosity of the model membrane and discuss the affecting dissipative mechanisms.

Influence of Surface Ligand Molecular Structure on Phospholipid Membrane Disruption by Cationic Nanoparticles

Cationic nanoparticles are known to interact with biological membranes and often cause serious membrane damage. Therefore, it is important to understand the molecular mechanism for such interactions and the factors that impact the degree of membrane damage. Previously, we have demonstrated that spatial distribution of molecular charge at cationic nanoparticle surfaces plays an important role in determining the cellular uptake and membrane damage of these nanoparticles. In this work, using diamond nanoparticles (DNPs) functionalized with five different amine-based surface ligands and small phospholipid unilamellar vesicles (SUVs), we further investigate how chemical features and conformational flexibility of surface ligands impact nanoparticle/membrane interactions. ³¹P-NMR T₂ relaxation measurements quantify the mobility changes in lipid dynamics upon exposing the SUVs to functional DNPs, and coarse-grained molecular dynamics simulations further elucidate molecular details for the different modes of DNP-SUV interactions depending on the surface ligands. Collectively, our results show that the length of the hydrophobic segment and conformational flexibility of surface ligands are two key factors that dictate the degree of membrane damage by the DNP, while the amount of surface charge alone is not predictive of the strength of interaction.

Extracellular mutation induces an allosteric effect across the membrane and hampers the activity of MRP1 (ABCC1)

Dynamic conformational changes play a major role in the function of proteins, including the ATP-Binding Cassette (ABC) transporters. Multidrug Resistance Protein 1 (MRP1) is an ABC exporter that protects cells from toxic molecules. Overexpression of MRP1 has been shown to confer Multidrug Resistance (MDR), a phenomenon in which cancer cells are capable to defend themselves against a broad variety of drugs. In this study, we used varied computational techniques to explore the unique F583A mutation that is known to essentially lock the transporter in a low-affinity solute binding state. We demonstrate how macro-scale conformational changes affect MRP1’s stability and dynamics, and how these changes correspond to micro-scale structural perturbations in helices 10–11 and the nucleotide-binding domains (NBDs) of the protein in regions known to be crucial for its ATPase activity. We demonstrate how a single substitution of an outward-facing aromatic amino acid causes a long-range allosteric effect that propagates across the membrane, ranging from the extracellular ECL5 loop to the cytoplasmic NBD2 over a distance of nearly 75 Å, leaving the protein in a non-functional state, and provide the putative allosteric pathway. The identified allosteric structural pathway is not only in agreement with experimental data but enhances our mechanical understanding of MRP1, thereby facilitating the rational design of chemosensitizers toward the success of chemotherapy treatments.

Tylopeptin B peptide antibiotic in lipid membranes at low concentrations: Self-assembling, mutual repulsion and localization

The medium-length peptide Tylopeptin B possesses activity against Gram-positive bacteria. It binds to bacterial membranes altering their mechanical properties and increasing their permeability. This action is commonly related with peptide self-assembling, resulting in the formation of membrane channels. Here, pulsed double electron-electron resonance (DEER) data for spin-labeled Tylopeptin B in palmitoyl-oleoyl-glycero-phosphocholine (POPC) model membrane reveal that peptide self-assembling starts at concentration as low as 0.1 mol%; above 0.2 mol% it attains a saturation-like dependence with a mean number of peptides in the cluster <n> = 3.3. Using the electron spin echo envelope modulation (ESEEM) technique, Tylopeptin B molecules are found to possess a planar orientation in the membrane. In the peptide concentration range between 0.1 and 0.2 mol%, DEER data show that the peptide clusters have tendency of mutual repulsion, with a circle of inaccessibility of radius around 20 nm. It may be proposed that within this radius the peptides destabilize the membrane, providing so the peptide antimicrobial activity. Exploiting spin-labeled stearic acids as a model for free fatty acids (FFA), we found that at concentrations of 0.1-0.2 mol% the peptide promotes formation of lipid-mediated FFA clusters; further increase in peptide concentration results in dissipation of these clusters.

Molecular dynamics simulations and experimental studies reveal differential permeability of withaferin-A and withanone across the model cell membrane

Poor bioavailability due to the inability to cross the cell membrane is one of the major reasons for the failure of a drug in clinical trials. We have used molecular dynamics simulations to predict the membrane permeability of natural drugs-withanolides (withaferin-A and withanone) that have similar structures but remarkably differ in their cytotoxicity. We found that whereas withaferin-A, could proficiently transverse through the model membrane, withanone showed weak permeability. The free energy profiles for the interaction of withanolides with the model bilayer membrane revealed that whereas the polar head group of the membrane caused high resistance for the passage of withanone, the interior of the membrane behaves similarly for both withanolides. The solvation analysis further revealed that the high solvation of terminal O5 oxygen of withaferin-A was the major driving force for its high permeability; it interacted with the phosphate group of the membrane that led to its smooth passage across the bilayer. The computational predictions were tested by raising and recruiting unique antibodies that react to withaferin-A and withanone. The time-lapsed analyses of control and treated cells demonstrated higher permeation of withaferin-A as compared to withanone. The concurrence between the computation and experimental results thus re-emphasised the use of computational methods for predicting permeability and hence bioavailability of natural drug compounds in the drug development process.

Effects of lipid heterogeneity on model human brain lipid membranes

Cell membranes naturally contain a heterogeneous lipid distribution. However, homogeneous bilayers are commonly preferred and utilised in computer simulations due to their relative simplicity, and the availability of lipid force field parameters. Recently, experimental lipidomics data for the human brain cell membranes under healthy and Alzheimer's disease (AD) conditions were investigated, since disruption to the lipid composition has been implicated in neurodegenerative disorders, including AD [R. B. Chan et al., J. Biol. Chem., 2012, 287, 2678-2688]. In order to observe the effects of lipid complexity on the various bilayer properties, molecular dynamics simulations were used to study four membranes with increasing heterogeneity: a pure POPC membrane, a POPC and cholesterol membrane in a 1 : 1 ratio (POPC-CHOL), and to our knowledge, the first realistic models of a healthy brain membrane and an Alzheimer's diseased brain membrane. Numerous structural, interfacial, and dynamical properties, including the area per lipid, interdigitation, dipole potential, and lateral diffusion of the two simple models, POPC and POPC-CHOL, were analysed and compared to those of the complex brain models consisting of 27 lipid components. As the membranes gain heterogeneity, a number of alterations were found in the structural and dynamical properties, and more significant differences were observed in the lateral diffusion. Additionally, we observed snorkeling behaviour of the lipid tails that may play a role in the permeation of small molecules across biological membranes. In this work, atomistic description of realistic brain membrane models is provided, which can add insight towards the permeability and transport pathways of small molecules across these membrane barriers.

Beyond Shielding: The Roles of Glycans in the SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 28,000,000 infections and 900,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viral fusion proteins, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of the glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans and on the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift toward the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of the SARS-CoV-2 S protein, which may be exploited in the therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Membrane-Sugar Interactions Probed by Low-Frequency Raman Spectroscopy: The Monolayer Adsorption Model

Small sugars are known to stabilize biological membranes under extreme conditions of freezing and desiccation. The proposed mechanisms of stabilization suggest membrane-sugar interactions to be either attractive or repulsive. To obtain new insight into the problem, we use a recently developed low-frequency Raman scattering approach which allows detecting membrane mechanical vibrations. For model membranes of palmitoyl-oleoyl-glycero-phosphocholine (POPC) hydrated in aqueous sucrose and trehalose solutions, we studied the Raman peak between 12 and 15 cm^-1 that is attributed to an eigenmode of the normal mechanical vibrations of a lipid monolayer. For both sugars, similar results were obtained. With an increase in sugar concentration in solution, the frequency position of the peak was found to decrease by ∼13% which was interpreted as a consequence of the membrane thickening due sugar monolayer adsorption on the membrane surface. The concentration dependence of the peak frequency position was satisfactorily described by a Langmuir monolayer adsorption model. It is concluded that, at small sugar concentrations (less than 0.2 M), the membrane-sugar interactions are attractive, while at higher concentrations (more than 0.4 M) the attraction disappears. The data obtained show that one sugar molecule on the surface interacts with approximately 3-4 polar lipid heads.

Cryostabilization of the Cell Membrane of a Psychrotolerant Bacteria via Homeoviscous Adaptation

Homeoviscous adaptation (maintenance of a critical balance between the saturated and unsaturated lipids) of the cell membrane of psychrotolerant bacteria is essential to protect them against freeze-thaw cycle. But how does the homeoviscous adaptation protect the cell membrane during cold stress? In this Letter, we answer this question using a coarse-grained molecular dynamics simulation technique. On the basis of the reported fatty acid profiles of psychrotolerant bacteria at different temperatures, multiple lipid membranes are simulated at a wide range of temperatures between 250 and 300 K. We explicate how the homeoviscous adaptation minimizes the effect of cold stress on the structure and fluidity of the membrane. Partial freezing of the saturated lipid domain occurs with the self-aggregation of saturated and unsaturated lipids near the melting temperature of the unadapted lipid membrane. The gel-like phase provides necessary local packing density that can be sensed by sensor proteins responsible for the homeoviscous adaptation.

Beyond Shielding: The Roles of Glycans in SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 15,000,000 infections and 600,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates the host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viruses, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans, and the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift towards the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of SARS-CoV-2 S protein, which may be exploited by therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

The influence of membrane bilayer thickness on KcsA channel activity

Atomic resolution structures have provided significant insight into the gating and permeation mechanisms of various ion channels, including potassium channels. However, ion channels may also be regulated by numerous factors, including the physiochemical properties of the membrane in which they are embedded. For example, the matching of the bilayer's hydrophobic region to the hydrophobic external surface of the ion channel is thought to minimize the energetic penalty needed to solvate hydrophobic residues or exposed lipid tails. To understand the molecular basis of such regulation by hydrophobic matching requires examining channels in the presence of the lipid membrane. Here we examine the role of hydrophobic matching in regulating the activity of the model potassium channel, KcsA. ⁸⁶Rb⁺ influx assays and single-channel recordings indicate that the non-inactivating E71A KcsA channel is most active in thin bilayers (<diC18:1PC). Bilayer thickness affects the open probability of KcsA and not its unitary conductance. Molecular dynamics simulations indicate that the bilayer can sufficiently modify its dimensions to accommodate KcsA channels without major perturbations in the protein helical packing within the nanosecond timescale. Based on experimental results and MD simulations, we present a model in which bilayer thickness influences the stability of the open and closed conformations of the intracellular gate of KcsA, with minimal impact on the stability of the selectivity filter of the non-inactivating mutant, E71A.

Highlighting the effect of amyloid beta assemblies on the mechanical properties and conformational stability of cell membrane

Alzheimer disease (AD) is the most common cause of dementia, characterized by a progressive decline in cognitive function due to the abnormal aggregation and deposition of Amyloid beta (Aβ) fibrils in the brain of patients. In this context, the molecular mechanisms of protein misfolding and aggregation that are known to induce significant biophysical alterations in cells, including destabilization of plasma membranes, remain partially unclear. Physical interaction between the Aβ assemblies and the membrane leads to the disruption of the cell membrane in multiple ways including, surface carpeting, generation of transmembrane channels and detergent-like membrane dissolution. Understanding the impact of amyloidogenic protein in different stages of aggregation with the plasma membrane, plays a crucial role to fully elucidate the pathological mechanisms of AD. Within this framework, computer simulations represent a powerful tool able to shed lights on the interactions governing the structural influence of Aβ proteins on biological membrane. In this study, molecular dynamics (MD) simulations have been performed in order to characterize how POPC bilayer conformational and mechanical properties are affected by the interaction with Aβ_11-42 peptide, oligomer and fibril.

Complexity in Lipid Membrane Composition Induces Resilience to Aβ42 Aggregation

The molecular origins of Alzheimer's disease are associated with the aggregation of the amyloid-β peptide (Aβ). This process is controlled by a complex cellular homeostasis system, which involves a variety of components, including proteins, metabolites, and lipids. It has been shown in particular that certain components of lipid membranes can speed up Aβ aggregation. This observation prompts the question of whether there are protective cellular mechanisms to counterbalance this effect. Here, to address this issue, we investigate the role of the composition of lipid membranes in modulating the aggregation process of Aβ. By adopting a chemical kinetics approach, we first identify a panel of lipids that affect the aggregation of the 42-residue form of Aβ (Aβ₄₂), ranging from enhancement to inhibition. We then show that these effects tend to average out in mixtures of these lipids, as such mixtures buffer extreme aggregation behaviors as the number of components increases. These results indicate that a degree of quality control on protein aggregation can be achieved through a mechanism by which an increase in the molecular complexity of lipid membranes balances opposite effects and creates resilience to aggregation.

Collective absorption of 2,4,6-trinitrotoluene into lipid membranes and its effects on bilayer properties. A computational study

The widely used explosive, 2,4,6-trinitrotoluene (TNT), is a highly toxic chemical, which can cause hepatitis, cataracts, jaundice and so on, in humans. The interaction between TNT and biological membranes is crucial for understanding its toxic effects. Here, we mainly focused on molecular-level mechanisms for the collective adsorption of TNT into lipid membranes and the corresponding effects on bilayer properties by all-atom molecular dynamics simulations. We revealed that TNT can readily form an aggregate in the aqueous phase and quickly approach the surface of the membrane. At low concentrations of TNT (7 mol%), the aggregate is unstable and breaks up after several nanoseconds, and then the dispersed TNT molecules enter the membrane alone. At high concentrations (14 mol%), the aggregate is adsorbed as a whole and remains stable inside the membrane. After some of the TNT is absorbed by the membrane, the remaining TNT across the membrane would have greater permeability, i.e., the calculated permeability coefficient (P) is increased from 1.7 × 10^-2 to 18.3 cm s^-1. Correspondingly, a higher bioconcentration factor (BCF) was also observed. The increased level is more pronounced in the presence of TNT aggregates (i.e., high concentrations). This phenomenon is closely related to the strong interaction between TNT molecules. The results suggested that TNT molecules that have entered into the membrane can facilitate the membrane uptake, permeation and bioaccumulation of subsequent TNT molecules, exhibiting a synergistic effect. This work has a certain significance for understanding the toxicity of TNT.

Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation

Membrane reorientation of oncogenic RAS proteins is emerging as an important modulator of their functions. Previous studies have shown that the most common orientations include those with either the three C-terminal α-helices (OS1) or N-terminal β-strands (OS2) of the catalytic domain facing the membrane. OS1 and OS2 differ by the degree to which the effector-interacting surface is occluded by the membrane. However, the relative stability of these states and the rates of transition between them remained undetermined. How mutations might modulate preferences for specific orientation states is also far from clear. The current work attempted to address these questions through a comprehensive analysis of two 20 μs-long atomistic molecular dynamics simulations. The simulations were conducted on the oncogenic G12D and Q61H KRAS mutants bound to an anionic lipid bilayer. G12D and Q61H are among the most prevalent cancer-causing mutations at the P-loop and switch 2 regions of KRAS, respectively. We found that both mutants fluctuate in a similar manner between OS1 and OS2 via an intermediate orientation OS0, and both favor the signaling competent OS1 and OS0 over the occluded OS2. However, they differ in the details, such as in the extent to which they sample OS1. Analysis of the orientation free-energy landscapes estimated from the simulations indicate that OS1 and OS2 are the most stable states. However, the overall free energy surface is rugged, indicating a large diversity of conformations including at least two substates in each orientation state that differ in stability only by about 0.5-1.0 kcal/mol. Reversible transitions between OS1 and OS2 occur via two well-defined pathways that traverse OS0. In the minimum energy path, helix 4 remains close to the membrane as the angle of the catalytic domain from the membrane plane changes, resulting in a barrier of ∼1 kcal/mol for OS1/OS2 interconversions. Estimation of the rates of the various transitions based on survival probabilities yielded two rate constants in the order of 10⁷ and 10⁶ s^-1, which we attribute to intrinsic protein conformational dynamics and transient protein-lipid interactions, respectively. The faster process dominates every transition, confirming a previous suggestion that RAS membrane reorientation is driven by conformational fluctuations rather than protein-lipid interactions.

Mechanisms by Which Lipids Influence Conformational Dynamics of the GlpG Intramembrane Protease

Rhomboid intramembrane proteases are bound to lipid membranes, where they dock and cleave other transmembrane substrates. How the lipid membrane surrounding the protease impacts the conformational dynamics of the protease is essential to understand because it informs on the reaction coordinate of substrate binding. Atomistic molecular dynamics simulations allow us to probe protein motions and characterize the coupling between protein and lipids. Simulations performed here on GlpG, the rhomboid protease from Escherichia coli, indicate that the thickness of the lipid membrane close to GlpG depends on both the composition of the lipid membrane and the conformation of GlpG. Transient binding of a lipid headgroup at the active site of the protease, as observed in some of the simulations reported here, suggests that a lipid headgroup might compete with the substrate for access to the GlpG active site. Interactions identified between lipid headgroups and the protein influence the dynamics of lipid interactions close to the substrate-binding site. These observations suggest that the lipid membrane environment shapes the energy profile of the substrate-docking region of the enzyme reaction coordinate.

Elucidating the Lipid Binding Properties of Membrane-Active Peptides Using Cyclised Nanodiscs

The lipid composition of the cellular membrane plays an important role in a number of biological processes including the binding of membrane-active peptides. Characterization of membrane binding remains challenging, due to the technical limitations associated with the use of standard biophysical techniques and available membrane models. Here, we investigate the lipid binding properties of two membrane-active peptides, VSTx1, a well characterized ion-channel inhibitor, identified from spider venom, that preferentially binds to anionic lipid mixtures, and AA139 an antimicrobial β-hairpin peptide with uncharacterised lipid binding properties, currently in pre-clinical development. The lipid binding properties of these peptides are elucidated using nanodiscs formed by both linear and circularized (sortase-mediated) forms of a membrane scaffold protein (MSP1D1ΔH5). We find that nanodiscs formed by circularized MSPs—in contrast to those formed by linear MSPs—are sufficiently stable under sample conditions typically used for biophysical measurements (including lipid composition, a range of buffers, temperatures and concentrations). Using these circularized nanodiscs, we are able to extract detailed thermodynamic data using isothermal titration calorimetry (ITC) as well as atomic resolution mapping of the lipid binding interfaces of our isotope labeled peptides using solution-state, heteronuclear, nuclear magnetic resonance (NMR) spectroscopy. This represents a novel and general approach for elucidating the thermodynamics and molecular interface of membrane-active peptides toward flat lipid bilayers of variable composition. Our approach is validated by first determining the thermodynamic parameters and binding interface of VSTx1 toward the lipid bilayer, which shows good agreement with previous studies using lipid micelles and liposomes. The method is then applied to AA139, where the membrane binding properties are unknown. This characterization, involved solving the high-resolution structure of AA139 in solution using NMR spectroscopy and the development of a suitable expression system for isotope labeling. AA139 was found to bind exclusively to anionic membranes with moderate affinity (K_d~low μM), and was found to have a lipid binding interface involving the termini of the β-hairpin structure. The preference of AA139 for anionic lipids supports a role for membrane binding in the mode-of-action of this peptide, which is also consistent with its higher inhibitory activity against bacterial cells compared to mammalian cells. The described approach is a powerful method for investigation of the membrane binding properties of this important class of molecules.

Superior Compatibility of C2 N with Human Red Blood Cell Membranes and the Underlying Mechanism

The widespread use of nanomaterials, such as carbon based 2D nanomaterials, in biomedical applications, has been accompanied by a growing concern on their biocompatibility, and in particular, on how they may affect the integrity of cell membranes. Herein, the interactions between C₂ N, a novel 2D nanomaterial, and human red blood cell membranes are explored using a combined experimental and theoretical approach. The experimental microscopies show that C₂ N exerts a negligible hemolysis effect on the blood cells with a superior compatibility to their cell membranes, when compared with the control system, reduced graphene oxide (rGO), which is found to be highly hemolytic. The molecular dynamics simulations further reveal the underlying molecular mechanisms, which indicate that C₂ N prefers to be adsorbed flat on the water-membrane interface. Interaction energy analyses demonstrate the crucial role of Coulombic contributions, originating from the unique electrostatic potential surface of C₂ N, in preventing C₂ N from penetrating into cell membranes. These findings indicate a high compatibility of C₂ N with cell membranes, which may provide useful foundation for the future exploration of this 2D nanomaterial in related biomedical applications.

The Impact of O-Glycosylation on Cyanidin Interaction with POPC Membranes: Structure-Activity Relationship

Cyanidin and its O-glycosides have many important physiological functions in plants and beneficial effects on human health. Their biological activity is not entirely clear and depends on the structure of the molecule, in particular, on the number and type of sugar substituents. Therefore, in this study the detailed structure-activity relationship (SARs) of the anthocyanins/anthocyanidins in relation to their interactions with lipid bilayer was determined. On the basis of their antioxidant activity and the changes induced by them in size and Zeta potential of lipid vesicles, and mobility and order of lipid acyl chains, the impact of the number and type of sugar substituents on the biological activity of the compounds was evaluated. The obtained results have shown, that 3-O-glycosylation changes the interaction of cyanidin with lipid bilayer entirely. The 3-O-glycosides containing a monosaccharide induces greater changes in physical properties of the lipid membrane than those containing disaccharides. The presence of additional sugar significantly reduces glycoside interaction with model lipid membrane. Furthermore, O-glycosylation alters the ability of cyanidin to scavenge free radicals. This alteration depends on the type of free radicals and the sensitivity of the method used for their determination.

Membrane potential and dynamics in a ternary lipid mixture: insights from molecular dynamics simulations

Transmembrane potential (Vm) plays critical roles in cell signaling and other functions. However, the impact of Vm on the structure and dynamics of membrane lipids and proteins, which are critical for the regulation of signaling, is still an open question. All-atom molecular dynamics (MD) simulation is emerging as a useful technique to address this issue. Previous atomistic MD simulations of pure or binary model membranes indicated that both ion imbalance and electric field can be used to generate Vm, but both approaches failed to yield structural changes in lipids with statistical significance. We hypothesized that a possible reason for this could be oversimplified membrane composition or limited sampling. In this work, we tested if and how Vm modulates the structure and dynamics of lipids in a physiologically relevant model membrane. Using a detailed side-by-side comparison, we first show that while both ion imbalance and electric field generate Vm in our complex membranes, only the latter could produce physiologically relevant Vm. We further show that double bonds in lipid acyl chains have a relatively large sensitivity to Vm. A single-bilayer model with an electric field showed the highest sensitivity in simulations under the isothermal-isobaric (NPT) ensemble, reproducing expected responses of head-group dipoles to Vm and suggesting that this approach may be more suitable for studying the structural effects of Vm. Our findings also shed light on the relationship between the macroscopic Vm and its atomic-level underpinnings.

Particle-based membrane model for mesoscopic simulation of cellular dynamics

We present a simple and computationally efficient coarse-grained and solvent-free model for simulating lipid bilayer membranes. In order to be used in concert with particle-based reaction-diffusion simulations, the model is purely based on interacting and reacting particles, each representing a coarse patch of a lipid monolayer. Particle interactions include nearest-neighbor bond-stretching and angle-bending and are parameterized so as to reproduce the local membrane mechanics given by the Helfrich energy density over a range of relevant curvatures. In-plane fluidity is implemented with Monte Carlo bond-flipping moves. The physical accuracy of the model is verified by five tests: (i) Power spectrum analysis of equilibrium thermal undulations is used to verify that the particle-based representation correctly captures the dynamics predicted by the continuum model of fluid membranes. (ii) It is verified that the input bending stiffness, against which the potential parameters are optimized, is accurately recovered. (iii) Isothermal area compressibility modulus of the membrane is calculated and is shown to be tunable to reproduce available values for different lipid bilayers, independent of the bending rigidity. (iv) Simulation of two-dimensional shear flow under a gravity force is employed to measure the effective in-plane viscosity of the membrane model and show the possibility of modeling membranes with specified viscosities. (v) Interaction of the bilayer membrane with a spherical nanoparticle is modeled as a test case for large membrane deformations and budding involved in cellular processes such as endocytosis. The results are shown to coincide well with the predicted behavior of continuum models, and the membrane model successfully mimics the expected budding behavior. We expect our model to be of high practical usability for ultra coarse-grained molecular dynamics or particle-based reaction-diffusion simulations of biological systems.

Faster Simulations with a 5 fs Time Step for Lipids in the CHARMM Force Field

The performance of all-atom molecular dynamics simulations is limited by an integration time step of 2 fs, which is needed to resolve the fastest degrees of freedom in the system, namely, the vibration of bonds and angles involving hydrogen atoms. The virtual interaction sites (VIS) method replaces hydrogen atoms by massless virtual interaction sites to eliminate these degrees of freedom while keeping intact nonbonded interactions and the explicit treatment of hydrogen atoms. We have modified the existing VIS algorithm for most lipids in the popular CHARMM36 force field by increasing the hydrogen atom masses at regular intervals in the lipid acyl chains and obtained lipid properties and pore formation free energies in very good agreement with those calculated in simulations without VIS. Our modified VIS scheme enables a 5 fs time step resulting in a significant performance gain for all-atom simulations of membranes. The method has the potential to make longer time and length scales accessible in all-atom simulations of membrane-protein complexes.

Effect of 5-trans Isomer of Arachidonic Acid on Model Liposomal Membranes Studied by a Combined Simulation and Experimental Approach

Unsaturated fatty acids are found in humans predominantly in the cis configuration. Fatty acids in the trans configuration are primarily the result of human processing (trans fats), but can also be formed endogenously by radical stress. The cis-trans isomerization of fatty acids by free radicals could be connected to several pathologies. Trans fats have been linked to an increased risk of coronary artery disease; however, the reasons for the resulting pathogenesis remain unclear. Here, we investigate the effect of a mono-trans isomer of arachidonic acid (C20:4-5trans, 8cis, 11cis, 14cis) produced by free radicals in physiological concentration on a model erythrocyte membrane using a combined experimental and theoretical approach. Molecular Dynamics (MD) simulations of two model lipid bilayers containing arachidonic acid and its 5-trans isomer in 3 mol% were carried out for this purpose. The 5-trans isomer formation in the phospholipids was catalyzed by HOCH₂CH₂S· radicals, generated from the corresponding thiol by γ-irradiation, in multilamellar vesicles of SAPC. Large unilamellar vesicles were made by the extrusion method (LUVET) as a biomimetic model for cis-trans isomerization. Atomic Force Microscopy and Dynamic Light Scattering were used to measure the average size, morphology, and the z-potential of the liposomes. Both results from MD simulations and experiments are in agreement and indicate that the two model membranes display different physicochemical properties in that the bilayers containing the trans fatty acids were more ordered and more rigid than those containing solely the cis arachidonic acid. Correspondingly, the average size of the liposomes containing trans isomers was smaller than the ones without.

Effect of ergosterol on the fungal membrane properties. All-atom and coarse-grained molecular dynamics study

Cell membranes are complex multicomponent systems consisting of thousands of different lipids with numerous embedded membrane proteins and many types of sterols. We used all-atom and coarse-grained molecular dynamics simulations to study the structural and dynamical properties of phospholipid bilayers containing four types of phospholipids and different amount of ergosterol, main sterol component in the fungal membranes. To characterize the influence of ergosterol on the membrane properties we analyzed the surface area per lipid, bilayer thickness, area compressibility modulus, mass density profiles, deuterium order parameters, and lateral diffusion coefficients. The presence of ergosterol induces the ordering of lipids leading to their denser packing, to reducing the lateral diffusion of lipids and lipid surface area, to increasing the thickness of bilayer and compressibility modulus. In addition, we evaluated each calculated property between the two simulation methods.

Molecular Dynamics Simulations of Ceramide and Ceramide-Phosphatidylcholine Bilayers

Recent studies in lipid raft formation and stratum corneum permeability have focused on the role of ceramides (CER). In this study, we use the all-atom CHARMM36 (C36) force field to simulate bilayers using N-palmitoylsphingosine (CER16) or α-hydroxy-N-stearoyl phytosphingosine (CER[AP]) in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) or 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), which serve as general membrane models. Conditions are replicated from experimental studies for comparison purposes, and concentration (X_CER) is varied to probe the effect of CER on these systems. Comparisons with experiment based on deuterium order parameters and bilayer thickness demonstrate good agreement, thus supporting further use of the C36 force field. CER concentration is shown to have a profound effect on nearly all membrane properties including surface area per lipid, chain order and tilt, area compressibility moduli, bilayer thickness, hydrogen bonding, and lipid clustering. Hydrogen bonding in particular can significantly affect other membrane properties and can even encourage transition to a gel phase. Despite CER's tendency to condense the membrane, an expansion of CER lipids with increasing X_CER is possible depending on how the balance between various hydrogen-bond pairs and lipid clustering is perturbed. Based on gel phase transitions, support is given for phytosphingosine's role as a hydrogen-bond bridge between sphingosine ordered domains in the stratum corneum.

Probing the disparate effects of arginine and lysine residues on antimicrobial peptide/bilayer association

Antimicrobial peptides (AMPs) are key components of the innate immune response and represent promising templates for the development of broad-spectrum alternatives to conventional antibiotics. Most AMPs are short, cationic peptides that interact more strongly with negatively charged prokaryotic membranes than net neutral eukaryotic ones. Both AMPs and synthetic analogues with arginine-like side chains are more active against bacteria than those with lysine-like amine groups, though the atomistic mechanism for this increase in potency remains unclear. To examine this, we conducted comparative molecular dynamics simulations of a model negatively-charged membrane system interacting with two mutants of the AMP KR-12: one with lysine residues mutated to arginines (R-KR12) and one with arginine residues mutated to lysine (K-KR12). Simulations show that both partition analogously to the bilayer and display similar preferences for hydrogen bonding with the anionic POPGs. However, R-KR12 binds stronger to the bilayer than K-KR12 and forms significantly more hydrogen bonds, leading to considerably longer interaction times. Additional simulations with methylated R-KR12 and charge-modified K-KR12 mutants show that the extensive interaction seen in the R-KR12 system is partly due to arginine's strong atomic charge distribution, rather than being purely an effect of the greater number of hydrogen bond donors. Finally, free energy simulations reveal that both peptides are disordered in solution but form an amphipathic α-helix when inserted into the bilayer headgroup region. Overall, these results highlight the role of charge and hydrogen bond strength in peptide bilayer insertion, and offer potential insights for designing more potent analogues in the future.

Molecular Understanding of Aβ-hIAPP Cross-Seeding Assemblies on Lipid Membranes

Amyloid-β (Aβ) and human islet polypeptide (hIAPP) are the causative agents responsible for Alzheimer's disease (AD) and type II diabetes (T2D), respectively. While numerous studies have reported the cross-seeding behavior of Aβ and hIAPP in solution, little effort has been made to examine the cross-seeding of Aβ and hIAPP in the presence of cell membranes, which is more biologically relevant to the pathological link between AD and T2D. In this work, we computationally study the cross-seeding and adsorption behaviors of Aβ and hIAPP on zwitterionic POPC and anionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG) mixed bilayers using all-atom molecular dynamics (MD) simulations, particularly aiming to the effects of the initial orientation of the Aβ-hIAPP assembly and the lipid composition of cell membranes on mutual structural and interaction changes in both Aβ-hIAPP assembly and lipid bilayers at the atomic level. Aβ-hIAPP cross-seeding assembly always preferred to adopt a specific orientation and interface to associate with both lipid bilayers strongly via the N-terminal strands of Aβ. Such membrane-bound orientation explains experimental observation that hybrid Aβ-hIAPP fibrils on cell membranes showed similar morphologies to pure hIAPP fibrils. Moreover, Aβ-hIAPP assembly, regardless of its initial orientations, interacted more strongly with POPC/POPG bilayer than POPC bilayer, indicating that electrostatic interactions are the major forces governing peptide-lipid interactions. Strong electrostatic interactions were also attributed to the formation of Ca²⁺ bridges connecting both negatively charged Glu of Aβ and PO₄ head groups of lipids, which facilitate the association of Aβ-hIAPP with the POPC/POPG bilayer. It was also found that the strong peptide-lipid binding reduced lipid fluidity. Both facts imply that Aβ-hIAPP assembly may induce cell damage by altering calcium homeostasis and cell membrane phase. This work provides a better fundamental understanding of cross-seeding of Aβ and hIAPP on cell membranes and a potential pathological link between AD and T2D.

The chain order of binary unsaturated lipid bilayers modulated by aromatic-residue-containing peptides: an ATR-FTIR spectroscopy study

It highlights the importance of aromatic residues in influencing peptide binding to the membrane, demonstrates that the stability of the membranes depends on the lipid composition and the sequence, structural context, and orientation of the peptides.

Kinetic Analysis of the Methyl-β-cyclodextrin-Mediated Intervesicular Transfer of Pyrene-Labeled Phospholipids

Methyl-β-cyclodextrin (MβCD) can transfer phospholipids between vesicles, and its transfer ability has been utilized for the preparation of asymmetric vesicle and lipid incorporation into culture cells. Nevertheless, a detailed kinetic analysis of the MβCD-mediated phospholipid transfer has not yet been carried out. We performed real-time monitoring of intervesicular lipid transfer by means of the fluorescence of pyrene-labeled phospholipids. Intermolecular excimer formation of the pyrene-labeled lipids in a membrane strongly depends on the local concentration of the fluorophore and decreases when the pyrene-labeled lipids are transferred from donor (fluorophore-containing) vesicles to acceptor (fluorophore-free) vesicles. We monitored the fluorescence intensity of the pyrene monomer and excimer simultaneously and found that the excimer/monomer ratio decreased in the presence of MβCD, pointing to MβCD-mediated lipid transfer. The transfer rate depended on the MβCD concentration but not on the lipid concentration, suggesting that dissociation from the membrane via extraction by MβCD is the rate-limiting step of the lipid transfer. Calibration of the excimer/monomer ratio to the molar fraction of the pyrene-labeled lipids enabled us to evaluate the dissociation rate constant correctly. From the temperature dependence of the transfer, we obtained the thermodynamic activation parameters, which revealed that the extraction of phosphatidylcholine by MβCD from membranes is less enthalpically unfavorable than that of phosphatidylglycerol and phosphatidylinositol.