Mapping structural and dynamic divergence across the MBOAT family

Membrane-bound O-acyltransferases (MBOATs) are membrane-embedded enzymes that catalyze acyl chain transfer to a diverse group of substrates, including lipids, small molecules, and proteins. MBOATs share a conserved structural core, despite wide-ranging functional specificity across both prokaryotes and eukaryotes. The structural basis of catalytic specificity, regulation and interactions with the surrounding environment remain uncertain. Here, we combine comparative molecular dynamics (MD) simulations with bioinformatics to assess molecular and interactional divergence across the family. In simulations, MBOATs differentially distort the bilayer depending on their substrate type. Additionally, we identify lipid binding sites surrounding reactant gates in the surrounding membrane. Complementary bioinformatic analyses reveal a conserved role for re-entrant loop-2 in MBOAT fold stabilization and a key hydrogen bond bridging DGAT1 dimerization. Finally, we predict differences in MBOAT solvation and water gating properties. These data are pertinent to the design of MBOAT-specific inhibitors that encompass dynamic information within cellular mimetic environments.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

The energetics and ion coupling of cholesterol transport through Patched1

Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, we find that transition between "inward-open" and solvent "occluded" states is accompanied by Na+-induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion coupling stoichiometries of PTCH1 transport mechanisms, whereby one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Rocking the MBOAT: Structural insights into the membrane bound O-acyltransferase family

The membrane-bound O-acyltransferase (MBOAT) superfamily catalyses the transfer of acyl chains to substrates implicated in essential cellular functions. Aberrant function of MBOATs is associated with various diseases and MBOATs are promising drug targets. There has been recent progress in structural characterisation of MBOATs, advancing our understanding of their functional mechanism. Integrating information across the MBOAT family, we characterise a common MBOAT fold and provide a blueprint for substrate and inhibitor engagement. This work provides context for the diverse substrates, mechanisms, and evolutionary relationships of protein and small-molecule MBOATs. Further work should aim to characterise MBOATs, as inherently lipid-associated proteins, within their membrane environment.

The Energetics and Ion Coupling of Cholesterol Transport Through Patched1

Patched1 (PTCH1) is the principal tumour suppressor protein of the mammalian Hedgehog (HH) signalling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the Class F G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating accessible cholesterol levels within ciliary membranes. Using extensive molecular dynamics (MD) simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15-20 kJ mol -1 for cholesterol export. In simulations we identify cation binding sites within the PTCH1 transmembrane domain (TMD) which may provide the energetic impetus for cholesterol transport. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between transmembrane motions and PTCH1 activity. Using complementary simulations of Dispatched1 (DISP1) we find that transition between 'inward-open' and solvent 'occluded' states is accompanied by Na + induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion-coupling stoichiometries of PTCH1 transport mechanisms, whereby 1-3 Na + or 2-3 K + couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Hepatocyte cholesterol content modulates glucagon receptor signalling

OBJECTIVE

To determine whether glucagon receptor (GCGR) actions are modulated by cellular cholesterol levels.

METHODS

We determined the effects of experimental cholesterol depletion and loading on glucagon-mediated cAMP production, ligand internalisation and glucose production in human hepatoma cells, mouse and human hepatocytes. GCGR interactions with lipid bilayers were explored using coarse-grained molecular dynamic simulations. Glucagon responsiveness was measured in mice fed a high cholesterol diet with or without simvastatin to modulate hepatocyte cholesterol content.

RESULTS

GCGR cAMP signalling was reduced by higher cholesterol levels across different cellular models. Ex vivo glucagon-induced glucose output from mouse hepatocytes was enhanced by simvastatin treatment. Mice fed a high cholesterol diet had increased hepatic cholesterol and a blunted hyperglycaemic response to glucagon, both of which were partially reversed by simvastatin. Simulations identified likely membrane-exposed cholesterol binding sites on the GCGR, including a site where cholesterol is a putative negative allosteric modulator.

CONCLUSIONS

Our results indicate that cellular cholesterol content influences glucagon sensitivity and indicate a potential molecular basis for this phenomenon. This could be relevant to the pathogenesis of non-alcoholic fatty liver disease, which is associated with both hepatic cholesterol accumulation and glucagon resistance.

Patched 1 regulates Smoothened by controlling sterol binding to its extracellular cysteine-rich domain

Smoothened (SMO) transduces the Hedgehog (Hh) signal across the plasma membrane in response to accessible cholesterol. Cholesterol binds SMO at two sites: one in the extracellular cysteine-rich domain (CRD) and a second in the transmembrane domain (TMD). How these two sterol-binding sites mediate SMO activation in response to the ligand Sonic Hedgehog (SHH) remains unknown. We find that mutations in the CRD (but not the TMD) reduce the fold increase in SMO activity triggered by SHH. SHH also promotes the photocrosslinking of a sterol analog to the CRD in intact cells. In contrast, sterol binding to the TMD site boosts SMO activity regardless of SHH exposure. Mutational and computational analyses show that these sites are in allosteric communication despite being 45 angstroms apart. Hence, sterols function as both SHH-regulated orthosteric ligands at the CRD and allosteric ligands at the TMD to regulate SMO activity and Hh signaling.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

Structure, mechanism, and inhibition of Hedgehog acyltransferase

The Sonic Hedgehog (SHH) morphogen pathway is fundamental for embryonic development and stem cell maintenance and is implicated in various cancers. A key step in signaling is transfer of a palmitate group to the SHH N terminus, catalyzed by the multi-pass transmembrane enzyme Hedgehog acyltransferase (HHAT). We present the high-resolution cryo-EM structure of HHAT bound to substrate analog palmityl-coenzyme A and a SHH-mimetic megabody, revealing a heme group bound to HHAT that is essential for HHAT function. A structure of HHAT bound to potent small-molecule inhibitor IMP-1575 revealed conformational changes in the active site that occlude substrate binding. Our multidisciplinary analysis provides a detailed view of the mechanism by which HHAT adapts the membrane environment to transfer an acyl chain across the endoplasmic reticulum membrane. This structure of a membrane-bound O-acyltransferase (MBOAT) superfamily member provides a blueprint for other protein-substrate MBOATs and a template for future drug discovery.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

The morphogen Sonic hedgehog inhibits its receptor Patched by a pincer grasp mechanism

Hedgehog (HH) ligands, classical morphogens that pattern embryonic tissues in all animals, are covalently coupled to two lipids-a palmitoyl group at the N terminus and a cholesteroyl group at the C terminus. While the palmitoyl group binds and inactivates Patched 1 (PTCH1), the main receptor for HH ligands, the function of the cholesterol modification has remained mysterious. Using structural and biochemical studies, along with reassessment of previous cryo-electron microscopy structures, we find that the C-terminal cholesterol attached to Sonic hedgehog (Shh) binds the first extracellular domain of PTCH1 and promotes its inactivation, thus triggering HH signaling. Molecular dynamics simulations show that this interaction leads to the closure of a tunnel through PTCH1 that serves as the putative conduit for sterol transport. Thus, Shh inactivates PTCH1 by grasping its extracellular domain with two lipidic pincers, the N-terminal palmitate and the C-terminal cholesterol, which are both inserted into the PTCH1 protein core.