Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

From Membrane Composition to Antimicrobial Strategies: Experimental and Computational Approaches to AMP Design and Selectivity

The United Nations have committed to end the epidemics of communicable diseases by 2030 (SDG Target 3.3). In contrast with this ambition, the rise of Multi Drug Resistant (MDR) and Pan Drug Resistant (PDR) bacteria poses a threat of a return to the pre-antibiotic era. It is of high priority to find new therapies that target the ESKAPEE group of pathogens and their drug-resistant strains. Antimicrobial peptides (AMPs) are an emerging class of antibiotics that hold promises of overcoming bacterial resistance by using both novel mechanisms of action as well as targeting already known pathways. The chemical space of AMPs is potentially huge and methodologies allowing the rational exploration of novel structures are highly needed. This review focuses on case studies that give novel insights about the mechanisms of action, resistance and selectivity of some relevant AMPs, exemplifying the importance of microscopic, computational and experimental tools. Particular focus will be devoted to bacterial membranes and how AMPs can target them while sparing human plasma membranes, in order to become safer drugs. The lessons learned from the literature cases give directions toward the development of AMPs as drug products.

Design, synthesis and evaluation of lysine- and leucine-rich hydrocarbon-stapled peptides as antibacterial agents

To address the challenge of antimicrobial resistance, we investigated new antibacterial peptides based on lysine- and leucine-rich sequences. We stabilised their membrane-active secondary structures by applying hydrocarbon stapling at sequence positions i and i+4. Stapling improved peptide structural stability in both aqueous and lipid environments, regardless of the staple position. It also enhanced antibacterial efficiency against both gram-negative and gram-positive bacteria, including antibiotic-resistant strains, with minimum inhibitory concentrations (MICs) of 2-4 μM (2.5-5.5 μg/mL). The stapled peptides showed increased resistance to enzymatic degradation, particularly with staples incorporated near the N-terminus, and were not haemolytic or cytotoxic at their MICs. Molecular dynamics simulations revealed how stapling aids in (i) stabilising the membrane-active secondary structure of amphipathic peptides and (ii) accelerating their membrane insertion. Our results provide insight into peptide design for antimicrobial use. We show that hydrocarbon stapling of lysine- and leucine-rich short sequences may offer a pathway towards more stable and effective antibacterial agents.

Identifying and quantifying membrane interactions of the protein human cis-prenyltransferase

Prenyl chains come in multiple sizes, fulfilling unique and essential functions across all domains of life. Prenyl chains are synthesized by prenyltransferase proteins. Despite their structural similarity, prenyltransferases exhibit substantial functional diversity to create lipophilic products of varying lengths. Human cis-prenyltransferase (h-cisPT) is a tetrameric enzyme responsible for the synthesis of long prenyl chains, consisting of 20-prenyl-unit products that are essential to specific posttranslational modifications such as N-glycosylation upon downstream processing. These long products are hypothesized to transfer from h-cisPT to the ER membrane, but the mechanism of this transfer is not known. We use molecular dynamics simulations to identify a consistent membrane binding pose for h-cisPT. By quantifying protein-membrane contacts, we identify the aromatic amino acid residues in the conserved catalytic domain as critical to membrane binding. Determining relative protein-membrane binding free energies through free energy perturbation highlights the importance of these residues for membrane association, as mutations lower membrane affinity by as much as 27 kcal/mol. These results are validated using FRET to demonstrate decreased catalytic activity and membrane binding in response to mutation. Together, our results suggest a possible mechanism for prenyl substrate transfer, where key aromatic residues facilitate h-cisPT binding to the ER membrane in an orientation that holds the substrate-containing active site near the membrane surface. Molecular dynamics simulations of the mutant exhibiting lower FRET show greater orientational variability relative to wild type. This evidence for a specific orientation of h-cisPT provides a structural basis for isoprenoid association to the membrane during synthesis and prior to its release.

Serum-Stable, Cationic, α-Helical AMPs to Combat Infections of ESKAPE Pathogens and C. albicans

Expedition in the rate of development of antimicrobial resistance accompanied by the slowdown in the development of new antimicrobials has led to a dire necessity to develop an alternate class of antimicrobial agents. Antimicrobial peptides (AMPs), available in nature, are effective molecules that can combat microbial infections. However, due to several inherent shortcomings such as salt sensitivity of their potency, short systemic half-lives owing to protease and serum degradation, and cytotoxicity, their commercial success is limited. Inspired by α helical AMPs present in nature, here in this work, we have developed two short, cationic, helical AMPs RR-12 and FL-13. Both peptides exhibited high broad-spectrum antimicrobial activity, salt tolerance, prompt bactericidal activity, considerable serum stability, remaining non-cytotoxic and non-hemolytic at relevant microbicidal concentrations. The designed AMPs were membranolytic toward the microbial strains, though there were subtle differences in the mechanism owing to the variation in the composition of the cell membranes in different microbes. Rigorous experimental techniques and molecular dynamics (MD) simulations were performed to understand the structure, activity, and their mechanisms in detail. Positive charge, balanced hydrophobicity-hydrophilicity, and helical conformation were the different attributes that led to the development of the superior performance of the AMPs, making them valuable additions to the repertoire of therapeutically promising antimicrobials.

Benzimidazole-Derived B2 as a Fluorescent Probe for Bacterial Outer Membrane Vesicle (OMV) Labeling: Integrating DFT, Molecular Dynamics, Flow Cytometry, and Confocal Microscopy

Bacterial outer membrane vesicles (OMVs) are nanoscale extracellular structures produced by Gram-negative bacteria that are critical for microbial biology and host-pathogen interactions and have great potential in biotechnological applications. Despite the availability of fluorescent dyes for OMV studies, many are repurposed from eukaryotic extracellular vesicle research and are not explicitly optimized for OMVs, leading to challenges in achieving consistent labeling, minimizing background noise, and preserving vesicle integrity during analyses. This study evaluates B2, a benzimidazole-derived fluorophore, for OMV labeling in advanced techniques like flow cytometry and confocal microscopy. OMVs were isolated from Escherichia coli strains BL21 and O157, and their integrity was confirmed using transmission electron microscopy (TEM). B2 staining protocols were optimized for OMVs, and fluorescence analyses revealed specific interactions with the vesicle membranes, reducing aggregation and enhancing signal uniformity. Flow cytometry indicated near-complete labeling efficiency (98-100%) with minimal background interference. Confocal microscopy further validated B2's effectiveness, showing evident OMV internalization into epithelial HT-29 cells and compatibility with other fluorophores. Density functional theory (DFT) calculations, including Fukui function analysis, identified key electrophilic and nucleophilic regions in B2 that facilitate specific hydrogen bonding and polar interactions with membrane components. Non-covalent interaction (NCI) analysis revealed pronounced intramolecular hydrogen bonding along with discrete regions of weak van der Waals interactions. Molecular dynamics simulations suggest that B2 exhibits an affinity for both the hydrophobic core of the lipid bilayer and the core oligosaccharide region of the LPS layer, which collectively ensures sustained retention of the dye. The findings presented in this study position B2 as a valuable fluorophore for OMV research.

Evolutionary responses of Escherichia coli to phage pressure: insights into mucoidy and colanic acid overexpression

BACKGROUND

Antibiotic resistance is a major issue affecting all spheres of human activity, including agriculture. One significant example is the Avian Pathogenic Escherichia coli (APEC), a bacterium that infects poultry and leads to substantial economic losses in the farming industry. As antibiotics lose efficacity, bacteriophages (phages) -viruses that specifically target bacteria-are emerging as a promising alternative to antibiotics for treating and preventing bacterial infections. However, bacteria can develop resistance to phages through various mechanisms. Studying the coevolution between a phage and its host bacterium is important to gain insight into the phage's potential as a therapeutic agent. This study investigates the evolutionary responses of an APEC strain and a laboratory E. coli strain to a commercial phage originally isolated from APEC.

RESULTS

In most cases, phage resistance resulted in a significant increase in mucoidy. Genomic analysis revealed that this resistance consistently correlated with amino acid changes, particularly in proteins involved in colanic acid production, such as YrfF. Further investigation of a mutation found in the YrfF protein demonstrated that this mutation altered the protein's structure and its interaction with the membrane. Transcriptomic analysis confirmed that the genes involved in colanic acid production were significantly overexpressed. Although the strains possessed a CRISPR-Cas system, it did not contribute to phage resistance.

CONCLUSIONS

This study suggests that specific amino acid changes in key proteins may be a mechanism employed by E. coli, including APEC, to defend against phage infections.

Atomic Insights Into Self-Assembly of Zingibroside R1 and its Therapeutic Action Against Fungal Diseases

Natural products are a crucial resource for drug discovery, but poor understanding of the molecular-scale mechanisms of their self-assembly into soluble, bioavailable hydrogels limits their applications and therapeutic potential. It is demonstrated that Zingibroside R1 (ZR1), derived from Panax notoginseng, undergoes spontaneous self-assemble into a hydrogel comprising helical nanofibrils with potent antifungal activity lacking in its monomeric state. Cryogenic electron microscopy (cryo-EM) revealed an intricate hydrogen-bonding network that facilitates ZR1 nanofibril formation, characterized by a hydrophobic core and hydrophilic exterior architecture, which underpin its binding activity with cell wall in the vulvovaginal candidiasis (VVC) pathogen, C. albicans. The hydrogen-bonding interface between ZR1 gel and glucan compromises membrane integrity, inhibiting C. albicans proliferation in vitro and in VVC model mice in vivo. ZR1 gel could also deliver probiotic Lactobacillus, synergistically inhibiting VVC and restoring the vaginal microenvironment. This study advances the mechanistic understanding of ZR1's structure-function relationships, offering valuable insights into the rational design and therapeutic optimization of natural product-based hydrogels.

Ergosterol depletion by fish AMP analogs likely enhances fungal membrane permeability

The molecular interactions between a fungal membrane model and SJGAP, a 32-amino-acid antimicrobial peptide (AMP) derived from skipjack tuna GAPDH, as well as four analogs, were investigated using molecular dynamics simulations and Fourier transform infrared (FTIR) spectroscopy. In a previous study, Analog 7, modified by replacing three alanine residues with leucine residues, exhibited unique antifungal activity without any antibacterial effect. This contrasts with other analogs, which showed both antifungal and antibacterial effects. In this study, Analog 7 displayed the strongest interactions with the membrane's hydrophobic core, inserting more deeply and causing significantly greater membrane deformation and thinning compared with the other analogs. Its presence caused significant membrane deformation, evident from the displacement of both the phosphate groups and terminal methyls of the lipids. Notably, Analog 7 was the only analog to induce a marked depletion of ergosterol around the peptide insertion site. FTIR spectroscopy experiments further confirmed the distinctive impact of Analog 7 on a fungal membrane model. The combined results from molecular dynamics simulations and spectroscopy emphasize the critical role of leucine substitutions in Analog 7, particularly at residues 18 and 19 within the central α helix, in promoting membrane thinning and inducing ergosterol depletion, suggesting increased membrane permeabilization, which could explain its previously reported antifungal specificity. This study provides the first insights into the molecular interactions between a GAPDH-derived AMP and a fungal membrane model, offering valuable information about its antifungal mechanism of action.

Modeling the different conformations of the human mitochondrial ADP/ATP carrier using AlphaFold and molecular dynamics simulations of the protein-ligand complexes

The ADP/ATP Carrier (AAC), a member of the mitochondrial Solute Carrier Family 25 (SLC25), facilitates the exchange of cytosolic ADP for mitochondrial ATP across the inner mitochondrial membrane (IMM). It serves as a master regulator of the cellular ADP/ATP ratio and is involved in various pathologies, including cancer. Its transport mechanism involves a conformational transition that alternates the accessibility of the binding site between the cytoplasmic (c-state) and mitochondrial (m-state) sides of the IMM. In this study, the human AAC was used as a case study to evaluate the performance of AlphaFold2 (AF2) and AlphaFold3 (AF3) for structural modeling of members of the SLC25 family. The study also compared the AF3 approach for predicting protein-ligand complexes with the standard methodology of modeling followed by molecular docking. Both AF2 and AF3 display a bias toward the c-state conformation. On the other hand, ColabFold implementation of AF2 successfully generated the first ab initio structural model of the human AAC in the m-state conformation. Modeling of the complexes coupled to molecular dynamics (MD) simulations allowed to obtain structural insight into AAC's substrate binding and stabilization mechanisms, and the possible effects of pathogenic mutations on its conformational dynamics and functionality. These analyses provided a deeper understanding of AAC's alternating access mechanism and highlighted the potential of AF3 in modeling protein-ligand interactions, though only in the c-state. This work demonstrates the reliability of AlphaFold models when aligned with experimental data and provides further confirmation of their utility for investigating solute carriers and membrane proteins.

A quantitative analysis of ligand binding at the protein-lipid bilayer interface

The majority of drugs target membrane proteins, and many of these proteins contain ligand binding sites embedded within the lipid bilayer. However, targeting these therapeutically relevant sites is hindered by limited characterization of both the sites and the molecules that bind to them. Here, we introduce the Lipid-Interacting LigAnd Complexes Database (LILAC-DB), a curated dataset of 413 structures of ligands bound at the protein-bilayer interface. Analysis of these structures reveals that ligands binding to lipid-exposed sites exhibit distinct chemical properties, such as higher calculated partition coefficient (clogP), molecular weight, and a greater number of halogen atoms, compared to ligands that bind to soluble proteins. Additionally, we demonstrate that the atomic properties of these ligands vary significantly depending on their depth within and exposure to the lipid bilayer. We also find that ligand binding sites exposed to the bilayer have distinct amino acid compositions compared to other protein regions, which may aid in the identification of lipid-exposed binding sites. This analysis provides valuable guidelines for researchers pursuing structure-based drug discovery targeting underexploited ligand binding sites at the protein-lipid bilayer interface.

Anomalous membrane organization by omega-6 and omega-9 fatty acids

Omega fatty acids are currently being marketed as healthy food supplements as they have been implicated in multiple pathophysiological conditions, such as reducing plaque formation of Aβ peptide and inhibiting SARS-CoV-2 infection. Their mode of action has been hypothesized to be via membrane reorganization by the unsaturated acyl chains, leading to the modulation of lipid-protein cross-talk. However, the lack of molecular details led us to evaluate the molecular effect of omega-6 (linolenic acid) and omega-9 (oleic acid) fatty acids on membrane organization using a consolidated approach of fluorescence spectroscopy and all-atom molecular dynamics simulation. Our results show that the effect of these omega fatty acids is sensitive to their protonation states. Contrary to the accepted notion that chain unsaturation causes membrane disordering, both experimental and simulation results demonstrate that protonated linoleic acid promotes membrane ordering, despite having two unsaturations at the fatty acyl chain. However, protonated oleic fatty acid, with reduced unsaturation, disordered the acyl chain area of the lipid membranes. Equally surprisingly, deprotonated oleic acid orders, whereas deprotonated linoleic acid disorders, the membrane core region. Interestingly, while the lipid order parameter measurements from simulations did not capture these subtle differences, the calculated rotational autocorrelation function of a membrane dye was in line with experimentally measured apparent rotational correlation times. Our work provides a comprehensive revised molecular picture of the effect of omega fatty acids on membranes and highlights the importance of rigorous comparative approaches, as experimental and simulation studies in isolation can sometimes lead to inconsistent results.

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.

Examining the Biophysical Properties of the Inner Membrane of Gram-Negative ESKAPE Pathogens

The World Health Organization has identified multidrug-resistant bacteria as a serious global health threat. Gram-negative bacteria are particularly prone to antibiotic resistance, and their high rate of antibiotic resistance has been suggested to be related to the complex structure of their cell membrane. The outer membrane of Gram-negative bacteria contains lipopolysaccharides that protect the bacteria against threats such as antibiotics, while the inner membrane houses 20-30% of the bacterial cellular proteins. Given the cell membrane's critical role in bacterial survival, antibiotics targeting the cell membrane have been proposed to combat bacterial infections. However, a deeper understanding of the biophysical properties of the bacterial cell membrane is crucial to developing effective and specific antibiotics. In this study, Martini coarse-grain molecular dynamics simulations were used to investigate the interplay between membrane composition and biophysical properties of the inner membrane across four pathogenic bacterial species: Klebsiella pneumoniae, Pseudomonas aeruginosa, Enterobacter cloacae, and Escherichia coli. The simulations indicate the impact of species-specific membrane composition on the overall membrane properties. Specifically, the cardiolipin concentration in the inner membrane is a key factor influencing the membrane features. Model membranes with varying concentrations of bacterial lipids (phosphatidylglycerol, phosphatidylethanolamine, and cardiolipin) further support the significant role of cardiolipin in determining the membrane biophysical properties. The bacterial inner membrane models developed in this work pave the way for future simulations of bacterial membrane proteins and for simulations investigating novel strategies aimed at disrupting the bacterial membrane to treat antibiotic-resistant infections.

Characterization of Key Lipid Components in the Cell Membrane of Freeze-Drying Resistant Lacticaseibacillus paracasei Strains Using Nontargeted Lipidomics

Lactic acid bacteria (LAB) are usually freeze-dried into powder for transportation and storage, with the bacterial membrane playing a crucial role in this process. However, different strains exhibit different levels of freeze-drying resistance in their cell membranes. In this study, Lacticaseibacillus paracasei (L. paracasei) strains 1F20, K56, and J5, demonstrating survival rates of 59.51, 25.86, and 4.05% after freeze-drying, respectively, were selected. The membrane structure and composition of these strains were subsequently analyzed. Bacterial live/dead staining results indicated that strain 1F20 maintained the highest membrane integrity after drying. Nontargeted lipidomics analysis revealed six differential lipid species that differed in membrane lipid compositions. KEGG functional enrichment analysis revealed 13 significantly different pathways, with glycerophospholipid metabolism being the most critical. This study explored the membrane composition of L. paracasei at the cellular level and identified key lipid species associated with freeze-drying resistance, providing a reference for screening highly resistant strains.

Split Membrane: A New Model to Accelerate All-Atom MD Simulation of Phospholipid Bilayers

All-atom molecular dynamics simulations are powerful tools for studying cell membranes and their interactions with proteins and other molecules. However, these processes occur on time scales determined by the diffusion rate of phospholipids, which are challenging to achieve in all-atom models. Here, we present a new all-atom model that accelerates lipid diffusion by splitting phospholipid molecules into head and tail groups. The bilayer structure is maintained by using external lateral potentials, which compensate for the lipid split. This split model enhances lateral lipid diffusion more than ten times, allowing faster and cheaper equilibration of large systems with different phospholipid types. The current model has been tested on membranes containing PSM, POPC, POPS, POPE, POPA, and cholesterol. We have also evaluated the interaction of the split model membranes with the Disheveled DEP domain and amphiphilic helix motif of the transcriptional repressor Opi1 as representative of peripheral proteins as well as the dimeric fragment of the epidermal growth factor receptor transmembrane domain and the Human A2A Adenosine of G protein-coupled receptors as representative of transmembrane proteins. The split model can predict the interaction sites of proteins and their preferred phospholipid type. Thus, the model could be used to identify lipid binding sites and equilibrate large membranes at an affordable computational cost.

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.

Structural basis for membrane association and catalysis by phosphatidylserine synthase in Escherichia coli

Phosphatidylserine synthase (PssA) is essential in the biosynthesis of phosphatidylethanolamine, a major phospholipid of bacterial membranes. A peripheral membrane protein PssA can associate with the cellular membrane in its active state or exist in the cytosol in an inactive form. The membrane-bound enzyme acts on cytidine diphosphate diacylglycerol (CDP-DG) to form cytidine monophosphate and a covalent intermediate, which is subsequently targeted by serine to produce phosphatidylserine. Here, we present two crystal structures of Escherichia coli PssA, one complexed with CDP-DG and the other without. The lipid-bound structure mimics the Michaelis complex before the formation of a covalent intermediate, revealing key determinants for substrate recognition and catalysis. Notably, membrane-free PssA is in a monomer-dimer equilibrium, with only the monomer capable of associating with the membrane, suggesting a regulatory mechanism for phospholipid biosynthesis dependent on the oligomerization state of the enzyme.

Functional Roles of the Charged Residues of the C- and M-Gates in the Yeast Mitochondrial NAD+ Transporter Ndt1p

Mitochondrial carriers transport organic acids, amino acids, nucleotides and cofactors across the mitochondrial inner membrane. These transporters consist of a three-fold symmetric bundle of six transmembrane α-helices that encircle a pore with a central substrate binding site, whose alternating access is controlled by a cytoplasmic and a matrix gate (C- and M-gates). The C- and M-gates close by forming two different salt-bridge networks involving the conserved motifs [YF][DE]XX[KR] on the even-numbered and PX[DE]XX[KR] on the odd-numbered transmembrane α-helices, respectively. We have investigated the effects on transport of mutating the C-gate charged residues of the yeast NAD+ transporter Ndt1p and performed molecular docking with NAD+ and other substrates into structural models of Ndt1p. Double-cysteine substitutions and swapping the positions of the C-gate charged-pair residues showed that all of them contribute to the high transport rate of wild-type Ndt1p, although no single salt bridge is essential for activity. The in silico docking results strongly suggest that both the C-gate motif mutations and our previously reported M-gate mutations affect gate closing, whereas those of the M-gate also affect substrate binding, which is further supported by molecular dynamics. In particular, NAD+ most likely interferes with the cation-π interaction between R303-W198, which has been proposed to exist in the Ndt1p M-gate in the place of one of the salt bridges. These findings contribute to understanding the roles of the charged C- and M-gate residues in the transport mechanism of Ndt1p.

Acidic sphingomyelinase interactions with lysosomal membranes and cation amphiphilic drugs: A molecular dynamics investigation

Lysosomes are pivotal in cellular functions and disease, influencing cancer progression and therapy resistance with Acid Sphingomyelinase (ASM) governing their membrane integrity. Moreover, cation amphiphilic drugs (CADs) are known as ASM inhibitors and have anti-cancer activity, but the structural mechanisms of their interactions with the lysosomal membrane and ASM are poorly explored. Our study, leveraging all-atom explicit solvent molecular dynamics simulations, delves into the interaction of glycosylated ASM with the lysosomal membrane and the effects of CAD representatives, i.e., ebastine, hydroxyebastine and loratadine, on the membrane and ASM. Our results confirm the ASM association to the membrane through the saposin domain, previously only shown with coarse-grained models. Furthermore, we elucidated the role of specific residues and ASM-induced membrane curvature in lipid recruitment and orientation. CADs also interfere with the association of ASM with the membrane at the level of a loop in the catalytic domain engaging in membrane interactions. Our computational approach, applicable to various CADs or membrane compositions, provides insights into ASM and CAD interaction with the membrane, offering a valuable tool for future studies.

Screening and computational characterization of novel antimicrobial cathelicidins from amphibian transcriptomic data

As cold-blooded organisms living in damp and dark environments, amphibians have evolved robust defense mechanisms to protect themselves from predators and infections. Among the wide repertoire of bioactive compounds they produce are antimicrobial peptides (AMPs), which are required as part of innate immunity. One important class of AMPs is cathelicidins, known for their broad-spectrum activity against pathogens and their immunoregulatory roles. However, despite their promising biomedical potential and the increasing availability of omics data, few cathelicidins have been studied in amphibians, mostly through conventional experimental techniques. Here, we present 210 novel cathelicidin sequences from amphibian transcriptomes, identified through a comprehensive computational pipeline, which employed HMMER and BLAST tools to screen cathelicidin domains. These sequences reveal a typical tripartite domain architecture that was confirmed by SignalP and InterProScan analysis. Phylogenetic inference with IQ-TREE classified the sequences into six categories based on evolutionary relationships. Compared to cathelicidins from other vertebrates, amphibian mature peptides exhibit longer average lengths (around 50 amino acids), fewer aromatic and hydrophobic residues, and reduced thermal stability. Furthermore, these amphibian cathelicidins were characterized for their physicochemical and biological properties, revealing significant antimicrobial potential with lower hemolytic capability, especially in anurans, which suggests a balance between their antimicrobial and hemolytic activities predicted through AMPlify, ampir, AmpGram, and HemoPI. Secondary structure estimations, including three-dimensional modeling using AlphaFold2, indicate that amphibian cathelicidins predominantly feature α-helices and coils. Some representative models also display a high α-helix composition with amphipathic topology, facilitating interactions with simulated bacterial membranes as assessed by the PPM approach. Thus, these findings highlight the functional role of cathelicidins in amphibian immunity and their promising biomedical applicability, emphasizing the importance of applying computational methods to expand the scope and reveal the diverse landscape of cathelicidins across vertebrates.

Dynamic Mechanism for Subtype Selectivity of Endocannabinoids

Endocannabinoids are naturally occurring lipid-like molecules that bind to cannabinoid receptors (CB1 and CB2) and regulate many of human bodily functions via the endocannabinoid system. There is a tremendous interest in developing selective drugs that target the CB receptors. However, the biophysical mechanisms responsible for the subtype selectivity for endocannbinoids have not been established. Recent experimental structures of CB receptors show that endocannbinoids potentially bind via membrane using the lipid access channel in the transmembrane region of the receptors. Furthermore, the N-terminus of the receptor could move in and out of the binding pocket thereby modulating both the pocket volume and its residue composition. On the basis of these observations, we propose two hypothesis to explain the selectivity of the endocannabinoid, anandamide for CB1 receptor. First, the selectivity arises from distinct enthalpic ligand-protein interactions along the ligand binding pathway formed due to the movement of N-terminus and subsequent shifts in the binding pocket composition. Second, selectivity arises from the volumetric differences in the binding pocket allowing for differences in ligand conformational entropy. To quantitatively test these hypothesis, we perform extensive molecular dynamics simulations (∼0.9 milliseconds) along with Markov state modeling and deep learning-based VAMP-nets to provide an interpretable characterization of the anandamide binding process to cannabinoid receptors and explain its selectivity for CB1. Our findings reveal that the distinct N-terminus positions along lipid access channels between TM1 and TM7 lead to different binding mechanisms and interactions between anandamide and the binding pocket residues. To validate the critical stabilizing interactions along the binding pathway, relative free energy calculations of anandamide analogs are used. Moreover, the larger CB2 pocket volume increases the entropic effects of ligand binding by allowing higher ligand fluctuations but reduced stable interactions. Therefore, the opposing enthalpy and entropy effects between the receptors shape the endocannabinoid selectivity. Overall, the CB1 selectivity of anandamide is explained by the dominant enthalpy contributions due to ligand-protein interactions in stable binding poses. This study shed lights on potential selectivity mechanisms for endocannabinoids that would aid in the discovery of CB selective drugs.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

ASM variants in the spotlight: A structure-based atlas for unraveling pathogenic mechanisms in lysosomal acid sphingomyelinase

Lysosomal acid sphingomyelinase (ASM), a critical enzyme in lipid metabolism encoded by the SMPD1 gene, plays a crucial role in sphingomyelin hydrolysis in lysosomes. ASM deficiency leads to acid sphingomyelinase deficiency, a rare genetic disorder with diverse clinical manifestations, and the protein can be found mutated in other diseases. We employed a structure-based framework to comprehensively understand the functional implications of ASM variants, integrating pathogenicity predictions with molecular insights derived from a molecular dynamics simulation in a lysosomal membrane environment. Our analysis, encompassing over 400 variants, establishes a structural atlas of missense variants of lysosomal ASM, associating mechanistic indicators with pathogenic potential. Our study highlights variants that influence structural stability or exert local and long-range effects at functional sites. To validate our predictions, we compared them to available experimental data on residual catalytic activity in 135 ASM variants. Notably, our findings also suggest applications of the resulting data for identifying cases suited for enzyme replacement therapy. This comprehensive approach enhances the understanding of ASM variants and provides valuable insights for potential therapeutic interventions.

Extreme makeover: the incredible cell membrane adaptations of extremophiles to harsh environments

The existence of life beyond Earth has long captivated humanity, and the study of extremophiles-organisms surviving and thriving in extreme environments-provides crucial insights into this possibility. Extremophiles overcome severe challenges such as enzyme inactivity, protein denaturation, and damage of the cell membrane by adopting several strategies. This feature article focuses on the molecular strategies extremophiles use to maintain the cell membrane's structure and fluidity under external stress. Key strategies include homeoviscous adaptation (HVA), involving the regulation of lipid composition, and osmolyte-mediated adaptation (OMA), where small organic molecules protect the lipid membrane under stress. Proteins also have direct and indirect roles in protecting the lipid membrane. Examining the survival strategies of extremophiles provides scientists with crucial insights into how life can adapt and persist in harsh conditions, shedding light on the origins of life. This article examines HVA and OMA and their mechanisms in maintaining membrane stability, emphasizing our contributions to this field. It also provides a brief overview of the roles of proteins and concludes with recommendations for future research directions.

Computational insights into the mechanisms underlying structural destabilization and recovery in trafficking-deficient hERG mutants

Cardiovascular diseases are a major global health concern, responsible for a significant number of deaths each year, often linked to cardiac arrhythmias resulting from dysfunction in ion channels. Hereditary Long QT Syndrome (LQTS) is a condition characterized by a prolonged QT interval on ECG, increasing the risk of sudden cardiac death. The most common type of LQTS, LQT2, is caused by mutations in the hERG gene, affecting a potassium ion channel. The majority of these mutations disrupt the channel's trafficking to the cell membrane, leading to intracellular retention. Specific high-affinity hERG blockers (e.g., E-4031) can rescue this mutant phenotype, but the exact mechanism is unknown. This study used accelerated molecular dynamics simulations to investigate how these mutations affect the hERG channel's structure, folding, endoplasmic reticulum (ER) retention, and trafficking. We reveal that these mutations induce structural changes in the channel, narrowing its central pore and altering the conformation of the intracellular domains. These changes expose internalization signals that contribute to ER retention and degradation of the mutant hERG channels. Moreover, the study found that the trafficking rescue drug E-4031 can inhibit these structural changes, potentially rescuing the mutant channels. This research offers valuable insights into the structural issues responsible for the degradation of rescuable transmembrane trafficking mutants. Understanding the defective trafficking structure of the hERG channel could help identify binding sites for small molecules capable of restoring proper folding and facilitating channel trafficking. This knowledge has the potential to lead to mechanism-based therapies that address the condition at the cellular level, which may prove more effective than treating clinical symptoms, ultimately offering hope for individuals with hereditary Long QT Syndrome.

Discovery of Antibacterial Compounds with Potential Multi-Pharmacology against Staphylococcus Mur ligase Family Members by In Silico Structure-Based Drug Screening

Staphylococcus aureus (S. aureus) is a major bacterial infection in humans, leading to severe disease and causing death. The stagnation of antibiotic development in recent decades has made it difficult to combat drug-resistant infections. In this study, we performed an in silico structure-based drug screening (SBDS) targeting the S. aureus MurE (saMurE) enzyme involved in cell wall synthesis of S. aureus. saMurE is an enzyme that is essential for the survival of S. aureus but not present in humans. SBDS identified nine saMurE inhibitor candidates, Compounds 1-9, from a structural library of 154,118 compounds. Among them, Compound 2 showed strong antibacterial activity against Staphylococcus epidermidis (S. epidermidis) used as a model bacterium. Amino acid sequence homology between saMurE and S. epidermidis MurE is 87.4%, suggesting that Compound 2 has a similar inhibitory effect on S. aureus. Compound 2 showed an IC50 value of 301 nM for S. epidermidis in the dose-dependent growth inhibition assay. Molecular dynamics simulation showed that Compound 2 binds stably to both S. aureus MurD and S. aureus MurF, suggesting that it is a potential multi-pharmacological pharmacological inhibitor. The structural and bioactivity information of Compound 2, as well as its potential multiple-target activity, could contribute to developing new antimicrobial agents based on MurE inhibition.

Development of protease resistant and non-cytotoxic Jelleine analogs with enhanced broad spectrum antimicrobial efficacy

Short systemic half- life of Antimicrobial Peptides (AMP) is one of the major bottlenecks that limits their successful commercialization as therapeutics. In this work, we have designed analogs of the natural AMP Jelleine, obtained from royal jelly of apis mellifera. Among the designed peptides, J3 and J4 were the most potent with broad spectrum activities against a varied class of ESKAPE pathogens and fungus C. albicans. All the developed peptides were more effective against Gram-negative bacteria in comparison to the Gram-positive pathogens, and were especially effective against P. aeruginosa and C. albicans.J3 and J4 were completely trypsin resistant and serum stable, while retaining the non-cytotoxicity of the parent Jelleine, Jc. The designed peptides were membranolytic in their mode of action. CD and MD simulations in the presence of bilayers, established that J3 and J4 were non-structured even upon membrane binding and suggested that biological properties of the AMPs were innocent of any specific secondary structural requirements. Enhancement of charge to increase the antimicrobial potency, controlling the hydrophobic-hydrophilic balance to maintain non-cytotoxicity and induction of unnatural amino acid residues to impart protease resistance, remains some of the fundamental principles in the design of more effective antimicrobial therapeutics of the future, which may help combat the quickly rising menace of antimicrobial resistance in the microbes.

CHARMM-GUI Membrane Builder for Lipid Droplet Modeling and Simulation

Lipid droplets (LDs) are organelles that are necessary for eukaryotic and prokaryotic metabolism and energy storage. They have a unique structure consisting of a spherical phospholipid monolayer encasing neutral lipids such as triacylglycerol (TAG). LDs have garnered increased interest for their implications in disease and for drug delivery applications. Consequently, there is an increased need for tools to study their structure, composition, and dynamics in biological contexts. In this work, we utilize CHARMM-GUI Membrane Builder to simulate and analyze LDs with and without a plant LD protein, oleosin. The results show that Membrane Builder can generate biologically relevant all-atom LD systems with relatively short equilibration times using a new TAG library having optimized headgroup parameters. TAG molecules originally inserted into a lipid bilayer aggregate in the membrane center, forming a TAG-only core flanked by two monolayers. The TAG-only core thickness stably grows with increasing TAG mole fraction. A 70 % TAG system has a core that is thick enough to house oleosin without its interactions with the distal leaflet or disruption of its secondary structure. We hope that Membrane Builder can aid in the future study of LD systems, including their structure and dynamics with and without proteins.

A facile method for monitoring sphingomyelin synthase activity in HeLa cells using liquid chromatography/mass spectrometry

Sphingomyelin synthase (SMS) is a sphingolipid-metabolizing enzyme involved in the de novo synthesis of sphingomyelin (SM) from ceramide (Cer). Recent studies have indicated that SMS is a key therapeutic target for metabolic diseases such as fatty liver, type 2 diabetes, atherosclerosis, and colorectal cancer. However, very few SMS inhibitors have been identified because of the limited sensitivity and selectivity of the current fluorescence-based screening assay. In this study, we developed a simple cell-based assay coupled with liquid chromatography/tandem mass spectrometry (LC-MS/MS) to screen for SMS inhibitors. HeLa cells stably expressing SMS1 or SMS2 were used for the screening. A non-fluorescent unnatural C6-Cer was used as a substrate for SMS to produce C6-SM. C6-Cer and C6-SM levels in the cells were monitored and quantified using LC-MS/MS. The activity of ginkgolic acid C15:1 (GA), a known SMS inhibitor, was measured. GA had half-maximal inhibitory concentrations of 5.5 μM and 3.6 μM for SMS1 and SMS2, respectively. To validate these findings, hSMS1 and hSMS2 proteins were optimized for molecular docking studies. In silico analyses were conducted to assess the interaction of GA with SMS1 and SMS2, and its binding affinity. This study offers an analytical approach for screening novel SMS inhibitors and provides in silico support for the experimental findings.

Membrane specificity of the human cholesterol transfer protein STARD4

STARD4 regulates cholesterol homeostasis by transferring cholesterol between the plasma membrane and endoplasmic reticulum. The STARD4 structure features a helix-grip fold surrounding a large hydrophobic cavity holding the sterol. Its access is controlled by a gate formed by the Ω1 and Ω4 loops and the C-terminal α-helix. Little is known about the mechanisms by which STARD4 binds to membranes and extracts/releases cholesterol. All available structures of STARD4 are without a bound sterol and display the same closed conformation of the gate. The cholesterol transfer activity of the mouse STARD4 is enhanced in the presence of anionic lipids, and in particular of phosphatidylinositol biphosphates (PIP2) for which two binding sites were proposed on the mouse STARD4 surface. Yet only one of these sites is conserved in human STARD4. We here report the results of a liposome microarray-based assay and microseconds-long molecular dynamics simulations of human STARD4 with complex lipid bilayers mimicking the composition of the donor and acceptor membranes. We show that the binding of apo form of human STARD4 is sensitive to the presence of PIP2 through two specific binding sites, one of which was not identified on mouse STARD4. We report two novel conformations of the gate in holo-STARD4: a yet-unobserved close conformation and an open conformation of Ω4 shedding light on the opening/closure mechanism needed for cholesterol uptake/release. Overall, the modulation of human STARD4 membrane-binding by lipid composition, and by the presence of the cargo supports the capacity of human STARD4 to achieve directed transfer between specific organelle membranes.

Recognition and remodeling of endosomal zones by sorting nexins

The proteolipid code determines how cytosolic proteins find and remodel membrane surfaces. Here, we investigate how this process works with sorting nexins Snx1 and Snx3. Both proteins form sorting machines by recognizing membrane zones enriched in phosphatidylinositol 3-phosphate (PI3P), phosphatidylserine (PS) and cholesterol. This co-localized combination forms a unique "lipid codon" or lipidon that we propose is responsible for endosomal targeting, as revealed by structures and interactions of their PX domain-based readers. We outline a membrane recognition and remodeling mechanism for Snx1 and Snx3 involving this code element alongside transmembrane pH gradients, dipole moment-guided docking and specific protein-protein interactions. This generates an initial membrane-protein assembly (memtein) that then recruits retromer and additional PX proteins to recruit cell surface receptors for sorting to the trans-Golgi network (TGN), lysosome and plasma membranes. Post-translational modification (PTM) networks appear to regulate how the sorting machines form and operate at each level. The commonalities and differences between these sorting nexins show how the proteolipid code orchestrates parallel flows of molecular information from ribosome emergence to organelle genesis, and illuminates a universally applicable model of the membrane.

Modeling asymmetric cell membranes at all-atom resolution

Molecular dynamics (MD) simulations are a useful tool when studying the properties of membranes as they allow for a molecular view of lipid interactions with proteins, nucleic acids, or small molecules. While model membranes are usually symmetric in their lipid composition between leaflets and include a small number of lipid components, physiological membranes are highly complex and vary in the level of asymmetry. Simulation studies have shown that changes in leaflet asymmetry can alter the properties of a membrane. It is therefore necessary to carefully build asymmetric membranes to accurately simulate membranes. This chapter carefully describes the different methods for building asymmetric membranes and the advantages/disadvantages of each method. The simplest methods involve building a membrane with either an equal number of lipids per leaflet or an equal initial surface area (SA) estimated by the area per lipid. More detailed methods include combining two symmetric membranes of equal SA or altering an asymmetric membrane and adjusting the number of lipids after equilibration to minimize an observable such as differential stress (0-DS). More complex methods that require specific simulation software are also briefly described. The challenges and assumptions are listed for each method which should help guide the researcher to choose the best method for their unique MD simulation of an asymmetric membrane.

Rendering protein structures inside cells at the atomic level with Unreal Engine

While the recent development of cryogenic electron tomography (CryoET) makes it possible to identify various macromolecules inside cells and determine their structure at near-atomic resolution, it remains challenging to visualize the complex cellular environment at the atomic level. One of the main hurdles in cell visualization is to render the millions of molecules in real time computationally. Here, using a video game engine, we demonstrate the capability of rendering massive biological macromolecules at the atomic level within their native environment. To facilitate the visualization, we also provide tools that help the interactive navigation inside the cells, as well as software that converts protein structures identified using CryoET to a scene that can be explored with the game engine.

Promiscuity and Quantitative Contribution of UGT2B17 in Drug and Steroid Metabolism Determined by Experimental and Computational Approaches

Uridine 5'-diphospho-glulcuronosyltransferase 2B17 (UGT2B17) is important in the metabolism of steroids and orally administered drugs due to its high interindividual variability. However, the structural basis governing the substrate selectivity or inhibition of UGT2B17 remains poorly understood. This study investigated 76 FDA-approved drugs and 20 steroids known to undergo glucuronidation for their metabolism by UGT2B17. Specifically, we assessed the substrate selectivity for UGT2B17 over other UGT enzymes using recombinant human UGT2B17 (rUGT2B17), human intestinal microsomes, and human liver microsomes. The quantitative contribution of intestinal UGT2B17 in the glucuronidation of these compounds was characterized using intestinal microsomes isolated from UGT2B17 expressors and nonexpressors. In addition, a structure-based pharmacophore model for UGT2B17 substrates was built and validated using the studied pool of substrates and nonsubstrates. The results show that UGT2B17 could metabolize 23 out of 96 compounds from various chemical classes, including alcohols and carboxylic acids, particularly in the intestine. Interestingly, amines were less susceptible to UGT2B17 metabolism, though they could inhibit the enzyme. Three main pharmacophoric features of UGT2B17 substrates include (1) the presence of an accessible -OH or -COOH group near His35 residue, (2) a hydrophobic functional group at ∼4.5-5 Å from feature 1, and (3) an aromatic ring ∼5-7 Å from feature 2. Most of the studied compounds inhibited UGT2B17 activity irrespective of their substrate potential, indicating the possibility of multiple mechanisms. These data suggest that UGT2B17 is promiscuous in substrate selectivity and inhibition and has a high potential to produce significant variability in the absorption and disposition of orally administered drugs.

Dynamic conformational changes of a tardigrade group-3 late embryogenesis abundant protein modulate membrane biophysical properties

A number of intrinsically disordered proteins (IDPs) encoded in stress-tolerant organisms, such as tardigrade, can confer fitness advantage and abiotic stress tolerance when heterologously expressed. Tardigrade-specific disordered proteins including the cytosolic-abundant heat-soluble proteins are proposed to confer stress tolerance through vitrification or gelation, whereas evolutionarily conserved IDPs in tardigrades may contribute to stress tolerance through other biophysical mechanisms. In this study, we characterized the mechanism of action of an evolutionarily conserved, tardigrade IDP, HeLEA1, which belongs to the group-3 late embryogenesis abundant (LEA) protein family. HeLEA1 homologs are found across different kingdoms of life. HeLEA1 is intrinsically disordered in solution but shows a propensity for helical structure across its entire sequence. HeLEA1 interacts with negatively charged membranes via dynamic disorder-to-helical transition, mainly driven by electrostatic interactions. Membrane interaction of HeLEA1 is shown to ameliorate excess surface tension and lipid packing defects. HeLEA1 localizes to the mitochondrial matrix when expressed in yeast and interacts with model membranes mimicking inner mitochondrial membrane. Yeast expressing HeLEA1 shows enhanced tolerance to hyperosmotic stress under nonfermentative growth and increased mitochondrial membrane potential. Evolutionary analysis suggests that although HeLEA1 homologs have diverged their sequences to localize to different subcellular organelles, all homologs maintain a weak hydrophobic moment that is characteristic of weak and reversible membrane interaction. We suggest that such dynamic and weak protein-membrane interaction buffering alterations in lipid packing could be a conserved strategy for regulating membrane properties and represent a general biophysical solution for stress tolerance across the domains of life.

Molecular Modeling and Simulation of the Mycobacterial Cell Envelope: From Individual Components to Cell Envelope Assemblies

Unlike typical Gram-positive bacteria, the cell envelope of mycobacteria is unique and composed of a mycobacterial outer membrane, also known as the mycomembrane, a peptidoglycan layer, and a mycobacterial inner membrane, which is analogous to that of Gram-negative bacteria. Despite its importance, however, our understanding of this complex cell envelope is rudimentary at best. Thus, molecular modeling and simulation of such an envelope can benefit the scientific community by proposing new hypotheses about the biophysical properties of its different layers. In this Perspective, we present recent advances in molecular modeling and simulation of the mycobacterial cell envelope from individual components to cell envelope assemblies. We also show how modeling other types of cell envelopes, such as that of Escherichia coli, may help modeling part of the mycobacterial envelopes. We hope that the studies presented here are just the beginning of the road and more and more new modeling and simulation studies help us to understand crucial questions related to mycobacteria such as antibiotic resistance or bacterial survival in the host.

A Model of the Full-Length Cytokinin Receptor: New Insights and Prospects

Cytokinins (CK) are one of the most important classes of phytohormones that regulate a wide range of processes in plants. A CK receptor, a sensor hybrid histidine kinase, was discovered more than 20 years ago, but the structural basis for its signaling is still a challenge for plant biologists. To date, only two fragments of the CK receptor structure, the sensory module and the receiver domain, were experimentally resolved. Some other regions were built up by molecular modeling based on structures of proteins homologous to CK receptors. However, in the long term, these data have proven insufficient for solving the structure of the full-sized CK receptor. The functional unit of CK receptor is the receptor dimer. In this article, a molecular structure of the dimeric form of the full-length CK receptor based on AlphaFold Multimer and ColabFold modeling is presented for the first time. Structural changes of the receptor upon interacting with phosphotransfer protein are visualized. According to mathematical simulation and available data, both types of dimeric receptor complexes with hormones, either half- or fully liganded, appear to be active in triggering signals. In addition, the prospects of using this and similar models to address remaining fundamental problems of CK signaling were outlined.

Synthesis of Pyrrolo[3,4-b]pyridin-5-ones via Ugi-Zhu Reaction and In Vitro-In Silico Studies against Breast Carcinoma

An Ugi-Zhu three-component reaction (UZ-3CR) coupled in a one-pot manner to a cascade process (N-acylation/aza Diels-Alder cycloaddition/decarboxylation/dehydration) was performed to synthesize a series of pyrrolo[3,4-b]pyridin-5-ones in 20% to 92% overall yields using ytterbium triflate as a catalyst, toluene as a solvent, and microwaves as a heat source. The synthesized molecules were evaluated in vitro against breast cancer cell lines MDA-MB-231 and MCF-7, finding that compound 1f, at a concentration of 6.25 μM, exhibited a potential cytotoxic effect. Then, to understand the interactions between synthesized compounds and the main proteins related to the cancer cell lines, docking studies were performed on the serine/threonine kinase 1 (AKT1) and Orexetine type 2 receptor (Ox2R), finding moderate to strong binding energies, which matched accurately with the in vitro results. Additionally, molecular dynamics were performed between proteins related to the studied cell lines and the three best ligands.

Structural Modeling of Cytokine-Receptor-JAK2 Signaling Complexes Using AlphaFold Multimer

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). These cell-surface single-pass transmembrane (TM) glycoproteins regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains, and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all of the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated here by using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts, with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices, causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to the TM α-helices of the active TPOR dimer was proposed. The models also help elucidate the molecular basis of oncogenic mutations that may involve a noncanonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Key Role of Choline Head Groups in Large Unilamellar Phospholipid Vesicles for the Interaction with and Rupture by Silica Nanoparticles

For highly abundant silica nanomaterials, detrimental effects on proteins and phospholipids are postulated as critical molecular initiating events that involve hydrogen-bonding, hydrophobic, and/or hydrophilic interactions. Here, large unilamellar vesicles with various well-defined phospholipid compositions are used as biomimetic models to recapitulate membranolysis, a process known to be induced by silica nanoparticles in human cells. Differential analysis of the dominant phospholipids determined in membranes of alveolar lung epithelial cells demonstrates that the quaternary ammonium head groups of phosphatidylcholine and sphingomyelin play a critical and dose-dependent role in vesicle binding and rupture by amorphous colloidal silica nanoparticles. Surface modification by either protein adsorption or by covalent coupling of carboxyl groups suppresses the disintegration of these lipid vesicles, as well as membranolysis in human A549 lung epithelial cells by the silica nanoparticles. Furthermore, molecular modeling suggests a preferential affinity of silanol groups for choline head groups, which is also modulated by the pH value. Biomimetic lipid vesicles can thus be used to better understand specific phospholipid-nanoparticle interactions at the molecular level to support the rational design of safe advanced materials.

Marine diterpenoid targets STING palmitoylation in mammalian cells

Natural products are important sources of therapeutic agents and useful drug discovery tools. The fused macrocycles and multiple stereocenters of briarane-type diterpenoids pose a major challenge to total synthesis and efforts to characterize their biological activities. Harnessing a scalable source of excavatolide B (excB) from cultured soft coral Briareum stechei, we generated analogs by late-stage diversification and performed structure-activity analysis, which was critical for the development of functional excB probes. We further used these probes in a chemoproteomic strategy to identify Stimulator of Interferon Genes (STING) as a direct target of excB in mammalian cells. We showed that the epoxylactone warhead of excB is required to covalently engage STING at its membrane-proximal Cys91, inhibiting STING palmitoylation and signaling. This study reveals a possible mechanism-of-action of excB, and expands the repertoire of covalent STING inhibitors.

Machine Learning Methods for Small Data Challenges in Molecular Science

Small data are often used in scientific and engineering research due to the presence of various constraints, such as time, cost, ethics, privacy, security, and technical limitations in data acquisition. However, big data have been the focus for the past decade, small data and their challenges have received little attention, even though they are technically more severe in machine learning (ML) and deep learning (DL) studies. Overall, the small data challenge is often compounded by issues, such as data diversity, imputation, noise, imbalance, and high-dimensionality. Fortunately, the current big data era is characterized by technological breakthroughs in ML, DL, and artificial intelligence (AI), which enable data-driven scientific discovery, and many advanced ML and DL technologies developed for big data have inadvertently provided solutions for small data problems. As a result, significant progress has been made in ML and DL for small data challenges in the past decade. In this review, we summarize and analyze several emerging potential solutions to small data challenges in molecular science, including chemical and biological sciences. We review both basic machine learning algorithms, such as linear regression, logistic regression (LR), k-nearest neighbor (KNN), support vector machine (SVM), kernel learning (KL), random forest (RF), and gradient boosting trees (GBT), and more advanced techniques, including artificial neural network (ANN), convolutional neural network (CNN), U-Net, graph neural network (GNN), Generative Adversarial Network (GAN), long short-term memory (LSTM), autoencoder, transformer, transfer learning, active learning, graph-based semi-supervised learning, combining deep learning with traditional machine learning, and physical model-based data augmentation. We also briefly discuss the latest advances in these methods. Finally, we conclude the survey with a discussion of promising trends in small data challenges in molecular science.

Structural modeling of cytokine-receptor-JAK2 signaling complexes using AlphaFold Multimer

Homodimeric class 1 cytokine receptors include the erythropoietin (EPOR), thrombopoietin (TPOR), granulocyte colony-stimulating factor 3 (CSF3R), growth hormone (GHR), and prolactin receptors (PRLR). They are cell-surface single-pass transmembrane (TM) glycoproteins that regulate cell growth, proliferation, and differentiation and induce oncogenesis. An active TM signaling complex consists of a receptor homodimer, one or two ligands bound to the receptor extracellular domains and two molecules of Janus Kinase 2 (JAK2) constitutively associated with the receptor intracellular domains. Although crystal structures of soluble extracellular domains with ligands have been obtained for all the receptors except TPOR, little is known about the structure and dynamics of the complete TM complexes that activate the downstream JAK-STAT signaling pathway. Three-dimensional models of five human receptor complexes with cytokines and JAK2 were generated using AlphaFold Multimer. Given the large size of the complexes (from 3220 to 4074 residues), the modeling required a stepwise assembly from smaller parts with selection and validation of the models through comparisons with published experimental data. The modeling of active and inactive complexes supports a general activation mechanism that involves ligand binding to a monomeric receptor followed by receptor dimerization and rotational movement of the receptor TM α-helices causing proximity, dimerization, and activation of associated JAK2 subunits. The binding mode of two eltrombopag molecules to TM α-helices of the active TPOR dimer was proposed. The models also help elucidating the molecular basis of oncogenic mutations that may involve non-canonical activation route. Models equilibrated in explicit lipids of the plasma membrane are publicly available.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Phase Transition and Phase Separation in Realistic Thylakoid Lipid Membrane of Marine Algae in All-Atom Simulations

Thylakoid membranes are specialized membranes predominantly composed of uncommon galacto- and sulfolipids, having distinct roles in photosynthesis. Large acyl chain variety and richness in polyunsaturated fatty acid (PUFA) content of thylakoid lipids further add to the compositional complexity. The function of these membrane systems is intimately dependent on the fluidity of its lipid matrix, which is strongly modulated by the lipid composition and temperature. The present work, employing extensive atomistic simulations, provides the first atomistic view of the phase transition and domain coexistence in a model membrane composed of thylakoid lipids of a commercially important red alga Gracilaria corticata between 10 and 40 °C. The growth and photosynthetic activity of marine algae are greatly influenced by the seawater temperature. So far, little is known about the molecular organization of lipids in thylakoid membranes, in particular their adaptive arrangements under temperature stress. Our simulations show that the algal thylakoid membrane undergoes a transition from a gel-like phase at a low temperature, 10-15 °C, to a homogeneous liquid-crystalline phase at a high temperature, 40 °C. Clear evidence of spontaneous phase separation into coexisting nanoscale domains is detected at intermediate temperatures nearing the optimal growth temperature range. Particularly, at 25-30 °C, we identified the formation of a stable ripple phase, where the gel-like domains rich in saturated and nearly hexagonally packed lipids were separated from fluid-like domains enriched in lipids containing PUFA chains. The phase separation is driven by the spontaneous and preferential segregation of lipids into differentially ordered domains, mainly depending on the acyl chain types. Cholesterol impairs the phase transition and the emergence of domains and induces a fairly uniform liquid-ordered phase in the membrane over the temperatures studied. This work improves the understanding of the properties and reorganization of lipids in the thylakoid membrane in response to temperature variation.

Ceramide-1-phosphate transfer protein enhances lipid transport by disrupting hydrophobic lipid-membrane contacts

Cellular distributions of the sphingolipid ceramide-1-phosphate (C1P) impact essential biological processes. C1P levels are spatiotemporally regulated by ceramide-1-phosphate transfer protein (CPTP), which efficiently shuttles C1P between organelle membranes. Yet, how CPTP rapidly extracts and inserts C1P into a membrane remains unknown. Here, we devise a multiscale simulation approach to elucidate biophysical details of CPTP-mediated C1P transport. We find that CPTP binds a membrane poised to extract and insert C1P and that membrane binding promotes conformational changes in CPTP that facilitate C1P uptake and release. By significantly disrupting a lipid's local hydrophobic environment in the membrane, CPTP lowers the activation free energy barrier for passive C1P desorption and enhances C1P extraction from the membrane. Upon uptake of C1P, further conformational changes may aid membrane unbinding in a manner reminiscent of the electrostatic switching mechanism used by other lipid transfer proteins. Insertion of C1P into an acceptor membrane, eased by a decrease in membrane order by CPTP, restarts the transfer cycle. Most notably, we provide molecular evidence for CPTP's ability to catalyze C1P extraction by breaking hydrophobic C1P-membrane contacts with compensatory hydrophobic lipid-protein contacts. Our work, thus, provides biophysical insights into how CPTP efficiently traffics C1P between membranes to maintain sphingolipid homeostasis and, additionally, presents a simulation method aptly suited for uncovering the catalytic mechanisms of other lipid transfer proteins.

In silico modeling human VPS13 proteins associated with donor and target membranes suggests lipid transfer mechanisms

The VPS13 protein family constitutes a novel class of bridge-like lipid transferases. Autosomal recessive inheritance of mutations in VPS13 genes is associated with the development of neurodegenerative diseases in humans. Bioinformatic approaches previously recognized the domain architecture of these proteins. In this study, we model the first ever full-length structures of the four human homologs VPS13A, VPS13B, VPS13C, and VPS13D in association with model membranes, to investigate their lipid transfer ability and potential structural association with membrane leaflets. We analyze the evolutionary conservation and physicochemical properties of these proteins, focusing on conserved C-terminal amphipathic helices that disturb organelle surfaces and that, adjoined, resemble a traditional Venetian gondola. The gondola domains share significant structural homology with lipid droplet surface-binding proteins. We introduce in silico protein-membrane models displaying the mode of association of VPS13A, VPS13B, VPS13C, and VPS13D to donor and target membranes, and present potential models of action for protein-mediated lipid transfer.

In Silico Analysis of the Structural Dynamics and Substrate Recognition Determinants of the Human Mitochondrial Carnitine/Acylcarnitine SLC25A20 Transporter

The Carnitine-Acylcarnitine Carrier is a member of the mitochondrial Solute Carrier Family 25 (SLC25), known as SLC25A20, involved in the electroneutral exchange of acylcarnitine and carnitine across the inner mitochondrial membrane. It acts as a master regulator of fatty acids β-oxidation and is known to be involved in neonatal pathologies and cancer. The transport mechanism, also known as "alternating access", involves a conformational transition in which the binding site is accessible from one side of the membrane or the other. In this study, through a combination of state-of-the-art modelling techniques, molecular dynamics, and molecular docking, the structural dynamics of SLC25A20 and the early substrates recognition step have been analyzed. The results obtained demonstrated a significant asymmetry in the conformational changes leading to the transition from the c- to the m-state, confirming previous observations on other homologous transporters. Moreover, analysis of the MD simulations' trajectories of the apo-protein in the two conformational states allowed for a better understanding of the role of SLC25A20 Asp231His and Ala281Val pathogenic mutations, which are at the basis of Carnitine-Acylcarnitine Translocase Deficiency. Finally, molecular docking coupled to molecular dynamics simulations lend support to the multi-step substrates recognition and translocation mechanism already hypothesized for the ADP/ATP carrier.

γ-Core Guided Antibiotic Design Based on Human Enteric Defensin 5

An increase in the number of infections caused by resistant bacteria worldwide necessitates the development of alternatives to antibiotics. Human defensin (HD) 5 is an innate immune peptide with broad-spectrum antibacterial activity, but its complicated structure makes its preparation difficult. Herein, we truncated the HD5 structure by extracting the highly conserved γ-core motif. A structure-activity study showed that this motif was ineffective in killing bacteria in the absence of specific spatial conformation. Notably, after the introduction of two intramolecular disulfide bonds, its antibacterial activity was markedly improved. Glu and Ser residues were then replaced with Arg to create the derivative RC18, which exhibited stronger potency than HD5, particularly against methicillin-resistant S. aureus (MRSA). Mechanistically, RC18 bound to lipid A and lipoteichoic acid at higher affinities than HD5. Furthermore, RC18 was more efficient than HD5 in penetrating the bacterial membranes. Molecular dynamics simulation revealed that five Arg residues, Arg1, Arg7, Arg9, Arg15, and Arg18, mediated most of the polar interactions of RC18 with the phospholipid head groups during membrane penetration. In vivo experiments indicated that RC18 decreased MRSA colonization and dramatically improved the survival of infected mice, thus demonstrating that RC18 is a promising drug candidate to treat MRSA infections.