Cholesterol Interaction Sites on the Transmembrane Domain of the Hedgehog Signal Transducer and Class F G Protein-Coupled Receptor Smoothened

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

De novo variants in GABRA4 are associated with a neurological phenotype including developmental delay, behavioral abnormalities and epilepsy

Nine out of 19 genes encoding GABAA receptor subunits have been linked to monogenic syndromes characterized by seizures and developmental disorders. Previously, we reported the de novo variant p.(Thr300Ile) in GABRA4 in a patient with epilepsy and neurodevelopmental abnormalities. However, no new cases have been reported since then. Through an international collaboration, we collected molecular and phenotype data of individuals carrying de novo variants in GABRA4. Patients and their parents were investigated either by exome or genome sequencing, followed by targeted Sanger sequencing in some cases. All variants within the transmembrane domain, including the previously reported p.(Thr300Ile) variant, were characterized in silico and analyzed by molecular dynamics (MD) simulation studies. We identified three novel de novo missense variants in GABRA4 (NM_000809.4): c.797 C > T, p.(Pro266Leu), c.899 C > A, p.(Thr300Asn), and c.634 G > A, p.(Val212Ile). The p.(Thr300Asn) variant impacts the same codon as the previously reported variant p.(Thr300Ile) and likely arose post-zygotically as evidenced by sequencing oral mucosal cells. Overlapping phenotypes among affected individuals included developmental delay (4/4), epileptiform EEG abnormalities (3/4), attention deficits (3/4), seizures (2/4), autistic features (2/4) and structural brain abnormalities (2/4). MD simulations of the three variants within the transmembrane domain of the receptor indicate that sub-microsecond scale dynamics differ between wild-type and mutated subunits. Taken together, our findings further corroborate an association between GABRA4 and a neurological phenotype including variable neurodevelopmental, behavioral and epileptic abnormalities.

Estimating the Cholesterol Affinity of Integral Membrane Proteins from Experimental Data

The cholesterol affinities of many integral plasma membrane proteins have been estimated by molecular computation. However, these values lack experimental confirmation. We therefore developed a simple mathematical model to extract sterol affinity constants and stoichiometries from published isotherms for the dependence of the activity of such proteins on the membrane cholesterol concentration. The binding curves for these proteins are sigmoidal, with strongly lagged thresholds attributable to competition for the cholesterol by bilayer phospholipids. The model provided isotherms that matched the experimental data using published values for the sterol association constants and stoichiometries of the phospholipids. Three oligomeric transporters were found to bind cholesterol without cooperativity, with dimensionless association constants of 35 for Kir3.4* and 100 for both Kir2 and a GAT transporter. (The corresponding ΔG° values were -8.8, -11.4, and -11.4 kJ/mol, respectively). These association constants are significantly lower than those for the phospholipids, which range from ∼100 to 6000. The BK channel, the nicotinic acetylcholine receptor, and the M192I mutant of Kir3.4* appear to bind multiple cholesterol molecules cooperatively (n = 2 or 4), with subunit affinities of 563, 950, and 700, respectively. The model predicts that the three less avid transporters are approximately half-saturated in their native plasma membranes; hence, they are sensitive to variations in cholesterol in vivo. The more avid proteins would be nearly saturated in vivo. The method can be applied to any integral protein or other ligands in any bilayer for which there are reasonable estimates of the sterol affinities and stoichiometries of the phospholipids.

Palmitoylation of the Glucagon-like Peptide-1 Receptor Modulates Cholesterol Interactions at the Receptor-Lipid Microenvironment

The G protein-coupled receptor (GPCR) superfamily of cell surface receptors has been shown to be functionally modulated by post-translational modifications. The glucagon-like peptide receptor-1 (GLP-1R), which is a drug target in diabetes and obesity, undergoes agonist-dependent palmitoyl tail conjugation. The palmitoylation in the C-terminal domain of GLP-1R has been suggested to modulate the receptor-lipid microenvironment. In this work, we have performed coarse-grain molecular dynamics simulations of palmitoylated and nonpalmitoylated GLP-1R to analyze the differential receptor-lipid interactions. Interestingly, the placement and dynamics of the C-terminal domain of GLP-1R are found to be directly dependent on the palmitoyl tail. We observe that both cholesterol and phospholipids interact with the receptor but display differential interactions in the presence and absence of the palmitoyl tail. We characterize important cholesterol-binding sites and validate sites that have been previously reported in experimentally resolved structures of the receptor. We show that the receptor acts like a conduit for cholesterol flip-flop by stabilizing cholesterol in the membrane core. Taken together, our work represents an important step in understanding the molecular effects of lipid modifications in GPCRs.

Fully activated structure of the sterol-bound Smoothened GPCR-Gi protein complex

Smoothened (SMO) is an oncoprotein and signal transducer in the Hedgehog signaling pathway that regulates cellular differentiation and embryogenesis. As a member of the Frizzled (Class F) family of G protein-coupled receptors (GPCRs), SMO biochemically and functionally interacts with Gi family proteins. However, key molecular features of fully activated, G protein-coupled SMO remain elusive. We present the atomistic structure of activated human SMO complexed with the heterotrimeric Gi protein and two sterol ligands, equilibrated at 310 K in a full lipid bilayer at physiological salt concentration and pH. In contrast to previous experimental structures, our equilibrated SMO complex exhibits complete breaking of the pi-cation interaction between R4516.32 and W5357.55, a hallmark of Class F receptor activation. The Gi protein couples to SMO at seven strong anchor points similar to those in Class A GPCRs: intracellular loop 1, intracellular loop 2, transmembrane helix 6, and helix 8. On the path to full activation, we find that the extracellular cysteine-rich domain (CRD) undergoes a dramatic tilt, following a trajectory suggested by positions of the CRD in active and inactive experimental SMO structures. Strikingly, a sterol ligand bound to a shallow transmembrane domain (TMD) site in the initial structure migrates to a deep TMD pocket found exclusively in activator-bound SMO complexes. Thus, our results indicate that SMO interacts with Gi prior to full activation to break the molecular lock, form anchors with Gi subunits, tilt the CRD, and facilitate migration of a sterol ligand in the TMD to an activated position.

Putative Role of Cholesterol in Shaping the Structural and Functional Dynamics of Smoothened (SMO)

The smoothened (SMO) receptor belongs to the superfamily of class F G protein-coupled receptors (GPCRs) and is a potential drug target in several types of cancer. It has two ligand binding sites, respectively, in the cysteine-rich domain (CRD) and the transmembrane domain (TMD). It has been shown that cholesterol is important for its activation and function. However, the molecular-level understanding of SMO dynamics in the presence of cholesterol has not been explored in sufficient detail. In this work, we have carried out atomistic molecular dynamics simulations totaling 3.6 μs to analyze the effect of cholesterol binding to TMD and/or CRD on the structure and dynamics of the SMO receptor. Our results show that the presence of cholesterol in the CRD and TMD, respectively, alters the conformational dynamics of SMO differently. We reported that the reorganization of the D-R-E network at the extracellular end of the TMD is important for the high activity of SMO. In general, the transmembrane helices 5, 6, and 7 and helix 8 are most affected, which, in turn, leads to changes in the CRD and intracellular cytoplasmic domain (ICD) dynamics patterns depending on the presence or absence of cholesterol in the CRD and/or the TMD. We have also reported that the interaction of membrane lipids with SMO is different in different SMO states. These results agree with the experimental structural observations and data of cholesterol-bound and unbound structures of SMO and add to our molecular understanding of the SMO-cholesterol interaction.

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.

Cholesterol catalyzes unfolding in membrane-inserted motifs of the pore forming protein cytolysin A

Plasma membrane-induced protein folding and conformational transitions play a central role in cellular homeostasis. Several transmembrane proteins are folded in the complex lipid milieu to acquire a specific structure and function. Bacterial pore forming toxins (PFTs) are proteins expressed by a large class of pathogenic bacteria that exploit the plasma membrane environment to efficiently undergo secondary structure changes, oligomerize, and form transmembrane pores. Unregulated pore formation causes ion imbalance, leading to cell death and infection. Determining the free energy landscape of these membrane-driven-driven transitions remains a challenging problem. Although cholesterol recognition is required for lytic activity of several proteins in the PFT family of toxins, the regulatory role of cholesterol for the α-PFT, cytolysin A expressed by Escherichia coli remains unexplained. In a recent free energy computation, we showed that the β tongue, a critical membrane-inserted motif of the ClyA toxin, has an on-pathway partially unfolded intermediate that refolds into the helix-turn-helix motif of the pore state. To understand the molecular role played by cholesterol, we carry out string-method-based computations in membranes devoid of cholesterol, which reveals an increase of ∼30 times in the free energy barrier for the loss of β sheet secondary structure when compared with membranes containing cholesterol. Specifically, the tyrosine-cholesterol interaction was found to be critical to creating the unfolded intermediate. Cholesterol also increases the packing and hydrophobicity of the bilayer, resulting in enhanced interactions of the bound protein before complete membrane insertion. Our study illustrates that cholesterol is critical to catalyzing and stabilizing the membrane-inserted unfolded state of the β tongue motif of ClyA, opening up fresh insights into cholesterol-assisted unfolding of membrane proteins.

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.

Dual mechanisms of cholesterol-GPCR interactions that depend on membrane phospholipid composition

Cholesterol is a critical component of mammalian cell membranes and an allosteric modulator of G protein-coupled receptors (GPCRs), but divergent views exist on the mechanisms by which cholesterol influences receptor functions. Leveraging the benefits of lipid nanodiscs, i.e., quantitative control of lipid composition, we observe distinct impacts of cholesterol in the presence and absence of anionic phospholipids on the function-related conformational dynamics of the human A2A adenosine receptor (A2AAR). Direct receptor-cholesterol interactions drive activation of agonist-bound A2AAR in membranes containing zwitterionic phospholipids. Intriguingly, the presence of anionic lipids attenuates cholesterol's impact through direct interactions with the receptor, highlighting a more complex role for cholesterol that depends on membrane phospholipid composition. Targeted amino acid replacements at two frequently predicted cholesterol interaction sites showed distinct impacts of cholesterol at different receptor locations, demonstrating the ability to delineate different roles of cholesterol in modulating receptor signaling and maintaining receptor structural integrity.

Activation mechanism of the human Smoothened receptor

Smoothened (SMO) is a membrane protein of the class F subfamily of G protein-coupled receptors (GPCRs) and maintains homeostasis of cellular differentiation. SMO undergoes conformational change during activation, transmitting the signal across the membrane, making it amenable to bind to its intracellular signaling partner. Receptor activation has been studied at length for class A receptors, but the mechanism of class F receptor activation remains unknown. Agonists and antagonists bound to SMO at sites in the transmembrane domain (TMD) and the cysteine-rich domain have been characterized, giving a static view of the various conformations SMO adopts. Although the structures of the inactive and active SMO outline the residue-level transitions, a kinetic view of the overall activation process remains unexplored for class F receptors. We describe SMO's activation process in atomistic detail by performing 300 μs of molecular dynamics simulations and combining it with Markov state model theory. A molecular switch, conserved across class F and analogous to the activation-mediating D-R-Y motif in class A receptors, is observed to break during activation. We also show that this transition occurs in a stage-wise movement of the transmembrane helices: TM6 first, followed by TM5. To see how modulators affect SMO activity, we simulated agonist and antagonist-bound SMO. We observed that agonist-bound SMO has an expanded hydrophobic tunnel in SMO's core TMD, whereas antagonist-bound SMO shrinks this tunnel, further supporting the hypothesis that cholesterol travels through a tunnel inside Smoothened to activate it. In summary, this study elucidates the distinct activation mechanism of class F GPCRs and shows that SMO's activation process rearranges the core TMD to open a hydrophobic conduit for cholesterol transport.

Cholesterol in Class C GPCRs: Role, Relevance, and Localization

G-protein coupled receptors (GPCRs), one of the largest superfamilies of cell-surface receptors, are heptahelical integral membrane proteins that play critical roles in virtually every organ system. G-protein-coupled receptors operate in membranes rich in cholesterol, with an imbalance in cholesterol level within the vicinity of GPCR transmembrane domains affecting the structure and/or function of many GPCRs, a phenomenon that has been linked to several diseases. These effects of cholesterol could result in indirect changes by altering the mechanical properties of the lipid environment or direct changes by binding to specific sites on the protein. There are a number of studies and reviews on how cholesterol modulates class A GPCRs; however, this area of study is yet to be explored for class C GPCRs, which are characterized by a large extracellular region and often form constitutive dimers. This review highlights specific sites of interaction, functions, and structural dynamics involved in the cholesterol recognition of the class C GPCRs. We summarize recent data from some typical family members to explain the effects of membrane cholesterol on the structural features and functions of class C GPCRs and speculate on their corresponding therapeutic potential.

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.

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.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations

We used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A_2A adenosine receptor (A_2AR) and the A₁ adenosine receptor (A₁R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP₂)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A_2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 μs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A_2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A₁R, three binding sites, and for the inactive A₁R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A_2AR and A₁R. The stability of the identified binding sites to A₁R or A_2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A_2AR have shallow binding free energy minima compared to the related binding sites in A₁R, suggesting a stronger binding for cholesterol to A₁R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A₁R and A_2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A₁R (i.e., generally higher association rates) compared to A_2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A_2AR and A₁R. Taken together, binding sites of active A₁R may be considered as promising allosteric targets.

The impact of bilayer composition on the dimerization properties of the Slg1 stress sensor TMD from a multiscale analysis

The impact of mutual interactions between the transmembrane domains of membrane proteins and lipids on bilayer properties has gained major attraction. Most simulation studies of membranes rely on the Martini force field, which has proven extremely helpful in providing molecular insights into realistic systems. Accordingly, an evaluation of the accuracy of the Martini force field is crucial to be able to correctly interpret the reported data. In this study, we combine atomistic and coarse-grained Martini simulations to investigate the properties of transmembrane domains (TMDs) in a model yeast membrane. The results show that the TMD binding state (monomeric and dimeric with positive or negative crossing angle) and the membrane composition significantly influence the properties around the TMDs and change TMD-TMD and TMD-lipid affinities. Furthermore, ergosterol (ERG) exhibits a strong affinity to TMD dimers. Importantly, the right-handed TMD dimer configuration is stabilized via TMD-TMD contacts by the addition of asymmetric anionic phosphatidylserine (PS). The coarse-grained simulations corroborate many of these findings, with two notable exceptions: a systematic overestimation of TMD-ERG interaction and lack of stabilization of the right-handed TMD dimers with the addition of PS.

Lipid Raft Facilitated Receptor Organization and Signaling: A Functional Rheostat in Embryonic Development, Stem Cell Biology and Cancer

Molecular views of plasma membrane organization and dynamics are gradually changing over the past fifty years. Dynamics of plasma membrane instigate several signaling nexuses in eukaryotic cells. The striking feature of plasma membrane dynamics is that, it is internally transfigured into various subdomains of clustered macromolecules. Lipid rafts are nanoscale subdomains, enriched with cholesterol and sphingolipids, reside as floating entity mostly on the exoplasmic leaflet of the lipid bilayer. In terms of functionality, lipid rafts are unique among other membrane subdomains. Herein, advances on the roles of lipid rafts in cellular physiology and homeostasis are discussed, precisely, on how rafts dynamically harbor signaling proteins, including GPCRs, catalytic receptors, and ionotropic receptors within it and orchestrate multiple signaling pathways. In the developmental proceedings signaling are designed for patterning of overall organism and they differ from the somatic cell physiology and signaling of fully developed organisms. Some of the developmental signals are characteristic in maintenance of stemness and activated during several types of tumor development and cancer progression. The harmony between extracellular signaling and lineage specific transcriptional programs are extremely important for embryonic development. The roles of plasma membrane lipid rafts mediated signaling in lineage specificity, early embryonic development, stem cell maintenance are emerging. In view of this, we have highlighted and analyzed the roles of lipid rafts in receptor organization, cell signaling, and gene expression during embryonic development; from pre-implantation through the post-implantation phase, in stem cell and cancer biology.

Synergistic and Competitive Lipid Interactions in the Serotonin1A Receptor Microenvironment

The interaction of lipids with G-protein-coupled receptors (GPCRs) has been shown to modulate and dictate several aspects of GPCR organization and function. Diverse lipid interaction sites have been identified from structural biology, bioinformatics, and molecular dynamics studies. For example, multiple cholesterol interaction sites have been identified in the serotonin_1A receptor, along with distinct and overlapping sphingolipid interaction sites. How these lipids interact with each other and what is the resultant effect on the receptor is still not clear. In this work, we have analyzed lipid-lipid crosstalk at the receptor of the serotonin_1A receptor embedded in a membrane bilayer that mimics the neuronal membrane composition by long coarse-grain simulations. Using a set of similarity coefficients, we classified lipids that bind at the receptor together as synergistic cobinding, and those that bind individually as competitive. Our results show that certain lipids interact with the serotonin_1A receptor in synergy with each other. Not surprisingly, the ganglioside GM1 and cholesterol show a synergistic cobinding, along with the relatively uncommon GM1-phosphatidylethanolamine (PE) and cholesterol-PE synergy. In contrast, certain lipid pairs such as cholesterol and sphingomyelin appear to be in competition at several sites, despite their coexistence in lipid nanodomains. In addition, we observed intralipid competition between two lipid tails, with the receptor exhibiting increased interactions with the unsaturated lipid tails. We believe our work represents an important step in understanding the diversity of GPCR-lipid interactions and exploring synergistic cobinding and competition in natural membranes.

Structural dynamics of Smoothened (SMO) in the ciliary membrane and its interaction with membrane lipids

The Smoothened receptor (SMO, a 7 pass transmembrane domain, Class F GPCR family protein) plays a crucial role in the Hedgehog (HH) signaling pathway, which is involved in embryonic development and is implicated in various types of cancer throughout the animal kingdom. In the absence of HH signaling, SMO is inhibited by Patched 1 (PTC1; a 12 pass transmembrane domain protein), which is localized in the primary cilia. HH binding leads to the dislocation of PTC1 from the cilia, thus making way for SMO to localize in the primary cilia, as an essential prerequisite for its activation. We have carried out MARTINI coarse-grained molecular dynamics simulations of SMO in POPC and in ciliary membrane models, respectively, to study the interactions of SMO with cholesterol and other lipid molecules in the ciliary membrane, and to gain molecular-level insights into the role of the primary cilia in shaping the functional dynamics of SMO. We are able to identify the interaction of membrane cholesterols with definite sites and domains within SMO and relate them with known cholesterol-binding sequence and structure motifs. We show that cholesterol interactions with the transmembrane domain TMD, unlike those with the cysteine-rich domain (CRD) and the intracellular domain (ICD), are through residues belonging to known cholesterol-binding motifs. Notably, a few persistent interactions of cholesterol with lower TM cholesterol-binding domains are governed by the presence of multiple cholesterol-binding motifs. These analyses have also helped to identify and define a strict cholesterol consensus motif (CCM), which may well steer cholesterol into the hitherto identified binding sites within the TMD of SMO. We have also reported the interaction of phosphatidylinositol 4-phosphate with the intracellular region of transmembrane (TM) helices (TM1, TM3, TM4, and TM5), intracellular loop1, helix8, and Arg/Lys clusters of the ICD. Structural analysis of SMO domains shows significant changes in the CRD and ICD, during the course of the simulation. Further detailed analysis of the dynamics of the TMD reveals the movements of TM5, TM6, and TM7, linked with the helix8, which are possibly involved in shaping the conformational disposition of the ICD. The movement of these TM helices could possibly be a consequence of interactions involving the extracellular domain and extracellular loops. In addition, our analysis also shows that phosphatidylinositol-4-phosphate (PI4P), along with some ICD cholesterols, are implicated in anchoring SMO in the membrane.

Diacylglycerols interact with the L2 lipidation site in TRPC3 to induce a sensitized channel state

Coordination of lipids within transient receptor potential canonical channels (TRPCs) is essential for their Ca2+ signaling function. Single particle cryo-EM studies identified two lipid interaction sites, designated L1 and L2, which are proposed to accommodate diacylglycerols (DAGs). To explore the role of L1 and L2 in TRPC3 function, we combined structure-guided mutagenesis and electrophysiological recording with molecular dynamics (MD) simulations. MD simulations indicate rapid DAG accumulation within both L1 and L2 upon its availability within the plasma membrane. Electrophysiological experiments using a photoswitchable DAG-probe reveal potentiation of TRPC3 currents during repetitive activation by DAG. Importantly, initial DAG exposure generates a subsequently sensitized channel state that is associated with significantly faster activation kinetics. TRPC3 sensitization is specifically promoted by mutations within L2, with G652A exhibiting sensitization at very low levels of active DAG. We demonstrate the ability of TRPC3 to adopt a closed state conformation that features partial lipidation of L2 sites by DAG and enables fast activation of the channel by the phospholipase C-DAG pathway.

The role of cholesterol binding in the control of cholesterol by the Scap–Insig system

Scap and Insig, two proteins embedded in the membrane of the endoplasmic reticulum (ER), regulate the synthesis of cholesterol in animal cells by forming a dimer in the presence of high concentrations of cholesterol. Cryo-electron microscopic structures for the Scap–Insig dimer show a sterol-binding site at the dimer interface, but none of the structures include cholesterol itself. Here, a molecular docking approach developed to characterise cholesterol binding to the transmembrane (TM) regions of membrane proteins is used to characterise cholesterol binding to sites on the TM surface of the dimer and to the interfacial binding site. Binding of cholesterol is also observed at sites on the extra-membranous luminal domains of Scap, but the properties of these sites suggest that they will be unoccupied in vivo. Comparing the structure of Scap in the dimer with that predicted by AlphaFold for monomeric Scap suggests that dimer formation could result in relocation of TM helix 7 of Scap and of the loop between TM6 and 7, and that this could be the key change on Scap that signals that there is a high concentration of cholesterol in the ER.

Exploring TRPC3 Interaction with Cholesterol through Coarse-Grained Molecular Dynamics Simulations

Transient receptor potential canonical 3 (TRPC3) channel belongs to the superfamily of transient receptor potential (TRP) channels which mediate Ca2+ influx into the cell. These channels constitute essential elements of cellular signalling and have been implicated in a wide range of diseases. TRPC3 is primarily gated by lipids and its surface expression has been shown to be dependent on cholesterol, yet a comprehensive exploration of its interaction with this lipid has thus far not emerged. Here, through 80 µs of coarse-grained molecular dynamics simulations, we show that cholesterol interacts with multiple elements of the transmembrane machinery of TRPC3. Through our approach, we identify an annular binding site for cholesterol on the pre-S1 helix and a non-annular site at the interface between the voltage-sensor-like domain and pore domains. Here, cholesterol interacts with exposed polar residues and possibly acts to stabilise the domain interface.

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.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

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.

Regulation of Cholesterol Binding to the Receptor Patched1 by its interactions With the Ligand Sonic Hedgehog (Shh)

The Hedgehog (Hh) signaling pathway is essential in cell development and regeneration, which is activated by the ligand Sonic hedgehog (Shh). The binding of Shh to its receptor Patched1 (PTCH1) releases the inhibitory effect on the downstream protein Smoothened (SMO), a G-protein-coupled-receptor (GPCR) protein. Cholesterol was supposed to function as a secondary messenger between PTCH1 and SMO. However, the molecular mechanism of this regulation process is still unclear. Therefore, microsecond coarse-grained molecular dynamics simulations were performed to investigate the protein-lipid interactions of the PTCH1 monomer and dimer-Shh complex. It was observed that the binding of cholesterols to the monomer is more stable than that to the dimer-Shh complex. It is regulated by the enrichment of Ganglioside lipids around proteins and the conformation of Y446, a residue in the sterol-sensing domain (SSD). The regulation of Shh on the dynamics of PTCH1 was further analyzed to explore the allosteric communication pathways between the Shh and the SSD. Our study provides structural and dynamic details of an additional perspective on the regulation of Hh signaling pathway through the lipid micro-environments of PTCH1.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

NanoBRET and NanoBiT/BRET-Based Ligand Binding Assays Permit Quantitative Assessment of Small Molecule Ligand Binding to Smoothened

Smoothened (SMO) is a G protein-coupled receptor (GPCR) that mediates Hedgehog (Hh) signaling. SMO activity is regulated following the binding of Hh to the transmembrane protein Patched. Overactive SMO signaling is oncogenic, and hence this receptor is a target for several marketed drugs. However, development of new SMO ligands has been hampered by the fact that current radioligand and fluorescence-based binding assays are not high-throughput scalable. Here, we demonstrate two Nanoluciferase (Nluc) bioluminescence resonance energy transfer-based ligand binding assays (NanoBRET and NanoBiT/BRET) which provide a sensitive and high-throughput-compatible tool in drug screening efforts. In the described assays, SMO is N-terminally tagged either with full-length nanoluciferase or the partial HiBiT sequence, and subsequently binding of BODIPY-cyclopamine is assessed by quantifying resonance energy transfer between the receptor and the fluorescent ligand. Additionally, the assay allows performing competition binding experiments using commercially available SMO ligands, such as the SMO agonist SAG1.3.

Mechanisms of Smoothened Regulation in Hedgehog Signaling

The seven-transmembrane protein, Smoothened (SMO), has shown to be critical for the hedgehog (HH) signal transduction on the cell membrane (and the cilium in vertebrates). SMO is subjected to multiple types of post-translational regulations, including phosphorylation, ubiquitination, and sumoylation, which alter SMO intracellular trafficking and cell surface accumulation. Recently, SMO is also shown to be regulated by small molecules, such as oxysterol, cholesterol, and phospholipid. The activity of SMO must be very well balanced by these different mechanisms in vivo because the malfunction of SMO will not only cause developmental defects in early stages, but also induce cancers in late stages. Here, we discuss the activation and inactivation of SMO by different mechanisms to better understand how SMO is regulated by the graded HH signaling activity that eventually governs distinct development outcomes.

ProLint: a web-based framework for the automated data analysis and visualization of lipid-protein interactions

The functional activity of membrane proteins is carried out in a complex lipid environment. Increasingly, it is becoming clear that lipids are an important player in regulating or generally modulating their activity. A routinely used method to gain insight into this interplay between lipids and proteins are Molecular Dynamics (MD) simulations, since they allow us to study interactions at atomic or near-atomic detail as a function of time. A major bottleneck, however, is analyzing and visualizing lipid-protein interactions, which, in practice, is a time-demanding task. Here, we present ProLint (www.prolint.ca), a webserver that completely automates analysis of MD generated files and visualization of lipid-protein interactions. Analysis is modular allowing users to select their preferred method, and visualization is entirely interactive through custom built applications that enable a detailed qualitative and quantitative exploration of lipid-protein interactions. ProLint also includes a database of published MD results that have been processed through the ProLint workflow and can be visualized by anyone regardless of their level of experience with MD. The automated analysis, feature-rich visualization, database integration, and open-source distribution with an easy to install process, will allow ProLint to become a routine workflow in lipid-protein interaction studies.

The structural biology of canonical Wnt signalling

The Wnt signalling pathways are of great importance in embryonic development and oncogenesis. Canonical and non-canonical Wnt signalling pathways are known, with the canonical (or β-catenin dependent) pathway being perhaps the best studied of these. While structural knowledge of proteins and interactions involved in canonical Wnt signalling has accumulated over the past 20 years, the pace of discovery has increased in recent years, with the structures of several key proteins and assemblies in the pathway being released. In this review, we provide a brief overview of canonical Wnt signalling, followed by a comprehensive overview of currently available X-ray, NMR and cryoEM data elaborating the structures of proteins and interactions involved in canonical Wnt signalling. While the volume of structures available is considerable, numerous gaps in knowledge remain, particularly a comprehensive understanding of the assembly of large multiprotein complexes mediating key aspects of pathway, as well as understanding the structure and activation of membrane receptors in the pathway. Nonetheless, the presently available data affords considerable opportunities for structure-based drug design efforts targeting canonical Wnt signalling.

Interfacial binding sites for cholesterol on GABAA receptors and competition with neurosteroids

γ-Aminobutyric acid type A (GABA_A) receptors in the brain are located in the outer membranes of brain cells where the concentration of cholesterol is high. Of the 25 available high-resolution structures available for GABA_A receptors, none were determined in the presence of cholesterol, but four include resolved molecules of cholesterol hemisuccinate (CHS). Here, a molecular docking procedure is used to sweep the transmembrane (TM) surfaces of the receptors for cholesterol binding sites. Cholesterol docking poses determined in this way match 89% of the resolved CHS when CHS molecules deemed unlikely to represent typical bound cholesterols are excluded. The receptors are pentameric, and their TM surfaces consist of a set of five facets, each including pairs of TM helices from two adjacent subunits. Each facet contains hydrophobic hollows running from one side of the membrane to the other, within which are six potential binding sites for cholesterol, three on each side of the membrane. High-resolution structures of GABA_A receptors with bound neurosteroids show that neurosteroids bind in these cholesterol binding sites, so the binding of neurosteroids and cholesterol will be competitive.

Statins repress hedgehog signaling in medulloblastoma with no bone toxicities

The Hedgehog (Hh) pathway plays an indispensable role in bone development and genetic activation of the pathway results in medulloblastoma (MB), the most common malignant brain tumor in children. Inhibitors of Hh pathway (such as vismodegib and sonedigib), which are used to treat MB, cause irreversible defects in bone growth in young children. Cholesterol is required for the activation of the Hh pathway, and statins, inhibitors of cholesterol biosynthesis, suppress MB growth by repressing Hh signaling in tumor cells. Here, we investigate the role of cholesterol biosynthesis in the proliferation and Hh signaling in chondrocytes, and examine the bone development in mice after statin treatment. Statins significantly inhibited MB growth in young mice, but caused no defects in bone development. Conditional deletion of NADP steroid dehydrogenase-like (NSDHL), an enzyme necessary for cholesterol biosynthesis, suppressed cholesterol synthesis in chondrocytes, and disrupted the growth plate in mouse femur and tibia, indicating the important function of intracellular cholesterol in bone development. Hh pathway activation and the proliferation of chondrocytes were inhibited by statin treatment in vitro; however, statins did not impair bone growth in vivo due to insufficient penetration into the bone. Our studies reveal a critical role of cholesterol in bone development, and support the utilization of statins for treatment of MB as well as other Hh pathway-associated malignancies.

The ligand-bound state of a G protein-coupled receptor stabilizes the interaction of functional cholesterol molecules

Cholesterol is a major component of mammalian plasma membranes that not only affects the physical properties of the lipid bilayer but also is the function of many membrane proteins including G protein-coupled receptors. The oxytocin receptor (OXTR) is involved in parturition and lactation of mammals and in their emotional and social behaviors. Cholesterol acts on OXTR as an allosteric modulator inducing a high-affinity state for orthosteric ligands through a molecular mechanism that has yet to be determined. Using the ion channel-coupled receptor technology, we developed a functional assay of cholesterol modulation of G protein-coupled receptors that is independent of intracellular signaling pathways and operational in living cells. Using this assay, we discovered a stable binding of cholesterol molecules to the receptor when it adopts an orthosteric ligand-bound state. This stable interaction preserves the cholesterol-dependent activity of the receptor in cholesterol-depleted membranes. This mechanism was confirmed using time-resolved FRET experiments on WT OXTR expressed in CHO cells. Consequently, a positive cross-regulation sequentially occurs in OXTR between cholesterol and orthosteric ligands.

Plasma Membrane Integrates Biophysical and Biochemical Regulation to Trigger Immune Receptor Functions

Plasma membrane provides a biophysical and biochemical platform for immune cells to trigger signaling cascades and immune responses against attacks from foreign pathogens or tumor cells. Mounting evidence suggests that the biophysical-chemical properties of this platform, including complex compositions of lipids and cholesterols, membrane tension, and electrical potential, could cooperatively regulate the immune receptor functions. However, the molecular mechanism is still unclear because of the tremendous compositional complexity and spatio-temporal dynamics of the plasma membrane. Here, we review the recent significant progress of dynamical regulation of plasma membrane on immune receptors, including T cell receptor, B cell receptor, Fc receptor, and other important immune receptors, to proceed mechano-chemical sensing and transmembrane signal transduction. We also discuss how biophysical-chemical cues couple together to dynamically tune the receptor’s structural conformation or orientation, distribution, and organization, thereby possibly impacting their in-situ ligand binding and related signal transduction. Moreover, we propose that electrical potential could potentially induce the biophysical-chemical coupling change, such as lipid distribution and membrane tension, to inevitably regulate immune receptor activation.

Understanding the complex pharmacology of cannabidiol: Mounting evidence suggests a common binding site with cholesterol

Cannabidiol is claimed to bind to a large number of protein targets based on in vitro assays. This suggests opportunities for a wide range of therapeutic applications. On the other hand, the existence of phytochemical 'nuisance compounds' suggests some measure of caution - these compounds are capable of altering membrane biophysical properties and changing protein function without directly contacting a binding site. Like cannabidiol, cholesterol alters membrane properties, but it also binds directly to membrane proteins through abundant cholesterol recognition sites. We present the evidence that cannabidiol and cholesterol may bind to the same site on some proteins. As a starting point for further research, we also used blind docking to show that cannabidiol binds to a cholesterol binding site on the CB1 receptor. Elucidation of the mechanism(s) of action of cannabidiol will assist the prioritisation of in vitro hits across targets, improve the success rate of medicinal chemistry campaigns, and ultimately benefit patient populations by focusing resources on programs with the most translational potential.

Activation barriers in Class F G protein-coupled receptors revealed by umbrella sampling simulations

The Class F G protein-coupled receptors (GPCRs) include Smoothened and the ten Frizzled receptors, which are major cell membrane receptors in the Hedgehog and Wnt signalling pathways respectively and of enormous interest in embryonic development and as therapeutic targets in cancer. Recent crystal structures of Smoothened provide the opportunity to investigate the structural biology of Class F GPCRs in more detail, in turn, informing the development of therapeutics. A key question in this area is how one receptor may trigger distinct pathways - particularly relevant for Wnt signalling, in which signals may be transduced from a Frizzled via Dishevelled or G proteins, depending on the context. In this study, we employ adiabatic biased molecular dynamics and umbrella sampling to investigate the activation of Smoothened and Frizzled-7 in both the native state and bound to endogenous ligands, as well as how the clinically used Smoothened antagonist vismodegib alters this signalling. The results highlight key energetic barriers in the activation of these receptors, and the molecular features of the receptors mediating these barriers, demonstrating our approach as a robust means of investigating signalling through these receptors.

Data analysis of molecular dynamics simulation trajectories of β-sitosterol, sonidegib and cholesterol in smoothened protein with the CHARMM36 force field

Inactivation of smoothened protein (SMO) by the antagonists in SHH-driven cancer types is essential for inhibition of cancer progression. This article presents molecular dynamics (MD) trajectories of water solution of three protein-ligand complexes smoothened-β-sitosterol (SMO-BST), smoothened-sonidegib (SMO-SNG) and smoothened-cholesterol (SMO-CLR) using CHARMM36 and SPC/E water model combination. Additionally, the work presents the topologies and trajectories of GROMACS files that were employed to analyse the protein-ligand interaction types (PyContact) and binding energy calculation (g_mmpbsa). The data demonstrated that equilibrated models of SMO-SNG and SMO-CLR complexes showed crucial residues that almost similar for interaction and contribution energy as previously reported in laboratory setup (in vitro). Initial simulations confirmed the role of ARG451 and TRP535 in the dynamic regulation of SMO. These data then were used as a reference for understanding the molecular dynamics of SMO-BST complex and thus predicted its mechanism of action.

Impact of Cholesterol Concentration and Lipid Phase on Structure and Fluctuation of Amyloid Precursor Protein

Elevated levels of cellular cholesterol have been identified as one factor contributing to the onset of Alzheimer's disease (AD). Specific interaction between cholesterol and the amyloid precursor protein (APP), investigated via NMR experiments and computational studies, has been proposed to play a critical role in the processing of APP by secretases and the biogenesis of amyloid-β (Aβ) protein. We present all-atom molecular dynamics simulations of the 40-residue congener of the C-terminal domain of APP, C9916-55 (C99), in cholesterol-enriched DMPC lipid bilayers. We investigated the effect of cholesterol concentration on the conformational ensemble of wild-type C99 and C99-cholesterol associations at the low pH of endosomal environments, at which residues E22 and D23 are neutral. C99 was also characterized in liquid ordered domains for Dutch (E22Q) and Iowa (D23N) Familial AD mutants at low pH and for the wild-type sequence using protonation states characteristic of neutral pH. Our results reproduce the equilibrium constant of past NMR characterizations of the C99-cholesterol interaction but are not consistent with the C99-cholesterol binding hypothesis. We find that the lifetimes of both DMPC and cholesterol complexed with C99 display a power-law distribution of residence lifetimes. Longer-lived C99-DMPC and C99-cholesterol complexes are primarily stabilized by salt bridges and hydrogen bonds of lysine amines to phosphate and hydroxyl groups. Nevertheless, specific interfaces for C99-cholesterol association which are not present for DMPC can be identified. Changes to C99-cholesterol interfaces are found to depend on C99 tilt angle and orientation of the juxtamembrane domain of C99 containing residues E22 and D23. These observations support a more nuanced view of the C99-cholesterol interaction than has previously been suggested. We propose that cholesterol modulates the conformation and activity of C99 and other small transmembrane proteins indirectly through induction of the liquid ordered phase and directly through hydrogen bonding. This suggests a critical role for membrane heterogeneity introduced by cholesterol in modulating the structural ensemble of C99 and the production of Aβ.

24-Dehydrocholesterol reductase promotes the growth of breast cancer stem-like cells through the Hedgehog pathway

Cholesterol is a risk factor for breast cancer. However, it is still unclear whether the cholesterol biosynthesis pathway plays any significant role in breast carcinogenesis. 24-Dehydrocholesterol reductase (DHCR24) is a key enzyme in the cholesterol synthesis pathway. Although DHCR24 is reported to have different functions in different cancers, it is not clear whether DHCR24 is involved in breast cancer. In this study, we found that DHCR24 expression was higher in breast cancer especially in luminal and HER2 positive breast cancer tissues compared with normal breast. Changes in DHCR24 expression altered cellular cholesterol content without affecting the adherent growth of breast cancer cells. However, DHCR24 knockdown reduced whereas DHCR24 overexpression enhanced breast cancer stem-like cell populations such as mammosphere and aldehyde dehydrogenase positive cell numbers. In addition, DHCR24 overexpression increased the expression of the Hedgehog pathway-regulated genes. Treating DHCR24 overexpressing breast cancer cell lines with the Hedgehog pathway inhibitor GANT61 blocked DHCR24-induced mammosphere growth and increased mRNA levels of the Hedgehog regulated genes. Furthermore, expression of a constitutively activated mutant of Smoothened, a key hedgehog signal transducer, rescued the decreases in mammosphere growth and Hedgehog regulated gene expression induced by knockdown of DHCR24. These results indicate that DHCR24 promotes the growth of breast cancer stem-like cells in part through enhancing the Hedgehog signaling pathway. Our data suggest that cholesterol contribute to breast carcinogenesis by enhancing Hedgehog signaling and cancer stem-like cell populations. Enzymes including DHCR24 involved in cholesterol biosynthesis should be considered as potential treatment targets for breast cancer.

Interfacial Binding Sites for Cholesterol on Kir, Kv, K2P, and Related Potassium Channels

Inwardly rectifying, voltage-gated, two-pore domain, and related K+ channels are located in eukaryotic membranes rich in cholesterol. Here, molecular docking is used to detect specific binding sites ("hot spots") for cholesterol on K+ channels with characteristics that match those of known cholesterol binding sites. The transmembrane surfaces of all available high-resolution structures for K+ channels were swept for potential binding sites. Cholesterol poses were found to be located largely in hollows between protein ridges. A comparison between cholesterol poses and resolved phospholipids suggests that not all cholesterol molecules binding to the transmembrane surface of a K+ channel will result in displacement of a phospholipid molecule from the surface. Competition between cholesterol binding and binding of anionic phospholipids essential for activity could explain some of the effects of cholesterol on channel function.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122–E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP₂) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP₂ compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP₂ on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.

Mechanistic Insights into the Generation and Transduction of Hedgehog Signaling

Cell differentiation and proliferation require Hedgehog (HH) signaling and aberrant HH signaling causes birth defects or cancers. In this signaling pathway, the N-terminally palmitoylated and C-terminally cholesterylated HH ligand is secreted into the extracellular space with help of the Dispatched-1 (DISP1) and Scube2 proteins. The Patched-1 (PTCH1) protein releases its inhibition of the oncoprotein Smoothened (SMO) after binding the HH ligand, triggering downstream signaling events. In this review, we discuss the recent structural and biochemical studies on four major components of the HH pathway: the HH ligand, DISP1, PTCH1, and SMO. This research provides mechanistic insights into how HH signaling is generated and transduced from the cell surface into the intercellular space and will aid in facilitating the treatment of HH-related diseases.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels

Protein-lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP₂) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP₂; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein-lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP₂ interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP₂ concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP₂ is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein-lipid interactions.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels

Ion channels form pores that allow for the selective transport of ions across cell membranes, generating electrical signals in response to a variety of signals. Inward rectifier potassium (Kir) channels in particular are regulated by direct interactions with the complex mixture of lipids that are present in eukaryotic cell membranes. However, the molecular details of these concurrent lipid interactions with Kir channels are not clear and difficult to access via experimental methods. Here, we simulate the Kir2.2 channel in a complex lipid mixture to explore how anionic phospholipids and cholesterol dynamically organize around the membrane protein. In particular we demonstrate a synergy between binding interactions of different anionic phospholipid species which are known to activate Kir channels.

Protein–lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP₂) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP₂; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein–lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP₂ interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP₂ concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP₂ is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein–lipid interactions.

Inserting Small Molecules across Membrane Mixtures: Insight from the Potential of Mean Force

Small solutes have been shown to alter the lateral organization of cell membranes and reconstituted phospholipid bilayers; however, the mechanisms by which these changes happen are still largely unknown. Traditionally, both experiment and simulation studies have been restricted to testing only a few compounds at a time, failing to identify general molecular descriptors or chemical properties that would allow extrapolating beyond the subset of considered solutes. In this work, we probe the competing energetics of inserting a solute in different membrane environments by means of the potential of mean force. We show that these calculations can be used as a computationally efficient proxy to establish whether a solute will stabilize or destabilize domain phase separation. Combined with umbrella-sampling simulations and coarse-grained molecular dynamics simulations, we are able to screen solutes across a wide range of chemistries and polarities. Our results indicate that for the system under consideration, preferential partitioning and therefore effectiveness in altering membrane phase separation are strictly linked to the location of insertion in the bilayer (i.e., midplane or interface). Our approach represents a fast and simple tool for obtaining structural and thermodynamic insight into the partitioning of small molecules between lipid domains and its relation to phase separation, ultimately providing a platform for identifying the key determinants of this process.

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.