CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans

Freezing-driven ionic charge imbalance leads to pore formation and osmotic injury of lipid membranes

In this work, we reveal the fundamental mechanism controlling osmotic injury of lipid membranes under low-temperature preservation using an in-vitro membrane system under freezing temperature and the molecular dynamic simulations. The freezing-driven ionic charge imbalance is the major factor affecting the membrane conformation and causing the osmotic injury. Under freezing temperature, the ionic charge imbalance, originating from the preferential incorporation of anions into the growing ice crystals, results in membrane poration with the directional penetration of ion molecules. Subsequently, the osmotic efflux of water molecules through the pore causes cell dehydration, eventually leading to the lethal osmotic injury of lipid membranes during freezing. Moreover, we find a stark difference in tolerance to freezing and the times required for pore formation in membranes with different lipid compositions. Membranes enriched with cholesterol and anionic lipids exhibit increased resistance to freezing-induced osmotic injury, as the addition of cholesterol and anionic lipids in membranes delays the pore formation under freezing temperature. These findings advance in depth the molecular-level understanding of freezing injury on lipid membranes and provide an opportunity to develop an alternative strategy to protect diverse cells during preservation at subzero temperatures by regulating the composition of lipid membranes.

Computational modeling of the anti-inflammatory complexes of IL37

Interleukin (IL) 37 is an anti-inflammatory cytokine belonging to the IL1 protein family. Owing to its pivotal role in modulating immune responses, elucidating the IL37 complex structures holds substantial therapeutic promise for various autoimmune disorders and cancers. However, none of the structures of IL37 complexes have been experimentally characterized. This computational study aims to address this gap through molecular modeling and classical molecular dynamics simulations. We modeled all protein-protein complexes of IL37 using a range of methods from homology modeling to AlphaFold2 multimer predictions. Models that successfully recapitulated experimental features underwent further analysis through molecular dynamics simulations. As positive controls, binary and ternary complexes of IL18 from PDB were included for comparison. Several key findings emerged from the comparative analysis of IL37 and IL18 complexes. IL37 complexes exhibited higher mobility than the IL18 complexes. Simulations of the IL37-IL18Rα complex revealed altered receptor conformations capable of accommodating a dimeric IL37, with the N-terminal loop of IL37 contributing significantly to complex mobility. Additionally, the glycosyl chain on N297 of IL18Rα, which contours one edge of the cytokine binding surface, acted as a steric block against the N-terminal loop of IL37. Further, investigations into interactions between IL37 and IL18BP suggested that a binding mode homologous to IL18 was unstable for IL37, indicating an alternative binding mechanism. Altogether, this study accesses to the structure and dynamics of IL37 complexes, revealing the structural underpinnings of the IL37's modulatory effect on the IL18 signaling pathway.

Mechanistic studies of mycobacterial glycolipid biosynthesis by the mannosyltransferase PimE

Tuberculosis (TB), a leading cause of death among infectious diseases globally, is caused by Mycobacterium tuberculosis (Mtb). The pathogenicity of Mtb is largely attributed to its complex cell envelope, which includes a class of glycolipids called phosphatidyl-myo-inositol mannosides (PIMs). These glycolipids maintain the integrity of the cell envelope, regulate permeability, and mediate host-pathogen interactions. PIMs comprise a phosphatidyl-myo-inositol core decorated with one to six mannose residues and up to four acyl chains. The mannosyltransferase PimE catalyzes the transfer of the fifth PIM mannose residue from a polyprenyl phosphate-mannose (PPM) donor. This step contributes to the proper assembly and function of the mycobacterial cell envelope; however, the structural basis for substrate recognition and the catalytic mechanism of PimE remain poorly understood. Here, we present the cryo-electron microscopy (cryo-EM) structures of PimE from Mycobacterium abscessus in its apo and product-bound form. The structures reveal a distinctive binding cavity that accommodates both donor and acceptor substrates/products. Key residues involved in substrate coordination and catalysis were identified and validated via in vitro assays and in vivo complementation, while molecular dynamics simulations delineated access pathways and binding dynamics. Our integrated approach provides comprehensive insights into PimE function and informs potential strategies for anti-TB therapeutics.

Ion-Induced Dipole Interactions Matter in Metadynamics Simulation of Transition Metal Ion Transporters

Metal transporters play crucial roles in the homeostasis and detoxification of beneficial and toxic metals in the human body. Due to experimental limitations in studying some metal transporters, numerous simulation studies have been conducted to understand the mechanisms of metal transport. However, studying the transport of divalent metal ions across the plasma membrane by metal transporters has been challenging with traditional molecular dynamics (MD) simulations. The metal ions often become trapped inside the transporter due to encountering high energy barriers during the transport process. In this study, we combined a recently developed metadynamics setup, known as well-tempered (WT) volume-based MTD, with the 12-6-4 Lennard-Jones (LJ) model representing transition metal-His/Asp/Glu side chain interactions. We used this approach to investigate the mechanism of action of a Zrt-/Irt-like protein (ZIP) transporter and compared the results with simulations using standard 12-6 LJ parameters for the transition metal-His/Asp/Glu side chain interactions. Our results show that the 12-6-4 LJ model for transition metal-His/Asp/Glu side chain interactions samples conformational space more broadly than the standard 12-6 LJ model for the same interactions in MTD simulations, facilitating the sampling of states that are hard to reach with the standard 12-6 model within the same time scale. This is even more remarkable given the fact that the model is dominated by 12-6 LJ interactions for the majority of the system, while the transition metal-His/Asp/Glu side chain interactions are the only interactions using the 12-6-4 LJ model. Hence, a small subset of interactions significantly modifies the states sampled by the entire protein leading to a more frequent observation of the transport of the transition metal ion. Overall, using 12-6-4 LJ to model the transition metal-His/Asp/Glu side chain interactions increases the potential for discovering additional metastable states by enabling metal ions to traverse more freely along the defined transport pathways.

Lanthanides Gd and Tm Can Inhibit Carnitine/Acylcarnitine Transporter: Insights from All-Atoms Simulations

Recent experimental evidence highlighted the inhibition of carnitine/acylcarnitine carrier (CAC), an important mitochondrial transmembrane protein for living organisms, by the early lanthanide Pr3+. A possible explanation of such a behaviour was found in the preference of the cation for amino acids like aspartate and glutamate containing a carboxylate in the side chain, laying in the inter-membrane space. Interaction of the cation with these residues can cause halt the transfer of the protein's substrates between the matrix and cytoplasm thus opening to new scenarios concerning the CAC-metal interactions and its relative inhibition. In the present work, the panel of metals binding the CAC protein is predictively expanded including Gd3+ and Tm3+, selected as representative species of middle and late lanthanides, respectively. A more realistic membrane-containing model of the protein was built and the comparative analysis of the molecular dynamics (MD) simulations of CAC apo-form with its complexed systems, named CAC-Pr, CAC-Gd and CAC-Tm, was performed. The analysis of the trajectories revealed that the inhibition is caused by the coordination of D132 and E179 to the cations and that such interactions generate a reorganization of important salt-bridges inside the framework of CAC. In detail, MD simulations highlighted that a spontaneous conformational change from cytoplasmatic-state (c-state) to matrix-state (m-state) induced by cations and that, in this condition, the protein channel is occluded, thus explaining the inhibition.

Analysis of the impact of protein conformational dynamics and intermolecular interactions on water flux through TIP3;1 aquaporins of Zea mays L

The discovery of aquaporins (AQPs) in 1992 had a profound impact on our understanding of the mechanisms underlying the transport of substances across cell membranes. To further understand water mobilization through AQPs, this study focuses on the characterization of water flux through the TIP3;1 aquaporins of Zea mays L. using molecular dynamics. The primary objective is to elucidate how protein-water intermolecular interactions and protein conformational dynamics impact water mobility across the cell membrane. To conduct this analysis, the three-dimensional structure of TIP3;1 was modeled using AlphaFold2, from which the complete system was constructed. This system consisted of a homotetramer of TIP3;1 immersed in a fragment of cell membrane and solvated with water molecules and ions. Subsequently, molecular dynamics simulations were conducted for 90 ns, resulting in the determination of an osmotic permeability coefficient (pf) of 0.8172 ± 0.146 × 10-14 cm3 s-1. In general, the mobility of water along the single-file water channel is influenced by the complex interplay of protein conformational dynamics and hydrogen bonding. The conformational dynamics of the protein channel modify the pore radius available for the passage of water, which affects the frequency of protein-water interactions and consequently influences the mobility of water in the channel. This study contributes to our understanding of the molecular mechanisms by which AQP activity is modulated without involving changes in protein chemical composition.

Long-chain cationic gemini surfactants as drug retention adjuvant on liposomes. A methodological approach with atorvastatin

This study delves into the development and characterization of dipalmitoyl phosphatidylcholine (DPPC) liposomes incorporated with gemini surfactant (tetradecamethylene-1,14 bis(dimethyl tetradecyl ammonium bromide); 14-14-14) and atorvastatin, aimed at enhancing drug delivery efficiency for cardiovascular diseases. The integration of gemini surfactants into liposomes is investigated for its potential to improve atorvastatin encapsulation and retention, addressing the drug's poor water solubility and the limitations of conventional liposomal systems. Through a combination of dynamic light scattering (DLS), differential scanning calorimetry (DSC), and molecular dynamics (MD) simulations, the study reveals that the presence of gemini surfactants significantly reduces liposome size and polydispersity, indicative of a more uniform and potentially unilamellar structure. DSC analysis highlights a decrease in transition temperatures and an alteration in transition symmetry, suggesting enhanced stability and a favourable drug release profile at physiological temperatures. MD simulations provide insight into the internalization mechanism of gemini surfactants and atorvastatin within the liposomal bilayer, demonstrating their mutual incorporation facilitated by polar interactions. Spectrophotometry-based retention studies further confirmed that liposomes containing gemini surfactants exhibit superior atorvastatin retention capabilities, nearly doubling the encapsulation efficiency compared to conventional liposomes. This research highlights the promising role of gemini surfactant-incorporated liposomes as an efficient drug delivery platform for cardiovascular therapeutics, offering insights into the molecular interactions and structural dynamics underlying their enhanced performance.

Assessing the interaction between the N-terminal region of the membrane protein magnesium transporter A and a lipid bilayer

This study investigates the interaction of KEIF, the intrinsically disordered N-terminal region of the magnesium transporter MgtA, with lipid bilayers mimicking cell membranes. Combining experimental techniques such as neutron reflectometry (NR), quartz-crystal microbalance with dissipation monitoring (QCM-D), synchrotron radiation circular dichroism (SRCD), and oriented circular dichroism (OCD), with molecular dynamics (MD) simulations, we characterized KEIF's adsorption behavior.

HYPOTHESIS

KEIF undergoes conformational changes upon interacting with lipid bilayers, potentially influencing MgtA's function within the plasma membrane.

EXPERIMENTS

The study assessed KEIF's structural transitions and position within lipid bilayers under various conditions, including zwitterionic versus anionic bilayers and different salt concentrations. The techniques analyzed adsorption-induced structural shifts and peptide localization within the bilayer.

FINDINGS

KEIF transitions from a disordered to a more structured state, notably increasing α-helical content upon adsorption to lipid bilayers. The peptide resides primarily in the hydrophobic tail region of the bilayer, where it may displace lipids. Electrostatic interactions, modulated by bilayer charge and ionic strength, play a critical role. These results suggest that KEIF's conformational changes and bilayer interactions can be integral to its potential modulatory role in MgtA function within the plasma membrane. This research highlights the importance of surface-induced structural transitions in intrinsically disordered proteins and their implications for membrane protein modulation.

Charge, Hydrophobicity, and Lipid Type Drive Antimicrobial Peptides' Unique Perturbation Ensembles

Antimicrobial peptides (AMPs) have emerged as a promising solution to the escalating public health threat caused by multidrug-resistant bacteria. Although ongoing research efforts have established AMP's role in membrane permeabilization and leakage, the precise mechanisms driving these disruption patterns remain unclear. We leverage molecular dynamics (MD) simulations enhanced by membrane mimetic (HMMM) to systematically investigate how the physiochemical properties of magainin (+3) and pexiganan (+9) affect their localization, insertion, curvature perturbation, and membrane binding ensemble. Building on existing microbiology, NMR, circular dichroism, and fluorescence data, our analysis reveals that the lipid makeup is a key determinant in the binding dynamics and structural conformation of AMPs. We find that phospholipid type is crucial for peptide localization, demonstrated through magainin's predominant interaction with lipid tails and pexiganan's with polar headgroups in POPC/POPS membranes. The membrane curvature changes induced by pexiganan relative to magainin suggest that AMPs with larger charges have more potential in modulating bilayer bending. These insights advance our understanding of AMP-membrane interactions at the molecular level, offering guidance for the design of targeted antimicrobial therapies.

Atomistic Insights into gp82 Binding: A Microsecond, Million-Atom Exploration of Trypanosoma cruzi Host-Cell Invasion

Chagas disease, caused by the protozoan Trypanosoma cruzi, affects millions globally, leading to severe cardiac and gastrointestinal complications in its chronic phase. The invasion of host cells by T. cruzi is mediated by the interaction between the parasite's glycoprotein gp82 and the human receptor lysosome-associated membrane protein 2 (LAMP2). While experimental studies have identified a few residues involved in this interaction, a comprehensive molecular-level understanding has been lacking. In this study, we present a 1.44-million-atom computational model of the gp82 complex, including over 3300 lipids, glycosylation sites, and full molecular representations of gp82 and LAMP2, making it the most complete model of a parasite-host interaction to date. Using microsecond-long molecular dynamics simulations and dynamic network analysis, we identified critical residue interactions, including novel regions of contact that were previously uncharacterized. Our findings also highlight the significance of the transmembrane domain of LAMP2 in stabilizing the complex. These insights extend beyond traditional hydrogen bond interactions, revealing a complex network of cooperative motions that facilitate T. cruzi invasion. This study not only confirms key experimental observations but also uncovers new molecular targets for therapeutic intervention, offering a potential pathway to disrupt T. cruzi infection and combat Chagas disease.

Cell membranes sustain phospholipid imbalance via cholesterol asymmetry

Membranes are molecular interfaces that compartmentalize cells to control the flow of nutrients and information. These functions are facilitated by diverse collections of lipids, nearly all of which are distributed asymmetrically between the two bilayer leaflets. Most models of biomembrane structure and function include the implicit assumption that these leaflets have similar abundances of phospholipids. Here, we show that this assumption is generally invalid and investigate the consequences of lipid abundance imbalances in mammalian plasma membranes (PMs). Using lipidomics, we report that cytoplasmic leaflets of human erythrocyte membranes have >50% overabundance of phospholipids compared with exoplasmic leaflets. This imbalance is enabled by an asymmetric interleaflet distribution of cholesterol, which regulates cellular cholesterol homeostasis. These features produce unique functional characteristics, including low PM permeability and resting tension in the cytoplasmic leaflet that regulates protein localization.

Long-term editing of brain circuits in mice using an engineered electrical synapse

Electrical signaling across distinct populations of brain cells underpins cognitive and emotional function; however, approaches that selectively regulate electrical signaling between two cellular components of a mammalian neural circuit remain sparse. Here, we engineered an electrical synapse composed of two connexin proteins found in Morone americana (white perch fish) - connexin34.7 and connexin35 - to accomplish mammalian circuit modulation. By exploiting protein mutagenesis, devising a new in vitro system for assaying connexin hemichannel docking, and performing computational modeling of hemichannel interactions, we uncovered a structural motif that contributes to electrical synapse formation. Targeting these motifs, we designed connexin34.7 and connexin35 hemichannels that dock with each other to form an electrical synapse, but not with other major connexins expressed in the mammalian central nervous system. We validated this electrical synapse in vivo using C. elegans and mice, demonstrating that it can strengthen communication across neural circuits composed of pairs of distinct cell types and modify behavior accordingly. Thus, we establish 'Long-term integration of Circuits using connexins' (LinCx) for precision circuit-editing in mammals.

Conformational dynamics and activation of membrane-associated human Group IVA cytosolic phospholipase A 2 (cPLA 2 )

Cytosolic phospholipase A 2 (cPLA 2 ) associates with membranes where it hydrolyzes phospholipids containing arachidonic acid to initiate an inflammatory cascade. All-atom molecular dynamics simulations were employed to understand the activation process when cPLA 2 associates with the endoplasmic reticulum (ER) membrane of macrophages where it acts. We found that membrane association causes the lid region of cPLA 2 to undergo a closed-to-open state transition that is accompanied by the sideways movement of loop 495-540, allowing the exposure of a cluster of lysine residues (K488, K541, K543, and K544), which binds the allosteric activator PIP 2 in the membrane. The active site of the open form of cPLA 2 , containing the catalytic dyad residues S228 and D549, exhibited a three-fold larger cavity than the closed form of cPLA 2 in aqueous solution. These findings provide mechanistic insight as to how cPLA 2 ER membrane association promotes major transitions between conformational states critical to allosteric activation and enzymatic phospholipid hydrolysis.

Computational analysis of the structural-functional dynamics of a Co-receptor proteoglycan

Nerve-Glial Antigen 2/Chondroitin Sulphate Proteoglycan 4 (NG2/CSPG4) is the largest membrane-intercalated cell surface component of the human proteome known to date. NG2/CSPG4 is endowed with the capability of engaging a myriad of molecular interactions and exert co-receptor functions, of which primary ones are sequestering of growth factors and the anchoring of cells to the extracellular matrix. However, the nature of the interactive dynamics of the proteoglycan remains veiled because of its conspicuous size and structural complexity. By leveraging on a multi-scale in silico approach, we have pioneered a comprehensive computational analysis of the structural-functional traits of the NG2/CSPG4 ectodomain. The modelling highlights an intricate assembly of β-sheet motifs linked together by flexible loops. Furthermore, our in silico predictions highlight that the previously delineated D1 domain may consistently remain more accessible for molecular interplays with respect to the D2 and D3 domains. Based on these findings, we have simulated the structural mechanism through the proteoglycan may serve as a co-receptor for growth factor FGF-2, showing that NG2/CSPG4 bends towards the receptor FGFR-1 for this growth factor and confirming the previously hypothesized trimeric complex formation promoted by FGF-2 dimers bridging the FGFR-1-proteoglycan interaction. The Chondroitin Sulphate Proteoglycan 4 is a large multi-domain transmembrane protein involved in several biological processes including pathological conditions. Despite its importance, it has never been studied at the atomistic level due to its large size. Here, we employed a multi-scale computer simulations approach to study its three-dimensional structure, its movements and co-receptor properties, showing that it can serve as mediator in the growth factor signaling process.

European Robin Cryptochrome-4a Associates with Lipid Bilayers in an Ordered Manner, Fulfilling a Molecular-Level Condition for Magnetoreception

Since the middle of the 20th century, long-distance avian migration has been known to rely partly on geomagnetic field. However, the underlying sensory mechanism is still not fully understood. Cryptochrome-4a (ErCry4a), found in European robin (Erithacus rubecula), a night-migratory songbird, has been suggested to be a magnetic sensory molecule. It is sensitive to external magnetic fields via the so-called radical-pair mechanism. ErCry4a is primarily located in the outer segments of the double-cone photoreceptor cells in the eye, which contain stacked and highly ordered membranes that could facilitate the anisotropic attachment of ErCry4a needed for magnetic compass sensing. Here, we investigate possible interactions of ErCry4a with a model membrane that mimics the lipid composition of outer segments of vertebrate photoreceptor cells using experimental and computational approaches. Experimental results show that the attachment of ErCry4a to the membrane could be controlled by the physical state of lipid molecules (average area per lipid) in the outer leaflet of the lipid bilayer. Furthermore, polarization modulation infrared reflection absorption spectroscopy allowed us to determine the conformation, motional freedom, and average orientation of the α-helices in ErCry4a in a membrane-associated state. Atomistic molecular dynamics studies supported the experimental results. A ∼ 1000 kcal mol-1 decrease in the interaction energy as a result of ErCry4a membrane binding was determined compared to cases where no protein binding to the membrane occurred. At the molecular level, the binding seems to involve negatively charged carboxylate groups of the phosphoserine lipids and the C-terminal residues of ErCry4a. Our study reveals a potential direct interaction of ErCry4a with the lipid membrane and discusses how this binding could be an essential step for ErCry4a to propagate a magnetic signal further and thus fulfill a role as a magnetoreceptor.

Narrowed pore conformations of aquaglyceroporins AQP3 and GlpF

Aquaglyceroporins such as aquaporin-3 (AQP3) and its bacterial homologue GlpF facilitate water and glycerol permeation across lipid bilayers. X-ray crystal structures of GlpF showed open pore conformations, and AQP3 has also been predicted to adopt this conformation. Here we present cryo-electron microscopy structures of rat AQP3 and GlpF in different narrowed pore conformations. In n-dodecyl-β-D-maltopyranoside detergent micelles, aromatic/arginine constriction filter residues of AQP3 containing Tyr212 form a 2.8-Å diameter pore, whereas in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) nanodiscs, Tyr212 inserts into the pore. Molecular dynamics simulation shows the Tyr212-in conformation is stable and largely suppresses water permeability. AQP3 reconstituted in POPC liposomes exhibits water and glycerol permeability, suggesting that the Tyr212-in conformation may be altered during permeation. AQP3 Y212F and Y212T mutant structures suggest that the aromatic residue drives the pore-inserted conformation. The aromatic residue is conserved in AQP7 and GlpF, but neither structure exhibits the AQP3-like conformation in POPC nanodiscs. Unexpectedly, the GlpF pore is covered by an intracellular loop, but the loop is flexible and not primarily related to the GlpF permeability. Our findings illuminate the unique AQP3 conformation and structural diversity of aquaglyceroporins.

Probing the Dynamics of Yersinia Adhesin A (YadA) in Outer Membranes Hints at Requirements for β-Barrel Membrane Insertion

The vast majority of cells are protected and functionalized by a dense surface layer of glycans, proteoglycans, and glycolipids. This surface represents an underexplored space in structural biology that is exceedingly challenging to recreate in vitro. Here, we investigate β-barrel protein dynamics within an asymmetric outer membrane environment, with the trimeric autotransporter Yersinia adhesin A (YadA) as an example. Magic-angle spinning NMR relaxation data and a model-free approach reveal increased mobility in the second half of strand β2 after the conserved G72, which is responsible for membrane insertion and autotransport, and in the subsequent loop toward β3. In contrast, the protomer-protomer interaction sites (β1i-β4i-1) are rigid. Intriguingly, the mobility in the β-strand section following G72 is substantially elevated in the outer membrane and less so in the detergent environment of microcrystals. A possible source is revealed by molecular dynamics simulations that show the formation of a salt bridge involving E79 and R76 in competition with a dynamic interplay of calcium binding by E79 and the phosphate groups of the lipids. An estimation of overall barrel motion in the outer membrane and detergent-containing crystals yields values of around 41 ns for both. The global motion of YadA in the outer membrane has a stronger rotational component orthogonal to the symmetry axis of the trimeric porin than in the detergent-containing crystal. In summary, our investigation shows that the mobility in the second half of β2 and the loop to β3 required for membrane insertion and autotransport is maintained in the final folded form of YadA.

Water occupancy in the Acinetobacter baumannii F-ATP synthase c-ring and its implications as a novel inhibitor target

The Acinetobacter baumannii F1FO-ATP synthase is essential for the opportunistic human pathogen. Its membrane-embedded FO domain consists of the c-ring and subunit a. The c-ring translocates protons via a conserved carboxylate across the membrane via two half-channels in subunit a, and its revolution enables the F1 domain to carry out ATP formation. Here, we used molecular dynamics simulations, free energy calculations, and in vivo mutational experiments to assess the likely existence of water molecules in the binding site of the A. baumannii c-ring. We first predicted its binding site structure in the ion-locked conformation and extrapolated the presence of two water molecules in the ion-binding site. Based on our predictions, amino acid point mutations confirmed the critical role of key residues involved in the water-binding site upon ATP synthesis ability and cell growth. We discuss the implications of our findings in the context of rational drug design to target the A. baumannii FO domain.

Computational Methods for Modeling Lipid-Mediated Active Pharmaceutical Ingredient Delivery

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

Virtual screening of potential inhibitors of the ATPase site in Acinetobacter baumannii DNA Gyrase

Bacterial resistance is a global public health problem because of the ineffectiveness of conventional antibiotics against super pathogens. To counter this situation, the search for or design of new molecules is essential to inhibit the key proteins involved in several stages of bacterial infection. One of these key proteins is DNA gyrase, which is responsible for packaging and unfolding of DNA chains during replication. Virtual screening calculations of 583,900 molecules against the ATPase site of DNA gyrase (PDB ID 7PQM) resulted in three promising molecules (Z927783420, Z4422201766, and Z2440107042) with significant binding modes at the active site, according to molecular docking studies. Additionally, they exhibited lower toxicological profiles than the previously reported 80S inhibitors. Molecular dynamics calculations revealed crucial interactions responsible for the inhibition process, with residues ASP87, GLU94, and ASN60 belonging to the ATPase site. On the other hand, the binding energy calculated using the MM/GBSA protocol highlighted Z2440107042 as the most promising inhibitor, with the best binding energy (-74.77 kcal/mol), suggesting that this molecule is a strong candidate for further biological studies.

Molecular mechanism of type ib MET inhibitors and their potential for CNS tumors

The emergence of targeted therapies for MET exon 14 (METex14) skipping mutations has significantly changed the treatment landscape for NSCLC and other solid tumors. The skipping of METex14 results in activating the MET-HGF pathway and promoting tumor cell proliferation, migration, and preventing apoptosis. Type Ib MET inhibitors, designed to selectively target the "DFG-in" conformation of MET, characteristically bind to the ATP-binding pocket of MET in a U-shaped conformation, extending into the solvent-accessible region and interact strongly with residue Y1230 through π-π interactions, have shown remarkable efficacy in treating METex14-altered NSCLC, including cases with brain metastases (BMs). Notably, vebreltinib and capmatinib have demonstrated superior blood-brain barrier (BBB) permeability in both computational and experimental models, highlighting their potential for treating the central nervous system (CNS) metastases. P-glycoprotein (P-gp) is highly expressed in the BBB, which limits the brain uptake of many highly lipophilic drugs. Despite challenges posed by P-gp mediated efflux, vebreltinib has emerged as a promising candidate for CNS treatment due to its favorable pharmacokinetic profile and minimal susceptibility to P-gp efflux. This study underscores the importance of molecular dynamics simulations in predicting drug efficacy and BBB penetration, providing valuable insights for the development of CNS-targeted metastases therapies.

Hidden complexity of α7 nicotinic acetylcholine receptor desensitization revealed by MD simulations and Markov state modeling

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that plays an important role in neuronal signaling throughout the nervous system. Its implication in neurological disorders and inflammation has spurred the development of numerous compounds that enhance channel activation. However, the therapeutic potential of these compounds has been limited by the characteristically fast desensitization of the α7 receptor. Using recent high-resolution structures from cryo-EM, and all-atom molecular dynamic simulations augmented by Markov state modeling, here we explore the mechanism of α7 receptor desensitization and its implication on allosteric modulation. The results provide a precise characterization of the desensitization gate and illuminate the mechanism of ion-pore opening/closing with an agonist bound. In addition, the simulations reveal the existence of a short-lived, open-channel intermediate between the activated and desensitized states that rationalizes the paradoxical pharmacology of the L247T mutant and may be relevant to type-II allosteric modulation. This analysis provides an interpretation of the signal transduction mechanism and its regulation in α7 receptors.

The Influence of Phosphoinositide Lipids in the Molecular Biology of Membrane Proteins: Recent Insights from Simulations

The phosphoinositide family of membrane lipids play diverse and critical roles in eukaryotic molecular biology. Much of this biological activity derives from interactions of phosphoinositide lipids with integral and peripheral membrane proteins, leading to modulation of protein structure, function, and cellular distribution. Since the discovery of phosphoinositides in the 1940s, combined molecular biology, biophysical, and structural approaches have made enormous progress in untangling this vast and diverse cellular network of interactions. More recently, in silico approaches such as molecular dynamics simulations have proven to be an asset in prospectively identifying, characterising, explaining the structural basis of these interactions, and in the best cases providing atomic level testable hypotheses on how such interactions control the function of a given membrane protein. This review details a number of recent seminal discoveries in phosphoinositide biology, enabled by advanced biomolecular simulation, and its integration with molecular biology, biophysical, and structural biology approaches. The results of the simulation studies agree well with experimental work, and in a number of notable cases have arrived at the key conclusion several years in advance of the experimental structures. SUMMARY: Hedger and Yen review developments in simulations of phosphoinositides and membrane proteins.

Mechanistic Insights into the Divergent Membrane Activities of a Viroporin from Chikungunya Virus and Its Transframe Variant

Alphaviruses, a genus of vector-borne viruses in the Togaviridae family, encode a small ion-channel-forming protein, 6K, and its transframe variant (TF) during infections. Although 6K/TF have vital roles in glycoprotein transport, virus assembly, and budding, there is no mechanistic explanation for these functions. We investigated the distinct biochemical functionalities of 6K and TF from the mosquito-borne alphavirus, Chikungunya Virus. We show that like 6K, TF is also capable of forming ion channels in bilayer membranes. The assemblies formed by 6K in membranes are structurally more complex and potentially more ion-restrictive than those formed by TF. Both 6K and TF show strong affinity toward the ER membranes, indicating that the localization of these components at the plasma membrane, as previously reported, is either linked to post-translational modification or mediated through interaction with binding partners. These structural and functional insights may elucidate the distinct roles of 6K and TF in the alphavirus life cycle.

Computational Insights into Membrane Disruption by Cell-Penetrating Peptides

Cell-penetrating peptides (CPPs) can translocate into cells without inducing cytotoxicity. The internalization process implies several steps at different time scales ranging from microseconds to minutes. We combine adaptive Steered Molecular Dynamics (aSMD) with conventional Molecular Dynamics (cMD) to observe nonequilibrium and equilibrium states to study the early mechanisms of peptide-bilayer interaction leading to CPPs internalization. We define three membrane compositions representing bilayer sections, neutral lipids (i.e., upper leaflet), neutral lipids with cholesterol (i.e., hydrophobic core), and neutral/negatively charged lipids with cholesterol (i.e., lower leaflet) to study the energy barriers and disruption mechanisms of Arg9, MAP, and TP2, representing cationic, amphiphilic, and hydrophobic CPPs, respectively. Cholesterol and negatively charged lipids increase the energetic barriers for the peptide-bilayer crossing. TP2 interacts with the bilayer by hydrophobic insertion, while Arg9 disrupts the bilayer by forming transient or stable pores. MAP has shown both behaviors. Collectively, these findings underscore the significance of innovative computational approaches in studying membrane-disruptive peptides and, more specifically, in harnessing their potential for cell penetration.

Role of Divalent Ions in Membrane Models of Polymyxin-Sensitive and Resistant Gram-Negative Bacteria

Polymyxins, critical last-resort antibiotics, impact the distribution of membrane-bound divalent cations in the outer membrane of Gram-negative bacteria. We employed atomistic molecular dynamics simulations to model the effect of displacing these ions. Two polymyxin-sensitive and two polymyxin-resistant models of the outer membrane of Salmonella enterica were investigated. First, we found that the removal of all calcium ions induces global stress on the model membranes, leading to substantial membrane restructuring. Next, we used enhanced sampling methods to explore the effects of localized stress by displacing membrane-bound ions. Our findings indicate that creating defects in the membrane-bound ion network facilitates polymyxin permeation. Additionally, our study of polymyxin-resistant mutations revealed that divalent ions in resistant model membranes are less likely to be displaced, potentially contributing to the increased resistance associated with these mutations. Lastly, we compared results from all-atom molecular dynamics simulations with coarse-grained simulations, demonstrating that the choice of force field significantly influences the behavior of membrane-bound ions under stress.

In Slico Screening and In Vitro Identification of Hyperuricemia-Inhibiting Peptides from Trachurus japonicus

Hyperuricemia arises from imbalanced uric acid metabolism, contributing to gout and related chronic diseases. When traditional drugs are used to treat hyperuricemia, side effects are inevitable, which promotes the exploration of new bioactive compounds. Protein hydrolysates and peptides are gradually showing potential in the treatment of hyperuricemia. This study investigated the uric acid inhibitory activity of peptides extracted from Trachurus japonicus using in silico and in vitro methods. We employed in silico virtual enzymolysis and experimental validation to identify bioactive peptides from Trachurus japonicus proteins. Four peptides (DF, AGF, QPSF, and AGDDAPR) were comprehensively screened by molecular docking and database analysis. After solid-phase synthesis, the inhibitory effects of these peptides on hyperuricemia were further verified in vitro and at the cellular level. The results showed that all four peptides have good hyperuricemia-inhibiting activities. Molecular docking and molecular dynamics revealed that peptides DF and AGDDAPR affect the production of uric acid by binding to the active sites of urate transporter 1 (URAT1), glucose transporter 9 (GLUT9), and xanthine oxidase (XOD), while peptides QPSF and AGF mainly influence the XOD active site, confirming that it is feasible to rapidly screen hyperuricemia-inhibiting peptides by molecular docking.

Investigation of Muscarinic Acetylcholine Receptor M3 Activation in Atomistic Detail: A Chemist's Viewpoint

We analyzed the precise ligand:receptor interactions required for activation of the muscarinic acetylcholine receptor M3, a prototypical G protein-coupled receptor and potential diabetes target. Starting from literature-known compounds and docking solutions, ligands were tailored for the modulation of this receptor's activation. Several aspects of the structure-activity relationship of agonists were investigated in atomistic detail, in order to delineate how the receptor can be activated via the orthosteric site. Such exquisitely precise knowledge is instrumental for designing potent and effiacious ligands. We put this strategy into practice and acquired or synthesized and measured a diverse set of 55 ligands ranging from small fragment-like amines coordinating D3.32 to bigger molecules extending towards helices 5 and 6 with diphenyl moieties. In the course of these investigations, we showed that the polarizability of the amine nitrogen and the rigidity and size of the moieties in the space delimited by helices 5 and 6 are the two key elements distinguishing potent and efficacious ligands from those that are not. The resulting data set will be highly useful in drug design and molecular machine learning alike.

Lipid droplet targeting of the lipase coactivator ABHD5 and the fatty liver disease-causing variant PNPLA3 I148M is required to promote liver steatosis

The storage and release of triacylglycerol (TAG) in lipid droplets (LDs) is regulated by dynamic protein interactions. α/β Hydrolase domain-containing protein 5 (ABHD5; also known as CGI-58) is a membrane/LD-bound protein that functions as a co-activator of patatin-like phospholipase domain-containing 2 (PNPLA2; also known as adipose triglyceride lipase) the rate-limiting enzyme for TAG hydrolysis. The dysregulation of TAG hydrolysis is involved in various metabolic diseases such as metabolic dysfunction-associated steatotic liver disease (MASLD). We previously demonstrated that ABHD5 interacted with PNPLA3, a closely related family member to PNPLA2. Importantly, a common missense variant in PNPLA3 (I148M) is the greatest genetic risk factor for MASLD. PNPLA3 148M functions to sequester ABHD5 and prevent coactivation of PNPLA2, which has implications for initiating MASLD; however, the exact mechanisms involved are not understood. Here, we demonstrate that LD targeting of both ABHD5 and PNPLA3 I148M is required for the interaction. Molecular modeling demonstrates important residues in the C terminus of PNPLA3 for LD binding and fluorescence cross-correlation spectroscopy demonstrates that PNPLA3 I148M has greater association with ABHD5 than WT PNPLA3. Moreover, the C terminus of PNPLA3 is sufficient for functional targeting of PNPLAs to LD and the interaction with ABHD5. In addition, ABHD5 is a general binding partner of LD-bound PNPLAs. Finally, PNPLA3 I148M targeting to LD is required to promote steatosis in vitro and in the liver. Overall results suggest that the interaction of PNPLA3 I148M with ABHD5 on LD is required to promote liver steatosis.

In silico evaluation of 4-thiazolidinone-based inhibitors against the receptor for advanced glycation end products (RAGE)

Non-enzymatic glycation of biomolecules by reducing sugars led to several products, including the advanced glycation end products (AGEs), the accumulation of which has been linked to various life-threatening diseases. The binding of AGEs to their respective protein receptors for advanced glycation end products (RAGE) can initiate a cascade of reactions, which may alter physiological conditions. The present work investigates the potential of 4-thiazolidinones as RAGE inhibitors. We performed an extensive computational study to identify the structural requirements needed to act as RAGE inhibitors. To achieve this goal, 4-thiazolidinone-based compounds available in PubChem, ZINC15, ChEMBL, and ChEBI databases were screened against RAGE (PDB: 4LP5), leading to the identification of top five drug-like candidates with a high binding affinity to RAGE V-domain catalytic region. Drug likeness, absorption, distribution, metabolism, excretion, and toxicity (ADMET) of the top-scoring compounds have been studied and discussed. Global molecular descriptors, chemical reactivity, hardness, softness, etc., have been estimated. Finally, molecular dynamics (MD) simulations at 100 ns were carried out to check the stability and other properties. Overall, we believe that the identified compounds can potentially attenuate RAGE-AGE interactions.Communicated by Ramaswamy H. Sarma.

Fluoxetine Alters the Biophysics of DPPC and DPPG Bilayers through Phase-Dependent and Electrostatic Interactions

Lipid membranes can control the permeability of a pharmaceutical drug, whereas the drug can induce changes in the structural and biophysical properties of the membranes. Understanding this interplay of drug-lipid membrane interactions can be of great importance in drug design. Here, we present a molecular dynamics study to provide insights into the interactions between the antidepressant fluoxetine and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) or 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG) bilayers. It was found that, due to the electrostatic interaction, the headgroup of the zwitterionic DPPC lipid is more stable than that of the negatively charged DPPG lipid, allowing the gel phase to persist even at the elevated temperature. At 25 °C, fluoxetine cannot penetrate into the gel-phase DPPC bilayer, while the electrostatic interaction between positively charged fluoxetine and negatively charged DPPG bilayer retains the drug within the lipid headgroup domain. When the temperature is increased to 45 °C, both neutral and charged forms of fluoxetine can partition into the DPPC and DPPG bilayers spontaneously. Analysis of the biophysical and structural changes in both DPPC and DPPG bilayers in the presence of fluoxetine revealed a phase-dependent effect. The binding of fluoxetine to the lipid bilayers limits the movement and orientation of the drug. These findings shed light on the interactions between a commonly prescribed antidepressant and lipid membranes, and such information can be beneficial to the development of potential therapeutic agents.

OPRLM: A Web Tool and a Database for Positioning and Simulations of Proteins in Realistic Lipid Membranes

Molecular dynamics (MD) simulations in explicit lipid bilayers enable modeling of protein-lipid interactions essential for membrane protein functions and regulation. The newly developed computational web tool, OPRLM (Orientations of Proteins in Realistic Lipid Membranes), automates the assembly of membrane protein structures with explicit lipids corresponding to 18 biological membrane types with symmetric or asymmetric lipid distributions, as well as 5 types of two-component lipid bilayers with varying cholesterol content. Built upon the CHARMM-GUI toolset and the PPM method, OPRLM simplifies the setup of complex simulation system involving integral and/or peripheral membrane proteins with explicit lipid mixtures and generates all necessary files for subsequent all-atom MD simulations. OPRLM has successfully generated protein-membrane systems for 286 tested protein structures in various biomembranes, including 138 structures containing ligands. The OPRLM database, an advanced successor of the OPM database, includes explicit protein-lipid systems for tested proteins in their native biomembranes. It provides coordinates of integral and peripheral membrane proteins from the Protein Data Bank embedded in planar or curved implicit lipid bilayers. Additionally, it includes the classification of proteins into types, superfamilies, and families, along with the information on intracellular localizations and membrane topology and visualization options. The OPRLM web tool and the database are publicly accessible at https://oprlm.org.

Exploring the impact of the stargazin V143L mutation on the dynamics of the AMPA receptor: stargazin complex

Stargazin, a transmembrane AMPAR regulatory protein (TARP), plays a crucial role in facilitating the transport of AMPA receptors to the cell surface, stabilising their localisation at synapses and influencing their gating properties. The primary objective of this study was to investigate the effect of the V143L mutation in stargazin, previously linked to intellectual disability, on the interaction between stargazin and AMPA receptors. To achieve this, we conducted a thorough examination of eight distinct molecular dynamics simulations of AMPA receptor-stargazin complexes, each associated with different conductance levels. Through extensive analysis of complex interface structures and dynamics, we revealed that the stargazin V143L mutation had a more pronounced destabilising effect on complexes with lower conductance levels than on the conductive states of the receptor, suggesting a potential association with impaired synaptic transmission in individuals with this mutation.

In silico drug repurposing at the cytoplasmic surface of human aquaporin 1

Aquaporin 1 (AQP1) is a key channel for water transport in peritoneal dialysis. Inhibition of AQP1 could therefore impair water transport during peritoneal dialysis. It is not known whether inhibition of AQP1 occurs unintentionally due to off-target interactions of administered medications. A high-throughput virtual screening study has been performed to investigate the possible binding of licensed medications to the water pore of human AQP1. A complete model of human AQP1 based on its canonical sequence was assembled using I-TASSER and MODELLER. The model was refined via the incorporation of pore water molecules from a high-resolution yeast aquaporin structure. Docking studies were conducted for the cytoplasmic domain of the AQP1 monomer against a library of all compounds listed in the British National Formulary (BNF), using the PLANTS software with the ChemPLP scoring function. The stability of the best docked conformations within the intrinsic water pore was assessed via short 15 nanosecond molecular dynamics (MD) simulations using the GROMACS-on-Colab utility. Of the 1512 compounds tested, 1002 docking results were obtained, and 198 of these conformations occupied a position within the intrinsic water pore. 30 compounds with promising docking scores were assessed by MD. The docked conformations for dopamine, gabapentin, pregabalin, and methyldopa were stable in these short MD studies. For furosemide and pravastatin, the MD trajectory suggested a binding mode different to the docking result. A small set of compounds which could impede water transport through human AQP1 have been identified in this computational screening study.

Complexes of vertebrate TMC1/2 and CIB2/3 proteins form hair-cell mechanotransduction cation channels

Calcium and integrin-binding protein 2 (CIB2) and CIB3 bind to transmembrane channel-like 1 (TMC1) and TMC2, the pore-forming subunits of the inner-ear mechano-electrical transduction (MET) apparatus. These interactions have been proposed to be functionally relevant across mechanosensory organs and vertebrate species. Here, we show that both CIB2 and CIB3 can form heteromeric complexes with TMC1 and TMC2 and are integral for MET function in mouse cochlea and vestibular end organs as well as in zebrafish inner ear and lateral line. Our AlphaFold 2 models suggest that vertebrate CIB proteins can simultaneously interact with at least two cytoplasmic domains of TMC1 and TMC2 as validated using nuclear magnetic resonance spectroscopy of TMC1 fragments interacting with CIB2 and CIB3. Molecular dynamics simulations of TMC1/2 complexes with CIB2/3 predict that TMCs are structurally stabilized by CIB proteins to form cation channels. Overall, our work demonstrates that intact CIB2/3 and TMC1/2 complexes are integral to hair-cell MET function in vertebrate mechanosensory epithelia.

Isoform-specific N-linked glycosylation of NaV channel α-subunits alters β-subunit binding sites

Voltage-gated sodium channel α-subunits (NaV1.1-1.9) initiate and propagate action potentials in neurons and myocytes. The NaV β-subunits (β1-4) have been shown to modulate α-subunit properties. Homo-oligomerization of β-subunits on neighboring or opposing plasma membranes has been suggested to facilitate cis or trans interactions, respectively. The interactions between several NaV channel isoforms and β-subunits have been determined using cryogenic electron microscopy (cryo-EM). Interestingly, the NaV cryo-EM structures reveal the presence of N-linked glycosylation sites. However, only the first glycan moieties are typically resolved at each site due to the flexibility of mature glycan trees. Thus, existing cryo-EM structures may risk de-emphasizing the structural implications of glycans on the NaV channels. Herein, molecular modeling and all-atom molecular dynamics simulations were applied to investigate the conformational landscape of N-linked glycans on NaV channel surfaces. The simulations revealed that negatively charged sialic acid residues of two glycan sites may interact with voltage-sensing domains. Notably, two NaV1.5 isoform-specific glycans extensively cover the α-subunit region that, in other NaV channel α-subunit isoforms, corresponds to the binding site for the β1- (and likely β3-) subunit immunoglobulin (Ig) domain. NaV1.8 contains a unique N-linked glycosylation site that likely prevents its interaction with the β2 and β4-subunit Ig-domain. These isoform-specific glycans may have evolved to facilitate specific functional interactions, for example, by redirecting β-subunit Ig-domains outward to permit cis or trans supraclustering within specialized cellular compartments such as the cardiomyocyte perinexal space. Further experimental work is necessary to validate these predictions.

Differential Inhibition by Cenobamate of Canonical Human Nav1.5 Ion Channels and Several Point Mutants

Cenobamate is a new and highly effective antiseizure compound used for the treatment of adults with focal onset seizures and particularly for epilepsy resistant to other antiepileptic drugs. It acts on multiple targets, as it is a positive allosteric activator of γ-aminobutyric acid type A (GABAA) receptors and an inhibitor of neuronal sodium channels, particularly of the late or persistent Na+ current. We recently evidenced the inhibitory effects of cenobamate on the peak and late current component of the human cardiac isoform hNav1.5. The determined apparent IC50 values of 87.6 µM (peak) and 46.5 µM (late current) are within a clinically relevant range of concentrations (the maximal plasma therapeutic effective concentration for a daily dose of 400 mg in humans is 170 µM). In this study, we built a 3D model of the canonical hNav1.5 channel (UniProt Q14524-1) in open conformation using AlphaFold2, embedded it in a DPPC lipid bilayer, corrected the residue protonation state (pH 7.2) with H++, and added 2 Na+ ions in the selectivity filter. By molecular docking, we found the cenobamate binding site in the central cavity. We identified 10-point mutant variants in the binding site region and explored them via docking and MD. Mutants N1462K/Y (rs1064795922, rs199473614) and M1765R (rs752476527) (by docking) and N932S (rs2061582195) (by MD) featured higher predicted affinity than wild-type.

Atomistic Simulations and Analysis of Peripheral Membrane Proteins with Model Lipid Bilayers

All-atom molecular dynamics (AAMD) is a computational technique that predicts the movement of particles based on the intermolecular forces acting on the system. It enables the study of biological systems at atomic detail, complements observations from experiments, and can help the selection of experimental targets. Here, we describe the applications of MD simulations to study the interaction between peripheral membrane proteins and lipid bilayers. Specifically, we provide step-by-step instructions to set up MD simulations to study the binding and interaction of ALPS, the amphipathic helix of the lipid transport protein Osh4, and Thanatin, an antimicrobial peptide with model membranes. We describe examples of systems built with fully atomistic lipid tails and those truncated with the highly-mobile-membrane-mimetic method to enhance conformational sampling. We also comment on the importance of lipid diversity, molecular resolution, and best practices for constructing, running, and analyzing protein-lipid simulation systems. In this second edition, we include a brief discussion on alternative approaches and software to construct protein-membrane coordinate systems, as well as analysis tools and practices that have become relevant to examining protein-lipid interactions since the first edition of this chapter.

Role of sequence length and functionalization in interactions of bioconjugated peptides with mitomembranes

Cell-penetrating peptides are efficient tools for intracellular delivery of a variety of cargoes. In this study, we explored the effect of chain length, side chain chemistry, and the locations of conjugated molecules on the interaction between iron-chelating peptides and a mitochondrial-mimicking membrane. We report that a longer chain length enhanced peptide/membrane interactions, and conjugation at the N-terminus lowered the free-energy barrier for peptide translocation across the membrane. Peptides containing Phe side chains and those containing modified Phe (cyclohexane) side chains showed comparable peptide/membrane energetics and translocation energy barriers. Using steered molecular dynamics (SMD) simulations, we further probed the mechanistic details of translocation of each N-terminated peptide across the membrane and compared their metastable states. At a higher steering velocity, the peptide adopted a compact structure due to frequent π-π interactions among conjugated molecules, but at lower steering velocities, each N-terminated peptide adopted an extended structure. This structure allowed cationic residues to maximize their interactions with phosphate headgroups in the mitomembrane. The hydrophobic residues also formed interactions with the lipid acyl tails, facilitating the passage of peptides across the membrane with decreased free energy barriers. Our results highlight the significance of peptide chain length and conjugation in facilitating peptide transport across the membrane.

Structural Dynamics of Neutral Amino Acid Transporter SLC6A19 in Simple and Complex Lipid Bilayers

B0AT1 (SLC6A19) is a major sodium-coupled neutral amino acid transporter that relies on angiotensin converting enzyme 2 (ACE2) or collectrin for membrane trafficking. Despite its significant role in disorders associated with amino acid metabolism, there is a deficit of comprehensive structure-function understanding of B0AT1 in lipid environment. Herein, we have employed molecular dynamics (MD) simulations to explore the architectural characteristics of B0AT1 in two distinct environments: a simplified POPC bilayer and a complex lipid system replicating the native membrane composition. Notably, our B0AT1 analysis in terms of structural stability and regions of maximum flexibility shows consistency in both the systems with enhanced structural features in the case of complex lipid system. Our findings suggest that diacylglycerol phospholipids significantly alter the pore radius, hydrophobic index, and surface charge distribution of B0AT1, thereby affecting the flexibility of transmembrane helices TM7, TM12, and loop connecting TM7-TM8, crucial for ACE2-B0AT1 interaction. Pro41, Ser190, Arg214, Arg240, Ser413, Pro414, Cys463, and Val582 are among the most prominent lipid binding residues that might influence B0AT1 functionality. We also perceive notable lipid mediated deviation in the degree of tilt and loss of helicity in TM1 and TM6 which might affect the substrate binding sites S1 and S2 in B0AT1. Considerably, destabilization in the structure of B0AT1 in lipid environment was evident upon mutation in TM domain, associated with Hartnup disorder through various structure-based protein stability tools. Our two-tiered approach allowed us to validate the use of POPC as a baseline for initial analyses of SLC transporters. Altogether, our all-atoms MD study provides a platform for future investigations into the structure-function mechanism of B0AT1 in realistic lipid mimetic bilayers and offers a framework for developing new therapeutic agents targeting this transporter.

The physicochemical properties of lipopolysaccharide chemotypes regulate activation of the contact pathway of blood coagulation

Lipopolysaccharide (LPS) is the primary pathogenic factor in Gram-negative sepsis. While the presence of LPS in the bloodstream during infection is associated with disseminated intravascular coagulation, the mechanistic link between LPS and blood coagulation activation remains ill-defined. The contact pathway of coagulation-a series of biochemical reactions that initiates blood clotting when plasma factors XII (FXII) and XI (FXI), prekallikrein (PK), and high molecular weight kininogen interact with anionic surfaces-has been shown to be activated in Gram-negative septic patients. In this study, using an in vivo baboon model of Gram-negative Escherichia coli sepsis, we observed activation of the contact pathway including FXII, FXI, and PK. We examined whether E.coli LPS molecules could bind and activate contact pathway members by quantifying the interaction and activation of either FXII, FXI, or PK with each of the three chemotypes of LPS: O111:B4, O26:B6, or Rd2. The LPS chemotypes exhibited distinct physicochemical properties as aggregates and formed complexes with FXII, FXI, and PK. The LPS chemotype O26:B6 uniquely promoted the autoactivation of FXII to FXIIa and, in complex with FXIIa, promoted the cleavage of FXI and PK to generate FXIa and plasma kallikrein, respectively. Furthermore, in complex with the active forms of FXI or PK, LPS chemotypes were able to regulate the catalytic activity of FXIa and plasma kallikrein, respectively, despite the inability to promote the autoactivation of either zymogen. These data suggest that the procoagulant phenotype of E.coli is influenced by bacterial strain and the physicochemical properties of the LPS chemotypes.

Insertion and Anchoring of the HIV-1 Fusion Peptide into a Complex Membrane Mimicking the Human T-Cell

A fundamental understanding of how the HIV-1 envelope (Env) protein facilitates fusion is still lacking. The HIV-1 fusion peptide, consisting of 15 to 22 residues, is the N-terminus of the gp41 subunit of the Env protein. Further, this peptide, a promising vaccine candidate, initiates viral entry into target cells by inserting and anchoring into human immune cells. The influence of membrane lipid reorganization and the conformational changes of the fusion peptide during the membrane insertion and anchoring processes, which can significantly affect HIV-1 cell entry, remains largely unexplored due to the limitations of experimental measurements. In this work, we investigate the insertion of the fusion peptide into an immune cell membrane mimic through multiscale molecular dynamics simulations. We mimic the native T-cell by constructing a nine-lipid asymmetric membrane, along with geometrical restraints accounting for insertion in the context of gp41. To account for the slow time scale of lipid mixing while enabling conformational changes, we implement a protocol to go back and forth between atomistic and coarse-grained simulations. Our study provides a molecular understanding of the interactions between the HIV-1 fusion peptide and the T-cell membrane, highlighting the importance of the conformational flexibility of fusion peptides and local lipid reorganization in stabilizing the anchoring of gp41 into the targeted host membrane during the early events of HIV-1 cell entry. Importantly, we identify a motif within the fusion peptide critical for fusion that can be further manipulated in future immunological studies.

Plasmalogens Reversed Oxidative Stress and Inflammatory Response Exacerbated by Damage to Cell Membrane Properties in Acute Liver Injury

BACKGROUND

In acute liver injury (ALI), cell membrane damage could induce an inflammatory response and oxidative stress. As a membrane glycerophospholipid, plasmalogens (PLS) are crucial in regulating the cell membrane properties and exhibit beneficial effects in various liver diseases. However, the specific regulatory effects of PLS in the ALI remain unknown.

METHODS

We utilized CCl4 to induce ALI in AML12 hepatocytes and C57BL/6J mice and examined oxidative stress indicators and inflammatory cytokine levels. Our study further validated the effect of PLS on cell membrane integrity by Lactate Dehydrogenase (LDH) release assay and Dil/Calcein assay, and molecular dynamics (MD) simulations were employed to elucidate the molecular mechanisms by which PLS affected cell membranes.

RESULTS

PLS attenuated hepatocyte damage both in vivo and in vitro. Moreover, PLS increased levels of SOD, GSH, and CAT and inhibited the production of malondialdehyde. PLS succeeded in decreasing proinflammatory cytokines (TNF-α, IL-1β, and IL-6) while increasing anti-inflammatory cytokines (IL-10). Furthermore, PLS effectively maintained the cell membrane integrity. The MD simulations well explained the molecular mechanisms: a high level of PLS modulated the cell membrane properties, enabling them to be more flexible, elastic, and less prone to rupture.

CONCLUSIONS

Our study illustrated the effect and molecular mechanisms of PLS against ALI, potentially broadening its application in liver diseases.

In Silico Modeling and Characterization of Epstein-Barr Virus Latent Membrane Protein 1 Protein

Latent membrane protein 1 (LMP1) plays a crucial role in Epstein-Barr virus (EBV)'s ability to establish latency and is involved in developing and progressing EBV-associated cancers. Additionally, EBV-infected cells affect the immune responses, making it challenging for the immune system to eliminate them. Due to the aforementioned reasons, it is crucial to understand the structural features of LMP1, which are essential for the development of novel cancer therapies that target its signaling pathways. To date, there is yet to be a complete LMP1 protein structure; therefore, in our work, we modeled the full-length LMP1 containing the short cytoplasmic N-terminus, six transmembrane domains (TMDs), and a long-simulated C-terminus. Our model showed good stability and protein compactness evaluated through accelerated-molecular dynamics, where the conformational ensemble exhibited compact folds, particularly in the TMDs. Our results suggest that specific domains or motifs, predominantly in the C-terminal domain of LMP1, show promise as potential drug targets. As a whole, our work provides insights into key structural features of LMP1 that will allow the development of novel LMP1 therapies.

Crosstalk of Nucleic Acid Mimics with Lipid Membranes: A Multifaceted Computational and Experimental Study

Nucleic acid mimics (NAMs) have demonstrated high potential as antibacterial drugs. However, very few studies have assessed their possible diffusion across the bacterial envelope. In this work, we studied NAMs' diffusion in lipid bilayer systems that mimic the bacterial outer membrane using molecular dynamics (MD) simulations. Additionally, we examined the interactions of a NAM sequence with lipid membranes and ascertained the partition constants (Kp) through MD and spectroscopic investigations. The NAM sequences were composed of locked nucleic acid (LNA) and 2'-O-methyl (2'-OMe) residues, whereas the membrane models were composed of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) phospholipids. The parametrization protocol followed was validated against literature data and demonstrated the reliability of our approach for simulating NAM sequences. Investigation into the interaction of the sequences with zwitterionic and anionic membranes revealed a preference for hydrogen bond formation with the anionic model over the zwitterionic one. Additionally, potential of mean force (PMF) calculations unveiled a lower free energy barrier for translocation across the zwitterionic bilayer model. Contrarily, the partition constants derived suggested a slightly higher partitioning within the anionic membrane, emphasizing a nuanced interplay of factors. Finally, spectroscopic partition measurements with liposomes presented challenges in quantifying the partition of NAMs due to minimal signal variations. However, a tendency for quenching in anionic vesicles suggested a potential, albeit small, partitioning effect that warrants further investigation. In summary, our study revealed that NAMs will not, in principle, be able to cross an intact bacterial outer membrane by passive diffusion.

Dynamic interplay between a TonB-dependent heme transporter and a TonB protein in a membrane environment

The envelope of Gram-negative bacteria is composed of two membranes separated by the periplasmic space. This organization imposes geometrical and distance constraints that are key for the mechanism of action of multicomponent systems spanning the envelope. However, consideration of all three compartments by experimental approaches is still elusive. Here, we have used the state-of-the-art molecular dynamics simulation in an Escherichia coli envelope model to obtain a dynamic view of molecular interactions between the outer membrane heme transporter HasR and the inner membrane TonB-like protein HasB. Their interaction allows the transfer of the inner membrane proton-motive force derived energy to the transporter for heme internalization. The simulations that incorporate both membranes show the key role of periplasmic domains of both proteins and their dynamics in complex formation and stability. They revealed a previously unidentified network of HasR-HasB protein-protein interactions in the periplasm. Experimental validation (mutations, in vivo phenotypic and biophysical assays) provides support for the simulation-predicted interactions. Based on structural and sequence conservation, the network of interaction revealed in this study is expected to occur in other nutrient import systems.

IMPORTANCE

Gram-negative bacteria import scarce nutrients such as metals and vitamins by an energized mechanism involving a multicomponent protein system that spans the cell envelope. It consists of an outer membrane TonB-dependent transporter (TBDT) and a TonB complex in the inner membrane that provides the proton motive force energy for the nutrient entry. Despite the intense research efforts focused on this system (a) from structural and fundamental microbiology perspectives and (b) for the interest in the development of new antibacterial strategies, the molecular mechanism of the system is not at all well understood. The lack of understanding comes from incomplete structural data and the experimental difficulties of studying an inherently flexible multicomponent complex that resides within the heterogeneous environment of the double membrane bacterial cell envelope. To address these challenges and obtain a comprehensive view of the molecular interactions at atomic level, here, we have used the combined power of advanced molecular simulations and complementary microbiology and biochemical experiments. Our results represent a significant step forward in understanding the structural and molecular bases of this vital mechanism.

A High-Throughput Method for Screening Peptide Activators of G-Protein-Coupled Receptors

Here, we describe an innovative and efficient method for screening peptide activators of G-protein-coupled receptors (GPCRs) utilizing a protein-protein interaction (PPI) approach. We designed a library of 92,918 peptides fused with transmembrane domains of glycosylphosphatidylinositol-anchored proteins (GPI-APs). We employed a pooled lentiviral system to promote the expression of these proteins at the cellular membrane and evaluate their ability to activate GPCRs. We then used fluorescence-activated cell sorting (FACS) to screen the GPI-AP-peptide library and identify novel peptide activators of the glucagon-like peptide-1 receptor (GLP-1R). We discovered one peptide PepA3 derived from the Frizzled-like (FZ) domain of human Carboxypeptidase Z (CPZ), a regulated secreted metallocarboxypeptidase. Notably, PepA3 and its two related variants, PepA and PepA2, activated the GLP-1R receptor with less potency but comparable efficacy to that of GLP-1. We then hypothesized that all of these peptides will bind differently to the GLP-1R than the normal ligand. Our technology could identify novel GPCR-activating peptides for structure-function or drug discovery research.

Cell membranes sustain phospholipid imbalance via cholesterol asymmetry

Membranes are molecular interfaces that compartmentalize cells to control the flow of nutrients and information. These functions are facilitated by diverse collections of lipids, nearly all of which are distributed asymmetrically between the two bilayer leaflets. Most models of biomembrane structure and function often include the implicit assumption that these leaflets have similar abundances of phospholipids. Here, we show that this assumption is generally invalid and investigate the consequences of lipid abundance imbalances in mammalian plasma membranes (PM). Using quantitative lipidomics, we discovered that cytoplasmic leaflets of human erythrocyte membranes have >50% overabundance of phospholipids compared to exoplasmic leaflets. This imbalance is enabled by an asymmetric interleaflet distribution of cholesterol, which regulates cellular cholesterol homeostasis. These features produce unique functional characteristics, including low PM permeability and resting tension in the cytoplasmic leaflet that regulates protein localization. These largely overlooked aspects of membrane asymmetry represent an evolution of classic paradigms of biomembrane structure and physiology.

An explainable few-shot learning model for the directed evolution of antimicrobial peptides

Due to the persistent threat of antibiotic resistance posed by Gram-negative pathogens, the discovery of new antimicrobial agents is of critical importance. In this study, we employed deep learning-guided directed evolution to explore the chemical space of antimicrobial peptides (AMPs), which present promising alternatives to traditional small-molecule antibiotics. Utilizing a fine-tuned protein language model tailored for small dataset learning, we achieved structural modifications of the lipopolysaccharide-binding domain (LBD) derived from Marsupenaeus japonicus, a prawn species of considerable value in aquaculture and commercial fisheries. The engineered LBDs demonstrated exceptional activity against a range of Gram-negative pathogens. Drawing inspiration from evolutionary principles, we elucidated the bactericidal mechanism through molecular dynamics simulations and mapped the directed evolution pathways using a ladderpath framework. This work highlights the efficacy of explainable few-shot learning in the rational design of AMPs through directed evolution.

Unraveling the Membrane Topology of TMEM151A: A Step Towards Understanding its Cellular Role

Transmembrane protein 151A (TMEM151A) has been identified as a causative gene for paroxysmal kinesigenic dyskinesia, though its molecular function remains almost completely unknown. Understanding the membrane topology of transmembrane proteins is crucial for elucidating their functions and possible interacting partners. In this study, we utilized molecular dynamics simulations, immunocytochemistry, and electron microscopy to define the topology of TMEM151A. Our results validate a starting AlphaFold model of TMEM151A and reveal that it comprises a transmembrane domain with two membrane-spanning alpha helices connected by a short extracellular loop and an intramembrane helix-hinge-helix structure. Notably, most of the protein is oriented towards the intracellular side of the membranes with a large cytosolic domain featuring a combination of alpha-helix and beta-sheet structures, as well as the protein N- and C-termini. These insights into TMEM151A's topology and orientation of its domains with respect of the cell membranes provide essential information for future functional studies and represent a first fundamental step for understanding its role in the pathogenesis of paroxysmal kinesigenic dyskinesia.