Leukocyte integrin αLβ2 transmembrane association dynamics revealed by coarse-grained molecular dynamics simulations

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Lipid-dependent conformational landscape of the ErbB2 growth factor receptor dimers

Altered lipid metabolism has been linked to cancer development and progression. Several roles have been attributed to the increased saturation and length of lipid acyl tails observed in tumors, but its effect on signaling receptors is still emerging. In this work, we have analyzed the lipid dependence of the ErbB2 growth factor receptor dimerization that plays an important role in the pathogenesis of breast cancer. We have performed coarse-grain ensemble molecular dynamics simulations to comprehensively sample the ErbB2 monomer-dimer association. Our results indicate a dynamic dimer state with a complex conformational landscape that is modulated with increasing lipid tail length. We resolve the native N-terminal "active" and C-terminal "inactive" conformations in all membrane compositions. However, the relative population of the N-terminal and C-terminal conformers is dependent on length of the saturated lipid tails. In short-tail membranes, additional non-specific dimers are observed which are reduced or absent in long-tailed bilayers. Our results indicate that the relative population as well as the structure of the dimer state is modulated by membrane composition. We have correlated these differences to local perturbations of the membrane around the receptor. Our work is an important step in characterizing ErbB dimers in healthy and diseased states and emphasize the importance of sampling lipid dynamics in understanding receptor association.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

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.

On the contributing role of the transmembrane domain for subunit-specific sensitivity of integrin activation

Integrins are α/β heterodimeric transmembrane adhesion receptors. Evidence exists that their transmembrane domain (TMD) separates upon activation. Subunit-specific differences in activation sensitivity of integrins were reported. However, whether sequence variations in the TMD lead to differential TMD association has remained elusive. Here, we show by molecular dynamics simulations and association free energy calculations on TMDs of integrin αIIbβ3, αvβ3, and α5β1 that αIIbβ3 TMD is most stably associated; this difference is related to interaction differences across the TMDs. The order of TMD association stability is paralleled by the basal activity of these integrins, which suggests that TMD differences can have a decisive effect on integrin conformational free energies. We also identified a specific order of clasp disintegration upon TMD dissociation, which suggests that the closed state of integrins may comprise several microstates. Our results provide unprecedented insights into a possibly contributing role of TMD towards subunit-specific sensitivity of integrin activation.

Self-assembling study of sarcolipin and its mutants in multiple molecular dynamic simulations

The Sarcolipin (SLN) is a single trans-membrane protein that can self-assembly to dimer and oligomer for playing importantphysiological function. In this work, we addressed the dimerization of wild type SLN (wSLN) and its mutants (mSLNs) - I17A and I20A, using both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations. Our results demonstrated that wSLN homodimer assembled as a left-handed helical complex, while mSLNs heterodimers assembled as right-handed complexes. Analysis of residue-residue contacts map indicated that isoleucine (Ile)-leucione (Leu) zipper domain played an important role in dimerization. The potential of mean force (PMF) demonstrated that wSLN homodimer was more stable than mSLNs heterodimers. Meanwhile, the mSLNs heterodimers preferred right-handed rather than left-handed helix. AT-MD simulations for wSLN and mSLNs were also in line with CG-MD simulations. These results provided the insights for understanding the mechanisms of SLNs self-assembling. Proteins 2017; 85:1065-1077. © 2017 Wiley Periodicals, Inc.

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2-F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2-F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes.

Homo-oligomerization of transmembrane α-domain of integrin

Integrins contribute to form focal adhesions complex. Therefore, simulation of integrin interactions can be helpful in clarifying the mechanism of focal adhesion formation. Interactions of integrins can also initiate signal transduction in the focal adhesions. Since integrins contain α and β subunits that are separated in an active state, studying both subunits separately is crucial, since, in the active state of integrins, the distance between these subunits is long enough that they do not influence one another significantly. Thus, this study aims to investigate the tendency of α subunits of integrins to form homodimers. All simulations were carried out via MARTINI coarse grain (CG) molecular dynamics technique. α subunits were placed in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer at a distance of 5 nm, and they were allowed to diffuse in the lipid bilayer. All simulations showed that α subunits have a tendency to form stable dimers.

Cholesterol modulates the dimer interface of the β₂-adrenergic receptor via cholesterol occupancy sites

The β2-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β2-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β2-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.

Molecular dynamic simulation of the self-assembly of DAP12-NKG2C activating immunoreceptor complex

The DAP12-NKG2C activating immunoreceptor complex is one of the multisubunit transmembrane protein complexes in which ligand-binding receptor chains assemble with dimeric signal-transducing modules through non-covalent associations in their transmembrane (TM) domains. In this work, both coarse grained and atomistic molecular dynamic simulation methods were applied to investigate the self-assembly dynamics of the transmembrane domains of the DAP12-NKG2C activating immunoreceptor complex. Through simulating the dynamics of DAP12-NKG2C TM heterotrimer and point mutations, we demonstrated that a five-polar-residue motif including: 2 Asps and 2 Thrs in DAP12 dimer, as well as 1 Lys in NKG2C TM plays an important role in the assembly structure of the DAP12-NKG2C TM heterotrimer. Furthermore, we provided clear evidences to exclude the possibility that another NKG2C could stably associate with the DAP12-NKG2C heterotrimer. Based on the simulation results, we proposed a revised model for the self-assembly of DAP12-NKG2C activating immunoreceptor complex, along with a plausible explanation for the association of only one NKG2C with a DAP12 dimer.

The free energy landscape of dimerization of a membrane protein, NanC

Membrane proteins are frequently present in crowded environments, which favour lateral association and, on occasions, two-dimensional crystallization. To better understand the non-specific lateral association of a membrane protein we have characterized the free energy landscape for the dimerization of a bacterial outer membrane protein, NanC, in a phospholipid bilayer membrane. NanC is a member of the KdgM-family of bacterial outer membrane proteins and is responsible for sialic acid transport in E. coli. Umbrella sampling and coarse-grained molecular dynamics were employed to calculate the potentials of mean force (PMF) for a variety of restrained relative orientations of two NanC proteins as the separation of their centres of mass was varied. We found the free energy of dimerization for NanC to be in the range of to . Differences in the depths of the PMFs for the various orientations are related to the shape of the proteins. This was quantified by calculating the lipid-inaccessible buried surface area of the proteins in the region around the minimum of each PMF. The depth of the potential well of the PMF was shown to depend approximately linearly on the buried surface area. We were able to resolve local minima in the restrained PMFs that would not be revealed using conventional umbrella sampling. In particular, these features reflected the local organization of the intervening lipids between the two interacting proteins. Through a comparison with the distribution of lipids around a single freely-diffusing NanC, we were able to predict the location of these restrained local minima for the orientational configuration in which they were most pronounced. Our ability to make this prediction highlights the important role that lipid organization plays in the association of two NanCs in a bilayer.

Cells are surrounded by selectively-permeable bilayer membranes, enabling the cell to control its internal environment. Embedded within these membranes are a variety of membrane proteins, many of which facilitate this environmental control and are integral to numerous metabolic processes. Their location within the membrane and their mutual association are controlled by many factors. We use molecular dynamics simulations to investigate the free energy of association for a pair of relatively simple membrane proteins. By doing so, we are able to characterize the effect that the geometrical properties of the protein have on their mutual association in a bilayer environment, showing that there is a correlation between the buried surface area of two proteins when in contact and the strength of their interaction. We also observe the effect of protein-lipid-protein interactions in this free energy characterization. Such interactions are related to the preferential distribution of lipids around proteins in the membrane.

Coarse-grain modelling of protein-protein interactions

Here, we review recent advances towards the modelling of protein-protein interactions (PPI) at the coarse-grained (CG) level, a technique that is now widely used to understand protein affinity, aggregation and self-assembly behaviour. PPI models of soluble proteins and membrane proteins are separately described, but we note the parallel development that is present in both research fields with three important themes: firstly, combining CG modelling with knowledge-based approaches to predict and refine protein-protein complexes; secondly, using physics-based CG models for de novo prediction of protein-protein complexes; and thirdly modelling of large scale protein aggregates.

The association of polar residues in the DAP12 homodimer: TOXCAT and molecular dynamics simulation studies

Dimerization of the transmembrane (TM) adaptor protein DAP12 plays a key role in mediating activation signals through TM-TM association with cell-surface receptors. Herein, we apply the TOXCAT assay and molecular dynamics simulation to analyze dynamics and dimerization of the TM helix of DAP12 in the membrane bilayer. In the TOXCAT assay, we performed site-specific mutagenesis of potential dimerization motifs in the DAP12 TM domain. Instead of the common GxxxG dimerization motif, mutating either of the polar residues Asp-50 and Thr-54 significantly decreased the TOXCAT signal for the dimerization of DAP12 TM domain. Furthermore, through the conformational difference between wild-type and mutant DAP12 TM homodimers, a combined coarse-grained and atomistic molecular dynamics simulation has identified both Asp-50 and Thr-54 at the dimerization interface. The experimental and computational results of the DAP12 TM dimer are in excellent agreement with the previously reported NMR structure obtained in detergent micelles. Such a combination of dynamics simulation and cell-based experiments can be applied to produce insights at the molecular level into the TM-TM association of many other transmembrane proteins.

Transmembrane and Juxtamembrane Structure of αL Integrin in Bicelles

The accepted model for the interaction of α and β integrins in the transmembrane (TM) domain is based on the pair αIIbβ3. This involves the so-called outer and inner membrane association clasps (OMC and IMC, respectively). In the α chain, the OMC involves a GxxxG-like motif, whereas in the IMC a conserved juxtamembrane GFFKR motif experiences a backbone reversal that partially fills the void generated by TM separation towards the cytoplasmic half. However, the GFFKR motif of several α integrin cytoplasmic tails in non-bicelle environments has been shown to adopt an α-helical structure that is not membrane-embedded and which was shown to bind a variety of cytoplasmic proteins. Thus it is not known if a membrane-embedded backbone reversal is a conserved structural feature in α integrins. We have studied the system αLβ2 because of its importance in leukocytes, where integrin deactivation is particularly important. Herein we show that the backbone reversal feature is not only present in αIIb but also in αL-TM when reconstituted in bicelles. Additionally, titration with β2 TM showed eight residues clustering along one side of αL-TM, forming a plausible interacting face with β2. The latter orientation is consistent with a previously predicted reported polar interaction between αL Ser-1071 and β2 Thr-686.

The leucocyte β2 (CD18) integrins: the structure, functional regulation and signalling properties

Leucocytes are highly motile cells. Their ability to migrate into tissues and organs is dependent on cell adhesion molecules. The integrins are a family of heterodimeric transmembrane cell adhesion molecules that are also signalling receptors. They are involved in many biological processes, including the development of metazoans, immunity, haemostasis, wound healing and cell survival, proliferation and differentiation. The leucocyte-restricted β2 integrins comprise four members, namely αLβ2, αMβ2, αXβ2 and αDβ2, which are required for a functional immune system. In this paper, the structure, functional regulation and signalling properties of these integrins are reviewed.

Transmembrane domains interactions within the membrane milieu: principles, advances and challenges

Protein-protein interactions within the membrane are involved in many vital cellular processes. Consequently, deficient oligomerization is associated with known diseases. The interactions can be partially or fully mediated by transmembrane domains (TMD). However, in contrast to soluble regions, our knowledge of the factors that control oligomerization and recognition between the membrane-embedded domains is very limited. Due to the unique chemical and physical properties of the membrane environment, rules that apply to interactions between soluble segments are not necessarily valid within the membrane. This review summarizes our knowledge on the sequences mediating TMD-TMD interactions which include conserved motifs such as the GxxxG, QxxS, glycine and leucine zippers, and others. The review discusses the specific role of polar, charged and aromatic amino acids in the interface of the interacting TMD helices. Strategies to determine the strength, dynamics and specificities of these interactions by experimental (ToxR, TOXCAT, GALLEX and FRET) or various computational approaches (molecular dynamic simulation and bioinformatics) are summarized. Importantly, the contribution of the membrane environment to the TMD-TMD interaction is also presented. Studies utilizing exogenously added TMD peptides have been shown to influence in vivo the dimerization of intact membrane proteins involved in various diseases. The chirality independent TMD-TMD interactions allows for the design of novel short d- and l-amino acids containing TMD peptides with advanced properties. Overall these studies shed light on the role of specific amino acids in mediating the assembly of the TMDs within the membrane environment and their contribution to protein function. This article is part of a Special Issue entitled: Protein Folding in Membranes.

A helix heterodimer in a lipid bilayer: prediction of the structure of an integrin transmembrane domain via multiscale simulations

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling.