mPPases create a conserved anionic membrane fingerprint as identified via multiscale simulations

Publisher Copyright: Copyright: © 2022 Holmes et al.Membrane-integral pyrophosphatases (mPPases) are membrane-bound enzymes responsible for hydrolysing inorganic pyrophosphate and translocating a cation across the membrane. Their function is essential for the infectivity of clinically relevant protozoan parasites and plant maturation. Recent developments have indicated that their mechanism is more complicated than previously thought and that the membrane environment may be important for their function. In this work, we use multiscale molecular dynamics simulations to demonstrate for the first time that mPPases form specific anionic lipid interactions at 4 sites at the distal and interfacial regions of the protein. These interactions are conserved in simulations of the mPPases from Thermotoga maritima, Vigna radiata and Clostridium leptum and characterised by interactions with positive residues on helices 1, 2, 3 and 4 for the distal site, or 9, 10, 13 and 14 for the interfacial site. Due to the importance of these helices in protein stability and function, these lipid interactions may play a crucial role in the mPPase mechanism and enable future structural and functional studies.Peer reviewe

The Function of Membrane Integral Pyrophosphatases From Whole Organism to Single Molecule

Membrane integral pyrophosphatases (mPPases) are responsible for the hydrolysis of pyrophosphate. This enzymatic mechanism is coupled to the pumping of H+ or Na+ across membranes in a process that can be K+ dependent or independent. Understanding the movements and dynamics throughout the mPPase catalytic cycle is important, as this knowledge is essential for improving or impeding protein function. mPPases have been shown to play a crucial role in plant maturation and abiotic stress tolerance, and so have the potential to be engineered to improve plant survival, with implications for global food security. mPPases are also selectively toxic drug targets, which could be pharmacologically modulated to reduce the virulence of common human pathogens. The last few years have seen the publication of many new insights into the function and structure of mPPases. In particular, there is a new body of evidence that the catalytic cycle is more complex than originally proposed. There are structural and functional data supporting a mechanism involving half-of-the-sites reactivity, inter-subunit communication, and exit channel motions. A more advanced and in-depth understanding of mPPases has begun to be uncovered, leaving the field of research with multiple interesting avenues for further exploration and investigation.Peer reviewe

Association of Peripheral Membrane Proteins with Membranes: Free Energy of Binding of GRP1 PH Domain with Phosphatidylinositol Phosphate-Containing Model Bilayers

Understanding the energetics of peripheral protein-membrane interactions is important to many areas of biophysical chemistry and cell biology. Estimating free-energy landscapes by molecular dynamics (MD) simulation is challenging for such systems, especially when membrane recognition involves complex lipids, e.g., phosphatidylinositol phosphates (PIPs). We combined coarse-grained MD simulations with umbrella sampling to quantify the binding of the well-explored GRP1 pleckstrin homology (PH) domain to model membranes containing PIP molecules. The experimentally observed preference of GRP1-PH for PIP3 over PIP2 was reproduced. Mutation of a key residue (K273A) within the canonical PIP-binding site significantly reduced the free energy of PIP binding. The presence of a noncanonical PIP-interaction site, observed experimentally in other PH domains but not previously in GRP1-PH, was also revealed. These studies demonstrate how combining coarse-grained simulations and umbrella sampling can unmask the molecular basis of the energetics of interactions between peripheral membrane proteins and complex cellular membranes

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

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK

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

Structure and lipid binding properties of the kindlin-3 pleckstrin homology domain

Kindlins co-activate integrins alongside talin. They possess, like talin, a FERM domain comprising F0-F3 subdomains, but with a pleckstrin homology (PH) domain inserted in the F2 subdomain that enables membrane association. We present the crystal structure of murine kindlin-3 PH domain determined at 2.23Å resolution and characterise its lipid binding using biophysical and computational approaches. Molecular dynamics (MD) simulations suggest flexibility in the PH domain loops connecting β-strands forming the putative phosphatidylinositol phosphate (PtdInsP) binding site. Simulations with PtdInsP-containing bilayers reveal that the PH domain associates with PtdInsP molecules mainly via the positively charged surface presented by the β1-β2 loop and that it binds with somewhat higher affinity to PtdIns(3,4,5)P3 compared to PtdIns(4,5)P2. Surface plasmon resonance (SPR) with lipid headgroups immobilised and the PH domain as analyte indicate affinities of 300 μM for PtdIns(3,4,5)P3 and 1mM for PtdIns(4,5)P2. In contrast, SPR studies with immobilised PH domain and lipid nanodiscs as analyte show affinities of 0.40 µM for PtdIns(3,4,5)P3 and no affinity for PtdIns(4,5)P2 when the inositol phosphate constitutes 5% of the total lipids (~5 molecules per nanodisc). Reducing the PtdIns(3,4,5)P3 composition to 1% abolishes nanodisc binding to the PH domain, as does site-directed mutagenesis of two lysines within the β1-β2 loop. Binding of PtdIns(3,4,5)P3 by a canonical PH domain, Grp1, is not similarly influenced by SPR experimental design. These data suggest a role for PtdIns(3,4,5)P3 clustering in the binding of some PH domains and not others, highlighting the importance lipid mobility and clustering for the biophysical assessment of protein-membrane interactions

Effects of Periplasmic Chaperones and Membrane Thickness on BamA-Catalyzed Outer-Membrane Protein Folding

The biogenesis of outer-membrane proteins (OMPs) in gram-negative bacteria involves delivery by periplasmic chaperones to the β-barrel assembly machinery (BAM), which catalyzes OMP insertion into the outer membrane. Here, we examine the effects of membrane thickness, the Escherichia coli periplasmic chaperones Skp and SurA, and BamA, the central subunit of the BAM complex, on the folding kinetics of a model OMP (tOmpA) using fluorescence spectroscopy, native mass spectrometry, and molecular dynamics simulations. We show that prefolded BamA promotes the release of tOmpA from Skp despite the nM affinity of the Skp:tOmpA complex. This activity is located in the BamA β-barrel domain, but is greater when full-length BamA is present, indicating that both the β-barrel and polypeptide transport-associated (POTRA) domains are required for maximal activity. By contrast, SurA is unable to release tOmpA from Skp, providing direct evidence against a sequential chaperone model. By varying lipid acyl chain length in synthetic liposomes we show that BamA has a greater catalytic effect on tOmpA folding in thicker bilayers, suggesting that BAM catalysis involves lowering of the kinetic barrier imposed by the hydrophobic thickness of the membrane. Consistent with this, molecular dynamics simulations reveal that increases in membrane thinning/disorder by the transmembrane domain of BamA is greatest in thicker bilayers. Finally, we demonstrate that cross-linking of the BamA barrel does not affect tOmpA folding kinetics in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes, suggesting that lateral gating of the BamA barrel and/or hybrid barrel formation is not required, at least for the assembly of a small 8-stranded OMP in vitro

Computational studies of talin-mediated integrin activation

Integrins are large heterodimeric (αβ) cell surface receptors that play a key role in the formation of focal adhesion complexes and are involved in various signal transduction pathways. They are ‘activated’ to a high affinity state by the formation of an intracellular complex between the membrane, the integrin β-subunit tail and talin, a process known as ‘inside-out activation’. The head domain of talin, a FERM domain homologue, plays a vital role in the formation of this complex. Recent studies also suggest that kindlins act in synergy with talin to induce integrin activation. Despite much available structural and functional data, details of how talin activates integrins remain elusive. In this thesis a multiscale simulation approach (using a combination of coarse-grained and atomistic molecular dynamics simulations) together with NMR experiments were employed to study talin-mediated integrin inside-out activation. A number of novel insights emerged from these studies including: (i) the crucial role of negatively charged lipids in talin/membrane association; (ii) a novel V-shape conformation of the talin head domain which optimizes its interactions with negatively charged lipids; (iii) that interactions of talin with negatively charged moieties in the membrane orient talin to optimize interactions with the β cytoplasmic tail; (iv) that binding of talin to the β cytoplasmic tail promotes rearrangement of the integrin TM helices and weakens the integrin α/β association; and (v) that an increase in the tilt angle of the β integrin TM helix initiates a scissoring movement of the two integrin TM helices. These results, combined with experimental data, provide new insights into the mechanism of integrin inside-out activation. The same multiscale approach was used to demonstrate the crucial role of the Gx3G motif in the packing of the integrin transmembrane helices. Using recent structural data the integrin/talin complex was modelled and inserted in bilayers which resemble the biological plasma membrane. The results demonstrate the dynamic nature of the integrin receptor and suggest that the integrin/talin complex alters the lipid organization and motion in the outer and inner bilayer leaflets in an asymmetric way and that diffusion of lipids in the vicinity of the protein is slowed down. The work in this thesis demonstrates that multiscale simulations have considerable strengths when applied to investigations of structure/function relationships in membrane proteins.</p