Influence of Monovalent Cation Size on Nanodomain Formation in Anionic–Zwitterionic Mixed Bilayers

Molecular Dynamics Study of the Effect of Charge and Glycosyl on Superoxide Anion Distribution near Lipid Membrane

To examine the effects of membrane charge, the electrolyte species and glycosyl on the distribution of negatively charged radical of superoxide anion (·O₂^-) around the cell membrane, different phospholipid bilayer systems containing ·O₂^- radicals, different electrolytes and phospholipid bilayers were constructed through Charmm-GUI and Amber16. These systems were equilibrated with molecular dynamics by using Gromacs 5.0.2 to analyze the statistical behaviors of ·O₂^- near the lipid membrane under different conditions. It was found that in the presence of potassium rather than sodium, the negative charge of the phospholipid membrane is more likely to rarefy the superoxide anion distribution near the membrane surface. Further, the presence of glycosyl significantly reduced the density of ·O₂^- near the phospholipid bilayer by 78.3% compared with that of the neutral lipid membrane, which may have a significant contribution to reducing the lipid peroxidation from decreasing the ·O₂^- density near the membrane.

Transferable and Polarizable Coarse Grained Model for Proteins─ProMPT

The application of classical molecular dynamics (MD) simulations at atomic resolution (fine-grained level, FG), to most biomolecular processes, remains limited because of the associated computational complexity of representing all the atoms. This problem is magnified in the presence of protein-based biomolecular systems that have a very large conformational space, and MD simulations with fine-grained resolution have slow dynamics to explore this space. Current transferable coarse grained (CG) force fields in literature are either limited to only peptides with the environment encoded in an implicit form or cannot capture transitions into secondary/tertiary peptide structures from a primary sequence of amino acids. In this work, we present a transferable CG force field with an explicit representation of the environment for accurate simulations with proteins. The force field consists of a set of pseudoatoms representing different chemical groups that can be joined/associated together to create different biomolecular systems. This preserves the transferability of the force field to multiple environments and simulation conditions. We have added electronic polarization that can respond to environmental heterogeneity/fluctuations and couple it to protein's structural transitions. The nonbonded interactions are parametrized with physics-based features such as solvation and partitioning free energies determined by thermodynamic calculations and matched with experiments and/or atomistic simulations. The bonded potentials are inferred from corresponding distributions in nonredundant protein structure databases. We present validations of the CG model with simulations of well-studied aqueous protein systems with specific protein fold types─Trp-cage, Trpzip4, villin, WW-domain, and β-α-β. We also explore the applications of the force field to study aqueous aggregation of Aβ 16-22 peptides.

Effects of applied surface-tension on membrane-assisted Aβ aggregation

Accumulation of protein-based (Aβ) aggregates on cellular membranes with varying structural properties is commonly recognized as the key step in Alzheimer's pathogenesis. But experimental and computational challenges have made this biophysical characterization difficult. In particular, studies connecting biological membrane organization and Aβ aggregation are limited. While experiments have suggested that an increased membrane curvature results in faster Aβ peptide aggregation in the context of Alzheimer's disease, a mechanistic explanation for this relation is missing. In this work, we are leveraging molecular simulations with a physics-based coarse grained model to address and understand the relationships between curved cellular membranes and aggregation of a model template peptide Aβ 16-22. In agreement with experimental results, our simulations also suggest a positive correlation between increased peptide aggregation and membrane curvature. More curved membranes have higher lipid packing defects that engage peptide hydrophobic groups and promote faster diffusion leading to peptide fibrillar structures. In addition, we curated the effects of peptide aggregation on the membrane's structure and organization. Interfacial peptide aggregation results in heterogeneous headgroup-peptide interactions and an induced crowding effect at the lipid headgroup region, leading to a more ordered headgroup region and disordered lipid-tails at the membrane core. This work presents a mechanistic and morphological overview of the relationships between the biomembrane local structure and organization, and Aβ peptide aggregation.

Water Pores in Planar Lipid Bilayers at Fast and Slow Rise of Transmembrane Voltage

Basic understanding of the barrier properties of biological membranes can be obtained by studying model systems, such as planar lipid bilayers. Here, we study water pores in planar lipid bilayers in the presence of transmembrane voltage. Planar lipid bilayers were exposed to fast and slow linearly increasing voltage and current signals. We measured the capacitance, breakdown voltage, and rupture time of planar lipid bilayers composed of 1-pamitoyl 2-oleoyl phosphatidylcholine (POPC), 1-pamitoyl 2-oleoyl phosphatidylserine (POPS), and a mixture of both lipids in a 1:1 ratio. Based on the measurements, we evaluated the change in the capacitance of the planar lipid bilayer corresponding to water pores, the radius of water pores at membrane rupture, and the fraction of the area of the planar lipid bilayer occupied by water pores.planar lipid bilayer capacitance, which corresponds to water pores, water pore radius at the membrane rupture, and a fraction of the planar lipid bilayer area occupied by water pores. The estimated pore radii determining the rupture of the planar lipid bilayer upon fast build-up of transmembrane voltage are 0.101 nm, 0.110 nm, and 0.106 nm for membranes composed of POPC, POPS, and POPC:POPS, respectively. The fraction of the surface occupied by water pores at the moment of rupture of the planar lipid bilayer The fraction of an area that is occupied by water pores at the moment of planar lipid bilayer rupture is in the range of 0.1–1.8%.

Structural information and membrane binding of truncated RGS9-1 Anchor Protein and its C-terminal hydrophobic segment

Visual phototransduction takes place in photoreceptor cells. Light absorption by rhodopsin leads to the activation of transducin as a result of the exchange of its GDP for GTP. The GTP-bound ⍺-subunit of transducin then activates phosphodiesterase (PDE), which in turn hydrolyzes cGMP leading to photoreceptor hyperpolarization. Photoreceptors return to the dark state upon inactivation of these proteins. In particular, PDE is inactivated by the protein complex R9AP/RGS9-1/Gβ5. R9AP (RGS9-1 anchor protein) is responsible for the membrane anchoring of this protein complex to photoreceptor outer segment disk membranes most likely by the combined involvement of its C-terminal hydrophobic domain as well as other types of interactions. This study thus aimed to gather information on the structure and membrane binding of the C-terminal hydrophobic segment of R9AP as well as of truncated R9AP (without its C-terminal domain, R9AP∆TM). Circular dichroism and infrared spectroscopic measurements revealed that the secondary structure of R9AP∆TM mainly includes ⍺-helical structural elements. Moreover, intrinsic fluorescence measurements of native R9AP∆TM and individual mutants lacking one tryptophan demonstrated that W79 is more buried than W173 but that they are both located in a hydrophobic environment. This method also revealed that membrane binding of R9AP∆TM does not involve regions near its tryptophan residues, while infrared spectroscopy validated its binding to lipid vesicles. Additional fluorescence measurements showed that the C-terminal segment of R9AP is membrane embedded. Maximum insertion pressure and synergy data using Langmuir monolayers suggest that interactions with specific phospholipids could be involved in the membrane binding of R9AP∆TM.

Microscopic Picture of Calcium-Assisted Lipid Demixing and Membrane Remodeling Using Multiscale Simulations

The specificity of anionic phospholipids-calcium ion interaction and lipid demixing has been established as a key regulatory mechanism in several cellular signaling processes. The mechanism and implications of this calcium-assisted demixing have not been elucidated from a microscopic point of view. Here, we present an overview of atomic interactions between calcium and phospholipids that can drive nonideal mixing of lipid molecules in a model lipid bilayer composed of zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) lipids with computer simulations at multiple resolutions. Lipid nanodomain formation and growth were driven by calcium-enabled lipid bridging of the charged phosphatidylserine (PS) headgroups, which were favored against inter-POPS dipole interactions. Consistent with several experimental studies of calcium-associated membrane sculpting, our analyses also suggest modifications in local membrane curvature and cross-leaflet couplings as a response to such induced lateral heterogeneity. In addition, reverse mapping to a complementary atomistic description revealed structural insights in the presence of anionic nanodomains, at timescales not accessed by previous computational studies. This work bridges information across multiple scales to reveal a mechanistic picture of calcium ion's impact on membrane biophysics.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Computational insights into lipid assisted peptide misfolding and aggregation in neurodegeneration

An overview of recent advances in computational investigation of peptide–lipid interactions in neurodegeneration – Alzheimer's, Parkinson's and Huntington's disease.

Pathways of amyloid-beta absorption and aggregation in a membranous environment

Aggregation of misfolded oligomeric amyloid-beta (Aβ) peptides on lipid membranes has been identified as a primary event in Alzheimer's pathogenesis. However, the structural and dynamical features of this membrane assisted Aβ aggregation have not been well characterized. The microscopic characterization of dynamic molecular-level interactions in peptide aggregation pathways has been challenging both computationally and experimentally. In this work, we explore differential patterns of membrane-induced Aβ 16-22 (K-L-V-F-F-A-E) aggregation from the microscopic perspective of molecular interactions. Physics-based coarse-grained molecular dynamics (CG-MD) simulations were employed to investigate the effect of lipid headgroup charge - zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPC) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine: POPS) - on Aβ 16-22 peptide aggregation. Our analyses present an extensive overview of multiple pathways for peptide absorption and biomechanical forces governing peptide folding and aggregation. In agreement with experimental observations, anionic POPS molecules promote extended configurations in Aβ peptides that contribute towards faster emergence of ordered β-sheet-rich peptide assemblies compared to POPC, suggesting faster fibrillation. In addition, lower cumulative rates of peptide aggregation in POPS due to higher peptide-lipid interactions and slower lipid diffusion result in multiple distinct ordered peptide aggregates that can serve as nucleation seeds for subsequent Aβ aggregation. This study provides an in-silico assessment of experimentally observed aggregation patterns, presents new morphological insights and highlights the importance of lipid headgroup chemistry in modulating the peptide absorption and aggregation process.

Effect of Curcumin Addition on the Adsorption and Transport of a Cationic Dye across DPPG-POPG Liposomes Probed by Second Harmonic Spectroscopy

The effect of addition of curcumin on the adsorption and transport characteristics of a cationic dye, LDS⁺, across negatively charged bilayers composed of POPG and DPPG lipids were investigated by the interface selective second harmonic (SH) spectroscopic technique. Curcumin induced changes in the SH electric field signal of the LDS⁺ ions (E_2ω (LDS⁺)) were observed to depend critically on the bilayer acyl chain saturation/unsaturation ratio (S/U). Following earlier works, the increase in the E_2ω (LDS⁺) signal is attributed to the release of the Na⁺ counterions present in the head group region of the bilayer by curcumin and the decay of the E_2ω (LDS⁺) signal is attributed to the bilayer intercalated state of curcumin. While the changes observed in the E_2ω (LDS⁺) signal in the presence of POPG liposomes were consistent with our earlier study ( Varshney, G. K. et al. Langmuir , 2016 , 32 , 10415 - 10421 ), they were significantly different for DPPG liposomes, following curcumin addition. While the increase in the E_2ω (LDS⁺) signal in the presence of POPG liposomes, is marginal (∼10-20%) and instantaneous (<1 s) followed by a rapid decay (completed within ∼100 s), in the presence of DPPG liposomes it was observed to increase slowly and at saturation shows a substantial increase (100-200%), following curcumin addition. When liposomes consisting of a mixture of POPG and DPPG lipids are used, curcumin induced kinetic characteristics of the E_2ω (LDS⁺) signal showed a mixture of the individual kinetic characteristics observed for the unsaturated (POPG) and saturated (DPPG) liposomes. The observed kinetic trends of the E_2ω (LDS⁺) signal following curcumin addition are explained on the basis of the relative strength of the Na⁺-POPG and Na⁺-DPPG interaction. Higher ordering of the lipid acyl chain region in DPPG liposome makes the Na⁺-DPPG interaction much stronger than the Na⁺-POPG interaction. Further, it is proposed that, in POPG-DPPG liposomes, individual domains of POPG and DPPG lipids exist at low temperature as suggested by the observed temperature dependent kinetic characteristics of the E_2ω (LDS⁺) signal following curcumin addition. These domains are dependent on the S/U ratio and phase state of the bilayer. The gel phase was observed to be more conducive for individual domain formation. Results presented in this work not only support the notion that biological activity of curcumin is associated with its bilayer altering properties, but more interestingly it provides a qualitative insight about how bilayer phase separation can be achieved by modulating the hydrophobic interactions between the lipid acyl chains.