Change In Binding Affinity Research Articles

Designing of nanobodies against Dengue virus Capsid: a computational affinity maturation approach

Dengue virus, an arbovirus, is one of the most prevalent diseases in the tropical environment and leads to huge number of casualties every year. No therapeutics are available till date against the viral disease and the only medications provide symptomatic relief. In this study, we have focused on utilizing conventional nanobodies and repurposing them for Dengue. Computationally affinity matured, best binding nanobodies tagged with constant antibody regions, could be proposed as therapeutics. These could also be applied for drug delivery purposes due to their high specificity against the viral Capsid. Another application of these nanobodies has been thought to utilize them for diagnostic purposes, to use the nanobodies for viral detection from patient samples at the earliest stage using ELISA. This study may open a new avenue for immunologic study in foreseeable future with the usage of the same molecules for multiple purposes. Highlights Natural nanobodies against viruses were modified for use against Dengue virus Capsid conserved regions. Computational affinity maturation was performed making use of change in binding affinities upon mutating various residues in the complementary determining regions. Docking studies performed to inspect the docking groove, interface analysis and energy calculations. MM/GBSA calculations done to calculate binding free energy of the complex to determine stability of the complex. Communicated by Ramaswamy H. Sarma

Engineering an Enhanced EGFR Engager: Humanization of Cetuximab for Improved Developability.

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase whose proliferative effects can contribute to the development of many types of solid tumors when overexpressed. For this reason, EGFR inhibitors such as cetuximab can play an important role in treating cancers such as colorectal cancer and head and neck cancer. Cetuximab is a chimeric monoclonal antibody containing mouse variable regions that bind to EGFR and prevent it from signaling. Although cetuximab has been used clinically since 2004 to successfully control solid tumors, advances in protein engineering have created the opportunity to address some of its shortcomings. In particular, the presence of mouse sequences could contribute to immunogenicity in the form of anti-cetuximab antibodies, and an occupied glycosylation site in FR3 can contribute to hypersensitivity reactions and product heterogeneity. Using simple framework graft or sequence-/structure-guided approaches, cetuximab was humanized onto 11 new frameworks. In addition to increasing humanness and removing the VH glycosylation site, dynamic light scattering revealed increases in stability, and bio-layer interferometry confirmed minimal changes in binding affinity, with patterns emerging across the humanization method. This work demonstrates the potential to improve the biophysical and clinical properties of first-generation protein therapeutics and highlights the advantages of computationally guided engineering.

Evaluation of protein descriptors in computer-aided rational protein engineering tasks and its application in property prediction in SARS-CoV-2 spike glycoprotein

The importance of protein engineering in the research and development of biopharmaceuticals and biomaterials has increased. Machine learning in computer-aided protein engineering can markedly reduce the experimental effort in identifying optimal sequences that satisfy the desired properties from a large number of possible protein sequences. To develop general protein descriptors for computer-aided protein engineering tasks, we devised new protein descriptors, one sequence-based descriptor (PCgrades), and three structure-based descriptors (PCspairs, 3D-SPIEs_5.4 Å, and 3D-SPIEs_8Å). While the PCgrades and PCspairs include general and statistical information in physicochemical properties in single and pairwise amino acids respectively, the 3D-SPIEs include specific and quantum–mechanical information with parameterized quantum mechanical calculations (FMO2-DFTB3/D/PCM). To evaluate the protein descriptors, we made prediction models with the new descriptors and previously developed descriptors for diverse protein datasets including protein expression and binding affinity change in SARS-CoV-2 spike glycoprotein. As a result, the newly devised descriptors showed a good performance in diverse datasets, in which the PCspairs showed the best performance (R2=0.783 for protein expression and R2=0.711 for binding affinity). As a result, the newly devised descriptors showed a good performance in diverse datasets, in which the PCspairs showed the best performance. Similar approaches with those descriptors would be promising and useful if the prediction models are trained with sufficient quantitative experimental data from high-throughput assays for industrial enzymes or protein drugs.

Structural exploration with AlphaFold2-generated STAT3\u03b1 structure reveals selective elements in STAT3\u03b1-GRIM-19 interactions involved in negative regulation

STAT3, an important transcription factor constitutively activated in cancers, is bound specifically by GRIM-19 and this interaction inhibits STAT3-dependent gene expression. GRIM-19 is therefore, considered as an inhibitor of STAT3 and may be an effective anti-cancer therapeutic target. While STAT3 exists in a dimeric form in the cytoplasm and nucleus, it is mostly present in a monomeric form in the mitochondria. Although GRIM-19-binding domains of STAT3 have been identified in independent experiments, yet the identified domains are not the same, and hence, discrepancies exist. Human STAT3-GRIM-19 complex has not been crystallised yet. Dictated by fundamental biophysical principles, the binding region, interactions and effects of hotspot mutations can provide us a clue to the negative regulatory mechanisms of GRIM-19. Prompted by the very nature of STAT3 being a challenging molecule, and to understand the structural basis of binding and interactions in STAT3α-GRIM-19 complex, we performed homology modelling and ab-initio modelling with evolutionary information using I-TASSER and avant-garde AlphaFold2, respectively, to generate monomeric, and subsequently, dimeric STAT3α structures. The dimeric form of STAT3α structure was observed to potentially exist in an anti-parallel orientation of monomers. We demonstrate that during the interactions with both unphosphorylated and phosphorylated STAT3α, the NTD of GRIM-19 binds most strongly to the NTD of STAT3α, in direct contrast to the earlier works. Key arginine residues at positions 57, 58 and 68 of GRIM-19 are mainly involved in the hydrogen-bonded interactions. An intriguing feature of these arginine residues is that these display a consistent interaction pattern across unphosphorylated and phosphorylated monomers as well as unphosphorylated dimers in STAT3α-GRIM-19 complexes. MD studies verified the stability of these complexes. Analysing the binding affinity and stability through free energy changes upon mutation, we found GRIM-19 mutations Y33P and Q61L and among GRIM-19 arginines, R68P and R57M, to be one of the top-most major and minor disruptors of binding, respectively. The proportionate increase in average change in binding affinity upon mutation was inclined more towards GRIM-19 mutants, leading to the surmise that GRIM-19 may play a greater role in the complex formation. These studies propound a novel structural perspective of STAT3α-GRIM-19 binding and inhibitory mechanisms in both the monomeric and dimeric forms of STAT3α as compared to that observed from the earlier experiments, these experimental observations being inconsistent among each other.

A 3D structural SARS-CoV-2-human interactome to explore genetic and drug perturbations.

Emergence of new viral agents is driven by evolution of interactions between viral proteins and host targets. For instance, increased infectivity of SARS-CoV-2 compared to SARS-CoV-1 arose in part through rapid evolution along the interface between the Spike protein and its human receptor ACE2, leading to increased binding affinity. To facilitate broader exploration of how pathogen-host interactions might impact transmission and virulence in the ongoing COVID-19 pandemic, we performed state-of-the-art interface prediction followed by molecular docking to construct a 3D structural interactome between SARS-CoV-2 and human. We additionally carried out downstream meta-analyses to investigate enrichment of sequence divergence between SARS-CoV-1 and SARS-CoV-2 or human population variants along viral-human protein interaction interfaces, predict changes in binding affinity by these mutations/variants, and further prioritize drug repurposing candidates predicted to competitively bind human targets. We believe this resource (http://3D-SARS2.yulab.org) will aid in development and testing of informed hypotheses for SARS-CoV-2 etiology and treatments.

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.

Elucidating important structural features for the binding affinity of spike - SARS-CoV-2 neutralizing antibody complexes.

The coronavirus disease 2019 (COVID‐19) has affected the lives of millions of people around the world. In an effort to develop therapeutic interventions and control the pandemic, scientists have isolated several neutralizing antibodies against SARS‐CoV‐2 from the vaccinated and convalescent individuals. These antibodies can be explored further to understand SARS‐CoV‐2 specific antigen–antibody interactions and biophysical parameters related to binding affinity, which can be utilized to engineer more potent antibodies for current and emerging SARS‐CoV‐2 variants. In the present study, we have analyzed the interface between spike protein of SARS‐CoV‐2 and neutralizing antibodies in terms of amino acid residue propensity, pair preference, and atomic interaction energy. We observed that Tyr residues containing contacts are highly preferred and energetically favorable at the interface of spike protein–antibody complexes. We have also developed a regression model to relate the experimental binding affinity for antibodies using structural features, which showed a correlation of 0.93. Moreover, several mutations at the spike protein–antibody interface were identified, which may lead to immune escape (epitope residues) and improved affinity (paratope residues) in current/emerging variants. Overall, the work provides insights into spike protein–antibody interactions, structural parameters related to binding affinity and mutational effects on binding affinity change, which can be helpful to develop better therapeutics against COVID‐19.

Rapid Assessment of Binding Affinity of SARS-COV-2 Spike Protein to the Human Angiotensin-Converting Enzyme 2 Receptor and to Neutralizing Biomolecules Based on Computer Simulations.

SARS-CoV-2 infects humans and causes Coronavirus disease 2019 (COVID-19). The S1 domain of the spike glycoprotein of SARS-CoV-2 binds to human angiotensin-converting enzyme 2 (hACE2) via its receptor-binding domain, while the S2 domain facilitates fusion between the virus and the host cell membrane for entry. The spike glycoprotein of circulating SARS-CoV-2 genomes is a mutation hotspot. Some mutations may affect the binding affinity for hACE2, while others may modulate S-glycoprotein expression, or they could result in a virus that can escape from antibodies generated by infection with the original variant or by vaccination. Since a large number of variants are emerging, it is of vital importance to be able to rapidly assess their characteristics: while changes of binding affinity alone do not always cause direct advantages for the virus, they still can provide important insights on where the evolutionary pressure is directed. Here, we propose a simple and cost-effective computational protocol based on Molecular Dynamics simulations to rapidly screen the ability of mutated spike protein to bind to the hACE2 receptor and selected neutralizing biomolecules. Our results show that it is possible to achieve rapid and reliable predictions of binding affinities. A similar approach can be used to perform preliminary screenings of the potential effects of S-RBD mutations, helping to prioritize the more time-consuming and expensive experimental work.

SCoV2-MD: a database for the dynamics of the SARS-CoV-2 proteome and variant impact predictions.

SCoV2-MD (www.scov2-md.org) is a new online resource that systematically organizes atomistic simulations of the SARS-CoV-2 proteome. The database includes simulations produced by leading groups using molecular dynamics (MD) methods to investigate the structure-dynamics-function relationships of viral proteins. SCoV2-MD cross-references the molecular data with the pandemic evolution by tracking all available variants sequenced during the pandemic and deposited in the GISAID resource. SCoV2-MD enables the interactive analysis of the deposited trajectories through a web interface, which enables users to search by viral protein, isolate, phylogenetic attributes, or specific point mutation. Each mutation can then be analyzed interactively combining static (e.g. a variety of amino acid substitution penalties) and dynamic (time-dependent data derived from the dynamics of the local geometry) scores. Dynamic scores can be computed on the basis of nine non-covalent interaction types, including steric properties, solvent accessibility, hydrogen bonding, and other types of chemical interactions. Where available, experimental data such as antibody escape and change in binding affinities from deep mutational scanning experiments are also made available. All metrics can be combined to build predefined or custom scores to interrogate the impact of evolving variants on protein structure and function.

High Temperature Requirement A (HtrA) protease of Listeria monocytogenes and its interaction with extracellular matrix molecules.

High Temperature Requirement A (HtrA) was identified as a secreted virulence factor in many pathogenic bacteria, including Listeria monocytogenes. Recently, it was discovered that Helicobacter pylori and Campylobacter jejuni HtrAs can directly cleave the human cell-adhesion molecule E-cadherin, which facilitates bacterial transmigration. HtrAs also interact with extracellular matrix (ECM) molecules. However, only a limited number of studies have been carried out in this regard. In the present study, the protease and ECM binding properties of L. monocytogenes HtrA (LmHtrA) were studied using native rLmHtrA, catalytically inactive rLmHtrA(S343A) and rLmHtrA lacking the PDZ domain (∆PDZ) to gain more insights into HtrA-ECM molecule interaction. The results show that (1) native rLmHtrA cleaves fibrinogen, fibronectin, plasminogen and casein in a time and temperature dependent manner, (2) interaction of rLmHtrA with various host proteins was found in the micromolar to nanomolar range, (3) in the absence of PDZ domain, rLmHtrA exhibits no drastic change in binding affinity toward the host molecules when compared with native rLmHtrA and (4) the PDZ domain plays an important role in the substrate cleavage as rLmHtrA1-394∆PDZ cleaves the substrates only under certain conditions. The proteolysis of various ECM molecules by rLmHtrA possibly highlights the role of HtrA in L. monocytogenes pathogenesis involving ECM degradation.

DsRBPBind: modeling the effect of RNA secondary structure on double-stranded RNA-protein binding.

RNA-binding proteins are fundamental to many cellular processes. Double-stranded RNA-binding proteins (dsRBPs) in particular are crucial for RNA interference, mRNA elongation, A-to-I editing, host defense, splicing and a multitude of other important mechanisms. Since dsRBPs require double-stranded RNA to bind, their binding affinity depends on the competition among all possible secondary structures of the target RNA molecule. Here, we introduce a quantitative model that allows calculation of the effective affinity of dsRBPs to any RNA given a principal affinity and the sequence of the RNA, while fully taking into account the entire secondary structure ensemble of the RNA. We implement our model within the ViennaRNA folding package while maintaining its O(N3) time complexity. We validate our quantitative model by comparing with experimentally determined binding affinities and stoichiometries for transactivation response element RNA-binding protein (TRBP). We also find that the change in dsRBP binding affinity purely due to the presence of alternative RNA structures can be many orders of magnitude and that the predicted affinity of TRBP for pre-miRNA-like constructs correlates with experimentally measured processing rates. Our modified version of the ViennaRNA package is available for download at http://bioserv.mps.ohio-state.edu/dsRBPBind, is free to use for research and educational purposes, and utilizes simple get/set methods for footprint size, concentration, cooperativity, principal dissociation constant and overlap.

Molecular insights into receptor binding of recent emerging SARS-CoV-2 variants

Multiple SARS-CoV-2 variants of concern (VOCs) have been emerging and some have been linked to an increase in case numbers globally. However, there is yet a lack of understanding of the molecular basis for the interactions between the human ACE2 (hACE2) receptor and these VOCs. Here we examined several VOCs including Alpha, Beta, and Gamma, and demonstrate that five variants receptor-binding domain (RBD) increased binding affinity for hACE2, and four variants pseudoviruses increased entry into susceptible cells. Crystal structures of hACE2-RBD complexes help identify the key residues facilitating changes in hACE2 binding affinity. Additionally, soluble hACE2 protein efficiently prevent most of the variants pseudoviruses. Our findings provide important molecular information and may help the development of novel therapeutic and prophylactic agents targeting these emerging mutants.

Predicting Mutational Effects on Receptor Binding of the Spike Protein of SARS-CoV-2 Variants.

The pandemic caused by SARS-CoV-2 has cost millions of lives and tremendous social/financial loss. The virus continues to evolve and mutate. In particular, the recently emerged “UK”, “South Africa”, and Delta variants show higher infectivity and spreading speed. Thus, the relationship between the mutations of certain amino acids and the spreading speed of the virus is a problem of great importance. In this respect, understanding the mutational mechanism is crucial for surveillance and prediction of future mutations as well as antibody/vaccine development. In this work, we used a coarse-grained model (that was used previously in predicting the importance of mutations of N501) to calculate the free energy change of various types of single-site or combined-site mutations. This was done for the UK, South Africa, and Delta mutants. We investigated the underlying mechanisms of the binding affinity changes for mutations at different spike protein domains of SARS-CoV-2 and provided the energy basis for the resistance of the E484 mutant to the antibody m396. Other potential mutation sites were also predicted. Furthermore, the in silico predictions were assessed by functional experiments. The results establish that the faster spreading of recently observed mutants is strongly correlated with the binding-affinity enhancement between virus and human receptor as well as with the reduction of the binding to the m396 antibody. Significantly, the current approach offers a way to predict new variants and to assess the effectiveness of different antibodies toward such variants.

ProNAB: database for binding affinities of protein-nucleic acid complexes and their mutants.

Protein–nucleic acid interactions are involved in various biological processes such as gene expression, replication, transcription, translation and packaging. The binding affinities of protein–DNA and protein–RNA complexes are important for elucidating the mechanism of protein–nucleic acid recognition. Although experimental data on binding affinity are reported abundantly in the literature, no well-curated database is currently available for protein–nucleic acid binding affinity. We have developed a database, ProNAB, which contains more than 20 000 experimental data for the binding affinities of protein–DNA and protein–RNA complexes. Each entry provides comprehensive information on sequence and structural features of a protein, nucleic acid and its complex, experimental conditions, thermodynamic parameters such as dissociation constant (Kd), binding free energy (ΔG) and change in binding free energy upon mutation (ΔΔG), and literature information. ProNAB is cross-linked with GenBank, UniProt, PDB, ProThermDB, PROSITE, DisProt and Pubmed. It provides a user-friendly web interface with options for search, display, sorting, visualization, download and upload the data. ProNAB is freely available at https://web.iitm.ac.in/bioinfo2/pronab/ and it has potential applications such as understanding the factors influencing the affinity, development of prediction tools, binding affinity change upon mutation and design complexes with the desired affinity.

WinBinVec: Cancer-Associated Protein-Protein Interaction Extraction and Identification of 20 Various Cancer Types and Metastasis Using Different Deep Learning Models.

Biophysical protein-protein interactions perform dominant roles in the initiation and progression of many cancer-related pathways. A protein-protein interaction might play different roles in diverse cancer types. Hence, prioritizing the PPIs in each cancer type would help detect cancer-associated pathways, find a better understanding of cancer biology, and facilitate drug discovery. Several studies to date have proposed computational methods for extracting the PPI essentiality of different cancer types based on the PPI network. The main drawback of these studies is not using a rich source such as genomics variant data. An amino acid sequence encodes useful information about protein structure and behavior. We represent each amino acid sequence based on its variants/mutations in seven different ways: binary vectors, pathogenicity scores, binding affinity changes upon mutations, gene expression-based network of the interactions, biophysicochemical properties, g-gap dipeptide, and one-hot vectors. Based on these representations, we design and consider seven different deep learning models. Then, we compare the accuracy of these models in predicting 20 different cancer types from the TCGA cohort. WinBinVec is a window-based model that outperforms the other models. Moreover, WinBinVec contains a PPI essentiality module that helps extract the essentiality probability of each PPI for every cancer type. Source code and Data: https://github.com/sabdollahi/WinBinVec.

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.

Determination of intracellular protein-ligand binding affinity by competition binding in-cell NMR.

Structure-based drug development suffers from high attrition rates due to the poor activity of lead compounds in cellular and animal models caused by low cell penetrance, off-target binding or changes in the conformation of the target protein in the cellular environment. The latter two effects cause a change in the apparent binding affinity of the compound, which is indirectly assessed by cellular activity assays. To date, direct measurement of the intracellular binding affinity remains a challenging task. In this work, in-cell NMR spectroscopy was applied to measure intracellular dissociation constants in the nanomolar range by means of protein-observed competition binding experiments. Competition binding curves relative to a reference compound could be retrieved either from a series of independent cell samples or from a single real-time NMR bioreactor run. The method was validated using a set of sulfonamide-based inhibitors of human carbonic anhydrase II with known activity in the subnanomolar to submicromolar range. The intracellular affinities were similar to those obtained in vitro, indicating that these compounds selectively bind to the intracellular target. In principle, the approach can be applied to any soluble intracellular target that gives rise to measurable chemical shift changes upon ligand binding.

Systematic Evaluation of DNA Sequence Variations on in vivo Transcription Factor Binding Affinity.

The majority of the single nucleotide variants (SNVs) identified by genome-wide association studies (GWAS) fall outside of the protein-coding regions. Elucidating the functional implications of these variants has been a major challenge. A possible mechanism for functional non-coding variants is that they disrupted the canonical transcription factor (TF) binding sites that affect the in vivo binding of the TF. However, their impact varies since many positions within a TF binding motif are not well conserved. Therefore, simply annotating all variants located in putative TF binding sites may overestimate the functional impact of these SNVs. We conducted a comprehensive survey to study the effect of SNVs on the TF binding affinity. A sequence-based machine learning method was used to estimate the change in binding affinity for each SNV located inside a putative motif site. From the results obtained on 18 TF binding motifs, we found that there is a substantial variation in terms of a SNV’s impact on TF binding affinity. We found that only about 20% of SNVs located inside putative TF binding sites would likely to have significant impact on the TF-DNA binding.

Computational studies of protein–drug binding affinity changes upon mutations in the drug target

AbstractMutations that lead to drug resistance limit the efficacy of antibiotics, antiviral drugs, targeted cancer therapies, and other treatments. Accurately calculating protein–drug binding affinity changes upon mutations in the drug target is of high interest as this can yield a better understanding into how such mutations drive drug‐resistance, especially when the mutation in question does not directly interfere with binding of the drug. The main aim of this article is to provide an up‐to‐date reference on the computational tools that are available for the calculation of Gibbs energy (free energy) changes upon mutation, their strengths, and limitations. The methods that are discussed include free energy calculations (free energy perturbation, thermodynamic integration, multistate Bennett acceptance ratio), analysis of molecular dynamics simulations (linear interaction energy, molecular mechanics [MM]/Poisson–Boltzmann solvated area, and MM/generalized Born solvated area), and methods that involve quantum mechanical calculations (including QM/MM). The possibility to use machine learning is also introduced. Given that the benefit of accurately calculating binding affinity changes upon mutation depends on comparing calculated values with experimental measurements, a brief survey on experimental methods and observables is provided. Examples of computational studies that go beyond calculating the Gibbs energy changes are given. Factors that need to be addressed by the computational chemist and potential pitfalls are discussed at length.This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Molecular and Statistical Mechanics > Free Energy Methods Molecular and Statistical Mechanics > Molecular Interactions

Nabe: an energetic database of amino acid mutations in protein-nucleic acid binding interfaces.

Protein–nucleic acid complexes play essential roles in regulating transcription, translation, DNA replication, repair and recombination, RNA processing and translocation. Site-directed mutagenesis has been extremely useful in understanding the principles of protein–DNA and protein–RNA interactions, and experimentally determined mutagenesis data are prerequisites for designing effective algorithms for predicting the binding affinity change upon mutation. However, a vital challenge in this area is the lack of sufficient public experimentally recognized mutation data, which leads to difficulties in developing computational prediction methods. In this article, we present Nabe, an integrated database of amino acid mutations and their effects on the binding free energy in protein–DNA and protein–RNA interactions for which binding affinities have been experimentally determined. Compared with existing databases and data sets, Nabe is the largest protein–nucleic acid mutation database, containing 2506 mutations in 473 protein–DNA and protein–RNA complexes, and of that 1751 are alanine mutations in 405 protein–nucleic acid complexes. For researchers to conveniently utilize the data, Nabe assembles protein–DNA and protein–RNA benchmark databases by adopting the data-processing procedures in the majority of models. To further facilitate users to query data, Nabe provides a searchable and graphical web page.Database URL: http://nabe.denglab.org