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

Mechanism of substrate binding and transport in BASS transporters

The Bile Acid Sodium Symporter (BASS) family transports a wide array of molecules across membranes, including bile acids in humans, and small metabolites in plants. These transporters, many of which are sodium-coupled, have been shown to use an elevator mechanism of transport, but exactly how substrate binding is coupled to sodium ion binding and transport is not clear. Here we solve the crystal structure at 2.3 Å of a transporter from Neisseria Meningitidis (ASBTNM) in complex with pantoate, a potential substrate of ASBTNM. The BASS family is characterised by two helices that cross-over in the centre of the protein in an arrangement that is intricately held together by two sodium ions. We observe that the pantoate binds, specifically, between the N-termini of two of the opposing helices in this cross-over region. During molecular dynamics simulations the pantoate remains in this position when sodium ions are present but is more mobile in their absence. Comparison of structures in the presence and absence of pantoate demonstrates that pantoate elicits a conformational change in one of the cross-over helices. This modifies the interface between the two domains that move relative to one another to elicit the elevator mechanism. These results have implications, not only for ASBTNM but for the BASS family as a whole and indeed other transporters that work through the elevator mechanism

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

Building a community-driven ecosystem for fast, reproducible, and reusable molecular simulation analysis using mdanalysis

A poster given at the 2023 BPS conference on MDAnalysis and the MDAKits ecosystem