Exploring GPCR signaling pathway networks as cancer therapeutic targets

GPCR signaling can contribute to establishing the tumor microenvironment and influence the progression and metabolism of tumors. Arora et al.1 describe a systems-level approach to investigate the patterns of co-expression of GPCR signaling pathway networks across diverse tumors and identify network components that correlate with patient-survival data across different cancer types.

Intrinsic disorder: A term to define the specific physicochemical characteristic of protein conformational heterogeneity

In his commentary in this issue of Molecular Cell,1 Struhl reasons that the term "intrinsically disordered regions" represents a vague and confusing concept for protein function. However, the term "intrinsically disordered" highlights the important physicochemical characteristic of conformational heterogeneity. Thus, "intrinsically disordered" is the counterpart to the term "folded, " with neither term having specific functional implications.

Modular chimeric cytokine receptors with leucine zippers enhance the antitumour activity of CAR T cells via JAK/STAT signalling

The limited availability of cytokines in solid tumours hinders maintenance of the antitumour activity of chimeric antigen receptor (CAR) T cells. Cytokine receptor signalling pathways in CAR T cells can be activated by transgenic expression or injection of cytokines in the tumour, or by engineering the activation of cognate cytokine receptors. However, these strategies are constrained by toxicity arising from the activation of bystander cells, by the suboptimal biodistribution of the cytokines and by downregulation of the cognate receptor. Here we show that replacement of the extracellular domains of heterodimeric cytokine receptors in T cells with two leucine zipper motifs provides optimal Janus kinase/signal transducer and activator of transcription signalling. Such chimeric cytokine receptors, which can be generated for common γ-chain receptors, interleukin-10 and -12 receptors, enabled T cells to survive cytokine starvation without induction of autonomous cell growth, and augmented the effector function of CAR T cells in vitro in the setting of chronic antigen exposure and in human tumour xenografts in mice. As a modular design, leucine zippers can be used to generate constitutively active cytokine receptors in effector immune cells.

Peptide-scFv antigen recognition domains effectively confer CAR T cell multiantigen specificity

The emergence of immune escape is a significant roadblock to developing effective chimeric antigen receptor (CAR) T cell therapies against hematological malignancies, including acute myeloid leukemia (AML). Here, we demonstrate feasibility of targeting two antigens simultaneously by combining a GRP78-specific peptide antigen recognition domain with a CD123-specific scFv to generate a peptide-scFv bispecific antigen recognition domain (78.123). To achieve this, we test linkers with varying length and flexibility and perform immunophenotypic and functional characterization. We demonstrate that bispecific CAR T cells successfully recognize and kill tumor cells that express GRP78, CD123, or both antigens and have improved antitumor activity compared to their monospecific counterparts when both antigens are expressed. Protein structure prediction suggests that linker length and compactness influence the functionality of the generated bispecific CARs. Thus, we present a bispecific CAR design strategy to prevent immune escape in AML that can be extended to other peptide-scFv combinations.

Dynamic conformational changes of a tardigrade group-3 late embryogenesis abundant protein modulate membrane biophysical properties

A number of intrinsically disordered proteins (IDPs) encoded in stress-tolerant organisms, such as tardigrade, can confer fitness advantage and abiotic stress tolerance when heterologously expressed. Tardigrade-specific disordered proteins including the cytosolic-abundant heat-soluble proteins are proposed to confer stress tolerance through vitrification or gelation, whereas evolutionarily conserved IDPs in tardigrades may contribute to stress tolerance through other biophysical mechanisms. In this study, we characterized the mechanism of action of an evolutionarily conserved, tardigrade IDP, HeLEA1, which belongs to the group-3 late embryogenesis abundant (LEA) protein family. HeLEA1 homologs are found across different kingdoms of life. HeLEA1 is intrinsically disordered in solution but shows a propensity for helical structure across its entire sequence. HeLEA1 interacts with negatively charged membranes via dynamic disorder-to-helical transition, mainly driven by electrostatic interactions. Membrane interaction of HeLEA1 is shown to ameliorate excess surface tension and lipid packing defects. HeLEA1 localizes to the mitochondrial matrix when expressed in yeast and interacts with model membranes mimicking inner mitochondrial membrane. Yeast expressing HeLEA1 shows enhanced tolerance to hyperosmotic stress under nonfermentative growth and increased mitochondrial membrane potential. Evolutionary analysis suggests that although HeLEA1 homologs have diverged their sequences to localize to different subcellular organelles, all homologs maintain a weak hydrophobic moment that is characteristic of weak and reversible membrane interaction. We suggest that such dynamic and weak protein-membrane interaction buffering alterations in lipid packing could be a conserved strategy for regulating membrane properties and represent a general biophysical solution for stress tolerance across the domains of life.

Molecular determinants of ligand efficacy and potency in GPCR signaling

Heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs) bind to extracellular ligands and drugs and modulate intracellular responses through conformational changes. Despite their importance as drug targets, the molecular origins of pharmacological properties such as efficacy (maximum signaling response) and potency (the ligand concentration at half-maximal response) remain poorly understood for any ligand-receptor-signaling system. We used the prototypical adrenaline-β2 adrenergic receptor-G protein system to reveal how specific receptor residues decode and translate the information encoded in a ligand to mediate a signaling response. We present a data science framework to integrate pharmacological and structural data to uncover structural changes and allosteric networks relevant for ligand pharmacology. These methods can be tailored to study any ligand-receptor-signaling system, and the principles open possibilities for designing orthosteric and allosteric compounds with defined signaling properties.

Defining the condensate landscape of fusion oncoproteins

Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.

Predicting evolutionary outcomes through the probability of accessing sequence variants

Natural selection can only operate on available genetic variation. Thus, determining the probability of accessing different sequence variants from a starting sequence can help predict evolutionary trajectories and outcomes. We define the concept of "variant accessibility" as the probability that a set of genotypes encoding a particular protein function will arise through mutations before subject to natural selection. This probability is shaped by the mutational biases of nucleotides and the structure of the genetic code. Using the influenza A virus as a model, we discuss how a more accessible but less fit variant can emerge as an adaptation rather than a more fit variant. We describe a genotype-accessibility landscape, complementary to the genotype-fitness landscape, that informs the likelihood of a starting sequence reaching different parts of genotype space. The proposed framework lays the foundation for predicting the emergence of adaptive genotypes in evolving systems such as viruses and tumors.

Recent breakthroughs in computational structural biology harnessing the power of sequences and structures

Recent advances in computational approaches and their integration into structural biology enable tackling increasingly complex questions. Here, we discuss several key areas, highlighting breakthroughs and remaining challenges. Theoretical modeling has provided tools to accurately predict and design protein structures on a scale currently difficult to achieve using experimental approaches. Molecular Dynamics simulations have become faster and more precise, delivering actionable information inaccessible by current experimental methods. Virtual screening workflows allow a high-throughput approach to discover ligands that bind and modulate protein function, while Machine Learning methods enable the design of proteins with new functionalities. Integrative structural biology combines several of these approaches, pushing the frontiers of structural and functional characterization to ever larger systems, advancing towards a complete understanding of the living cell. These breakthroughs will accelerate and significantly impact diverse areas of science.

Data-driven design of orthogonal protein-protein interactions

Engineering protein-protein interactions to generate new functions presents a challenge with great potential for many applications, ranging from therapeutics to synthetic biology. To avoid unwanted cross-talk with preexisting protein interaction networks in a cell, the specificity and selectivity of newly engineered proteins must be controlled. Here, we developed a computational strategy that mimics gene duplication and the divergence of preexisting interacting protein pairs to design new interactions. We used the bacterial PhoQ-PhoP two-component system as a model system to demonstrate the feasibility of this strategy and validated the approach with known experimental results. The designed protein pairs are predicted to exclusively interact with each other and to be insulated from potential cross-talk with their native partners. Thus, our approach enables exploration of uncharted regions of the protein sequence space and the design of new interacting protein pairs.

Dynamic spatiotemporal determinants modulate GPCR:G protein coupling selectivity and promiscuity

Recent studies have shown that G protein coupled receptors (GPCRs) show selective and promiscuous coupling to different Gα protein subfamilies and yet the mechanisms of the range of coupling preferences remain unclear. Here, we use Molecular Dynamics (MD) simulations on ten GPCR:G protein complexes and show that the location (spatial) and duration (temporal) of intermolecular contacts at the GPCR:Gα protein interface play a critical role in how GPCRs selectively interact with G proteins. We identify that some GPCR:G protein interface contacts are common across Gα subfamilies and others specific to Gα subfamilies. Using large scale data analysis techniques on the MD simulation snapshots we derive a spatio-temporal code for contacts that confer G protein selective coupling and validated these contacts using G protein activation BRET assays. Our results demonstrate that promiscuous GPCRs show persistent sampling of the common contacts more than G protein specific contacts. These findings suggest that GPCRs maintain contact with G proteins through a common central interface, while the selectivity comes from G protein specific contacts at the periphery of the interface.

Discovery and Genetic Code Expansion of a Polyethylene Terephthalate (PET) Hydrolase from the Human Saliva Metagenome for the Degradation and Bio‐Functionalization of PET

We report a bioinformatic workflow and subsequent discovery of a new polyethylene terephthalate (PET) hydrolase, which we named MG8, from the human saliva metagenome. MG8 has robust PET plastic degradation activities under different temperature and salinity conditions, outperforming several naturally occurring and engineered hydrolases in degrading PET. Moreover, we genetically encoded 2,3-diaminopropionic acid (DAP) in place of the catalytic serine residue of MG8, thereby converting a PET hydrolase into a covalent binder for bio-functionalization of PET. We show that MG8(DAP), in conjunction with a split green fluorescent protein system, can be used to attach protein cargos to PET as well as other polyester plastics. The discovery of a highly active PET hydrolase from the human metagenome—currently an underexplored resource for industrial enzyme discovery—as well as the repurposing of such an enzyme into a plastic functionalization tool, should facilitate ongoing efforts to degrade and maximize reusability of PET.

Conformational selection guides β-arrestin recruitment at a biased G protein-coupled receptor

G protein-coupled receptors (GPCRs) recruit β-arrestins to coordinate diverse cellular processes, but the structural dynamics driving this process are poorly understood. Atypical chemokine receptors (ACKRs) are intrinsically biased GPCRs that engage β-arrestins but not G proteins, making them a model system for investigating the structural basis of β-arrestin recruitment. Here, we performed nuclear magnetic resonance (NMR) experiments on ¹³CH₃-ε-methionine-labeled ACKR3, revealing that β-arrestin recruitment is associated with conformational exchange at key regions of the extracellular ligand-binding pocket and intracellular β-arrestin-coupling region. NMR studies of ACKR3 mutants defective in β-arrestin recruitment identified an allosteric hub in the receptor core that coordinates transitions among heterogeneously populated and selected conformational states. Our data suggest that conformational selection guides β-arrestin recruitment by tuning receptor dynamics at intracellular and extracellular regions.

Integrated Genomic Analysis Identifies UBTF Tandem Duplications as a Recurrent Lesion in Pediatric Acute Myeloid Leukemia

The genetics of relapsed pediatric acute myeloid leukemia (AML) has yet to be comprehensively defined. Here, we present the spectrum of genomic alterations in 136 relapsed pediatric AMLs. We identified recurrent exon 13 tandem duplications (TD) in upstream binding transcription factor (UBTF) in 9% of relapsed AML cases. UBTF-TD AMLs commonly have normal karyotype or trisomy 8 with cooccurring WT1 mutations or FLT3-ITD but not other known oncogenic fusions. These UBTF-TD events are stable during disease progression and are present in the founding clone. In addition, we observed that UBTF-TD AMLs account for approximately 4% of all de novo pediatric AMLs, are less common in adults, and are associated with poor outcomes and MRD positivity. Expression of UBTF-TD in primary hematopoietic cells is sufficient to enhance serial clonogenic activity and to drive a similar transcriptional program to UBTF-TD AMLs. Collectively, these clinical, genomic, and functional data establish UBTF-TD as a new recurrent mutation in AML.

SIGNIFICANCE

We defined the spectrum of mutations in relapsed pediatric AML and identified UBTF-TDs as a new recurrent genetic alteration. These duplications are more common in children and define a group of AMLs with intermediate-risk cytogenetic abnormalities, FLT3-ITD and WT1 alterations, and are associated with poor outcomes. See related commentary by Hasserjian and Nardi, p. 173. This article is highlighted in the In This Issue feature, p. 171.

The chemotherapeutic CX-5461 primarily targets TOP2B and exhibits selective activity in high-risk neuroblastoma

Survival in high-risk pediatric neuroblastoma has remained around 50% for the last 20 years, with immunotherapies and targeted therapies having had minimal impact. Here, we identify the small molecule CX-5461 as selectively cytotoxic to high-risk neuroblastoma and synergistic with low picomolar concentrations of topoisomerase I inhibitors in improving survival in vivo in orthotopic patient-derived xenograft neuroblastoma mouse models. CX-5461 recently progressed through phase I clinical trial as a first-in-human inhibitor of RNA-POL I. However, we also use a comprehensive panel of in vitro and in vivo assays to demonstrate that CX-5461 has been mischaracterized and that its primary target at pharmacologically relevant concentrations, is in fact topoisomerase II beta (TOP2B), not RNA-POL I. This is important because existing clinically approved chemotherapeutics have well-documented off-target interactions with TOP2B, which have previously been shown to cause both therapy-induced leukemia and cardiotoxicity-often-fatal adverse events, which can emerge several years after treatment. Thus, while we show that combination therapies involving CX-5461 have promising anti-tumor activity in vivo in neuroblastoma, our identification of TOP2B as the primary target of CX-5461 indicates unexpected safety concerns that should be examined in ongoing phase II clinical trials in adult patients before pursuing clinical studies in children.

GPCR activation mechanisms across classes and macro/microscales

Two-thirds of human hormones and one-third of clinical drugs activate ~350 G-protein-coupled receptors (GPCR) belonging to four classes: A, B1, C and F. Whereas a model of activation has been described for class A, very little is known about the activation of the other classes, which differ by being activated by endogenous ligands bound mainly or entirely extracellularly. Here we show that, although they use the same structural scaffold and share several 'helix macroswitches', the GPCR classes differ in their 'residue microswitch' positions and contacts. We present molecular mechanistic maps of activation for each GPCR class and methods for contact analysis applicable for any functional determinants. This provides a superfamily residue-level rationale for conformational selection and allosteric communication by ligands and G proteins, laying the foundation for receptor-function studies and drugs with the desired modality.

Lipolysis drives expression of the constitutively active receptor GPR3 to induce adipose thermogenesis

Thermogenic adipocytes possess a therapeutically appealing, energy-expending capacity, which is canonically cold-induced by ligand-dependent activation of β-adrenergic G protein-coupled receptors (GPCRs). Here, we uncover an alternate paradigm of GPCR-mediated adipose thermogenesis through the constitutively active receptor, GPR3. We show that the N terminus of GPR3 confers intrinsic signaling activity, resulting in continuous Gs-coupling and cAMP production without an exogenous ligand. Thus, transcriptional induction of Gpr3 represents the regulatory parallel to ligand-binding of conventional GPCRs. Consequently, increasing Gpr3 expression in thermogenic adipocytes is alone sufficient to drive energy expenditure and counteract metabolic disease in mice. Gpr3 transcription is cold-stimulated by a lipolytic signal, and dietary fat potentiates GPR3-dependent thermogenesis to amplify the response to caloric excess. Moreover, we find GPR3 to be an essential, adrenergic-independent regulator of human brown adipocytes. Taken together, our findings reveal a noncanonical mechanism of GPCR control and thermogenic activation through the lipolysis-induced expression of constitutively active GPR3.

Biological impact of mutually exclusive exon switching

Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.

Human MC4R variants affect endocytosis, trafficking and dimerization revealing multiple cellular mechanisms involved in weight regulation

The Melanocortin-4 Receptor (MC4R) plays a pivotal role in energy homeostasis. We used human MC4R mutations associated with an increased or decreased risk of obesity to dissect mechanisms that regulate MC4R function. Most obesity-associated mutations impair trafficking to the plasma membrane (PM), whereas obesity-protecting mutations either accelerate recycling to the PM or decrease internalization, resulting in enhanced signaling. MC4R mutations that do not affect canonical Gα_s protein-mediated signaling, previously considered to be non-pathogenic, nonetheless disrupt agonist-induced internalization, β-arrestin recruitment, and/or coupling to Gα_s, establishing their causal role in severe obesity. Structural mapping reveals ligand-accessible sites by which MC4R couples to effectors and residues involved in the homodimerization of MC4R, which is disrupted by multiple obesity-associated mutations. Human genetic studies reveal that endocytosis, intracellular trafficking, and homodimerization regulate MC4R function to a level that is physiologically relevant, supporting the development of chaperones, agonists, and allosteric modulators of MC4R for weight loss therapy.

Evolution of fold switching in a metamorphic protein

Metamorphic proteins switch between different folds, defying the protein folding paradigm. It is unclear how fold switching arises during evolution. With ancestral reconstruction and nuclear magnetic resonance, we studied the evolution of the metamorphic human protein XCL1, which has two distinct folds with different functions, making it an unusual member of the chemokine family, whose members generally adopt one conserved fold. XCL1 evolved from an ancestor with the chemokine fold. Evolution of a dimer interface, changes in structural constraints and molecular strain, and alteration of intramolecular protein contacts drove the evolution of metamorphosis. Then, XCL1 likely evolved to preferentially populate the noncanonical fold before reaching its modern-day near-equal population of folds. These discoveries illuminate how one sequence has evolved to encode multiple structures, revealing principles for protein design and engineering.

Combinatorial multivalent interactions drive cooperative assembly of the COPII coat

Protein secretion is initiated at the endoplasmic reticulum by the COPII coat, which self-assembles to form vesicles. Here, we examine the mechanisms by which a cargo-bound inner coat layer recruits and is organized by an outer scaffolding layer to drive local assembly of a stable structure rigid enough to enforce membrane curvature. An intrinsically disordered region in the outer coat protein, Sec31, drives binding with an inner coat layer via multiple distinct interfaces, including a newly defined charge-based interaction. These interfaces combinatorially reinforce each other, suggesting coat oligomerization is driven by the cumulative effects of multivalent interactions. The Sec31 disordered region could be replaced by evolutionarily distant sequences, suggesting plasticity in the binding interfaces. Such a multimodal assembly platform provides an explanation for how cells build a powerful yet transient scaffold to direct vesicle traffic.

Combinatorial expression of GPCR isoforms affects signalling and drug responses

G-protein-coupled receptors (GPCRs) are membrane proteins that modulate physiology across human tissues in response to extracellular signals. GPCR-mediated signalling can differ because of changes in the sequence1,2 or expression3 of the receptors, leading to signalling bias when comparing diverse physiological systems4. An underexplored source of such bias is the generation of functionally diverse GPCR isoforms with different patterns of expression across different tissues. Here we integrate data from human tissue-level transcriptomes, GPCR sequences and structures, proteomics, single-cell transcriptomics, population-wide genetic association studies and pharmacological experiments. We show how a single GPCR gene can diversify into several isoforms with distinct signalling properties, and how unique isoform combinations expressed in different tissues can generate distinct signalling states. Depending on their structural changes and expression patterns, some of the detected isoforms may influence cellular responses to drugs and represent new targets for developing drugs with improved tissue selectivity. Our findings highlight the need to move from a canonical to a context-specific view of GPCR signalling that considers how combinatorial expression of isoforms in a particular cell type, tissue or organism collectively influences receptor signalling and drug responses.

Structural and functional characterization of G protein-coupled receptors with deep mutational scanning

The >800 human G protein–coupled receptors (GPCRs) are responsible for transducing diverse chemical stimuli to alter cell state- and are the largest class of drug targets. Their myriad structural conformations and various modes of signaling make it challenging to understand their structure and function. Here, we developed a platform to characterize large libraries of GPCR variants in human cell lines with a barcoded transcriptional reporter of G protein signal transduction. We tested 7800 of 7828 possible single amino acid substitutions to the beta-2 adrenergic receptor (β₂AR) at four concentrations of the agonist isoproterenol. We identified residues specifically important for β₂AR signaling, mutations in the human population that are potentially loss of function, and residues that modulate basal activity. Using unsupervised learning, we identify residues critical for signaling, including all major structural motifs and molecular interfaces. We also find a previously uncharacterized structural latch spanning the first two extracellular loops that is highly conserved across Class A GPCRs and is conformationally rigid in both the inactive and active states of the receptor. More broadly, by linking deep mutational scanning with engineered transcriptional reporters, we establish a generalizable method for exploring pharmacogenomics, structure and function across broad classes of drug receptors.

Context-Specific Striatal Astrocyte Molecular Responses Are Phenotypically Exploitable

Astrocytes tile the central nervous system and are widely implicated in brain diseases, but the molecular mechanisms by which astrocytes contribute to brain disorders remain incompletely explored. By performing astrocyte gene expression analyses following 14 experimental perturbations of relevance to the striatum, we discovered that striatal astrocytes mount context-specific molecular responses at the level of genes, pathways, and upstream regulators. Through data mining, we also identified astrocyte pathways in Huntington's disease (HD) that were reciprocally altered with respect to the activation of striatal astrocyte G protein-coupled receptor (GPCR) signaling. Furthermore, selective striatal astrocyte stimulation of the G_i-GPCR pathway in vivo corrected several HD-associated astrocytic, synaptic, and behavioral phenotypes, with accompanying improvement of HD-associated astrocyte signaling pathways, including those related to synaptogenesis and neuroimmune functions. Overall, our data show that astrocytes are malleable, using context-specific responses that can be dissected molecularly and used for phenotypic benefit in brain disorders.

A Global Map of G Protein Signaling Regulation by RGS Proteins

The control over the extent and timing of G protein signaling is provided by the regulator of G protein signaling (RGS) proteins that deactivate G protein α subunits (Gα). Mammalian genomes encode 20 canonical RGS and 16 Gα genes with key roles in physiology and disease. To understand the principles governing the selectivity of Gα regulation by RGS, we examine the catalytic activity of all canonical human RGS proteins and their selectivity for a complete set of Gα substrates using real-time kinetic measurements in living cells. The data reveal rules governing RGS-Gα recognition, the structural basis of its selectivity, and provide principles for engineering RGS proteins with defined selectivity. The study also explores the evolution of RGS-Gα selectivity through ancestral reconstruction and demonstrates how naturally occurring non-synonymous variants in RGS alter signaling. These results provide a blueprint for decoding signaling selectivity and advance our understanding of molecular recognition principles.

•

Systematic analysis reveals G protein selectivity of all canonical RGS proteins

•

RGS proteins rely on selectivity bar codes for selective G protein recognition

•

Transplantation of bar codes across RGS proteins switches their G protein preferences

•

Natural variants, mutations, and evolution shape RGS selectivity

Amino acid homorepeats in proteins

Amino acid homorepeats, or homorepeats, are polypeptide segments found in proteins that contain stretches of identical amino acid residues. Although abnormal homorepeat expansions are linked to pathologies such as neurodegenerative diseases, homorepeats are prevalent in eukaryotic proteomes, suggesting that they are important for normal physiology. In this Review, we discuss recent advances in our understanding of the biological functions of homorepeats, which range from facilitating subcellular protein localization to mediating interactions between proteins across diverse cellular pathways. We explore how the functional diversity of homorepeat-containing proteins could be linked to the ability of homorepeats to adopt different structural conformations, an ability influenced by repeat composition, repeat length and the nature of flanking sequences. We conclude by highlighting how an understanding of homorepeats will help us better characterize and develop therapeutics against the human diseases to which they contribute.

Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method

G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.

Pathway perturbations in signaling networks: Linking genotype to phenotype

Genes and gene products interact with each other to form signal transduction networks in the cell. The interactome networks are under intricate regulation in physiological conditions, but could go awry upon genome instability caused by genetic mutations. In the past decade with next-generation sequencing technologies, an increasing number of genomic mutations have been identified in a variety of disease patients and healthy individuals. As functional and systematic studies on these mutations leap forward, they begin to reveal insights into cellular homeostasis and disease mechanisms. In this review, we discuss recent advances in the field of network biology and signaling pathway perturbations upon genomic changes, and highlight the success of various omics datasets in unraveling genotype-to-phenotype relationships.

Common activation mechanism of class A GPCRs

Class A G-protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology. Understanding receptor activation mechanism is critical for discovering novel therapeutics since about one-third of all marketed drugs target members of this family. GPCR activation is an allosteric process that couples agonist binding to G-protein recruitment, with the hallmark outward movement of transmembrane helix 6 (TM6). However, what leads to TM6 movement and the key residue level changes of this movement remain less well understood. Here, we report a framework to quantify conformational changes. By analyzing the conformational changes in 234 structures from 45 class A GPCRs, we discovered a common GPCR activation pathway comprising of 34 residue pairs and 35 residues. The pathway unifies previous findings into a common activation mechanism and strings together the scattered key motifs such as CWxP, DRY, Na+ pocket, NPxxY and PIF, thereby directly linking the bottom of ligand-binding pocket with G-protein coupling region. Site-directed mutagenesis experiments support this proposition and reveal that rational mutations of residues in this pathway can be used to obtain receptors that are constitutively active or inactive. The common activation pathway provides the mechanistic interpretation of constitutively activating, inactivating and disease mutations. As a module responsible for activation, the common pathway allows for decoupling of the evolution of the ligand binding site and G-protein-binding region. Such an architecture might have facilitated GPCRs to emerge as a highly successful family of proteins for signal transduction in nature.

Molecular Deconvolution Platform to Establish Disease Mechanisms by Surveying GPCR Signaling

Despite the wealth of genetic information available, mechanisms underlying pathological effects of disease-associated mutations in components of G protein-coupled receptor (GPCR) signaling cascades remain elusive. In this study, we developed a scalable approach for the functional analysis of clinical variants in GPCR pathways along with a complete analytical framework. We applied the strategy to evaluate an extensive set of dystonia-causing mutations in G protein Gαolf. Our quantitative analysis revealed diverse mechanisms by which pathogenic variants disrupt GPCR signaling, leading to a mechanism-based classification of dystonia. In light of significant clinical heterogeneity, the mechanistic analysis of individual disease-associated variants permits tailoring personalized intervention strategies, which makes it superior to the current phenotype-based approach. We propose that the platform developed in this study can be universally applied to evaluate disease mechanisms for conditions associated with genetic variation in all components of GPCR signaling.

•

A scalable platform allows multidimensional analysis of GPCR signaling

•

The approach is applied to dystonia-causing mutations in G protein Gαolf

•

Pathogenic variants in Gαolf disrupt GPCR signaling by diverse mechanisms

•

Mechanism-based disease classification could allow targeted therapies

Madan Babu

Molecular Systems Biology warmly welcomes its new academic Chief Editor, M. Madan Babu. Madan shared his thoughts on the evolution of the field and the importance of bridging disciplines to enable next generation systems biology.

An intrinsically disordered proteins community for ELIXIR

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) are now recognised as major determinants in cellular regulation. This white paper presents a roadmap for future e-infrastructure developments in the field of IDP research within the ELIXIR framework. The goal of these developments is to drive the creation of high-quality tools and resources to support the identification, analysis and functional characterisation of IDPs. The roadmap is the result of a workshop titled “An intrinsically disordered protein user community proposal for ELIXIR” held at the University of Padua. The workshop, and further consultation with the members of the wider IDP community, identified the key priority areas for the roadmap including the development of standards for data annotation, storage and dissemination; integration of IDP data into the ELIXIR Core Data Resources; and the creation of benchmarking criteria for IDP-related software. Here, we discuss these areas of priority, how they can be implemented in cooperation with the ELIXIR platforms, and their connections to existing ELIXIR Communities and international consortia. The article provides a preliminary blueprint for an IDP Community in ELIXIR and is an appeal to identify and involve new stakeholders.

The fitness cost and benefit of phase-separated protein deposits

Phase separation of soluble proteins into insoluble deposits is associated with numerous diseases. However, protein deposits can also function as membrane-less compartments for many cellular processes. What are the fitness costs and benefits of forming such deposits in different conditions? Using a model protein that phase-separates into deposits, we distinguish and quantify the fitness contribution due to the loss or gain of protein function and deposit formation in yeast. The environmental condition and the cellular demand for the protein function emerge as key determinants of fitness. Protein deposit formation can influence cell-to-cell variation in free protein abundance between individuals of a cell population (i.e., gene expression noise). This results in variable manifestation of protein function and a continuous range of phenotypes in a cell population, favoring survival of some individuals in certain environments. Thus, protein deposit formation by phase separation might be a mechanism to sense protein concentration in cells and to generate phenotypic variability. The selectable phenotypic variability, previously described for prions, could be a general property of proteins that can form phase-separated assemblies and may influence cell fitness.

Molecular Signatures of Fusion Proteins in Cancer

Although gene fusions are recognized as driver mutations in a wide variety of cancers, the general molecular mechanisms underlying oncogenic fusion proteins are insufficiently understood. Here, we employ large-scale data integration and machine learning and (1) identify three functionally distinct subgroups of gene fusions and their molecular signatures; (2) characterize the cellular pathways rewired by fusion events across different cancers; and (3) analyze the relative importance of over 100 structural, functional, and regulatory features of ∼2200 gene fusions. We report subgroups of fusions that likely act as driver mutations and find that gene fusions disproportionately affect pathways regulating cellular shape and movement. Although fusion proteins are similar across different cancer types, they affect cancer type-specific pathways. Key indicators of fusion-forming proteins include high and nontissue specific expression, numerous splice sites, and higher centrality in protein-interaction networks. Together, these findings provide unifying and cancer type-specific trends across diverse oncogenic fusion proteins.

A conserved molecular switch in Class F receptors regulates receptor activation and pathway selection

Class F receptors are considered valuable therapeutic targets due to their role in human disease, but structural changes accompanying receptor activation remain unexplored. Employing population and cancer genomics data, structural analyses, molecular dynamics simulations, resonance energy transfer-based approaches and mutagenesis, we identify a conserved basic amino acid in TM6 in Class F receptors that acts as a molecular switch to mediate receptor activation. Across all tested Class F receptors (FZD_4,5,6,7, SMO), mutation of the molecular switch confers an increased potency of agonists by stabilizing an active conformation as assessed by engineered mini G proteins as conformational sensors. Disruption of the switch abrogates the functional interaction between FZDs and the phosphoprotein Dishevelled, supporting conformational selection as a prerequisite for functional selectivity. Our studies reveal the molecular basis of a common activation mechanism conserved in all Class F receptors, which facilitates assay development and future discovery of Class F receptor-targeting drugs.

Class F receptors are therapeutic targets in human disease and understanding their structural changes during receptor activation may provide important pharmacological insight. Here, the authors combine computational and experimental methods to identify a molecular switch in TM6/7 of Class F receptors that mediates receptor activation.

From prioritisation to understanding: mechanistic predictions of variant effects

The widespread application of sequencing technologies, used for example to obtain data from healthy individuals or patient cohorts, has led to the identification of numerous mutations, the effect of which remains largely unclear. Therefore, developing approaches allowing accurate in‐silico prediction of mutation effects is becoming increasingly important. In their recent study, Beltrao and colleagues (Wagih et al , 2018 ) describe an integrative approach for determining the effects of mutations from the perspective of protein structure, conservation and transcription factor binding. This allows for predicting the mechanisms underlying the most impactful variants rather than just identifying these variants.

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

Revealing the Determinants of Widespread Alternative Splicing Perturbation in Cancer

It is increasingly appreciated that alternative splicing plays a key role in generating functional specificity and diversity in cancer. However, the mechanisms by which cancer mutations perturb splicing remain unknown. Here, we developed a network-based strategy, DrAS-Net, to investigate more than 2.5 million variants across cancer types and link somatic mutations with cancer-specific splicing events. We identified more than 40,000 driver variant candidates and their 80,000 putative splicing targets deregulated in 33 cancer types and inferred their functional impact. Strikingly, tumors with splicing perturbations show reduced expression of immune system-related genes and increased expression of cell proliferation markers. Tumors harboring different mutations in the same gene often exhibit distinct splicing perturbations. Further stratification of 10,000 patients based on their mutation-splicing relationships identifies subtypes with distinct clinical features, including survival rates. Our work reveals how single-nucleotide changes can alter the repertoires of splicing isoforms, providing insights into oncogenic mechanisms for precision medicine.

Molecular mechanism of modulating arrestin conformation by GPCR phosphorylation

Arrestins regulate the signaling of ligand-activated, phosphorylated G-protein-coupled receptors (GPCRs). Different patterns of receptor phosphorylation (phosphorylation barcode) can modulate arrestin conformations, resulting in distinct functional outcomes (for example, desensitization, internalization, and downstream signaling). However, the mechanism of arrestin activation and how distinct receptor phosphorylation patterns could induce different conformational changes on arrestin are not fully understood. We analyzed how each arrestin amino acid contributes to its different conformational states. We identified a conserved structural motif that restricts the mobility of the arrestin finger loop in the inactive state and appears to be regulated by receptor phosphorylation. Distal and proximal receptor phosphorylation sites appear to selectively engage with distinct arrestin structural motifs (that is, micro-locks) to induce different arrestin conformations. These observations suggest a model in which different phosphorylation patterns of the GPCR C terminus can combinatorially modulate the conformation of the finger loop and other phosphorylation-sensitive structural elements to drive distinct arrestin conformation and functional outcomes.

High-throughput discovery of functional disordered regions: investigation of transactivation domains

Over 40% of proteins in any eukaryotic genome encode intrinsically disordered regions (IDRs) that do not adopt defined tertiary structures. Certain IDRs perform critical functions, but discovering them is non‐trivial as the biological context determines their function. We present IDR‐Screen, a framework to discover functional IDRs in a high‐throughput manner by simultaneously assaying large numbers of DNA sequences that code for short disordered sequences. Functionality‐conferring patterns in their protein sequence are inferred through statistical learning. Using yeast HSF1 transcription factor‐based assay, we discovered IDRs that function as transactivation domains (TADs) by screening a random sequence library and a designed library consisting of variants of 13 diverse TADs. Using machine learning, we find that segments devoid of positively charged residues but with redundant short sequence patterns of negatively charged and aromatic residues are a generic feature for TAD functionality. We anticipate that investigating defined sequence libraries using IDR‐Screen for specific functions can facilitate discovering novel and functional regions of the disordered proteome as well as understand the impact of natural and disease variants in disordered segments.

Function and Regulation of Phase-Separated Biological Condensates

Achieving functional specificity while minimizing cost to fitness is a key constraint during evolution. Formation of biological condensates by liquid-liquid phase separation (LLPS) appears to serve as an important regulatory mechanism to generate moderate specificity in molecular recognition while maintaining a reasonable cost for fitness in terms of design complexity. Formation of biological condensates serves as a unique mechanism of molecular recognition achieving some level of specificity without a huge cost to fitness. Rapid formation of biological condensates in vivo induced by specific cellular or environmental triggers has been shown to be an important mechanism for increasing cellular fitness. Here we discuss the functions and regulation of biological condensates, especially those formed by LLPS, involving interactions between proteins and nucleic acids. These condensates are spatially isolated within the cytosol or nucleus and can facilitate specific biochemical functions under conditions such as stress. The misregulation of biological condensates resulting in nondynamic aggregates has been implicated in a number of diseases. Understanding the functional importance of biological condensates and their regulation opens doors for development of therapies targeting dysfunctional biological condensates, as well as spatiotemporal engineering of functions in cells.