SNARE complex alters the interactions of the Ca2+ sensor synaptotagmin 1 with lipid bilayers

Release of neuronal transmitters from nerve terminals is triggered by the molecular Ca²⁺ sensor synaptotagmin 1 (Syt1). Syt1 is a transmembrane protein attached to the synaptic vesicle (SV), and its cytosolic region comprises two domains, C2A and C2B, which are thought to penetrate into lipid bilayers upon Ca²⁺ binding. Before fusion, SVs become attached to the presynaptic membrane (PM) by the four-helical SNARE complex, which is thought to bind the C2B domain in vivo. To understand how the interactions of Syt1 with lipid bilayers and the SNARE complex trigger fusion, we performed molecular dynamics (MD) simulations at a microsecond scale. We investigated how the isolated C2 modules and the C2AB tandem of Syt1 interact with membranes mimicking either SV or PM. The simulations showed that the C2AB tandem can either bridge SV and PM or insert into PM with its Ca²⁺-bound tips and that the latter configuration is more favorable. Surprisingly, C2 domains did not cooperate in penetrating into PM but instead mutually hindered their insertion into the bilayer. To test whether the interaction of Syt1 with lipid bilayers could be affected by the C2B-SNARE attachment, we performed systematic conformational analysis of the C2AB-SNARE complex. Notably, we found that the C2B-SNARE interface precludes the coupling of C2 domains and promotes their insertion into PM. We performed the MD simulations of the prefusion protein complex positioned between the lipid bilayers mimicking PM and SV, and our results demonstrated in silico that the presence of the Ca²⁺ bound C2AB tandem promotes lipid merging. Altogether, our MD simulations elucidated the role of the Syt1-SNARE interactions in the fusion process and produced the dynamic all-atom model of the prefusion protein-lipid complex.     

Computational study of colipase interaction with lipid droplets and bile salt micelles

Colipase is a key element in the lipase-catalyzed hydrolysis of dietary lipids. Although devoid of enzymatic activity, colipase promotes the pancreatic lipase activity in physiological intestinal conditions by anchoring the enzyme at the surface of lipid droplets. Analysis of structures of NMR colipase models and simulations of their interactions with various lipid aggregates, lipid droplet, and bile salt micelle, were carried out to determine and to map the lipid binding sites on colipase. We show that the micelle and the oil droplet bind to the same side of colipase 3D structure, mainly the hydrophobic fingers. Moreover, it appears that, although colipase has a single direction of interaction with a lipid interface, it does not bind in a specific way but rather oscillates between different positions. Indeed, different NMR models of colipase insert different fragments of sequence in the interface, either simultaneously or independently. This supports the idea that colipase finger plasticity may be crucial to adapt the lipase activity to different lipid aggregates.     

Computer modeling of the membrane interaction of FYVE domains

FYVE domains are membrane targeting domains that are found in proteins involved in endosomal trafficking and signal transduction pathways. Most FYVE domains bind specifically to phosphatidylinositol 3-phosphate (PI(3)P), a lipid that resides mainly in endosomal membranes. Though the specific interactions between FYVE domains and the headgroup of PI(3)P have been well characterized, principally through structural studies, the available experimental structures suggest several different models for FYVE/membrane association. Thus, the manner in which FYVE domains adsorb to the membrane surface remains to be elucidated. Towards this end, recent experiments have shown that FYVE domains bind PI(3)P in the context of phospholipid bilayers and that hydrophobic residues on a conserved loop are able to penetrate the membrane interface in a PI(3)P-dependent manner.Here, the finite difference Poisson-Boltzmann (FDPB) method has been used to calculate the energetic interactions of FYVE domains with phospholipid membranes. Based on the computational analysis, it is found that (1) recruitment to membranes is facilitated by non-specific electrostatic interactions that occur between basic residues on the domains and acidic phospholipids in the membrane, (2) the energetic analysis can quantitatively differentiate among the modes of membrane association proposed by the experimentally determined structures, (3) FDPB calculations predict energetically feasible models for the membrane-associated states of FYVE domains, (4) these models are consistent with the observation that conserved hydrophobic residues insert into the membrane interface, and (5) the calculations provide a molecular model for the hydrophobic partitioning: binding of PI(3)P significantly neutralizes positive potential in the region of the hydrophobic residues, which acts as an "electrostatic switch" by reducing the energetic barrier for membrane penetration. Finally, the computational results are extended to FYVE domains of unknown structure through the construction of high quality homology models for human FYVE sequences.     

Hydrophobic Compounds Reshape Membrane Domains

Cell membranes have a complex lateral organization featuring domains with distinct composition, also known as rafts, which play an essential role in cellular processes such as signal transduction and protein trafficking. In vivo, perturbations of membrane domains (e.g., by drugs or lipophilic compounds) have major effects on the activity of raft-associated proteins and on signaling pathways, but they are difficult to characterize because of the small size of the domains, typically below optical resolution. Model membranes, instead, can show macroscopic phase separation between liquid-ordered and liquid-disordered domains, and they are often used to investigate the driving forces of membrane lateral organization. Studies in model membranes have shown that some lipophilic compounds perturb membrane domains, but it is not clear which chemical and physical properties determine domain perturbation. The mechanisms of domain stabilization and destabilization are also unknown. Here we describe the effect of six simple hydrophobic compounds on the lateral organization of phase-separated model membranes consisting of saturated and unsaturated phospholipids and cholesterol. Using molecular simulations, we identify two groups of molecules with distinct behavior: aliphatic compounds promote lipid mixing by distributing at the interface between liquid-ordered and liquid-disordered domains; aromatic compounds, instead, stabilize phase separation by partitioning into liquid-disordered domains and excluding cholesterol from the disordered domains. We predict that relatively small concentrations of hydrophobic species can have a broad impact on domain stability in model systems, which suggests possible mechanisms of action for hydrophobic compounds in vivo.     

Achieving Peptide Binding Specificity and Promiscuity by Loops: Case of the Forkhead-Associated Domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.     

Phase separation and bistability in a three-dimensional model for protein domain formation at biomembranes

Proteins in living cells interact with membranes. They may bind to or unbind from the membrane to the cytosol depending on the lipid composition of the membrane and their interaction with cytosolic enzymes. Moreover, proteins can accumulate at the membrane and assemble in spatial domains. Here, a simple model of protein cycling at biomembranes is studied, when the total number of proteins is conserved. Specifically, we consider the spatio-temporal dynamics of MARCKS proteins and their interactions with enzymes facilitating translocation from and rebinding to the membrane. The model exhibits two qualitatively different mechanisms of protein domain formation: phase separation related to a long-wave instability of a membrane state with homogeneous protein coverage and stable coexistence of two states with different homogeneous protein coverage in bistable media. We evaluate the impact of the cytosolic volume on the occurrence of protein pattern formation by simulations in a three-dimensional model. We show that the explicit treatment of the volume in the model leads to an effective rescaling of the reaction rates. For a simplified model of protein cycling, we can derive analytical expressions for the rescaling coefficients and verify them by direct simulations with the complete three-dimensional model.     

BAR domains and membrane curvature: bringing your curves to the BAR

BAR (bin, amphiphysin and Rvs161/167) domains are a unique class of dimerization domains, whose dimerization interface is edged by a membrane-binding surface. In its dimeric form, the membrane-binding interface is concave, and this gives the ability to bind better to curved membranes, i.e. to sense membrane curvature. When present at higher concentrations, the domain can stabilize membrane curvature, generating lipid tubules. This domain is found in many contexts in a wide variety of proteins, where the dimerization and membrane-binding function of this domain is likely to have a profound effect on protein activity. If these proteins function as predicted, then there will be membrane subdomains based on curvature, and thus there is an additional layer of compartmentalization on membranes. These and other possible functions of the BAR domain are discussed.     

Lipid asymmetry in DLPC/DSPC-supported lipid bilayers: a combined AFM and fluorescence microscopy study

A fundamental attribute of cell membranes is transmembrane asymmetry, specifically the formation of ordered phase domains in one leaflet that are compositionally different from the opposing leaflet of the bilayer. Using model membrane systems, many previous studies have demonstrated the formation of ordered phase domains that display complete transmembrane symmetry; but there have been few reports on the more biologically relevant asymmetric membrane structures. Here we report on a combined atomic force microscopy and fluorescence microscopy study whereby we observe three different states of transmembrane symmetry in phase-separated supported lipid bilayers formed by vesicle fusion. We find that if the leaflets differ in gel-phase area fraction, then the smaller domains in one leaflet are in registry with the larger domains in the other leaflet and the system is dynamic. In a presumed lipid flip-flop process similar to Ostwald ripening, the smaller domains in one leaflet erode away whereas the large domains in the other leaflet grow until complete compositional asymmetry is reached and remains stable. We have quantified this evolution and determined that the lipid flip-flop event happens most frequently at the interface between symmetric and asymmetric DSPC domains. If both leaflets have identical area fraction of gel-phase, gel-phase domains are in registry and are static in comparison to the first state. The stability of these three DSPC domain distributions, the degree of registry observed, and the domain immobility have biological significance with regards to maintenance of lipid asymmetry in living cell membranes, communication between inner leaflet and outer leaflet, membrane adhesion, and raft mobility.     

State-dependent protein-lipid interactions of a pentameric ligand-gated ion channel in a neuronal membrane

Pentameric ligand-gated ion channels (pLGICs) are receptor proteins that are sensitive to their membrane environment, but the mechanism for how lipids modulate function under physiological conditions in a state dependent manner is not known. The glycine receptor is a pLGIC whose structure has been resolved in different functional states. Using a realistic model of a neuronal membrane coupled with coarse-grained molecular dynamics simulations, we demonstrate that some key lipid-protein interactions are dependent on the receptor state, suggesting that lipids may regulate the receptor’s conformational dynamics. Comparison with existing structural data confirms known lipid binding sites, but we also predict further protein-lipid interactions including a site at the communication interface between the extracellular and transmembrane domain. Moreover, in the active state, cholesterol can bind to the binding site of the positive allosteric modulator ivermectin. These protein-lipid interaction sites could in future be exploited for the rational design of lipid-like allosteric drugs.     

Domain formation in a fluid mixed lipid bilayer modulated through binding of the C2 protein motif

The role and mechanism of formation of lipid domains in a functional membrane have generally received limited attention. Our approach, based on the hypothesis that thermodynamic coupling between lipid-lipid and protein-lipid interactions can lead to domain formation, uses a combination of an experimental lipid bilayer model system and Monte Carlo computer simulations of a simple model of that system. The experimental system is a fluid bilayer composed of a binary mixture of phosphatidylcholine (PC) and phosphatidylserine (PS), containing 4% of a pyrene-labeled anionic phospholipid. Addition of the C2 protein motif (a structural domain found in proteins implicated in eukaryotic signal transduction and cellular trafficking processes) to the bilayer first increases and then decreases the excimer/monomer ratio of the pyrene fluorescence. We interpret this to mean that protein binding induces anionic lipid domain formation until the anionic lipid becomes saturated with protein. Monte Carlo simulations were performed on a lattice representing the lipid bilayer to which proteins were added. The important parameters are an unlike lipid-lipid interaction term and an experimentally derived preferential protein-lipid interaction term. The simulations support the experimental conclusion and indicate the existence of a maximum in PS domain size as a function of protein concentration. Thus, lipid-protein coupling is a possible mechanism for both lipid and protein clustering on a fluid bilayer. Such domains could be precursors of larger lipid-protein clusters ('rafts'), which could be important in various biological processes such as signal transduction at the level of the cell membrane.     

FSCS Reveals the Complexity of Lipid Domain Dynamics in the Plasma Membrane of Live Cells

The coexistence of lipid domains with different degrees of lipid packing in the plasma membrane of mammalian cells has been postulated, but direct evidence has so far been challenging to obtain because of the small size and short lifetime of these domains in live cells. Here, we use fluorescence spectral correlation spectroscopy in conjunction with a probe sensitive to the membrane environment to quantify spectral fluctuations associated with dynamics of membrane domains in live cells. With this method, we show that membrane domains are present in live COS-7 cells and have a lifetime lower bound of 5.90 and 14.69 ms for the ordered and disordered phases, respectively. Comparisons to simulations indicate that the underlying mechanism of these fluctuations is complex but qualitatively described by a combination of dye diffusion between membrane domains as well as the motion of domains within the membrane.     

Molecular dynamics simulations of inwardly rectifying (Kir) potassium channels: a comparative study

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.     

Molecular modelling study of antigen binding to oxazolone-specific antibodies: the Ox1 idiotypic IgG and its mature variant with increased affinity to 2-phenyloxazolone

Structural models of the variable domains of the murine anti-2-phenyloxazolone IgG (Ox1 idiotype) and its somatic variant, which has higher affinity to the hapten 2-phenyloxazolone, were constructed by computer-aided model building using known structures of highly homologous immunoglobulins as templates. Molecular dynamics simulations were used to dock the hapten between the VL and VH domains. The hapten is predicted to bind to slightly different sites in the two models. Hypotheses concerning the role of a number of preferred mutations in anti-oxazolone variants are presented. These can be tested by mutagenesis and crystallography. In particular, the higher binding affinities of the different antibody variants are shown to correlate with better complementarity of electrostatics. The molecular dynamic simulations also suggest that two mobile tryptophans at the mouth of the pocket may play an important role in the binding of hapten.     

Insights into the behavior of unsaturated diacylglycerols in mixed lipid bilayers in relation to protein kinase C activation—A molecular dynamics simulation study

The lipid second messenger diacylglycerol (DAG) is known for its involvement in many types of cellular signaling, especially as an endogenous agonist for protein kinase C (PKC). Evidence has emerged that the degree of saturation of the DAG molecules can affect PKC activation. DAG molecules with different acyl chain saturation have not only been observed to induce varying extents of PKC activation, but also to express selectivity towards different PKC isozymes. Both qualities are important for precise therapeutic activation of PKC; understanding DAG behavior at the molecular level in different environments has much potential in the development of drugs to target PKC. We used molecular dynamics simulations to study the behavior of two different unsaturated DAG species in lipid environments with varying degrees of unsaturation. We focus on phosphatidylethanolamine (PE) instead of phosphatidylcholine (PC) to more accurately model the relevant biomembranes. The effect of cholesterol (CHOL) on these systems was also explored. We found that both high level of unsaturation in the acyl chains of the DAG species and presence of CHOL in the surrounding membrane increase DAG molecule availability at the lipid–water interface. This can partially explain the previously observed differences in PKC activation strength and specificity, the complete mechanism is, however, likely to be more complex. Our simulations coupled with the current understanding of lipids highlight the need for more simulations of biologically accurate lipid environments in order to determine the correct correlations between molecular mechanisms and biological behavior when studying PKC activation.     

Position of synaptotagmin I at the membrane interface: cooperative interactions of tandem C2 domains

Synaptotagmin I is a synaptic vesicle associated membrane protein that appears to regulate Ca(2+)-mediated exocytosis. Here, the Ca(2+)-dependent membrane interactions of a water soluble fragment of synaptotagmin I (C2AB) that contains its two C2 domains (C2A and C2B) were determined using site-directed spin labeling. Membrane depth parameters were obtained for 19 spin-labeled mutants of C2AB when bound to phosphatidylcholine and phosphatidylserine membranes, and these distance constraints were used in combination with the high-resolution structures of C2A and C2B to generate a model for the membrane orientation and position of synaptotagmin at the bilayer interface. Both C2A and C2B bind to the membrane interface with their first and third Ca(2+) binding loops penetrating the membrane interface. The polybasic face of C2B does not interact with the membrane lipid but is available for electrostatic interaction with other components of the fusion machinery. When compared to positions determined previously for the isolated domains, both C2A and C2B have similar orientations; however, the two domains are positioned deeper into the bilayer interior when present in the tandem construct. These data indicate that C2A and C2B do not act independently but influence their mutual membrane penetration. This may explain the occurrence of multiple C2 domains in proteins that function in membrane trafficking and repair.     

The common phospholipid-binding activity of the N-terminal domains of PEX1 and VCP/p97

PEX1 is a type II AAA-ATPase that is indispensable for biogenesis and maintenance of the peroxisome, an organelle responsible for the primary metabolism of lipids, such as beta-oxidation and lipid biosynthesis. Recently, we demonstrated a striking structural similarity between its N-terminal domain and those of other membrane-related AAA-ATPases, such as valosine-containing protein (p97). The N-terminal domain of valosine-containing protein serves as an interface to its adaptor proteins p47 and Ufd1, whereas the physiologic interaction partner of the N-terminal domain of PEX1 remains unknown. Here we found that N-terminal domains isolated from valosine-containing protein, as well as from PEX1, bind phosphoinositides. The N-terminal domain of PEX1 appears to preferentially bind phosphatidylinositol 3-monophosphate and phosphatidylinositol 4-monophosphate, whereas the N-terminal domain of valosine-containing protein displays broad and nonspecific lipid binding. Although N-ethylmaleimide-sensitive fusion protein, CDC48 and Ufd1 have structures similar to that of valosine-containing protein, they displayed lipid specificity similar to that of the N-terminal domain of PEX1 in the assays. By mutational analysis, we demonstrate that a conserved arginine surrounded by hydrophobic residues is essential for lipid binding, despite very low sequence similarity between PEX1 and valosine-containing protein.     

Uncoupling of catalysis and colipase binding in pancreatic lipase by limited proteolysis.

In the intestine, the hydrolysis of triglycerides by pancreatic lipase is performed only in the presence of colipase, whose function is to anchor lipase to the bile-salt-coated lipid interface. Biochemical and crystallographic data on porcine and human lipases have shown that the molecule is made of two well-delimited domains. In order to get more information on the role of the domains in catalysis and colipase binding, we performed limited proteolysis on lipase from various species and obtained different patterns of cleavage. In the case of porcine and human lipases, only the C-terminal domain (12 kDa) could be obtained after chymotryptic attack, whereas in the horse enzyme the cleavage of the Leu410-Thr411 bond gave rise to a large N-terminal (45 kDa) and a small C-terminal (4 kDa) fragment. The isolated porcine and human C-terminal domains were completely inactive towards emulsified tributyrin, though were able to bind colipase. Conversely, the horse 45 kDa fragment retained the lipase activity but failed to correctly bind colipase. This work definitely proves that catalysis and colipase binding are separate events involving topographically distinct regions of the molecule and focuses attention on the role of the C-terminal domain in colipase binding.     

A Monte Carlo simulation study of protein-induced heat capacity changes and lipid-induced protein clustering.

Monte Carlo simulations were used to describe the interaction of peripheral and integral proteins with lipids in terms of heat capacity profiles and protein distribution. The simulations were based on a two-state model for the lipid, representing the lipid state as being either gel or fluid. The interaction between neighboring lipids has been taken into account through an unlike nearest neighbor free energy term delta omega, which is a measure of the cooperativity of the lipid transition. Lipid/protein interaction was considered using the experimental observation that the transition midpoints of lipid membranes are shifted upon protein binding, a thermodynamic consequence of different binding constants of protein with fluid or gel lipids. The difference of the binding free energies was used as an additional parameter to describe lipid-protein interaction. The heat capacity profiles of lipid/protein complexes could be well described for both peripheral and integral proteins. Binding of proteins results in a shift and an asymmetric broadening of the melting profile. The model results in a coexistence of gel and fluid lipid domains in the proximity of the thermotropic transition. As a consequence, bound peripheral proteins aggregate in the temperature range of the lipid transition. Integral proteins induce calorimetric melting curves that are qualitatively different from that of peripheral proteins and aggregate in either gel or liquid crystalline lipid phase. The results presented here are in good agreement with calorimetric experiments on lipid-protein complexes and have implementations for the functional control of proteins.     

The C1B domains of novel PKCε and PKCη have a higher membrane binding affinity than those of the also novel PKCδ and PKCθ

The C1 domains of novel PKCs mediate the diacylglycerol-dependent translocation of these enzymes. The four different C1B domains of novel PKCs (δ, ε, θ and η) were studied, together with different lipid mixtures containing acidic phospholipids and diacylglycerol or phorbol ester. The results show that either in the presence or in the absence of diacylglycerol, C1Bε and C1Bη exhibit a substantially higher propensity to bind to vesicles containing negatively charged phospholipids than C1Bδ and C1Bθ. The observed differences between the C1B domains of novel PKCs (in two groups of two each) were also evident in RBL-2H3 cells and it was found that, as with model membranes, in which C1Bε and C1Bη could be translocated to membranes by the addition of a soluble phosphatidic acid without diacylglycerol or phorbol ester, C1Bδ and C1Bθ were not translocated when soluble phosphatidic acid was added, and diacylglycerol was required to achieve a detectable binding to cell membranes. It is concluded that two different subfamilies of novel PKCs can be established with respect to their propensity to bind to the cell membrane and that these peculiarities in recognizing lipids may explain why these isoenzymes are specialized in responding to different triggering signals and bind to different cell membranes.     

Structural Conservation and Effects of Alterations in T Cell Receptor Transmembrane Interfaces

T cell receptors (TCRs) are octameric assemblies of type-I membrane proteins in which a receptor heterodimer (αβ, δγ, or pre-Tαβ) is associated with three dimeric signaling modules (CD3δε, CD3γε, and ζζ) at the T cell or pre-T cell surface. In the human αβTCR, the α and β transmembrane (TM) domains form a specific structure that acts as a hub for assembly with the signaling modules inside the lipid bilayer. Conservation of key polar contacts across the C-terminal half of this TM interface suggests that the structure is a common feature of all TCR types. In this study, using molecular dynamics simulations in explicit lipid bilayers, we show that human δγ and pre-Tαβ TM domains also adopt stable αβ-like interfaces, yet each displays unique features that modulate the stability of the interaction and are related to sequences that are conserved within TCR types, but are distinct from the αβ sequences. We also performed simulations probing effects of previously reported mutations in the human αβ TM interface, and observed that the most disruptive mutations caused substantial departures from the wild-type TM structure and increased dynamics. These simulations show a strong correlation between structural instability, increased conformational variation, and the severity of structural defects in whole-TCR complexes measured in our previous biochemical assays. These results thus support the view that the stability of the core TM structure is a key determinant of TCR structural integrity and suggest that the interface has been evolutionarily optimized for different forms of TCRs.