Ionic Liquid Gating Induces Anomalous Permeation through Membrane Channel Proteins.

Membrane channel proteins (MCPs) play key roles in matter transport through cell membranes and act as major targets for vaccines and drugs. For emerging ionic liquid (IL) drugs, a rational understanding of how ILs affect the structure and transport function of MCP is crucial to their design. In this work, GPU-accelerated microsecond-long molecular dynamics simulations were employed to investigate the modulating mechanism of ILs on MCP. Interestingly, ILs prefer to insert into the lipid bilayer and channel of aquaporin-2 (AQP2) but adsorb on the entrance of voltage-gated sodium channels (Nav). Molecular trajectory and free energy analysis reflect that ILs have a minimal impact on the structure of MCPs but significantly influence MCP functions. It demonstrates that ILs can decrease the overall energy barrier for water through AQP2 by 1.88 kcal/mol, whereas that for Na+ through Nav is increased by 1.70 kcal/mol. Consequently, the permeation rates of water and Na+ can be enhanced and reduced by at least 1 order of magnitude, respectively. Furthermore, an abnormal IL gating mechanism was proposed by combining the hydrophobic nature of MCP and confined water/ion coordination effects. More importantly, we performed experiments to confirm the influence of ILs on AQP2 in human cells and found that treatment with ILs significantly accelerated the changes in cell volume in response to altered external osmotic pressure. Overall, these quantitative results will not only deepen the understanding of IL-cell interactions but may also shed light on the rational design of drugs and disease diagnosis.

Which Coverages of Arc-Shaped Proteins Are Required for Membrane Tubulation?

Arc-shaped BIN/Amphiphysin/Rvs (BAR) domain proteins generate curvature by binding to membranes and induce membrane tubulation at sufficiently large protein coverages. For the amphiphysin N-BAR domain, Le Roux et al., Nat. Commun. 2021, 12, 6550, measured a threshold coverage of 0.44 ± 0.097 for nanotubules emerging from the supported lipid bilayer. In this article, we systematically investigate membrane tubulation induced by arc-shaped protein-like particles with coarse-grained modeling and simulations and determine the threshold coverages at different particle-particle interaction strengths and membrane spontaneous curvatures. In our simulations, the binding of arc-shaped particles induces a membrane shape transition from spherical vesicles to tubules at a particle threshold coverage of about 0.5, which is rather robust to variations of the direct attractive particle interactions or spontaneous membrane curvature in the coarse-grained model. Our study suggests that threshold coverages of around or slightly below 0.5 are a general requirement for membrane tubulation by arc-shaped BAR domain proteins.

Reconstitution and characterization of BRAF in complex with 14-3-3 and KRAS4B on nanodiscs.

RAF kinases are key components of the RAS-MAPK signaling pathway, which drives cell growth and is frequently overactivated in cancer. Upstream signaling activates the small GTPase RAS, which recruits RAF to the cell membrane, driving a transition of the latter from an auto-inhibited monomeric conformation to an active dimer. Despite recent progress, mechanistic details underlying RAF activation remain unclear, particularly the role of RAS and the membrane in mediating this conformational rearrangement of RAF together with 14-3-3 to permit RAF kinase domain dimerization. Here, we reconstituted an active complex of dimeric BRAF, a 14-3-3 dimer and two KRAS4B on a nanodisc bilayer and verified that its assembly is GTP-dependent. Biolayer interferometry (BLI) was used to compare the binding affinities of monomeric versus dimeric full-length BRAF:14-3-3 complexes for KRAS4B-conjugated nanodiscs (RAS-ND) and to investigate the effects of membrane lipid composition and spatial density of KRAS4B on binding. 1,2-Dioleoyl-sn-glycero-3-phospho-L-serine (DOPS) and higher KRAS4B density enhanced the interaction of BRAF:14-3-3 with RAS-ND to different degrees depending on BRAF oligomeric state. We utilized our reconstituted system to dissect the effects of KRAS4B and the membrane on the kinase activity of monomeric and dimeric BRAF:14-3-3 complexes, finding that KRAS4B or nanodiscs alone were insufficient to stimulate activity, whereas RAS-ND increased activity of both states of BRAF. The reconstituted assembly of full-length BRAF with 14-3-3 and KRAS on a cell-free, defined lipid bilayer offers a more holistic biophysical perspective to probe regulation of this multimeric signaling complex at the membrane surface.

In Bacterial Membranes Lipid II Changes the Stability of Pores Formed by the Antimicrobial Peptide Nisin.

Resistance to available antibiotics poses a growing challenge to modern medicine, as this often disallows infections to be controlled. This problem can only be alleviated by the development of new drugs. Nisin, a natural lantibiotic with broad antimicrobial activity, has shown promise as a potential candidate for combating antibiotic-resistant bacteria. However, nisin is poorly soluble and barely stable at physiological pH, which despite attempts to address these issues through mutant design has restricted its use as an antibacterial drug. Therefore, gaining a deeper understanding of the antimicrobial effectiveness, which relies in part on its ability to form pores, is crucial for finding innovative ways to manage infections caused by resistant bacteria. Using large-scale molecular dynamics simulations, we find that the bacterial membrane-specific lipid II increases the stability of pores formed by nisin and that the interplay of nisin and lipid II reduces the overall integrity of bacterial membranes by changing the local thickness and viscosity.

Protein Interface Regulating the Inserting Process of Imidazole Ionic Liquids into the Cell Membrane.

Ionic liquids (ILs) have shown promising potential in membrane protein extraction; however, the underlying mechanism remains unclear. Herein, we employed GPU-accelerated molecular dynamics (MD) simulations to investigate the dynamic insertion process of ILs into cell membranes containing membrane proteins. Our findings reveal that ILs spontaneously insert into the membrane, and the presence of membrane proteins significantly decelerates the rate of IL insertion into the membrane. Specifically, the relationship between the insertion rate and inserting free energy exhibits non-monotonic changes, which can be attributed to interfacial effects. The protein-water interface acts as trap for free ions and ionic clusters, while free ions preferentially insert into the membrane from the protein-lipid interface, which limits the insertion rate due to its narrowness. Thus, the insertion rate is governed by a combination of the free energy and interfacial effects. These findings provide valuable insights into the interfacial effects of protein-lipid bilayers and have implications for various biochemical-related applications.

Multiwaves, breathers, lump and other solutions for the Heimburg model in biomembranes and nerves.

In this manuscript, a mathematical model known as the Heimburg model is investigated analytically to get the soliton solutions. Both biomembranes and nerves can be studied using this model. The cell membrane's lipid bilayer is regarded by the model as a substance that experiences phase transitions. It implies that the membrane responds to electrical disruptions in a nonlinear way. The importance of ionic conductance in nerve impulse propagation is shown by Heimburg's model. The dynamics of the electromechanical pulse in a nerve are analytically investigated using the Hirota Bilinear method. The various types of solitons are investigates, such as homoclinic breather waves, interaction via double exponents, lump waves, multi-wave, mixed type solutions, and periodic cross kink solutions. The electromechanical pulse's ensuing three-dimensional and contour shapes offer crucial insight into how nerves function and may one day be used in medicine and the biological sciences. Our grasp of soliton dynamics is improved by this research, which also opens up new directions for biomedical investigation and medical developments. A few 3D and contour profiles have also been created for new solutions, and interaction behaviors have also been shown.

Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study.

The construction of the stratum corneum (SC) is crucial to the problems of transdermal drug delivery. SC consists of the keratinocyte layers and the lipid matrix surrounding it. Among them, the lipid matrix is the barrier for many exogenous molecules, mainly composed of ceramides (CERs), free fatty acids (FFA), and cholesterol (CHOL). In this work, we developed single-component (CERs, CER-NS, and CER-EOS) and six three-component models, and each model was simulated by using the GROMOS-54A7 force field. Short-period phase (SPP) and long-period phase (LPP) systems were established separately, and area per lipid (APL), thickness, order of carbon chain (SCD), and density distribution were analyzed. The transition of CER-NS and CER-EOS in LPP was observed. The results of hydrogen bonds in the lipid systems indicated that a strong hydrogen-bond network was formed between the skin-lipid bilayers. Umbrella sampling method simulations were performed to calculate the free energy change of ethanol moving into the skin-lipid bilayer. The results revealed that ethanol molecules pulled some water molecules into the membrane when they passed through SPP-1. Our findings provided some insights and models of the stratum corneum that could be used for the subsequent mechanism of macromolecule permeation through membranes in drugs, cosmetics, and so on.

Scaling the Functional Nanopore (FuN) Screen: Systematic Evaluation of Self-Assembling Membrane Peptides and Extension with a K+-Responsive Fluorescent Protein Sensor.

The functional analysis of protein nanopores is typically conducted in planar lipid bilayers or liposomes exploiting high-resolution but low-throughput electrical and optical read-outs. Yet, the reconstitution of protein nanopores in vitro still constitutes an empiric and low-throughput process. Addressing these limitations, nanopores can now be analyzed using the functional nanopore (FuN) screen exploiting genetically encoded fluorescent protein sensors that resolve distinct nanopore-dependent Ca2+ in- and efflux patterns across the inner membrane of Escherichia coli. With a primary proof-of-concept established for the S2168 holin, and thereof based recombinant nanopore assemblies, the question arises to what extent alternative nanopores can be analyzed with the FuN screen and to what extent alternative fluorescent protein sensors can be adapted. Focusing on self-assembling membrane peptides, three sets of 13 different nanopores are assessed for their capacity to form nanopores in the context of the FuN screen. Nanopores tested comprise both natural and computationally designed nanopores. Further, the FuN screen is extended to K+-specific fluorescent protein sensors and now provides a capacity to assess the specificity of a nanopore or ion channel. Finally, a comparison to high-resolution biophysical and electrophysiological studies in planar lipid bilayers provides an experimental benchmark for future studies.

Diffusion and Oligomerization States of the Muscarinic M1 Receptor in Live Cells─The Impact of Ligands and Membrane Disruptors.

G protein-coupled receptors (GPCRs) are a major gateway to cellular signaling, which respond to ligands binding at extracellular sites through allosteric conformational changes that modulate their interactions with G proteins and arrestins at intracellular sites. High-resolution structures in different ligand states, together with spectroscopic studies and molecular dynamics simulations, have revealed a rich conformational landscape of GPCRs. However, their supramolecular structure and spatiotemporal distribution is also thought to play a significant role in receptor activation and signaling bias within the native cell membrane environment. Here, we applied single-molecule fluorescence techniques, including single-particle tracking, single-molecule photobleaching, and fluorescence correlation spectroscopy, to characterize the diffusion and oligomerization behavior of the muscarinic M1 receptor (M1R) in live cells. Control samples included the monomeric protein CD86 and fixed cells, and experiments performed in the presence of different orthosteric M1R ligands and of several compounds known to change the fluidity and organization of the lipid bilayer. M1 receptors exhibit Brownian diffusion characterized by three diffusion constants: confined/immobile (∼0.01 µm2/s), slow (∼0.04 µm2/s), and fast (∼0.14 µm2/s), whose populations were found to be modulated by both orthosteric ligands and membrane disruptors. The lipid raft disruptor C6 ceramide led to significant changes for CD86, while the diffusion of M1R remained unchanged, indicating that M1 receptors do not partition in lipid rafts. The extent of receptor oligomerization was found to be promoted by increasing the level of expression and the binding of orthosteric ligands; in particular, the agonist carbachol elicited a large increase in the fraction of M1R oligomers. This study provides new insights into the balance between conformational and environmental factors that define the movement and oligomerization states of GPCRs in live cells under close-to-native conditions.

Effect of GM1 concentration change on plasma membrane: molecular dynamics simulation and analysis.

Ganglioside GM1 is a class of glycolipids predominantly located in the nervous system. Comprising a ceramide anchor and an oligosaccharide chain containing sialic acid, GM1 plays a pivotal role in various cellular processes, including signal transduction, cell adhesion, and membrane organization. Moreover, GM1 has been implicated in the pathogenesis of several neurological disorders, such as Parkinson's disease, Alzheimer's disease, and stroke. In this study, by creating a neural cell model membrane simulation system and employing rigorous molecular models, we utilize a coarse-grained molecular dynamics approach to explore the structural and dynamic characteristics of multi-component neuronal plasma membranes at varying GM1 ganglioside concentrations. The simulation results reveal that as GM1 concentration increases, a greater number of hydrogen bonds form between GM1 molecules, resulting in the formation of larger clusters, which leads to reduced membrane fluidity, increased lipid ordering, decreased membrane thickness and surface area and higher levels of GM1 dissociation. Through a meticulous analysis, while considering GM1's structural attributes, we offer valuable insights into the structural and dynamic traits of the cell membrane. This study provides a robust methodology for exploring membrane characteristics and enhances our comprehension of GM1 molecules, serving as a resource for both experimental and computational researchers in this field.

Probing the stability and interdomain interactions in the ABC transporter OpuA using single-molecule optical tweezers.

Transmembrane transporter proteins are essential for maintaining cellular homeostasis and, as such, are key drug targets. Many transmembrane transporter proteins are known to undergo large structural rearrangements during their functional cycles. Despite the wealth of detailed structural and functional data available for these systems, our understanding of their dynamics and, consequently, how they function is generally limited. We introduce an innovative approach that enables us to directly measure the dynamics and stability of interdomain interactions of transmembrane proteins using optical tweezers. Focusing on the osmoregulatory ATP-binding cassette transporter OpuA from Lactococcus lactis, we examine the mechanical properties and potential interactions of its substrate-binding domains. Our measurements are performed in lipid nanodiscs, providing a native-mimicking environment for the transmembrane protein. The technique provides high spatial and temporal resolution and allows us to study the functionally relevant motions and interdomain interactions of individual transmembrane transporter proteins in real time in a lipid bilayer.

Water Hydrogen-Bond Mediated Layer by Layer Alignment of Lipid Rafts as a Precursor of Intermembrane Processes.

Lipid rafts, which are dynamic nanodomains in the plasma membrane, play a crucial role in intermembrane processes by clustering together and growing in size within the plane of the membrane while also aligning with each other across different membranes. However, the physical origin of layer by layer alignment of lipid rafts remains to be elucidated. Here, by using fluorescence imaging and synchrotron X-ray reflectivity in a phase-separated multilayer system, we find that the alignment of raft-mimicking Lo domains is regulated by the distance between bilayers. Molecular dynamics simulations reveal that the aligned state is energetically preferred when the intermembrane distance is small due to its ability to minimize the volume of surface water, which has fewer water hydrogen bonds (HBs) compared to bulk water. Our results suggest that water HB-driven alignment of lipid rafts plays a role as a precursor of intermembrane processes such as cell-cell fusion, virus entry, and signaling.

Thermodynamics of oligomerization and Helix-to-sheet structural transition of amyloid ß-protein on anionic phospholipid vesicles.

Understanding oligomerization and aggregation of the amyloid-ß protein is important to elucidate the pathological mechanisms of Alzheimer's disease, and lipid membranes play critical roles in this process. In addition to studies reported by other groups, our group has also reported that the negatively-charged lipid bilayers with a high positive curvature induced α-helix-to-ß-sheet conformational transitions of amyloid-ß-(1-40) upon increase in protein density on the membrane surface and promoted amyloid fibril formation of the protein. Herein, we investigated detailed mechanisms of the conformational transition and oligomer formation of the amyloid-ß protein on the membrane surface. Changes in the fractions of the three protein conformers (free monomer, membrane-bound α-helix-rich conformation, and ß-sheet-rich conformation) were determined from the fluorescent spectral changes of the tryptophan probe in the protein. The helix-to-sheet structural transition on the surface was described by a thermodynamic model of octamer formation driven by entropic forces including hydrophobic interactions. These findings provide useful information for understanding the self-assembly of amyloidogenic proteins on lipid membrane surfaces.

CHARMM GUI Membrane Builder for oxidized phospholipid membrane modeling and simulation.

Oxidative stress leads to the production of oxidized phospholipids (oxPLs) that modulate the biophysical properties of phospholipid monolayers and bilayers. As many immune cells are responsible for surveilling cells and tissues for the presence of oxPLs, oxPL-dependent mechanisms have been suggested as targets for treating chronic kidney disease, atherosclerosis, diabetes, and cancer metastasis. This review details recent experimental and computational studies that characterize oxPLs' behaviors in various monolayers and bilayers. These studies investigate how the tail length and polar functional groups of OxPLs impact membrane properties, how oxidized membranes can be stabilized, and how membrane integrity is generally affected by oxidized lipids. In addition, for oxPL-containing membrane modeling and simulation, CHARMM-GUI Membrane Builder has been extended to support a variety of oxPLs, accelerating the simulation system building process for these biologically relevant lipid bilayers.

NMR methods to detect fluoride binding and transport by membrane proteins.

Solid-state nuclear magnetic resonance (NMR) methods can probe the motions of membrane proteins in liposomes at the atomic level, and propel the understanding of biomolecular processes for which static structures cannot provide a satisfactory description. High-resolution crystallography snapshots have provided a structural basis for fluoride channels. NMR is a powerful tool to build upon these snapshots and depict a dynamic picture of fluoride channels in native-like lipid bilayers. In this contribution, we discuss solid-state and solution NMR experiments to detect fluoride binding and transport by fluoride channels. Ongoing developments in membrane protein sample preparation and ssNMR methodology, particularly in using 1H, 19F and 13C-detection schemes, offer additional opportunities to study structure and functional aspects of fluoride channels.

High Charge Density in Peptide Dendrimers is Required to Destabilize Membranes: Insights into Endosome Evasion.

Peptide dendrimers are a type of branched, symmetric, and topologically well-defined molecule that have already been used as delivery systems for nucleic acid transfection. Several of the most promising sequences showed high efficiency in many key steps of transfection, namely, binding siRNA, entering cells, and evading the endosome. However, small changes to the peptide dendrimers, such as in the hydrophobic core, the amino acid chirality, or the total available charges, led to significantly different experimental results with unclear mechanistic insights. In this work, we built a computational model of several of those peptide dendrimers (MH18, MH13, and MH47) and some of their variants to study the molecular details of the structure and function of these molecules. We performed CpHMD simulations in the aqueous phase and in interaction with a lipid bilayer to assess how conformation and protonation are affected by pH in different environments. We found that while the different peptide dendrimer sequences lead to no substantial structural differences in the aqueous phase, the total charge and, more importantly, the total charge density are key for the capacity of the dendrimer to interact and destabilize the membrane. These dendrimers become highly charged when the pH changes from 7.5 to 4.5, and the presence of a high charge density, which is decreased for MH47 that has four fewer titratable lysines, is essential to trigger membrane destabilization. These findings are in excellent agreement with the experimental data and help us to understand the high efficiency of some dendrimers and why the dendrimer MH47 is unable to complete the transfection process. This evidence provides further understanding of the mode of action of these peptide dendrimers and will be pivotal for the future design of new sequences with improved transfection capabilities.

Membrane damage by MBP-1 is mediated by pore formation and amplified by mtDNA.

Eosinophils play a crucial role in host defense while also contributing to immunopathology through the release of inflammatory mediators. Characterized by distinctive cytoplasmic granules, eosinophils securely store and rapidly release various proteins exhibiting high toxicity upon extracellular release. Among these, major basic protein 1 (MBP-1) emerges as an important mediator in eosinophil function against pathogens and in eosinophil-associated diseases. While MBP-1 targets both microorganisms and host cells, its precise mechanism remains elusive. We demonstrate that formation of small pores by MBP-1 in lipid bilayers induces membrane permeabilization and disrupts potassium balance. Additionally, we reveal that mitochondrial DNA (mtDNA) present in eosinophil extracellular traps (EETs) amplifies MBP-1 toxic effects, underscoring the pivotal role of mtDNA in EETs. Furthermore, we present evidence indicating that absence of CpG methylation in mtDNA contributes to the regulation of MBP-1-mediated toxicity. Taken together, our data suggest that the mtDNA scaffold within extracellular traps promotes MBP-1 toxicity.

Vesicle protrusion induced by antimicrobial peptides suggests common carpet mechanism for short antimicrobial peptides.

Short-cationic alpha-helical antimicrobial peptides (SCHAMPs) are promising candidates to combat the growing global threat of antimicrobial resistance. They are short-sequenced, selective against bacteria, and have rapid action by destroying membranes. A full understanding of their mechanism of action will provide key information to design more potent and selective SCHAMPs. Molecular Dynamics (MD) simulations are invaluable tools that provide detailed insights into the peptide-membrane interaction at the atomic- and meso-scale level. We use atomistic and coarse-grained MD to look into the exact steps that four promising SCHAMPs-BP100, Decoralin, Neurokinin-1, and Temporin L-take when they interact with membranes. Following experimental set-ups, we explored the effects of SCHAMPs on anionic membranes and vesicles at multiple peptide concentrations. Our results showed all four peptides shared similar binding steps, initially binding to the membrane through electrostatic interactions and then flipping on their axes, dehydrating, and inserting their hydrophobic moieties into the membrane core. At higher concentrations, fully alpha-helical peptides induced membrane budding and protrusions. Our results suggest the carpet mode of action is fit for the description of SCHAMPs lysis activity and discuss the importance of large hydrophobic residues in SCHAMPs design and activity.

Mechanisms of the Drug Penetration Enhancer Propylene Glycol Interacting with Skin Lipid Membranes.

Very few drugs have the necessary physicochemical properties to cross the skin's main permeability barrier, the stratum corneum (SC), in sufficient amounts. Propylene glycol (PG) is a chemical penetration enhancer that could be included in topical formulations in order to overcome the barrier properties of the skin and facilitate the transport of drugs across it. Experiments have demonstrated that PG increases the mobility and disorder of SC lipids and may extract cholesterol from the SC, but little is known about the molecular mechanisms of drug permeation enhancement by PG. In this work, we have performed molecular dynamics (MD) simulations to investigate the molecular-level effects of PG on the structure and properties of model SC lipid bilayers. The model bilayers were simulated in the presence of PG concentrations over the range of 0-100% w/w PG, using both an all-atom and a united atom force field. PG was found to localize in the hydrophilic headgroup regions at the bilayer interface, to occupy the lipid-water hydrogen-bonding sites, and to slightly increase lipid tail disorder in a concentration-dependent manner. We showed with MD simulation that PG enhances the permeation of small molecules such as water by interacting with the bilayer interface; the results of our study may be used to guide the design of formulations for transdermal drug delivery with enhanced skin permeation, as well as topical formulations and cosmetic products.

Computational Insights into the Conformational Dynamics of HIV-1 Vpr in a Lipid Bilayer for Ion Channel Modeling.

HIV-1 Vpr is a multifunctional accessory protein consisting of 96 amino acids that play a critical role in viral pathogenesis. Among its diverse range of activities, Vpr can create a cation-selective ion channel within the plasma membrane. However, the oligomeric state of this channel has not yet been elucidated. In this study, we investigated the conformational dynamics of Vpr helices to model the ion channel topology. First, we employed a series of multiscale simulations to investigate the specific structure of monomeric Vpr in a membrane model. During the lipid bilayer self-assembly coarse grain simulation, the C-terminal helix (residues 56-77) effectively formed the transmembrane region, while the N-terminal helix exhibited an amphipathic nature by associating horizontally with a single leaflet. All-atom molecular dynamics (MD) simulations of full-length Vpr inside a phospholipid bilayer show that the C-terminal helix remains very stable inside the bilayer core in a vertical orientation. Subsequently, using the predicted C-terminal helix orientation and conformation, various oligomeric states (ranging from tetramer to heptamer) possibly forming the Vpr ion channel were built and further evaluated. Among these models, the pentameric form exhibited consistent stability in MD simulations and displayed a compatible conformation for a water-assisted ion transport mechanism. This study provides structural insights into the ion channel activity of the Vpr protein and the foundation for developing therapeutics against HIV-1 Vpr-related conditions.