SAAMBE-3D: Predicting Effect of Mutations on Protein-Protein Interactions

Robust Prediction of Relative Binding Energies for Protein-Protein Complex Mutations Using Free Energy Perturbation Calculations

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how Protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

A benchmark for evaluation of structure-based online tools for antibody-antigen binding affinity

The prediction of binding affinity changes caused by missense mutations can elucidate antigen-antibody interactions. A few accessible structure-based online computational tools have been proposed. However, selecting suitable software for particular research is challenging, especially research on the SARS-CoV-2 spike protein with antibodies. Therefore, benchmarking of the mutation-diverse SARS-CoV-2 datasets is critical. Here, we collected the datasets including 1216 variants about the changes in binding affinity of antigens from 22 complexes for SARS-CoV-2 S proteins and 22 monoclonal antibodies as well as applied them to evaluate the performance of seven binding affinity prediction tools. The tested tools' Pearson correlations between predicted and measured changes in binding affinity were between -0.158 and 0.657, while accuracy in classification tasks on predicting increasing or decreasing affinity ranged from 0.444 to 0.834. These tools performed relatively better on predicting single mutations, especially at epitope sites, whereas poor performance on extremely decreasing affinity. The tested tools were relatively insensitive to the experimental techniques used to obtain structures of complexes. In summary, we constructed a list of datasets and evaluated a range of structure-based online prediction tools that will explicate relevant processes of antigen-antibody interactions and enhance the computational design of therapeutic monoclonal antibodies.

MuToN Quantifies Binding Affinity Changes upon Protein Mutations by Geometric Deep Learning

Assessing changes in protein-protein binding affinity due to mutations helps understanding a wide range of crucial biological processes within cells. Despite significant efforts to create accurate computational models, predicting how mutations affect affinity remains challenging due to the complexity of the biological mechanisms involved. In the present work, a geometric deep learning framework called MuToN is introduced for quantifying protein binding affinity change upon residue mutations. The method, designed with geometric attention networks, is mechanism-aware. It captures changes in the protein binding interfaces of mutated complexes and assesses the allosteric effects of amino acids. Experimental results highlight MuToN's superiority compared to existing methods. Additionally, MuToN's flexibility and effectiveness are illustrated by its precise predictions of binding affinity changes between SARS-CoV-2 variants and the ACE2 complex.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years 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 25 years, 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 186 VIPs, resulting in a total of 403 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 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

ProBASS - a language model with sequence and structural features for predicting the effect of mutations on binding affinity

Protein-protein interactions (PPIs) govern virtually all cellular processes. Even a single mutation within PPI can significantly influence overall protein functionality and potentially lead to various types of diseases. To date, numerous approaches have emerged for predicting the change in free energy of binding (ΔΔGbind) resulting from mutations, yet the majority of these methods lack precision. In recent years, protein language models (PLMs) have been developed and shown powerful predictive capabilities by leveraging both sequence and structural data from protein-protein complexes. Yet, PLMs have not been optimized specifically for predicting ΔΔGbind. We developed an approach to predict effects of mutations on PPI binding affinity based on two most advanced protein language models ESM2 and ESM-IF1 that incorporate PPI sequence and structural features, respectively. We used the two models to generate embeddings for each PPI mutant and subsequently fine-tuned our model by training on a large dataset of experimental ΔΔGbind values. Our model, ProBASS (Protein Binding Affinity from Structure and Sequence) achieved a correlation with experimental ΔΔGbind values of 0.83 ± 0.05 for single mutations and 0.69 ± 0.04 for double mutations when model training and testing was done on the same PDB. Moreover, ProBASS exhibited very high correlation (0.81 ± 0.02) between prediction and experiment when training and testing was performed on a dataset containing 2325 single mutations in 132 PPIs. ProBASS surpasses the state-of-the-art methods in correlation with experimental data and could be further trained as more experimental data becomes available. Our results demonstrate that the integration of extensive datasets containing ΔΔGbind values across multiple PPIs to refine the pre-trained PLMs represents a successful approach for achieving a precise and broadly applicable model for ΔΔGbind prediction, greatly facilitating future protein engineering and design studies.

DeepPPAPredMut: deep ensemble method for predicting the binding affinity change in protein-protein complexes upon mutation

MOTIVATION

Protein-protein interactions underpin many cellular processes and their disruption due to mutations can lead to diseases. With the evolution of protein structure prediction methods like AlphaFold2 and the availability of extensive experimental affinity data, there is a pressing need for updated computational tools that can efficiently predict changes in binding affinity caused by mutations in protein-protein complexes.

RESULTS

We developed a deep ensemble model that leverages protein sequences, predicted structure-based features, and protein functional classes to accurately predict the change in binding affinity due to mutations. The model achieved a correlation of 0.97 and a mean absolute error (MAE) of 0.35 kcal/mol on the training dataset, and maintained robust performance on the test set with a correlation of 0.72 and a MAE of 0.83 kcal/mol. Further validation using Leave-One-Out Complex (LOOC) cross-validation exhibited a correlation of 0.83 and a MAE of 0.51 kcal/mol, indicating consistent performance.

AVAILABILITY AND IMPLEMENTATION

https://web.iitm.ac.in/bioinfo2/DeepPPAPredMut/index.html.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Predicting the immune escape of SARS-CoV-2 neutralizing antibodies upon mutation

COVID-19 has resulted in millions of deaths and severe impact on economies worldwide. Moreover, the emergence of SARS-CoV-2 variants presented significant challenges in controlling the pandemic, particularly their potential to avoid the immune system and evade vaccine immunity. This has led to a growing need for research to predict how mutations in SARS-CoV-2 reduces the ability of antibodies to neutralize the virus. In this study, we assembled a set of 1813 mutations from the interface of SARS-CoV-2 spike protein's receptor binding domain (RBD) and neutralizing antibody complexes and developed a machine learning model to classify high or low escape mutations using interaction energy, inter-residue contacts and predicted binding free energy change. Our approach achieved an Area under the Receiver Operating Characteristics (ROC) Curve (AUC) of 0.91 using the Random Forest classifier on the test dataset with 217 mutations. The model was further utilized to predict the escape mutations on a dataset of 29,165 mutations located at the interface of 83 RBD-neutralizing antibody complexes. A small subset of this dataset was also validated based on available experimental data. We found that top 10 % high escape mutations were dominated by charged to nonpolar mutations whereas low escape mutations were dominated by polar to nonpolar mutations. We believe that the present method will allow prioritization of high/low escape mutations in the context of neutralizing antibodies targeting SARS-CoV-2 RBD region and assist antibody design for current and emerging variants.

Deciphering principles of nucleosome interactions and impact of cancer-associated mutations from comprehensive interaction network analysis

Nucleosomes represent hubs in chromatin organization and gene regulation and interact with a plethora of chromatin factors through different modes. In addition, alterations in histone proteins such as cancer mutations and post-translational modifications have profound effects on histone/nucleosome interactions. To elucidate the principles of histone interactions and the effects of those alterations, we developed histone interactomes for comprehensive mapping of histone-histone interactions (HHIs), histone-DNA interactions (HDIs), histone-partner interactions (HPIs) and DNA-partner interactions (DPIs) of 37 organisms, which contains a total of 3808 HPIs from 2544 binding proteins and 339 HHIs, 100 HDIs and 142 DPIs across 110 histone variants. With the developed networks, we explored histone interactions at different levels of granularities (protein-, domain- and residue-level) and performed systematic analysis on histone interactions at a large scale. Our analyses have characterized the preferred binding hotspots on both nucleosomal/linker DNA and histone octamer and unraveled diverse binding modes between nucleosome and different classes of binding partners. Last, to understand the impact of histone cancer-associated mutations on histone/nucleosome interactions, we complied one comprehensive cancer mutation dataset including 7940 cancer-associated histone mutations and further mapped those mutations onto 419,125 histone interactions at the residue level. Our quantitative analyses point to histone cancer-associated mutations' strongly disruptive effects on HHIs, HDIs and HPIs. We have further predicted 57 recurrent histone cancer mutations that have large effects on histone/nucleosome interactions and may have driver status in oncogenesis.

VarEPS-Influ:an risk evaluation system of occurred and virtual variations of influenza virus genomes

Influenza viruses undergo frequent genomic mutations, leading to potential cross-species transmission, phenotypic changes, and challenges in diagnostic reagents and vaccines. Accurately evaluating and predicting the risk of such variations remain significant challenges. To address this, we developed the VarEPS-Influ database, an influenza virus variations risk evaluation system (VarEPS-Influ). This database employs a 'multi-dimensional evaluation of mutations' strategy, utilizing various tools to assess the physical and chemical properties, primary, secondary, and tertiary structures, receptor affinity, antibody binding capacity, antigen epitopes, and other aspects of the variation's impact. Additionally, we consider space-time distribution, host species distribution, pedigree analysis, and frequency of mutations to provide a comprehensive risk evaluation of mutations and viruses. The VarEPS-Influ database evaluates both observed variations and virtual variations (variations that have not yet occurred), thereby addressing the time-lag issue in risk predictions. Our current one-stop evaluation system for influenza virus genomic variation integrates 1065290 sequences from 224 927 Influenza A, B and C isolates retrieved from public resources. Researchers can freely access the data at https://nmdc.cn/influvar/.

Software for Predicting Binding Free Energy of Protein-Protein Complexes and Their Mutants

Protein-protein binding affinity prediction is important for understanding complex biochemical pathways and to uncover protein interaction networks. Quantitative estimation of the binding affinity changes caused by mutations can provide critical information for protein function annotation and genetic disease diagnoses. The binding free energies of protein-protein complexes can be predicted using several computational tools. This chapter is a summary of software developed for the prediction of binding free energies for protein-protein complexes and their mutants.

What Strengthens Protein-Protein Interactions: Analysis and Applications of Residue Correlation Networks

Identifying residues critical to protein-protein binding and efficient design of stable and specific protein binders are challenging tasks. Extending beyond the direct contacts in a protein-protein binding interface, our study employs computational modeling to reveal the essential network of residue interactions and dihedral angle correlations critical in protein-protein recognition. We hypothesized that mutating residues exhibiting highly correlated dynamic motion within the interaction network could efficiently optimize protein-protein interactions to create tight and selective protein binders. We tested this hypothesis using the ubiquitin (Ub) and MERS coronaviral papain-like protease (PLpro) complex, since Ub is a central player in multiple cellular functions and PLpro is an antiviral drug target. Our designed ubiquitin variant (UbV) hosting three mutated residues displayed a ∼3,500-fold increase in functional inhibition relative to wild-type Ub. Further optimization of two C-terminal residues within the Ub network resulted in a KD of 1.5 nM and IC50 of 9.7 nM for the five-point Ub mutant, eliciting 27,500-fold and 5,500-fold enhancements in affinity and potency, respectively, as well as improved selectivity, without destabilizing the UbV structure. Our study highlights residue correlation and interaction networks in protein-protein interactions, and introduces an effective approach to design high-affinity protein binders for cell biology research and future therapeutics.

Quantification of biases in predictions of protein-protein binding affinity changes upon mutations

Understanding the impact of mutations on protein-protein binding affinity is a key objective for a wide range of biotechnological applications and for shedding light on disease-causing mutations, which are often located at protein-protein interfaces. Over the past decade, many computational methods using physics-based and/or machine learning approaches have been developed to predict how protein binding affinity changes upon mutations. They all claim to achieve astonishing accuracy on both training and test sets, with performances on standard benchmarks such as SKEMPI 2.0 that seem overly optimistic. Here we benchmarked eight well-known and well-used predictors and identified their biases and dataset dependencies, using not only SKEMPI 2.0 as a test set but also deep mutagenesis data on the severe acute respiratory syndrome coronavirus 2 spike protein in complex with the human angiotensin-converting enzyme 2. We showed that, even though most of the tested methods reach a significant degree of robustness and accuracy, they suffer from limited generalizability properties and struggle to predict unseen mutations. Interestingly, the generalizability problems are more severe for pure machine learning approaches, while physics-based methods are less affected by this issue. Moreover, undesirable prediction biases toward specific mutation properties, the most marked being toward destabilizing mutations, are also observed and should be carefully considered by method developers. We conclude from our analyses that there is room for improvement in the prediction models and suggest ways to check, assess and improve their generalizability and robustness.

Inferring the Effects of Protein Variants on Protein-Protein Interactions with Interpretable Transformer Representations

Identifying pathogenetic variants and inferring their impact on protein-protein interactions sheds light on their functional consequences on diseases. Limited by the availability of experimental data on the consequences of protein interaction, most existing methods focus on building models to predict changes in protein binding affinity. Here, we introduced MIPPI, an end-to-end, interpretable transformer-based deep learning model that learns features directly from sequences by leveraging the interaction data from IMEx. MIPPI was specifically trained to determine the types of variant impact (increasing, decreasing, disrupting, and no effect) on protein-protein interactions. We demonstrate the accuracy of MIPPI and provide interpretation through the analysis of learned attention weights, which exhibit correlations with the amino acids interacting with the variant. Moreover, we showed the practicality of MIPPI in prioritizing de novo mutations associated with complex neurodevelopmental disorders and the potential to determine the pathogenic and driving mutations. Finally, we experimentally validated the functional impact of several variants identified in patients with such disorders. Overall, MIPPI emerges as a versatile, robust, and interpretable model, capable of effectively predicting mutation impacts on protein-protein interactions and facilitating the discovery of clinically actionable variants.

Analysis of RBD-ACE2 interactions in livestock species as a factor in the spread of SARS-CoV-2 among animals

The high mutation rate of SARS-CoV-2, which has led to the emergence of a number of virus variants, creates risks of transmission from humans to animal species and the emergence of new animal reservoirs of COVID-19. This study aimed to identify animal species among livestock susceptible to infection and develop an approach that would be possible to use for assessing the hazards caused by new SARS-CoV-2 variants for animals. Bioinformatic analysis was used to evaluate the ability of receptor-binding domains (RBDs) of different SARS-CoV-2 variants to interact with ACE2 receptors of livestock species. The results indicated that the stability of RBD-ACE2 complexes depends on both amino acid residues in the ACE2 sequences of animal species and on mutations in the RBDs of SARS-CoV-2 variants, with the residues in the interface of the RBD-ACE2 complex being the most important. All studied SARS-CoV-2 variants had high affinity for ferret and American mink receptors, while the affinity for horse, donkey, and bird species' receptors significantly increased in the highly mutated Omicron variant. Hazards that future SARS-CoV-2 variants may acquire specificity to new animal species remain high given the mutability of the virus. The continued use and expansion of the bioinformatic approach presented in this study may be relevant for monitoring transmission risks and preventing the emergence of new reservoirs of COVID-19 among animals.

Predator: Predicting the Impact of Cancer Somatic Mutations on Protein-Protein Interactions

Since many biological processes are governed by protein-protein interactions, understanding which mutations lead to a disruption in these interactions is profoundly important for cancer research. Most of the existing methods focus on the stability of the protein without considering the specific effects of a mutation on its interactions with other proteins. Here, we focus on somatic mutations that appear on the interface regions of the protein and predict the interactions that would be affected by a mutation of interest. We build an ensemble model, Predator, that classifies the interface mutations as disruptive or nondisruptive based on the predicted effects of mutations on specific protein-protein interactions. We show that Predator outperforms existing approaches in literature in terms of prediction accuracy. We then apply Predator on various TCGA cancer cohorts and perform comprehensive analysis at cohort level, patient level, and gene level in determining the genes whose interface mutations tend to yield a disruption in its interactions. The predictions obtained by Predator shed light on interesting patterns on several genes for each cohort regarding their potential as cancer drivers. Our analyses further reveal that the identified genes and their frequently disrupted partners exhibit patterns of mutually exclusivity across cancer cohorts under study.

Insight into de-regulation of amino acid feedback inhibition: a focus on structure analysis method

Regulation of amino acid's biosynthetic pathway is of significant importance to maintain homeostasis and cell functions. Amino acids regulate their biosynthetic pathway by end-product feedback inhibition of enzymes catalyzing committed steps of a pathway. Discovery of new feedback resistant enzyme variants to enhance industrial production of amino acids is a key objective in industrial biotechnology. Deregulation of feedback inhibition has been achieved for various enzymes using in vitro and in silico mutagenesis techniques. As enzyme's function, its substrate binding capacity, catalysis activity, regulation and stability are dependent on its structural characteristics, here, we provide detailed structural analysis of all feedback sensitive enzyme targets in amino acid biosynthetic pathways. Current review summarizes information regarding structural characteristics of various enzyme targets and effect of mutations on their structures and functions especially in terms of deregulation of feedback inhibition. Furthermore, applicability of various experimental as well as computational mutagenesis techniques to accomplish feedback resistance has also been discussed in detail to have an insight into various aspects of research work reported in this particular field of study.

Predicting the Effect of Single Mutations on Protein Stability and Binding with Respect to Types of Mutations

The development of methods and algorithms to predict the effect of mutations on protein stability, protein-protein interaction, and protein-DNA/RNA binding is necessitated by the needs of protein engineering and for understanding the molecular mechanism of disease-causing variants. The vast majority of the leading methods require a database of experimentally measured folding and binding free energy changes for training. These databases are collections of experimental data taken from scientific investigations typically aimed at probing the role of particular residues on the above-mentioned thermodynamic characteristics, i.e., the mutations are not introduced at random and do not necessarily represent mutations originating from single nucleotide variants (SNV). Thus, the reported performance of the leading algorithms assessed on these databases or other limited cases may not be applicable for predicting the effect of SNVs seen in the human population. Indeed, we demonstrate that the SNVs and non-SNVs are not equally presented in the corresponding databases, and the distribution of the free energy changes is not the same. It is shown that the Pearson correlation coefficients (PCCs) of folding and binding free energy changes obtained in cases involving SNVs are smaller than for non-SNVs, indicating that caution should be used in applying them to reveal the effect of human SNVs. Furthermore, it is demonstrated that some methods are sensitive to the chemical nature of the mutations, resulting in PCCs that differ by a factor of four across chemically different mutations. All methods are found to underestimate the energy changes by roughly a factor of 2.

Combining MOE Bioinformatics Analysis and In Vitro Pseudovirus Neutralization Assays to Predict the Neutralizing Ability of CV30 Monoclonal Antibody on SARS-CoV-2 Variants

Combining bioinformatics and in vitro cytology assays, a predictive method was established to quickly evaluate the protective effect of immunity acquired through SARS-CoV-2 infection against variants. Bioinformatics software was first used to predict the changes in the affinity of variant antigens to the CV30 monoclonal antibody by integrating bioinformatics and cytology assays. Then, the ability of the antibody to neutralize the variant antigen was further verified, and the ability of the CV30 to neutralize the new variant strain was predicted through pseudovirus neutralization experiments. The current study has demonstrated that when the Molecular Operating Environment (MOE) predicts |ΔBFE| ≤ 3.0003, it suggests that the CV30 monoclonal antibody exhibits some affinity toward the variant strain and can potentially neutralize it. However, if |ΔBFE| ≥ 4.1539, the CV30 monoclonal antibody does not display any affinity for the variant strain and cannot neutralize it. In contrast, if 3.0003 < |ΔBFE| < 4.1539, it is necessary to conduct a series of neutralization tests promptly with the CV30 monoclonal antibody and the variant pseudovirus to obtain results and supplement the existing method, which is faster than the typical procedures. This approach allows for a rapid assessment of the protective efficacy of natural immunity gained through SARS-CoV-2 infection against variants.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

3DVizSNP: a tool for rapidly visualizing missense mutations identified in high throughput experiments in iCn3D

BACKGROUND

High throughput experiments in cancer and other areas of genomic research identify large numbers of sequence variants that need to be evaluated for phenotypic impact. While many tools exist to score the likely impact of single nucleotide polymorphisms (SNPs) based on sequence alone, the three-dimensional structural environment is essential for understanding the biological impact of a nonsynonymous mutation.

RESULTS

We present a program, 3DVizSNP, that enables the rapid visualization of nonsynonymous missense mutations extracted from a variant caller format file using the web-based iCn3D visualization platform. The program, written in Python, leverages REST APIs and can be run locally without installing any other software or databases, or from a webserver hosted by the National Cancer Institute. It automatically selects the appropriate experimental structure from the Protein Data Bank, if available, or the predicted structure from the AlphaFold database, enabling users to rapidly screen SNPs based on their local structural environment. 3DVizSNP leverages iCn3D annotations and its structural analysis functions to assess changes in structural contacts associated with mutations.

CONCLUSIONS

This tool enables researchers to efficiently make use of 3D structural information to prioritize mutations for further computational and experimental impact assessment. The program is available as a webserver at https://analysistools.cancer.gov/3dvizsnp or as a standalone python program at https://github.com/CBIIT-CGBB/3DVizSNP .

On the linkage of thermodynamics and pathogenicity

This review outlines the effect of disease-causing mutations on proteins' thermodynamics. Two major thermodynamics quantities, which are essential for structural integrity, the folding and binding free energy changes caused by missense mutations, are considered. It is emphasized that disease effects in case of complex diseases may originate from several mutations over several genes, while monogenic diseases are caused by mutation is a single gene. Nevertheless, in both cases it is shown that pathogenic mutations cause larger perturbations of the above-mentioned thermodynamics quantities as compared with the benign mutations. Recent works demonstrating the effect of pathogenic mutations on the above-mentioned thermodynamics quantities, as well as on structural dynamics and allosteric pathways, are reviewed.

An optimized thermodynamics integration protocol for identifying beneficial mutations in antibody design

Accurate identification of beneficial mutations is central to antibody design. Many knowledge-based (KB) computational approaches have been developed to predict beneficial mutations, but their accuracy leaves room for improvement. Thermodynamic integration (TI) is an alchemical free energy algorithm that offers an alternative technique for identifying beneficial mutations, but its performance has not been evaluated. In this study, we developed an efficient TI protocol with high accuracy for predicting binding free energy changes of antibody mutations. The improved TI method outperforms KB methods at identifying both beneficial and deleterious mutations. We observed that KB methods have higher accuracies in predicting deleterious mutations than beneficial mutations. A pipeline using KB methods to efficiently exclude deleterious mutations and TI to accurately identify beneficial mutations was developed for high-throughput mutation scanning. The pipeline was applied to optimize the binding affinity of a broadly sarbecovirus neutralizing antibody 10-40 against the circulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) omicron variant. Three identified beneficial mutations show strong synergy and improve both binding affinity and neutralization potency of antibody 10-40. Molecular dynamics simulation revealed that the three mutations improve the binding affinity of antibody 10-40 through the stabilization of an altered binding mode with increased polar and hydrophobic interactions. Above all, this study presents an accurate and efficient TI-based approach for optimizing antibodies and other biomolecules.

The impact of mutation sets in receptor-binding domain of SARS-CoV-2 variants on the stability of RBD–ACE2 complex

Aim: Bioinformatic analysis of mutation sets in receptor-binding domain (RBD) of currently and previously circulating SARS-CoV-2 variants of concern (VOCs) and interest (VOIs) to assess their ability to bind the ACE2 receptor. Methods: In silico sequence and structure-oriented approaches were used to evaluate the impact of single and multiple mutations. Results: Mutations detected in VOCs and VOIs led to the reduction of binding free energy of the RBD–ACE2 complex, forming additional chemical bonds with ACE2, and to an increase of RBD–ACE2 complex stability. Conclusion: Mutation sets characteristic of SARS-CoV-2 variants have complex effects on the ACE2 receptor-binding affinity associated with amino acid interactions at mutation sites, as well as on the acquisition of other viral adaptive advantages.

Structural basis for the mutation-induced dysfunction of the human IL-15/IL-15α receptor complex

In silico strategies offer a reliable, fast, and inexpensive, way compared to the clumsy in vitro approaches to boost understanding of the effect of amino acid substitution on the structure and consequently the associated function of proteins. In the present work, we report an atomistic-based, reliable in silico structural and energetic framework of the interactions between the receptor-binding domain of the Interleukin-15 (IL-15) protein and its receptor Interleukin-15α (IL-15α), consequently, providing qualitative and quantitative details of the key molecular determinants in ligand/receptor recognition. Molecular dynamics simulations were used to investigate the dynamic behavior of the specific binding between IL-15 and IL-15α followed by estimation of the free energies via molecular mechanics/generalized Born surface area (MM/GBSA). In particular, residues Y26, E46, E53, and E89 of the IL-15 protein receptor-binding domain are identified as main hot spots, shaping and governing the stability of the assembly. These results can be used for the development of neutralizing antibodies and the effective structure-based design of protein-protein interaction inhibitors against the so-called orphan disease, vitiligo.

Predicting mutational function using machine learning

Genetic variations are one of the major causes of phenotypic variations between human individuals. Although beneficial as being the substrate of evolution, germline mutations may cause diseases, including Mendelian diseases and complex diseases such as diabetes and heart diseases. Mutations occurring in somatic cells are a main cause of cancer and likely cause age-related phenotypes and other age-related diseases. Because of the high abundance of genetic variations in the human genome, i.e., millions of germline variations per human subject and thousands of additional somatic mutations per cell, it is technically challenging to experimentally verify the function of every possible mutation and their interactions. Significant progress has been made to solve this problem using computational approaches, especially machine learning (ML). Here, we review the progress and achievements made in recent years in this field of research. We classify the computational models in two ways: one according to their prediction goals including protein structural alterations, gene expression changes, and disease risks, and the other according to their methodologies, including non-machine learning methods, classical machine learning methods, and deep neural network methods. For models in each category, we discuss their architecture, prediction accuracy, and potential limitations. This review provides new insights into the applications and future directions of computational approaches in understanding the role of mutations in aging and disease.

Toward reporting standards for the pathogenicity of variant combinations involved in multilocus/oligogenic diseases

Although standards and guidelines for the interpretation of variants identified in genes that cause Mendelian disorders have been developed, this is not the case for more complex genetic models including variant combinations in multiple genes. During a large curation process conducted on 318 research articles presenting oligogenic variant combinations, we encountered several recurring issues concerning their proper reporting and pathogenicity assessment. These mainly concern the absence of strong evidence that refutes a monogenic model and the lack of a proper genetic and functional assessment of the joint effect of the involved variants. With the increasing accumulation of such cases, it has become essential to develop standards and guidelines on how these oligogenic/multilocus variant combinations should be interpreted, validated, and reported in order to provide high-quality data and supporting evidence to the scientific community.

Genetic landscape in coagulation factor XIII associated defects – Advances in coagulation and beyond

Coagulation factor XIII (FXIII) acts as a fine fulcrum in blood plasma that maintains the balance between bleeding and thrombosis by covalently crosslinking the pre-formed fibrin clot into an insoluble one that is resistant to premature fibrinolysis. In plasma, FXIII circulates as a pro-transglutaminase complex composed of the dimeric catalytic FXIII-A encoded by the F13A1 gene and dimeric carrier/regulatory FXIII-B subunits encoded by the F13B gene. Growing evidence accumulated over decades of exhaustive research shows that not only does FXIII play major roles in both pathological extremes of hemostasis i.e. bleeding and thrombosis, but that it is, in fact, a pleiotropic protein with physiological roles beyond coagulation. However, the current FXIII genetic-epidemiological literature is overwhelmingly derived from the bleeding pathology associated with its deficiency. In this article we review the current clinical, functional, and molecular understanding of this fascinating multifaceted protein, especially putting into the same perspective its genetic landscape.

Hereditary cerebral small vessel disease: Assessment of a HTRA1 variant using protein stability predictors and 3D modelling

Cerebral autosomal recessive arteriopathy with subcortical infarcts and leukoencephalopathy (CARASIL) is an autosomal recessive vascular disorder caused by biallellic variants in HTRA1. Recently, it has been reported that several heterozygous mutations in HTRA1 are responsible for a milder late-onset cerebral small vessel disease (CSVD) with an autosomal dominant pattern of inheritance. The majority of them are missense that affects the Htr1A protease activity due to a dominant-negative effect caused by defective trimerization or monomer activation. The molecular mechanism related to the structural destabilization of the protein supports the practical utility of integrating computational stability predictors to prioritize candidate variants in this gene. In this work, we report a family with several members diagnosed with subcortical ischemic events and progressive cognitive impairment caused by the novel c.820C > G, p.(Arg274Gly) heterozygous variant in HTRA1 segregating in an autosomal dominant manner and propose its molecular mechanism by a three-dimensional model of the protein's structure.

A comprehensive in silico exploration of the impacts of missense variants on two different conformations of human pirin protein

BACKGROUND

Pirin, a member of the cupin superfamily, is an iron-binding non-heme protein. It acts as a coregulator of several transcription factors, especially the members of NFκB transcription factor family. Based on the redox state of its iron cofactor, it can assume two different conformations and thereby act as a redox sensor inside the nucleus. Previous studies suggested that pirin may be associated with cancer, inflammatory diseases as well as COVID-19 severities. Hence, it is important to explore the pathogenicity of its missense variants. In this study, we used a number of in silico tools to investigate the effects of missense variants of pirin on its structure, stability, metal cofactor binding affinity and interactions with partner proteins. In addition, we used protein dynamics simulation to elucidate the effects of selected variants on its dynamics. Furthermore, we calculated the frequencies of haplotypes containing pirin missense variants across five major super-populations (African, Admixed American, East Asian, European and South Asian).

RESULTS

Among a total of 153 missense variants of pirin, 45 were uniformly predicted to be pathogenic. Of these, seven variants can be considered for further experimental studies. Variants R59P and L116P were predicted to significantly destabilize and damage pirin structure, substantially reduce its affinity to its binding partners and alter pirin residue fluctuation profile via changing the flexibility of several key residues. Additionally, variants R59Q, F78V, G98D, V151D and L220P were found to impact pirin structure and function in multiple ways. As no haplotype was identified to be harboring more than one missense variant, further interrogation of the individual effects of these seven missense variants is highly recommended.

CONCLUSIONS

Pirin is involved in the transcriptional regulation of several genes and can play an important role in inflammatory responses. The variants predicted to be pathogenic in this study may thus contribute to a better understanding of the underlying molecular mechanisms of various inflammatory diseases. Future studies should be focused on clarifying if any of these variants can be used as disease biomarkers.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1186/s42269-022-00917-7.

Comparative molecular dynamics study of the receptor-binding domains in SARS-CoV-2 and SARS-CoV and the effects of mutations on the binding affinity

Here, we report on a computational comparison of the receptor-binding domains (RBDs) on the spike proteins of severe respiratory syndrome coronavirus-2 (SARS-CoV-2) and SARS-CoV in free forms and as complexes with angiotensin-converting enzyme 2 (ACE2) as their receptor in humans. The impact of 42 mutations discovered so far on the structure and thermodynamics of SARS-CoV-2 RBD was also assessed. The binding affinity of SARS-CoV-2 RBD for ACE2 is higher than that of SARS-CoV RBD. The binding of COVA2-04 antibody to SARS-CoV-2 RBD is more energetically favorable than the binding of COVA2-39, but also less favorable than the formation of SARS-CoV-2 RBD-ACE2 complex. The net charge, the dipole moment and hydrophilicity of SARS-CoV-2 RBD are higher than those of SARS-CoV RBD, producing lower solvation and surface free energies and thus lower stability. The structure of SARS-CoV-2 RBD is also more flexible and more open, with a larger solvent-accessible surface area than that of SARS-CoV RBD. Single-point mutations have a dramatic effect on distribution of charges, most prominently at the site of substitution and its immediate vicinity. These charge alterations alter the free energy landscape, while X→F mutations exhibit a stabilizing effect on the RBD structure through π stacking. F456 and W436 emerge as two key residues governing the stability and affinity of the spike protein for its ACE2 receptor. These analyses of the structural differences and the impact of mutations on different viral strains and members of the coronavirus genera are an essential aid in the development of effective therapeutic strategies. Communicated by Ramaswamy H. Sarma.

Substituted Syndecan-2-Derived Mimetic Peptides Show Improved Antitumor Activity over the Parent Syndecan-2-Derived Peptide

We previously showed that a synthetic peptide (S2-P) corresponding to a portion of the human syndecan-2 (SDC2) sequence can bind to the pro-domain of matrix metalloproteinase-7 (MMP-7) to inhibit colon cancer activities. Since S2-P had a relatively weak binding affinity for the MMP-7 pro-domain, we herein modified the amino acid sequence of S2-P to improve the anticancer potential. On the basis of the interaction structure of S2-P and MMP-7, four peptides were generated by replacing amino acids near Tyr 51, which is critical for the interaction. The SDC2-mimetic peptides harboring an Ala-to-Asp substitution at the C-terminal side of Tyr 51 (S2-D) or with an Ala-to-Phe substitution at the N-terminal side of Tyr 51 and an Ala-to-Asp substitution at the C-terminal side of Tyr 51 (S2-FE) showed improved interaction affinities for the MMP-7 pro-domain. Compared to S2-P, S2-FE was better able to inhibit the SDC2-MMP-7 interaction, the cell surface localization of MMP-7, the gelatin degradation activity of MMP-7, and the cancer activities (cell migration, invasion, and colony-forming activity) of human HCT116 colon cancer cells in vitro. In vivo, S2-FE inhibited the primary tumor growth and lung metastasis of CT26 mouse colon cancer cells in a xenograft mouse model. Together, these data suggest that S2-FE could be useful therapeutic anticancer peptides for colon cancer.

PLXNB1 mutations in the etiology of idiopathic hypogonadotropic hypogonadism

Idiopathic hypogonadotropic hypogonadism (IHH) comprises a group of rare genetic disorders characterized by pubertal failure caused by gonadotropin-releasing hormone (GnRH) deficiency. Genetic factors involved in semaphorin/plexin signaling have been identified in patients with IHH. PlexinB1, a member of the plexin family receptors, serves as the receptor for semaphorin 4D (Sema4D). In mice, perturbations in Sema4D/PlexinB1 signaling leads to improper GnRH development, highlighting the importance of investigating PlexinB1 mutations in IHH families. In total, 336 IHH patients (normosmic IHH, n = 293 and Kallmann syndrome, n = 43) from 290 independent families were included in the present study. Six PLXNB1 rare sequence variants (p.N361S, p.V608A, p.R636C, p.V672A, p.R1031H, and p.C1318R) are described in eight normosmic IHH patients from seven independent families. These variants were examined using bioinformatic modeling and compared to mutants reported in PLXNA1. Based on these analyses, the variant p.R1031H was assayed for alterations in cell morphology, PlexinB1 expression, and migration using a GnRH cell line and Boyden chambers. Experiments showed reduced membrane expression and impaired migration in cells expressing this variant compared to the wild-type. Our results provide clinical, genetic, molecular/cellular, and modeling evidence to implicate variants in PLXNB1 in the etiology of IHH.

Dynamic disequilibrium-based pathogenicity model in mutated pyrin’s B30.2 domain—Casp1/p20 complex

Background

The B30.2 variants lead to most relevant severity forms of familial Mediterranean fever (FMF) manifestations. The B30.2 domain plays a key role in protein-protein interaction (PPI) of pyrin with other apoptosis proteins and in regulation the cascade of inflammatory reactions. Pyrin-casp1 interaction is mainly responsible for the dysregulation of the inflammatory responses in FMF. Lower binding affinity was observed between the mutant B30.2 pyrin and casp1 without the release of the complete pathogenicity mechanism. The aim of this study was to identify the possible effects of the interface pocked residues in B30.2/SPRY-Casp1/p20 complex using molecular mechanics simulation and in silico analysis.

Results

It was found that Lys671Met, Ser703Ile, and Ala744Ser variants led mainly to shift of the binding affinity (∆G), dissociation constant (K_d), and root mean square deviation (RMSD) in B30.2/SPRY-Casp1/p20 complex leading to dynamic disequilibrium of the p20-B30.2/SPRY complex toward its complex form. The current pathogenicity model and its predicted implementation in the relevant colchicine dosage were delineated.

Conclusion

The molecular mechanics analysis of B30.2/SPRY-p20 complex harboring Lys671Met, Ser703Ile, and Ala744Ser variants showed dynamic disequilibrium of B30.2/SPRY-casp1/p20complex in context of the studied variants that could be a new computational model for FMF pathogenicity. This study also highlighted the specific biochemical markers that could be useful to adjust the colchicine dose in FMF patients.

Missense variants in the TNFA epitopes and their effects on interaction with therapeutic antibodies-in silico analysis

BACKGROUND

Tumor necrosis factor alpha (TNFA) is an important cytokine that influences multiple biological processes. It is one of the key mediators of acute and chronic systemic inflammatory reactions and plays a central role in several autoimmune diseases. A number of approved monoclonal antibodies (mAbs) are widely used to subside these autoimmune diseases. However, there is a high rate of non-responsiveness to treatments with these mAbs. Therefore, it is important to be able to predict responses of the patients in an individualistic manner to these therapeutic antibodies before administration. In the present study, we used in silico tools to explore the effects of missense variants in the respective epitopes of four therapeutic anti-TNFA mAbs-adalimumab (ADA), certolizumab pegol (CZP), golimumab (GLM), and infliximab (IFX)-on their interactions with TNFA.

RESULTS

The binding affinities of CZP and ADA to corresponding epitopes appear to be reduced by four (TNFAR131Q, TNFAE135G, TNFAR138Q, and TNFAR138W) and two (TNFAG66C and TNFAG66S) variants, respectively. The binding of GLM and IFX appears to be affected by TNFAR141S and TNFAR138W, respectively. TNFAG66C and TNFAG66S may be associated with autoimmune diseases, whereas TNFAE135G, TNFAR138W, and TNFAR141S may be pathogenic per se.

CONCLUSION

These variants may contribute to the observed inter-individual variability in response to anti-TNFA mAbs treatments and be used as markers to predict responses, and thus optimize therapeutic benefits to the patients.

Agonist antibody discovery: Experimental, computational, and rational engineering approaches

Agonist antibodies that activate cellular signaling have emerged as promising therapeutics for treating myriad pathologies. Unfortunately, the discovery of rare antibodies with the desired agonist functions is a major bottleneck during drug development. Nevertheless, there has been important recent progress in discovering and optimizing agonist antibodies against a variety of therapeutic targets that are activated by diverse signaling mechanisms. Herein, we review emerging high-throughput experimental and computational methods for agonist antibody discovery as well as rational molecular engineering methods for optimizing their agonist activity.

Structure-function relationships of the disease-linked A218T oxytocin receptor variant

Various single nucleotide polymorphisms (SNPs) in the oxytocin receptor (OXTR) gene have been associated with behavioral traits, autism spectrum disorder (ASD) and other diseases. The non-synonymous SNP rs4686302 results in the OXTR variant A218T and has been linked to core characteristics of ASD, trait empathy and preterm birth. However, the molecular and intracellular mechanisms underlying those associations are still elusive. Here, we uncovered the molecular and intracellular consequences of this mutation that may affect the psychological or behavioral outcome of oxytocin (OXT)-treatment regimens in clinical studies, and provide a mechanistic explanation for an altered receptor function. We created two monoclonal HEK293 cell lines, stably expressing either the wild-type or A218T OXTR. We detected an increased OXTR protein stability, accompanied by a shift in Ca²⁺ dynamics and reduced MAPK pathway activation in the A218T cells. Combined whole-genome and RNA sequencing analyses in OXT-treated cells revealed 7823 differentially regulated genes in A218T compared to wild-type cells, including 429 genes being associated with ASD. Furthermore, computational modeling provided a molecular basis for the observed change in OXTR stability suggesting that the OXTR mutation affects downstream events by altering receptor activation and signaling, in agreement with our in vitro results. In summary, our study provides the cellular mechanism that links the OXTR rs4686302 SNP with genetic dysregulations associated with aspects of ASD.

Cluster learning-assisted directed evolution

Directed evolution, a strategy for protein engineering, optimizes protein properties (i.e., fitness) by expensive and time-consuming screening or selection of large mutational sequence space. Machine learning-assisted directed evolution (MLDE), which screens sequence properties in silico, can accelerate the optimization and reduce the experimental burden. This work introduces a MLDE framework, cluster learning-assisted directed evolution (CLADE), that combines hierarchical unsupervised clustering sampling and supervised learning to guide protein engineering. The clustering sampling selectively picks and screens variants in targeted subspaces, which guides the subsequent generation of diverse training sets. In the last stage, accurate predictions via supervised learning models improve final outcomes. By sequentially screening 480 sequences out of 160,000 in a four-site combinatorial library with five equal experimental batches, CLADE achieves the global maximal fitness hit rate up to 91.0% and 34.0% for GB1 and PhoQ datasets, respectively, improved from 18.6% and 7.2% obtained by random-sampling-based MLDE.

PremPLI: a machine learning model for predicting the effects of missense mutations on protein-ligand interactions

Resistance to small-molecule drugs is the main cause of the failure of therapeutic drugs in clinical practice. Missense mutations altering the binding of ligands to proteins are one of the critical mechanisms that result in genetic disease and drug resistance. Computational methods have made a lot of progress for predicting binding affinity changes and identifying resistance mutations, but their prediction accuracy and speed are still not satisfied and need to be further improved. To address these issues, we introduce a structure-based machine learning method for quantitatively estimating the effects of single mutations on ligand binding affinity changes (named as PremPLI). A comprehensive comparison of the predictive performance of PremPLI with other available methods on two benchmark datasets confirms that our approach performs robustly and presents similar or even higher predictive accuracy than the approaches relying on first-principle statistical mechanics and mixed physics- and knowledge-based potentials while requires much less computational resources. PremPLI can be used for guiding the design of ligand-binding proteins, identifying and understanding disease driver mutations, and finding potential resistance mutations for different drugs. PremPLI is freely available at https://lilab.jysw.suda.edu.cn/research/PremPLI/ and allows to do large-scale mutational scanning.

Sun et al. present PremPLI, a machine learning approach and web tool to predict how missense mutations in a drug’s target protein will affect the drug’s binding affinity. PremPLI can be applied to identify drug resistance mechanisms in cancer cells and microorganisms and develop drugs to combat drug resistance.

Identification of Hotspot Residues in Binding of SARS-CoV-2 Spike and Human ACE2 Proteins

Coronaviruses are a large family of viruses that can cause respiratory infections with varying severity from common cold to severe diseases such as novel coronavirus disease (COVID-19). COVID-19 has been declared as a global pandemic by the World Health Organization on March 11, 2020 and with the development of vaccines it slowed down as of this date. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses its spike glycoprotein (Sgp) to bind human angiotensin-converting enzyme 2 (hACE2) receptor, and mediates membrane fusion and virus entry. The recognition of Sgp to human ACE2 and its high affinity for it has been of great importance since this provides the first step in viral entry to human cells. Therefore, it is crucial to identify key residues (hotspots) in this process. In this study, computational alanine scanning has been performed for Sgp and hACE2. The residues identified with significance in binding and other residues in close proximity were studied further through molecular mechanics-based protein binding free energy change prediction methods. Moreover, the interfacial residues in both proteins were investigated for their cooperative binding. Additionally, folding free energy changes upon mutation to Ala were calculated to assess their effect on stability of Sgp and hACE2. Our results taken together with findings from previous studies revealed the residues that are most significant and are relatively significant in binding of Sgp to human ACE2 protein.

Computational prediction of the effect of amino acid changes on the binding affinity between SARS-CoV-2 spike RBD and human ACE2

The association of the receptor binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein with human angiotensin-converting enzyme 2 (hACE2) represents the first required step for cellular entry. SARS-CoV-2 has continued to evolve with the emergence of several novel variants, and amino acid changes in the RBD have been implicated with increased fitness and potential for immune evasion. Reliably predicting the effect of amino acid changes on the ability of the RBD to interact more strongly with the hACE2 can help assess the implications for public health and the potential for spillover and adaptation into other animals. Here, we introduce a two-step framework that first relies on 48 independent 4-ns molecular dynamics (MD) trajectories of RBD-hACE2 variants to collect binding energy terms decomposed into Coulombic, covalent, van der Waals, lipophilic, generalized Born solvation, hydrogen bonding, π-π packing, and self-contact correction terms. The second step implements a neural network to classify and quantitatively predict binding affinity changes using the decomposed energy terms as descriptors. The computational base achieves a validation accuracy of 82.8% for classifying single-amino acid substitution variants of the RBD as worsening or improving binding affinity for hACE2 and a correlation coefficient of 0.73 between predicted and experimentally calculated changes in binding affinities. Both metrics are calculated using a fivefold cross-validation test. Our method thus sets up a framework for screening binding affinity changes caused by unknown single- and multiple-amino acid changes offering a valuable tool to predict host adaptation of SARS-CoV-2 variants toward tighter hACE2 binding.

Molecular insights into the role of genetic determinants of congenital hypothyroidism

In our previous studies, we have demonstrated the association of certain variants of the thyroid-stimulating hormone receptor (TSHR), thyroid peroxidase (TPO), and thyroglobulin (TG) genes with congenital hypothyroidism. Herein, we explored the mechanistic basis for this association using different in silico tools. The mRNA 3'-untranslated region (3'-UTR) plays key roles in gene expression at the post-transcriptional level. In TSHR variants (rs2268477, rs7144481, and rs17630128), the binding affinity of microRNAs (miRs) (hsa-miR-154-5p, hsa-miR-376a-2-5p, hsa-miR-3935, hsa-miR-4280, and hsa-miR-6858-3p) to the 3'-UTR is disrupted, affecting post-transcriptional gene regulation. TPO and TG are the two key proteins necessary for the biosynthesis of thyroid hormones in the presence of iodide and H₂O₂. Reduced stability of these proteins leads to aberrant biosynthesis of thyroid hormones. Compared to the wild-type TPO protein, the p.S398T variant was found to exhibit less stability and significant rearrangements of intra-atomic bonds affecting the stoichiometry and substrate binding (binding energies, ΔG of wild-type vs. mutant: ‒15 vs. ‒13.8 kcal/mol; and dissociation constant, K_d of wild-type vs. mutant: 7.2E^-12 vs. 7.0E^-11 M). The missense mutations p.G653D and p.R1999W on the TG protein showed altered ΔG (0.24 kcal/mol and 0.79 kcal/mol, respectively). In conclusion, an in silico analysis of TSHR genetic variants in the 3'-UTR showed that they alter the binding affinities of different miRs. The TPO protein structure and mutant protein complex (p.S398T) are less stable, with potentially deleterious effects. A structural and energy analysis showed that TG mutations (p.G653D and p.R1999W) reduce the stability of the TG protein and affect its structure-functional relationship.

VarEPS: an evaluation and prewarning system of known and virtual variations of SARS-CoV-2 genomes

The genomic variations of SARS-CoV-2 continue to emerge and spread worldwide. Some mutant strains show increased transmissibility and virulence, which may cause reduced protection provided by vaccines. Thus, it is necessary to continuously monitor and analyze the genomic variations of SARS-COV-2 genomes. We established an evaluation and prewarning system, SARS-CoV-2 variations evaluation and prewarning system (VarEPS), including known and virtual mutations of SARS-CoV-2 genomes to achieve rapid evaluation of the risks posed by mutant strains. From the perspective of genomics and structural biology, the database comprehensively analyzes the effects of known variations and virtual variations on physicochemical properties, translation efficiency, secondary structure, and binding capacity of ACE2 and neutralizing antibodies. An AI-based algorithm was used to verify the effectiveness of these genomics and structural biology characteristic quantities for risk prediction. This classifier could be further used to group viral strains by their transmissibility and affinity to neutralizing antibodies. This unique resource makes it possible to quickly evaluate the variation risks of key sites, and guide the research and development of vaccines and drugs. The database is freely accessible at www.nmdc.cn/ncovn.

Opioid Addiction and Opioid Receptor Dimerization: Structural Modeling of the OPRD1 and OPRM1 Heterodimer and Its Signaling Pathways

Opioid addiction is a complex phenomenon with genetic, social, and other components. Due to such complexity, it is difficult to interpret the outcome of clinical studies, and thus, mutations found in individuals with these addictions are still not indisputably classified as opioid addiction-causing variants. Here, we computationally investigated two such mutations, A6V and N40D, found in the mu opioid receptor gene OPRM1. The mutations are located in the extracellular domain of the corresponding protein, which is important to the hetero-dimerization of OPRM1 with the delta opioid receptor protein (OPRD1). The hetero-dimerization of OPRD1-OPRM1 affects the signaling pathways activated by opioids and natural peptides and, thus, could be considered a factor contributing to addiction. In this study, we built four 3D structures of molecular pathways, including the G-protein signaling pathway and the β-arrestin signaling pathway of the heterodimer of OPRD1-OPRM1. We also analyzed the effect of mutations of A6V and N40D on the stability of individual OPRM1/OPRD1 molecules and the OPRD1-OPRM1 heterodimer with the goal of inferring their plausible linkage with opioid addiction. It was found that both mutations slightly destabilize OPRM1/OPRD1 monomers and weaken their association. Since hetero-dimerization is a key step for signaling processes, it is anticipated that both mutations may be causing increased addiction risk.

In-silico identification of frequently mutated genes and their co-enriched metabolic pathways associated with Prostate cancer progression

Prostate cancer (PCa) has emerged as a significant health burden in men globally. Several genetic anomalies such as mutations and also epigenetic aberrations are responsible for the heterogeneity of this disease. This study identified the 20 most frequently mutated genes reported in PCa based on literature and database survey. Further gene ontology and functional enrichment analysis were conducted to determine their co-modulated molecular and biological pathways. A protein-protein interaction network was used for the identification of hub genes. These hub genes identified were then subjected to survival analysis. The prognostic values of these identified genes were investigated using GEPIA and HPA. Gene Ontology analysis of the identified genes depicted that these genes significantly contributed to the cell cycle, apoptosis, angiogenesis and TGF-β receptor signalling. Further, the research showed that high expressions of identified mutated genes led to a reduction in the long-term survival of PCa patients, which was supported by immunohistochemical and mRNA expression level data. Our results suggest that identified panel of mutated genes viz., CTNNB1, TP53, ATM, AR and KMT2D play crucial roles in the onset and progression of PCa, thereby providing candidate diagnostic markers for PCa for individualised treatment in the future.

SAMPDI-3D: predicting the effects of protein and DNA mutations on protein-DNA interactions

MOTIVATION

Mutations that alter protein-DNA interactions may be pathogenic and cause diseases. Therefore, it is extremely important to quantify the effect of mutations on protein-DNA binding free energy to reveal the molecular origin of diseases and to assist the development of treatments. Although several methods that predict the change of protein-DNA binding affinity upon mutations in the binding protein were developed, the effect of DNA mutations was not considered yet.

RESULTS

Here, we report a new version of SAMPDI, the SAMPDI-3D, which is a gradient boosting decision tree machine learning method to predict the change of the protein-DNA binding free energy caused by mutations in both the binding protein and the bases of the corresponding DNA. The method is shown to achieve Pearson correlation coefficient of 0.76 and 0.80 in a benchmarking test against experimentally determined change of the binding free energy caused by mutations in the binding protein or DNA, respectively. Furthermore, three datasets collected from literature were used to do blind benchmark for SAMPDI-3D and it is shown that it outperforms all existing state-of-the-art methods. The method is very fast allowing for genome-scale investigations.

AVAILABILITY

It is available as a web server and a stand-code at http://compbio.clemson.edu/SAMPDI-3D/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Non-RBM Mutations Impaired SARS-CoV-2 Spike Protein Regulated to the ACE2 Receptor Based on Molecular Dynamic Simulation

The emergence of novel coronavirus mutants is a main factor behind the deterioration of the epidemic situation. Further studies into the pathogenicity of these mutants are thus urgently needed. Binding of the spinous protein receptor binding domain (RBD) of SARS-CoV-2 to the angiotensin-converting enzyme 2 (ACE2) receptor was shown to initiate coronavirus entry into host cells and lead to their infection. The receptor-binding motif (RBM, 438-506) is a region that directly interacts with ACE2 receptor in the RBD and plays a crucial role in determining affinity. To unravel how mutations in the non-RBM regions impact the interaction between RBD and ACE2, we selected three non-RBM mutant systems (N354D, D364Y, and V367F) from the documented clinical cases, and the Q498A mutant system located in the RBM region served as the control. Molecular dynamics simulation was conducted on the mutant systems and the wild-type (WT) system, and verified experiments also performed. Non-RBM mutations have been shown not only to change conformation of the RBM region but also to significantly influence its hydrogen bonding and hydrophobic interactions. In particular, the D364Y and V367F systems showed a higher affinity for ACE2 owing to their electrostatic interactions and polar solvation energy changes. In addition, although the binding free energy at this point increased after the mutation of N354D, the conformation of the random coil (Pro384-Asp389) was looser than that of other systems, and the combined effect weakened the binding free energy between RBD and ACE2. Interestingly, we also found a random coil (Ala475-Gly485). This random coil is very sensitive to mutations, and both types of mutations increase the binding free energy of residues in this region. We found that the binding loop (Tyr495-Tyr505) in the RBD domain strongly binds to Lys353, an important residue of the ACE2 domain previously identified. The binding free energy of the non-RBM mutant group at the binding loop had positive and negative changes, and these changes were more obvious than that of the Q498A system. The results of this study elucidate the effect of non-RBM mutation on ACE2-RBD binding, and provide new insights for SARS-CoV-2 mutation research.

Application of an integrated computational antibody engineering platform to design SARS-CoV-2 neutralizers

As the COVID-19 pandemic continues to spread, hundreds of new initiatives including studies on existing medicines are running to fight the disease. To deliver a potentially immediate and lasting treatment to current and emerging SARS-CoV-2 variants, new collaborations and ways of sharing are required to create as many paths forward as possible. Here, we leverage our expertise in computational antibody engineering to rationally design/engineer three previously reported SARS-CoV neutralizing antibodies and share our proposal towards anti-SARS-CoV-2 biologics therapeutics. SARS-CoV neutralizing antibodies, m396, 80R and CR-3022 were chosen as templates due to their diversified epitopes and confirmed neutralization potency against SARS-CoV (but not SARS-CoV-2 except for CR3022). Structures of variable fragment (Fv) in complex with receptor binding domain (RBD) from SARS-CoV or SARS-CoV-2 were subjected to our established in silico antibody engineering platform to improve their binding affinity to SARS-CoV-2 and developability profiles. The selected top mutations were ensembled into a focused library for each antibody for further screening. In addition, we convert the selected binders with different epitopes into the trispecific format, aiming to increase potency and to prevent mutational escape. Lastly, to avoid antibody-induced virus activation or enhancement, we suggest application of NNAS and DQ mutations to the Fc region to eliminate effector functions and extend half-life.

Rational optimization of a human neutralizing antibody of SARS-CoV-2

SARS-CoV-2 has caused a worldwide epidemic of coronavirus disease 19 (COVID-19). Antibody drugs present an effective weapon for tens of millions of COVID-19 patients. Antibodies disrupting the interactions between the receptor-binding domain (RBD) of SARS-CoV-2 S protein and the angiotensin converting enzyme 2 (ACE2) effectively block SARS-CoV-2 cell entry into host cells. In order to rapidly develop more potent neutralizing antibodies, we utilized virtual scanning mutageneses and molecular dynamics simulations to optimize the antibody of P2B-2F6 isolated from single B cells of SARS-CoV-2 infected patients. Two potent P2B-2F6 mutants, namely H:V106R and H:V106R/H:P107Y, were found to possess higher binding affinities with the RBD domain of SARS-CoV-2 than others. Polar interactions are preferred near 106 and 107 paratope residues of the heavy chain. The mutations also increase the hydrogen-bonding network formed between the antibody and the RBD. Notably, the optimized antibodies possess potential neutralizing activity against the alarming SARS-CoV-2 variant of N501Y. This study provides insights into structure-based optimization of antibodies with higher affinity to the antigen. We hope that our proposed antibody mutants could contribute to the development of improved therapies against COVID-19.