This study employs molecular dynamics (MD) simulations to investigate the adsorption and aggregation behavior of simple polyarginine cell-penetrating peptides (CPPs), specifically modeled as R9 peptides, at zwitterionic phosphocholine POPC membranes under varying ionic strengths of two peptide concentrations and two concentrations of NaCl and CaCl2. The results reveal an intriguing phenomenon of R9 aggregation at the membrane, which is dependent on the ionic strength, indicating a salting-out effect. As the peptide concentration and ionic strength increase, peptide aggregation also increases, with aggregate lifetimes and sizes showing a corresponding rise, accompanied by the total decrease of adsorbed peptides at the membrane surface. Notably, in high ionic strength environments, large R9 aggregates, such as octamers, are also observed occasionally. The salting-out, typically uncommon for short positively charged peptides, is attributed to the unique properties of arginine amino acid, specifically by its side chain containing amphiphilic guanidinium (Gdm+) ion which makes both intermolecular hydrophobic like-charge Gdm+ - Gdm+ and salt-bridge Gdm+ - C-terminus interactions, where the former are increased with the ionic strength, and the latter decreased due to electrostatic screening. The aggregation behavior of R9 peptides at membranes can also be linked to their CPP translocation properties, suggesting that aggregation may aid in translocation across cellular membranes.

• Keywords
• Ionic strength, Molecular dynamics simulations, Peptide aggregation, Phosphocholine lipid bilayers, Polyarginines, Salting-out,
• Publication type
• Journal Article MeSH

Perforation of the outer mitochondrial membrane triggered by BAX and facilitated by its main activator cBID is a fundamental process in cell apoptosis. Here, we employ a newly designed correlative approach based on a combination of a fluorescence cross correlation binding with a calcein permeabilization assay to understand the involvement of BAX in pore formation under oxidative stress conditions. To mimic the oxidative stress, we enriched liposomal membranes by phosphatidylcholines with truncated sn-2 acyl chains terminated by a carboxyl or aldehyde moiety. Our observations revealed that oxidative stress enhances proapoptotic conditions involving accelerated pore-opening kinetics. This enhancement is achieved through increased recruitment of BAX to the membrane and facilitation of BAX membrane insertion. Despite these effects, the fundamental mechanism of pore formation remained unchanged, suggesting an all-or-none mechanism. In line with this mechanism, we demonstrated that the minimal number of BAX molecules at the membrane necessary for pore formation remains constant regardless of BAX activation by cBID or the presence of oxidized lipids. Overall, our findings give a comprehensive picture of the molecular mechanisms underlying apoptotic pore formation and highlight the selective amplifying role of oxidized lipids in triggering formation of membrane pores.

• MeSH
• Apoptosis MeSH
• Fluoresceins chemistry   metabolism   MeSH
• Phosphatidylcholines chemistry   metabolism   MeSH
• Humans MeSH
• Liposomes chemistry   metabolism   MeSH
• Mitochondrial Membranes metabolism   MeSH
• Oxidative Stress * MeSH
• Porosity MeSH
• bcl-2-Associated X Protein * metabolism   chemistry   MeSH
• Check Tag
• Humans MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Fluoresceins MeSH
• fluorexon MeSH Browser
• Phosphatidylcholines MeSH
• Liposomes MeSH
• bcl-2-Associated X Protein * MeSH

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

• MeSH
• Glycosaminoglycans * chemistry   MeSH
• Sulfonic Acids chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Sulfates * chemistry   MeSH
• Static Electricity * MeSH
• Calcium chemistry   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Glycosaminoglycans * MeSH
• Sulfonic Acids MeSH
• Sulfates * MeSH
• sulfamic acid MeSH Browser
• Calcium MeSH

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

• MeSH
• Lipid Bilayers * chemistry   MeSH
• Solvents MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Lipid Bilayers * MeSH
• Solvents MeSH

Ethylenediaminetetraacetic acid (EDTA) is frequently used in lipid experiments to remove redundant ions, such as Ca2+, from the sample solution. In this work, combining molecular dynamics (MD) simulations and Langmuir monolayer experiments, we show that on top of the expected Ca2+ depletion, EDTA anions themselves bind to phosphatidylcholine (PC) monolayers. This binding, originating from EDTA interaction with choline groups of PC lipids, leads to the adsorption of EDTA anions at the monolayer surface and concentration-dependent changes in surface pressure as measured by monolayer experiments and explained by MD simulations. This surprising observation emphasizes that lipid experiments carried out using EDTA-containing solutions, especially of high concentrations, must be interpreted very carefully due to potential interfering interactions of EDTA with lipids and other biomolecules involved in the experiment, e.g., cationic peptides, that may alter membrane-binding affinities of studied compounds.

• MeSH
• Edetic Acid MeSH
• Phosphatidylcholines * chemistry   MeSH
• Ions MeSH
• Membranes, Artificial * MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Edetic Acid MeSH
• Phosphatidylcholines * MeSH
• Ions MeSH
• Membranes, Artificial * MeSH