CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model

Atomistic characterization of β2-glycoprotein I domain V interaction with anionic membranes

BACKGROUND

Interaction of β2-glycoprotein I (β2GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of its membrane-binding domain, domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood.

OBJECTIVES

This study aimed to elucidate the molecular details of β2GPI DV binding to anionic membranes.

METHODS

We utilized molecular dynamics simulations to investigate the structural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the highly mobile membrane mimetic simulations, we conducted additional simulations using a full membrane model.

RESULTS

The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, which are essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β2GPI's membrane interaction. Simulations revealed 2 distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β2GPI by antiphospholipid antibodies.

CONCLUSION

The study advances our understanding of the molecular interactions between β2GPI's DV and anionic membranes, which are crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant binding mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

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.

Salt-Induced Hydrophobic C-Terminal Region of α-Synuclein Triggers Its Fibrillation under the Mimic Physiologic Condition

α-Synuclein (αS) causes Parkinson's disease due to the structural alteration into amyloid fibrils that form after the interaction with synaptic membranes in neurons. To understand the alternation mechanism, the effect of salt (NaCl) on the interaction of αS with synaptic mimic membrane was characterized at the molecular level because salt triggered the amyloid fibril formation. The membrane-bound conformation (or the initial conformation before fibrillation) showed that NaCl decreased the number of helical structures and Tyr residues interacting with the membrane surface compared to when NaCl was absent, implying an increase in solvent-exposed regions. The N-terminal region of αS interacted with the membrane, forming the helical structures regardless of NaCl, while the C-terminal region formed a random structure with weak membrane interaction, but NaCl inhibited the interaction of its hydrophobic area, suggesting that salt promoted amyloid fibril formations by exposing the hydrophobic C-terminal region, which can intermolecularly interact with free αS.

Virtual screening, molecular docking and dynamics simulation studies to identify potential agonists of orphan receptor GPR78 targeting CNS disorders

G protein-coupled receptors (GPCRs) are important targets in drug discovery because of their roles in physiological and pathological processes. Orphan GPCRs are GPCR proteins for which endogenous ligands have not yet been identified and they present interesting avenues for therapeutic intervention. This study focuses on GPR78, an orphan GPCR that is expressed in the central nervous system and linked to neurological disorders. GPR78 has no reported crystal structure and there is limited research. In this study, we have predicted the three dimensional model of GPR78 and its probable binding pocket. Structure-based virtual screening was carried out using the ChemDiv and Enamine REAL databases, followed by induced-fit docking studies to identify potential lead compounds with favorable interactions. These lead compounds were then embedded into a POPC lipid bilayer for a 200 ns molecular dynamics simulation. Free energy landscapes and MM-PBSA analyses were performed to assess the binding energies and conformational dynamics. The results highlight the dynamic nature of GPR78 in the presence of lead compounds and show favorable binding interactions. This study aims to predict a reliable three dimensional model of GPR78 and identify novel lead compounds through a comprehensive in silico approach. The identification of these potential GPR78 agonists represents a significant step in the development of new therapeutics for neurological disorders, highlighting the therapeutic potential of orphan GPR78 in CNS disorders.

Faa1 membrane binding drives positive feedback in autophagosome biogenesis via fatty acid activation

Autophagy serves as a stress response pathway by mediating the degradation of cellular material within lysosomes. In autophagy, this material is encapsulated in double-membrane vesicles termed autophagosomes, which form from precursors referred to as phagophores. Phagophores grow by lipid influx from the endoplasmic reticulum into Atg9-positive compartments and local lipid synthesis provides lipids for their expansion. How phagophore nucleation and expansion are coordinated with lipid synthesis is unclear. Here, we show that Faa1, an enzyme activating fatty acids, is recruited to Atg9 vesicles by directly binding to negatively charged membranes with a preference for phosphoinositides such as PI3P and PI4P. We define the membrane-binding surface of Faa1 and show that its direct interaction with the membrane is required for its recruitment to phagophores. Furthermore, the physiological localization of Faa1 is key for its efficient catalysis and promotes phagophore expansion. Our results suggest a positive feedback loop coupling phagophore nucleation and expansion to lipid synthesis.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

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.

Atomistic Characterization of Beta-2-Glycoprotein I Domain V Interaction with Anionic Membranes

Background

Interaction of beta-2-glycoprotein I ( β 2 GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of it's membrane bind-ing domain, Domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood.

Objectives

This study aims to elucidate the mechanism by which DV of β 2 GPI binds to anionic membranes.

Methods

We utilized molecular dynamics (MD) simulations to investigate the struc-tural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the HMMM simulations, we conducted additional simulations using a full mem-brane model.

Results

The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β 2 GPI's membrane interaction. Simulations revealed two distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β 2 GPI by antiphospholipid antibodies.

Conclusion

The study advances our understanding of the molecular interactions be-tween β 2 GPI's DV and anionic membranes, crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant bind-ing mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

Structural basis for the oligomerization-facilitated NLRP3 activation

The NACHT-, leucine-rich-repeat-, and pyrin domain-containing protein 3 (NLRP3) is a critical intracellular inflammasome sensor and an important clinical target against inflammation-driven human diseases. Recent studies have elucidated its transition from a closed cage to an activated disk-like inflammasome, but the intermediate activation mechanism remains elusive. Here we report the cryo-electron microscopy structure of NLRP3, which forms an open octamer and undergoes a ~ 90° hinge rotation at the NACHT domain. Mutations on open octamer's interfaces reduce IL-1β signaling, highlighting its essential role in NLRP3 activation/inflammasome assembly. The centrosomal NIMA-related kinase 7 (NEK7) disrupts large NLRP3 oligomers and forms NEK7/NLRP3 monomers/dimers which is a critical step preceding the assembly of the disk-like inflammasome. These data demonstrate an oligomeric cooperative activation of NLRP3 and provide insight into its inflammasome assembly mechanism.

The carboxy terminus causes interfacial assembly of oleate hydratase on a membrane bilayer

The soluble flavoprotein oleate hydratase (OhyA) hydrates the 9-cis double bond of unsaturated fatty acids. OhyA substrates are embedded in membrane bilayers; OhyA must remove the fatty acid from the bilayer and enclose it in the active site. Here, we show that the positively charged helix-turn-helix motif in the carboxy terminus (CTD) is responsible for interacting with the negatively charged phosphatidylglycerol (PG) bilayer. Super-resolution microscopy of Staphylococcus aureus cells expressing green fluorescent protein fused to OhyA or the CTD sequence shows subcellular localization along the cellular boundary, indicating OhyA is membrane-associated and the CTD sequence is sufficient for membrane recruitment. Using cryo-electron microscopy, we solved the OhyA dimer structure and conducted 3D variability analysis of the reconstructions to assess CTD flexibility. Our surface plasmon resonance experiments corroborated that OhyA binds the PG bilayer with nanomolar affinity and we found the CTD sequence has intrinsic PG binding properties. We determined that the nuclear magnetic resonance structure of a peptide containing the CTD sequence resembles the OhyA crystal structure. We observed intermolecular NOE from PG liposome protons next to the phosphate group to the CTD peptide. The addition of paramagnetic MnCl2 indicated the CTD peptide binds the PG surface but does not insert into the bilayer. Molecular dynamics simulations, supported by site-directed mutagenesis experiments, identify key residues in the helix-turn-helix that drive membrane association. The data show that the OhyA CTD binds the phosphate layer of the PG surface to obtain bilayer-embedded unsaturated fatty acids.

Alchemical Free Energy Workflows for the Computation of Protein-Ligand Binding Affinities

Alchemical free energy methods can be used for the efficient computation of relative binding free energies during preclinical drug discovery stages. In recent years, this has been facilitated further by the implementation of workflows that enable non-experts to quickly and consistently set up the required simulations. Given the correct input structures, workflows handle the difficult aspects of setting up perturbations, including consistently defining the perturbable molecule, its atom mapping and topology generation, perturbation network generation, running of the simulations via different sampling methods, and analysis of the results. Different academic and commercial workflows are discussed, including FEW, FESetup, FEPrepare, CHARMM-GUI, Transformato, PMX, QLigFEP, TIES, ProFESSA, PyAutoFEP, BioSimSpace, FEP+, Flare, and Orion. These workflows differ in various aspects, such as mapping algorithms or enhanced sampling methods. Some workflows can accommodate more than one molecular dynamics (MD) engine and use external libraries for tasks. Differences between workflows can present advantages for different use cases, however a lack of interoperability of the workflows' components hinders systematic comparisons.

Myr-Arf1 conformational flexibility at the membrane surface sheds light on the interactions with ArfGAP ASAP1

ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Characterization of membrane-interaction mechanisms of proteins using vacuum-ultraviolet circular dichroism spectroscopy

Protein-membrane interactions play an important role in various biological phenomena, such as material transport, demyelinating diseases, and antimicrobial activity. We combined vacuum-ultraviolet circular dichroism (VUVCD) spectroscopy with theoretical (e.g., molecular dynamics and neural networks) and polarization experimental (e.g., linear dichroism and fluorescence anisotropy) methods to characterize the membrane interaction mechanisms of three soluble proteins (or peptides). α1 -Acid glycoprotein has the drug-binding ability, but the combination of VUVCD and neural-network method revealed that the membrane interaction causes the extension of helix in the N-terminal region, which reduces the binding ability. Myelin basic protein (MBP) is an essential component of the myelin sheath with a multi-layered structure. Molecular dynamics simulations using a VUVCD-guided system showed that MBP forms two amphiphilic and three non-amphiphilic helices as membrane interaction sites. These multivalent interactions may allow MBP to interact with two opposing membrane leaflets, contributing to the formation of a multi-layered myelin structure. The antimicrobial peptide magainin 2 interacts with the bacterial membrane, causing damage to its structure. VUVCD analysis revealed that the M2 peptides assemble in the membrane and turn into oligomers with a β-strand structure. Linear dichroism and fluorescence anisotropy suggested that the oligomers are inserted into the hydrophobic core of the membrane, disrupting the bacterial membrane. Overall, our findings demonstrate that VUVCD and its combination with theoretical and polarization experimental methods pave the way for unraveling the molecular mechanisms of biological phenomena related to protein-membrane interactions.

Lactadherin's Multistate Binding Predicts Stable Membrane-Bound Conformations of Factors V and VIII's C Domains

Protein binding to negatively charged lipids is essential for maintaining numerous vital cellular processes where its dysfunction can lead to various diseases. One such protein that plays a crucial role in this process is lactadherin, which competes with coagulation factors for membrane binding sites to regulate blood clotting. Despite identifying key binding regions of these proteins through structural and biochemical studies, models incorporating membrane dynamics are still lacking. In this study, we report on the multimodal binding of lactadherin and use it to gain insight into the binding mechanisms of its C domain homologs, factor V and factor VIII. Molecular dynamics simulations enhanced with the highly mobile mimetic model enabled the determination of lactadherin's multimodal binding on membranes that revealed critical interacting residues consistent with prior NMR and mutagenesis data. The binding occurred primarily via two dynamic structural ensembles: an inserted state and an unreported, highly conserved side-lying state driven by a cationic patch. We utilized these findings to analyze the membrane binding domains of coagulation factors V and VIII and identified their preferred membrane-bound conformations. Specifically, factor V's C domains maintained an inserted state, while factor VIII preferred a tilted, side-lying state that permitted antibody binding. Insight into lactadherin's atomistically resolved membrane interactions from a multistate perspective can guide new therapeutic opportunities in treating diseases related to blood coagulation.

Simulation study of membrane bending by protein crowding: a case study with the epsin N-terminal homology domain

The mechanisms by which peripheral membrane proteins generate curvature is currently an active area of research. One of the proposed mechanisms is amphipathic insertion or the 'wedge' mechanism in which the protein shallowly inserts an amphipathic helix inside the membrane to drive the curvature. However, recent experimental studies have challenged the efficiency of the 'wedge' mechanism as it requires unusual protein densities. These studies proposed an alternative mechanism, namely 'protein-crowding', in which the lateral pressure generated by the random collisions among the membrane bound proteins drives the bending. In this study, we employ atomistic and coarse-grained molecular dynamics simulations to investigate the effects of amphipathic insertion and protein crowding on the membrane surface. Considering epsin N-terminal homology (ENTH) domain as a model protein, we show that amphipathic insertion is not essential for membrane bending. Our results suggest that ENTH domains can aggregate on the membrane surface by employing another structured region (H₃ helix). And this protein crowding decreases the cohesive energy of the lipid tails which causes a significant decrease in the membrane bending rigidity. The ENTH domain can generate a similar degree of membrane curvature irrespective of the activity of its H₀ helix. Our results are consistent with the recent experimental results.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

A novel machine learning-based screening identifies statins as inhibitors of the calcium pump SERCA

We report a novel small-molecule screening approach that combines data augmentation and machine learning to identify Food and Drug Administration (FDA)-approved drugs interacting with the calcium pump (Sarcoplasmic reticulum Ca2+-ATPase, SERCA) from skeletal (SERCA1a) and cardiac (SERCA2a) muscle. This approach uses information about small-molecule effectors to map and probe the chemical space of pharmacological targets, thus allowing to screen with high precision large databases of small molecules, including approved and investigational drugs. We chose SERCA because it plays a major role in the excitation-contraction-relaxation cycle in muscle and it represents a major target in both skeletal and cardiac muscle. The machine learning model predicted that SERCA1a and SERCA2a are pharmacological targets for seven statins, a group of FDA-approved 3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitors used in the clinic as lipid-lowering medications. We validated the machine learning predictions by using in vitro ATPase assays to show that several FDA-approved statins are partial inhibitors of SERCA1a and SERCA2a. Complementary atomistic simulations predict that these drugs bind to two different allosteric sites of the pump. Our findings suggest that SERCA-mediated Ca2+ transport may be targeted by some statins (e.g., atorvastatin), thus providing a molecular pathway to explain statin-associated toxicity reported in the literature. These studies show the applicability of data augmentation and machine learning-based screening as a general platform for the identification of off-target interactions and the applicability of this approach extends to drug discovery.

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.

Effect of phosphorylation on the structural dynamics, thermal stability of human dopamine transporter: A simulation study using normal modes, molecular dynamics and Markov State Model

The Human Dopamine Transporter (hDAT) plays an essential role in modulating the Influx/Efflux of dopamine, and it is involved in the mechanism of certain neurodegenerative diseases such as Parkinson's disease. Several studies have reported important states for Dopamine transport: outward-facing open state (OFo), the outward-facing closed state (OFc), the holo-occluded state closed (holo), and the inward-facing open state (IFo). Furthermore, experimental assays have shown that different phosphorylation conditions in hDAT can affect the rate of dopamine absorption. We present a protocol using hybrid simulation methods to study the conformational dynamics and stability of states of hDAT under different phosphorylation sites. With this protocol, we explored the conformational space of hDAT, identified the states, and evaluated the free energy differences and the transition probabilities between them in each of the phosphorylation cases. We also presented the conformational changes and correlated them with those described in the literature. There is a thesis/hypothesis that the phosphorylation condition corresponding to NP-333 system (where all sites Ser/Thr from residue 2 to 62 and 254 to 613 are phosphorylated, except residue 333) would decrease the rate of dopamine transport from the extracellular medium to the intracellular medium by hDAT as previously described in the literature by Lin et al., 2003. Our results corroborated this thesis/hypothesis and the data reported. It is probably due to the affectation/changes/alteration of the conformational dynamics of this system that makes the intermediate states more likely and makes it difficult to initial states associated with the uptake of dopamine in the extracellular medium, corroborating the experimental results. Furthermore, our results showed that just single phosphorylation/dephosphorylation could alter intrinsic protein motions affecting the sampling of one or more states necessary for dopamine transport. In this sense, the modification of phosphorylation influences protein movements and conformational preferences, affecting the stability of states and the transition between them and, therefore, the transport.

Nanomaterials-Mediated Co-Stimulation of Toll-Like Receptors and CD40 for Antitumor Immunity

Toll-like receptors (TLRs) and CD40-related signaling pathways represent critical bridges between innate and adaptive immune responses. Here, an immunotherapy regimen that enables co-stimulation of TLR7/8- and CD40-mediated pathways is developed. TLR7/8 agonist resiquimod (R848) derived amino lipids, RAL1 and RAL2, are synthesized and formulated into RAL-derived lipid nanoparticles (RAL-LNPs). The RAL2-LNPs show efficient CD40 mRNA delivery to DCs both in vitro (90.8 ± 2.7%) and in vivo (61.3 ± 16.4%). When combined with agonistic anti-CD40 antibody, this approach can produce effective antitumor activities in mouse melanoma tumor models, thereby suppressing tumor growth, prolonging mouse survival, and establishing antitumor memory immunity. Overall, RAL2-LNPs provide a novel platform toward cancer immunotherapy by integrating innate and adaptive immunity.

Molecular View into Preferential Binding of the Factor VII Gla Domain to Phosphatidic Acid

Factor VII (FVII) is a serine protease with a key role in initiating the coagulation cascade. It is part of a family of vitamin K-dependent clotting proteins, which require vitamin K for formation of their specialized membrane-binding domains (Gla domains). Membrane binding of the FVII Gla domain is critical to the activity of FVII, mediating the formation of its complex with other clotting factors. While Gla domains among coagulation factors are highly conserved in terms of amino acid sequence and structure, they demonstrate differential binding specificity toward anionic lipids. Although most Gla domain-containing clotting proteins display a strong preference for phosphatidylserine (PS), it has been demonstrated that FVII and protein C instead bind preferentially to phosphatidic acid (PA). We have developed the first model of the FVII Gla domain bound to PA lipids in membranes containing PA, the highly mobile membrane mimetic model, which accelerates slow diffusion of lipids in molecular dynamics simulations and therefore facilitates the membrane binding process and enhances sampling of lipid interactions. Simulations were performed using atomic level molecular dynamics, requiring a fixed charge to all atoms. The overall charge assigned to each PA lipid for this study was -1. We also developed an additional model of the FVII Gla domain bound to a 1:1 PS/PC membrane and compared the modes of binding of PS and PA lipids to FVII, allowing us to identify potential PA-specific binding sites.

Characterization of folic acid-functionalized PLA-PEG nanomicelle to deliver Letrozole: A nanoinformatics study

Effective cancer treatment is currently the number one challenge to human health. To date, several treatment methods have been introduced for cancer cell targeting. Among the proposed new methods for attacking cancer cells, nanotechnology has attracted much attention. Hence, various nanocarriers have been developed for targeted delivery of available drugs and improve their effectiveness against malignant cells. The PLA-PEG functionalised with folic acid (PLA-PEG-FA) is one of the nanocarriers with a limited range of applications for targeting cancer cells. In this investigation, different types of in-silico methods such as molecular docking approach, molecular dynamics simulation and free energy calculations are employed to characterise the carriers studied. The effectiveness of PLA-PEG-FA and PLA-PEG in delivering Letrozole as an aromatase inhibitor in cancer cells is examined. It is found that in the presence of folic acid, the stability and cell membrane permeability of nanomicelle are increased. Therefore, PLA-PEG-FA can be considered as a versatile carrier that can increase the effectiveness of aromatase inhibitors (such as Letrozole) and reduce their side effects.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Investigation of the molecular and mechanistic basis for the regioselective metabolism of midazolam by cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is the most important P450 enzyme for drug metabolism and drug-drug interaction, due to it being responsible for the biotransformation of approximately 50% of clinically used drugs. Advance knowledge of the molecular and mechanistic basis of CYP3A4 regioselective metabolism is beneficial for understanding the production of metabolites, and may allow personalized metabolic pathways or designing pathway-specific therapeutics. In this work, we focus on investigating the ligand-receptor interactions, substrate conformational transition, and key factors regulating the specificity of metabolic pathways using midazolam (MDZ) as a probe. Here, three types of substrate-binding conformations related to the diversity of MDZ metabolites are identified. The results also suggest that an allosteric site for MDZ is located near the F'-helix, A-anchor, and C-terminal loop of CYP3A4. The presence of an effector in the allosteric site can accelerate the conformational transition of the substrate via modulating a "sandwich" structure, and may affect the proportion of metabolites at high substrate concentration. We hope that the results can improve the understanding of the CYP3A4 structure and function, and provide a new perspective for drug development.

Modeling the Membrane Binding Mechanism of a Lipid Transport Protein Osh4 to Single Membranes

All atom molecular dynamics simulations are utilized to unravel the binding mechanism of yeast oxysterol binding protein (Osh4) to model membranes with varying anionic lipid concentration using all-atom (AA) and the highly mobile membrane mimetic (HMMM) representations. For certain protein-lipid interactions, an improved force field description is used (CUFIX) to accurately describe lipid-protein electrostatic interactions. Our detailed computational studies have identified a single, β-crease orientated, membrane bound conformation of Osh4 for all anionic membranes. The penetration of the PHE-239 residue below the membrane phosphate plane is the characteristic signature of the membrane-bound state of Osh4. As the phenylalanine loop anchors itself deeply in the membrane; the other regions of the Osh4, namely, ALPS motif, ß6- ß7 loop, ß14- ß15 loop and ß16- ß17 loop, maximize their contact with the membrane. Furthermore, loose lipid packing and higher mobility of HMMM enables stronger association of ALPS motif with the membrane lipids through its hydrophobic surface and after the HMMM is converted to AA and equilibrated the binding is 2-3 times stronger compared to simulations started with the AA representation yielding the major importance of the ALPS motif to binding. Quantitative estimation of binding energy revealed that the phenylalanine loop plays a crucial role in stable membrane attachment of Osh4 and contributes significantly toward overall binding process. The CUFIX parameters provide a more balanced picture of hydrophobic and electrostatic interactions between the protein and the membrane which differs from our past work that showed salt bridges alone stabilized Osh4-membrane contact. Our study provides a comprehensive picture of the binding mechanism of Osh4 with model single membranes and thus, understanding of the initial interactions is important for elucidating the biological function of this protein to shuttle lipids between organelles.

Variability in the Spatial Structure of the Central Loop in Cobra Cytotoxins Revealed by X-ray Analysis and Molecular Modeling

Cobra cytotoxins (CTs) belong to the three-fingered protein family and possess membrane activity. Here, we studied cytotoxin 13 from Naja naja cobra venom (CT13Nn). For the first time, a spatial model of CT13Nn with both “water” and “membrane” conformations of the central loop (loop-2) were determined by X-ray crystallography. The “water” conformation of the loop was frequently observed. It was similar to the structure of loop-2 of numerous CTs, determined by either NMR spectroscopy in aqueous solution, or the X-ray method. The “membrane” conformation is rare one and, to date has only been observed by NMR for a single cytotoxin 1 from N. oxiana (CT1No) in detergent micelle. Both CT13Nn and CT1No are S-type CTs. Membrane-binding of these CTs probably involves an additional step—the conformational transformation of the loop-2. To confirm this suggestion, we conducted molecular dynamics simulations of both CT1No and CT13Nn in the Highly Mimetic Membrane Model of palmitoiloleoylphosphatidylglycerol, starting with their “water” NMR models. We found that the both toxins transform their “water” conformation of loop-2 into the “membrane” one during the insertion process. This supports the hypothesis that the S-type CTs, unlike their P-type counterparts, require conformational adaptation of loop-2 during interaction with lipid membranes.

Ataxia-linked SLC1A3 mutations alter EAAT1 chloride channel activity and glial regulation of CNS function

Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system (CNS). Excitatory Amino Acid Transporters (EAATs) regulate extracellular glutamate by transporting it into cells, mostly glia, to terminate neurotransmission and to avoid neurotoxicity. EAATs are also chloride (Cl-) channels, but the physiological role of Cl- conductance through EAATs is poorly understood. Mutations of human EAAT1 (hEAAT1) have been identified in patients with episodic ataxia type 6 (EA6). One mutation showed increased Cl- channel activity and decreased glutamate transport, but the relative contributions of each function of hEAAT1 to mechanisms underlying the pathology of EA6 remain unclear. Here we investigated the effects of five additional EA6-related mutations on hEAAT1 function in Xenopus laevis oocytes, and on CNS function in a Drosophila melanogaster model of locomotor behavior. Our results indicate that mutations resulting in decreased hEAAT1 Cl- channel activity but with functional glutamate transport can also contribute to the pathology of EA6, highlighting the importance of Cl- homeostasis in glial cells for proper CNS function. We also identified a novel mechanism involving an ectopic sodium (Na+) leak conductance in glial cells. Together, these results strongly support the idea that EA6 is primarily an ion channelopathy of CNS glia.

Identification of Structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

Cholesterol-Mediated Clustering of the HIV Fusion Protein gp41 in Lipid Bilayers

The envelope glycoprotein (Env) of the human immunodeficient virus (HIV-1) is known to cluster on the viral membrane surface to attach to target cells and cause membrane fusion for HIV-1 infection. However, the molecular structural mechanisms that drive Env clustering remain opaque. Here, we use solid-state NMR spectroscopy and molecular dynamics (MD) simulations to investigate nanometer-scale clustering of the membrane-proximal external region (MPER) and transmembrane domain (TMD) of gp41, the fusion protein component of Env. Using ¹⁹F solid-state NMR experiments of mixed fluorinated peptides, we show that MPER-TMD trimers form clusters with interdigitated MPER helices in cholesterol-containing membranes. Inter-trimer ¹⁹F-¹⁹F cross peaks, which are indicative of spatial contacts within ∼2 nm, are observed in cholesterol-rich virus-mimetic membranes but are suppressed in cholesterol-free model membranes. Water-peptide and lipid-peptide cross peaks in 2D ¹H-¹⁹F correlation spectra indicate that the MPER is well embedded in model phosphocholine membranes but is more exposed to the surface of the virus-mimetic membrane. These experimental results are reproduced in coarse-grained and atomistic molecular dynamics simulations, which suggest that the effects of cholesterol on gp41 clustering is likely via indirect modulation of the MPER orientation. Cholesterol binding to the helix-turn-helix region of the MPER-TMD causes a parallel orientation of the MPER with the membrane surface, thus allowing MPERs of neighboring trimers to interact with each other to cause clustering. These solid-state NMR data and molecular dynamics simulations suggest that MPER and cholesterol cooperatively govern the clustering of gp41 trimers during virus-cell membrane fusion.

Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Subtype-selective mechanisms of negative allosteric modulators binding to group I metabotropic glutamate receptors

Group I metabotropic glutamate receptors (mGlu₁ and mGlu₅) are promising targets for multiple psychiatric and neurodegenerative disorders. Understanding the subtype selectivity of mGlu₁ and mGlu₅ allosteric sites is essential for the rational design of novel modulators with single- or dual-target mechanism of action. In this study, starting from the deposited mGlu₁ and mGlu₅ crystal structures, we utilized computational modeling approaches integrating docking, molecular dynamics simulation, and efficient post-trajectory analysis to reveal the subtype-selective mechanism of mGlu₁ and mGlu₅ to 10 diverse drug scaffolds representing known negative allosteric modulators (NAMs) in the literature. The results of modeling identified six pairs of non-conserved residues and four pairs of conserved ones as critical features to distinguish the selective NAMs binding to the corresponding receptors. In addition, nine pairs of residues are beneficial to the development of novel dual-target NAMs of group I metabotropic glutamate receptors. Furthermore, the binding modes of a reported dual-target NAM (VU0467558) in mGlu₁ and mGlu₅ were predicted to verify the identified residues that play key roles in the receptor selectivity and the dual-target binding. The results of this study can guide rational structure-based design of novel NAMs, and the approach can be generally applicable to characterize the features of selectivity for other G-protein-coupled receptors.

Conformation of myelin basic protein bound to phosphatidylinositol membrane characterized by vacuum-ultraviolet circular-dichroism spectroscopy and molecular-dynamics simulations

The 18.5-kDa isoform of myelin basic protein (MBP) interacts with the membrane surface of the myelin sheath to construct its compact multilamellar structure. This study characterized the conformation of MBP in the membrane by measuring the vacuum-ultraviolet circular-dichroism (VUVCD) spectra of MBP in the bilayer liposome comprising the following essential lipid constituents of the myelin sheath: phosphatidylinositol (PI), phosphatidylinositol-4-phosphate (PIP), and phosphatidylinositol-4,5-bisphosphate (PIP2). The spectra of MBP exhibited the characteristic peaks of the helix structure in the presence of PI liposome, and the intensity increased markedly in the presence of PIP and PIP2 liposomes to show an isodichroic point. This suggests that the amount of the membrane-bound conformation of MBP enhanced due to the increased number of negative net charges on the liposome surfaces. Secondary-structure analysis revealed that MBP in the membrane comprised approximately 40% helix contents and eight helix segments. Molecular-dynamics (MD) simulations of the eight segments were conducted for 250 ns in the presence of PI membrane, which predicted two amphiphilic and three nonamphiphilic helices as the membrane-interaction sites. Further analysis of the distances of the amino-acid residues in each segment from the phosphate group suggested that the nonamphiphilic helices interact with the membrane surface electrostatically, while the amphiphilic ones invade the inside of the membrane to produce electrostatic and hydrophobic interactions. These results show that MBP can interact with the PI membrane via amphiphilic and nonamphiphilic helices under the control of a delicate balance between electrostatic and hydrophobic interactions.

Structural basis of complex formation between mitochondrial anion channel VDAC1 and Hexokinase-II

Complex formation between hexokinase-II (HKII) and the mitochondrial VDAC1 is crucial to cell growth and survival. We hypothesize that HKII first inserts into the outer membrane of mitochondria (OMM) and then interacts with VDAC1 on the cytosolic leaflet of OMM to form a binary complex. To systematically investigate this process, we devised a hybrid approach. First, we describe membrane binding of HKII with molecular dynamics (MD) simulations employing a membrane mimetic model with enhanced lipid diffusion capturing membrane insertion of its H-anchor. The insertion depth of the H-anchor was then used to derive positional restraints in subsequent millisecond-scale Brownian dynamics (BD) simulations to preserve the membrane-bound pose of HKII during the formation of the HKII/VDAC1 binary complex. Multiple BD-derived structural models for the complex were further refined and their structural stability probed with additional MD simulations, resulting in one stable complex. A major feature in the complex is the partial (not complete) blockade of VDAC1's permeation pathway, a result supported by our comparative electrophysiological measurements of the channel in the presence and absence of HKII. We also show how VDAC1 phosphorylation disrupts HKII binding, a feature that is verified by our electrophysiology recordings and has implications in mitochondria-mediated cell death.

Impact of PIP2 Lipids, Force Field Parameters, and Mutational Analysis on the Binding of the Osh4's α6-α7 Domain

All-atom molecular dynamics simulations are used with the highly mobile membrane mimetic method to study the α₆-α₇ peptide of the critical yeast Osh4 peripheral membrane protein. This research focuses on the impact of 1-palmitoyl-2-oleoyl-sn-glycero-phosphatidylinoside 4,5-bisphosphate (PIP2) lipids and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine on the protein's ability to bind to the membrane. Details of the binding mechanism are described qualitatively and quantitatively by measuring the position of the deepest residues, angle of the peptide during binding, root mean square deviation of the atomic positions within the peptide, and interaction energy, while changing variables, such as the force field used and the presence of the PIP2 lipids. The negatively charged PIP2 has a large head group that is a few Ångstroms above the main membrane phosphates enabling the PIP2 lipids to interact with the peptide before it binds deeper into the membrane. The PIP2 lipids can alter the position of the peptide during binding by recruiting charged residues on the α₇ helix, such as R344 and R347. Residues R347 and R344 are unusual because they are slightly out of the reach of the main membrane phosphates but optimally positioned to interact with the PIP2 lipids. The salt-bridge interactions can also typically occur between cationic peptide residues such as R314, K325, and K336. The force field interaction effect on peptide binding was also investigated by changing the standard CHARMM36m to an improved description between some amino acids and lipid moieties (Phys. Chem. Chem. Phys. 20, 8432-8449). This resulted in the total number of salt bridges and hydrogen bonds being drastically reduced, the interaction energy was also reduced, and there was more balance between electrostatic and nonpolar interactions, but the general bound structure is maintained. This work is an important initial step to understand the effect of the Osh4 protein on the membrane binding and to quantify the effect of PIP2 lipids on this domain.

Molecular dynamics simulation for mechanism revelation of the safety and nutrition of lipids and derivatives in food: State of the art

Molecular dynamics (MD) simulation has proved to be a powerful tool in the study of proteins, nucleic acids, lipids, and carbohydrates et al. in fields of health, nutrition, and food science. In particular, MD simulation has been employed in the investigation of various lipid systems such as triglycerides, phospholipid membranes, etc. Due to the continuous updating of computing resources and the development of new MD simulation methods and force field parameters, the simulation's time and size scale of lipids system has increased by several orders of magnitude. However, MD simulation cannot be used for systems invovle chemical reactions. These greatly limit its further application in the field of lipid research. This paper reviews the progress and development of MD simulation, especially for the application of MD simulation in different lipid systems. In this paper, MD simulation and its general workflow was briefly introduced firstly. Subsequently, the application of MD simulation in various lipid systems was reviewed in-depth. Finally, the limitation and future prospects of MD simulation in lipid research were also discussed. This review provided new insights into the investigation of MD simulation, and a novel thought for lipid study. We believe that MD simulation will exhibit more and more great advantages in the investigation of lipids in the future due to the development of novlel methods.

Hedgehog proteins create a dynamic cholesterol interface

During formation of the Hedgehog (Hh) signaling proteins, cooperative activities of the Hedgehog INTein (Hint) fold and Sterol Recognition Region (SRR) couple autoproteolysis to cholesterol ligation. The cholesteroylated Hh morphogens play essential roles in embryogenesis, tissue regeneration, and tumorigenesis. Despite the centrality of cholesterol in Hh function, the full structure of the Hint-SRR ("Hog") domain that attaches cholesterol to the last residue of the active Hh morphogen remains enigmatic. In this work, we combine molecular dynamics simulations, photoaffinity crosslinking, and mutagenesis assays to model cholesterolysis intermediates in the human Sonic Hedgehog (hSHH) protein. Our results provide evidence for a hydrophobic Hint-SRR interface that forms a dynamic, non-covalent cholesterol-Hog complex. Using these models, we suggest a unified mechanism by which Hh proteins can recruit, sequester, and orient cholesterol, and offer a molecular basis for the effects of disease-causing hSHH mutations.

Identification of a Potential Membrane-Targeting Sequence in the C-Terminus of the F Plasmid Segregation Protein SopA

Stable maintenance and partitioning of the 'Fertility' plasmid or the F plasmid in its host Escherichia coli require the function of a ParA superfamily of proteins known as SopA. The mechanism by which SopA mediates plasmid segregation is well studied. SopA is a nucleoid-binding protein and binds DNA in an ATP-dependent but sequence non-specific manner. ATP hydrolysis stimulated by the binding of the SopBC complex mediates the release of SopA from the nucleoid. Cycles of ATP-binding and hydrolysis generate an ATPase gradient that moves the plasmid through a chemophoresis force. Nucleoid binding of SopA thus assumes a central role in its plasmid-partitioning function. However, earlier work also suggests that the F plasmid can be partitioned into anucleate cells, thus implicating nucleoid independent partitioning. Interestingly, SopA is also reported to be associated with the inner membrane of the bacteria. Here, we report the identification of a possible membrane-targeting sequence, a predicted amphipathic helix, at the C-terminus of SopA. Molecular dynamics simulations indicate that the predicted amphipathic helical motif of SopA has weak affinity for membranes. Moreover, we experimentally show that SopA can associate with bacterial membranes, is detectable in the membrane fractions of bacterial lysates, and is sensitive to the membrane potential. Further, unlike the wild-type SopA, a deletion of the C-terminal 29 amino acids results in the loss of F plasmids from bacterial cells.

Molecular Dynamics Simulation of T10609C and C10676G Mutations of Mitochondrial ND4L Gene Associated With Proton Translocation in Type 2 Diabetes Mellitus and Cataract Patients

The mutation rate of mitochondrial DNA (mtDNA) is 17 times higher than nuclear DNA, and these mutations can cause mitochondrial disease in 1 of 10.000 people. The T10609C mutation was identified in type 2 diabetes mellitus (T2DM) patients and the C10676G mutation in cataract patients, with both mutations occurring in the ND4L gene of mtDNA that encodes ND4L protein. ND4L protein, a subunit of complex I in the respiratory complex, has been shown to play a role in the proton translocation process. The purpose of this study was to investigate the effect of both mutations on the proton translocation mechanism. Mutation mapping showed changes in amino acids M47T (T10609C) and C69W (C10676G). The 100 ns molecular dynamics (MD) simulations performed on native and mutants of ND4L-ND6 subunits. It is revealed that the native model had a similar proton translocation pathway to that of complex I from other organisms. Interestingly, the mutant M47T and C69W showed the interruption of the translocation pathway by a hydrogen bond formation between Glu34 and Tyr157. It is observed that the mutations were restricting the passage of water molecules through the transmembrane region. These results could help to develop the computational assay for the validation of a specific genetic biomarker for T2DM and cataracts.

Membrane surface recognition by the ASAP1 PH domain and consequences for interactions with the small GTPase Arf1

Adenosine diphosphate-ribosylation factor (Arf) guanosine triphosphatase-activating proteins (GAPs) are enzymes that need to bind to membranes to catalyze the hydrolysis of guanosine triphosphate (GTP) bound to the small GTP-binding protein Arf. Binding of the pleckstrin homology (PH) domain of the ArfGAP With SH3 domain, ankyrin repeat and PH domain 1 (ASAP1) to membranes containing phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂] is key for maximum GTP hydrolysis but not fully understood. By combining nuclear magnetic resonance, neutron reflectometry, and molecular dynamics simulation, we show that binding of multiple PI(4,5)P₂ molecules to the ASAP1 PH domain (i) triggers a functionally relevant allosteric conformational switch and (ii) maintains the PH domain in a well-defined orientation, allowing critical contacts with an Arf1 mimic to occur. Our model provides a framework to understand how binding of the ASAP1 PH domain to PI(4,5)P₂ at the membrane may play a role in the regulation of ASAP1.

Conserved Luminal C-Terminal Domain Dynamically Controls Interdomain Communication in Sarcolipin

Sarcolipin (SLN) mediates Ca²⁺ transport and metabolism in muscle by regulating the activity of the Ca²⁺ pump SERCA. SLN has a conserved luminal C-terminal domain that contributes to its functional divergence among homologous SERCA regulators, but the precise mechanistic role of this domain remains poorly understood. We used all-atom molecular dynamics (MD) simulations of SLN totaling 77.5 μs to show that the N- (NT) and C-terminal (CT) domains function in concert. Analysis of the MD simulations showed that serial deletions of the SLN C-terminus do not affect the stability of the peptide nor induce dissociation of SLN from the membrane but promote a gradual decrease in both the tilt angle of the transmembrane helix and the local thickness of the lipid bilayer. Mutual information analysis showed that the NT and CT domains communicate with each other in SLN and that interdomain communication is partially or completely abolished upon deletion of the conserved segment Tyr29-Tyr31 as well as by serial deletions beyond this domain. Phosphorylation of SLN at residue Thr5 also induces changes in the communication between the CT and NT domains, which thus provides additional evidence for interdomain communication within SLN. We found that interdomain communication is independent of the force field used and lipid composition, which thus demonstrates that communication between the NT and CT domains is an intrinsic functional feature of SLN. We propose the novel hypothesis that the conserved C-terminus is an essential element required for dynamic control of SLN regulatory function.

Membrane Interactions of Cy3 and Cy5 Fluorophores and Their Effects on Membrane-Protein Dynamics

Organic fluorophores, such as Cy3 and Cy5, have been widely used as chemical labels to probe the structure and dynamics of membrane proteins. Although a number of previous studies have reported on the possibility of some of the water-soluble fluorophores to interact with lipid bilayers, detailed fluorophore-lipid interactions and, more importantly, the potential effect of such interactions on the natural dynamics of the labeled membrane proteins have not been well studied. We have performed a large set of all-atom molecular dynamics simulations employing the highly mobile membrane mimetic model to describe spontaneous partitioning of the fluorophores into lipid bilayers with different lipid compositions. Spontaneous membrane partitioning of Cy3 and Cy5 fluorophores captured in these simulations proceeds in two steps. Electrostatic interaction between the fluorophores and the lipid headgroups facilitates the initial, fast membrane association of the fluorophores, followed by slow insertion of hydrophobic moieties into the lipid bilayer core. After the conversion of the resulting membrane-bound systems to full-membrane representations, biased-exchange umbrella sampling simulations are performed for free energy calculations, revealing a higher energy barrier for partitioning into negatively charged (phosphatidylserine or phosphatidylcholine) membranes than purely zwitterionic (phosphatidylcholine or phosphatidylethanolamine) ones. Furthermore, the potential effect of fluorophore-lipid interactions on membrane proteins has been examined by covalently linking Cy5 to single- and multipass transmembrane helical proteins. Equilibrium simulations show strong position-dependent effects of Cy5-tagging on the structure and natural dynamics of membrane proteins. Interactions between the tagged protein and Cy5 were also observed. Our results suggest that fluorophore-lipid interactions can affect the structure and dynamics of membrane proteins to various extents, especially in systems with higher structural flexibility.

Delivery of Alpha-Mangostin Using Cyclodextrins through a Biological Membrane: Molecular Dynamics Simulation

α-Mangostin (MGS) exhibits various pharmacological activities, including antioxidant, anticancer, antibacterial, and anti-inflammatory properties. However, its low water solubility is the major obstacle for its use in pharmaceutical applications. To increase the water solubility of MGS, complex formation with beta-cyclodextrins (βCDs), particularly with the native βCD and/or its derivative 2,6-dimethyl-β-CD (DMβCD) is a promising technique. Although there have been several reports on the adsorption of βCDs on the lipid bilayer, the release of the MGS/βCDs inclusion complex through the biological membrane remains unclear. In this present study, the release the MGS from the two different βCDs (βCD and DMβCD) across the lipid bilayer was investigated. Firstly, the adsorption of the free MGS, free βCDs, and inclusion complex formation was studied by conventional molecular dynamics simulation. The MGS in complex with those two βCDs was able to spontaneously release free MGS into the inner membrane. However, both MGS and DMβCD molecules potentially permeated into the deeper region of the interior membrane, whereas βCD only adsorbed at the outer membrane surface. The interaction between secondary rim of βCD and the 1-palmitoeyl-2-oleoyl-glycero-3-phosphocholine (POPC) phosphate groups showed the highest number of hydrogen bonds (up to 14) corresponding to the favorable location of βCD on the POPC membrane. Additionally, the findings suggested that electrostatic energy was the main driving force for βCD adsorption on the POPC membrane, while van der Waals interactions played a predominant role in DMβCD adsorption. The release profile of MGS from the βCDs pocket across the lipid bilayer exhibited two energy minima along the reaction coordinate associated with the permeation of the MGS molecule into the deeper region of the POPC membrane.

Molecular Simulation of Mechanical Properties and Membrane Activities of the ESCRT-III Complexes

The endosomal sorting complex required for transport (ESCRT) machinery carries out the membrane scission reactions that are required for many biological processes throughout cells. How ESCRTs bind and deform cellular membranes and ultimately produce vesicles has been a matter of active research in recent years. In this study, we use fully atomistic molecular dynamics simulations to scrutinize the structural details of a filament composed of Vps32 protomers, a major component of ESCRT-III complexes. The simulations show that both hydrophobic and electrostatic interactions between monomers help maintain the structural stability of the filament, which exhibits an intrinsic bend and twist. Our findings suggest that the accumulation of bending and twisting stresses as the filament elongates on the membrane surface likely contributes to the driving force for membrane invagination. The filament exposes a large cationic surface that senses the negatively charged lipids in the membrane, and the N-terminal amphipathic helix of the monomers not only acts as a membrane anchor but also generates significant positive membrane curvature. Taking all results together, we discuss a plausible mechanism for membrane invagination driven by ESCRT-III.

Microscopic Characterization of GRP1 PH Domain Interaction with Anionic Membranes

The pleckstrin homology (PH) domain of general receptor for phosphoionositides 1 (GRP1-PHD) binds specifically to phosphatidylinositol (3,4,5)-triphosphate (PIP3 ), and acts as a second messenger. Using an extensive array of molecular dynamics (MD) simulations employing highly mobile membrane mimetic (HMMM) model as well as complementary full membrane simulations, we capture differentiable binding and dynamics of GRP1-PHD in the presence of membranes containing PC, PS, and PIP3 lipids in varying compositions. While GRP1-PHD forms only transient interactions with pure PC membranes, incorporation of anionic lipids resulted in stable membrane-bound configurations. We report the first observation of two distinct PIP3 binding modes on GRP1-PHD, involving PIP3 interactions at a "canonical" and at an "alternate" site, suggesting the possibility of simultaneous binding of multiple anionic lipids. The full membrane simulations confirmed the stability of the membrane bound pose of GRP1-PHD as captured from our HMMM membrane binding simulations. By performing additional steered membrane unbinding simulations and calculating nonequilibrium work associated with the process, as well as metadynamics simulations, on the protein bound to full membranes, allowing for more quantitative examination of the binding strength of the GRP1-PHD to the membrane, we demonstrate that along with the bound PIP3 , surrounding anionic PS lipids increase the energetic cost of unbinding of GRP1-PHD from the canonical mode, causing them to dissociate more slowly than the alternate mode. Our results demonstrate that concurrent binding of multiple anionic lipids by GRP1-PHD contributes to its membrane affinity, which in turn control its signaling activity. © 2019 Wiley Periodicals, Inc.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Coaction of Electrostatic and Hydrophobic Interactions: Dynamic Constraints on Disordered TrkA Juxtamembrane Domain

In the receptor tyrosine kinase family, conformational change induced by ligand binding is transmitted across the membrane via a single transmembrane helix and a flexible juxtamembrane domain (JMD). Membrane dynamics makes it challenging to study the structural mechanism of receptor activation experimentally. In this study, we employ all-atom molecular dynamics with highly mobile membrane mimetic (HMMM) to capture the native conformation of the JMD in tropomyosin receptor kinase A (TrkA). We find that phosphatidylinositol 4,5-bisphosphate (PIP₂) lipids engage in stable binding with multiple basic residues. Anionic lipids can compete with salt bridges within the peptide and alter TrkA-JMD conformation. We discover three-residue insertion into the membrane and are able to either enhance or reduce the level of insertion through computationally-designed point mutations. The vesicle-binding experiment supports computational results and indicates that hydrophobic insertion is comparable to electrostatic binding for membrane anchoring. Biochemical assays on cell lines with mutated TrkA show that enhanced TrkA-JMD insertion promotes receptor degradation but does not affect the short-term signaling capacity. Our joint work points to a scenario where lipid headgroups and tails interact with basic and hydrophobic residues on disordered domain, respectively, to restrain flexibility and potentially modulate protein function.

A network of phosphatidylinositol 4,5-bisphosphate binding sites regulates gating of the Ca2+-activated Cl- channel ANO1 (TMEM16A)

ANO1 (TMEM16A) is a Ca²⁺-activated Cl^- channel that regulates diverse cellular functions including fluid secretion, neuronal excitability, and smooth muscle contraction. ANO1 is activated by elevation of cytosolic Ca²⁺ and modulated by phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂]. Here, we describe a closely concerted experimental and computational study, including electrophysiology, mutagenesis, functional assays, and extended sampling of lipid-protein interactions with molecular dynamics (MD) to characterize PI(4,5)P₂ binding modes and sites on ANO1. ANO1 currents in excised inside-out patches activated by 270 nM Ca²⁺ at +100 mV are increased by exogenous PI(4,5)P₂ with an EC₅₀ = 1.24 µM. The effect of PI(4,5)P₂ is dependent on membrane voltage and Ca²⁺ and is explained by a stabilization of the ANO1 Ca²⁺-bound open state. Unbiased atomistic MD simulations with 1.4 mol% PI(4,5)P₂ in a phosphatidylcholine bilayer identified 8 binding sites with significant probability of binding PI(4,5)P₂ Three of these sites captured 85% of all ANO1-PI(4,5)P₂ interactions. Mutagenesis of basic amino acids near the membrane-cytosol interface found 3 regions of ANO1 critical for PI(4,5)P₂ regulation that correspond to the same 3 sites identified by MD. PI(4,5)P₂ is stabilized by hydrogen bonding between amino acid side chains and phosphate/hydroxyl groups on PI(4,5)P₂ Binding of PI(4,5)P₂ alters the position of the cytoplasmic extension of TM6, which plays a crucial role in ANO1 channel gating, and increases the accessibility of the inner vestibule to Cl^- ions. We propose a model consisting of a network of 3 PI(4,5)P₂ binding sites at the cytoplasmic face of the membrane allosterically regulating ANO1 channel gating.