Efficient Exploration of Membrane-Associated Phenomena at Atomic Resolution

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.

Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

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.

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.

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.

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.

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.

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.

Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane

Infection of human cells by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) relies on its binding to a specific receptor and subsequent fusion of the viral and host cell membranes. The fusion peptide (FP), a short peptide segment in the spike protein, plays a central role in the initial penetration of the virus into the host cell membrane, followed by the fusion of the two membranes. Here, we use an array of molecular dynamics simulations that take advantage of the highly mobile membrane mimetic model to investigate the interaction of the SARS-CoV2 FP with a lipid bilayer representing mammalian cellular membranes at an atomic level and to characterize the membrane-bound form of the peptide. Six independent systems were generated by changing the initial positioning and orientation of the FP with respect to the membrane, and each system was simulated in five independent replicas, each for 300 ns. In 73% of the simulations, the FP reaches a stable, membrane-bound configuration, in which the peptide deeply penetrated into the membrane. Clustering of the results reveals three major membrane-binding modes (binding modes 1-3), in which binding mode 1 populates over half of the data points. Taking into account the sequence conservation among the viral FPs and the results of mutagenesis studies establishing the role of specific residues in the helical portion of the FP in membrane association, the significant depth of penetration of the whole peptide, and the dense population of the respective cluster, we propose that the most deeply inserted membrane-bound form (binding mode 1) represents more closely the biologically relevant form. Analysis of FP-lipid interactions shows the involvement of specific residues, previously described as the "fusion-active core residues," in membrane binding. Taken together, the results shed light on a key step involved in SARS-CoV2 infection, with potential implications in designing novel inhibitors.

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.

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.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Calcium-Induced Lipid Nanocluster Structures: Sculpturing of the Plasma Membrane

The plasma membrane of the cell is a complex, tightly regulated, heterogeneous environment shaped by proteins, lipids, and small molecules. Ca²⁺ ions are important cellular messengers, spatially separated from anionic lipids. After cell injury, disease, or apoptotic events, anionic lipids are externalized to the outer leaflet of the plasma membrane and encounter Ca²⁺, resulting in dramatic changes in the plasma membrane structure and initiation of signaling cascades. Despite the high chemical and biological significance, the structures of lipid-Ca²⁺ nanoclusters are still not known. Previously, we demonstrated by solid-state nuclear magnetic resonance (NMR) spectroscopy that upon binding to Ca²⁺, individual phosphatidylserine lipids populate two distinct yet-to-be-characterized structural environments. Here, we concurrently employ extensive all-atom molecular dynamics (MD) simulations with our accelerated membrane mimetic and detailed NMR measurements to identify lipid-Ca²⁺ nanocluster conformations. We find that major structural characteristics of these nanoclusters, including interlipid pair distances and chemical shifts, agree with observable NMR parameters. Simulations reveal that lipid-ion nanoclusters are shaped by two characteristic, long-lived lipid structures induced by divalent Ca²⁺. Using ab initio quantum mechanical calculations of chemical shifts on MD-captured lipid-ion complexes, we show that computationally observed conformations are validated by experimental NMR data. Both NMR measurements of diluted specifically labeled lipids and MD simulations reveal that the basic structural unit that reshapes the membrane is a Ca²⁺-coordinated phosphatidylserine tetramer. Our combined computational and experimental approach presented here can be applied to other complex systems in which charged membrane-active molecular agents leave structural signatures on lipids.

Molecular details of spontaneous insertion and interaction of HCV non-structure 3 protease protein domain with PIP2-containing membrane

Hepatitis C virus (HCV), known as the leading cause of liver cirrhosis, viral hepatitis, and hepatocellular carcinoma, has been affecting more than 150 million people globally. The HCV non-structure 3 (NS3) protease protein domain plays a key role in HCV replication and pathogenesis; and is currently a primary target for HCV antiviral therapy. Through unbiased molecular dynamics simulations which take advantage of the novel highly mobile membrane mimetic model, we constructed the membrane-bound state of the protein domain at the atomic level. Our results indicated that protease domain of HCV NS3 protein can spontaneously bind and penetrate to an endoplasmic reticulum complex membrane containing phosphatidylinositol 4,5-bisphosphate (PIP2). An amphipathic helix α₀ and loop S1 show their anchoring role to keep the protein on the membrane surface. Proper orientation of the protein domain at membrane surface was identified through measuring tilt angles of two specific vectors, wherein residue R161 plays a crucial role in its final orientation. Remarkably, PIP2 molecules were observed to bind to three main sites of the protease domain via specific electrostatic contacts and hydrogen bonds. PIP2-interaction determines the protein orientation at the membrane while both hydrophobic interplay and PIP2-interaction can stabilize the NS3 - membrane complex. Simulated results provide us with a detailed characterization of insertion, orientation and PIP2-interaction of NS3 protease domain at membrane environment, thus enhancing our understanding of structural functions and mechanism for the association of HCV non-structure 3 protein with respect to ER membranes.

Coupling X-Ray Reflectivity and In Silico Binding to Yield Dynamics of Membrane Recognition by Tim1

The dynamic nature of lipid membranes presents significant challenges with respect to understanding the molecular basis of protein/membrane interactions. Consequently, there is relatively little known about the structural mechanisms by which membrane-binding proteins might distinguish subtle variations in lipid membrane composition and/or structure. We have previously developed a multidisciplinary approach that combines molecular dynamics simulation with interfacial x-ray scattering experiments to produce an atomistic model for phosphatidylserine recognition by the immune receptor Tim4. However, this approach requires a previously determined protein crystal structure in a membrane-bound conformation. Tim1, a Tim4 homolog with distinct differences in both immunological function and sensitivity to membrane composition, was crystalized in a closed-loop conformation that is unlikely to support membrane binding. Here we have used a previously described highly mobile membrane mimetic membrane in combination with a conventional lipid bilayer model to generate a membrane-bound configuration of Tim1 in silico. This refined structure provided a significantly improved fit of experimental x-ray reflectivity data. Moreover, the coupling of the x-ray reflectivity analysis with both highly mobile membrane mimetic membranes and conventional lipid bilayer molecular dynamics simulations yielded a dynamic model of phosphatidylserine membrane recognition by Tim1 with atomic-level detail. In addition to providing, to our knowledge, new insights into the molecular mechanisms that distinguish the various Tim receptors, these results demonstrate that in silico membrane-binding simulations can remove the requirement that the existing crystal structure be in the membrane-bound conformation for effective x-ray reflectivity analysis. Consequently, this refined methodology has the potential for much broader applicability with respect to defining the atomistic details of membrane-binding proteins.

Lipid specificity of the membrane binding domain of coagulation factor X

Essentials Membrane-binding GLA domains of coagulation factors are essential for proper clot formation. Factor X (FX) is specific to phosphatidylserine (PS) lipids through unknown atomic-level interactions. Molecular dynamics simulations were used to develop the first membrane-bound model of FX-GLA. PS binding modes of FX-GLA were described, and potential PS-specific binding sites identified.

SUMMARY

Background Factor X (FX) binds to cell membranes in a highly phospholipid-dependent manner and, in complex with tissue factor and factor VIIa (FVIIa), initiates the clotting cascade. Experimental information concerning the membrane-bound structure of FX with atomic resolution has remained elusive because of the fluid nature of cellular membranes. FX is known to bind preferentially to phosphatidylserine (PS). Objectives To develop the first membrane-bound model of the FX-GLA domain to PS at atomic level, and to identify PS-specific binding sites of the FX-GLA domain. Methods Molecular dynamics (MD) simulations were performed to develop an atomic-level model for the FX-GLA domain bound to PS bilayers. We utilized a membrane representation with enhanced lipid mobility, termed the highly mobile membrane mimetic (HMMM), permitting spontaneous membrane binding and insertion by FX-GLA in multiple 100-ns simulations. In 14 independent simulations, FX-GLA bound spontaneously to the membrane. The resulting membrane-bound models were converted from HMMM to conventional membrane and simulated for an additional 100 ns. Results The final membrane-bound FX-GLA model allowed for detailed characterization of the orientation, insertion depth and lipid interactions of the domain, providing insight into the molecular basis of its PS specificity. All binding simulations converged to the same configuration despite differing initial orientations. Conclusions Analysis of interactions between residues in FX-GLA and lipid-charged groups allowed for potential PS-specific binding sites to be identified. This new structural and dynamic information provides an additional step towards a full understanding of the role of atomic-level lipid-protein interactions in regulating the critical and complex clotting cascade.

Extension of the Highly Mobile Membrane Mimetic to Transmembrane Systems through Customized in Silico Solvents

The mechanics of the protein-lipid interactions of transmembrane proteins are difficult to capture with conventional atomic molecular dynamics, due to the slow lateral diffusion of lipids restricting sampling to states near the initial membrane configuration. The highly mobile membrane mimetic (HMMM) model accelerates lipid dynamics by modeling the acyl tails nearest to the membrane center as a fluid organic solvent while maintaining an atomic description of the lipid headgroups and short acyl tails. The HMMM has been applied to many peripheral protein systems; however, the organic solvent used to date caused deformations in transmembrane proteins by intercalating into the protein and disrupting interactions between individual side chains. We ameliorate the effect of the solvent on transmembrane protein structure through the development of two new in silico Lennard-Jones solvents. The parameters for the new solvents were determined through an extensive parameter search in order to match the bulk properties of alkanes in a highly simplified model. Using these new solvents, we substantially improve the insertion free energy profiles of 10 protein side chain analogues across the entire bilayer. In addition, we reduce the intercalation of solvent into transmembrane systems, resulting in native-like transmembrane protein structures from five different topological classes within a HMMM bilayer. The parametrization of the solvents, in addition to their computed physical properties, is discussed. By combining high lipid lateral diffusion with intact transmembrane proteins, we foresee the developed solvents being useful to efficiently identify membrane composition inhomogeneities and lipid binding caused by the presence of membrane proteins.

Nanodiscs in Membrane Biochemistry and Biophysics

Membrane proteins play a most important part in metabolism, signaling, cell motility, transport, development, and many other biochemical and biophysical processes which constitute fundamentals of life on the molecular level. Detailed understanding of these processes is necessary for the progress of life sciences and biomedical applications. Nanodiscs provide a new and powerful tool for a broad spectrum of biochemical and biophysical studies of membrane proteins and are commonly acknowledged as an optimal membrane mimetic system that provides control over size, composition, and specific functional modifications on the nanometer scale. In this review we attempted to combine a comprehensive list of various applications of nanodisc technology with systematic analysis of the most attractive features of this system and advantages provided by nanodiscs for structural and mechanistic studies of membrane proteins.

Differential Membrane Binding Mechanics of Synaptotagmin Isoforms Observed in Atomic Detail

Synaptotagmin (Syt) is a membrane-associated protein involved in vesicle fusion through the SNARE complex that is found throughout the human body in 17 different isoforms. These isoforms have two membrane-binding C2 domains, which sense Ca²⁺ and thereby promote anionic membrane binding and lead to vesicle fusion. Through molecular dynamics simulations using the highly mobile membrane mimetic acclerated bilayer model, we have investigated how small protein sequence changes in the Ca²⁺-binding loops of the C2 domains may give rise to the experimentally determined difference in binding kinetics between Syt-1 and Syt-7 isoforms. Syt-7 C2 domains are found to form more close contacts with anionic phospholipid headgroups, particularly in loop 1, where an additional positive charge in Syt-7 draws the loop closer to the membrane and causes the anchoring residue F167 to insert deeper into the bilayer than the corresponding methionine in Syt-1 (M173). By performing additional replica exchange umbrella sampling calculations, we demonstrate that these additional contacts increase the energetic cost of unbinding the Syt-7 C2 domains from the bilayer, causing them to unbind more slowly than their counterparts in Syt-1.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

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

Slow diffusion of the lipids in conventional all-atom simulations of membrane systems makes it difficult to sample large rearrangements of lipids and protein-lipid interactions. Recently, Tajkhorshid and co-workers developed the highly mobile membrane-mimetic (HMMM) model with accelerated lipid motion by replacing the lipid tails with small organic molecules. The HMMM model provides accelerated lipid diffusion by one to two orders of magnitude, and is particularly useful in studying membrane-protein associations. However, building an HMMM simulation system is not easy, as it requires sophisticated treatment of the lipid tails. In this study, we have developed CHARMM-GUI HMMM Builder (http://www.charmm-gui.org/input/hmmm) to provide users with ready-to-go input files for simulating HMMM membrane systems with/without proteins. Various lipid-only and protein-lipid systems are simulated to validate the qualities of the systems generated by HMMM Builder with focus on the basic properties and advantages of the HMMM model. HMMM Builder supports all lipid types available in CHARMM-GUI and also provides a module to convert back and forth between an HMMM membrane and a full-length membrane. We expect HMMM Builder to be a useful tool in studying membrane systems with enhanced lipid diffusion.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

Interactions of Pleckstrin Homology Domains with Membranes: Adding Back the Bilayer via High-Throughput Molecular Dynamics

A molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Atomic-level description of protein-lipid interactions using an accelerated membrane model

Peripheral membrane proteins are structurally diverse proteins that are involved in fundamental cellular processes. Their activity of these proteins is frequently modulated through their interaction with cellular membranes, and as a result techniques to study the interfacial interaction between peripheral proteins and the membrane are in high demand. Due to the fluid nature of the membrane and the reversibility of protein-membrane interactions, the experimental study of these systems remains a challenging task. Molecular dynamics simulations offer a suitable approach to study protein-lipid interactions; however, the slow dynamics of the lipids often prevents sufficient sampling of specific membrane-protein interactions in atomistic simulations. To increase lipid dynamics while preserving the atomistic detail of protein-lipid interactions, in the highly mobile membrane-mimetic (HMMM) model the membrane core is replaced by an organic solvent, while short-tailed lipids provide a nearly complete representation of natural lipids at the organic solvent/water interface. Here, we present a brief introduction and a summary of recent applications of the HMMM to study different membrane proteins, complementing the experimental characterization of the presented systems, and we offer a perspective of future applications of the HMMM to study other classes of membrane proteins. This article is part of a Special Issue entitled: Membrane proteins edited by J.C. Gumbart and Sergei Noskov.