Predicting folding free energy changes upon single point mutations

Enhancing thermostability of tryptophan hydroxylase via protein engineering and its application in 5-hydroxytryptophan production

5-Hydroxytryptophan (5-HTP), as the precursor of serotonin and melatonin in animals, can regulate mood, sleep, and behavior, which is widely used in pharmaceutical and health products industry. The enzymatic production of 5-hydroxytryptophan (5-HTP) from L-tryptophan (L-Trp) using tryptophan hydroxylase (TPH) show huge potential in application due to its advantages, such as mild reaction conditions, avoidance of protection/deprotection processes, excellent regioselectivity and considerable catalytic efficiency, compared with chemical synthesis and natural extraction. However, the low thermostability of TPH restricted its hydroxylation efficiency toward L-Trp. In this study, we aimed to improve the thermostability of TPH via semi-rational design guided by (folding free energy) ΔΔG fold calculation. After two rounds of evolution, two beneficial mutants M1 (S422V) and M30 (V275L/I412K) were obtained. Thermostability evaluation showed that M1 and M30 possessed 5.66-fold and 6.32-fold half-lives (t1/2) at 37 °C, and 4.2 °C and 6.0 °C higher melting temperature (Tm) than the WT, respectively. The mechanism behind thermostability improvement was elucidated with molecular dynamics simulation. Furthermore, biotransformation of 5-HTP from L-Trp was performed, M1 and M30 displayed 1.80-fold and 2.30-fold than that of WT, respectively. This work provides important insights into the thermostability enhancement of TPH and generate key mutants that could be robust candidates for practical production of 5-HTP.

Exploring Free Energies of Specific Protein Conformations Using the Martini Force Field

Coarse-grained (CG) level molecular dynamics simulations are routinely used to study various biomolecular processes. The Martini force field is currently the most widely adopted parameter set for such simulations. The functional form of this and several other CG force fields enforces secondary protein structure support by employing a variety of harmonic potentials or restraints that favor the protein's native conformation. We propose a straightforward method to calculate the energetic consequences of transitions between predefined conformational states in systems in which multiple factors can affect protein conformational equilibria. This method is designed for use within the Martini force field and involves imposing conformational transitions by linking a Martini-inherent elastic network to the coupling parameter λ. We demonstrate the applicability of our method using the example of five biomolecular systems that undergo experimentally characterized conformational transitions between well-defined structures (Staphylococcal nuclease, C-terminal segment of surfactant protein B, LAH4 peptide, and β2-adrenergic receptor) as well as between folded and unfolded states (GCN4 leucine zipper protein). The results show that the relative free energy changes associated with protein conformational transitions, which are affected by various factors, such as pH, mutations, solvent, and lipid membrane composition, are correctly reproduced. The proposed method may be a valuable tool for understanding how different conditions and modifications affect conformational equilibria in proteins.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Bioinformatic analysis of the effect of SNPs in the pig TERT gene on the structural and functional characteristics of the enzyme to develop new genetic markers of productivity traits

BACKGROUND

Telomerase reverse transcriptase (TERT) plays a crucial role in synthesizing telomeric repeats that safeguard chromosomes from damage and fusion, thereby maintaining genome stability. Mutations in the TERT gene can lead to a deviation in gene expression, impaired enzyme activity, and, as a result, abnormal telomere shortening. Genetic markers of productivity traits in livestock can be developed based on the TERT gene polymorphism for use in marker-associated selection (MAS). In this study, a bioinformatic-based approach is proposed to evaluate the effect of missense single-nucleotide polymorphisms (SNPs) in the pig TERT gene on enzyme function and structure, with the prospect of developing genetic markers.

RESULTS

A comparative analysis of the coding and amino acid sequences of the pig TERT was performed with corresponding sequences of other species. The distribution of polymorphisms in the pig TERT gene, with respect to the enzyme's structural-functional domains, was established. A three-dimensional model of the pig TERT structure was obtained through homological modeling. The potential impact of each of the 23 missense SNPs in the pig TERT gene on telomerase function and stability was assessed using predictive bioinformatic tools utilizing data on the amino acid sequence and structure of pig TERT.

CONCLUSIONS

According to bioinformatic analysis of 23 missense SNPs of the pig TERT gene, a predictive effect of rs789641834 (TEN domain), rs706045634 (TEN domain), rs325294961 (TRBD domain) and rs705602819 (RTD domain) on the structural and functional parameters of the enzyme was established. These SNPs hold the potential to serve as genetic markers of productivity traits. Therefore, the possibility of their application in MAS should be further evaluated in associative analysis studies.

Characterizing the combined effects of cytochrome P450 missense variation within star allele definitions

Background: Cytochrome P450 (CYP) genetic variation largely impacts drug response. However, many CYP star alleles (haplotypes) lack functional annotation, impeding our understanding of drug metabolism mechanisms. We aimed to investigate the impact of missense variant combinations on CYP protein structures. Methods: Normal mode analysis was conducted on 261 missense variants within 91 CYP haplotypes. CYP2D6*2 and CYP2D6*17 were prioritized for molecular dynamics simulation. Results: Normal mode analysis and molecular dynamics highlight the effects of known CYP missense variants on protein stability and conformational dynamics. Missense variants within haplotypes may have intermodulating effects on protein structure and function. Conclusion: This study highlights the utility of multiscale modeling in interpreting CYP missense variants and particularly their combinations within various star alleles.

An efficient method to predict protein thermostability in alanine mutation

The relationship between protein sequence and its thermodynamic stability is a critical aspect of computational protein design. In this work, we present a new theoretical method to calculate the free energy change (ΔΔG) resulting from a single-point amino acid mutation to alanine in a protein sequence. The method is derived based on physical interactions and is very efficient in estimating the free energy changes caused by a series of alanine mutations from just a single molecular dynamics (MD) trajectory. Numerical calculations are carried out on a total of 547 alanine mutations in 19 diverse proteins whose experimental results are available. The comparison between the experimental ΔΔG_exp and the calculated values shows a generally good correlation with a correlation coefficient of 0.67. Both the advantages and limitations of this method are discussed. This method provides an efficient and valuable tool for protein design and engineering.

Stability profile of the neuronal SNARE complex reflects its potency to drive fast membrane fusion

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) form the SNARE complex to mediate most fusion events of the secretory pathway. The neuronal SNARE complex is featured by its high stability and half-zippered conformation required for driving robust and fast synaptic exocytosis. However, these two features seem to be thermodynamically mutually exclusive. In this study, we have employed temperature-dependent disassociation assays and single-molecule Förster resonance energy transfer (FRET) experiments to analyze the stability and conformation of the neuronal SNARE complex. We reclassified the amino acids of the SNARE motif into four sub-groups (core, core-side I and II, and non-contact). Our data showed that the core residues predominantly contribute to the complex stability to meet a basal requirement for SNARE-mediated membrane fusion, while the core-side residues exert an unbalanced effect on the N- and C-half bundle stability that determines the half-zippered conformation of the neuronal SNARE complex, which would accommodate essential regulations by complexins and synaptotagmins for fast Ca2+-triggered membrane fusion. Furthermore, our data confirmed a strong coupling of folding energy between the N- and C-half assembly of the neuronal SNARE complex, which rationalizes the strong potency of the half-zippered conformation to conduct robust and fast fusion. Overall, these results uncovered that the stability profile of the neuronal SNARE complex reflects its potency to drive fast and robust membrane fusion. Based on these results, we also developed a new parameter, the stability factor (Fs), to characterize the overall stability of the neuronal SNARE complex and resolved a linear correlation between the stability and inter-residue coulombic interactions of the neuronal SNARE complex, which would help rationally design artificial SNARE complexes and remold functional SNARE complexes with desirable stability.

Thermostability engineering of industrial enzymes through structure modification

Thermostability is an essential requirement of enzymes in the industrial processes to catalyze the reactions at high temperatures; thus, enzyme engineering through directed evolution, semi-rational design and rational design are commonly employed to construct desired thermostable mutants. Several strategies are implemented to fulfill enzymes' thermostability demand including decreasing the entropy of the unfolded state through substitutions Gly → Xxx or Xxx → Pro, hydrogen bond, salt bridge, introducing two different simultaneous interactions through single mutant, hydrophobic interaction, filling the hydrophobic cavity core, decreasing surface hydrophobicity, truncating loop, aromatic-aromatic interaction and introducing positively charged residues to enzyme surface. In the current review, horizons about compatibility between secondary structures and substitutions at preferable structural positions to generate the most desirable thermostability in industrial enzymes are broadened. KEY POINTS: • Protein engineering is a powerful tool for generating thermostable industrial enzymes. • Directed evolution and rational design are practical approaches in enzyme engineering. • Substitutions in preferable structural positions can increase thermostability.

Dissecting mutational allosteric effects in alkaline phosphatases associated with different Hypophosphatasia phenotypes: An integrative computational investigation

Hypophosphatasia (HPP) is a rare inherited disorder characterized by defective bone mineralization and is highly variable in its clinical phenotype. The disease occurs due to various loss-of-function mutations in ALPL, the gene encoding tissue-nonspecific alkaline phosphatase (TNSALP). In this work, a data-driven and biophysics-based approach is proposed for the large-scale analysis of ALPL mutations-from nonpathogenic to severe HPPs. By using a pipeline of synergistic approaches including sequence-structure analysis, network modeling, elastic network models and atomistic simulations, we characterized allosteric signatures and effects of the ALPL mutations on protein dynamics and function. Statistical analysis of molecular features computed for the ALPL mutations showed a significant difference between the control, mild and severe HPP phenotypes. Molecular dynamics simulations coupled with protein structure network analysis were employed to analyze the effect of single-residue variation on conformational dynamics of TNSALP dimers, and the developed machine learning model suggested that the topological network parameters could serve as a robust indicator of severe mutations. The results indicated that the severity of disease-associated mutations is often linked with mutation-induced modulation of allosteric communications in the protein. This study suggested that ALPL mutations associated with mild and more severe HPPs can exert markedly distinct effects on the protein stability and long-range network communications. By linking the disease phenotypes with dynamic and allosteric molecular signatures, the proposed integrative computational approach enabled to characterize and quantify the allosteric effects of ALPL mutations and role of allostery in the pathogenesis of HPPs.

Exploring the Effect of Enhanced Sampling on Protein Stability Prediction

Changes in protein stability due to side-chain mutations are evaluated by alchemical free-energy calculations based on classical molecular dynamics (MD) simulations in explicit solvent using the GROMOS force field. Three proteins of different complexity with a total number of 93 single-point mutations are analyzed, and the relative free-energy differences are discussed with respect to configurational sampling and (dis)agreement with experimental data. For the smallest protein studied, a 34-residue WW domain, the starting structure dependence of the alchemical free-energy changes, is discussed in detail. Deviations from previous simulations for the two other proteins are shown to result from insufficient sampling in the earlier studies. Hamiltonian replica exchange in combination with multiple starting structures and sufficient sampling time of more than 100 ns per intermediate alchemical state is required in some cases to reach convergence.

Critical assessment of structure-based approaches to improve protein resistance in aqueous ionic liquids by enzyme-wide saturation mutagenesis

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for green industrial processes involving biocatalysts, but often reduce enzyme activity. Experimental and computational methods are applied to predict favorable substitution sites and, most often, subsequent site-directed surface charge modifications are introduced to enhance enzyme resistance towards aIL. However, almost no studies evaluate the prediction precision with random mutagenesis or the application of simple data-driven filtering processes. Here, we systematically and rigorously evaluated the performance of 22 previously described structure-based approaches to increase enzyme resistance to aIL based on an experimental complete site-saturation mutagenesis library of Bacillus subtilis Lipase A (BsLipA) screened against four aIL. We show that, surprisingly, most of the approaches yield low gain-in-precision (GiP) values, particularly for predicting relevant positions: 14 approaches perform worse than random mutagenesis. Encouragingly, exploiting experimental information on the thermostability of BsLipA or structural weak spots of BsLipA predicted by rigidity theory yields GiP = 3.03 and 2.39 for relevant variants and GiP = 1.61 and 1.41 for relevant positions. Combining five simple-to-compute physicochemical and evolutionary properties substantially increases the precision of predicting relevant variants and positions, yielding GiP = 3.35 and 1.29. Finally, combining these properties with predictions of structural weak spots identified by rigidity theory additionally improves GiP for relevant variants up to 4-fold to ∼10 and sustains or increases GiP for relevant positions, resulting in a prediction precision of ∼90% compared to ∼9% in random mutagenesis. This combination should be applicable to other enzyme systems for guiding protein engineering approaches towards improved aIL resistance.

A guide to machine learning for biologists

The expanding scale and inherent complexity of biological data have encouraged a growing use of machine learning in biology to build informative and predictive models of the underlying biological processes. All machine learning techniques fit models to data; however, the specific methods are quite varied and can at first glance seem bewildering. In this Review, we aim to provide readers with a gentle introduction to a few key machine learning techniques, including the most recently developed and widely used techniques involving deep neural networks. We describe how different techniques may be suited to specific types of biological data, and also discuss some best practices and points to consider when one is embarking on experiments involving machine learning. Some emerging directions in machine learning methodology are also discussed.

LYRUS: a machine learning model for predicting the pathogenicity of missense variants

SUMMARY

Single amino acid variations (SAVs) are a primary contributor to variations in the human genome. Identifying pathogenic SAVs can provide insights to the genetic architecture of complex diseases. Most approaches for predicting the functional effects or pathogenicity of SAVs rely on either sequence or structural information. This study presents 〈Lai Yang Rubenstein Uzun Sarkar〉 (LYRUS), a machine learning method that uses an XGBoost classifier to predict the pathogenicity of SAVs. LYRUS incorporates five sequence-based, six structure-based and four dynamics-based features. Uniquely, LYRUS includes a newly proposed sequence co-evolution feature called the variation number. LYRUS was trained using a dataset that contains 4363 protein structures corresponding to 22 639 SAVs from the ClinVar database, and tested using the VariBench testing dataset. Performance analysis showed that LYRUS achieved comparable performance to current variant effect predictors. LYRUS's performance was also benchmarked against six Deep Mutational Scanning datasets for PTEN and TP53.

AVAILABILITY AND IMPLEMENTATION

LYRUS is freely available and the source code can be found at https://github.com/jiaying2508/LYRUS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Advances online.

Overview of strategies for developing high thermostability industrial enzymes: Discovery, mechanism, modification and challenges

Biocatalysts such as enzymes are environmentally friendly and have substrate specificity, which are preferred in the production of various industrial products. However, the strict reaction conditions in industry including high temperature, organic solvents, strong acids and bases and other harsh environments often destabilize enzymes, and thus substantially compromise their catalytic functions, and greatly restrict their applications in food, pharmaceutical, textile, bio-refining and feed industries. Therefore, developing industrial enzymes with high thermostability becomes very important in industry as thermozymes have more advantages under high temperature. Discovering new thermostable enzymes using genome sequencing, metagenomics and sample isolation from extreme environments, or performing molecular modification of the existing enzymes with poor thermostability using emerging protein engineering technology have become an effective means of obtaining thermozymes. Based on the thermozymes as biocatalytic chips in industry, this review systematically analyzes the ways to discover thermostable enzymes from extreme environment, clarifies various interaction forces that will affect thermal stability of enzymes, and proposes different strategies to improve enzymes' thermostability. Furthermore, latest development in the thermal stability modification of industrial enzymes through rational design strategies is comprehensively introduced from structure-activity relationship point of view. Challenges and future research perspectives are put forward as well.

Computational and synthetic biology approaches for the biosynthesis of antiviral and anticancer terpenoids from Bacillus subtilis

Recent advancements in medicinal research have identified several antiviral and anticancer terpenoids that are usually deployed as a source of flavor, fragrances and pharmaceuticals. Under the current COVID-19 pandemic conditions, natural therapeutics with least side effects are the need of the hour to save the patients, especially, which are pre-affected with other medical complications. Although, plants are the major sources of terpenoids; however, for the environmental concerns, the global interest has shifted to the biocatalytic production of molecules from microbial sources. The gram-positive bacterium Bacillus subtilis is a suitable host in this regard due to its GRAS (generally regarded as safe) status, ease in genetic manipulations and wide industrial acceptability. The B. subtilis synthesizes its terpenoid molecules from 1-deoxy-d-xylulose-5-phosphate (DXP) pathway, a common route in almost all microbial strains. Here, we summarize the computational and synthetic biology approaches to improve the production of terpenoid-based therapeutics from B. subtilis by utilizing DXP pathway. We focus on the in-silico approaches for screening the functionally improved enzyme-variants of the two crucial enzymes namely, the DXP synthase (DXS) and farnesyl pyrophosphate synthase (FPPS). The approaches for engineering the active sites are subsequently explained. It will be helpful to construct the functionally improved enzymes for the high-yield production of terpenoid-based anticancer and antiviral metabolites, which would help to reduce the cost and improve the availability of such therapeutics for the humankind.

Analysis and Interpretation of the Impact of Missense Variants in Cancer

Large scale genome sequencing allowed the identification of a massive number of genetic variations, whose impact on human health is still unknown. In this review we analyze, by an in silico-based strategy, the impact of missense variants on cancer-related genes, whose effect on protein stability and function was experimentally determined. We collected a set of 164 variants from 11 proteins to analyze the impact of missense mutations at structural and functional levels, and to assess the performance of state-of-the-art methods (FoldX and Meta-SNP) for predicting protein stability change and pathogenicity. The result of our analysis shows that a combination of experimental data on protein stability and in silico pathogenicity predictions allowed the identification of a subset of variants with a high probability of having a deleterious phenotypic effect, as confirmed by the significant enrichment of the subset in variants annotated in the COSMIC database as putative cancer-driving variants. Our analysis suggests that the integration of experimental and computational approaches may contribute to evaluate the risk for complex disorders and develop more effective treatment strategies.

Predicting rifampicin resistance mutations in bacterial RNA polymerase subunit beta based on majority consensus

BACKGROUND

Mutations in an enzyme target are one of the most common mechanisms whereby antibiotic resistance arises. Identification of the resistance mutations in bacteria is essential for understanding the structural basis of antibiotic resistance and design of new drugs. However, the traditionally used experimental approaches to identify resistance mutations were usually labor-intensive and costly.

RESULTS

We present a machine learning (ML)-based classifier for predicting rifampicin (Rif) resistance mutations in bacterial RNA Polymerase subunit β (RpoB). A total of 186 mutations were gathered from the literature for developing the classifier, using 80% of the data as the training set and the rest as the test set. The features of the mutated RpoB and their binding energies with Rif were calculated through computational methods, and used as the mutation attributes for modeling. Classifiers based on five ML algorithms, i.e. decision tree, k nearest neighbors, naïve Bayes, probabilistic neural network and support vector machine, were first built, and a majority consensus (MC) approach was then used to obtain a new classifier based on the classifications of the five individual ML algorithms. The MC classifier comprehensively improved the predictive performance, with accuracy, F-measure and AUC of 0.78, 0.83 and 0.81for training set whilst 0.84, 0.87 and 0.83 for test set, respectively.

CONCLUSION

The MC classifier provides an alternative methodology for rapid identification of resistance mutations in bacteria, which may help with early detection of antibiotic resistance and new drug discovery.

A Novel SEMA3G Mutation in Two Siblings Affected by Syndromic GnRH Deficiency

INTRODUCTION

Gonadotropin-releasing hormone (GnRH) deficiency causes hypogonadotropic hypogonadism (HH), a rare genetic disorder that impairs sexual reproduction. HH can be due to defective GnRH-secreting neuron development or function and may be associated with other clinical signs in overlapping genetic syndromes. With most of the cases being idiopathic, genetics underlying HH is still largely unknown.

OBJECTIVE

To assess the contribution of mutated Semaphorin 3G (SEMA3G) in the onset of a syndromic form of HH, characterized by intellectual disability and facial dysmorphic features.

METHOD

By combining homozygosity mapping with exome sequencing, we identified a novel variant in the SEMA3G gene. We then applied mouse as a model organism to examine SEMA3Gexpression and its functional requirement in vivo. Further, we applied homology modelling in silico and cell culture assays in vitro to validate the pathogenicity of the identified gene variant.

RESULTS

We found that (i) SEMA3G is expressed along the migratory route of GnRH neurons and in the developing pituitary, (ii) SEMA3G affects GnRH neuron development, but is redundant in the adult hypothalamic-pituitary-gonadal axis, and (iii) mutated SEMA3G alters binding properties in silico and in vitro to its PlexinA receptors and attenuates its effect on the migration of immortalized GnRH neurons.

CONCLUSION

In silico, in vitro, and in vivo models revealed that SEMA3G regulates GnRH neuron migration and that its mutation affecting receptor selectivity may be responsible for the HH-related defects.

PremPS: Predicting the impact of missense mutations on protein stability

Computational methods that predict protein stability changes induced by missense mutations have made a lot of progress over the past decades. Most of the available methods however have very limited accuracy in predicting stabilizing mutations because existing experimental sets are dominated by mutations reducing protein stability. Moreover, few approaches could consistently perform well across different test cases. To address these issues, we developed a new computational method PremPS to more accurately evaluate the effects of missense mutations on protein stability. The PremPS method is composed of only ten evolutionary- and structure-based features and parameterized on a balanced dataset with an equal number of stabilizing and destabilizing mutations. A comprehensive comparison of the predictive performance of PremPS with other available methods on nine benchmark datasets confirms that our approach consistently outperforms other methods and shows considerable improvement in estimating the impacts of stabilizing mutations. A protein could have multiple structures available, and if another structure of the same protein is used, the predicted change in stability for structure-based methods might be different. Thus, we further estimated the impact of using different structures on prediction accuracy, and demonstrate that our method performs well across different types of structures except for low-resolution structures and models built based on templates with low sequence identity. PremPS can be used for finding functionally important variants, revealing the molecular mechanisms of functional influences and protein design. PremPS is freely available at https://lilab.jysw.suda.edu.cn/research/PremPS/, which allows to do large-scale mutational scanning and takes about four minutes to perform calculations for a single mutation per protein with ~ 300 residues and requires ~ 0.4 seconds for each additional mutation.

The development of computational methods to accurately predict the impacts of amino acid substitutions on protein stability is of paramount importance for the field of protein design and understanding the roles of missense mutations in disease. However, most of the available methods have very limited predictive accuracy for mutations increasing stability and few could consistently perform well across different test cases. Here we present a new computational approach PremPS, which is capable of predicting the effects of single point mutations on protein stability. PremPS employs only ten evolutionary- and structure-based features and is trained on a symmetrical dataset consisting of the same number of cases of stabilizing and destabilizing mutations. Our method was tested against numerous blind datasets and shows a considerable improvement especially in evaluating the effects of stabilizing mutations, outperforming previously developed methods. PremPS is freely available as a user-friendly web server at http://lilab.jysw.suda.edu.cn/research/PremPS/, which is fast enough to handle the large number of cases.

Variation benchmark datasets: update, criteria, quality and applications

Development of new computational methods and testing their performance has to be carried out using experimental data. Only in comparison to existing knowledge can method performance be assessed. For that purpose, benchmark datasets with known and verified outcome are needed. High-quality benchmark datasets are valuable and may be difficult, laborious and time consuming to generate. VariBench and VariSNP are the two existing databases for sharing variation benchmark datasets used mainly for variation interpretation. They have been used for training and benchmarking predictors for various types of variations and their effects. VariBench was updated with 419 new datasets from 109 papers containing altogether 329 014 152 variants; however, there is plenty of redundancy between the datasets. VariBench is freely available at http://structure.bmc.lu.se/VariBench/. The contents of the datasets vary depending on information in the original source. The available datasets have been categorized into 20 groups and subgroups. There are datasets for insertions and deletions, substitutions in coding and non-coding region, structure mapped, synonymous and benign variants. Effect-specific datasets include DNA regulatory elements, RNA splicing, and protein property for aggregation, binding free energy, disorder and stability. Then there are several datasets for molecule-specific and disease-specific applications, as well as one dataset for variation phenotype effects. Variants are often described at three molecular levels (DNA, RNA and protein) and sometimes also at the protein structural level including relevant cross references and variant descriptions. The updated VariBench facilitates development and testing of new methods and comparison of obtained performances to previously published methods. We compared the performance of the pathogenicity/tolerance predictor PON-P2 to several benchmark studies, and show that such comparisons are feasible and useful, however, there may be limitations due to lack of provided details and shared data.

Database URL: http://structure.bmc.lu.se/VariBench

Whole Genome In-Silico Analysis of South African G1P[8] Rotavirus Strains Before and After Vaccine Introduction Over A Period of 14 Years

Rotavirus G1P[8] strains account for more than half of the group A rotavirus (RVA) infections in children under five years of age, globally. A total of 103 stool samples previously characterized as G1P[8] and collected seven years before and seven years after introducing the Rotarix® vaccine in South Africa were processed for whole-genome sequencing. All the strains analyzed had a Wa-like constellation (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1). South African pre- and post-vaccine G1 strains were clustered in G1 lineage-I and II while the majority (84.2%) of the P[8] strains were grouped in P[8] lineage-III. Several amino acid sites across ten gene segments with the exception of VP7 were under positive selective pressure. Except for the N147D substitution in the antigenic site of eight post-vaccine G1 strains when compared to both Rotarix® and pre-vaccine strains, most of the amino acid substitutions in the antigenic regions of post-vaccine G1P[8] strains were already present during the pre-vaccine period. Therefore, Rotarix® did not appear to have an impact on the amino acid differences in the antigenic regions of South African post-vaccine G1P[8] strains. However, continued whole-genome surveillance of RVA strains to decipher genetic changes in the post-vaccine period remains imperative.

Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior

Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZA^R) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZA^R mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 μg ml⁻¹). The PZA^R samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZA^R for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZA^R TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level.

We performed several analyses to unveil the pyrazinamide-resistance mechanism of R123P, T76P, G150A, and H71R mutants through molecular dynamics simulations.

Computational analysis of androgen receptor (AR) variants to decipher the relationship between protein stability and related-diseases

Although more than 1,000 androgen receptor (AR) mutations have been identified and these mutants are pathologically important, few theoretical studies have investigated the role of AR protein folding stability in disease and its relationship with the phenotype of the patients. Here, we extracted AR variant data from four databases: ARDB, HGMD, Cosmic, and 1,000 genome. 905 androgen insensitivity syndrome (AIS)-associated loss-of-function mutants and 168 prostate cancer-associated gain-of-function mutants in AR were found. We analyzed the effect of single-residue variation on the folding stability of AR by FoldX and guanidine hydrochloride denaturation experiment, and found that genetic disease-associated mutations tend to have a significantly greater effect on protein stability than gene polymorphisms. Moreover, AR mutants in complete androgen insensitivity syndrome (CAIS) tend to have a greater effect on protein stability than in partial androgen insensitive syndrome (PAIS). This study, by linking disease phenotypes to changes in AR stability, demonstrates the importance of protein stability in the pathogenesis of hereditary disease.

Mutations in transmembrane proteins: diseases, evolutionary insights, prediction and comparison with globular proteins

Membrane proteins are unique in that they interact with lipid bilayers, making them indispensable for transporting molecules and relaying signals between and across cells. Due to the significance of the protein's functions, mutations often have profound effects on the fitness of the host. This is apparent both from experimental studies, which implicated numerous missense variants in diseases, as well as from evolutionary signals that allow elucidating the physicochemical constraints that intermembrane and aqueous environments bring. In this review, we report on the current state of knowledge acquired on missense variants (referred to as to single amino acid variants) affecting membrane proteins as well as the insights that can be extrapolated from data already available. This includes an overview of the annotations for membrane protein variants that have been collated within databases dedicated to the topic, bioinformatics approaches that leverage evolutionary information in order to shed light on previously uncharacterized membrane protein structures or interaction interfaces, tools for predicting the effects of mutations tailored specifically towards the characteristics of membrane proteins as well as two clinically relevant case studies explaining the implications of mutated membrane proteins in cancer and cardiomyopathy.

Systematically Scrutinizing the Impact of Substitution Sites on Thermostability and Detergent Tolerance for Bacillus subtilis Lipase A

Improving an enzyme's (thermo-)stability or tolerance against solvents and detergents is highly relevant in protein engineering and biotechnology. Recent developments have tended toward data-driven approaches, where available knowledge about the protein is used to identify substitution sites with high potential to yield protein variants with improved stability, and subsequently, substitutions are engineered by site-directed or site-saturation (SSM) mutagenesis. However, the development and validation of algorithms for data-driven approaches have been hampered by the lack of availability of large-scale data measured in a uniform way and being unbiased with respect to substitution types and locations. Here, we extend our knowledge on guidelines for protein engineering following a data-driven approach by scrutinizing the impact of substitution sites on thermostability or/and detergent tolerance for Bacillus subtilis lipase A (BsLipA) at very large scale. We systematically analyze a complete experimental SSM library of BsLipA containing all 3439 possible single variants, which was evaluated as to thermostability and tolerances against four detergents under respectively uniform conditions. Our results provide systematic and unbiased reference data at unprecedented scale for a biotechnologically important protein, identify consistently defined hot spot types for evaluating the performance of data-driven protein-engineering approaches, and show that the rigidity theory and ensemble-based approach Constraint Network Analysis yields hot spot predictions with an up to ninefold gain in precision over random classification.

Approaching protein design with multisite λ dynamics: Accurate and scalable mutational folding free energies in T4 lysozyme

The estimation of changes in free energy upon mutation is central to the problem of protein design. Modern protein design methods have had remarkable success over a wide range of design targets, but are reaching their limits in ligand binding and enzyme design due to insufficient accuracy in mutational free energies. Alchemical free energy calculations have the potential to supplement modern design methods through more accurate molecular dynamics based prediction of free energy changes, but suffer from high computational cost. Multisite λ dynamics (MSλD) is a particularly efficient and scalable free energy method with potential to explore combinatorially large sequence spaces inaccessible with other free energy methods. This work aims to quantify the accuracy of MSλD and demonstrate its scalability. We apply MSλD to the classic problem of calculating folding free energies in T4 lysozyme, a system with a wealth of experimental measurements. Single site mutants considering 32 mutations show remarkable agreement with experiment with a Pearson correlation of 0.914 and mean unsigned error of 1.19 kcal/mol. Multisite mutants in systems with up to five concurrent mutations spanning 240 different sequences show comparable agreement with experiment. These results demonstrate the promise of MSλD in exploring large sequence spaces for protein design.

Structural Perspective on Revealing and Altering Molecular Functions of Genetic Variants Linked with Diseases

Structural information of biological macromolecules is crucial and necessary to deliver predictions about the effects of mutations-whether polymorphic or deleterious (i.e., disease causing), wherein, thermodynamic parameters, namely, folding and binding free energies potentially serve as effective biomarkers. It may be emphasized that the effect of a mutation depends on various factors, including the type of protein (globular, membrane or intrinsically disordered protein) and the structural context in which it occurs. Such information may positively aid drug-design. Furthermore, due to the intrinsic plasticity of proteins, even mutations involving radical change of the structural and physico⁻chemical properties of the amino acids (native vs. mutant) can still have minimal effects on protein thermodynamics. However, if a mutation causes significant perturbation by either folding or binding free energies, it is quite likely to be deleterious. Mitigating such effects is a promising alternative to the traditional approaches of designing inhibitors. This can be done by structure-based in silico screening of small molecules for which binding to the dysfunctional protein restores its wild type thermodynamics. In this review we emphasize the effects of mutations on two important biophysical properties, stability and binding affinity, and how structures can be used for structure-based drug design to mitigate the effects of disease-causing variants on the above biophysical properties.

In silico mutational analysis and identification of stability centers in human interleukin-4

Interleukin-4 (IL-4) is a multifunctional cytokine that plays a critical role in apoptosis, differentiation and proliferation. The intensity of IL4 response depends upon binding to its receptor, IL-4R. The therapeutic efficiency of interleukins can be increased by generating structural mutants having greater stability. In the present work, attempts were made to increase the stability of human IL-4 using in-silico site directed mutagenesis. Different orthologous sequences of IL4 from Pan troglodytes, Aotusnigriceps, Macacamulatta, Papiohamadryas, Chlorocebusaethiops, Vicugnapacos, Susscrofa and Homo sapiens were aligned using Clustal Omega that revealed the conserved and non-conserved positions. For each non-conserved position, possible favorable and stabilizing mutations were found using CUPSAT with predicted ΔΔG (kcal/mol). The one with highest ΔΔG (kcal/mol) among all possible mutations, for each non-conserved position was selected and introduced manually in human IL-4 sequence resulting in multiple mutants of IL-4. Mutant proteins were modeled using structure of IL4 (PDB ID: 2B8U) as a template by SWISS MODEL. The mutants A49L and Q106T were identified to have stability centre using SCide. Molecular dynamics and docking analysis also confirmed the mutants stability and binding respectively. Mutants A49L and Q106T had -7.580079 kcal/mol and -39.418124 kcal/mol respectively lesser energy value than the wild type IL4. The result suggested that, the stability of human IL-4 has been increased by mutation.

Computational Approaches to Prioritize Cancer Driver Missense Mutations

Cancer is a complex disease that is driven by genetic alterations. There has been a rapid development of genome-wide techniques during the last decade along with a significant lowering of the cost of gene sequencing, which has generated widely available cancer genomic data. However, the interpretation of genomic data and the prediction of the association of genetic variations with cancer and disease phenotypes still requires significant improvement. Missense mutations, which can render proteins non-functional and provide a selective growth advantage to cancer cells, are frequently detected in cancer. Effects caused by missense mutations can be pinpointed by in silico modeling, which makes it more feasible to find a treatment and reverse the effect. Specific human phenotypes are largely determined by stability, activity, and interactions between proteins and other biomolecules that work together to execute specific cellular functions. Therefore, analysis of missense mutations’ effects on proteins and their complexes would provide important clues for identifying functionally important missense mutations, understanding the molecular mechanisms of cancer progression and facilitating treatment and prevention. Herein, we summarize the major computational approaches and tools that provide not only the classification of missense mutations as cancer drivers or passengers but also the molecular mechanisms induced by driver mutations. This review focuses on the discussion of annotation and prediction methods based on structural and biophysical data, analysis of somatic cancer missense mutations in 3D structures of proteins and their complexes, predictions of the effects of missense mutations on protein stability, protein-protein and protein-nucleic acid interactions, and assessment of conformational changes in protein conformations induced by mutations.

Atomistic simulations and network-based modeling of the Hsp90-Cdc37 chaperone binding with Cdk4 client protein: A mechanism of chaperoning kinase clients by exploiting weak spots of intrinsically dynamic kinase domains

A fundamental role of the Hsp90 and Cdc37 chaperones in mediating conformational development and activation of diverse protein kinase clients is essential in signal transduction. There has been increasing evidence that the Hsp90-Cdc37 system executes its chaperoning duties by recognizing conformational instability of kinase clients and modulating their folding landscapes. The recent cryo-electron microscopy structure of the Hsp90-Cdc37-Cdk4 kinase complex has provided a framework for dissecting regulatory principles underlying differentiation and recruitment of protein kinase clients to the chaperone machinery. In this work, we have combined atomistic simulations with protein stability and network-based rigidity decomposition analyses to characterize dynamic factors underlying allosteric mechanism of the chaperone-kinase cycle and identify regulatory hotspots that control client recognition. Through comprehensive characterization of conformational dynamics and systematic identification of stabilization centers in the unbound and client- bound Hsp90 forms, we have simulated key stages of the allosteric mechanism, in which Hsp90 binding can induce instability and partial unfolding of Cdk4 client. Conformational landscapes of the Hsp90 and Cdk4 structures suggested that client binding can trigger coordinated dynamic changes and induce global rigidification of the Hsp90 inter-domain regions that is coupled with a concomitant increase in conformational flexibility of the kinase client. This process is allosteric in nature and can involve reciprocal dynamic exchanges that exert global effect on stability of the Hsp90 dimer, while promoting client instability. The network-based rigidity analysis and emulation of thermal unfolding of the Cdk4-cyclin D complex and Hsp90-Cdc37-Cdk4 complex revealed weak spots of kinase instability that are present in the native Cdk4 structure and are targeted by the chaperone during client recruitment. Our findings suggested that this mechanism may be exploited by the Hsp90-Cdc37 chaperone to recruit and protect intrinsically dynamic kinase clients from degradation. The results of this investigation are discussed and interpreted in the context of diverse experimental data, offering new insights into mechanisms of chaperone regulation and binding.

A component analysis of the free energies of folding of 35 proteins: A consensus view on the thermodynamics of folding at the molecular level

What factors favor protein folding? This is a textbook question. Parsing the experimental free energies of folding/unfolding into diverse enthalpic and entropic components of solute and solvent favoring or disfavoring folding is not an easy task. In this study, we present a computational protocol for estimating the free energy contributors to protein folding semi-quantitatively using ensembles of unfolded and native states generated via molecular dynamics simulations. We tested the methodology on 35 proteins with diverse structural motifs and sizes and found that the calculated free energies correlate well with experiment (correlation coefficient ∼ 0.85), enabling us to develop a consensus view of the energetics of folding. As a more sensitive test of the methodology, we also investigated the free energies of folding of an additional 33 single point mutants and obtained a correlation coefficient of 0.8. A notable observation is that the folding free energy components appear to carry signatures of the fold (SCOP classification) of the protein. © 2017 Wiley Periodicals, Inc.

Improvement thermostability of Pseudoalteromonas carrageenovora arylsulfatase by rational design

This study aimed to improve the thermostability of arylsulfatase from Pseudoalteromonas carrageenovora. A total of 10 single-site mutants were chosen using the PoPMuSiC program, and two mutants of K253N and P314T showed enhanced thermal stability. By saturation mutagenesis and thermostability analysis, K253H and P314T were the best mutants at the two sites. Combinational mutations of K253H, P314T and H260L were subsequently introduced, and the best mutant of K253H/H260L was selected. Thermal inactivation analysis showed the half-life (t_1/2) value at 55°C for K253H/H260L was 7.7-fold that of the wild-type enzyme (WT), meanwhile this mutant maintained the specific enzyme activity. Structure modeling demonstrated that the additional hydrogen bonds, optimization of surface charge-charge interactions, and increasing of hydrophobic interaction could account for the improved thermostability imparted by K253H/H260L.

Annotating Mutational Effects on Proteins and Protein Interactions: Designing Novel and Revisiting Existing Protocols

In this review we describe a protocol to annotate the effects of missense mutations on proteins, their functions, stability, and binding. For this purpose we present a collection of the most comprehensive databases which store different types of sequencing data on missense mutations, we discuss their relationships, possible intersections, and unique features. Next, we suggest an annotation workflow using the state-of-the art methods and highlight their usability, advantages, and limitations for different cases. Finally, we address a particularly difficult problem of deciphering the molecular mechanisms of mutations on proteins and protein complexes to understand the origins and mechanisms of diseases.

Localized structural frustration for evaluating the impact of sequence variants

Population-scale sequencing is increasingly uncovering large numbers of rare single-nucleotide variants (SNVs) in coding regions of the genome. The rarity of these variants makes it challenging to evaluate their deleteriousness with conventional phenotype-genotype associations. Protein structures provide a way of addressing this challenge. Previous efforts have focused on globally quantifying the impact of SNVs on protein stability. However, local perturbations may severely impact protein functionality without strongly disrupting global stability (e.g. in relation to catalysis or allostery). Here, we describe a workflow in which localized frustration, quantifying unfavorable local interactions, is employed as a metric to investigate such effects. Using this workflow on the Protein Databank, we find that frustration produces many immediately intuitive results: for instance, disease-related SNVs create stronger changes in localized frustration than non-disease related variants, and rare SNVs tend to disrupt local interactions to a larger extent than common variants. Less obviously, we observe that somatic SNVs associated with oncogenes and tumor suppressor genes (TSGs) induce very different changes in frustration. In particular, those associated with TSGs change the frustration more in the core than the surface (by introducing loss-of-function events), whereas those associated with oncogenes manifest the opposite pattern, creating gain-of-function events.

A pH-dependent computational approach to the effect of mutations on protein stability

This article describes a novel software implementation for high-throughput scanning mutagenesis with a focus on protein stability. The approach combines molecular mechanics calculations with calculations of protein ionization and a Gaussian-chain model of electrostatic interactions in unfolded state. Comprehensive testing demonstrates a state-of-the-art accuracy for predicted free energy differences on single, double, and triple mutations with a correlation coefficient R above 0.7, which takes about 1.5 min per mutation on a single CPU. Unlike most of existing in silico methods for fast mutagenesis, the stability changes are reported as a continuous function of solution pH for wide pH intervals. We also propose a novel in silico strategy for searching stabilized protein variants that is based on combinatorial scanning mutagenesis using representative amino acid types. Our in silico predictions are in excellent agreement with the hyper-stabilized variants of mesophilic cold shock protein found using the Proside method of direct evolution. © 2016 Wiley Periodicals, Inc.

Electrostatic effects on the folding stability of FKBP12

The roles of electrostatic interactions in protein folding stability have been a matter of debate, largely due to the complexity in the theoretical treatment of these interactions. We have developed computational methods for calculating electrostatic effects on protein folding stability. To rigorously test and further refine these methods, here we carried out experimental studies into electrostatic effects on the folding stability of the human 12-kD FK506 binding protein (FKBP12). This protein has a close homologue, FKBP12.6, with amino acid substitutions in only 18 of their 107 residues. Of the 18 substitutions, 8 involve charged residues. Upon mutating FKBP12 residues at these 8 positions individually into the counterparts in FKBP12.6, the unfolding free energy (ΔGu) of FKBP12 changed by -0.3 to 0.7 kcal/mol. Accumulating stabilizing substitutions resulted in a mutant with a 0.9 kcal/mol increase in stability. Additional charge mutations were grafted from a thermophilic homologue, MtFKBP17, which aligns to FKBP12 with 31% sequence identity over 89 positions. Eleven such charge mutations were studied, with ΔΔGu varying from -2.9 to 0.1 kcal/mol. The predicted electrostatic effects by our computational methods with refinements herein had a root-mean-square deviation of 0.9 kcal/mol from the experimental ΔΔGu values on 16 single mutations of FKBP12. The difference in ΔΔGu between mutations grafted from FKBP12.6 and those from MtFKBP17 suggests that more distant homologues are less able to provide guidance for enhancing folding stability.

MutaBind estimates and interprets the effects of sequence variants on protein-protein interactions

Proteins engage in highly selective interactions with their macromolecular partners. Sequence variants that alter protein binding affinity may cause significant perturbations or complete abolishment of function, potentially leading to diseases. There exists a persistent need to develop a mechanistic understanding of impacts of variants on proteins. To address this need we introduce a new computational method MutaBind to evaluate the effects of sequence variants and disease mutations on protein interactions and calculate the quantitative changes in binding affinity. The MutaBind method uses molecular mechanics force fields, statistical potentials and fast side-chain optimization algorithms. The MutaBind server maps mutations on a structural protein complex, calculates the associated changes in binding affinity, determines the deleterious effect of a mutation, estimates the confidence of this prediction and produces a mutant structural model for download. MutaBind can be applied to a large number of problems, including determination of potential driver mutations in cancer and other diseases, elucidation of the effects of sequence variants on protein fitness in evolution and protein design. MutaBind is available at http://www.ncbi.nlm.nih.gov/projects/mutabind/.

SAAFEC: Predicting the Effect of Single Point Mutations on Protein Folding Free Energy Using a Knowledge-Modified MM/PBSA Approach

Folding free energy is an important biophysical characteristic of proteins that reflects the overall stability of the 3D structure of macromolecules. Changes in the amino acid sequence, naturally occurring or made in vitro, may affect the stability of the corresponding protein and thus could be associated with disease. Several approaches that predict the changes of the folding free energy caused by mutations have been proposed, but there is no method that is clearly superior to the others. The optimal goal is not only to accurately predict the folding free energy changes, but also to characterize the structural changes induced by mutations and the physical nature of the predicted folding free energy changes. Here we report a new method to predict the Single Amino Acid Folding free Energy Changes (SAAFEC) based on a knowledge-modified Molecular Mechanics Poisson-Boltzmann (MM/PBSA) approach. The method is comprised of two main components: a MM/PBSA component and a set of knowledge based terms delivered from a statistical study of the biophysical characteristics of proteins. The predictor utilizes a multiple linear regression model with weighted coefficients of various terms optimized against a set of experimental data. The aforementioned approach yields a correlation coefficient of 0.65 when benchmarked against 983 cases from 42 proteins in the ProTherm database. Availability: the webserver can be accessed via http://compbio.clemson.edu/SAAFEC/.

Investigating the linkage between disease-causing amino acid variants and their effect on protein stability and binding

Single amino acid variations (SAV) occurring in human population result in natural differences between individuals or cause diseases. It is well understood that the molecular effect of SAV can be manifested as changes of the wild type characteristics of the corresponding protein, among which are the protein stability and protein interactions. Typically the effect of SAV on protein stability and interactions was assessed via the changes of the wild type folding and binding free energies. However, in terms of SAV affecting protein functionally and disease susceptibility, one wants to know to what extend the wild type function is perturbed by the SAV. Here it is demonstrated that relative, rather than the absolute, change of the folding and binding free energy serves as a good indicator for SAV association with disease. Using HumVar as a source for disease-causing SAV and experimentally determined free energy changes from ProTherm and SKEMPI databases, correlation coefficients (CC) between the disease index (Pd) and relative folding (Ppr,f) and binding (Ppr,b) probability indexes, respectively, was achieved. The obtained CCs demonstrated the applicability of the proposed approach and it served as good indicator for SAV association with disease.

Applications of Protein Thermodynamic Database for Understanding Protein Mutant Stability and Designing Stable Mutants

Protein stability is the free energy difference between unfolded and folded states of a protein, which lies in the range of 5-25 kcal/mol. Experimentally, protein stability is measured with circular dichroism, differential scanning calorimetry, and fluorescence spectroscopy using thermal and denaturant denaturation methods. These experimental data have been accumulated in the form of a database, ProTherm, thermodynamic database for proteins and mutants. It also contains sequence and structure information of a protein, experimental methods and conditions, and literature information. Different features such as search, display, and sorting options and visualization tools have been incorporated in the database. ProTherm is a valuable resource for understanding/predicting the stability of proteins and it can be accessed at http://www.abren.net/protherm/ . ProTherm has been effectively used to examine the relationship among thermodynamics, structure, and function of proteins. We describe the recent progress on the development of methods for understanding/predicting protein stability, such as (1) general trends on mutational effects on stability, (2) relationship between the stability of protein mutants and amino acid properties, (3) applications of protein three-dimensional structures for predicting their stability upon point mutations, (4) prediction of protein stability upon single mutations from amino acid sequence, and (5) prediction methods for addressing double mutants. A list of online resources for predicting has also been provided.

Integration of structural dynamics and molecular evolution via protein interaction networks: a new era in genomic medicine

Sequencing technologies are revealing many new non-synonymous single nucleotide variants (nsSNVs) in each personal exome. To assess their functional impacts, comparative genomics is frequently employed to predict if they are benign or not. However, evolutionary analysis alone is insufficient, because it misdiagnoses many disease-associated nsSNVs, such as those at positions involved in protein interfaces, and because evolutionary predictions do not provide mechanistic insights into functional change or loss. Structural analyses can aid in overcoming both of these problems by incorporating conformational dynamics and allostery in nSNV diagnosis. Finally, protein-protein interaction networks using systems-level methodologies shed light onto disease etiology and pathogenesis. Bridging these network approaches with structurally resolved protein interactions and dynamics will advance genomic medicine.

Mutations in the KDM5C ARID Domain and Their Plausible Association with Syndromic Claes-Jensen-Type Disease

Mutations in KDM5C gene are linked to X-linked mental retardation, the syndromic Claes-Jensen-type disease. This study focuses on non-synonymous mutations in the KDM5C ARID domain and evaluates the effects of two disease-associated missense mutations (A77T and D87G) and three not-yet-classified missense mutations (R108W, N142S, and R179H). We predict the ARID domain’s folding and binding free energy changes due to mutations, and also study the effects of mutations on protein dynamics. Our computational results indicate that A77T and D87G mutants have minimal effect on the KDM5C ARID domain stability and DNA binding. In parallel, the change in the free energy unfolding caused by the mutants A77T and D87G were experimentally measured by urea-induced unfolding experiments and were shown to be similar to the in silico predictions. The evolutionary conservation analysis shows that the disease-associated mutations are located in a highly-conserved part of the ARID structure (N-terminal domain), indicating their importance for the KDM5C function. N-terminal residues’ high conservation suggests that either the ARID domain utilizes the N-terminal to interact with other KDM5C domains or the N-terminal is involved in some yet unknown function. The analysis indicates that, among the non-classified mutations, R108W is possibly a disease-associated mutation, while N142S and R179H are probably harmless.