Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91.

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

AMBRA1 regulates cyclin D to guard S-phase entry and genomic integrity.

Mammalian development, adult tissue homeostasis and the avoidance of severe diseases including cancer require a properly orchestrated cell cycle, as well as error-free genome maintenance. The key cell-fate decision to replicate the genome is controlled by two major signalling pathways that act in parallel-the MYC pathway and the cyclin D-cyclin-dependent kinase (CDK)-retinoblastoma protein (RB) pathway1,2. Both MYC and the cyclin D-CDK-RB axis are commonly deregulated in cancer, and this is associated with increased genomic instability. The autophagic tumour-suppressor protein AMBRA1 has been linked to the control of cell proliferation, but the underlying molecular mechanisms remain poorly understood. Here we show that AMBRA1 is an upstream master regulator of the transition from G1 to S phase and thereby prevents replication stress. Using a combination of cell and molecular approaches and in vivo models, we reveal that AMBRA1 regulates the abundance of D-type cyclins by mediating their degradation. Furthermore, by controlling the transition from G1 to S phase, AMBRA1 helps to maintain genomic integrity during DNA replication, which counteracts developmental abnormalities and tumour growth. Finally, we identify the CHK1 kinase as a potential therapeutic target in AMBRA1-deficient tumours. These results advance our understanding of the control of replication-phase entry and genomic integrity, and identify the AMBRA1-cyclin D pathway as a crucial cell-cycle-regulatory mechanism that is deeply interconnected with genomic stability in embryonic development and tumorigenesis.

Minimum information guidelines for experiments structurally characterizing intrinsically disordered protein regions.

An unambiguous description of an experiment, and the subsequent biological observation, is vital for accurate data interpretation. Minimum information guidelines define the fundamental complement of data that can support an unambiguous conclusion based on experimental observations. We present the Minimum Information About Disorder Experiments (MIADE) guidelines to define the parameters required for the wider scientific community to understand the findings of an experiment studying the structural properties of intrinsically disordered regions (IDRs). MIADE guidelines provide recommendations for data producers to describe the results of their experiments at source, for curators to annotate experimental data to community resources and for database developers maintaining community resources to disseminate the data. The MIADE guidelines will improve the interpretability of experimental results for data consumers, facilitate direct data submission, simplify data curation, improve data exchange among repositories and standardize the dissemination of the key metadata on an IDR experiment by IDR data sources.

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks.

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

TFG binds LC3C to regulate ULK1 localization and autophagosome formation.

The early secretory pathway and autophagy are two essential and evolutionarily conserved endomembrane processes that are finely interlinked. Although growing evidence suggests that intracellular trafficking is important for autophagosome biogenesis, the molecular regulatory network involved is still not fully defined. In this study, we demonstrate a crucial effect of the COPII vesicle-related protein TFG (Trk-fused gene) on ULK1 puncta number and localization during autophagy induction. This, in turn, affects formation of the isolation membrane, as well as the correct dynamics of association between LC3B and early ATG proteins, leading to the proper formation of both omegasomes and autophagosomes. Consistently, fibroblasts derived from a hereditary spastic paraparesis (HSP) patient carrying mutated TFG (R106C) show defects in both autophagy and ULK1 puncta accumulation. In addition, we demonstrate that TFG activity in autophagy depends on its interaction with the ATG8 protein LC3C through a canonical LIR motif, thereby favouring LC3C-ULK1 binding. Altogether, our results uncover a link between TFG and autophagy and identify TFG as a molecular scaffold linking the early secretion pathway to autophagy.

A workflow to study mechanistic indicators for driver gene prediction with Moonlight.

Prediction of driver genes (tumor suppressors and oncogenes) is an essential step in understanding cancer development and discovering potential novel treatments. We recently proposed Moonlight as a bioinformatics framework to predict driver genes and analyze them in a system-biology-oriented manner based on -omics integration. Moonlight uses gene expression as a primary data source and combines it with patterns related to cancer hallmarks and regulatory networks to identify oncogenic mediators. Once the oncogenic mediators are identified, it is important to include extra levels of evidence, called mechanistic indicators, to identify driver genes and to link the observed gene expression changes to the underlying alteration that promotes them. Such a mechanistic indicator could be for example a mutation in the regulatory regions for the candidate gene. Here, we developed new functionalities and released Moonlight2 to provide the user with a mutation-based mechanistic indicator as a second layer of evidence. These functionalities analyze mutations in a cancer cohort to classify them into driver and passenger mutations. Those oncogenic mediators with at least one driver mutation are retained as the final set of driver genes. We applied Moonlight2 to the basal-like breast cancer subtype, lung adenocarcinoma and thyroid carcinoma using data from The Cancer Genome Atlas. For example, in basal-like breast cancer, we found four oncogenes (COPZ2, SF3B4, KRTCAP2 and POLR2J) and nine tumor suppressor genes (KIR2DL4, KIF26B, ARL15, ARHGAP25, EMCN, GMFG, TPK1, NR5A2 and TEK) containing a driver mutation in their promoter region, possibly explaining their deregulation. Moonlight2R is available at https://github.com/ELELAB/Moonlight2R.

MutateX: an automated pipeline for in silico saturation mutagenesis of protein structures and structural ensembles.

Mutations, which result in amino acid substitutions, influence the stability of proteins and their binding to biomolecules. A molecular understanding of the effects of protein mutations is both of biotechnological and medical relevance. Empirical free energy functions that quickly estimate the free energy change upon mutation (ΔΔG) can be exploited for systematic screenings of proteins and protein complexes. In silico saturation mutagenesis can guide the design of new experiments or rationalize the consequences of known mutations. Often software such as FoldX, while fast and reliable, lack the necessary automation features to apply them in a high-throughput manner. We introduce MutateX, a software to automate the prediction of ΔΔGs associated with the systematic mutation of each residue within a protein, or protein complex to all other possible residue types, using the FoldX energy function. MutateX also supports ΔΔG calculations over protein ensembles, upon post-translational modifications and in multimeric assemblies. At the heart of MutateX lies an automated pipeline engine that handles input preparation, parallelization and outputs publication-ready figures. We illustrate the MutateX protocol applied to different case studies. The results of the high-throughput scan provided by our tools can help in different applications, such as the analysis of disease-associated mutations, to complement experimental deep mutational scans, or assist the design of variants for industrial applications. MutateX is a collection of Python tools that relies on open-source libraries. It is available free of charge under the GNU General Public License from https://github.com/ELELAB/mutatex.

DisProt in 2022: improved quality and accessibility of protein intrinsic disorder annotation.

The Database of Intrinsically Disordered Proteins (DisProt, URL: https://disprot.org) is the major repository of manually curated annotations of intrinsically disordered proteins and regions from the literature. We report here recent updates of DisProt version 9, including a restyled web interface, refactored Intrinsically Disordered Proteins Ontology (IDPO), improvements in the curation process and significant content growth of around 30%. Higher quality and consistency of annotations is provided by a newly implemented reviewing process and training of curators. The increased curation capacity is fostered by the integration of DisProt with APICURON, a dedicated resource for the proper attribution and recognition of biocuration efforts. Better interoperability is provided through the adoption of the Minimum Information About Disorder (MIADE) standard, an active collaboration with the Gene Ontology (GO) and Evidence and Conclusion Ontology (ECO) consortia and the support of the ELIXIR infrastructure.

PDBminer to Find and Annotate Protein Structures for Computational Analysis.

Computational methods relying on protein structure strongly depend on the structure selected for investigation. Typical sources of protein structures include experimental structures available at the Protein Data Bank (PDB) and high-quality in silico model structures, such as those available at the AlphaFold Protein Structure Database. Either option has significant advantages and drawbacks, and exploring the wealth of available structures to identify the most suitable ones for specific applications can be a daunting task. We provide an open-source software package, PDBminer, with the purpose of making structure identification and selection easier, faster, and less error prone. PDBminer searches the AlphaFold Database and the PDB for available structures of interest and provides an up-to-date, quality-ranked table of structures applicable for further use. PDBminer provides an overview of the available protein structures to one or more input proteins, parallelizing the runs if multiple cores are specified. The output table reports the coverage of the protein structures aligned to the UniProt sequence, overcoming numbering differences in PDB structures and providing information regarding model quality, protein complexes, ligands, and nucleic acid chain binding. The PDBminer2coverage and PDBminer2network tools assist in visualizing the results. PDBminer can be applied to overcome the tedious task of choosing a PDB structure without losing the wealth of additional information available in the PDB. Here, we showcase the main functionalities of the package on the p53 tumor suppressor protein. The package is available at http://github.com/ELELAB/PDBminer.

PyInteraph2 and PyInKnife2 to Analyze Networks in Protein Structural Ensembles.

Due to the complex nature of noncovalent interactions and their long-range effects, analyzing protein conformations using network theory can be enlightening. Protein Structure Networks (PSNs) provide a convenient formalism to study protein structures in relation to essential properties such as key residues for structural stability, allosteric communication, and the effects of modifications of the protein. PSNs can be defined according to very different principles, and the available tools have limitations in input formats, supported models, and version control. Other outstanding problems are related to the definition of network cutoffs and the assessment of the stability of the network properties. The protein science community could benefit from a common framework to carry out these analyses and make them easier to reproduce, reuse, and evaluate. We here provide two open-source software packages, PyInteraph2 and PyInKnife2, to implement and analyze PSNs in a reproducible and documented manner. PyInteraph2 interfaces with multiple formats for protein ensembles and incorporates different network models with the possibility of integrating them into a macronetwork and performing various downstream analyses, including hubs, connected components, and several other centrality measures, and visualizes the networks or further analyzes them thanks to compatibility with Cytoscape.PyInKnife2 that supports the network models implemented in PyInteraph2. It employs a jackknife resampling approach to estimate the convergence of network properties and streamline the selection of distance cutoffs. We foresee that the modular structure of the code and the supported version control system will promote the transition to a community-driven effort, boost reproducibility, and establish common protocols in the PSN field. As developers, we will guarantee the introduction of new functionalities and maintenance, assistance, and training of new contributors.

Disease-associated mutations impacting BC-loop flexibility trigger long-range transthyretin tetramer destabilization and aggregation.

Hereditary transthyretin amyloidosis (ATTR) is an autosomal dominant disease characterized by the extracellular deposition of the transport protein transthyretin (TTR) as amyloid fibrils. Despite the progress achieved in recent years, understanding why different TTR residue substitutions lead to different clinical manifestations remains elusive. Here, we studied the molecular basis of disease-causing missense mutations affecting residues R34 and K35. R34G and K35T variants cause vitreous amyloidosis, whereas R34T and K35N mutations result in amyloid polyneuropathy and restrictive cardiomyopathy. All variants are more sensitive to pH-induced dissociation and amyloid formation than the wild-type (WT)-TTR counterpart, specifically in the variants deposited in the eyes amyloid formation occurs close to physiological pHs. Chemical denaturation experiments indicate that all the mutants are less stable than WT-TTR, with the vitreous amyloidosis variants, R34G and K35T, being highly destabilized. Sequence-induced stabilization of the dimer-dimer interface with T119M rendered tetramers containing R34G or K35T mutations resistant to pH-induced aggregation. Because R34 and K35 are among the residues more distant to the TTR interface, their impact in this region is therefore theorized to occur at long range. The crystal structures of double mutants, R34G/T119M and K35T/T119M, together with molecular dynamics simulations indicate that their strong destabilizing effect is initiated locally at the BC loop, increasing its flexibility in a mutation-dependent manner. Overall, the present findings help us to understand the sequence-dynamic-structural mechanistic details of TTR amyloid aggregation triggered by R34 and K35 variants and to link the degree of mutation-induced conformational flexibility to protein aggregation propensity.

DisProt: intrinsic protein disorder annotation in 2020.

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the 'dark' proteome.

Molecular dynamics-guided discovery of an ago-allosteric modulator for GPR40/FFAR1.

The long-chain fatty acid receptor FFAR1/GPR40 binds agonists in both an interhelical site between the extracellular segments of transmembrane helix (TM)-III and TM-IV and a lipid-exposed groove between the intracellular segments of these helices. Molecular dynamics simulations of FFAR1 with agonist removed demonstrated a major rearrangement of the polar and charged anchor point residues for the carboxylic acid moiety of the agonist in the interhelical site, which was associated with closure of a neighboring, solvent-exposed pocket between the extracellular poles of TM-I, TM-II, and TM-VII. A synthetic compound designed to bind in this pocket, and thereby prevent its closure, was identified through structure-based virtual screening and shown to function both as an agonist and as an allosteric modulator of receptor activation. This discovery of an allosteric agonist for a previously unexploited, dynamic pocket in FFAR1 demonstrates both the power of including molecular dynamics in the drug discovery process and that this specific, clinically proven, but difficult, antidiabetes target can be addressed by chemotypes different from existing ligands.

Targeting cancer stem cells in medulloblastoma by inhibiting AMBRA1 dual function in autophagy and STAT3 signalling.

Medulloblastoma (MB) is a childhood malignant brain tumour comprising four main subgroups characterized by different genetic alterations and rate of mortality. Among MB subgroups, patients with enhanced levels of the c-MYC oncogene (MBGroup3) have the poorest prognosis. Here we identify a previously unrecognized role of the pro-autophagy factor AMBRA1 in regulating MB. We demonstrate that AMBRA1 expression depends on c-MYC levels and correlates with Group 3 patient poor prognosis; also, knockdown of AMBRA1 reduces MB stem potential, growth and migration of MBGroup3 stem cells. At a molecular level, AMBRA1 mediates these effects by suppressing SOCS3, an inhibitor of STAT3 activation. Importantly, pharmacological inhibition of autophagy profoundly affects both stem and invasion potential of MBGroup3 stem cells, and a combined anti-autophagy and anti-STAT3 approach impacts the MBGroup3 outcome. Taken together, our data support the c-MYC/AMBRA1/STAT3 axis as a strong oncogenic signalling pathway with significance for both patient stratification strategies and targeted treatments of MBGroup3.

CAncer bioMarker Prediction Pipeline (CAMPP)-A standardized framework for the analysis of quantitative biological data.

With the improvement of -omics and next-generation sequencing (NGS) methodologies, along with the lowered cost of generating these types of data, the analysis of high-throughput biological data has become standard both for forming and testing biomedical hypotheses. Our knowledge of how to normalize datasets to remove latent undesirable variances has grown extensively, making for standardized data that are easily compared between studies. Here we present the CAncer bioMarker Prediction Pipeline (CAMPP), an open-source R-based wrapper (https://github.com/ELELAB/CAncer-bioMarker-Prediction-Pipeline -CAMPP) intended to aid bioinformatic software-users with data analyses. CAMPP is called from a terminal command line and is supported by a user-friendly manual. The pipeline may be run on a local computer and requires little or no knowledge of programming. To avoid issues relating to R-package updates, a renv .lock file is provided to ensure R-package stability. Data-management includes missing value imputation, data normalization, and distributional checks. CAMPP performs (I) k-means clustering, (II) differential expression/abundance analysis, (III) elastic-net regression, (IV) correlation and co-expression network analyses, (V) survival analysis, and (VI) protein-protein/miRNA-gene interaction networks. The pipeline returns tabular files and graphical representations of the results. We hope that CAMPP will assist in streamlining bioinformatic analysis of quantitative biological data, whilst ensuring an appropriate bio-statistical framework.

Secreted breast tumor interstitial fluid microRNAs and their target genes are associated with triple-negative breast cancer, tumor grade, and immune infiltration.

BACKGROUND: Studies on tumor-secreted microRNAs point to a functional role of these in cellular communication and reprogramming of the tumor microenvironment. Uptake of tumor-secreted microRNAs by neighboring cells may result in the silencing of mRNA targets and, in turn, modulation of the transcriptome. Studying miRNAs externalized from tumors could improve cancer patient diagnosis and disease monitoring and help to pinpoint which miRNA-gene interactions are central for tumor properties such as invasiveness and metastasis. METHODS: Using a bioinformatics approach, we analyzed the profiles of secreted tumor and normal interstitial fluid (IF) microRNAs, from women with breast cancer (BC). We carried out differential abundance analysis (DAA), to obtain miRNAs, which were enriched or depleted in IFs, from patients with different clinical traits. Subsequently, miRNA family enrichment analysis was performed to assess whether any families were over-represented in the specific sets. We identified dysregulated genes in tumor tissues from the same cohort of patients and constructed weighted gene co-expression networks, to extract sets of co-expressed genes and co-abundant miRNAs. Lastly, we integrated miRNAs and mRNAs to obtain interaction networks and supported our findings using prediction tools and cancer gene databases. RESULTS: Network analysis showed co-expressed genes and miRNA regulators, associated with tumor lymphocyte infiltration. All of the genes were involved in immune system processes, and many had previously been associated with cancer immunity. A subset of these, BTLA, CXCL13, IL7R, LAMP3, and LTB, was linked to the presence of tertiary lymphoid structures and high endothelial venules within tumors. Co-abundant tumor interstitial fluid miRNAs within this network, including miR-146a and miR-494, were annotated as negative regulators of immune-stimulatory responses. One co-expression network encompassed differences between BC subtypes. Genes differentially co-expressed between luminal B and triple-negative breast cancer (TNBC) were connected with sphingolipid metabolism and predicted to be co-regulated by miR-23a. Co-expressed genes and TIF miRNAs associated with tumor grade were BTRC, CHST1, miR-10a/b, miR-107, miR-301a, and miR-454. CONCLUSION: Integration of IF miRNAs and mRNAs unveiled networks associated with patient clinicopathological traits, and underlined molecular mechanisms, specific to BC sub-groups. Our results highlight the benefits of an integrative approach to biomarker discovery, placing secreted miRNAs within a biological context.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level.

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site.

New functionalities in the TCGAbiolinks package for the study and integration of cancer data from GDC and GTEx.

The advent of Next-Generation Sequencing (NGS) technologies has opened new perspectives in deciphering the genetic mechanisms underlying complex diseases. Nowadays, the amount of genomic data is massive and substantial efforts and new tools are required to unveil the information hidden in the data. The Genomic Data Commons (GDC) Data Portal is a platform that contains different genomic studies including the ones from The Cancer Genome Atlas (TCGA) and the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiatives, accounting for more than 40 tumor types originating from nearly 30000 patients. Such platforms, although very attractive, must make sure the stored data are easily accessible and adequately harmonized. Moreover, they have the primary focus on the data storage in a unique place, and they do not provide a comprehensive toolkit for analyses and interpretation of the data. To fulfill this urgent need, comprehensive but easily accessible computational methods for integrative analyses of genomic data that do not renounce a robust statistical and theoretical framework are required. In this context, the R/Bioconductor package TCGAbiolinks was developed, offering a variety of bioinformatics functionalities. Here we introduce new features and enhancements of TCGAbiolinks in terms of i) more accurate and flexible pipelines for differential expression analyses, ii) different methods for tumor purity estimation and filtering, iii) integration of normal samples from other platforms iv) support for other genomics datasets, exemplified here by the TARGET data. Evidence has shown that accounting for tumor purity is essential in the study of tumorigenesis, as these factors promote confounding behavior regarding differential expression analysis. With this in mind, we implemented these filtering procedures in TCGAbiolinks. Moreover, a limitation of some of the TCGA datasets is the unavailability or paucity of corresponding normal samples. We thus integrated into TCGAbiolinks the possibility to use normal samples from the Genotype-Tissue Expression (GTEx) project, which is another large-scale repository cataloging gene expression from healthy individuals. The new functionalities are available in the TCGAbiolinks version 2.8 and higher released in Bioconductor version 3.7.

Predicting the impact of Lynch syndrome-causing missense mutations from structural calculations.

Accurate methods to assess the pathogenicity of mutations are needed to fully leverage the possibilities of genome sequencing in diagnosis. Current data-driven and bioinformatics approaches are, however, limited by the large number of new variations found in each newly sequenced genome, and often do not provide direct mechanistic insight. Here we demonstrate, for the first time, that saturation mutagenesis, biophysical modeling and co-variation analysis, performed in silico, can predict the abundance, metabolic stability, and function of proteins inside living cells. As a model system, we selected the human mismatch repair protein, MSH2, where missense variants are known to cause the hereditary cancer predisposition disease, known as Lynch syndrome. We show that the majority of disease-causing MSH2 mutations give rise to folding defects and proteasome-dependent degradation rather than inherent loss of function, and accordingly our in silico modeling data accurately identifies disease-causing mutations and outperforms the traditionally used genetic disease predictors. Thus, in conclusion, in silico biophysical modeling should be considered for making genotype-phenotype predictions and for diagnosis of Lynch syndrome, and perhaps other hereditary diseases.

Toward mechanistic models for genotype-phenotype correlations in phenylketonuria using protein stability calculations.

Phenylketonuria (PKU) is a genetic disorder caused by variants in the gene encoding phenylalanine hydroxylase (PAH), resulting in accumulation of phenylalanine to neurotoxic levels. Here, we analyzed the cellular stability, localization, and interaction with wild-type PAH of 20 selected PKU-linked PAH protein missense variants. Several were present at reduced levels in human cells, and the levels increased in the presence of a proteasome inhibitor, indicating that proteins are proteasome targets. We found that all the tested PAH variants retained their ability to associate with wild-type PAH, and none formed aggregates, suggesting that they are only mildly destabilized in structure. In all cases, PAH variants were stabilized by the cofactor tetrahydrobiopterin (BH4 ), a molecule known to alleviate symptoms in certain PKU patients. Biophysical calculations on all possible single-site missense variants using the full-length structure of PAH revealed a strong correlation between the predicted protein stability and the observed stability in cells. This observation rationalizes previously observed correlations between predicted loss of protein destabilization and disease severity, a correlation that we also observed using new calculations. We thus propose that many disease-linked PAH variants are structurally destabilized, which in turn leads to proteasomal degradation and insufficient amounts of cellular PAH protein.