Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature

Absorption of Polypropylene in Dipalmitoylphosphatidylcholine Membranes: The Role of Molecular Weight and Initial Configuration of Polymer Chains

We study by molecular dynamics simulations the absorption of polypropylene (PP) chains within a dipalmitoylphosphatidylcholine (DPPC) lipid membrane in aqueous solvent. DPPC represents the most abundant phospholipid in biological membranes, while PP is one of the most common synthetic polymers diffused in the anthropic environment. By following in detail the absorption process, and the corresponding structural modification undergone by the membrane, we show how the initial configuration and the PP molecular weight determine the overall behavior of the system. Specifically, if PP chains initially lie on the DPPC surface, they are fully absorbed; likewise, polymers initially included within the membrane cannot escape from. On the other hand, if polymers are placed sufficiently apart from the membrane, they have time to join together and coalesce into a few nanoparticles. At contact, such nanoparticles may completely dissolve (for low molecular weight) and then be absorbed. For high molecular weight, not all of them dissolve, and therefore the system attains a condition in which some of the chains are absorbed, while others form a residual nanoparticle staying outside (but in contact with) the membrane. Such a state─albeit energetically unfavorable with respect to a condition in which all PP chains are absorbed─remains stable, at the least over a substantial simulation time, extending in our study up to 1.6 μs. The tendency for polymers to spontaneously form aggregates, which then prefer to stay in contact with the membrane, is further corroborated by calculation of the DPPC-nanoparticle potential of mean force.

Changes in free energy barrier for water permeation by stretch-induced phase transitions in phospholipid/cholesterol bilayers

Water permeation through phospholipid/cholesterol bilayers is the key to understanding tension-induced rupture of biological cell membranes. We performed molecular dynamics simulations of stretched phospholipid/cholesterol bilayers to investigate changes in the free energy profile of water molecules across the bilayer and the lipid structure responsible for water permeation. We modeled stretching of the bilayer by applying areal strain. In stretched phospholipid/cholesterol bilayers, the hydrophobic tail of the phospholipids became disordered and the free energy barrier to water permeation decreased. Upon exceeding the critical areal strain, a phase transition to an interdigitated gel phase occurred before rupture, and the hydrophobic tail ordering as well as the free energy barrier were restored. In pure phospholipid bilayers, we did not observe such recoveries. These transient recoveries in the phospholipid/cholesterol bilayer suppressed water permeation and membrane rupture, followed by an increase in the critical areal strain at which the bilayer ruptured. This result agrees with experimental results and provides a reasonable molecular mechanism for the toughness of phospholipid/cholesterol bilayers under tension.Communicated by Ramaswamy H. Sarma.

Pathways to a Shiny Future: Building the Foundation for Computational Physical Chemistry and Biophysics in 2050

In the last quarter-century, the field of molecular dynamics (MD) has undergone a remarkable transformation, propelled by substantial enhancements in software, hardware, and underlying methodologies. In this Perspective, we contemplate the future trajectory of MD simulations and their possible look at the year 2050. We spotlight the pivotal role of artificial intelligence (AI) in shaping the future of MD and the broader field of computational physical chemistry. We outline critical strategies and initiatives that are essential for the seamless integration of such technologies. Our discussion delves into topics like multiscale modeling, adept management of ever-increasing data deluge, the establishment of centralized simulation databases, and the autonomous refinement, cross-validation, and self-expansion of these repositories. The successful implementation of these advancements requires scientific transparency, a cautiously optimistic approach to interpreting AI-driven simulations and their analysis, and a mindset that prioritizes knowledge-motivated research alongside AI-enhanced big data exploration. While history reminds us that the trajectory of technological progress can be unpredictable, this Perspective offers guidance on preparedness and proactive measures, aiming to steer future advancements in the most beneficial and successful direction.

Hydration- and Temperature-Dependent Fluorescence Spectra of Laurdan Conformers in a DPPC Membrane

The widely used Laurdan probe has two conformers, resulting in different optical properties when embedded in a lipid bilayer membrane, as demonstrated by our previous simulations. Up to now, the two conformers' optical responses have, however, not been investigated when the temperature and the phase of the membrane change. Since Laurdan is known to be both a molecular rotor and a solvatochromic probe, it is subject to a profound interaction with both neighboring lipids and water molecules. In the current study, molecular dynamics simulations and hybrid Quantum Mechanics/Molecular Mechanics calculations are performed for a DPPC membrane at eight temperatures between 270K and 320K, while the position, orientation, fluorescence lifetime and fluorescence anisotropy of the embedded probes are monitored. The importance of both conformers is proven through a stringent comparison with experiments, which corroborates the theoretical findings. It is seen that for Conf-I, the excited state lifetime is longer than the relaxation of the environment, while for Conf-II, the surroundings are not yet adapted when the probe returns to the ground state. Throughout the temperature range, the lifetime and anisotropy decay curves can be used to identify the different membrane phases. The current work might, therefore, be of importance for biomedical studies on diseases, which are associated with cell membrane transformations.

Lateral diffusion of ions near membrane surface

Biological membranes isolate living cells from their environment, while allowing selective molecular transport between the inner and outer realms. For example, Na+ and K+ permeability through ionic channels contributes to neural conduction. Whether the ionic currents arise directly from cations in the bulk, or from the interface, is currently unclear. There are only scant results concerning lateral diffusion of ions on aquated membrane surfaces (and strong belief that this occurs through binding to a diffusing lipid). We performed classical molecular dynamics (MD) simulations of monovalent ions, Na+, K+, and Cl-, near the surface of the zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) membrane. Realistic force-fields for lipids (Amber's Lipid17 and Lipid21) and water (TIP4P-Ew) are tested for the mass and charge densities and the electrostatic potential across the membrane. These calculations reveal that the chloride can bind to the choline moiety through an intervening water molecule by forming a CH⋯OH hydrogen bond, while cations bind to both the phosphatic and carbonyl oxygens of phosphatidylcholine moieties. Upon transitioning from the bulk to the interface, a cation sheds some of its hydration water, which are replaced by headgroup atoms. Notably, an interfacial cation can bind 1-4 headgroup atoms, which is a key to understanding its surface hopping mechanism. We find that cation binding to three headgroup atoms immobilizes it, while binding to four energizes it. Consequently, the lateral cation diffusion rate is only 15-25 times slower than in the bulk, and 4-5 times faster than lipid self-diffusion. K+ diffusion is notably more anomalous than Na+, switching from sub- to super-diffusion after about 2 ns.

Learning Protein-Ligand Unbinding Pathways via Single-Parameter Community Detection

Understanding the dynamics of biomolecular complexes, e.g., of protein-ligand (un)binding, requires the comprehension of paths such systems take between metastable states. In MD simulations, paths are usually not observable per se, but they need to be inferred from simulation trajectories. Here, we present a novel approach to cluster trajectories based on a community detection algorithm that necessitates only the definition of a single parameter. The unbinding of the streptavidin-biotin complex is used as a benchmark system and the A2a adenosine receptor in complex with the inhibitor ZM241385 as an elaborate application. We demonstrate how such clusters of trajectories correspond to pathways and how the approach helps in the identification of reaction coordinates for a considered (un)binding process.

Understanding Aβ Peptide Binding to Lipid Membranes: A Biophysical Perspective

Aβ peptides are known to bind neural plasma membranes in a process leading to the deposit of Aβ-enriched plaques. These extracellular structures are characteristic of Alzheimer's disease, the major cause of late-age dementia. The mechanisms of Aβ plaque formation and deposition are far from being understood. A vast number of studies in the literature describe the efforts to analyze those mechanisms using a variety of tools. The present review focuses on biophysical studies mostly carried out with model membranes or with computational tools. This review starts by describing basic physical aspects of lipid phases and commonly used model membranes (monolayers and bilayers). This is followed by a discussion of the biophysical techniques applied to these systems, mainly but not exclusively Langmuir monolayers, isothermal calorimetry, density-gradient ultracentrifugation, and molecular dynamics. The Methodological Section is followed by the core of the review, which includes a summary of important results obtained with each technique. The last section is devoted to an overall reflection and an effort to understand Aβ-bilayer binding. Concepts such as Aβ peptide membrane binding, adsorption, and insertion are defined and differentiated. The roles of membrane lipid order, nanodomain formation, and electrostatic forces in Aβ-membrane interaction are separately identified and discussed.

Quantum dipole interactions and transient hydrogen bond orientation order in cells, cellular membranes and myelin sheath: Implications for MRI signal relaxation, anisotropy, and T1 magnetic field dependence

PURPOSE

Despite significant impact on the study of human brain, MRI lacks a theory of signal formation that integrates quantum interactions involving proton dipoles (a primary MRI signal source) with brain intricate cellular environment. The purpose of the present study is developing such a theory.

METHODS

We introduce the Transient Hydrogen Bond (THB) model, where THB-mediated quantum dipole interactions between water and protons of hydrophilic heads of amphipathic biomolecules forming cells, cellular membranes and myelin sheath serve as a major source of MR signal relaxation.

RESULTS

The THB theory predicts the existence of a hydrogen-bond-driven structural order of dipole-dipole connections within THBs as a primary factor for the anisotropy observed in MRI signal relaxation. We have also demonstrated that the conventional Lorentzian spectral density function decreases too fast at high frequencies to adequately capture the field dependence of brain MRI signal relaxation. To bridge this gap, we introduced a stretched spectral density function that surpasses the limitations of Lorentzian dispersion. In human brain, our findings reveal that at any time point only about 4% to 7% of water protons are engaged in quantum encounters within THBs. These ultra-short (2 to 3 ns), but frequent quantum spin exchanges lead to gradual recovery of magnetization toward thermodynamic equilibrium, that is, relaxation of MRI signal.

CONCLUSION

By incorporating quantum proton interactions involved in brain imaging, the THB approach introduces new insights on the complex relationship between brain tissue cellular structure and MRI measurements, thus offering a promising new tool for better understanding of brain microstructure in health and disease.

Atomistic mechanism of coupling between cytosolic sensor domain and selectivity filter in TREK K2P channels

The two-pore domain potassium (K2P) channels TREK-1 and TREK-2 link neuronal excitability to a variety of stimuli including mechanical force, lipids, temperature and phosphorylation. This regulation involves the C-terminus as a polymodal stimulus sensor and the selectivity filter (SF) as channel gate. Using crystallographic up- and down-state structures of TREK-2 as a template for full atomistic molecular dynamics (MD) simulations, we reveal that the SF in down-state undergoes inactivation via conformational changes, while the up-state structure maintains a stable and conductive SF. This suggests an atomistic mechanism for the low channel activity previously assigned to the down state, but not evident from the crystal structure. Furthermore, experimentally by using (de-)phosphorylation mimics and chemically attaching lipid tethers to the proximal C-terminus (pCt), we confirm the hypothesis that moving the pCt towards the membrane induces the up-state. Based on MD simulations, we propose two gating pathways by which movement of the pCt controls the stability (i.e., conductivity) of the filter gate. Together, these findings provide atomistic insights into the SF gating mechanism and the physiological regulation of TREK channels by phosphorylation.

Discovery of Salifungin as a Repurposed Antibiotic against Methicillin-Resistant Staphylococcus aureus with Limited Resistance Development

Exploring novel antimicrobial drugs and strategies has become essential to the fight MRSA-associated infections. Herein, we found that membrane-disrupted repurposed antibiotic salifungin had excellent bactericidal activity against MRSA, with limited development of drug resistance. Furthermore, adding salifungin effectively decreased the minimum inhibitory concentrations of clinical antibiotics against Staphylococcus aureus. Evaluations of the mechanism demonstrated that salifungin disrupted the level of H+ and K+ ions using hydrophilic and lipophilic groups to interact with bacterial membranes, causing the disruption of bacterial proton motive force followed by impacting on bacterial the function of the respiratory chain and adenosine 5'-triphosphate, thereby inhibiting phosphatidic acid biosynthesis. Moreover, salifungin also significantly inhibited the formation of bacterial biofilms and eliminated established bacterial biofilms by interfering with bacterial membrane potential and inhibiting biofilm-associated gene expression, which was even better than clinical antibiotics. Finally, salifungin exhibited efficacy comparable to or even better than that of vancomycin in the MRSA-infected animal models. In conclusion, these results indicate that salifungin can be a potential drug for treating MRSA-associated infections.

Molecular dynamics simulations of intracellular lipid droplets: a new tool in the toolbox

Lipid droplets (LDs) are ubiquitous intracellular organelles with a central role in multiple lipid metabolic pathways. However, identifying correlations between their structural properties and their biological activity has proved challenging, owing to their unique physicochemical properties as compared with other cellular membranes. In recent years, molecular dynamics (MD) simulations, a computational methodology allowing the accurate description of molecular assemblies down to their individual components, have been demonstrated to be a useful and powerful approach for studying LD structural and dynamical properties. In this short review, we attempt to highlight, as comprehensively as possible, how MD simulations have contributed to our current understanding of multiple molecular mechanisms involved in LD biology.

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.

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Molecular Simulations Reveal the Free Energy Landscape and Transition State of Membrane Electroporation

The formation of pores over lipid membranes by the application of electric fields, termed membrane electroporation, is widely used in biotechnology and medicine to deliver drugs, vaccines, or genes into living cells. Continuum models for describing the free energy landscape of membrane electroporation were proposed decades ago, but they have never been tested against spatially detailed atomistic models. Using molecular dynamics (MD) simulations with a recently proposed reaction coordinate, we computed potentials of mean force of pore nucleation and pore expansion in lipid membranes at various transmembrane potentials. Whereas the free energies of pore expansion are compatible with previous continuum models, the experimentally important free energy barrier of pore nucleation is at variance with established models. The discrepancy originates from different geometries of the transition state; previous continuum models assumed the presence of a membrane-spanning defect throughout the process, whereas, according to the MD simulations, the transition state of pore nucleation is typically passed before a transmembrane defect has formed. A modified continuum model is presented that qualitatively agrees with the MD simulations. Using kinetics of pore opening together with transition state theory, our free energies of pore nucleation are in excellent agreement with previous experimental data.

Computational Insights into the Interaction of the Conserved Cysteine-Noose Domain of the Human Respiratory Syncytial Virus G Protein with the Canonical Fractalkine Binding site of Transmembrane Receptor CX3CR1 Isoforms

The human Respiratory Syncytial Virus (hRSV) stands as one of the most common causes of acute respiratory diseases. The infectivity of this virus is intricately linked to its membrane proteins, notably the attachment glycoprotein (G protein). The latter plays a key role in facilitating the attachment of hRSV to respiratory tract epithelial cells, thereby initiating the infection process. The present study aimed to characterize the interaction of the conserved cysteine-noose domain of hRSV G protein (cndG) with the transmembrane CX3C motif chemokine receptor 1 (CX3CR1) isoforms using computational tools of molecular modeling, docking, molecular dynamics simulations, and binding free energy calculations. From MD simulations of the molecular system embedded in the POPC lipid bilayer, we showed a stable interaction of cndG with the canonical fractalkine binding site in the N-terminal cavity of the CX3CR1 isoforms and identified that residues in the extracellular loop 2 (ECL2) region and Glu279 of this receptor are pivotal for the stabilization of CX3CR1/cndG binding, corroborating what was reported for the interaction of the chemokine fractalkine with CX3CR1 and its structure homolog US28. Therefore, the results presented here contribute by revealing key structural points for the CX3CR1/G interaction, allowing us to better understand the biology of hRSV from its attachment process and to develop new strategies to combat it.

Exploring Biginelli-based scaffolds as A2B adenosine receptor antagonists: Unveiling novel structure-activity relationship trends, lead compounds, and potent colorectal anticancer agents

Antagonists of the A2B adenosine receptor have recently emerged as targeted anticancer agents and immune checkpoint inhibitors within the realm of cancer immunotherapy. This study presents a comprehensive evaluation of novel Biginelli-assembled pyrimidine chemotypes, including mono-, bi-, and tricyclic derivatives, as A2BAR antagonists. We conducted a comprehensive examination of the adenosinergic profile (both binding and functional) of a large compound library consisting of 168 compounds. This approach unveiled original lead compounds and enabled the identification of novel structure-activity relationship (SAR) trends, which were supported by extensive computational studies, including quantum mechanical calculations and free energy perturbation (FEP) analysis. In total, 25 molecules showed attractive affinity (Ki < 100 nM) and outstanding selectivity for A2BAR. From these, five molecules corresponding to the new benzothiazole scaffold were below the Ki < 10 nM threshold, in addition to a novel dual A2A/A2B antagonist. The most potent compounds, and the dual antagonist, showed enantiospecific recognition in the A2BAR. Two A2BAR selective antagonists and the dual A2AAR/A2BAR antagonist reported in this study were assessed for their impact on colorectal cancer cell lines. The results revealed a significant and dose-dependent reduction in cell proliferation. Notably, the A2BAR antagonists exhibited remarkable specificity, as they did not impede the proliferation of non-tumoral cell lines. These findings support the efficacy and potential that A2BAR antagonists as valuable candidates for cancer therapy, but also that they can effectively complement strategies involving A2AAR antagonism in the context of immune checkpoint inhibition.

A Methylazanediyl Bisacetamide Derivative Sensitizes Staphylococcus aureus Persisters to a Combination of Gentamicin And Daptomycin

Infections caused by Staphylococcus aureus, notably methicillin-resistant S. aureus (MRSA), pose treatment challenges due to its ability to tolerate antibiotics and develop antibiotic resistance. The former, a mechanism independent of genetic changes, allows bacteria to withstand antibiotics by altering metabolic processes. Here, a potent methylazanediyl bisacetamide derivative, MB6, is described, which selectively targets MRSA membranes over mammalian membranes without observable resistance development. Although MB6 is effective against growing MRSA cells, its antimicrobial activity against MRSA persisters is limited. Nevertheless, MB6 significantly potentiates the bactericidal activity of gentamicin against MRSA persisters by facilitating gentamicin uptake. In addition, MB6 in combination with daptomycin exhibits enhanced anti-persister activity through mutual reinforcement of their membrane-disrupting activities. Crucially, the "triple" combination of MB6, gentamicin, and daptomycin exhibits a marked enhancement in the killing of MRSA persisters compared to individual components or any double combinations. These findings underscore the potential of MB6 to function as a potent and selective membrane-active antimicrobial adjuvant to enhance the efficacy of existing antibiotics against persister cells. The molecular mechanisms of MB6 elucidated in this study provide valuable insights for designing anti-persister adjuvants and for developing new antimicrobial combination strategies to overcome the current limitations of antibiotic treatments.

Mathematical and computational modeling of fats and triacylglycerides

Fats and oils are found in many food products; however, their macroscopic properties are difficult to predict, especially when blending different fats or oils together. With difficulties in sourcing specific fats or oils, whether due to availability or pricing, food companies may be required to find alternative sources for these ingredients, with possible differences in ingredient performance. Mathematical and computational modeling of these ingredients can provide a quick way to predict their properties, avoiding costly trials or manufacturing problems, while, most importantly, keeping the consumers happy. This review covers a range of mathematical models for triacylglycerides (TAGs) and fats, namely, models for the prediction of melting point, solid fat content, and crystallization temperature and composition. There are a number of models that have been designed for both TAGs and fats and which have been shown to agree very well with empirical measurements, using both kinetic and thermodynamic approaches, with models for TAGs being used to, in turn, predict fat properties. The last section describes computational models to simulate the behavior of TAGs using molecular dynamics (MD). Simulation of TAGs using MD, however, is still at an early stage, although the most recent papers on this topic are bringing this area up to speed.

Cholesterols Induced Distinctive Entry of the Graphene Nanosheet into the Cell Membrane

Graphene nanosheets are highly valued in the biomedical field due to their potential applications in drug delivery, biological imaging, and biosensors. Their biological effects on mammalian cells may be influenced by cholesterols, which are crucial components in cell membranes that take part in many vital processes. Therefore, it is particularly important to investigate the effect of cholesterols on the transport mechanism of graphene nanosheets in the cell membrane as well as the final stable configuration of graphene, which may have an impact on cytotoxicity. In this paper, the molecular details of a graphene nanosheet interacting with a 1,2-dipalmitoyl-sn-glycero-3-phosphorylcholine (DPPC) membrane with cholesterols were studied using molecular dynamics simulations. Results showed that the structure of the graphene nanosheet transits from the cut-in state in a pure DPPC membrane to being sandwiched between two DPPC leaflets when cholesterols reach a certain concentration. The underlying mechanism showed that cholesterols are preferentially adsorbed on the graphene nanosheet, which causes a larger disturbance to the nearby DPPC tails and thus guides the graphene nanosheet into the core of lipid bilayers to form a sandwiched structure. Our results are helpful for understanding the fundamental interaction mechanism between the graphene nanosheet and cell membrane and to explore the potential applications of the graphene nanosheet in biomedical sciences.

N-(Heteroaryl)thiophene sulfonamides as angiotensin AT2 receptor ligands

Two series of N-(heteroaryl)thiophene sulfonamides, encompassing either a methylene imidazole group or a tert-butylimidazolylacetyl group in the meta position of the benzene ring, have been synthesized. An AT2R selective ligand with a Ki of 42 nM was identified in the first series and in the second series, six AT2R selective ligands with significantly improved binding affinities and Ki values of <5 nM were discovered. The binding modes to AT2R were explored by docking calculations combined with molecular dynamics simulations. Although some of the high affinity ligands exhibited fair stability in human liver microsomes, comparable to that observed with C21 undergoing clinical trials, most ligands displayed a very low metabolic stability with t½ of less than 10 min in human liver microsomes. The most promising ligand, with an AT2R Ki value of 4.9 nM and with intermediate stability in human hepatocytes (t½ = 77 min) caused a concentration-dependent vasorelaxation of pre-contracted mouse aorta.

Discovery of lipid binding sites in a ligand-gated ion channel by integrating simulations and cryo-EM

Ligand-gated ion channels transduce electrochemical signals in neurons and other excitable cells. Aside from canonical ligands, phospholipids are thought to bind specifically to the transmembrane domain of several ion channels. However, structural details of such lipid contacts remain elusive, partly due to limited resolution of these regions in experimental structures. Here, we discovered multiple lipid interactions in the channel GLIC by integrating cryo-electron microscopy and large-scale molecular simulations. We identified 25 bound lipids in the GLIC closed state, a conformation where none, to our knowledge, were previously known. Three lipids were associated with each subunit in the inner leaflet, including a buried interaction disrupted in mutant simulations. In the outer leaflet, two intrasubunit sites were evident in both closed and open states, while a putative intersubunit site was preferred in open-state simulations. This work offers molecular details of GLIC-lipid contacts particularly in the ill-characterized closed state, testable hypotheses for state-dependent binding, and a multidisciplinary strategy for modeling protein-lipid interactions.

The rise of FTIR spectroscopy in the characterization of asymmetric lipid membranes

In contrast to symmetric unilamellar liposomes (sLUVs) prepared from a mixture of different lipids, asymmetric ones (aLUVs) with different lipid composition in the inner and outer membrane leaflets are more suitable model systems of eukaryotic plasma membranes. However, apart from the challenging preparation of asymmetric liposomes and small amounts of obtained asymmetric unilamellar liposomes (aLUVs), a major drawback is the qualitative characterization of asymmetry, as each of the techniques used so far has certain limitations. In this regard, we prepared aLUVs composed dominantly of DPPC(out)/DPPS(in) lipids and, along with 1H NMR and DSC characterization, we showed for the first time how FTIR spectroscopy can be used in the presence of (a)symmetry between DPPC/DPPS lipid bilayers. Using second derivative FTIR spectra we demonstrated not only that the hydration of lipids glycerol backbone and choline moiety of DPPC differs in s/aLUVs, but in addition that the lateral interactions between hydrocarbon chains during the phase change display different trend in s/aLUVs. Molecular dynamics simulations confirmed different chain ordering and packing between s/a bilayers, with a significant influence of temperature, i.e. membrane phase.

Unraveling the Molecular Mechanisms of Alcohol-Mediated Skin Permeation Enhancement: Insights from Molecular Dynamics Simulations

The application of alcohols as permeation enhancers in pharmaceutical and cosmetic formulations has attracted considerable attention, owing to their skin permeation-enhancing effect. Nonetheless, the elucidation of the fundamental mechanisms underlying the skin permeation-enhancing effect remains elusive. In this study, molecular dynamics (MD) simulations were employed to investigate the effect of 1,2-propanediol (1,2-PDO), 1,2-butanediol (1,2-BDO), and ethanol (EtOH) on the stratum corneum (SC) model membrane. The results showed that the effect of alcohols on the SC model membrane displayed a concentration-dependent nature. The alcohols can interact with SC lipids and exhibit a remarkable ability to selectively extract free fatty acid (FFA) molecules from the SC model membrane and make the SC looser. Meanwhile, 1,2-BDO and EtOH can penetrate into SC lipid bilayers at higher concentrations, leading to the formation of continuous hydrophilic defects in SC. The FFA extraction and the formation of continuous hydrophilic defects induced ceramide (CER) tail chains to become more disordered and fluid and also weakened the hydrogen bonding (H-bonding) network among SC lipids. Both the FFA extraction and the continuous hydrophilic defect formation endowed alcohols with the permeation-enhancing effect. The constrained simulations revealed that the free energy barriers decreased for the permeation of the hydrophilic model molecule (COL) across the SC model membranes containing alcohols, particularly for 1,2-BDO and EtOH. The possible permeation-enhancing mechanisms of alcohols were proposed correspondingly. This work not only provided a deep understanding of the transdermal permeation-enhancing behavior of alcohols at the molecular level but also provided necessary reference information for designing effective transdermal drug delivery systems in applications.

A molecular dynamics study of cell-penetrating peptide transportan-10 (TP10): Binding, folding and insertion to transmembrane state in zwitterionic membrane

Transportan 10 (TP10) is a 21-residue, cationic, α-helical cell-penetrating peptide that can be used as a delivery vector for various bioactive molecules. Based on recent confocal microscopy studies, it is believed that TP10 can translocate across neutral lipid membrane passively, possibly as a monomer, without the formation of permanent pore. Here, we performed extensive molecular dynamics (MD) simulations of TP10W (Y3W variant of TP10) to find the microscopic details of binding, folding and insertion of TP10W to transmembrane state in POPC bilayer. Binding study with CHARMM36 force field showed that TP10W initially binds to the membrane surface in unstructured configuration, but it spontaneously folds into α-helical conformation under the lipid head groups. Further insertion of TP10W, changing from a surface bound state to a vertically oriented transmembrane state, was investigated via umbrella simulations. The resulting free energy profile shows a relatively small barrier between two states, suggesting a possible translocation pathway as a monomer. In fact, unbiased simulation of transmembrane TP10W revealed how a charged Lys side chain can move from one leaflet to the other without a significant free energy cost. Finally, we compared the results of TP10W simulations with those of point mutated variants (TP10W-K12A18 and TP10W-K19L) to understand the effect of charge distribution on the peptide. It was observed that such a conservative mutation can cause noticeable changes in the conformations of both surface bound and transmembrane states. The results of present study will be discussed in relation to the experimentally observed activities of TP10W against neutral membrane.

Generalized polarization and time-resolved fluorescence provide evidence for different populations of Laurdan in lipid vesicles

The solvatochromic dye Laurdan is widely used in sensing the lipid packing of both model and biological membranes. The fluorescence emission maximum shifts from about 440 nm (blue channel) in condensed membranes (So) to about 490 nm (green channel) in the liquid-crystalline phase (Lα). Although the fluorescence intensity based generalized polarization (GP) is widely used to characterize lipid membranes, the fluorescence lifetime of Laurdan, in the blue and the green channel, is less used for that purpose. Here we explore the correlation between GP and fluorescence lifetimes by spectroscopic measurements on the So and Lα phases of large unilamellar vesicles of DMPC and DPPC. A positive correlation between GP and the lifetimes is observed in each of the optical channels for the two lipid phases. Microfluorimetric determinations on giant unilamellar vesicles of DPPC and DOPC at room temperature are performed under linearly polarized two-photon excitation to disentangle possible subpopulations of Laurdan at a scale below the optical resolution. Fluorescence intensities, GP and fluorescence lifetimes depend on the angle between the orientation of the linear polarization of the excitation light and the local normal to the membrane of the optical cross-section. This angular variation depends on the lipid phase and the emission channel. GP and fluorescence intensities in the blue and green channel in So and in the blue channel in Lα exhibit a minimum near 90o. Surprisingly, the intensity in the green channel in Lα reaches a maximum near 90o. The fluorescence lifetimes in the two optical channels also reach a pronounced minimum near 90o in So and Lα, apart from the lifetime in the blue channel in Lα where the lifetime is short with minimal angular variation. To our knowledge, these experimental observations are the first to demonstrate the existence of a bent conformation of Laurdan in lipid membranes, as previously suggested by molecular dynamics calculations.

Computational approaches to lipid-based nucleic acid delivery systems

Nucleic acid-based therapies have shown enormous effectiveness as vaccines against the recent COVID19 pandemics and hold great promises in the fight of a broad spectrum of diseases ranging from viral infections to cancer up to genetically transmitted pathologies. Due to their highly degradable polyanionic nature, nucleic acids need to be packed in sophisticate delivery vehicles which compact them up, protect them from early degradation and help delivery them to the right tissue/cells. Lipid-based nanoparticles (LNP) represent, at present, the main solution for nucleic acid delivery. They are made of a mixture of lipids whose key ingredient is an ionizable cationic lipid. Indeed, the interactions between the polyanionic nucleic acids and the ionizable cationic lipids, and their pH-dependent regulation in the life cycle of the nanoparticle, from production to cargo delivery, mostly determine the effectiveness of the therapeutic approach. Notwithstanding the large improvements in the delivery efficiency of LNPs in the last two decades, it is estimated that only a small fraction of the cargo is actually delivered, stimulating further research for the design of more effective LNP formulations. A rationally driven design would profit from the knowledge of the precise molecular structure of these materials, which is however still either missing or characterized by poor spatial resolution. Computational approaches have often been used as a molecular microscope either to enrich the available experimental data and provide a molecular-level picture of the LNPs or even simulate specific processes involving the formation and/or the molecular mechanisms of action of the LNP. Here, I review the recent literature in the field.

Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants

BACKGROUND

Many organisms rely on mineral nutrients taken directly from the soil or aquatic environment, and therefore, developed mechanisms to cope with the limitation of a given essential nutrient. For example, photosynthetic cells have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. Under phosphate starvation, phospholipids in extraplastidial membranes are replaced by betaine lipids in microalgae. In higher plants, the synthesis of betaine lipid is lost, driving plants to other strategies to cope with phosphate starvation where they replace their phospholipids by glycolipids.

RESULTS

The aim of this work was to evaluate to what extent betaine lipids and PC lipids share physicochemical properties and could substitute for each other. By neutron diffraction experiments and dynamic molecular simulation of two synthetic lipids, the dipalmitoylphosphatidylcholine (DPPC) and the dipalmitoyl-diacylglyceryl-N,N,N-trimethylhomoserine (DP-DGTS), we found that DP-DGTS bilayers are thicker than DPPC bilayers and therefore are more rigid. Furthermore, DP-DGTS bilayers are more repulsive, especially at long range, maybe due to unexpected unscreened electrostatic contribution. Finally, DP-DGTS bilayers could coexist in the gel and fluid phases.

CONCLUSION

The different properties and hydration responses of PC and DGTS provide an explanation for the diversity of betaine lipids observed in marine organisms and for their disappearance in seed plants.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Effect of quaternary ammonium surfactants on biomembranes using molecular dynamics simulation

Research conducted both prior to and after the emergence of the COVID-19 pandemic reveals a notable rise in human exposure to cleaning products, hand sanitizers, and personal care items. Moreover, there has been a corresponding increase in the environmental release of these chemicals. Cleaning and disinfecting products often contain quaternary ammonium compounds (QACs) with alkyl chains as long as 8-12 carbon atoms. The attachment of quaternary ammonium surfactants to the membrane resulted in the deformation of the bilayer and membrane disruption. Before interactions with cell membranes, these surfactant molecules may form different aggregates depending on their architecture. Interaction of surfactant monomers or clusters with the cell membrane changes the physiochemical properties of the biomembranes. To investigate this interaction and its influence on membrane properties, we conducted molecular dynamics simulations of cationic quaternary ammonium surfactants interacting with dipalmitoylphosphatidylcholine (DPPC) membranes. Our simulations revealed significant interactions between the surfactants and the phospholipids, leading to substantial alterations in the structure of the bilayer. The results are compared with the simulated anionic (SDS) and nonionic surfactants/bilayer systems. Various aspects were considered, including the aggregation process, migration behavior, and eventual equilibrium of these molecules at the interface between the membrane and water. This analysis used various techniques such as density profiles, distribution functions, cluster analysis, order parameters, hydrogen bonding (H-bonding), and mean-square displacements. The results indicate that while surfactants with shorter alkyl tails (N-(2-hydroxyethyl)-N,N-dimethyloctan-1-aminium chloride (HEDMOAC)) make strong hydrogen bonds with the phosphate group and ester oxygen of the phosphatidylcholine bilayer and enter toward the bilayer in the monomer form, surfactants of longer alkyl tails aggregated on the membrane head-water interface and interact minimally with the head groups of the DPPC bilayer. For DDEDMEAC, a quaternary ammonium surfactant with a hydrophobic alkyl chain consisting of two decanoate groups, alteration of the structural and dynamical properties of the bilayer is expected to be governed by two different factors. First, the structural order of DPPC increases as surfactant aggregates interact with the membrane head group. Second, the decrease in the order of the bilayer occurs due to the insertion of surfactant monomers within the hydrophobic region of the bilayer. Strong interactions between constituents of tetraoctylammonium bromide (TOABr) and lipid head groups lead to a reduction in interlipid interactions and order, which further results in increased porosity of cellular membranes. Understanding the extent of these interactions plays a pivotal role in the toxicological assessment of these surfactants.

Head-to-head comparison of in vitro and in vivo efficacy of pHLIP-conjugated anti-seed gamma peptide nucleic acids

Gamma peptide nucleic acids (γPNAs) have recently garnered attention in diverse therapeutic and diagnostic applications. Serine and diethylene-glycol-containing γPNAs have been tested for numerous RNA-targeting purposes. Here, we comprehensively evaluated the in vitro and in vivo efficacy of pH-low insertion peptide (pHLIP)-conjugated serine and diethylene-based γPNAs. pHLIP targets only the acidic tumor microenvironment and not the normal cells. We synthesized and parallelly tested pHLIP-serine γPNAs and pHLIP-diethylene glycol γPNAs that target the seed region of microRNA-155, a microRNA that is upregulated in various cancers. We performed an all-atom molecular dynamics simulation-based computational study to elucidate the interaction of pHLIP-γPNA constructs with the lipid bilayer. We also determined the biodistribution and efficacy of the pHLIP constructs in the U2932-derived xenograft model. Overall, we established that the pHLIP-serine γPNAs show superior results in vivo compared with the pHLIP-diethylene glycol-based γPNA.

Alkali metal cations modulate the geometry of different binding sites in HCN4 selectivity filter for permeation or block

Hyperpolarization-activated cyclic-nucleotide gated (HCN) channels are important for timing biological processes like heartbeat and neuronal firing. Their weak cation selectivity is determined by a filter domain with only two binding sites for K+ and one for Na+. The latter acts as a weak blocker, which is released in combination with a dynamic widening of the filter by K+ ions, giving rise to a mixed K+/Na+ current. Here, we apply molecular dynamics simulations to systematically investigate the interactions of five alkali metal cations with the filter of the open HCN4 pore. Simulations recapitulate experimental data like a low Li+ permeability, considerable Rb+ conductance, a block by Cs+ as well as a punch through of Cs+ ions at high negative voltages. Differential binding of the cation species in specific filter sites is associated with structural adaptations of filter residues. This gives rise to ion coordination by a cation-characteristic number of oxygen atoms from the filter backbone and solvent. This ion/protein interplay prevents Li+, but not Na+, from entry into and further passage through the filter. The site equivalent to S3 in K+ channels emerges as a preferential binding and presumably blocking site for Cs+. Collectively, the data suggest that the weak cation selectivity of HCN channels and their block by Cs+ are determined by restrained cation-generated rearrangements of flexible filter residues.

Structural insights on the deformations induced by various mutations on cholesteryl ester transfer protein

Cholesteryl Ester Transfer Protein (CETP) is a plasma glycoprotein that intervenes the reverse cholesterol transport (RCT) by equimolar exchange of Cholesteryl esters (CE) and Triglycerides (TGs) between anti-atherogenic High-Density Lipoproteins (HDLs) and pro-atherogenic Low-Density Lipoproteins (LDLs) resulting in the increased concentration of CEs in LDL. This is a potential cause for the formation of atherosclerotic plaques in blood vessels leading to fatality. Therefore, blocking the function of CETP has emerged as a novel strategy for suppressing atherosclerotic plaques. The crystal structure of CETP revealed two Cholesteryl esters (CEs) in the hydrophobic tunnel and two phospholipids (PLs) plugged on the concave surface. Previous lipid transfer assay experimental studies have shown a substantial reduction in the neutral lipid transfer in [R201S] and [I443W, V198W] mutants. However, the protein conformational arrangements due to the mutations present in the CETP system leading to a decrease in the transfer rate of neutral lipids is not explored. Thus, I explored the reason behind the decreased transfer rate in mutants using molecular dynamics (MD) simulations and free energy calculations. Resulting evidences show that R201S mutant induces unfavorable bending angle to CETP with a decreased binding efficiency between N-terminal phospholipid of CETP with S201. Also, an unfavorable conformation state of TGs is formed which makes them difficult to transfer across CETP. Likewise, [I443W, V198W] mutant induces unfavorable CE, TG, and bending angle conformation to CETP impeding neutral lipid transfer. Thus, my results provide sufficient insights on the causation for a decreased transfer rate as reported earlier. The detailed understanding obtained here could help in developing a new strategy in preventing the function of CETP by blocking the role of potential hot spot residues.

The Membrane-Mediated Interaction of Liquid-Ordered Lipid Domains in the Presence of Amphipathic Peptides

The lipid membranes of living cells are composed of a large number of lipid types and can undergo phase separation with the formation of nanometer-scale liquid-ordered lipid domains, also called rafts. Raft coalescence, i.e., the fusion of lipid domains, is involved in important cell processes, such as signaling and trafficking. In this work, within the framework of the theory of elasticity of lipid membranes, we explore how amphipathic peptides adsorbed on lipid membranes may affect the domain-domain fusion processes. We show that the elastic deformations of lipid membranes drive amphipathic peptides to the boundary of lipid domains, which leads to an increase in the average energy barrier of the domain-domain fusion, even if the surface concentration of amphipathic peptides is low and the domain boundaries are only partially occupied by the peptides. This inhibition of the fusion of lipid domains may lead to negative side effects of using amphipathic peptides as antimicrobial agents.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

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

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

Viroporins of Mpox Virus

Mpox or monkeypox virus (MPXV) belongs to the subclass of Poxviridae and has emerged recently as a global threat. With a limited number of anti-viral drugs available for this new virus species, it is challenging to thwart the illness it begets. Therefore, characterizing new drug targets in the virus may prove advantageous to curbing the disease. Since channels as a family are excellent drug targets, we have sought to identify viral ion channels for this virus, which are instrumental in formulating channel-blocking anti-viral drugs. Bioinformatics analyses yielded eight transmembranous proteins smaller or equal to 100 amino acids in length. Subsequently, three independent bacteria-based assays have pointed to five of the eight proteins that exhibit ion channel activity. Finally, we propose a tentative structure of four ion channels from their primary amino acid sequences, employing AlphaFold2 and molecular dynamic simulation methods. These results may represent the first steps in characterizing MPXV viroporins en route to developing blockers that inhibit their function.

Structural and thermodynamic properties of bulk triglycerides and triglyceride/water mixtures reproduced using a polarizable coarse-grained model

Triglycerides (TGs) play important roles in renewable energies, food production, medicine, and metabolism in organisms. Here, we developed a novel coarse-grained (CG) force field (FF) for triglycerides to reproduce both the structural and thermodynamic properties of bulk TGs, TG/air interfaces, and TG/water mixtures using molecular dynamics (MD) simulations. We rigorously optimized the bonded and nonbonded force parameters between the CG beads of TGs and nonbonded force parameters between TG beads and polarizable CG water beads by employing an efficient meta-multilinear interpolation parameterization algorithm recently developed by us. This CG FF performs very well in reproducing the percolating network of the TG bulk phase self-assembled in water and a variety of molecular conformations predicted by all-atom MD simulations. More importantly, it also correctly reproduces multiple experimentally measurable macroscopic thermodynamic properties, including the density and surface tensions of both the TG/air and TG/water interfaces. This paves the way for studying more complicated systems involving TGs on a large scale.

The opening dynamics of the lateral gate regulates the activity of rhomboid proteases

Rhomboid proteases hydrolyze substrate helices within the lipid bilayer to release soluble domains from the membrane. Here, we investigate the mechanism of activity regulation for this unique but wide-spread protein family. In the model rhomboid GlpG, a lateral gate formed by transmembrane helices TM2 and TM5 was previously proposed to allow access of the hydrophobic substrate to the shielded hydrophilic active site. In our study, we modified the gate region and either immobilized the gate by introducing a maleimide-maleimide (M2M) crosslink or weakened the TM2/TM5 interaction network through mutations. We used solid-state nuclear magnetic resonance (NMR), molecular dynamics (MD) simulations, and molecular docking to investigate the resulting effects on structure and dynamics on the atomic level. We find that variants with increased dynamics at TM5 also exhibit enhanced activity, whereas introduction of a crosslink close to the active site strongly reduces activity. Our study therefore establishes a strong link between the opening dynamics of the lateral gate in rhomboid proteases and their enzymatic activity.

Eukaryotic Cell Membranes: Structure, Composition, Research Methods and Computational Modelling

This paper deals with the problems encountered in the study of eukaryotic cell membranes. A discussion on the structure and composition of membranes, lateral heterogeneity of membranes, lipid raft formation, and involvement of actin and cytoskeleton networks in the maintenance of membrane structure is included. Modern methods for the study of membranes and their constituent domains are discussed. Various simplified models of biomembranes and lipid rafts are presented. Computer modelling is considered as one of the most important methods. This is stated that from the study of the plasma membrane structure, it is desirable to proceed to the diverse membranes of all organelles of the cell. The qualitative composition and molar content of individual classes of polar lipids, free sterols and proteins in each of these membranes must be considered. A program to create an open access electronic database including results obtained from the membrane modelling of individual cell organelles and the key sites of the membranes, as well as models of individual molecules composing the membranes, has been proposed.

Permeation of a Homologous Series of NBD-Labeled Fatty Amines through Lipid Bilayers: A Molecular Dynamics Study

Permeation through biomembranes is ubiquitous for drugs to reach their active sites. Asymmetry of the cell plasma membrane (PM) has been described as having an important role in this process. Here we describe the interaction of a homologous series of 7-nitrobenz-2-oxa-1,3-diazol-4-yl (NBD)-labeled amphiphiles (NBD-Cn, n = 4 to 16) with lipid bilayers of different compositions (1-palmitoyl, 2-oleoyl-sn-glycero-3-phosphocholine (POPC):cholesterol (1:1) and palmitoylated sphingomyelin (SpM):cholesterol (6:4)), including an asymmetric bilayer. Both unrestrained and umbrella sampling (US) simulations (at varying distances to the bilayer center) were carried out. The free energy profile of NBD-Cn at different depths in the membrane was obtained from the US simulations. The behavior of the amphiphiles during the permeation process was described regarding their orientation, chain elongation, and H-bonding to lipid and water molecules. Permeability coefficients were also calculated for the different amphiphiles of the series, using the inhomogeneous solubility-diffusion model (ISDM). Quantitative agreement with values obtained from kinetic modeling of the permeation process could not be obtained. However, for the longer, and more hydrophobic amphiphiles, the variation trend along the homologous series was qualitatively better matched by the ISDM when the equilibrium location of each amphiphile was taken as reference (ΔG = 0), compared to the usual choice of bulk water.

Biophysical Analysis to Assess the Interaction of CRAC and CARC Motif Peptides of Alpha Hemolysin of Escherichia coli with Membranes

Alpha hemolysin of Escherichia coli (HlyA) is a pore-forming protein, which is a prototype of the "Repeat in Toxins" (RTX) family. It was demonstrated that HlyA-cholesterol interaction facilitates the insertion of the toxin into membranes. Putative cholesterol-binding sites, called cholesterol recognition/amino acid consensus (CRAC), and CARC (analogous to CRAC but with the opposite orientation) were identified in the HlyA sequence. In this context, two peptides were synthesized, one derived from a CARC site from the insertion domain of the toxin (residues 341-353) (PEP 1) and the other one from a CRAC site from the domain between the acylated lysines (residues 639-644) (PEP 2), to study their role in the interaction of HlyA with membranes. The interaction of peptides with membranes of different lipid compositions (pure POPC and POPC/Cho of 4:1 and 2:1 molar ratios) was analyzed by surface plasmon resonance and molecular dynamics simulations. Results demonstrate that both peptides interact preferentially with Cho-containing membranes, although PEP 2 presents a lower K_D than PEP 1. Molecular dynamics simulation results indicate that the insertion and interaction of PEP 2 with Cho-containing membranes are more prominent than those caused by PEP 1. The hemolytic activity of HlyA in the presence of peptides indicates that PEP 2 was the only one that inhibits HlyA activity, interfering in the binding between the toxin and cholesterol.

Bound Phospholipids Assist Cholesteryl Ester Transfer in the Cholesteryl Ester Transfer Protein

Cholesteryl ester transfer protein (CETP) is a plasma glycoprotein that assists the transfer of cholesteryl esters (CEs) from antiatherogenic high-density lipoproteins (HDLs) to proatherogenic low-density lipoproteins (LDLs), initiating cholesterol plaques in the arteries. Consequently, inhibiting the activity of CETP is therefore being pursued as a novel strategy to reduce the risk of cardiovascular diseases (CVDs). The crystal structure of CETP has revealed the presence of two CEs running in the hydrophobic tunnel and two plugged-in phospholipids (PLs) near the concave surface. Other than previous animal models that rule out the PL transfer by CETP and PLs in providing the structural stability, the functional importance of bound phospholipids in CETP is not fully explored. Here, we employ a series of molecular dynamics (MD) simulations, steered molecular dynamics (SMD) simulations, and free energy calculations to unravel the effect of PLs on the functionality of the protein. Our results suggest that PLs play an important role in the transfer of neutral lipids by transforming the unfavorable bent conformation of CEs into a favorable linear conformation to facilitate the smooth transfer. The results also suggest that the making and breaking interactions of the hydrophobic tunnel residues with CEs with a combined effort from PLs are responsible for the transfer of CEs. Further, the findings demonstrate that the N-PL has a more pronounced effort on CE transfer than C-PL but efforts from both PLs are essential in the transfer. Thus, we propose that the functionally important PLs can be considered with potential research interest in targeting cardiovascular diseases.

Fifty Years of Biophysics at the Membrane Frontier

The author first describes his childhood in the South and the ways in which it fostered the values he has espoused throughout his life, his development of a keen fascination with science, and the influences that supported his progress toward higher education. His experiences in ROTC as a student, followed by two years in the US Army during the Vietnam War, honed his leadership skills. The bulk of the autobiography is a chronological journey through his scientific career, beginning with arrival at the University of California, Irvine in 1972, with an emphasis on the postdoctoral students and colleagues who have contributed substantially to each phase of his lab's progress. White's fundamental findings played a key role in the development of membrane biophysics, helping establish it as fertile ground for research. A story gradually unfolds that reveals the deeply collaborative and painstakingly executed work necessary for a successful career in science.

Equivalence of Charge Imbalance and External Electric Fields during Free Energy Calculations of Membrane Electroporation

Electric fields across lipid membranes play important roles in physiology, medicine, and biotechnology, rationalizing the wide interest in modeling transmembrane potentials in molecular dynamics simulations. Transmembrane potentials have been implemented with external electric fields or by imposing charge imbalance between the two water compartments of a stacked double-membrane system. We compare the two methods in the context of membrane electroporation, which involves a large change of membrane structure and capacitance. We show that, given that Ewald electrostatics are defined with tinfoil boundary conditions, the two methods lead to (i) identical potentials of mean force (PMFs) of pore formation and expansion at various potentials, demonstrating that the two methods impose equivalent driving forces for large-scale transitions at membranes, and (ii) to identical polarization of water within thin water wires or open pores, suggesting that the two methods furthermore impose equivalent local electric fields. Without tinfoil boundary conditions, effects from external fields on pore formation are spuriously suppressed or even removed. Together, our study shows that both methods, external fields and charge imbalance, are well suitable for studying large-scale transitions of lipid membranes that involve changes of membrane capacitance. However, using charge imbalance is technically more challenging for maintaining a constant transmembrane potential since it requires updating of the charge imbalance as the membrane capacitance changes.

The effects of adsorbed benzo(a)pyrene on dynamic behavior of polystyrene nanoplastics through phospholipid membrane: A molecular simulation study

Nanoplastics (NPs) are mainly generated from the decomposition of plastic waste and industrial production, which have attracted much attention due to the potential risk for humans. The ability of NPs to penetrate different biological barriers has been proved, but the understanding of molecular details is very limited, especially for organic pollutant-NP combinations. Here, we investigated the uptake process of polystyrene NPs (PSNPs) combined with benzo(a)pyrene (BAP) molecules by dipalmitoylphosphatidylcholine (DPPC) bilayers by molecular dynamics (MD) simulations. The results showed that the PSNPs can adsorb and accumulate BAP molecules in water phase and then carried BAP molecules to enter DPPC bilayers. At the same time, the adsorbed BAP promoted the penetration of PSNPs into DPPC bilayers effectively by hydrophobic effect. The process of BAP-PSNP combinations penetrating into DPPC bilayers can be summarized into four steps including adhesion on the DPPC bilayer surface, uptake by the DPPC bilayer, BAP molecules detached from the PSNPs, and the PSNPs depolymerized in the bilayer interior. Furthermore, the amount of adsorbed BAP on PSNPs affected the properties of DPPC bilayers directly, especially the fluidity of DPPC bilayers that determine the physiologic function. Obviously, the combined effect of PSNPs and BAP enhanced the cytotoxicity. This work not only presented a vivid transmembrane process of BAP-PSNP combinations and revealed the nature of the effects of adsorbed benzo(a)pyrene on the dynamic behavior of polystyrene nanoplastics through phospholipid membrane, but also provide some necessary information of the potential damage for organic pollutant-nanoplastic combinations on human health at a molecular level.

Interaction of guanidinium and ammonium cations with phosphatidylcholine and phosphatidylserine lipid bilayers - Calorimetric, spectroscopic and molecular dynamics simulations study

The ability of arginine-rich peptides to cross the lipid bilayer and enter cytoplasm, unlike their lysine-based analogues, is intensively studied in the context of cell-penetrating peptides. Although the experiments have not yet reconstructed their internalization mechanism, the computational studies have shown that the type or charge of lipid polar groups is one of the crucial factors in their translocation. In order to gain more detailed insight into the interaction of guanidinium (Gdm⁺) and ammonium (NH₄⁺) cations, as important building blocks in arginine and lysine amino acids, with lipid bilayers, we conducted the experimental and computational study that tackles this phenomenon. The adsorption of Gdm⁺ and NH₄⁺ on lipid bilayers prepared from a zwitterionic (DPPC) and an anionic (DPPS) lipid was examined by thermoanalytic and spectroscopic techniques. Using temperature-dependent UV-Vis spectroscopy and DSC calorimetry we determined the impact of Gdm⁺ and NH₄⁺ on the thermotropic properties of lipid bilayers. FTIR data, along with molecular dynamics simulations, unraveled the molecular-level details on the nature of their interactions, showing the proton transfer between NH₄⁺ and DPPS, but not between Gdm⁺ and DPPS. The findings originated from this work imply that Gdm⁺ and NH₄⁺ form qualitatively different interactions with lipids of different charge which is reflected in the physico-chemical interactions that arginine-and lysine-based peptides establish at a complex and chemically heterogeneous environment such as the biological membrane.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

In Silico Prediction of Stratum Corneum Partition Coefficients via COSMOmic and Molecular Dynamics Simulations

Stratum corneum (SC) is the main barrier of human skin where the inter-corneocytes lipids provide the main pathway for transdermal permeation of functional actives of skin care and health. Molecular dynamics (MD) has been increasingly used to simulate the SC lipid bilayer structure so that the barrier property and its affecting factors can be elucidated. Among reported MD simulation studies, solute partition in the SC lipids, an important parameter affecting SC permeability, has received limited attention. In this work, we combine MD simulation with COSMOmic to predict the partition coefficients of dermatologically relevant solutes in SC lipid bilayer. Firstly, we run MD simulations to obtain equilibrated SC lipid bilayers with different lipid types, compositions, and structures. Then, the simulated SC lipid bilayer structures are fed to COSMOmic to calculate the partition coefficients of the solutes. The results show that lipid types and bilayer geometries play a minor role in the predicted partition coefficients. For the more lipophilic solutes, the predicted results of solute partition in SC lipid bilayers agree well with reported experimental values of solute partition in extracted SC lipids. For the more hydrophilic molecules, there is a systematical underprediction. Nevertheless, the MD/COSMOmic approach correctly reproduces the phenomenological correlation between the SC lipid/water partition coefficients and the octanol/water partition coefficients. Overall, the results show that the MD/COSMOmic approach is a fast and valid method for predicting solute partitioning into SC lipids and hence supporting the assessment of percutaneous absorption of skin care ingredients, dermatological drugs as well as environmental pollutants.

Subunit gating resulting from individual protonation events in Kir2 channels

Inwardly rectifying potassium (Kir) channels play a critical role in stabilizing the membrane potential, thus controlling numerous physiological phenomena in multiple tissues. Channel conductance is activated by cytoplasmic modulators that open the channel at the 'helix bundle crossing' (HBC), formed by the coming together of the M2 helices from each of the four subunits, at the cytoplasmic end of the transmembrane pore. We introduced a negative charge at the bundle crossing region (G178D) in classical inward rectifier Kir2.2 channel subunits that forces channel opening, allowing pore wetting and free movement of permeant ions between the cytoplasm and the inner cavity. Single-channel recordings reveal a striking pH-dependent subconductance behavior in G178D (or G178E and equivalent Kir2.1[G177E]) mutant channels that reflects individual subunit events. These subconductance levels are well resolved temporally and occur independently, with no evidence of cooperativity. Decreasing cytoplasmic pH shifts the probability towards lower conductance levels, and molecular dynamics simulations show how protonation of Kir2.2[G178D] and, additionally, the rectification controller (D173) pore-lining residues leads to changes in pore solvation, K+ ion occupancy, and ultimately K+ conductance. While subconductance gating has long been discussed, resolution and explanation have been lacking. The present data reveals how individual protonation events change the electrostatic microenvironment of the pore, resulting in distinct, uncoordinated, and relatively long-lasting conductance states, which depend on levels of ion pooling in the pore and the maintenance of pore wetting. Gating and conductance are classically understood as separate processes in ion channels. The remarkable sub-state gating behavior of these channels reveals how intimately connected 'gating' and 'conductance' are in reality.