A family of evolution-entropy hybrid methods for ranking protein residues by importance

Unraveling the potential effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on the Protein structure and function of the human SLC30A8 gene on type 2 diabetes and colorectal cancer: An In silico approach

Background and aims

The single nucleotide polymorphisms (SNPs) in SLC30A8 gene have been recognized as contributing to type 2 diabetes (T2D) susceptibility and colorectal cancer. This study aims to predict the structural stability, and functional impacts on variations in non-synonymous SNPs (nsSNPs) in the human SLC30A8 gene using various computational techniques.

Materials and methods

Several in silico tools, including SIFT, Predict-SNP, SNPs&GO, MAPP, SNAP2, PhD-SNP, PANTHER, PolyPhen-1,PolyPhen-2, I-Mutant 2.0, and MUpro, have been used in our study.

Results

After data analysis, out of 336 missenses, the eight nsSNPs, namely R138Q, I141N, W136G, I349N, L303R, E140A, W306C, and L308Q, were discovered by ConSurf to be in highly conserved regions, which could affect the stability of their proteins. Project HOPE determines any significant molecular effects on the structure and function of eight mutated proteins and the three-dimensional (3D) structures of these proteins. The two pharmacologically significant compounds, Luzonoid B and Roseoside demonstrate strong binding affinity to the mutant proteins, and they are more efficient in inhibiting them than the typical SLC30A8 protein using Autodock Vina and Chimera. Increased binding affinity to mutant SLC30A8 proteins has been determined not to influence drug resistance. Ultimately, the Kaplan-Meier plotter study revealed that alterations in SLC30A8 gene expression notably affect the survival rates of patients with various cancer types.

Conclusion

Finally, the study found eight highly deleterious missense nsSNPs in the SLC30A8 gene that can be helpful for further proteomic and genomic studies for T2D and colorectal cancer diagnosis. These findings also pave the way for personalized treatments using biomarkers and more effective healthcare strategies.

Conformational Dynamics, Energetics, and the Divergent Evolution of Allosteric Regulation: The Case of the Yeast MAPK Family

Allosteric mechanisms provide finely-tuned control over signalling proteins. Proteins of the same family may share high sequence identity and structural similarity but show distinct traits of allosteric control and evolutionary divergent regulation. Revealing the determinants of such properties may be important to understand the molecular bases of different regulatory pathways. Herein, we investigate whether and how evolutionarily-divergent traits of allosteric regulation in homologous proteins can be decoded in terms of internal dynamics and interaction networks that support functionally oriented conformations. In this framework, we start from the comparative analysis of the dynamics and energetics of the yeast MAP Kinases (MAPKs) Fus3 and Kss1 in their native basins. Importantly, distinctive dynamic and energetic stabilization features emerge, which can be related to the two proteins' differential ability to be phosphorylated and engage with the allosteric activator Ste5. We then expanded our study to other evolutionarily-related MAPKs. We show that the dynamical and energetical traits defining the distinct regulatory profiles of Fus3 and Kss1 can be traced along their evolutionary tree. Overall, our approach is able to reconnect (latent) allostery with the principal elements of protein structural stabilization and dynamics, showing how allosteric regulation was encrypted in MAPKs structure well before Ste5 appearance.

Variants in ZFX are associated with an X-linked neurodevelopmental disorder with recurrent facial gestalt

Pathogenic variants in multiple genes on the X chromosome have been implicated in syndromic and non-syndromic intellectual disability disorders. ZFX on Xp22.11 encodes a transcription factor that has been linked to diverse processes including oncogenesis and development, but germline variants have not been characterized in association with disease. Here, we present clinical and molecular characterization of 18 individuals with germline ZFX variants. Exome or genome sequencing revealed 11 variants in 18 subjects (14 males and 4 females) from 16 unrelated families. Four missense variants were identified in 11 subjects, with seven truncation variants in the remaining individuals. Clinical findings included developmental delay/intellectual disability, behavioral abnormalities, hypotonia, and congenital anomalies. Overlapping and recurrent facial features were identified in all subjects, including thickening and medial broadening of eyebrows, variations in the shape of the face, external eye abnormalities, smooth and/or long philtrum, and ear abnormalities. Hyperparathyroidism was found in four families with missense variants, and enrichment of different tumor types was observed. In molecular studies, DNA-binding domain variants elicited differential expression of a small set of target genes relative to wild-type ZFX in cultured cells, suggesting a gain or loss of transcriptional activity. Additionally, a zebrafish model of ZFX loss displayed an altered behavioral phenotype, providing additional evidence for the functional significance of ZFX. Our clinical and experimental data support that variants in ZFX are associated with an X-linked intellectual disability syndrome characterized by a recurrent facial gestalt, neurocognitive and behavioral abnormalities, and an increased risk for congenital anomalies and hyperparathyroidism.

An Update on Protein Kinases as Therapeutic Targets-Part II: Peptides as Allosteric Protein Kinase C Modulators Targeting Protein-Protein Interactions

Human protein kinases are highly-sought-after drug targets, historically harnessed for treating cancer, cardiovascular disease, and an increasing number of autoimmune and inflammatory conditions. Most current treatments involve small molecule protein kinase inhibitors that interact orthosterically with the protein kinase ATP-binding pocket. As a result, these compounds are often poorly selective and highly toxic. Part I of this series reviews the role of PKC isoforms in various human diseases, featuring cancer and cardiovascular disease, as well as translational examples of PKC modulation applied to human health and disease. In the present Part II, we discuss alternative allosteric binding mechanisms for targeting PKC, as well as novel drug platforms, such as modified peptides. A major goal is to design protein kinase modulators with enhanced selectivity and improved pharmacological properties. To this end, we use molecular docking analysis to predict the mechanisms of action for inhibitor-kinase interactions that can facilitate the development of next-generation PKC modulators.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

H2S-Generating Cytosolic L-Cysteine Desulfhydrase and Mitochondrial D-Cysteine Desulfhydrase from Sweet Pepper (Capsicum annuum L.) Are Regulated During Fruit Ripening and by Nitric Oxide

Aims: Pepper fruit is a horticultural product worldwide consumed that has great nutritional and economic relevance. Besides the phenotypical changes that undergo pepper fruit during ripening, there are many associated modifications at transcriptomic, proteomic, biochemical, and metabolic levels. Nitric oxide (NO) and hydrogen sulfide (H2S) are recognized signal molecules that can exert regulatory functions in diverse plant processes. This study aims at analyzing the interrelationship between NO and H2S during fruit ripening. Results: Our data indicate that the H2S-generating cytosolic L-cysteine desulfhydrase (LCD) and the mitochondrial D-cysteine desulfhydrase (DCD) activities are downregulated during ripening but this effect was reverted after NO treatment of fruits. Innovation and Conclusion: Using as a model the non-climacteric pepper fruits at different ripening stages and under an NO-enriched atmosphere, the activity of the H2S-generating LCD and DCD was analyzed. LCD and DCD activities were downregulated during ripening, but this effect was reverted after NO treatment of fruits. The analysis of LCD activity by non-denaturing polyacrylamide gel electrophoresis (PAGE) allowed identifying three isozymes designated CaLCD I to CaLCD III, which were differentially modulated by NO and strictly dependent on pyridoxal 5'-phosphate (PLP). In vitro analyses of green fruit samples in the presence of different compounds including NO donors, peroxynitrite (ONOO-), and reducing agents such as reduced glutathione (GSH) and L-cysteine (L-Cys) triggered an almost 100% inhibition of CaLCD II and CaLCD III. This redox adaptation process of both enzymes could be cataloged as a hormesis phenomenon. The protein tyrosine (Tyr) nitration (an NO-promoted post-translational modification) of the recombinant LCD was corroborated by immunoblot and by mass spectrometry (MS) analyses. Among the 11 Tyr residues present in this enzyme, MS of the recombinant LCD enabled us to identify that Tyr82 and Tyr254 were nitrated by ONOO-, this occurring near the active center on the enzyme, where His237 and Lys260 together with the cofactor PLP are involved. These data support the relationship between NO and H2S during pepper fruit ripening, since LCD and DCD are regulated by NO during this physiological event, and this could also be extrapolated to other plant species.

CovET: A covariation-evolutionary trace method that identifies protein structure-function modules

Measuring the relative effect that any two sequence positions have on each other may improve protein design or help better interpret coding variants. Current approaches use statistics and machine learning but rarely consider phylogenetic divergences which, as shown by Evolutionary Trace studies, provide insight into the functional impact of sequence perturbations. Here, we reframe covariation analyses in the Evolutionary Trace framework to measure the relative tolerance to perturbation of each residue pair during evolution. This approach (CovET) systematically accounts for phylogenetic divergences: at each divergence event, we penalize covariation patterns that belie evolutionary coupling. We find that while CovET approximates the performance of existing methods to predict individual structural contacts, it performs significantly better at finding structural clusters of coupled residues and ligand binding sites. For example, CovET found more functionally critical residues when we examined the RNA recognition motif and WW domains. It correlates better with large-scale epistasis screen data. In the dopamine D2 receptor, top CovET residue pairs recovered accurately the allosteric activation pathway characterized for Class A G protein-coupled receptors. These data suggest that CovET ranks highest the sequence position pairs that play critical functional roles through epistatic and allosteric interactions in evolutionarily relevant structure-function motifs. CovET complements current methods and may shed light on fundamental molecular mechanisms of protein structure and function.

The ribosome-inactivating proteins MAP30 and Momordin inhibit SARS-CoV-2

The continuing emergence of SARS-CoV-2 variants has highlighted the need to identify additional points for viral inhibition. Ribosome inactivating proteins (RIPs), such as MAP30 and Momordin which are derived from bitter melon (Momordica charantia), have been found to inhibit a broad range of viruses. MAP30 has been shown to potently inhibit HIV-1 with minimal cytotoxicity. Here we show that MAP30 and Momordin potently inhibit SARS-CoV-2 replication in A549 human lung cells (IC50 ~ 0.2 μM) with little concomitant cytotoxicity (CC50 ~ 2 μM). Both viral inhibition and cytotoxicity remain unaltered by appending a C-terminal Tat cell-penetration peptide to either protein. Mutation of tyrosine 70, a key residue in the active site of MAP30, to alanine completely abrogates both viral inhibition and cytotoxicity, indicating the involvement of its RNA N-glycosylase activity. Mutation of lysine 171 and lysine 215, residues corresponding to those in Ricin which when mutated prevented ribosome binding and inactivation, to alanine in MAP30 decreased cytotoxicity (CC50 ~ 10 μM) but also the viral inhibition (IC50 ~ 1 μM). Unlike with HIV-1, neither Dexamethasone nor Indomethacin exhibited synergy with MAP30 in the inhibition of SARS-CoV-2. From a structural comparison of the two proteins, one can explain their similar activities despite differences in both their active-sites and ribosome-binding regions. We also note points on the viral genome for potential inhibition by these proteins.

Evolutionary conservation of sequence motifs at sites of protein modification

Gene duplications are common in biology and are likely to be an important source of functional diversification and specialization. The yeast Saccharomyces cerevisiae underwent a whole-genome duplication event early in evolution, and a substantial number of duplicated genes have been retained. We identified more than 3500 instances where only one of two paralogous proteins undergoes posttranslational modification despite having retained the same amino acid residue in both. We also developed a web-based search algorithm (CoSMoS.c.) that scores conservation of amino acid sequences based on 1011 wild and domesticated yeast isolates and used it to compare differentially modified pairs of paralogous proteins. We found that the most common modifications-phosphorylation, ubiquitylation, and acylation but not N-glycosylation-occur in regions of high sequence conservation. Such conservation is evident even for ubiquitylation and succinylation, where there is no established 'consensus site' for modification. Differences in phosphorylation were not associated with predicted secondary structure or solvent accessibility but did mirror known differences in kinase-substrate interactions. Thus, differences in posttranslational modification likely result from differences in adjoining amino acids and their interactions with modifying enzymes. By integrating data from large-scale proteomics and genomics analysis, in a system with such substantial genetic diversity, we obtained a more comprehensive understanding of the functional basis for genetic redundancies that have persisted for 100 million years.

Substrate-binding glycine residues are major determinants for hydrolase and ligase activity of plant legumains

Peptide asparaginyl ligases (PALs) are useful tools for precision modifications of proteins and live-cell surfaces by ligating peptides after Asn/Asp (Asx). They share high sequence and structural similarity to plant legumains that are generally known as asparaginyl endopeptidases (AEPs), thus making it challenging to identify PALs from AEPs. In this study, we investigate 875 plant species from algae to seed plants with available sequence data in public databases to identify new PALs. We conducted evolutionary trace analysis on 1500 plant legumains, including eight known PALs, to identify key residues that could differentiate ligases and proteases, followed by recombinant expression and functional validation of 16 novel legumains. Previously, we showed that the substrate-binding sequences flanking the catalytic site can strongly influence the enzymatic direction of a legumain and which we named as ligase-activity determinants (LADs). Here, we show that two conserved substrate-binding Gly residues of LADs are critical, but negative determinants for ligase activity. Our results suggest that specific glycine residues are molecular determinants to identify PALs and AEPs as two different legumain subfamilies, accounting for c. 1% and 88%, respectively.

Soybean (Glycine max L.) Lipoxygenase 1 (LOX 1) Is Modulated by Nitric Oxide and Hydrogen Sulfide: An In Vitro Approach

Hydrogen sulfide (H₂S) and nitric oxide (NO) are two relevant signal molecules that can affect protein function throughout post-translational modifications (PTMs) such as persulfidation, S-nitrosation, metal-nitrosylation, and nitration. Lipoxygenases (LOXs) are a group of non-heme iron enzymes involved in a wide range of plant physiological functions including seed germination, plant growth and development, and fruit ripening and senescence. Likewise, LOXs are also involved in the mechanisms of response to diverse environmental stresses. Using purified soybean (Glycine max L.) lipoxygenase type 1 (LOX 1) and nitrosocysteine (CysNO) and sodium hydrosulfide (NaHS) as NO and H₂S donors, respectively, the present study reveals that both compounds negatively affect LOX activity, suggesting that S-nitrosation and persulfidation are involved. Mass spectrometric analysis of nitrated soybean LOX 1 using a peroxynitrite (ONOO^-) donor enabled us to identify that, among the thirty-five tyrosine residues present in this enzyme, only Y214 was exclusively nitrated by ONOO^-. The nitration of Y214 seems to affect its interaction with W500, a residue involved in the substrate binding site. The analysis of the structure 3PZW demonstrates the existence of several tunnels that directly communicate the surface of the protein with different internal cysteines, thus making feasible their potential persulfidation, especially C429 and C127. On the other hand, the CysNO molecule, which is hydrophilic and bulkier than H₂S, can somehow be accommodated throughout the tunnel until it reaches C127, thus facilitating its nitrosation. Overall, a large number of potential persulfidation targets and the ease by which H₂S can reach them through the diffuse tunneling network could be behind their efficient inhibition.

Coevolutionary signals in metabotropic glutamate receptors capture residue contacts and long-range functional interactions

Upon ligand binding to a G protein-coupled receptor, extracellular signals are transmitted into a cell through sets of residue interactions that translate ligand binding into structural rearrangements. These interactions needed for functions impose evolutionary constraints so that, on occasion, mutations in one position may be compensated by other mutations at functionally coupled positions. To quantify the impact of amino acid substitutions in the context of major evolutionary divergence in the G protein-coupled receptor subfamily of metabotropic glutamate receptors (mGluRs), we combined two phylogenetic-based algorithms, Evolutionary Trace and covariation Evolutionary Trace, to infer potential structure-function couplings and roles in mGluRs. We found a subset of evolutionarily important residues at known functional sites and evidence of coupling among distinct structural clusters in mGluR. In addition, experimental mutagenesis and functional assays confirmed that some highly covariant residues are coupled, revealing their synergy. Collectively, these findings inform a critical step toward understanding the molecular and structural basis of amino acid variation patterns within mGluRs and provide insight for drug development, protein engineering, and analysis of naturally occurring variants.

The Solvation of the E. coli CheY Phosphorylation Site Mapped by XFMS

The Escherichia coli CheY protein belongs to a large bacterial response regulator superfamily. X-ray hydroxy radical foot-printing with mass spectroscopy (XFMS) has shown that allosteric activation of CheY by its motor target triggers a concerted internalization of aromatic sidechains. We reanalyzed the XFMS data to compare polar versus non-polar CheY residue positions. The polar residues around and including the 57D phosphorylated site had an elevated hydroxy radical reactivity. Bioinformatic measures revealed that a water-mediated hydrogen bond network connected this ring of residues with the central 57D. These residues solvated 57D to energetically stabilize the apo-CheY fold. The abundance of these reactive residues was reduced upon activation. This result was supported by the bioinformatics and consistent with the previously reported activation-induced increase in core hydrophobicity. It further illustrated XFMS detection of structural waters. Direct contacts between the ring residues and the phosphorylation site would stabilize the aspartyl phosphate. In addition, we report that the ring residue, 18R, is a constant central node in the 57D solvation network and that 18R non-polar substitutions determine CheY diversity as assessed by its evolutionary trace in bacteria with well-studied chemotaxis. These results showcase the importance of structured water dynamics for phosphorylation-mediated signal transduction.

Genome interpretation using in silico predictors of variant impact

Estimating the effects of variants found in disease driver genes opens the door to personalized therapeutic opportunities. Clinical associations and laboratory experiments can only characterize a tiny fraction of all the available variants, leaving the majority as variants of unknown significance (VUS). In silico methods bridge this gap by providing instant estimates on a large scale, most often based on the numerous genetic differences between species. Despite concerns that these methods may lack reliability in individual subjects, their numerous practical applications over cohorts suggest they are already helpful and have a role to play in genome interpretation when used at the proper scale and context. In this review, we aim to gain insights into the training and validation of these variant effect predicting methods and illustrate representative types of experimental and clinical applications. Objective performance assessments using various datasets that are not yet published indicate the strengths and limitations of each method. These show that cautious use of in silico variant impact predictors is essential for addressing genome interpretation challenges.

Impact of TP53 mutations in Triple Negative Breast Cancer

Identifying triple negative breast cancer (TNBC) patients expected to have poor outcomes provides an opportunity to enhance clinical management. We applied an Evolutionary Action Score to functionally characterize TP53 mutations (EAp53) in 96 TNBC patients and observed that EAp53 stratification may identify TP53 mutations associated with worse outcomes. These findings merit further exploration in larger TNBC cohorts and in patients treated with neoadjuvant chemotherapy regimens.

Evolutionary action of mutations reveals antimicrobial resistance genes in Escherichia coli

Since antibiotic development lags, we search for potential drug targets through directed evolution experiments. A challenge is that many resistance genes hide in a noisy mutational background as mutator clones emerge in the adaptive population. Here, to overcome this noise, we quantify the impact of mutations through evolutionary action (EA). After sequencing ciprofloxacin or colistin resistance strains grown under different mutational regimes, we find that an elevated sum of the evolutionary action of mutations in a gene identifies known resistance drivers. This EA integration approach also suggests new antibiotic resistance genes which are then shown to provide a fitness advantage in competition experiments. Moreover, EA integration analysis of clinical and environmental isolates of antibiotic resistant of E. coli identifies gene drivers of resistance where a standard approach fails. Together these results inform the genetic basis of de novo colistin resistance and support the robust discovery of phenotype-driving genes via the evolutionary action of genetic perturbations in fitness landscapes.

Conformational spread drives the evolution of the calcium-calmodulin protein kinase II

The calcium calmodulin (Ca²⁺/CaM) dependent protein kinase II (CaMKII) decodes Ca²⁺ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca²⁺ frequency response range for complex learning functions.

Computational prediction of protein functional sites-Applications in biotechnology and biomedicine

There are many computational approaches for predicting protein functional sites based on different sequence and structural features. These methods are essential to cope with the sequence deluge that is filling databases with uncharacterized protein sequences. They complement the more expensive and time-consuming experimental approaches by pointing them to possible candidate positions. In many cases they are jointly used to characterize the functional sites in proteins of biotechnological and biomedical interest and eventually modify them for different purposes. There is a clear trend towards approaches based on machine learning and those using structural information, due to the recent developments in these areas. Nevertheless, "classic" methods based on sequence and evolutionary features are still playing an important role as these features are strongly related to functionality. In this review, the main approaches for predicting general functional sites in a protein are discussed, with a focus on sequence-based approaches.

Decoding Cancer Variants of Unknown Significance for Helicase-Nuclease-RPA Complexes Orchestrating DNA Repair During Transcription and Replication

All tumors have DNA mutations, and a predictive understanding of those mutations could inform clinical treatments. However, 40% of the mutations are variants of unknown significance (VUS), with the challenge being to objectively predict whether a VUS is pathogenic and supports the tumor or whether it is benign. To objectively decode VUS, we mapped cancer sequence data and evolutionary trace (ET) scores onto crystallography and cryo-electron microscopy structures with variant impacts quantitated by evolutionary action (EA) measures. As tumors depend on helicases and nucleases to deal with transcription/replication stress, we targeted helicase–nuclease–RPA complexes: (1) XPB-XPD (within TFIIH), XPF-ERCC1, XPG, and RPA for transcription and nucleotide excision repair pathways and (2) BLM, EXO5, and RPA plus DNA2 for stalled replication fork restart. As validation, EA scoring predicts severe effects for most disease mutations, but disease mutants with low ET scores not only are likely destabilizing but also disrupt sophisticated allosteric mechanisms. For sites of disease mutations and VUS predicted to be severe, we found strong co-localization to ordered regions. Rare discrepancies highlighted the different survival requirements between disease and tumor mutations, as well as the value of examining proteins within complexes. In a genome-wide analysis of 33 cancer types, we found correlation between the number of mutations in each tumor and which pathways or functional processes in which the mutations occur, revealing different mutagenic routes to tumorigenesis. We also found upregulation of ancient genes including BLM, which supports a non-random and concerted cancer process: reversion to a unicellular, proliferation-uncontrolled, status by breaking multicellular constraints on cell division. Together, these genes and global analyses challenge the binary “driver” and “passenger” mutation paradigm, support a gradient impact as revealed by EA scoring from moderate to severe at a single gene level, and indicate reduced regulation as well as activity. The objective quantitative assessment of VUS scoring and gene overexpression in the context of functional interactions and pathways provides insights for biology, oncology, and precision medicine.

Identification of evolutionarily stable functional and immunogenic sites across the SARS-CoV-2 proteome and greater coronavirus family

MOTIVATION

Since the first recognized case of COVID-19, more than 100 million people have been infected worldwide. Global efforts in drug and vaccine development to fight the disease have yielded vaccines and drug candidates to cure COVID-19. However, the spread of SARS-CoV-2 variants threatens the continued efficacy of these treatments. In order to address this, we interrogate the evolutionary history of the entire SARS-CoV-2 proteome to identify evolutionarily conserved functional sites that can inform the search for treatments with broader coverage across the coronavirus family.

RESULTS

Combining coronavirus family sequence information with the mutations observed in the current COVID-19 outbreak, we systematically and comprehensively define evolutionarily stable sites that may provide useful drug and vaccine targets and which are less likely to be compromised by the emergence of new virus strains. Several experimentally validated effective drugs interact with these proposed target sites. In addition, the same evolutionary information can prioritize cross reactive antigens that are useful in directing multi-epitope vaccine strategies to illicit broadly neutralizing immune responses to the betacoronavirus family. Although the results are focused on SARS-CoV-2, these approaches stem from evolutionary principles that are agnostic to the organism or infective agent.

AVAILABILITY AND IMPLEMENTATION

The results of this work are made interactively available at http://cov.lichtargelab.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Uncovering DNA-PKcs ancient phylogeny, unique sequence motifs and insights for human disease

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

Targeting SARS-CoV-2 Nsp3 macrodomain structure with insights from human poly(ADP-ribose) glycohydrolase (PARG) structures with inhibitors

Arrival of the novel SARS-CoV-2 has launched a worldwide effort to identify both pre-approved and novel therapeutics targeting the viral proteome, highlighting the urgent need for efficient drug discovery strategies. Even with effective vaccines, infection is possible, and at-risk populations would benefit from effective drug compounds that reduce the lethality and lasting damage of COVID-19 infection. The CoV-2 MacroD-like macrodomain (Mac1) is implicated in viral pathogenicity by disrupting host innate immunity through its mono (ADP-ribosyl) hydrolase activity, making it a prime target for antiviral therapy. We therefore solved the structure of CoV-2 Mac1 from non-structural protein 3 (Nsp3) and applied structural and sequence-based genetic tracing, including newly determined A. pompejana MacroD2 and GDAP2 amino acid sequences, to compare and contrast CoV-2 Mac1 with the functionally related human DNA-damage signaling factor poly (ADP-ribose) glycohydrolase (PARG). Previously, identified targetable features of the PARG active site allowed us to develop a pharmacologically useful PARG inhibitor (PARGi). Here, we developed a focused chemical library and determined 6 novel PARGi X-ray crystal structures for comparative analysis. We applied this knowledge to discovery of CoV-2 Mac1 inhibitors by combining computation and structural analysis to identify PARGi fragments with potential to bind the distal-ribose and adenosyl pockets of the CoV-2 Mac1 active site. Scaffold development of these PARGi fragments has yielded two novel compounds, PARG-345 and PARG-329, that crystallize within the Mac1 active site, providing critical structure-activity data and a pathway for inhibitor optimization. The reported structural findings demonstrate ways to harness our PARGi synthesis and characterization pipeline to develop CoV-2 Mac1 inhibitors targeting the ADP-ribose active site. Together, these structural and computational analyses reveal a path for accelerating development of antiviral therapeutics from pre-existing drug optimization pipelines.

A method to delineate de novo missense variants across pathways prioritizes genes linked to autism

Genotype-phenotype relationships shape health and population fitness but remain difficult to predict and interpret. Here, we apply an evolutionary action method to de novo missense variants in whole-exome sequences of individuals with autism spectrum disorder (ASD) to unravel genes and pathways connected to ASD. Evolutionary action predicts the impact of missense variants on protein function by measuring the fitness effect based on phylogenetic distances and substitution odds in homologous gene sequences. By examining de novo missense variants in 2384 individuals with ASD (probands) compared to matched siblings without ASD, we found missense variants in 398 genes representing 23 pathways that were biased toward higher evolutionary action scores than expected by random chance; these pathways were involved in axonogenesis, synaptic transmission, and neurodevelopment. The predicted fitness impact of de novo and inherited missense variants in candidate genes correlated with the IQ of individuals with ASD, even for new gene candidates. Taking an evolutionary action method, we detected those missense variants most likely to contribute to ASD pathogenesis and elucidated their phenotypic impact. This approach could be applied to integrate missense variants across a patient cohort to identify genes contributing to a shared phenotype in other complex diseases.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations

MOTIVATION

Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated.

RESULTS

Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance.

AVAILABILITY AND IMPLEMENTATION

All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structure and evolutionary trace-assisted screening of a residue swapping the substrate ambiguity and chiral specificity in an esterase

Our understanding of enzymes with high substrate ambiguity remains limited because their large active sites allow substrate docking freedom to an extent that seems incompatible with stereospecificity. One possibility is that some of these enzymes evolved a set of evolutionarily fitted sequence positions that stringently allow switching substrate ambiguity and chiral specificity. To explore this hypothesis, we targeted for mutation a serine ester hydrolase (EH3) that exhibits an impressive 71-substrate repertoire but is not stereospecific (e.e. 50%). We used structural actions and the computational evolutionary trace method to explore specificity-swapping sequence positions and hypothesized that position I244 was critical. Driven by evolutionary action analysis, this position was substituted to leucine, which together with isoleucine appears to be the amino acid most commonly present in the closest homologous sequences (max. identity, ca. 67.1%), and to phenylalanine, which appears in distant homologues. While the I244L mutation did not have any functional consequences, the I244F mutation allowed the esterase to maintain a remarkable 53-substrate range while gaining stereospecificity properties (e.e. 99.99%). These data support the possibility that some enzymes evolve sequence positions that control the substrate scope and stereospecificity. Such residues, which can be evolutionarily screened, may serve as starting points for further designing substrate-ambiguous, yet chiral-specific, enzymes that are greatly appreciated in biotechnology and synthetic chemistry.

A Conserved Motif in Intracellular Loop 1 Stabilizes the Outward-Facing Conformation of TmrAB

The ATP binding cassette (ABC) family of transporters moves small molecules (lipids, sugars, peptides, drugs, nutrients) across membranes in nearly all organisms. Transport activity requires conformational switching between inward-facing and outward-facing states driven by ATP-dependent dimerization of two nucleotide binding domains (NBDs). The mechanism that connects ATP binding and hydrolysis in the NBDs to conformational changes in a substrate binding site in the transmembrane domains (TMDs) is currently an outstanding question. Here we use sequence coevolution analyses together with biochemical characterization to investigate the role of a highly conserved region in intracellular loop 1 we define as the GRD motif in coordinating domain rearrangements in the heterodimeric peptide exporter from Thermus thermophilus, TmrAB. Mutations in the GRD motif alter ATPase activity as well as transport. Disulfide crosslinking, evolutionary trace, and evolutionary coupling analysis reveal that these effects are likely due to the destabilization of a network in which the GRD motif in TmrA bridges residues of the Q-loop, X-loop, and ABC motif in the NBDs to residues in the TmrAB peptide substrate binding site, thus providing an avenue for conformational coupling. We further find that disruption of this network in TmrA versus TmrB has different functional consequences, hinting at an intrinsic asymmetry in heterodimeric ABC transporters extending beyond that of the NBDs. These results support a mechanism in which the GRD motifs help coordinate a transition to an outward open conformation, and each half of the transporter likely plays a different role in the conformational cycle of TmrAB.

An efficient chemical screening method for structure-based inhibitors to nucleic acid enzymes targeting the DNA repair-replication interface and SARS CoV-2

We present a Chemistry and Structure Screen Integrated Efficiently (CASSIE) approach (named for Greek prophet Cassandra) to design inhibitors for cancer biology and pathogenesis. CASSIE provides an effective path to target master keys to control the repair-replication interface for cancer cells and SARS CoV-2 pathogenesis as exemplified here by specific targeting of Poly(ADP-ribose) glycohydrolase (PARG) and ADP-ribose glycohydrolase ARH3 macrodomains plus SARS CoV-2 nonstructural protein 3 (Nsp3) Macrodomain 1 (Mac1) and Nsp15 nuclease. As opposed to the classical massive effort employing libraries with large numbers of compounds against single proteins, we make inhibitor design for multiple targets efficient. Our compact, chemically diverse, 5000 compound Goldilocks (GL) library has an intermediate number of compounds sized between fragments and drugs with predicted favorable ADME (absorption, distribution, metabolism, and excretion) and toxicological profiles. Amalgamating our core GL library with an approved drug (AD) library, we employ a combined GLAD library virtual screen, enabling an effective and efficient design cycle of ranked computer docking, top hit biophysical and cell validations, and defined bound structures using human proteins or their avatars. As new drug design is increasingly pathway directed as well as molecular and mechanism based, our CASSIE approach facilitates testing multiple related targets by efficiently turning a set of interacting drug discovery problems into a tractable medicinal chemistry engineering problem of optimizing affinity and ADME properties based upon early co-crystal structures. Optimization efforts are made efficient by a computationally-focused iterative chemistry and structure screen. Thus, we herein describe and apply CASSIE to define prototypic, specific inhibitors for PARG vs distinct inhibitors for the related macrodomains of ARH3 and SARS CoV-2 Nsp3 plus the SARS CoV-2 Nsp15 RNA nuclease.

Understanding substrate binding and the role of gatekeeping residues in PigC access tunnels

Prodigiosin ligase PigC has been engineered by semi-rational design to accept short chain-pyrroles, providing molecular understanding of access tunnels and the substrate-binding pocket.

Refinement of pairwise potentials via logistic regression to score protein-protein interactions

Protein-protein interactions (PPIs) are ubiquitous and functionally of great importance in biological systems. Hence, the accurate prediction of PPIs by protein-protein docking and scoring tools is highly desirable in order to characterize their structure and biological function. Ab initio docking protocols are divided into the sampling of docking poses to produce at least one near-native structure, and then to evaluate the vast candidate structures by scoring. Concurrent development in both sampling and scoring is crucial for the deployment of protein-protein docking software. In the present work, we apply a machine learning model on pairwise potentials to refine the task of protein quaternary structure native structure detection among decoys. A decoy set was featurized using the Knowledge and Empirical Combined Scoring Algorithm 2 (KECSA2) pairwise potential. The highly unbalanced decoy set was then balanced using a comparison concept between native and decoy structures. The resultant comparison descriptors were used to train a logistic regression (LR) classifier. The LR model yielded the optimal performance for native detection among decoys compared with conventional scoring functions, while exhibiting lesser performance for the detection of low root mean square deviation decoy structures. Its deployment on an independent benchmark set confirms that the scoring function performs competitively relative to other scoring functions. The scripts used are available at https://github.com/TanemuraKiyoto/PPI-native-detection-via-LR.

An Evolutionary Marker of the Ribokinase Superfamily Is Responsible for Zinc-Mediated Regulation of Human Pyridoxal Kinase

The ribokinase superfamily catalyzes the phosphorylation of a vast diversity of substrates, and its members are characterized by the conservation of a common structural fold along with highly conserved sequence motifs responsible for phosphoryl transfer (GXGD) and stabilization of the metal-nucleotide complex (NXXE). Recently, a third motif (HXE) exclusive from ADP-dependent enzymes was identified, with its glutamic acid participating in water-mediated interactions with the metal-nucleotide complex and in stabilization of the ternary complex during catalysis. In this work, we bioinformatically determine that the aspartic acid of another motif (DPV), exclusively found in hydroxyethyl thiazole (THZK), hydroxymethyl pyrimidine (HMPK) and pyridoxal kinases (PLK), is structurally equivalent to the acidic residue in the HXE motif. Moreover, this residue is highly conserved among all ribokinase superfamily members. To determine whether the functional role of the DPV motif is similar to the HXE motif, we employed molecular dynamics simulations using crystal structures of phosphoryl donor substrate-complexed THZK and PLK, showing that its aspartic acid participated in water-mediated or direct interactions with the divalent metal of the metal-nucleotide complex. Lastly, enzyme kinetic assays on human PLK, an enzyme that utilizes zinc, showed that site-directed mutagenesis of the aspartic acid from the DPV motif abolishes the inhibition of this enzyme by increasing free zinc concentrations. Altogether, our results highlight that the DPV and HXE motifs are evolutionary markers of the functional and structural divergence of the ribokinase superfamily and evidence the role of the DPV motif in the interaction with both free and nucleotide-complexed divalent metals in the binding site of these enzymes.

Sequence-structure-function relationships in class I MHC: A local frustration perspective

Class I Major Histocompatibility Complex (MHC) binds short antigenic peptides with the help of Peptide Loading Complex (PLC), and presents them to T-cell Receptors (TCRs) of cytotoxic T-cells and Killer-cell Immunglobulin-like Receptors (KIRs) of Natural Killer (NK) cells. With more than 10000 alleles, human MHC (Human Leukocyte Antigen, HLA) is the most polymorphic protein in humans. This allelic diversity provides a wide coverage of peptide sequence space, yet does not affect the three-dimensional structure of the complex. Moreover, TCRs mostly interact with HLA in a common diagonal binding mode, and KIR-HLA interaction is allele-dependent. With the aim of establishing a framework for understanding the relationships between polymorphism (sequence), structure (conserved fold) and function (protein interactions) of the human MHC, we performed here a local frustration analysis on pMHC homology models covering 1436 HLA I alleles. An analysis of local frustration profiles indicated that (1) variations in MHC fold are unlikely due to minimally-frustrated and relatively conserved residues within the HLA peptide-binding groove, (2) high frustration patches on HLA helices are either involved in or near interaction sites of MHC with the TCR, KIR, or tapasin of the PLC, and (3) peptide ligands mainly stabilize the F-pocket of HLA binding groove.

An Evolutionary Trace method defines functionally important bases and sites common to RNA families

Functional non-coding (fnc)RNAs are nucleotide sequences of varied lengths, structures, and mechanisms that ubiquitously influence gene expression and translation, genome stability and dynamics, and human health and disease. Here, to shed light on their functional determinants, we seek to exploit the evolutionary record of variation and divergence read from sequence comparisons. The approach follows the phylogenetic Evolutionary Trace (ET) paradigm, first developed and extensively validated on proteins. We assigned a relative rank of importance to every base in a study of 1070 functional RNAs, including the ribosome, and observed evolutionary patterns strikingly similar to those seen in proteins, namely, (1) the top-ranked bases clustered in secondary and tertiary structures. (2) In turn, these clusters mapped functional regions for catalysis, binding proteins and drugs, post-transcriptional modification, and deleterious mutations. (3) Moreover, the quantitative quality of these clusters correlated with the identification of functional regions. (4) As a result of this correlation, smoother structural distributions of evolutionary important nucleotides improved functional site predictions. Thus, in practice, phylogenetic analysis can broadly identify functional determinants in RNA sequences and functional sites in RNA structures, and reveal details on the basis of RNA molecular functions. As example of application, we report several previously undocumented and potentially functional ET nucleotide clusters in the ribosome. This work is broadly relevant to studies of structure-function in ribonucleic acids. Additionally, this generalization of ET shows that evolutionary constraints among sequence, structure, and function are similar in structured RNA and proteins. RNA ET is currently available as part of the ET command-line package, and will be available as a web-server.

Traditionally, RNA has been delegated to the role of an intermediate between DNA and proteins. However, we now recognize that RNAs are broadly functional beyond their role in translation, and that a number of diverse classes exist. Because functional, non-coding RNAs are prevalent in biology and impact human health, it is important to better understand their functional determinants. However, the classical solution to this problem, targeted mutagenesis, is time-consuming and scales poorly. We propose an alternative computational approach to this problem, the Evolutionary Trace method. Previously developed and validated for proteins, Evolutionary Trace examines evolutionary history of a molecule and predicts evolutionarily important residues in the sequence. We apply Evolutionary Trace to a set of diverse RNAs, and find that the evolutionarily important nucleotides cluster on the three-dimensional structure, and that these clusters closely overlap functional sites. We also find that the clustering property can be used to refine and improve predictions. These findings are in close agreement with our observations of Evolutionary Trace in proteins, and suggest that structured functional RNAs and proteins evolve under similar constraints. In practice, the approach is to be used by RNA researches seeking insight into their molecule of interest, and the Evolutionary Trace program, along with a working example, is available at https://github.com/LichtargeLab/RNA_ET_ms.

Deep Analysis of Residue Constraints (DARC): identifying determinants of protein functional specificity

Protein functional constraints are manifest as superfamily and functional-subgroup conserved residues, and as pairwise correlations. Deep Analysis of Residue Constraints (DARC) aids the visualization of these constraints, characterizes how they correlate with each other and with structure, and estimates statistical significance. This can identify determinants of protein functional specificity, as we illustrate for bacterial DNA clamp loader ATPases. These load ring-shaped sliding clamps onto DNA to keep polymerase attached during replication and contain one δ, three γ, and one δ’ AAA+ subunits semi-circularly arranged in the order δ-γ₁-γ₂-γ₃-δ’. Only γ is active, though both γ and δ’ functionally influence an adjacent γ subunit. DARC identifies, as functionally-congruent features linking allosterically the ATP, DNA, and clamp binding sites: residues distinctive of γ and of γ/δ’ that mutually interact in trans, centered on the catalytic base; several γ/δ’-residues and six γ/δ’-covariant residue pairs within the DNA binding N-termini of helices α2 and α3; and γ/δ’-residues associated with the α2 C-terminus and the clamp-binding loop. Most notable is a trans-acting γ/δ’ hydroxyl group that 99% of other AAA+ proteins lack. Mutation of this hydroxyl to a methyl group impedes clamp binding and opening, DNA binding, and ATP hydrolysis—implying a remarkably clamp-loader-specific function.

Multiple protein-DNA interfaces unravelled by evolutionary information, physico-chemical and geometrical properties

Interactions between proteins and nucleic acids are at the heart of many essential biological processes. Despite increasing structural information about how these interactions may take place, our understanding of the usage made of protein surfaces by nucleic acids is still very limited. This is in part due to the inherent complexity associated to protein surface deformability and evolution. In this work, we present a method that contributes to decipher such complexity by predicting protein-DNA interfaces and characterizing their properties. It relies on three biologically and physically meaningful descriptors, namely evolutionary conservation, physico-chemical properties and surface geometry. We carefully assessed its performance on several hundreds of protein structures and compared it to several machine-learning state-of-the-art methods. Our approach achieves a higher sensitivity compared to the other methods, with a similar precision. Importantly, we show that it is able to unravel ‘hidden’ binding sites by applying it to unbound protein structures and to proteins binding to DNA via multiple sites and in different conformations. It is also applicable to the detection of RNA-binding sites, without significant loss of performance. This confirms that DNA and RNA-binding sites share similar properties. Our method is implemented as a fully automated tool, JETDNA2, freely accessible at: http://www.lcqb.upmc.fr/JET2DNA. We also provide a new dataset of 187 protein-DNA complex structures, along with a subset of 82 associated unbound structures. The set represents the largest body of high-resolution crystallographic structures of protein-DNA complexes, use biological protein assemblies as DNA-binding units, and covers all major types of protein-DNA interactions. It is available at: http://www.lcqb.upmc.fr/PDNAbenchmarks.

Protein-DNA interactions are essential to living organisms and their impairment is associated to many diseases. For these reasons, they have become increasingly important therapeutic targets. Experimental structure determination has revealed different binding motifs and modes, associated to different functions. Yet, the available structural data gives us only a glimpse of the multiplicity and complexity of protein surface usage by DNA. In this work, we use a three-layer model to describe and predict DNA-binding sites at protein surfaces. Given a protein, we consider the way its residues are conserved through evolution, their physico-chemical properties and geometrical shapes to decrypt its surface. We are able to detect a large portion of interacting residues with good precision, even when they are ‘hidden’ by conformational changes. We highlight cases where one protein binds DNA via distinct regions to perform different functions. We are able to uncover the alternative binding sites and relate their properties with their specific roles. Our work can help guiding mutagenesis experiments and the development of new drugs specifically targeting one site while limiting possible side effects.

Deleterious Non-Synonymous Single Nucleotide Polymorphism Predictions on Human Transcription Factors

Transcription factors (TFs) are the major components of human gene regulation. In particular, they bind onto specific DNA sequences and regulate neighborhood genes in different tissues at different developmental stages. Non-synonymous single nucleotide polymorphisms on its protein-coding sequences could result in undesired consequences in human. Therefore, it is necessary to develop methods for predicting any abnormality among those non-synonymous single nucleotide polymorphisms. To address it, we have developed and compared different strategies to predict deleterious non-synonymous single nucleotide polymorphisms (also known as missense mutations) on the protein-coding sequences of human TFs. Taking advantage of evolutionary conservation signals, we have developed and compared different classifiers with different feature sets as computed from different evolutionarily related sequence collections. The results indicate that the classic ensemble algorithm, Adaboost with decision stumps, with orthologous sequence collection, has performed the best (namely, TFmedic). We have further compared TFmedic with other state-of-the-arts methods (i.e., PolyPhen-2 and SIFT) on PolyPhen-2's own datasets, demonstrating that TFmedic can outperform the others. As applications, we have further applied TFmedic to all possible missense mutations on all human transcription factors; the proteome-wide results reveal interesting insights, consistent with the existing physiochemical knowledge. A case study with the actual 3D structure is conducted, revealing how TFmedic can be contributed to protein-DNA binding complex studies.

ConSurf-DB: An accessible repository for the evolutionary conservation patterns of the majority of PDB proteins

Patterns observed by examining the evolutionary relationships among proteins of common origin can reveal the structural and functional importance of specific residue positions. In particular, amino acids that are highly conserved (i.e., their positions evolve at a slower rate than other positions) are particularly likely to be of biological importance, for example, for ligand binding. ConSurf is a bioinformatics tool for accurately estimating the evolutionary rate of each position in a protein family. Here we introduce a new release of ConSurf‐DB, a database of precalculated ConSurf evolutionary conservation profiles for proteins of known structure. ConSurf‐DB provides high‐accuracy estimates of the evolutionary rates of the amino acids in each protein. A reliable estimate of a query protein's evolutionary rates depends on having a sufficiently large number of effective homologues (i.e., nonredundant yet sufficiently similar). With current sequence data, ConSurf‐DB covers 82% of the PDB proteins. It will be updated on a regular basis to ensure that coverage remains high—and that it might even increase. Much effort was dedicated to improving the user experience. The repository is available at https://consurfdb.tau.ac.il/.

Broader audience

By comparing a protein to other proteins of similar origin, it is possible to determine the extent to which each amino acid position in the protein evolved slowly or rapidly. A protein's evolutionary profile can provide valuable insights: For example, amino acid positions that are highly conserved (i.e., evolved slowly) are particularly likely to be of structural and/or functional importance, for example, for ligand binding and catalysis. We introduce here a new and improved version of ConSurf‐DB, a continually updated database that provides precalculated evolutionary profiles of proteins with known structure.

Mutational Landscape of the BAP1 Locus Reveals an Intrinsic Control to Regulate the miRNA Network and the Binding of Protein Complexes in Uveal Melanoma

The BAP1 (BRCA1-associated protein 1) gene is associated with a variety of human cancers. With its gene product being a nuclear ubiquitin carboxy-terminal hydrolase with deubiquitinase activity, BAP1 acts as a tumor suppressor gene with potential pleiotropic effects in multiple tumor types. Herein, we focused specifically on uveal melanoma (UM) in which BAP1 mutations are associated with a metastasizing phenotype and decreased survival rates. We identified the ubiquitin carboxyl hydrolase (UCH) domain as a major hotspot region for the pathogenic mutations with a high evolutionary action (EA) score. This also includes the mutations at conserved catalytic sites and the ones overlapping with the phosphorylation residues. Computational protein interaction studies revealed that distant BAP1-associated protein complexes (FOXK2, ASXL1, BARD1, BRCA1) could be directly impacted by this mutation paradigm. We also described the conformational transition related to BAP1-BRCA-BARD1 complex, which may pose critical implications for mutations, especially at the docking interfaces of these three proteins. The mutations affect - independent of being somatic or germline - the binding affinity of miRNAs embedded within the BAP1 locus, thereby altering the unique regulatory network. Apart from UM, BAP1 gene expression and survival associations were found to be predictive for the prognosis in several (n = 29) other cancer types. Herein, we suggest that although BAP1 is conceptually a driver gene in UM, it might contribute through its interaction partners and its regulatory miRNA network to various aspects of cancer. Taken together, these findings will pave the way to evaluate BAP1 in a variety of other human cancers with a shared mutational spectrum.

Short-Term Low Temperature Induces Nitro-Oxidative Stress that Deregulates the NADP-Malic Enzyme Function by Tyrosine Nitration in Arabidopsis thaliana

Low temperature (LT) negatively affects plant growth and development via the alteration of the metabolism of reactive oxygen and nitrogen species (ROS and RNS). Among RNS, tyrosine nitration, the addition of an NO₂ group to a tyrosine residue, can modulate reduced nicotinamide-dinucleotide phosphate (NADPH)-generating systems and, therefore, can alter the levels of NADPH, a key cofactor in cellular redox homeostasis. NADPH also acts as an indispensable electron donor within a wide range of enzymatic reactions, biosynthetic pathways, and detoxification processes, which could affect plant viability. To extend our knowledge about the regulation of this key cofactor by this nitric oxide (NO)-related post-translational modification, we analyzed the effect of tyrosine nitration on another NADPH-generating enzyme, the NADP-malic enzyme (NADP-ME), under LT stress. In Arabidopsis thaliana seedlings exposed to short-term LT (4 °C for 48 h), a 50% growth reduction accompanied by an increase in the content of superoxide, nitric oxide, and peroxynitrite, in addition to diminished cytosolic NADP-ME activity, were found. In vitro assays confirmed that peroxynitrite inhibits cytosolic NADP-ME2 activity due to tyrosine nitration. The mass spectrometric analysis of nitrated NADP-ME2 enabled us to determine that Tyr-73 was exclusively nitrated to 3-nitrotyrosine by peroxynitrite. The in silico analysis of the Arabidopsis NADP-ME2 protein sequence suggests that Tyr73 nitration could disrupt the interactions between the specific amino acids responsible for protein structure stability. In conclusion, the present data show that short-term LT stress affects the metabolism of ROS and RNS, which appears to negatively modulate the activity of cytosolic NADP-ME through the tyrosine nitration process.

Assessment of blind predictions of the clinical significance of BRCA1 and BRCA2 variants

Testing for variation in BRCA1 and BRCA2 (commonly referred to as BRCA1/2), has emerged as a standard clinical practice and is helping countless women better understand and manage their heritable risk of breast and ovarian cancer. Yet the increased rate of BRCA1/2 testing has led to an increasing number of Variants of Uncertain Significance (VUS), and the rate of VUS discovery currently outpaces the rate of clinical variant interpretation. Computational prediction is a key component of the variant interpretation pipeline. In the CAGI5 ENIGMA Challenge, six prediction teams submitted predictions on 326 newly-interpreted variants from the ENIGMA Consortium. By evaluating these predictions against the new interpretations, we have gained a number of insights on the state of the art of variant prediction and specific steps to further advance this state of the art.

GEMME: a simple and fast global epistatic model predicting mutational effects

The systematic and accurate description of protein mutational landscapes is a question of utmost importance in biology, bioengineering, and medicine. Recent progress has been achieved by leveraging on the increasing wealth of genomic data and by modeling intersite dependencies within biological sequences. However, state-of-the-art methods remain time consuming. Here, we present Global Epistatic Model for predicting Mutational Effects (GEMME) (www.lcqb.upmc.fr/GEMME), an original and fast method that predicts mutational outcomes by explicitly modeling the evolutionary history of natural sequences. This allows accounting for all positions in a sequence when estimating the effect of a given mutation. GEMME uses only a few biologically meaningful and interpretable parameters. Assessed against 50 high- and low-throughput mutational experiments, it overall performs similarly or better than existing methods. It accurately predicts the mutational landscapes of a wide range of protein families, including viral ones and, more generally, of much conserved families. Given an input alignment, it generates the full mutational landscape of a protein in a matter of minutes. It is freely available as a package and a webserver at www.lcqb.upmc.fr/GEMME/.

CAGI5: Objective performance assessments of predictions based on the Evolutionary Action equation

Many computational approaches estimate the effect of coding variants, but their predictions often disagree with each other. These contradictions confound users and raise questions regarding reliability. Performance assessments can indicate the expected accuracy for each method and highlight advantages and limitations. The Critical Assessment of Genome Interpretation (CAGI) community aims to organize objective and systematic assessments: They challenge predictors on unpublished experimental and clinical data and assign independent assessors to evaluate the submissions. We participated in CAGI experiments as predictors, using the Evolutionary Action (EA) method to estimate the fitness effect of coding mutations. EA is untrained, uses homology information, and relies on a formal equation: The fitness effect equals the functional sensitivity to residue changes multiplied by the magnitude of the substitution. In previous CAGI experiments (between 2011 and 2016), our submissions aimed to predict the protein activity of single mutants. In 2018 (CAGI5), we also submitted predictions regarding clinical associations, folding stability, and matching genomic data with phenotype. For all these diverse challenges, we used EA to predict the fitness effect of variants, adjusted to specifically address each question. Our submissions had consistently good performance, suggesting that EA predicts reliably the effects of genetic variants.

Evolutionary Stability of Salt Bridges Hints Its Contribution to Stability of Proteins

The contribution of newly designed salt bridges to protein stabilization remains controversial even today. In order to solve this problem, we investigated salt bridges from two aspects: spatial distribution and evolutionary characteristics of salt bridges. Firstly, we analyzed spatial distribution of salt bridges in proteins, elucidating the basic requirements of forming salt bridges. Then, from an evolutionary point of view, the evolutionary characteristics of salt bridges as well as their neighboring residues were investigated in our study. The results demonstrate that charged residues appear more frequently than other neutral residues at certain positions of sequence even under evolutionary pressure, which are able to form electrostatic interactions that could increase the evolutionary stability of corresponding amino acid regions, enhancing their importance to stability of proteins. As a corollary, we conjectured that the newly designed salt bridges with more contribution to proteins, not only, are qualified spatial distribution of salt bridges, but also, are needed to further increase the evolutionary stability of corresponding amino acid regions. Based on analysis, the 8 mutations were accordingly constructed in the 1,4-α-glucan branching enzyme (EC 2.4.1.18, GBE) from Geobacillus thermoglucosidans STB02, of which 7 mutations improved thermostability of GBE. The enhanced thermostability of 7 mutations might be a result of additional salt bridges on residue positions that at least one of amino acids positions is conservative, improving their contribution of stabilization to proteins.

Decrypting protein surfaces by combining evolution, geometry, and molecular docking

The growing body of experimental and computational data describing how proteins interact with each other has emphasized the multiplicity of protein interactions and the complexity underlying protein surface usage and deformability. In this work, we propose new concepts and methods toward deciphering such complexity. We introduce the notion of interacting region to account for the multiple usage of a protein's surface residues by several partners and for the variability of protein interfaces coming from molecular flexibility. We predict interacting patches by crossing evolutionary, physicochemical and geometrical properties of the protein surface with information coming from complete cross-docking (CC-D) simulations. We show that our predictions match well interacting regions and that the different sources of information are complementary. We further propose an indicator of whether a protein has a few or many partners. Our prediction strategies are implemented in the dynJET² algorithm and assessed on a new dataset of 262 protein on which we performed CC-D. The code and the data are available at: http://www.lcqb.upmc.fr/dynJET2/.

The Subtle Trade-Off between Evolutionary and Energetic Constraints in Protein-Protein Interactions

Life machinery, although overwhelmingly complex, is rooted on a rather limited number of molecular processes. One of the most important is protein-protein interaction. Metabolic regulation, protein folding control, and cellular motility are examples of processes based on the fine-tuned interaction of several protein partners. The region on the protein surface devoted to the recognition of a specific partner is essential for the function of the protein and is, therefore, likely to be conserved during evolution. On the other hand, the physical chemistry of amino acids underlies the mechanism of interactions. Both evolutionary and energetic constraints can then be used to build scoring functions capable of recognizing interaction sites. Our working hypothesis is that residues within the interaction interface tend at the same time to be evolutionarily conserved (to preserve their function) and to provide little contribution to the internal stabilization of the structure of their cognate protein, to facilitate conformational adaptation to the partner. Here, we show that for some classes of protein partners (for example, those involved in signal transduction and in enzymes) evolutionary constraints play the key role in defining the interaction surface. In contrast, energetic constraints emerge as more important in protein partners involved in immune response, in inhibitor proteins, and in structural proteins. Our results indicate that a general-purpose scoring function for protein-protein interaction should not be agnostic of the biological function of the partners.

Learning protein constitutive motifs from sequence data

Statistical analysis of evolutionary-related protein sequences provides information about their structure, function, and history. We show that Restricted Boltzmann Machines (RBM), designed to learn complex high-dimensional data and their statistical features, can efficiently model protein families from sequence information. We here apply RBM to 20 protein families, and present detailed results for two short protein domains (Kunitz and WW), one long chaperone protein (Hsp70), and synthetic lattice proteins for benchmarking. The features inferred by the RBM are biologically interpretable: they are related to structure (residue-residue tertiary contacts, extended secondary motifs (α-helixes and β-sheets) and intrinsically disordered regions), to function (activity and ligand specificity), or to phylogenetic identity. In addition, we use RBM to design new protein sequences with putative properties by composing and 'turning up' or 'turning down' the different modes at will. Our work therefore shows that RBM are versatile and practical tools that can be used to unveil and exploit the genotype-phenotype relationship for protein families.

FARS2 deficiency; new cases, review of clinical, biochemical, and molecular spectra, and variants interpretation based on structural, functional, and evolutionary significance

An increasing number of mitochondrial diseases are found to be caused by pathogenic variants in nuclear encoded mitochondrial aminoacyl-tRNA synthetases. FARS2 encodes mitochondrial phenylalanyl-tRNA synthetase (mtPheRS) which transfers phenylalanine to its cognate tRNA in mitochondria. Since the first case was reported in 2012, a total of 21 subjects with FARS2 deficiency have been reported to date with a spectrum of disease severity that falls between two phenotypes; early onset epileptic encephalopathy and a less severe phenotype characterized by spastic paraplegia. In this report, we present an additional 15 individuals from 12 families who are mostly Arabs homozygous for the pathogenic variant Y144C, which is associated with the more severe early onset phenotype. The total number of unique pathogenic FARS2 variants known to date is 21 including three different partial gene deletions reported in four individuals. Except for the large deletions, all variants but two (one in-frame deletion of one amino acid and one splice-site variant) are missense. All large deletions and the single splice-site variant are in trans with a missense variant. This suggests that complete loss of function may be incompatible with life. In this report, we also review structural, functional, and evolutionary significance of select FARS2 pathogenic variants reported here.

Endogenous hydrogen sulfide (H2S) is up-regulated during sweet pepper (Capsicum annuum L.) fruit ripening. In vitro analysis shows that NADP-dependent isocitrate dehydrogenase (ICDH) activity is inhibited by H2S and NO

Like nitric oxide (NO), hydrogen sulfide (H₂S) has been recognized as a new gasotransmitter which plays an important role as a signaling molecule in many physiological processes in higher plants. Although fruit ripening is a complex process associated with the metabolism of reactive oxygen species (ROS) and nitrogen oxygen species (RNS), little is known about the potential involvement of endogenous H₂S. Using sweet pepper (Capsicum annuum L.) as a model non-climacteric fruit during the green and red ripening stages, we studied endogenous H₂S content and cytosolic l-cysteine desulfhydrase (L-DES) activity which increased by 14% and 28%, respectively, in red pepper fruits. NADPH is a redox compound and key cofactor required for cell growth, proliferation and detoxification. We studied the NADPH-regenerating enzyme, NADP-isocitrate dehydrogenase (NADP-ICDH), whose activity decreased by 34% during ripening. To gain a better understanding of its potential regulation by H₂S, we obtained a 50-75% ammonium sulfate-enriched protein fraction containing the NADP-ICDH protein; with the aid of in vitro assays in the presence of H₂S, we observed that 2 and 10 mM NaHS used as H₂S donors resulted in a decrease of up to 36% and 45%, respectively, in NADP-ICDH activity, which was unaffected by reduced glutathione (GSH). On the other hand, peroxynitrite (ONOO^-), S-nitrosocyteine (CysNO) and DETA-NONOate, with the last two acting as NO donors, also inhibited NADP-ICDH activity. In silico analysis of the tertiary structure of sweet pepper NADP-ICDH activity (UniProtKB ID A0A2G2Y555) suggests that residues Cys133 and Tyr450 are the most likely potential targets for S-nitrosation and nitration, respectively. Taken together, the data reveal that the increase in the H₂S production capacity of red fruits is due to higher L-DES activity during non-climacteric pepper fruit ripening. In vitro assays appear to show that H₂S inhibits NADP-ICDH activity, thus suggesting that this enzyme may be regulated by persulfidation, as well as by S-nitrosation and nitration. NO and H₂S may therefore regulate NADPH production and consequently cellular redox status during pepper fruit ripening.

Outcome analysis of Phase I trial patients with metastatic KRAS and/or TP53 mutant non-small cell lung cancer

KRAS and TP53 mutations, which are the most common genetic drivers of tumorigenesis, are still considered undruggable targets. Therefore, we analyzed these genetic aberrations in metastatic non-small cell lung cancer (NSCLC) for the development of potential therapeutics. One hundred eighty-five consecutive patients with metastatic NSCLC in a phase 1 trial center were included. Their genomic aberrations, clinical characteristics, survivals, and phase 1 trial therapies were analyzed. About 10%, 18%, 36%, and 36% of the patients had metastatic KRAS+/TP53+, KRAS+/TP53-,KRAS-/TP53+, and KRAS-/TP53- NSCLC, respectively. The most common concurrent genetic aberrations beside KRAS and/or TP53 (>5%) were KIT, epidermal growth factor receptor, PIK3CA, c-MET, BRAF, STK11, ATM, CDKN2A, and APC. KRAS+/TP53+ NSCLC did not respond well to the phase 1 trial therapy and was associated with markedly worse progression-free (PFS) and overall (OS) survivals than the other three groups together. KRAS hotspot mutations at locations other than codon G12 were associated with considerably worse OS than those at this codon. Gene aberration-matched therapy produced prolonged PFS and so was anti-angiogenesis in patients with TP53 mutations. Introduction of the evolutionary action score system of TP53 missense mutations enabled us to identify a subgroup of NSCLC patients with low-risk mutant p53 proteins having a median OS duration of 64.5 months after initial diagnosis of metastasis. These data suggested that patients with metastatic dual KRAS+/TP53+ hotspot-mutant NSCLC had poor clinical outcomes. Further analysis identified remarkably prolonged survival in patients with low-risk mutant p53 proteins, which warrants confirmatory studies.

Co-evolution networks of HIV/HCV are modular with direct association to structure and function

Mutational correlation patterns found in population-level sequence data for the Human Immunodeficiency Virus (HIV) and the Hepatitis C Virus (HCV) have been demonstrated to be informative of viral fitness. Such patterns can be seen as footprints of the intrinsic functional constraints placed on viral evolution under diverse selective pressures. Here, considering multiple HIV and HCV proteins, we demonstrate that these mutational correlations encode a modular co-evolutionary structure that is tightly linked to the structural and functional properties of the respective proteins. Specifically, by introducing a robust statistical method based on sparse principal component analysis, we identify near-disjoint sets of collectively-correlated residues (sectors) having mostly a one-to-one association to largely distinct structural or functional domains. This suggests that the distinct phenotypic properties of HIV/HCV proteins often give rise to quasi-independent modes of evolution, with each mode involving a sparse and localized network of mutational interactions. Moreover, individual inferred sectors of HIV are shown to carry immunological significance, providing insight for guiding targeted vaccine strategies.

BALCONY: an R package for MSA and functional compartments of protein variability analysis

BACKGROUND

Here, we present an R package for entropy/variability analysis that facilitates prompt and convenient data extraction, manipulation and visualization of protein features from multiple sequence alignments. BALCONY can work with residues dispersed across a protein sequence and map them on the corresponding alignment of homologous protein sequences. Additionally, it provides several entropy and variability scores that indicate the conservation of each residue.

RESULTS

Our package allows the user to visualize evolutionary variability by locating the positions most likely to vary and to assess mutation candidates in protein engineering.

CONCLUSION

In comparison to other R packages BALCONY allows conservation/variability analysis in context of protein structure with linkage of the appropriate metrics with physicochemical features of user choice.

AVAILABILITY

CRAN project page: https://cran.r-project.org/package=BALCONY and our website: http://www.tunnelinggroup.pl/software/ for major platforms: Linux/Unix, Windows and Mac OS X.