Computational Methods for Controlling Binding Specificity

Engineered TIMP2 with narrow MMP-9 specificity is an effective inhibitor of invasion and proliferation of triple negative breast cancer cells

Matrix metalloproteinases (MMPs) are a family of endopeptidases that degrade extracellular matrix (ECM) proteins, functioning in various physiological processes such as tissue remodeling, embryogenesis, and morphogenesis. Dysregulation of these enzymes is linked to multiple diseases. Specific inhibition of particular MMPs is crucial for anti-MMP drug development as some MMPs have shown anti-disease properties. In this study, we aimed to design a highly specific inhibitor of MMP-9, that plays a crucial role in cell invasion and metastasis, using tissue inhibitor of metalloproteinases 2 (TIMP2), an endogenous broad-family MMP inhibitor, as a prototype. In our earlier work, we were able to narrow down the specificity of the N-terminal domain of TIMP2 (N-TIMP2) toward MMP-9, yet at the expense of lowering its affinity to MMP-9. In this study, a library of N-TIMP2 mutants based on previous design with randomized additional positions was sorted for binding to MMP-9 using yeast surface display. Two selected N-TIMP2 mutants were expressed, purified and their inhibitory activity against a panel of MMPs was measured. The best engineered N-TIMP2 mutant (REY) exhibited a 2-fold higher affinity to MMP-9 compared to that of the WT N-TIMP2, and 6- to 1.1x104-fold increase in binding specificity toward MMP-9 compared to five alternative MMPs. Moreover, REY demonstrated a significant increase in inhibition of cell invasion and proliferation compared to the WT N-TIMP2 in MDA-MB-231 breast cancer cells. Therefore, our engineered N-TIMP2 mutant emerges as a promising candidate for future therapeutic development, offering precise targeting of MMP-9 in MMP-9-driven diseases.

Cold Spot SCANNER: Colab Notebook for predicting cold spots in protein-protein interfaces

BACKGROUND

Protein-protein interactions (PPIs) are conveyed through binding interfaces or surface patches on proteins that become buried upon binding. Structural and biophysical analysis of many protein-protein interfaces revealed certain unique features of these surfaces that determine the energetics of interactions and play a critical role in protein evolution. One of the significant aspects of binding interfaces is the presence of binding hot spots, where mutations are highly deleterious for binding. Conversely, binding cold spots are positions occupied by suboptimal amino acids and several mutations in such positions could lead to affinity enhancement. While there are many software programs for identification of hot spot positions, there is currently a lack of software for cold spot detection.

RESULTS

In this paper, we present Cold Spot SCANNER, a Colab Notebook, which scans a PPI binding interface and identifies cold spots resulting from cavities, unfavorable charge-charge, and unfavorable charge-hydrophobic interactions. The software offers a Py3DMOL-based interface that allows users to visualize cold spots in the context of the protein structure and generates a zip file containing the results for easy download.

CONCLUSIONS

Cold spot identification is of great importance to protein engineering studies and provides a useful insight into protein evolution. Cold Spot SCANNER is open to all users without login requirements and can be accessible at: https://colab.

RESEARCH

google.com/github/sagagugit/Cold-Spot-Scanner/blob/main/Cold_Spot_Scanner.ipynb .

Engineered variants of the Ras effector protein RASSF5 (NORE1A) promote anticancer activities in lung adenocarcinoma

Within the superfamily of small GTPases, Ras appears to be the master regulator of such processes as cell cycle progression, cell division, and apoptosis. Several oncogenic Ras mutations at amino acid positions 12, 13, and 61 have been identified that lose their ability to hydrolyze GTP, giving rise to constitutive signaling and eventually development of cancer. While disruption of the Ras/effector interface is an attractive strategy for drug design to prevent this constitutive activity, inhibition of this interaction using small molecules is impractical due to the absence of a cavity to which such molecules could bind. However, proteins and especially natural Ras effectors that bind to the Ras/effector interface with high affinity could disrupt Ras/effector interactions and abolish procancer pathways initiated by Ras oncogene. Using a combination of computational design and in vitro evolution, we engineered high-affinity Ras-binding proteins starting from a natural Ras effector, RASSF5 (NORE1A), which is encoded by a tumor suppressor gene. Unlike previously reported Ras oncogene inhibitors, the proteins we designed not only inhibit Ras-regulated procancer pathways, but also stimulate anticancer pathways initiated by RASSF5. We show that upon introduction into A549 lung carcinoma cells, the engineered RASSF5 mutants decreased cell viability and mobility to a significantly greater extent than WT RASSF5. In addition, these mutant proteins induce cellular senescence by increasing acetylation and decreasing phosphorylation of p53. In conclusion, engineered RASSF5 variants provide an attractive therapeutic strategy able to oppose cancer development by means of inhibiting of procancer pathways and stimulating anticancer processes.

Optimizing the Specificity Window of Biomolecular Receptors Using Structure-Switching and Allostery

To ensure maximum specificity (i.e., minimize cross-reactivity with structurally similar analogues of the desired target), most bioassays invoke "stringency", the careful tuning of the conditions employed (e.g., pH, ionic strength, or temperature). Willingness to control assay conditions will fall, however, as quantitative, single-step biosensors begin to replace multistep analytical processes. This is especially true for sensors deployed in vivo, where the tuning of such parameters is not just inconvenient but impossible. In response, we describe here the rational adaptation of two strategies employed by nature to tune the affinity of biomolecular receptors so as to optimize the placement of their specificity "windows" without the need to alter measurement conditions: structure-switching and allosteric control. We quantitatively validate these approaches using two distinct, DNA-based receptors: a simple, linear-chain DNA suitable for detecting a complementary DNA strand and a structurally complex DNA aptamer used for the detection of a small-molecule drug. Using these models, we show that, without altering assay conditions, structure-switching and allostery can tune the concentration range over which a receptor achieves optimal specificity over orders of magnitude, thus optimally matching the specificity window with the range of target concentrations expected to be seen in a given application.

Increasing the Affinity of an O-Antigen Polysaccharide Binding Site in Shigella flexneri Bacteriophage Sf6 Tailspike Protein

Broad and unspecific use of antibiotics accelerates spread of resistances. Sensitive and robust pathogen detection is thus important for a more targeted application. Bacteriophages contain a large repertoire of pathogen-binding proteins. These tailspike proteins (TSP) often bind surface glycans and represent a promising design platform for specific pathogen sensors. We analysed bacteriophage Sf6 TSP that recognizes the O-polysaccharide of dysentery-causing Shigella flexneri to develop variants with increased sensitivity for sensor applications. Ligand polyrhamnose backbone conformations were obtained from 2D 1 H,1 H-trNOESY NMR utilizing methine-methine and methine-methyl correlations. They agreed well with conformations obtained from molecular dynamics (MD), validating the method for further predictions. In a set of mutants, MD predicted ligand flexibilities that were in good correlation with binding strength as confirmed on immobilized S. flexneri O-polysaccharide (PS) with surface plasmon resonance. In silico approaches combined with rapid screening on PS surfaces hence provide valuable strategies for TSP-based pathogen sensor design.

Generating quantitative binding landscapes through fractional binding selections combined with deep sequencing and data normalization

Quantifying the effects of various mutations on binding free energy is crucial for understanding the evolution of protein-protein interactions and would greatly facilitate protein engineering studies. Yet, measuring changes in binding free energy (ΔΔG_bind) remains a tedious task that requires expression of each mutant, its purification, and affinity measurements. We developed an attractive approach that allows us to quantify ΔΔG_bind for thousands of protein mutants in one experiment. Our protocol combines protein randomization, Yeast Surface Display technology, deep sequencing, and a few experimental ΔΔG_bind data points on purified proteins to generate ΔΔG_bind values for the remaining numerous mutants of the same protein complex. Using this methodology, we comprehensively map the single-mutant binding landscape of one of the highest-affinity interaction between BPTI and Bovine Trypsin (BT). We show that ΔΔG_bind for this interaction could be quantified with high accuracy over the range of 12 kcal mol⁻¹ displayed by various BPTI single mutants.

Quantifying the effect of mutations on binding free energy is important to understand protein-protein interaction (PPI). Here the authors develop a method based on yeast display and next-generation sequencing to generate quantitative binding landscapes for any PPI regardless of their Kd value.

Converting a broad matrix metalloproteinase family inhibitor into a specific inhibitor of MMP-9 and MMP-14

MMP-14 and MMP-9 are two well-established cancer targets for which no specific clinically relevant inhibitor is available. Using a powerful combination of computational design and yeast surface display technology, we engineered such an inhibitor starting from a nonspecific MMP inhibitor, N-TIMP2. The engineered purified N-TIMP2 variants showed enhanced specificity toward MMP-14 and MMP-9 relative to a panel of off-target MMPs. MMP-specific N-TIMP2 sequence signatures were obtained that could be understood from the structural perspective of MMP/N-TIMP2 interactions. Our MMP-9 inhibitor exhibited 1000-fold preference for MMP-9 vs. MMP-14, which is likely to translate into significant differences under physiological conditions. Our results provide new insights regarding evolution of promiscuous proteins and optimization strategies for design of inhibitors with single-target specificities.

Development of High Affinity and High Specificity Inhibitors of Matrix Metalloproteinase 14 through Computational Design and Directed Evolution

Degradation of the extracellular matrices in the human body is controlled by matrix metalloproteinases (MMPs), a family of more than 20 homologous enzymes. Imbalance in MMP activity can result in many diseases, such as arthritis, cardiovascular diseases, neurological disorders, fibrosis, and cancers. Thus, MMPs present attractive targets for drug design and have been a focus for inhibitor design for as long as 3 decades. Yet, to date, all MMP inhibitors have failed in clinical trials because of their broad activity against numerous MMP family members and the serious side effects of the proposed treatment. In this study, we integrated a computational method and a yeast surface display technique to obtain highly specific inhibitors of MMP-14 by modifying the natural non-specific broad MMP inhibitor protein N-TIMP2 to interact optimally with MMP-14. We identified an N-TIMP2 mutant, with five mutations in its interface, that has an MMP-14 inhibition constant (K_i ) of 0.9 pm, the strongest MMP-14 inhibitor reported so far. Compared with wild-type N-TIMP2, this variant displays ∼900-fold improved affinity toward MMP-14 and up to 16,000-fold greater specificity toward MMP-14 relative to other MMPs. In an in vitro and cell-based model of MMP-dependent breast cancer cellular invasiveness, this N-TIMP2 mutant acted as a functional inhibitor. Thus, our study demonstrates the enormous potential of a combined computational/directed evolution approach to protein engineering. Furthermore, it offers fundamental clues into the molecular basis of MMP regulation by N-TIMP2 and identifies a promising MMP-14 inhibitor as a starting point for the development of protein-based anticancer therapeutics.

The Framework of Computational Protein Design

Computational protein design (CPD) has established itself as a leading field in basic and applied science with a strong coupling between the two. Proteins are computationally designed from the level of amino acids to the level of a functional protein complex. Design targets range from increased thermo- (or other) stability to specific requested reactions such as protein-protein binding, enzymatic reactions, or nanotechnology applications. The design scheme may encompass small regions of the proteins or the entire protein. In either case, the design may aim at the side-chains or at the full backbone conformation. Herein, the main framework for the process is outlined highlighting key elements in the CPD iterative cycle. These include the very definition of CPD, the diverse goals of CPD, components of the CPD protocol, methods for searching sequence and structure space, scoring functions, and augmenting the CPD with other optimization tools. Taken together, this chapter aims to introduce the framework of CPD.

Identifying Residues that Determine SCF Molecular-Level Interactions through a Combination of Experimental and In silico Analyses

The stem cell factor (SCF)/c-Kit receptor tyrosine kinase complex-with its significant roles in hematopoiesis and angiogenesis-is an attractive target for rational drug design. There is thus a need to map, in detail, the SCF/c-Kit interaction sites and the mechanisms that modulate this interaction. While most residues in the direct SCF/c-Kit binding interface can be identified from the existing crystal structure of the complex, other residues that affect binding through protein unfolding, intermolecular interactions, allosteric or long-distance electrostatic effects cannot be directly inferred. Here, we describe an efficient method for protein-wide epitope mapping using yeast surface display. A library of single SCF mutants that span the SCF sequence was screened for decreased affinity to soluble c-Kit. Sequencing of selected clones allowed the identification of mutations that reduce SCF binding affinity to c-Kit. Moreover, the screening of these SCF clones for binding to a structural antibody helped identify mutations that result in small or large conformational changes in SCF. Computational modeling of the experimentally identified mutations showed that these mutations reduced the binding affinity through one of the three scenarios: through SCF destabilization, through elimination of favorable SCF/c-Kit intermolecular interactions, or through allosteric changes. Eight SCF variants were expressed and purified. Experimentally measured in vitro binding affinities of these mutants to c-Kit confirmed both the yeast surface display selection results and the computational predictions. This study has thus identified the residues crucial for c-Kit/SCF binding and has demonstrated the advantages of using a combination of computational and combinatorial methods for epitope mapping.

Cold Spots in Protein Binding

Understanding the energetics and architecture of protein-binding interfaces is important for basic research and could potentially facilitate the design of novel binding domains for biotechnological applications. It is well accepted that a few key residues at binding interfaces (binding hot spots) are responsible for contributing most to the free energy of binding. In this opinion article, we introduce a new concept of 'binding cold spots', or interface positions occupied by suboptimal amino acids. Such positions exhibit a potential for affinity enhancement through various mutations. We give several examples of cold spots from different protein-engineering studies and argue that identification of such positions is crucial for studies of protein evolution and protein design.

Saturation scanning of ubiquitin variants reveals a common hot spot for binding to USP2 and USP21

A detailed understanding of the molecular mechanisms whereby ubiquitin (Ub) recognizes enzymes in the Ub proteasome system is crucial for understanding the biological function of Ub. Many structures of Ub complexes have been solved and, in most cases, reveal a large structural epitope on a common face of the Ub molecule. However, owing to the generally weak nature of these interactions, it has been difficult to map in detail the functional contributions of individual Ub side chains to affinity and specificity. Here we took advantage of Ub variants (Ubvs) that bind tightly to particular Ub-specific proteases (USPs) and used phage display and saturation scanning mutagenesis to comprehensively map functional epitopes within the structural epitopes. We found that Ubvs that bind to USP2 or USP21 contain a remarkably similar core functional epitope, or "hot spot," consisting mainly of positions that are conserved as the wild type sequence, but also some positions that prefer mutant sequences. The Ubv core functional epitope contacts residues that are conserved in the human USP family, and thus it is likely important for the interactions of Ub across many family members.

Protein Engineering by Combined Computational and In Vitro Evolution Approaches

Two alternative strategies are commonly used to study protein-protein interactions (PPIs) and to engineer protein-based inhibitors. In one approach, binders are selected experimentally from combinatorial libraries of protein mutants that are displayed on a cell surface. In the other approach, computational modeling is used to explore an astronomically large number of protein sequences to select a small number of sequences for experimental testing. While both approaches have some limitations, their combination produces superior results in various protein engineering applications. Such applications include the design of novel binders and inhibitors, the enhancement of affinity and specificity, and the mapping of binding epitopes. The combination of these approaches also aids in the understanding of the specificity profiles of various PPIs.

Combinatorial and Computational Approaches to Identify Interactions of Macrophage Colony-stimulating Factor (M-CSF) and Its Receptor c-FMS

The molecular interactions between macrophage colony-stimulating factor (M-CSF) and the tyrosine kinase receptor c-FMS play a key role in the immune response, bone metabolism, and the development of some cancers. Because no x-ray structure is available for the human M-CSF · c-FMS complex, the binding epitope for this complex is largely unknown. Our goal was to identify the residues that are essential for binding of the human M-CSF to c-FMS. For this purpose, we used a yeast surface display (YSD) approach. We expressed a combinatorial library of monomeric M-CSF (M-CSFM) single mutants and screened this library to isolate variants with reduced affinity for c-FMS using FACS. Sequencing yielded a number of single M-CSFM variants with mutations both in the direct binding interface and distant from the binding site. In addition, we used computational modeling to map the identified mutations onto the M-CSFM structure and to classify the mutations into three groups as follows: those that significantly decrease protein stability; those that destroy favorable intermolecular interactions; and those that decrease affinity through allosteric effects. To validate the YSD and computational data, M-CSFM and three variants were produced as soluble proteins; their affinity and structure were analyzed; and very good correlations with both YSD data and computational predictions were obtained. By identifying the M-CSFM residues critical for M-CSF · c-FMS interactions, we have laid down the basis for a deeper understanding of the M-CSF · c-FMS signaling mechanism and for the development of target-specific therapeutic agents with the ability to sterically occlude the M-CSF·c-FMS binding interface.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

Affinity- and specificity-enhancing mutations are frequent in multispecific interactions between TIMP2 and MMPs

Multispecific proteins play a major role in controlling various functions such as signaling, regulation of transcription/translation, and immune response. Hence, a thorough understanding of the atomic-level principles governing multispecific interactions is important not only for the advancement of basic science but also for applied research such as drug design. Here, we study evolution of an exemplary multispecific protein, a Tissue Inhibitor of Matrix Metalloproteinases 2 (TIMP2) that binds with comparable affinities to more than twenty-six members of the Matrix Metalloproteinase (MMP) and the related ADAMs families. We postulate that due to its multispecific nature, TIMP2 is not optimized to bind to any individual MMP type, but rather embodies a compromise required for interactions with all MMPs. To explore this hypothesis, we perform computational saturation mutagenesis of the TIMP2 binding interface and predict changes in free energy of binding to eight MMP targets. Computational results reveal the non-optimality of the TIMP2 binding interface for all studied proteins, identifying many affinity-enhancing mutations at multiple positions. Several TIMP2 point mutants predicted to enhance binding affinity and/or binding specificity towards MMP14 were selected for experimental verification. Experimental results show high abundance of affinity-enhancing mutations in TIMP2, with some point mutations producing more than ten-fold improvement in affinity to MMP14. Our computational and experimental results collaboratively demonstrate that the TIMP2 sequence lies far from the fitness maximum when interacting with its target enzymes. This non-optimality of the binding interface and high potential for improvement might characterize all proteins evolved for binding to multiple targets.

Mapping of the binding landscape for a picomolar protein-protein complex through computation and experiment

Our understanding of protein evolution would greatly benefit from mapping of binding landscapes, i.e., changes in protein-protein binding affinity due to all single mutations. However, experimental generation of such landscapes is a tedious task due to a large number of possible mutations. Here, we use a simple computational protocol to map the binding landscape for two homologous high-affinity complexes, involving a snake toxin fasciculin and acetylcholinesterase from two different species. To verify our computational predictions, we experimentally measure binding between 25 Fas mutants and the 2 enzymes. Both computational and experimental results demonstrate that the Fas sequence is close to the optimum when interacting with its targets, yet a few mutations could further improve Kd, kon, and koff. Our computational predictions agree well with experimental results and generate distributions similar to those observed in other high-affinity PPIs, demonstrating the potential of simple computational protocols in capturing realistic binding landscapes.