A saturation-mutagenesis analysis of the interplay between stability and activation in Ras

Comprehensive structure-function analysis reveals gain- and loss-of-function mechanisms impacting oncogenic KRAS activity

To dissect variant-function relationships in the KRAS oncoprotein, we performed deep mutational scanning (DMS) screens for both wild-type and KRAS G12D mutant alleles. We defined the spectrum of oncogenic potential for nearly all possible KRAS variants, identifying several novel transforming alleles and elucidating a model to describe the frequency of KRAS mutations in human cancer as a function of transforming potential, mutational probability, and tissue-specific mutational signatures. Biochemical and structural analyses of variants identified in a KRAS G12D second-site suppressor DMS screen revealed that attenuation of oncogenic KRAS can be mediated by protein instability and conformational rigidity, resulting in reduced binding affinity to effector proteins, such as RAF and PI3-kinases, or reduced SOS-mediated nucleotide exchange activity. These studies define the landscape of single amino acid alterations that modulate the function of KRAS, providing a resource for the clinical interpretation of KRAS variants and elucidating mechanisms of oncogenic KRAS inactivation for therapeutic exploitation.

Functional distinction in oncogenic Ras variant activity in Caenorhabditis elegans

Ras genes are important oncogenes that are frequently mutated in cancer. Human oncogenic variants exhibit functional distinctions in terms of their representation in different cancer types, impact on cellular targets and sensitivity to pharmacological treatments. However, how these distinct variants influence and respond to the cellular networks in which they are embedded is poorly understood. To identify novel participants in the complex interplay between Ras genotype and cell interaction networks in vivo, we have developed and tested an experimental framework using a simple vulva-development assay in the nematode C. elegans. Using this system, we evaluated a set of Ras oncogenic substitution changes at G12, G13 and Q61. We found that these variants fall into distinct groups based on phenotypic differences, sensitivity to gene dosage and inhibition of the downstream kinase MEK and their response to genetic modulators that influence Ras activity in a non-autonomous manner. Together, our results demonstrated that oncogenic C. elegans Ras variants exhibit clear distinctions in how they interface with the vulva-development network and showed that extracellular modulators yield variant-restricted effects in vivo.

Structure-based network analysis predicts pathogenic variants in human proteins associated with inherited retinal disease

Advances in gene sequencing technologies have accelerated the identification of genetic variants, but better tools are needed to understand which are causal of disease. This would be particularly useful in fields where gene therapy is a potential therapeutic modality for a disease-causing variant such as inherited retinal disease (IRD). Here, we apply structure-based network analysis (SBNA), which has been successfully utilized to identify variant-constrained amino acid residues in viral proteins, to identify residues that may cause IRD if subject to missense mutation. SBNA is based entirely on structural first principles and is not fit to specific outcome data, which makes it distinct from other contemporary missense prediction tools. In 4 well-studied human disease-associated proteins (BRCA1, HRAS, PTEN, and ERK2) with high-quality structural data, we find that SBNA scores correlate strongly with deep mutagenesis data. When applied to 47 IRD genes with available high-quality crystal structure data, SBNA scores reliably identified disease-causing variants according to phenotype definitions from the ClinVar database. Finally, we applied this approach to 63 patients at Massachusetts Eye and Ear (MEE) with IRD but for whom no genetic cause had been identified. Untrained models built using SBNA scores and BLOSUM62 scores for IRD-associated genes successfully predicted the pathogenicity of novel variants (AUC = 0.851), allowing us to identify likely causative disease variants in 40 IRD patients. Model performance was further augmented by incorporating orthogonal data from EVE scores (AUC = 0.927), which are based on evolutionary multiple sequence alignments. In conclusion, SBNA can used to successfully identify variants as causal of disease in human proteins and may help predict variants causative of IRD in an unbiased fashion.

Mutational scanning pinpoints distinct binding sites of key ATGL regulators in lipolysis

ATGL is a key enzyme in intracellular lipolysis and plays an important role in metabolic and cardiovascular diseases. ATGL is tightly regulated by a known set of protein-protein interaction partners with activating or inhibiting functions in the control of lipolysis. Here, we use deep mutational protein interaction perturbation scanning and generate comprehensive profiles of single amino acid variants that affect the interactions of ATGL with its regulatory partners: CGI-58, G0S2, PLIN1, PLIN5 and CIDEC. Twenty-three ATGL amino acid variants yield a specific interaction perturbation pattern when validated in co-immunoprecipitation experiments in mammalian cells. We identify and characterize eleven highly selective ATGL switch mutations which affect the interaction of one of the five partners without affecting the others. Switch mutations thus provide distinct interaction determinants for ATGL's key regulatory proteins at an amino acid resolution. When we test triglyceride hydrolase activity in vitro and lipolysis in cells, the activity patterns of the ATGL switch variants trace to their protein interaction profile. In the context of structural data, the integration of variant binding and activity profiles provides insights into the regulation of lipolysis and the impact of mutations in human disease.

The energetic and allosteric landscape for KRAS inhibition

Thousands of proteins have been validated genetically as therapeutic targets for human diseases1. However, very few have been successfully targeted, and many are considered 'undruggable'. This is particularly true for proteins that function via protein-protein interactions-direct inhibition of binding interfaces is difficult and requires the identification of allosteric sites. However, most proteins have no known allosteric sites, and a comprehensive allosteric map does not exist for any protein. Here we address this shortcoming by charting multiple global atlases of inhibitory allosteric communication in KRAS. We quantified the effects of more than 26,000 mutations on the folding of KRAS and its binding to six interaction partners. Genetic interactions in double mutants enabled us to perform biophysical measurements at scale, inferring more than 22,000 causal free energy changes. These energy landscapes quantify how mutations tune the binding specificity of a signalling protein and map the inhibitory allosteric sites for an important therapeutic target. Allosteric propagation is particularly effective across the central β-sheet of KRAS, and multiple surface pockets are genetically validated as allosterically active, including a distal pocket in the C-terminal lobe of the protein. Allosteric mutations typically inhibit binding to all tested effectors, but they can also change the binding specificity, revealing the regulatory, evolutionary and therapeutic potential to tune pathway activation. Using the approach described here, it should be possible to rapidly and comprehensively identify allosteric target sites in many proteins.

RAS signaling and immune cells: a sinister crosstalk in the tumor microenvironment

The rat sarcoma virus (RAS) gene is the most commonly mutated oncogene in cancer, with about 19% of cancer patients carrying RAS mutations. Studies on the interaction between RAS mutation and tumor immune microenvironment (TIM) have been flourishing in recent years. More and more evidence has proved that RAS signals regulate immune cells' recruitment, activation, and differentiation while assisting tumor cells to evade immune surveillance. This review concluded the direct and indirect treatment strategies for RAS mutations. In addition, we updated the underlying mechanisms by which RAS signaling modulated immune infiltration and immune escape. Finally, we discussed advances in RAS-targeted immunotherapies, including cancer vaccines and adoptive cell therapies, with a particular focus on combination strategies with personalized therapy and great potential to achieve lasting clinical benefits.

Structure-based network analysis predicts mutations associated with inherited retinal disease

With continued advances in gene sequencing technologies comes the need to develop better tools to understand which mutations cause disease. Here we validate structure-based network analysis (SBNA)1,2 in well-studied human proteins and report results of using SBNA to identify critical amino acids that may cause retinal disease if subject to missense mutation. We computed SBNA scores for genes with high-quality structural data, starting with validating the method using 4 well-studied human disease-associated proteins. We then analyzed 47 inherited retinal disease (IRD) genes. We compared SBNA scores to phenotype data from the ClinVar database and found a significant difference between benign and pathogenic mutations with respect to network score. Finally, we applied this approach to 65 patients at Massachusetts Eye and Ear (MEE) who were diagnosed with IRD but for whom no genetic cause was found. Multivariable logistic regression models built using SBNA scores for IRD-associated genes successfully predicted pathogenicity of novel mutations, allowing us to identify likely causative disease variants in 37 patients with IRD from our clinic. In conclusion, SBNA can be meaningfully applied to human proteins and may help predict mutations causative of IRD.

Emerging maps of allosteric regulation in cellular networks

Allosteric regulation is classically defined as action at a distance, where a perturbation outside of a protein active site affects function. While this definition has motivated many studies of allosteric mechanisms at the level of protein structure, translating these insights to the allosteric regulation of entire cellular processes - and their crosstalk - has received less attention, despite the broad importance of allostery for cellular regulation foreseen by Jacob and Monod. Here, we revisit an evolutionary model for the widespread emergence of allosteric regulation in colocalized proteins, describe supporting evidence, and discuss emerging advances in mapping allostery in cellular networks that link precise and often allosteric perturbations at the molecular level to functional changes at the pathway and systems levels.

A complete allosteric map of a GTPase switch in its native cellular network

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

Saturation mutagenesis of a predicted ancestral Syk-family kinase

Many tyrosine kinases cannot be expressed readily in Escherichia coli, limiting facile production of these proteins for biochemical experiments. We used ancestral sequence reconstruction to generate a spleen tyrosine kinase (Syk) variant that can be expressed in bacteria and purified in soluble form, unlike the human members of this family (Syk and zeta-chain-associated protein kinase of 70 kDa [ZAP-70]). The catalytic activity, substrate specificity, and regulation by phosphorylation of this Syk variant are similar to the corresponding properties of human Syk and ZAP-70. Taking advantage of the ability to express this novel Syk-family kinase in bacteria, we developed a two-hybrid assay that couples the growth of E. coli in the presence of an antibiotic to successful phosphorylation of a bait peptide by the kinase. Using this assay, we screened a site-saturation mutagenesis library of the kinase domain of this reconstructed Syk-family kinase. Sites of loss-of-function mutations identified in the screen correlate well with residues established previously as critical to function and/or structure in protein kinases. We also identified activating mutations in the regulatory hydrophobic spine and activation loop, which are within key motifs involved in kinase regulation. Strikingly, one mutation in an ancestral Syk-family variant increases the soluble expression of the protein by 75-fold. Thus, through ancestral sequence reconstruction followed by deep mutational scanning, we have generated Syk-family kinase variants that can be expressed in bacteria with very high yield.