Protein-lipid interactions and protein anchoring modulate the modes of association of the globular domain of the Prion protein and Doppel protein to model membrane patches

The Prion protein is the molecular hallmark of the incurable prion diseases affecting mammals, including humans. The protein-only hypothesis states that the misfolding, accumulation, and deposition of the Prion protein play a critical role in toxicity. The cellular Prion protein (PrPC) anchors to the extracellular leaflet of the plasma membrane and prefers cholesterol- and sphingomyelin-rich membrane domains. Conformational Prion protein conversion into the pathological isoform happens on the cell surface. In vitro and in vivo experiments indicate that Prion protein misfolding, aggregation, and toxicity are sensitive to the lipid composition of plasma membranes and vesicles. A picture of the underlying biophysical driving forces that explain the effect of Prion protein - lipid interactions in physiological conditions is needed to develop a structural model of Prion protein conformational conversion. To this end, we use molecular dynamics simulations that mimic the interactions between the globular domain of PrPC anchored to model membrane patches. In addition, we also simulate the Doppel protein anchored to such membrane patches. The Doppel protein is the closest in the phylogenetic tree to PrPC, localizes in an extracellular milieu similar to that of PrPC, and exhibits a similar topology to PrPC even if the amino acid sequence is only 25% identical. Our simulations show that specific protein-lipid interactions and conformational constraints imposed by GPI anchoring together favor specific binding sites in globular PrPC but not in Doppel. Interestingly, the binding sites we found in PrPC correspond to prion protein loops, which are critical in aggregation and prion disease transmission barrier (β2-α2 loop) and in initial spontaneous misfolding (α2-α3 loop). We also found that the membrane re-arranges locally to accommodate protein residues inserted in the membrane surface as a response to protein binding.

Improved Protein Model in SPICA Force Field

The previous version of the SPICA coarse-grained (CG) force field (FF) protein model focused primarily on membrane proteins and successfully reproduced the dimerization free energies of several transmembrane helices and the stable structures of various membrane protein assemblies. However, that model had limited accuracy when applied to other proteins, such as intrinsically disordered proteins (IDPs) and peripheral proteins, because the dimensions of the IDPs in an aqueous solution were too compact, and protein binding on the lipid membrane surface was overstabilized. To improve the accuracy of the SPICA FF model for the simulation of such systems, in this study, we introduce protein secondary structure-dependent nonbonded interaction parameters to the backbone segments and reoptimize almost all nonbonded parameters for amino acids. The improved FF proposed here successfully reproduces the radii of gyration of various IDPs, the binding sensitivity of several peripheral membrane proteins, and the dimerization free energies of several transmembrane helices. The new model also shows improved agreement with experiments on the free energy of peptide association in water. In addition, an extensive library of nonbonded interactions between proteins and lipids, including various glycerophospholipids, sphingolipids, and cholesterol, allows the study of specific interactions between lipids and peripheral and transmembrane proteins. Hence, the new SPICA FF (version 2) proposed herein is applicable with high accuracy for simulating a wide range of protein systems.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Computational modeling of membrane trafficking processes: From large molecular assemblies to chemical specificity

In the last decade, molecular dynamics (MD) simulations have become an essential tool to investigate the molecular properties of membrane trafficking processes, often in conjunction with experimental approaches. The combination of MD simulations with recent developments in structural biology, such as cryo-electron microscopy and artificial intelligence-based structure determination, opens new, exciting possibilities for future investigations. However, the full potential of MD simulations to provide a molecular view of the complex and dynamic processes involving membrane trafficking can only be realized if certain limitations are addressed, and especially those concerning the quality of coarse-grain models, which, despite recent successes in describing large-scale systems, still suffer from far-from-ideal chemical accuracy. In this review, we will highlight recent success stories of MD simulations in the investigation of membrane trafficking processes, their implications for future research, and the challenges that lie ahead in this specific research domain.

Insights into membrane association of the SMP domain of extended synaptotagmin

The Synaptotagmin-like Mitochondrial-lipid-binding Protein (SMP) domain is a newly identified lipid transfer module present in proteins that regulate lipid homeostasis at membrane contact sites (MCSs). However, how the SMP domain associates with the membrane to extract and unload lipids is unclear. Here, we performed in vitro DNA brick-assisted lipid transfer assays and in silico molecular dynamics simulations to investigate the molecular basis of the membrane association by the SMP domain of extended synaptotagmin (E-Syt), which tethers the tubular endoplasmic reticulum (ER) to the plasma membrane (PM). We demonstrate that the SMP domain uses its tip region to recognize the extremely curved subdomain of tubular ER and the acidic-lipid-enriched PM for highly efficient lipid transfer. Supporting these findings, disruption of these mechanisms results in a defect in autophagosome biogenesis contributed by E-Syt. Our results suggest a model that provides a coherent picture of the action of the SMP domain at MCSs.

NS1 from Two Zika Virus Strains Differently Interact with a Membrane: Insights to Understand Their Differential Virulence

Zika virus (ZIKV) from Uganda (UG) expresses a phenotype related to fetal loss, whereas the variant from Brazil (BR) induces microcephaly in neonates. The differential virulence has a direct relation to biomolecular mechanisms that make one strain more aggressive than the other. The nonstructural protein 1 (NS1) is a key viral toxin to comprehend these viral discrepancies because of its versatility in many processes of the virus life cycle. Here, we aim to examine through coarse-grained models and molecular dynamics simulations the protein-membrane interactions for both NS1_ZIKV-UG and NS1_ZIKV-BR dimers. A first evaluation allowed us to establish that the NS1 proteins, in the membrane presence, explore new conformational spaces when compared to systems simulated without a lipid bilayer. These events derive from both differential coupling patterns and discrepant binding affinities to the membrane. The N-terminal domain, intertwined loop, and greasy finger proposed previously as binding membrane regions were also computationally confirmed by us. The anchoring sites have aromatic and ionizable residues that manage the assembly of NS1 toward the membrane, especially for the Ugandan variant. Furthermore, in the presence of the membrane, the difference in the dynamic cross-correlation of residues between the two strains is particularly pronounced in the putative epitope regions. This suggests that the protein-membrane interaction induces changes in the distal part related to putative epitopes. Taken together, these results open up new strategies for the treatment of flaviviruses that would specifically target these dynamic differences.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

Efficient Quantification of Lipid Packing Defect Sensing by Amphipathic Peptides: Comparing Martini 2 and 3 with CHARMM36

In biological systems, proteins can be attracted to curved or stretched regions of lipid bilayers by sensing hydrophobic defects in the lipid packing on the membrane surface. Here, we present an efficient end-state free energy calculation method to quantify such sensing in molecular dynamics simulations. We illustrate that lipid packing defect sensing can be defined as the difference in mechanical work required to stretch a membrane with and without a peptide bound to the surface. We also demonstrate that a peptide's ability to concurrently induce excess leaflet area (tension) and elastic softening─a property we call the "characteristic area of sensing" (CHAOS)─and lipid packing sensing behavior are in fact two sides of the same coin. In essence, defect sensing displays a peptide's propensity to generate tension. The here-proposed mechanical pathway is equally accurate yet, computationally, about 40 times less costly than the commonly used alchemical pathway (thermodynamic integration), allowing for more feasible free energy calculations in atomistic simulations. This enabled us to directly compare the Martini 2 and 3 coarse-grained and the CHARMM36 atomistic force fields in terms of relative binding free energies for six representative peptides including the curvature sensor ALPS and two antiviral amphipathic helices (AH). We observed that Martini 3 qualitatively reproduces experimental trends while producing substantially lower (relative) binding free energies and shallower membrane insertion depths compared to atomistic simulations. In contrast, Martini 2 tends to overestimate (relative) binding free energies. Finally, we offer a glimpse into how our end-state-based free energy method can enable the inverse design of optimal lipid packing defect sensing peptides when used in conjunction with our recently developed evolutionary molecular dynamics (Evo-MD) method. We argue that these optimized defect sensors─aside from their biomedical and biophysical relevance─can provide valuable targets for the development of lipid force fields.

SEED LIPID DROPLET PROTEIN1, SEED LIPID DROPLET PROTEIN2, and LIPID DROPLET PLASMA MEMBRANE ADAPTOR mediate lipid droplet-plasma membrane tethering

Membrane contact sites (MCSs) are interorganellar connections that allow for the direct exchange of molecules, such as lipids or Ca2+ between organelles, but can also serve to tether organelles at specific locations within cells. Here, we identified and characterized three proteins of Arabidopsis thaliana that form a lipid droplet (LD)-plasma membrane (PM) tethering complex in plant cells, namely LD-localized SEED LD PROTEIN (SLDP) 1 and SLDP2 and PM-localized LD-PLASMA MEMBRANE ADAPTOR (LIPA). Using proteomics and different protein-protein interaction assays, we show that both SLDPs associate with LIPA. Disruption of either SLDP1 and SLDP2 expression, or that of LIPA, leads to an aberrant clustering of LDs in Arabidopsis seedlings. Ectopic co-expression of one of the SLDPs with LIPA is sufficient to reconstitute LD-PM tethering in Nicotiana tabacum pollen tubes, a cell type characterized by dynamically moving LDs in the cytosolic streaming. Furthermore, confocal laser scanning microscopy revealed both SLDP2.1 and LIPA to be enriched at LD-PM contact sites in seedlings. These and other results suggest that SLDP and LIPA interact to form a tethering complex that anchors a subset of LDs to the PM during post-germinative seedling growth in Arabidopsis.

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.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

SPICA Force Field for Proteins and Peptides

A coarse-grained (CG) model for peptides and proteins was developed as an extension of the Surface Property fItting Coarse grAined (SPICA) force field (FF). The model was designed to examine membrane proteins that are fully compatible with the lipid membranes of the SPICA FF. A preliminary version of this protein model was created using thermodynamic properties, including the surface tension and density in the SPICA (formerly called SDK) FF. In this study, we improved the CG protein model to facilitate molecular dynamics (MD) simulations with a reproduction of multiple properties from both experiments and all-atom (AA) simulations. An elastic network model was adopted to maintain the secondary structure within a single chain. The side-chain analogues reproduced the transfer free energy profiles across the lipid membrane and demonstrated reasonable association free energy (potential of mean force) in water compared to those from AA MD. A series of peptides/proteins adsorbed onto or penetrated into the membrane simulated by the CG MD correctly predicted the penetration depths and tilt angles of peripheral and transmembrane peptides/proteins as comparable to those in the orientations of proteins in membranes (OPM) database. In addition, the dimerization free energies of several transmembrane helices within a lipid bilayer were comparable to those from experimental estimation. Application studies on a series of membrane protein assemblies, scramblases, and poliovirus capsids demonstrated the good performance of the SPICA FF.

A conserved function of Human DLC3 and Drosophila Cv-c in testis development

The identification of genes affecting gonad development is essential to understand the mechanisms causing Variations/Differences in Sex Development (DSD). Recently, a DLC3 mutation was associated with male gonadal dysgenesis in 46,XY DSD patients. We have studied the requirement of Cv-c, the Drosophila ortholog of DLC3, in Drosophila gonad development, as well as the functional capacity of DLC3 human variants to rescue cv-c gonad defects. We show that Cv-c is required to maintain testis integrity during fly development. We find that Cv-c and human DLC3 can perform the same function in fly embryos, as flies carrying wild type but not patient DLC3 variations can rescue gonadal dysgenesis, suggesting functional conservation. We also demonstrate that the StART domain mediates Cv-c's function in the male gonad independently from the GAP domain's activity. This work demonstrates a role for DLC3/Cv-c in male gonadogenesis and highlights a novel StART domain mediated function required to organize the gonadal mesoderm and maintain its interaction with the germ cells during testis development.

Synaptotagmin-1 C2B domains cooperatively stabilize the fusion stalk via a master-servant mechanism

Synaptotagmin-1 is a low-affinity Ca2+ sensor that triggers synchronous vesicle fusion. It contains two similar C2 domains (C2A and C2B) that cooperate in membrane binding, being the C2B domain the...

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.