News from the protein mutability landscape

Investigation of single nucleotide polymorphisms in differentially expressed genes and proteins reveals the genetic basis of skeletal muscle growth differences between Tibetan and Large White pigs

OBJECTIVE

Skeletal muscle growth is an important economic trait for meat production, with notable differences between Tibetan pigs (TIBPs, a slow-growing breed) and Large White pigs (LWPs, a fast-growing breed). However, the genetic underpinnings of this disparity remain unclear.

METHODS

In the current study, we integrated differentially expressed genes (DEGs) and proteins (DEPs) from 60-day-old embryonic muscle tissue, along with whole-genome single nucleotide polymorphisms (SNPs) displaying absolute allele frequency differences (ΔAF) of 0.5 or more between the TIBP and LWP breeds, to unravel the genetic factors influencing skeletal muscle growth.

RESULTS

Our analysis revealed 3,499 DEGs and 628 DEPs with SNPs having a ΔAF equal to or greater than 0.5. Further functional analysis identified 145 DEGs and 23 DEPs involved in biological processes related to skeletal muscle development, and 22 DEGs and 3 DEPs implicated in the mechanistic target of rapamycin kinase signaling pathway, which is known for positively regulating protein synthesis. Among these genes, several DEGs and DEPs, enriched with TIPB-specific SNPs in regulatory or/and coding regions, showed marked ΔAF between the TIBP and LWP breeds, including MYF5, MYOF, ASB2, PDE9A, SDC1, PDGFRA, MYOM2, ACVR1, ZIC3, COL11A1, TGFBR1, EDNRA, TGFB2, PDE4D, PGAM2, GRK2, SCN4B, CACNA1S, MYL4, IGF1, and FOXO1. Additionally, genes such as CAPN3, MYOM2, and PGAM2, identified as both DEPs and DEGs related to skeletal muscle development, contained multiple TIBP-specific and LWP-predominant SNPs in regulatory and/or coding regions, underscoring significant ΔAF differences between the two breeds.

CONCLUSION

This comprehensive investigation of SNPs in DEGs and DEPs identified a significant number of SNPs and genes related to skeletal muscle development during the prenatal stage. These findings not only shed light on potential causal genes for muscle divergence between the TIBP and LWP breeds but also offer valuable insights for pig breeding strategies aimed at enhancing meat production.

Determinants of Unfolding Cooperativity and Binding Are Decoupled in a DNA Binding Domain

The relative magnitudes of noncovalent stabilization energies or the coupling free energies in folded proteins are anisotropically distributed, uniquely influencing folding and functional behaviors. In this regard, the fructose repressor (FruR) DBD belonging to the LacR repressor family harbors a three-residue insertion─KQY─between the canonical second and third helices. This sequence insertion promotes a strong Tyr-Tyr stacking interaction that is not observed in related homologues. Combining experiments with simulations, we show that the Tyr-Tyr stacking contributes to a decoupled unfolding due to the localization of a large part of the stabilization energy in this specific structural region. This leads to melting temperatures from different probes spanning nearly 10 K, while concomitantly stabilizing a partially structured intermediate state. Disruption of the aromatic stacking interaction via an alanine mutation promotes a molten-globular state whose native ensemble is replete with non-native interactions while displaying enhanced thermodynamic fluctuations and minimal calorimetric cooperativity. Surprisingly, the molten-globular variant of FruR DBD binds to the operator site on DNA with an affinity similar to that of the wild-type but with altered secondary-structure characteristics in the bound state, underscoring the chaperone-like role of DNA through its large negative electrostatic potential. FruR DBD thus appears to be at the verge of disorder as expected of an entropically destabilizing three-residue insertion but is rescued by the aromatic stacking interaction that distinctly dictates the finer details of stability, cooperativity, and binding.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

Rendering protein mutation movies with MutAmore

BACKGROUND

The success of AlphaFold2 in reliable protein three-dimensional (3D) structure prediction, assists the move of structural biology toward studies of protein dynamics and mutational impact on structure and function. This transition needs tools that qualitatively assess alternative 3D conformations.

RESULTS

We introduce MutAmore, a bioinformatics tool that renders individual images of protein 3D structures for, e.g., sequence mutations into a visually intuitive movie format. MutAmore streamlines a pipeline casting single amino-acid variations (SAVs) into a dynamic 3D mutation movie providing a qualitative perspective on the mutational landscape of a protein. By default, the tool first generates all possible variants of the sequence reachable through SAVs (L*19 for proteins with L residues). Next, it predicts the structural conformation for all L*19 variants using state-of-the-art models. Finally, it visualizes the mutation matrix and produces a color-coded 3D animation. Alternatively, users can input other types of variants, e.g., from experimental structures.

CONCLUSION

MutAmore samples alternative protein configurations to study the dynamical space accessible from SAVs in the post-AlphaFold2 era of structural biology. As the field shifts towards the exploration of alternative conformations of proteins, MutAmore aids in the understanding of the structural impact of mutations by providing a flexible pipeline for the generation of protein mutation movies using current and future structure prediction models.

Protein engineering using mutability landscapes: Controlling site-selectivity of P450-catalyzed steroid hydroxylation

Directed evolution and rational design have been used widely in engineering enzymes for their application in synthetic organic chemistry and biotechnology. With stereoselectivity playing a crucial role in catalysis for the synthesis of valuable chemical and pharmaceutical compounds, rational design has not achieved such wide success in this specific area compared to directed evolution. Nevertheless, one bottleneck of directed evolution is the laborious screening efforts and the observed trade-offs in catalytic profiles. This has motivated researchers to develop more efficient protein engineering methods. As a prime approach, mutability landscaping avoids such trade-offs by providing more information of sequence-function relationships. Here, we describe an application of this efficient protein engineering method to improve the regio-/stereoselectivity and activity of P450BM3 for steroid hydroxylation, while keeping the mutagenesis libraries small so that they will require only minimal screening.

In Silico Examination of Single Nucleotide Missense Mutations in NHLH2, a Gene Linked to Infertility and Obesity

Continual advances in our understanding of the human genome have led to exponential increases in known single nucleotide variants. The characterization of each of the variants lags behind. For researchers needing to study a single gene, or multiple genes in a pathway, there must be ways to narrow down pathogenic variants from those that are silent or pose less pathogenicity. In this study, we use the NHLH2 gene which encodes the nescient helix-loop-helix 2 (Nhlh2) transcription factor in a systematic analysis of all missense mutations to date in the gene. The NHLH2 gene was first described in 1992. Knockout mice created in 1997 indicated a role for this protein in body weight control, puberty, and fertility, as well as the motivation for sex and exercise. Only recently have human carriers of NHLH2 missense variants been characterized. Over 300 missense variants for the NHLH2 gene are listed in the NCBI single nucleotide polymorphism database (dbSNP). Using in silico tools, predicted pathogenicity of the variants narrowed the missense variants to 37 which were predicted to affect NHLH2 function. These 37 variants cluster around the basic-helix-loop-helix and DNA binding domains of the transcription factor, and further analysis using in silico tools provided 21 SNV resulting in 22 amino acid changes for future wet lab analysis. The tools used, findings, and predictions for the variants are discussed considering the known function of the NHLH2 transcription factor. Overall use of these in silico tools and analysis of these data contribute to our knowledge of a protein which is both involved in the human genetic syndrome, Prader-Willi syndrome, and in controlling genes involved in body weight control, fertility, puberty, and behavior in the general population, and may provide a systematic methodology for others to characterize variants for their gene of interest.

Embeddings from protein language models predict conservation and variant effects

The emergence of SARS-CoV-2 variants stressed the demand for tools allowing to interpret the effect of single amino acid variants (SAVs) on protein function. While Deep Mutational Scanning (DMS) sets continue to expand our understanding of the mutational landscape of single proteins, the results continue to challenge analyses. Protein Language Models (pLMs) use the latest deep learning (DL) algorithms to leverage growing databases of protein sequences. These methods learn to predict missing or masked amino acids from the context of entire sequence regions. Here, we used pLM representations (embeddings) to predict sequence conservation and SAV effects without multiple sequence alignments (MSAs). Embeddings alone predicted residue conservation almost as accurately from single sequences as ConSeq using MSAs (two-state Matthews Correlation Coefficient-MCC-for ProtT5 embeddings of 0.596 ± 0.006 vs. 0.608 ± 0.006 for ConSeq). Inputting the conservation prediction along with BLOSUM62 substitution scores and pLM mask reconstruction probabilities into a simplistic logistic regression (LR) ensemble for Variant Effect Score Prediction without Alignments (VESPA) predicted SAV effect magnitude without any optimization on DMS data. Comparing predictions for a standard set of 39 DMS experiments to other methods (incl. ESM-1v, DeepSequence, and GEMME) revealed our approach as competitive with the state-of-the-art (SOTA) methods using MSA input. No method outperformed all others, neither consistently nor statistically significantly, independently of the performance measure applied (Spearman and Pearson correlation). Finally, we investigated binary effect predictions on DMS experiments for four human proteins. Overall, embedding-based methods have become competitive with methods relying on MSAs for SAV effect prediction at a fraction of the costs in computing/energy. Our method predicted SAV effects for the entire human proteome (~ 20 k proteins) within 40 min on one Nvidia Quadro RTX 8000. All methods and data sets are freely available for local and online execution through bioembeddings.com, https://github.com/Rostlab/VESPA , and PredictProtein.

Synthetic Biology: Bottom-Up Assembly of Molecular Systems

The bottom-up assembly of biological and chemical components opens exciting opportunities to engineer artificial vesicular systems for applications with previously unmet requirements. The modular combination of scaffolds and functional building blocks enables the engineering of complex systems with biomimetic or new-to-nature functionalities. Inspired by the compartmentalized organization of cells and organelles, lipid or polymer vesicles are widely used as model membrane systems to investigate the translocation of solutes and the transduction of signals by membrane proteins. The bottom-up assembly and functionalization of such artificial compartments enables full control over their composition and can thus provide specifically optimized environments for synthetic biological processes. This review aims to inspire future endeavors by providing a diverse toolbox of molecular modules, engineering methodologies, and different approaches to assemble artificial vesicular systems. Important technical and practical aspects are addressed and selected applications are presented, highlighting particular achievements and limitations of the bottom-up approach. Complementing the cutting-edge technological achievements, fundamental aspects are also discussed to cater to the inherently diverse background of the target audience, which results from the interdisciplinary nature of synthetic biology. The engineering of proteins as functional modules and the use of lipids and block copolymers as scaffold modules for the assembly of functionalized vesicular systems are explored in detail. Particular emphasis is placed on ensuring the controlled assembly of these components into increasingly complex vesicular systems. Finally, all descriptions are presented in the greater context of engineering valuable synthetic biological systems for applications in biocatalysis, biosensing, bioremediation, or targeted drug delivery.

Understanding the mutational frequency in SARS-CoV-2 proteome using structural features

The prolonged transmission of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus in the human population has led to demographic divergence and the emergence of several location-specific clusters of viral strains. Although the effect of mutation(s) on severity and survival of the virus is still unclear, it is evident that certain sites in the viral proteome are more/less prone to mutations. In fact, millions of SARS-CoV-2 sequences collected all over the world have provided us a unique opportunity to understand viral protein mutations and develop novel computational approaches to predict mutational patterns. In this study, we have classified the mutation sites into low and high mutability classes based on viral isolates count containing mutations. The physicochemical features and structural analysis of the SARS-CoV-2 proteins showed that features including residue type, surface accessibility, residue bulkiness, stability and sequence conservation at the mutation site were able to classify the low and high mutability sites. We further developed machine learning models using above-mentioned features, to predict low and high mutability sites at different selection thresholds (ranging 5-30% of topmost and bottommost mutated sites) and observed the improvement in performance as the selection threshold is reduced (prediction accuracy ranging from 65 to 77%). The analysis will be useful for early detection of variants of concern for the SARS-CoV-2, which can also be applied to other existing and emerging viruses for another pandemic prevention.

Integrative assessment of CIP2A overexpression and mutational effects in human malignancies identifies possible deleterious variants

KIAA1524 is the gene encoding the human cancerous inhibitor of PP2A (CIP2A) protein which is regarded as a novel target for cancer therapy. It is overexpressed in 65%-90% of tissues in almost all studied human cancers. CIP2A expression correlates with cancer progression, disease aggressivity in lung cancer besides poor survival and resistance to chemotherapy in breast cancer. Herein, a pan-cancer analysis of public gene expression datasets was conducted showing significant upregulation of CIP2A in cancerous and metastatic tissues. CIP2A overexpression also correlated with poor survival of cancer patients. To determine the non-coding variants associated with CIP2A overexpression, 5'UTR and 3'UTR variants were annotated and scored using RegulomeDB and Enformer deep learning model. The 5'UTR variants rs1239349555, rs1576326380, and rs1231839144 were predicted to be potential regulators of CIP2A overexpression scoring best on RegulomeDB annotations with a high "2a" rank of supporting experimental data. These variants also scored the highest on Enformer predictions. Analysis of the 3'UTR variants of CIP2A predicted rs56255137 and rs58758610 to alter binding sites of hsa-miR-500a-5 and (hsa-miR-3671, hsa-miR-5692a) respectively. Both variants were also found in linkage disequilibrium with rs11709183 and rs147863209 respectively at r² ≥ 0.8. The aforementioned variants were found to be eQTL hits significantly associated with CIP2A overexpression. Further, analysis of rs11709183 and rs147863209 revealed a high "2b" rank on RegulomeDB annotations indicating a probable effect on DNAse transcription factors binding. The MuTarget analysis indicated that somatic mutations in TP53 are significantly associated with upregulated CIP2A in human cancers. Analysis of missense SNPs on CIP2A solved structure predicted seven deleterious effects. Four of these variants were also predicted as structurally and functionally destabilizing to CIP2A including; rs375108755, rs147942716, rs368722879, and rs367941403. Variant rs1193091427 was predicted as a potential intronic splicing mutation that might be responsible for the novel CIP2A variant (NOCIVA) in multiple myeloma. Finally, Enrichment of the Wnt/β-catenin pathway within the CIP2A regulatory gene network suggested potential of therapeutic combinations between FTY720 with Wnt/β-catenin, Plk1 and/or HDAC inhibitors to downregulate CIP2A which has been shown to be essential for the survival of different cancer cell lines.

Protein engineering for natural product biosynthesis and synthetic biology applications

As protein engineering grows more salient, many strategies have emerged to alter protein structure and function, with the goal of redesigning and optimizing natural product biosynthesis. Computational tools, including machine learning and molecular dynamics simulations, have enabled the rational mutagenesis of key catalytic residues for enhanced or altered biocatalysis. Semi-rational, directed evolution and microenvironment engineering strategies have optimized catalysis for native substrates and increased enzyme promiscuity beyond the scope of traditional rational approaches. These advances are made possible using novel high-throughput screens, including designer protein-based biosensors with engineered ligand specificity. Herein, we detail the most recent of these advances, focusing on polyketides, non-ribosomal peptides and isoprenoids, including their native biosynthetic logic to provide clarity for future applications of these technologies for natural product synthetic biology.

Exploring mutable conserved sites and fatal non-conserved sites by random mutation of esterase from Sulfolobus tokodaii and subtilisin from Thermococcus kodakarensis

Homologous proteins differ in their amino acid sequences at several positions. Generally, conserved sites are recognized as not suitable for amino acid substitution, and thus in evolutionary protein engineering, non-conserved sites are often selected as mutation sites. However, there have also been reports of possible mutations in conserved sites. In this study, we explored mutable conserved sites and immutable non-conserved sites by testing random mutations of two thermostable proteins, an esterase from Sulfolobus tokodaii (Sto-Est) and a subtilisin from Thermococcus kodakarensis (Tko-Sub). The subtilisin domain of Tko-Sub needs Ca²⁺ ions and the propeptide domain for stability, folding and maturation. The results from the two proteins showed that about one-third of the mutable sites were detected in conserved sites and some non-conserved sites lost enzymatic activity at high temperatures due to mutation. Of the conserved sites in Sto-Est, the sites on the loop, on the surface, and far from the active site are more resistant to mutation. In Tko-Sub, the sites flanking Ca²⁺-binding sites and propeptide were undesirable for mutation. The results presented here serve as an index for selecting mutation sites and contribute to the expansion of available sequence range by introducing mutations at conserved sites.

When a foreign gene meets its native counterpart: computational biophysics analysis of two PgiC loci in the grass Festuca ovina

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the “native” PgiC1 and “foreign” PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

Variant effect predictions capture some aspects of deep mutational scanning experiments

BACKGROUND

Deep mutational scanning (DMS) studies exploit the mutational landscape of sequence variation by systematically and comprehensively assaying the effect of single amino acid variants (SAVs; also referred to as missense mutations, or non-synonymous Single Nucleotide Variants - missense SNVs or nsSNVs) for particular proteins. We assembled SAV annotations from 22 different DMS experiments and normalized the effect scores to evaluate variant effect prediction methods. Three trained on traditional variant effect data (PolyPhen-2, SIFT, SNAP2), a regression method optimized on DMS data (Envision), and a naïve prediction using conservation information from homologs.

RESULTS

On a set of 32,981 SAVs, all methods captured some aspects of the experimental effect scores, albeit not the same. Traditional methods such as SNAP2 correlated slightly more with measurements and better classified binary states (effect or neutral). Envision appeared to better estimate the precise degree of effect. Most surprising was that the simple naïve conservation approach using PSI-BLAST in many cases outperformed other methods. All methods captured beneficial effects (gain-of-function) significantly worse than deleterious (loss-of-function). For the few proteins with multiple independent experimental measurements, experiments differed substantially, but agreed more with each other than with predictions.

CONCLUSIONS

DMS provides a new powerful experimental means of understanding the dynamics of the protein sequence space. As always, promising new beginnings have to overcome challenges. While our results demonstrated that DMS will be crucial to improve variant effect prediction methods, data diversity hindered simplification and generalization.

Computer-Aided Saturation Mutagenesis of Arabidopsis thaliana Ent-Copalyl Diphosphate Synthase

Ent-copalyl diphosphate synthase controls the biosynthesis of gibberellin plant hormones, which in turn coordinate the expression of numerous enzymes. Some gibberellin-dependent genes encode enzymes coordinating the biosynthesis of tanshinones: diterpene derivatives with broad medical applications. New biotechnological approaches, such as metabolic engineering using naturally occurring or mutated enzymes, have been proposed to meet the growing demand for tanshinones which is currently met by the Chinese medicinal plant Salvia miltiorrhiza Bunge. These mutants may be prepared by directed evolution, saturation mutagenesis or rational enzyme design. In the presented paper, 15,257 non-synonymous variants of Arabidopsis thaliana ent-copalyl diphosphate synthase were obtained using the SNAP2 tool. The obtained forms were screened to isolate variants with potentially improved biological functions. A group of 455 mutants with potentially improved stability was isolated and subjected to further screening on the basis of ligand–substrate affinity, and both secondary structure and active site structure stability. Finally, a group of six single mutants was obtained, which were used to construct double mutants with potentially improved stability and ligand affinity. The potential influence of single mutations on protein stability and ligand affinity was evaluated by double mutant cycle analysis. Finally, the procedure was validated by in silico assessment of the experimentally verified enzyme mutants with reduced enzymatic activity.

Human collectin-11 (COLEC11) and its synergic genetic interaction with MASP2 are associated with the pathophysiology of Chagas Disease

Chagas Disease (CD) is an anthropozoonosis caused by Trypanosoma cruzi. With complex pathophysiology and variable clinical presentation, CD outcome can be influenced by parasite persistence and the host immune response. Complement activation is one of the primary defense mechanisms against pathogens, which can be initiated via pathogen recognition by pattern recognition molecules (PRMs). Collectin-11 is a multifunctional soluble PRM lectin, widely distributed throughout the body, with important participation in host defense, homeostasis, and embryogenesis. In complex with mannose-binding lectin-associated serine proteases (MASPs), collectin-11 may initiate the activation of complement, playing a role against pathogens, including T. cruzi. In this study, collectin-11 plasma levels and COLEC11 variants in exon 7 were assessed in a Brazilian cohort of 251 patients with chronic CD and 108 healthy controls. Gene-gene interactions between COLEC11 and MASP2 variants were analyzed. Collectin-11 levels were significantly decreased in CD patients compared to controls (p<0.0001). The allele rs7567833G, the genotypes rs7567833AG and rs7567833GG, and the COLEC11*GGC haplotype were related to T. cruzi infection and clinical progression towards symptomatic CD. COLEC11 and MASP2*CD risk genotypes were associated with cardiomyopathy (p = 0.014; OR 9.3, 95% CI 1.2–74) and with the cardiodigestive form of CD (p = 0.005; OR 15.2, 95% CI 1.7–137), suggesting that both loci act synergistically in immune modulation of the disease. The decreased levels of collectin-11 in CD patients may be associated with the disease process. The COLEC11 variant rs7567833G and also the COLEC11 and MASP2*CD risk genotype interaction were associated with the pathophysiology of CD.

The heterogeneity of clinical progression during chronic Trypanosoma cruzi infection and the mechanisms determining why some individuals develop symptoms whereas others remain asymptomatic are still poorly understood. The pathogenesis of chronic Chagas Disease (CD) has been attributed mainly to the persistence of the causing parasite and the character of individual host immune responses. Collectin-11 is a host immune response molecule with affinity for sugars found on the T. cruzi’s surface. Together with mannose-binding lectin-associated serine proteases (MASPs), it triggers the host defense response against pathogens. Genetic variants and protein levels of MASP-2 and the mannose-binding lectin (MBL), a molecule structurally similar to collectin-11, have been found to be associated with susceptibility to T. cruzi infection and clinical progression to cardiomyopathy. This prompted us to investigate collectin-11 genetic variants and protein levels in 251 patients with chronic CD and 108 healthy individuals, and to examine the effect of gene interaction between COLEC11 and MASP2 risk mutations. We found an association to CD infection with COLEC11 gene variants and reduced collectin-11 levels. The concomitant presence of these genetic variants and MASP2 risk mutations greatly increased the odds for cardiomyopathy. This is the first study to reveal a role for collectin-11 and COLEC11-MASP2 gene interaction in the pathogenesis of CD.

Large-scale in-silico statistical mutagenesis analysis sheds light on the deleteriousness landscape of the human proteome

Next generation sequencing technologies are providing increasing amounts of sequencing data, paving the way for improvements in clinical genetics and precision medicine. The interpretation of the observed genomic variants in the light of their phenotypic effects is thus emerging as a crucial task to solve in order to advance our understanding of how exomic variants affect proteins and how the proteins' functional changes affect human health. Since the experimental evaluation of the effects of every observed variant is unfeasible, Bioinformatics methods are being developed to address this challenge in-silico, by predicting the impact of millions of variants, thus providing insight into the deleteriousness landscape of entire proteomes. Here we show the feasibility of this approach by using the recently developed DEOGEN2 variant-effect predictor to perform the largest in-silico mutagenesis scan to date. We computed the deleteriousness score of 170 million variants over 15000 human proteins and we analysed the results, investigating how the predicted deleteriousness landscape of the proteins relates to known functionally and structurally relevant protein regions and biophysical properties. Moreover, we qualitatively validated our results by comparing them with two mutagenesis studies targeting two specific proteins, showing the consistency of DEOGEN2 predictions with respect to experimental data.

In Silico Analysis of Missense Mutations as a First Step in Functional Studies: Examples from Two Sphingolipidoses

In order to delineate a better approach to functional studies, we have selected 23 missense mutations distributed in different domains of two lysosomal enzymes, to be studied by in silico analysis. In silico analysis of mutations relies on computational modeling to predict their effects. Various computational platforms are currently available to check the probable causality of mutations encountered in patients at the protein and at the RNA levels. In this work we used four different platforms freely available online (Protein Variation Effect Analyzer- PROVEAN, PolyPhen-2, Swiss-model Expert Protein Analysis System—ExPASy, and SNAP2) to check amino acid substitutions and their effect at the protein level. The existence of functional studies, regarding the amino acid substitutions, led to the selection of the distinct protein mutants. Functional data were used to compare the results obtained with different bioinformatics tools. With the advent of next-generation sequencing, it is not feasible to carry out functional tests in all the variants detected. In silico analysis seems to be useful for the delineation of which mutants are worth studying through functional studies. Therefore, prediction of the mutation impact at the protein level, applying computational analysis, confers the means to rapidly provide a prognosis value to genotyping results, making it potentially valuable for patient care as well as research purposes. The present work points to the need to carry out functional studies in mutations that might look neutral. Moreover, it should be noted that single nucleotide polymorphisms (SNPs), occurring in coding and non-coding regions, may lead to RNA alterations and should be systematically verified. Functional studies can gain from a preliminary multi-step approach, such as the one proposed here.

Characterizing the relation of functional and Early Folding Residues in protein structures using the example of aminoacyl-tRNA synthetases

Proteins are chains of amino acids which adopt a three-dimensional structure and are then able to catalyze chemical reactions or propagate signals in organisms. Without external influence, many proteins fold into their native structure, and a small number of Early Folding Residues (EFR) have previously been shown to initiate the formation of secondary structure elements and guide their respective assembly. Using the two diverse superfamilies of aminoacyl-tRNA synthetases (aaRS), it is shown that the position of EFR is preserved over the course of evolution even when the corresponding sequence conservation is small. Folding initiation sites are positioned in the center of secondary structure elements, independent of aaRS class. In class I, the predicted position of EFR resembles an ancient structural packing motif present in many seemingly unrelated proteins. Furthermore, it is shown that EFR and functionally relevant residues in aaRS are almost entirely disjoint sets of residues. The Start2Fold database is used to investigate whether this separation of EFR and functional residues can be observed for other proteins. EFR are found to constitute crucial connectors of protein regions which are distant at sequence level. Especially, these residues exhibit a high number of non-covalent residue-residue contacts such as hydrogen bonds and hydrophobic interactions. This tendency also manifests as energetically stable local regions, as substantiated by a knowledge-based potential. Despite profound differences regarding how EFR and functional residues are embedded in protein structures, a strict separation of structurally and functionally relevant residues cannot be observed for a more general collection of proteins.

Sequence variant analysis of RNA sequences in severe equine asthma

Background

Severe equine asthma is a chronic inflammatory disease of the lung in horses similar to low-Th2 late-onset asthma in humans. This study aimed to determine the utility of RNA-Seq to call gene sequence variants, and to identify sequence variants of potential relevance to the pathogenesis of asthma.

Methods

RNA-Seq data were generated from endobronchial biopsies collected from six asthmatic and seven non-asthmatic horses before and after challenge (26 samples total). Sequences were aligned to the equine genome with Spliced Transcripts Alignment to Reference software. Read preparation for sequence variant calling was performed with Picard tools and Genome Analysis Toolkit (GATK). Sequence variants were called and filtered using GATK and Ensembl Variant Effect Predictor (VEP) tools, and two RNA-Seq predicted sequence variants were investigated with both PCR and Sanger sequencing. Supplementary analysis of novel sequence variant selection with VEP was based on a score of <0.01 predicted with Sorting Intolerant from Tolerant software, missense nature, location within the protein coding sequence and presence in all asthmatic individuals. For select variants, effect on protein function was assessed with Polymorphism Phenotyping 2 and screening for non-acceptable polymorphism 2 software. Sequences were aligned and 3D protein structures predicted with Geneious software. Difference in allele frequency between the groups was assessed using a Pearson’s Chi-squared test with Yates’ continuity correction, and difference in genotype frequency was calculated using the Fisher’s exact test for count data.

Results

RNA-Seq variant calling and filtering correctly identified substitution variants in PACRG and RTTN. Sanger sequencing confirmed that the PACRG substitution was appropriately identified in all 26 samples while the RTTN substitution was identified correctly in 24 of 26 samples. These variants of uncertain significance had substitutions that were predicted to result in loss of function and to be non-neutral. Amino acid substitutions projected no change of hydrophobicity and isoelectric point in PACRG, and a change in both for RTTN. For PACRG, no difference in allele frequency between the two groups was detected but a higher proportion of asthmatic horses had the altered RTTN allele compared to non-asthmatic animals.

Discussion

RNA-Seq was sensitive and specific for calling gene sequence variants in this disease model. Even moderate coverage (<10–20 counts per million) yielded correct identification in 92% of samples, suggesting RNA-Seq may be suitable to detect sequence variants in low coverage samples. The impact of amino acid alterations in PACRG and RTTN proteins, and possible association of the sequence variants with asthma, is of uncertain significance, but their role in ciliary function may be of future interest.

Rapid and Error-Free Site-Directed Mutagenesis by a PCR-Free In Vitro CRISPR/Cas9-Mediated Mutagenic System

The quality and efficiency of any PCR-based mutagenesis technique may not be optimal due to, among other things, amino acid bias, which means that the development of efficient PCR-free methods is desirable. Here, we present a highly efficient in vitro CRISPR/Cas9-mediated mutagenic (ICM) system that allows rapid construction of designed mutants in a PCR-free manner. First, it involves plasmid digestion by utilizing a complex of Cas9 with specific single guide RNA (sgRNA) followed by degradation with T5 exonuclease to generate a 15 nt homologous region. Second, primers containing the desired mutations are annealed to form the double-stranded DNA fragments, which are then ligated into the linearized plasmid. In theory, neither the size of the target plasmid nor the unavailable restriction enzyme site poses any problems that may arise in traditional techniques. In this study, single and multiple site-directed mutagenesis were successfully performed even for a large size plasmid (up to 9.0 kb). Moreover, a PCR-free site-saturation mutagenesis library on single site and two adjacent sites of a green fluorescent protein was also generated with promising results. This demonstrates the great potential of the ICM system for creating high-quality mutant libraries in directed evolution as an alternative to PCR-based saturation mutagenesis, thus facilitating research on synthetic biology.

Boosting the efficiency of site-saturation mutagenesis for a difficult-to-randomize gene by a two-step PCR strategy

Site-saturation mutagenesis (SSM) has been used in directed evolution of proteins for a long time. As a special form of saturation mutagenesis, it involves individual randomization at a given residue with formation of all 19 amino acids. To date, the most efficient embodiment of SSM is a one-step PCR-based approach using NNK codon degeneracy. However, in the case of difficult-to-randomize genes, SSM may not deliver all of the expected 19 mutants, which compels the user to invest further efforts by applying site-directed mutagenesis for the construction of the missing mutants. To solve this problem, we developed a two-step PCR-based technique in which a mutagenic primer and a non-mutagenic (silent) primer are used to generate a short DNA fragment, which is recovered and then employed as a megaprimer to amplify the whole plasmid. The present two-step and older one-step (partially overlapped primer approach) procedures were compared by utilizing cytochrome P450-BM3, which is a "difficult-to-randomize" gene. The results document the distinct superiority of the new method by checking the library quality on DNA level based on massive sequence data, but also at amino acid level. Various future applications in biotechnology can be expected, including the utilization when constructing mutability landscapes, which provide semi-rational information for identifying hot spots for protein engineering and directed evolution.

PhosphoPICK-SNP: quantifying the effect of amino acid variants on protein phosphorylation

Motivation

Genome-wide association studies are identifying single nucleotide variants (SNVs) linked to various diseases, however the functional effect caused by these variants is often unknown. One potential functional effect, the loss or gain of protein phosphorylation sites, can be induced through variations in key amino acids that disrupt or introduce valid kinase binding patterns. Current methods for predicting the effect of SNVs on phosphorylation operate on the sequence content of reference and variant proteins. However, consideration of the amino acid sequence alone is insufficient for predicting phosphorylation change, as context factors determine kinase-substrate selection.

Results

We present here a method for quantifying the effect of SNVs on protein phosphorylation through an integrated system of motif analysis and context-based assessment of kinase targets. By predicting the effect that known variants across the proteome have on phosphorylation, we are able to use this background of proteome-wide variant effects to quantify the significance of novel variants for modifying phosphorylation. We validate our method on a manually curated set of phosphorylation change-causing variants from the primary literature, showing that the method predicts known examples of phosphorylation change at high levels of specificity. We apply our approach to data-sets of variants in phosphorylation site regions, showing that variants causing predicted phosphorylation loss are over-represented among disease-associated variants.

Availability and Implementation

The method is freely available as a web-service at the website http://bioinf.scmb.uq.edu.au/phosphopick/snp.

Contact

m.boden@uq.edu.au.

Supplementary information

Supplementary data are available at Bioinformatics online.

Computational Analysis of High-Risk SNPs in Human DBY Gene Responsible for Male Infertility: A Functional and Structural Impact

BACKGROUND

DEAD-box helicase 3, Y-linked (DBY) is a candidate gene of the AZF region which is involved in spermatogenesis process. Mutations in the DBY gene may disrupt the spermatogenesis and lead to infertility in men. Identification of functionally neutral mutation from the disease-causing mutation is the biggest challenge in human genetic variation analysis. Owing to the importance of DBY in male infertility, functional analysis was carried out to reveal the association between genetic mutation and phenotypic variation through various in silico approaches.

METHODS

The present study analyzed the functional consequences of the nsSNPs in human DBY gene using SIFT, PolyPhen 2, PROVEAN, SNAP2, PMut, nsSNPAnalyzer, PhD-SNP and SNPs&GO along with stability analysis through I-Mutant2.0, MuPro and iPTREE-STAB. The conservational analysis of amino acid residues, biophysical properties and conserved domains of the DBY protein was analyzed using various computational tools. The 3D structure of the protein was generated using SPARKS-X and validated using RAMPAGE.

RESULTS

Out of 1130 SNPs reported in dbSNP, only one nsSNP (G300D) was found to have a functional effect on stability as well as the function of the DBY protein. The results showed the presence of G300 in the putative structure of DBY domain.

CONCLUSION

To the best of our knowledge, this is the first study to detect pathologically significant nsSNPs (G300D) through a computational approach in the DBY which can be useful for development in potent drug discovery studies.

Comprehensive Analysis of Non-Synonymous Natural Variants of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are the largest superfamily of transmembrane receptors and have vital signaling functions in various organs. Because of their critical roles in physiology and pathology, GPCRs are the most commonly used therapeutic target. It has been suggested that GPCRs undergo massive genetic variations such as genetic polymorphisms and DNA insertions or deletions. Among these genetic variations, non-synonymous natural variations change the amino acid sequence and could thus alter GPCR functions such as expression, localization, signaling, and ligand binding, which may be involved in disease development and altered responses to GPCR-targeting drugs. Despite the clinical importance of GPCRs, studies on the genotype-phenotype relationship of GPCR natural variants have been limited to a few GPCRs such as β-adrenergic receptors and opioid receptors. Comprehensive understanding of non-synonymous natural variations within GPCRs would help to predict the unknown genotype-phenotype relationship and yet-to-be-discovered natural variants. Here, we analyzed the non-synonymous natural variants of all non-olfactory GPCRs available from a public database, UniProt. The results suggest that non-synonymous natural variations occur extensively within the GPCR superfamily especially in the N-terminus and transmembrane domains. Within the transmembrane domains, natural variations observed more frequently in the conserved residues, which leads to disruption of the receptor function. Our analysis also suggests that only few non-synonymous natural variations have been studied in efforts to link the variations with functional consequences.

Directed Evolution of Proteins Based on Mutational Scanning

Directed evolution has emerged as one of the most effective protein engineering methods in basic research as well as in applications in synthetic organic chemistry and biotechnology. The successful engineering of protein activity, allostery, binding affinity, expression, folding, fluorescence, solubility, substrate scope, selectivity (enantio-, stereo-, and regioselectivity), and/or stability (temperature, organic solvents, pH) is just limited by the throughput of the genetic selection, display, or screening system that is available for a given protein. Sometimes it is possible to analyze millions of protein variants from combinatorial libraries per day. In other cases, however, only a few hundred variants can be screened in a single day, and thus the creation of smaller yet smarter libraries is needed. Different strategies have been developed to create these libraries. One approach is to perform mutational scanning or to construct "mutability landscapes" in order to understand sequence-function relationships that can guide the actual directed evolution process. Herein we provide a protocol for economically constructing scanning mutagenesis libraries using a cytochrome P450 enzyme in a high-throughput manner. The goal is to engineer activity, regioselectivity, and stereoselectivity in the oxidative hydroxylation of a steroid, a challenging reaction in synthetic organic chemistry. Libraries based on mutability landscapes can be used to engineer any fitness trait of interest. The protocol is also useful for constructing gene libraries for deep mutational scanning experiments.

Novel non-neutral mitochondrial DNA mutations found in childhood acute lymphoblastic leukemia

Mitochondria produce adenosine triphosphate (ATP) for energy requirements via the mitochondrial oxidative phosphorylation (OXPHOS) system. One of the hallmarks of cancer is the energy shift toward glycolysis. Low OXPHOS activity and increased glycolysis are associated with aggressive types of cancer. Mitochondria have their own genome (mitochondrial DNA [mtDNA]) encoding for 13 essential subunits of the OXPHOS enzyme complexes. We studied mtDNA in childhood acute lymphoblastic leukemia (ALL) to detect potential pathogenic mutations in OXPHOS complexes. The whole mtDNA from blood and bone marrow samples at diagnosis and follow-up from 36 ALL patients were analyzed. Novel or previously described pathogenic mtDNA mutations were identified in 8 out of 36 patients. Six out of these 8 patients had died from ALL. Five out of 36 patients had an identified poor prognosis genetic marker, and 4 of these patients had mtDNA mutations. Missense or nonsense mtDNA mutations were detected in the genes encoding subunits of OXPHOS complexes, as follows: MT-ND1, MT-ND2, MT-ND4L and MT-ND6 of complex I; MT-CO3 of complex IV; and MT-ATP6 and MT-ATP8 of complex V. We discovered mtDNA mutations in childhood ALL supporting the hypothesis that non-neutral variants in mtDNA affecting the OXPHOS function may be related to leukemic clones.

DYRK1B mutations associated with metabolic syndrome impair the chaperone-dependent maturation of the kinase domain

Two missense mutations of the DYRK1B gene have recently been found to co-segregate with a rare autosomal-dominant form of metabolic syndrome. This gene encodes a member of the DYRK family of protein kinases, which depend on tyrosine autophosphorylation to acquire the catalytically active conformation. The mutations (H90P and R102C) affect a structural element named DYRK homology (DH) box and did not directly interfere with the conformation of the catalytic domain in a structural model of DYRK1B. Cellular assays showed that the mutations did not alter the specific activity of mature kinase molecules. However, a significant part of the mutant DYRK1B protein accumulated in detergent-insoluble cytoplasmic aggregates and was underphosphorylated on tyrosine. The mutant DYRK1B variants were more vulnerable to the HSP90 inhibitor ganetespib and showed enhanced binding to the co-chaperone CDC37 as compared to wild type DYRK1B. These results support the hypothesis that the mutations in the DH box interfere with the maturation of DYRK1B by tyrosine autophosphorylation and compromise the conformational stability of the catalytic domain, which renders the kinase susceptible to misfolding.

Computational analysis for the determination of deleterious nsSNPs in human MTHFD1 gene

Single nucleotide polymorphisms (SNPs) are the most common genetic polymorphisms and play a major role in many inherited diseases. Methylenetetrahydrofolate dehydrogenase 1 (MTHFD1) is one of the enzymes involved in folate metabolism. In the present study, the functional and structural consequences of nsSNPs of human MTHFD1 gene was analyzed using various computational tools like SIFT, PolyPhen2, PANTHER, PROVEAN, SNAP2, nsSNPAnalyzer, PhD-SNP, SNPs&GO, I-Mutant, MuPro, ConSurf, InterPro, NCBI Conserved Domain Search tool, ModPred, SPARKS-X, RAMPAGE, FT Site and PyMol. Out of 327 nsSNPs form human MTHFD1 gene, total 45 SNPs were predicted as functionally most significant SNPs, among which 17 were highly conserved and functional, 17 were highly conserved and structural residues. Among 45 most significant SNPs, 15 were predicted to be involved in post translational modifications. The p.Gly165Arg may interfere in homodimer interface formation. The p.Asn439Lys and p.Asp445Asn may interfere in binding interactions of MTHFD1 protein with cesium cation and potassium. The two SNPs (p.Asp562Gly and p.Gly637Cys) might interfere in interactions of MTHFD1 with ligand.

The Generation and Exploitation of Protein Mutability Landscapes for Enzyme Engineering

The increasing number of enzyme applications in chemical synthesis calls for new engineering methods to develop the biocatalysts of the future. An interesting concept in enzyme engineering is the generation of large-scale mutational data in order to chart protein mutability landscapes. These landscapes allow the important discrimination between beneficial mutations and those that are neutral or detrimental, thus providing detailed insight into sequence-function relationships. As such, mutability landscapes are a powerful tool with which to identify functional hotspots at any place in the amino acid sequence of an enzyme. These hotspots can be used as targets for combinatorial mutagenesis to yield superior enzymes with improved catalytic properties, stability, or even new enzymatic activities. The generation of mutability landscapes for multiple properties of one enzyme provides the exciting opportunity to select mutations that are beneficial either for one or for several of these properties. This review presents an overview of the recent advances in the construction of mutability landscapes and discusses their importance for enzyme engineering.

Predicted Molecular Effects of Sequence Variants Link to System Level of Disease

Developments in experimental and computational biology are advancing our understanding of how protein sequence variation impacts molecular protein function. However, the leap from the micro level of molecular function to the macro level of the whole organism, e.g. disease, remains barred. Here, we present new results emphasizing earlier work that suggested some links from molecular function to disease. We focused on non-synonymous single nucleotide variants, also referred to as single amino acid variants (SAVs). Building upon OMIA (Online Mendelian Inheritance in Animals), we introduced a curated set of 117 disease-causing SAVs in animals. Methods optimized to capture effects upon molecular function often correctly predict human (OMIM) and animal (OMIA) Mendelian disease-causing variants. We also predicted effects of human disease-causing variants in the mouse model, i.e. we put OMIM SAVs into mouse orthologs. Overall, fewer variants were predicted with effect in the model organism than in the original organism. Our results, along with other recent studies, demonstrate that predictions of molecular effects capture some important aspects of disease. Thus, in silico methods focusing on the micro level of molecular function can help to understand the macro system level of disease.

The variations in the genetic sequence between individuals affect the gene-product, i.e. the protein differently. Some variants have no measurable effect (are neutral), while others affect protein function. Some of those effects are so severe they cause so called monogenic Mendelian diseases, i.e. diseases triggered by a single letter change. Some in silico methods predict the molecular impact of sequence variation. However, both experimental and computational analyses struggle to generalize from the effect upon molecular protein function to the effect upon the organism such as a disease. Here, we confirmed that methods predicting molecular effects correctly capture the type of effects causing Mendelian diseases in human and introduced a data set for animal diseases that was also captured by predictions methods. Predicted effects were less when in silico testing human variants in an animal model (here mouse). This is important to know because “mouse models” are common to study human diseases. Overall, we provided some evidence for a link between the molecular level and some type of disease.

Bioinformatics Approach for Prediction of Functional Coding/Noncoding Simple Polymorphisms (SNPs/Indels) in Human BRAF Gene

This study was carried out for Homo sapiens single variation (SNPs/Indels) in BRAF gene through coding/non-coding regions. Variants data was obtained from database of SNP even last update of November, 2015. Many bioinformatics tools were used to identify functional SNPs and indels in proteins functions, structures and expressions. Results shown, for coding polymorphisms, 111 SNPs predicted as highly damaging and six other were less. For UTRs, showed five SNPs and one indel were altered in micro RNAs binding sites (3' UTR), furthermore nil SNP or indel have functional altered in transcription factor binding sites (5' UTR). In addition for 5'/3' splice sites, analysis showed that one SNP within 5' splice site and one Indel in 3' splice site showed potential alteration of splicing. In conclude these previous functional identified SNPs and indels could lead to gene alteration, which may be directly or indirectly contribute to the occurrence of many diseases.

Detection and sequence/structure mapping of biophysical constraints to protein variation in saturated mutational libraries and protein sequence alignments with a dedicated server

Background

Protein variability can now be studied by measuring high-resolution tolerance-to-substitution maps and fitness landscapes in saturated mutational libraries. But these rich and expensive datasets are typically interpreted coarsely, restricting detailed analyses to positions of extremely high or low variability or dubbed important beforehand based on existing knowledge about active sites, interaction surfaces, (de)stabilizing mutations, etc.

Results

Our new webserver PsychoProt (freely available without registration at http://psychoprot.epfl.ch or at http://lucianoabriata.altervista.org/psychoprot/index.html) helps to detect, quantify, and sequence/structure map the biophysical and biochemical traits that shape amino acid preferences throughout a protein as determined by deep-sequencing of saturated mutational libraries or from large alignments of naturally occurring variants.

Discussion

We exemplify how PsychoProt helps to (i) unveil protein structure-function relationships from experiments and from alignments that are consistent with structures according to coevolution analysis, (ii) recall global information about structural and functional features and identify hitherto unknown constraints to variation in alignments, and (iii) point at different sources of variation among related experimental datasets or between experimental and alignment-based data. Remarkably, metabolic costs of the amino acids pose strong constraints to variability at protein surfaces in nature but not in the laboratory. This and other differences call for caution when extrapolating results from in vitro experiments to natural scenarios in, for example, studies of protein evolution.

Conclusion

We show through examples how PsychoProt can be a useful tool for the broad communities of structural biology and molecular evolution, particularly for studies about protein modeling, evolution and design.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1124-4) contains supplementary material, which is available to authorized users.

Transthyretin variants with improved inhibition of β-amyloid aggregation

Aggregation of β-amyloid (Aβ) is widely believed to cause neuronal dysfunction in Alzheimer's disease. Transthyretin (TTR) binds to Aβ and inhibits its aggregation and neurotoxicity. TTR is a homotetrameric protein, with each monomer containing a short α-helix and two anti-parallel β-sheets. Dimers pack into tetramers to form a hydrophobic cavity. Here we report the discovery of a TTR mutant, N98A, that was more effective at inhibiting Aβ aggregation than wild-type (WT) TTR, although N98A and WT bound Aβ equally. The N98A mutation is located on a flexible loop distant from the putative Aβ-binding sites and does not alter secondary and tertiary structures nor prevent correct assembly into tetramers. Under non-physiological conditions, N98A tetramers were kinetically and thermodynamically less stable than WT, suggesting a difference in the tetramer folded structure. In vivo, the lone cysteine in TTR is frequently modified by S-cysteinylation or S-sulfonation. Like the N98A mutation, S-cysteinylation of TTR modestly decreased tetramer stability and increased TTR's effectiveness at inhibiting Aβ aggregation. Collectively, these data indicate that a subtle change in TTR tetramer structure measurably increases TTR's ability to inhibit Aβ aggregation.

Maize peroxidase Px5 has a highly conserved sequence in inbreds resistant to mycotoxin producing fungi which enhances fungal and insect resistance

Mycotoxin presence in maize causes health and economic issues for humans and animals. Although many studies have investigated expression differences of genes putatively governing resistance to producing fungi, few have confirmed a resistance role, or examined putative resistance gene structure in more than a couple of inbreds. The pericarp expression of maize Px5 has previously been associated with resistance to Aspergillus flavus growth and insects in a set of inbreds. Genes from 14 different inbreds that included ones with resistance and susceptibility to A. flavus, Fusarium proliferatum, F. verticillioides and F. graminearum and/or mycotoxin production were cloned using high fidelity enzymes, and sequenced. The sequence of Px5 from all resistant inbreds was identical, except for a single base change in two inbreds, only one of which affected the amino acid sequence. Conversely, the Px5 sequence from several susceptible inbreds had several base variations, some of which affected amino acid sequence that would potentially alter secondary structure, and thus enzyme function. The sequence of the maize peroxidase Px5 common to inbreds resistant to mycotoxigenic fungi was overexpressed in maize callus. Callus transformants overexpressing the gene caused significant reductions in growth for fall armyworms, corn earworms, and F. graminearum compared to transformant callus with a β-glucuronidase gene. This study demonstrates rarer transcripts of potential resistance genes overlooked by expression screens can be identified by sequence comparisons. A role in pest resistance can be verified by callus expression of the candidate genes, which can thereby justify larger scale transformation and regeneration of transgenic plants expressing the resistance gene for further evaluation.

Antibody Binding Selectivity: Alternative Sets of Antigen Residues Entail High-Affinity Recognition

Understanding the relationship between protein sequence and molecular recognition selectivity remains a major challenge. The antibody fragment scFv1F4 recognizes with sub nM affinity a decapeptide (sequence 6TAMFQDPQER15) derived from the N-terminal end of human papilloma virus E6 oncoprotein. Using this decapeptide as antigen, we had previously shown that only the wild type amino-acid or conservative replacements were allowed at positions 9 to 12 and 15 of the peptide, indicating a strong binding selectivity. Nevertheless phenylalanine (F) was equally well tolerated as the wild type glutamine (Q) at position 13, while all other amino acids led to weaker scFv binding. The interfaces of complexes involving either Q or F are expected to diverge, due to the different physico-chemistry of these residues. This would imply that high-affinity binding can be achieved through distinct interfacial geometries. In order to investigate this point, we disrupted the scFv-peptide interface by modifying one or several peptide positions. We then analyzed the effect on binding of amino acid changes at the remaining positions, an altered susceptibility being indicative of an altered role in complex formation. The 23 starting variants analyzed contained replacements whose effects on scFv1F4 binding ranged from minor to drastic. A permutation analysis (effect of replacing each peptide position by all other amino acids except cysteine) was carried out on the 23 variants using the PEPperCHIP® Platform technology. A comparison of their permutation patterns with that of the wild type peptide indicated that starting replacements at position 11, 12 or 13 modified the tolerance to amino-acid changes at the other two positions. The interdependence between the three positions was confirmed by SPR (Biacore® technology). Our data demonstrate that binding selectivity does not preclude the existence of alternative high-affinity recognition modes.

Systematic Mapping of Protein Mutational Space by Prolonged Drift Reveals the Deleterious Effects of Seemingly Neutral Mutations

Systematic mappings of the effects of protein mutations are becoming increasingly popular. Unexpectedly, these experiments often find that proteins are tolerant to most amino acid substitutions, including substitutions in positions that are highly conserved in nature. To obtain a more realistic distribution of the effects of protein mutations, we applied a laboratory drift comprising 17 rounds of random mutagenesis and selection of M.HaeIII, a DNA methyltransferase. During this drift, multiple mutations gradually accumulated. Deep sequencing of the drifted gene ensembles allowed determination of the relative effects of all possible single nucleotide mutations. Despite being averaged across many different genetic backgrounds, about 67% of all nonsynonymous, missense mutations were evidently deleterious, and an additional 16% were likely to be deleterious. In the early generations, the frequency of most deleterious mutations remained high. However, by the 17th generation, their frequency was consistently reduced, and those remaining were accepted alongside compensatory mutations. The tolerance to mutations measured in this laboratory drift correlated with sequence exchanges seen in M.HaeIII’s natural orthologs. The biophysical constraints dictating purging in nature and in this laboratory drift also seemed to overlap. Our experiment therefore provides an improved method for measuring the effects of protein mutations that more closely replicates the natural evolutionary forces, and thereby a more realistic view of the mutational space of proteins.

Understanding and predicting the effects of single nucleotide polymorphisms (SNPs) is of fundamental importance in many fields. Systematic experimental mappings of the effects of such mutations within a given gene/protein comprise an essential experimental tool for determining protein function and for refining models of protein evolution, as well as an important resource for improving prediction algorithms. Here, we present the results of a laboratory system that mimics the manner by which protein sequences diverge in nature: a prolonged process of gradually accumulating random mutations that retain the protein’s structure and function. The change in frequencies of mutations over generations, as obtained by deep sequencing, enabled us to assess the relative effects of all possible SNPs at the background of an accumulating number of mutations. Compared to previous reports, we found that > 80% of all possible amino acid exchanges have potential deleterious effects, with 67% being clearly deleterious. Tolerance vs. purging of mutations in our prolonged drift also showed better correlation with natural diversity. Overall, our experimental setup provides a better understanding of how protein sequences diverge in nature, plus a new basis for improving the prediction accuracy of the effects of protein mutations, and specifically of SNPs.

Better prediction of functional effects for sequence variants

Elucidating the effects of naturally occurring genetic variation is one of the major challenges for personalized health and personalized medicine. Here, we introduce SNAP2, a novel neural network based classifier that improves over the state-of-the-art in distinguishing between effect and neutral variants. Our method's improved performance results from screening many potentially relevant protein features and from refining our development data sets. Cross-validated on >100k experimentally annotated variants, SNAP2 significantly outperformed other methods, attaining a two-state accuracy (effect/neutral) of 83%. SNAP2 also outperformed combinations of other methods. Performance increased for human variants but much more so for other organisms. Our method's carefully calibrated reliability index informs selection of variants for experimental follow up, with the most strongly predicted half of all effect variants predicted at over 96% accuracy. As expected, the evolutionary information from automatically generated multiple sequence alignments gave the strongest signal for the prediction. However, we also optimized our new method to perform surprisingly well even without alignments. This feature reduces prediction runtime by over two orders of magnitude, enables cross-genome comparisons, and renders our new method as the best solution for the 10-20% of sequence orphans. SNAP2 is available at: https://rostlab.org/services/snap2web

Structure-PPi: a module for the annotation of cancer-related single-nucleotide variants at protein-protein interfaces

Motivation: The interpretation of cancer-related single-nucleotide variants (SNVs) considering the protein features they affect, such as known functional sites, protein–protein interfaces, or relation with already annotated mutations, might complement the annotation of genetic variants in the analysis of NGS data. Current tools that annotate mutations fall short on several aspects, including the ability to use protein structure information or the interpretation of mutations in protein complexes.

Results: We present the Structure–PPi system for the comprehensive analysis of coding SNVs based on 3D protein structures of protein complexes. The 3D repository used, Interactome3D, includes experimental and modeled structures for proteins and protein–protein complexes. Structure–PPi annotates SNVs with features extracted from UniProt, InterPro, APPRIS, dbNSFP and COSMIC databases. We illustrate the usefulness of Structure–PPi with the interpretation of 1 027 122 non-synonymous SNVs from COSMIC and the 1000G Project that provides a collection of ∼172 700 SNVs mapped onto the protein 3D structure of 8726 human proteins (43.2% of the 20 214 SwissProt-curated proteins in UniProtKB release 2014_06) and protein–protein interfaces with potential functional implications.

Availability and implementation: Structure–PPi, along with a user manual and examples, isavailable at http://structureppi.bioinfo.cnio.es/Structure, the code for local installations at https://github.com/Rbbt-Workflows

Contact:tpons@cnio.es

Supplementary Information:Supplementary data are available at Bioinformatics online.

How structural and physicochemical determinants shape sequence constraints in a functional enzyme

The need for interfacing structural biology and biophysics to molecular evolution is being increasingly recognized. One part of the big problem is to understand how physics and chemistry shape the sequence space available to functional proteins, while satisfying the needs of biology. Here we present a quantitative, structure-based analysis of a high-resolution map describing the tolerance to all substitutions in all positions of a functional enzyme, namely a TEM lactamase previously studied through deep sequencing of mutants growing in competition experiments with selection against ampicillin. Substitutions are rarely observed within 7 Å of the active site, a stringency that is relaxed slowly and extends up to 15–20 Å, with buried residues being especially sensitive. Substitution patterns in over one third of the residues can be quantitatively modeled by monotonic dependencies on amino acid descriptors and predictions of changes in folding stability. Amino acid volume and steric hindrance shape constraints on the protein core; hydrophobicity and solubility shape constraints on hydrophobic clusters underneath the surface, and on salt bridges and polar networks at the protein surface together with charge and hydrogen bonding capacity. Amino acid solubility, flexibility and conformational descriptors also provide additional constraints at many locations. These findings provide fundamental insights into the chemistry underlying protein evolution and design, by quantitating links between sequence and different protein traits, illuminating subtle and unexpected sequence-trait relationships and pinpointing what traits are sacrificed upon gain-of-function mutation.