- Department of Cancer Biology, Mayo Clinic, Jacksonville, FL, United States
The study of protein-protein interactions (PPIs) and the engineering of protein-based inhibitors often employ two distinct strategies. One approach leverages the power of combinatorial libraries, displaying large ensembles of mutant proteins, for example, on the yeast cell surface, to select binders. Another approach harnesses computational modeling, sifting through an astronomically large number of protein sequences and attempting to predict the impact of mutations on PPI binding energy. Individually, each approach has inherent limitations, but when combined, they generate superior outcomes across diverse protein engineering endeavors. This synergistic integration of approaches aids in identifying novel binders and inhibitors, fine-tuning specificity and affinity for known binding partners, and detailed mapping of binding epitopes. It can also provide insight into the specificity profiles of varied PPIs. Here, we outline strategies for directing the evolution of tissue inhibitors of metalloproteinases (TIMPs), which act as natural inhibitors of matrix metalloproteinases (MMPs). We highlight examples wherein design of combinatorial TIMP libraries using structural and computational insights and screening these libraries of variants using yeast surface display (YSD), has successfully optimized for MMP binding and selectivity, and conferred insight into the PPIs involved.
1 Introduction
The matrix metalloproteinase (MMP) family plays a pivotal role in embryonic development and tissue regeneration. However, the unchecked activity of MMPs is linked to various diseases, including cancer (Gialeli et al., 2011), cardiovascular and pulmonary diseases (Hu et al., 2007; Murphy and Nagase, 2008), rheumatoid arthritis (Grillet et al., 2023), and encephalomyelitis (Clements et al., 1997). Although targeting MMPs presents an attractive therapeutic approach, prior efforts to develop therapeutic MMP inhibitors met with limited success. The primary issue is the similarity of catalytic site structure among MMPs, and consequent poor selectivity of small molecule inhibitors, frequently resulting in musculoskeletal toxicity (Overall and Kleifeld, 2006; Krzeski et al., 2007). This makes the development of highly selective MMP inhibitors a complex challenge.
The engineered TIMPs present a novel alternative to small molecule inhibitors for targeting MMPs in disease (Levin et al., 2017; Radisky et al., 2017). Natural TIMPs are strong protein-based inhibitors of MMPs, and engineered mutations of TIMPs can improve their affinity and selectivity for MMPs (Nagase and Brew, 2002). Furthermore, TIMP inhibitory profiles can be more comprehensively manipulated by using the powerful protein engineering platform, YSD. By leveraging structural knowledge and computational tools, we can design efficient focused libraries, trimming down the vast array of potential mutation combinations (Bonadio and Shifman, 2021). In this review, we highlight recent advances that have employed structure and computation-guided YSD approaches in combination as a cutting-edge method for designing and creating TIMP-based MMP inhibitors.
2 Integrating yeast surface display (YSD) with structural and computational design
Combinatorial protein libraries in YSD are presented on the surface of Saccharomyces cerevisiae cells; such libraries can contain up to 108 distinct protein variants (Gai and Wittrup, 2007). In comparison to other protein display technologies, the large size of yeast cells enables the use of flow cytometry and fluorescence-activated cell sorting (FACS) to identify antigen-binding cells. YSD typically employs a two-color fluorescence labeling method, where one fluorophore measures antigen binding, while the other measures protein expression (Cherf and Cochran, 2015). This simultaneous approach allows concurrent screening for both affinity and protein stability. Following the screening process, the binding affinity of individual protein variants may also be quantitatively assessed directly on the yeast surface. This efficient workflow allows for rapid prioritization of clones without the need for extensive production and purification of soluble protein for many individual variants (Zahradnik et al., 2021).
A limitation of YSD remains the number of distinct variants that can be generated and screened within a library. While a YSD library may encompass up to 108 variants, this number is dwarfed by the vast sequence space that proteins inhabit. Furthermore, mutations at disparate positions can interact in complex and sometimes unpredictable ways, further complicating the landscape. The limitation in library size thus necessitates efficient library design to make the most of the available screening capacity and ensure that libraries are relevant and enriched in potential solutions (Boder and Wittrup, 1997; Park et al., 2006).
Here, structural insights and computational strategies play a pivotal role. Structural data can highlight key interaction sites, pinpointing regions most amenable to mutations that are likely to affect binding characteristics. Furthermore, computational strategies can simulate and predict the effects of myriad mutations, allowing prioritization of specific high-value sites for diversification, as well as pairs or groups of sites likely to interact. By integrating these structural and computational insights, library design within YSD is elevated to a more strategic level, creating libraries enriched in sequences with the highest potential to achieve desired binding outcomes (Rosenfeld et al., 2016).
After YSD selection, structural and computational methods become essential in interpreting the results. Through next-generation sequencing (NGS), we can gain insights into enriched sequences, revealing preferred mutations and potential interaction domains. Structural analyses of selected variants, combined with molecular modeling, reveal the conformational adjustments and interactions driving enhanced binding (Rosenfeld et al., 2016). Importantly, these insights can be channeled back into refining library designs, paving the way for iterative improvements and increasingly effective protein engineering endeavors.
3 Matrix metalloproteinases (MMPs)
MMPs, or matrixins, belong to a larger family of metalloproteinases called the metzincins, which share characteristic structural features of the catalytic domain including a conserved zinc-binding motif at the catalytic site (Cerda-Costa and Gomis-Ruth, 2014; Radisky and Coban, 2021). MMPs degrade different types of extracellular matrix proteins, as well as other substrates including cytokines and cell-surface receptors (Nagase, 1996). These enzymes play a critical role in various biological processes, such as embryogenesis, morphogenesis, and wound healing. Moreover, they participate in pathological processes, such as tumor metastasis, angiogenesis, and inflammation (Malemud, 2006).
3.1 MMP classification
MMPs are commonly classified into several groups based on their substrate specificity toward extracellular matrix (ECM) proteins and similarities in their domain structures. Accordingly, MMPs are categorized into gelatinases, collagenases, stromelysins, matrilysins, membrane-type MMPs (MT-MMPs), and others (Table1) (Sternlicht and Werb, 2001; Overall, 2002; Visse and Nagase, 2003; Nagase et al., 2006). All possess N-terminal signal peptides, prodomains, and catalytic domains; other accessory domains vary among the groups.
TABLE 1. Members of the human MMP family, and their domains and extracellular substrates (Sternlicht and Werb, 2001; Overall, 2002; Visse and Nagase, 2003; Nagase et al., 2006).
Gelatinases (MMP-2 and MMP-9) can digest gelatin (denatured collagen), some collagens, elastin, fibronectin, and laminin. They possess fibronectin-like repeats as an insertion within in their catalytic domain, which assist in binding to gelatin; they also possess a C-terminal hemopexin domain connected to the catalytic domain by a hinge domain. Gelatinases play key roles in osteogenesis, embryogenesis, angiogenesis, and wound healing. Overexpression of gelatinase was reported in many invasive and metastatic tumors (Nagase et al., 2006; Klein and Bischoff, 2011).
True collagenases (those that efficiently digest intact triple-helical collagen) include MMP-1, MMP-8, and MMP-13, along with MT1-MMP (Amar et al., 2017). C-terminal to the catalytic domain, they possess a hinge domain and hemopexin domain. The collagenases can cleave fibrillar collagen, which is the main component of cartilage and bone (Fischer et al., 2019).
Stromelysins have the same domain organization as collagenases but do not cleave fibrillar collagen, instead digesting a wide variety of collagen fragments and other extracellular matrix components. MMP-3, MMP-10, and MMP-11 belong to this family (Nagase et al., 2006).
Matrilysins (MMP-7 and MMP-26) have the simplest domain structure, comprised of only the signal sequence, prodomain, and catalytic domain. They cleave various components of the extracellular matrix and other substrates. For example, MMP-7 secreted by epithelial cells of the intestine can cleave α-defensin precursors to bactericidal forms, and thus play an important role in innate immunity (Wilson et al., 1999).
Membrane-type metalloproteinases (MT-MMPs) include six MMPs anchored to the plasma membrane. This group can be subdivided into type-I transmembrane proteins including MT1-MMP, MT2-MMP, MT3-MMP, and MT5-MMP, and glycosylphosphatidylinositol (GPI)-anchored proteins including MT4-MMP and MT6-MMP. The type-I transmembrane group have cytoplasmic domains which can play essential roles in cellular signaling (Murphy and Nagase, 2008; Itoh, 2015).
3.2 MMP structure
The 3D structures of many MMPs have been revealed by NMR spectroscopy and X-ray crystallography (Maskos, 2005). The prodomain contains ∼88 amino acids and includes a conserved motif (PRCGxPD) containing a cysteine residue which coordinates to the active site Zn ion, maintaining the enzyme in an inactive form (Tallant et al., 2010). Upon proteolytic cleavage between the prodomain and the catalytic domain, the binding between Cys and zinc is disrupted, and the enzyme becomes activated (Ra and Parks, 2007).
The catalytic domain contains the substrate binding cleft and the active site and consists of approximately 170 amino acids (Lovejoy et al., 1994). A conserved sequence, HExxHxxGxxH, consists of three histidine residues interacting with zinc (Bode et al., 1993). There are two zinc atoms; one is catalytic, and the other is structural and additionally, MMPs bind calcium atoms for structural integrity (Massova et al., 1998). The hemopexin-like domain comprises approximately 190 amino acids forming a 4-propeller structure located at or near the C-terminus and is pivotal in interacting with various other proteins (Libson et al., 1995). The hinge region is a flexible linker that differs considerably in length among different MMPs and connects the hemopexin-like domain to the catalytic domain (Nagase et al., 2006). The gelatinases (MMP-2 and MMP-9) additionally possess three fibronectin type II domains inserted within the catalytic domain after the fifth β-sheet (Elkins et al., 2002).
4 Tissue inhibitors of metalloproteinases (TIMPs)
Imbalances in production of MMPs and their inhibitors dysregulate ECM turnover and contribute to pathology in various diseases including cancer, cardiovascular diseases, neurological conditions, and arthritis. In healthy tissues, appropriate levels of MMP activity are controlled at multiple levels, including transcriptional regulation, zymogen activation of precursor proteins, localization and binding to specific ECM proteins, and inhibition by endogenous protein inhibitors. TIMPs are the major endogenous inhibitors of MMPs and play a critical role in controlling the balance between the degradation and remodeling of ECM (Murphy and Nagase, 2008).
4.1 TIMP classification
Four members of the human TIMP family are encoded by paralogous genes on different chromosomes: TIMP-1 on chromosome Xp, TIMP-2 on chromosome 17p, TIMP-3 on chromosome 22q, and TIMP-4 on chromosome 3p (Murphy, 2011). Orthologs of all four TIMPs are conserved throughout mammals, while homologs are distributed across vertebrate and invertebrate animals (Brew and Nagase, 2010). All four TIMPs inhibit MMPs broadly, but with some differences in affinities toward specific MMPs. For example, MT1-MMP and MMP-19 are only very weakly inhibited by TIMP-1, but more potently inhibited by the other three TIMPs (Batra and Radisky, 2011). All TIMPs are soluble and secreted extracellularly; uniquely, TIMP-3 can also bind tightly to glycosaminoglycans and become incorporated within ECM (Yu et al., 2000; Murphy and Nagase, 2008).
4.2 TIMP structure
The TIMP structure is comprised of two domains situated side-by-side (Figure 1) (Batra and Radisky, 2011). The N-terminal domain, also known as the inhibitory domain, consists of approximately 125 amino acids and directly interacts with the zinc ion in the catalytic site, leading to inhibition of catalytic activity. The C-terminal domain has approximately 65 amino acids and for each TIMP may mediate unique interactions with binding proteins or domains distinct from the MMP catalytic domains. For example, TIMP-2 can serve as an adapter protein to assist in the activation of proMMP-2 by MT1-MMP; the TIMP-2 inhibitory domain docks with one catalytic domain of an MT1-MMP dimer, while the TIMP-2 C-terminal domain binds to the hemopexin domain of proMMP-2, bringing the zymogen into proximity with the activating catalytic domain of the second MT1-MMP molecule (Morgunova et al., 2002). TIMPs can also bind to cell-surface receptors and stimulate cell signaling; these activities vary among the four TIMPs. The structure of the N-terminal domain contains a five-stranded β-barrel and three α-helices, while the C-terminal domain has two pairs of β-strands (parallel and antiparallel), which are connected by an α-helix. These domains are about 40% conserved among the four human TIMPs (Batra et al., 2013; Raeeszadeh-Sarmazdeh et al., 2019).
FIGURE 1. Details of the MMP/TIMP interaction. The MMP-3 catalytic domain (light green) bound to TIMP-1, colored by the N-terminal domain (yellow) and C-terminal domain (pink), illustrates the broad protein-protein interaction interface (PDB ID 1UEA). Structural calciums (bright green), the structural zinc (grey) and catalytic zinc (brick red) are displayed for context. Detail panel shows the conserved zinc-chelating histidines and catalytic Glu202 of the MMP-3 active site, with Cys1 of TIMP-1 chelating the zinc ion.
X-ray crystallographic structures of MMP/TIMP complexes reveal that the most crucial protein-protein interactions are made with a conserved ridge on the N-terminal domain of TIMP, which consists of the N-terminal five amino acids Cys1-X-Cys3-X-Pro5, and the C-connector loop which connects β-strands sC to sD. A disulfide bond between Cys1 and Cys70 covalently fastens these two segments, which together fill both the nonprimed and primed subsites of the MMP substrate-binding cleft. The N-terminal amine of Cys1 is positioned at the catalytic site of the MMP, and coordinates to the catalytic zinc ion, displacing the hydrolytic water molecule (Gomis-Ruth et al., 1997). Additional interaction regions which extend the binding ridge are comprised of the TIMP AB-loop of the N-terminal domain and the GH-loop and multiple-turn-loop of the C-terminal domain; these regions have variable involvement in different MMP-TIMP complexes (Brew and Nagase, 2010; Batra and Radisky, 2011; Batra et al., 2013).
5 Applications of yeast surface display (YSD) for TIMPs
A primary reason for the failure of MMP inhibitors in clinical trials has been the lack of adequately selective drugs (Clements et al., 1997; Moore et al., 2003; Sparano et al., 2004; Bissett et al., 2005). Therefore, identifying highly selective MMP inhibitors has become a significant goal to improve the effectiveness of these agents and minimize off-target effects. Here, we overview recent work in the field, focusing on novel protein engineering strategies for developing selective TIMP-based therapeutics.
Raeeszadeh-Sarmazdeh et al. (2019) screened a large library of full-length TIMP-1 variants displayed on the yeast surface to optimize MMP-3 binding affinity. The library was designed to incorporate diversity across 8 positions in the N-terminal domain and 9 positions in the C-terminal domain, all of which participate in the interaction interface with the MMP catalytic domain (Figure 2). The large number of diversified positions did not allow for comprehensive coverage of sequence-space, but instead allowed broad sampling of variants with multiple mutations (∼4 on average) and enabled us to identify cooperative mutations that work together. Follow-up experiments probed this cooperativity by comparing binding profiles of single and double mutants, while structural studies analyzed the crystal structure of an improved variant to yield insight into the structural basis for binding improvements. Our results showed that the Leu34Gly mutation on the N-terminal domain AB-loop increases the affinity of TIMP-1 to MMP-3 by facilitating formation of a reciprocal clasp between Tyr153 of MMP-3 and Tyr35 of TIMP-1. Moreover, a Gly154Ala mutation on the C-terminal domain changes the conformation of the multiple-turn loop resulting in stronger interaction between two domains and improving MMP-3 binding in comparison with the wild-type TIMP-1 (Raeeszadeh-Sarmazdeh et al., 2019). This study revealed that, despite the common belief that the N-terminal domain of TIMPs is uniquely important for MMP binding and inhibition, the C-terminal domain and the interaction between the two domains can also be crucial in targeting an MMP with the highest affinity.
FIGURE 2. Schematic overview of yeast surface display for library screening and directed evolution. Directed protein evolution cycle utilizing YSD entails producing, expressing, displaying on the yeast surface, and screening libraries, and then sequencing and assessing the isolated variants. TIMP variants are displayed on the yeast surface via fusion at the N-terminus of yeast cell wall protein Aga2; the native N-terminus of the mature TIMP protein must be accessible to enable binding to MMPs during the screening process. TIMP expression is detected through immunolabeling of the c-myc epitope tag, and binding of an MMP is detected using a different fluorescent label (red highlight). (Illustration created with BioRender.com).
A subsequent study delved into the challenge of distinguishing between MMP-3 and MMP-10, two MMPs noted for their pronounced structural and functional conservation. Utilizing a counter-selective strategy to screen the library described above, we identified TIMP-1 variants with up to 23-fold improved selectivity for MMP-3 relative to MMP-10 (Raeeszadeh-Sarmazdeh et al., 2022). The two most improved variants each possessed the Leu34Gly mutation, as found in the previous study, along with additional mutations in the N-terminal domain C-connector loop and the C-terminal domain GH- and multiple turn loops. Crystal structures of these TIMP-1 variants in complex with MMP-3 provided insights into how the mutations were accommodated at the binding interface, while structural modeling revealed explanations for losses in affinity toward the counter-target MMP-10. For example, mutations at Met66 in the C-connector loop and Pro131 and Leu133 in the GH loop each disrupt predicted hydrophobic packing interactions with MMP-10, while having less impact on the interaction with MMP-3. Together, these mutations preserved TIMP-1 affinity toward MMP-3 while simultaneously diminishing affinity toward MMP-10 (Raeeszadeh-Sarmazdeh et al., 2022). This study underscores the power of the YSD platform in identifying rare variants that excel in discriminating between closely related targets. Moreover, it emphasizes the indispensable role of structural data and computational modeling in deciphering the underpinnings of enhanced selectivity.
Other efforts to engineer selective TIMPs have employed the more minimal scaffold of the N-terminal domain of TIMP-2 (N-TIMP-2), which retains inhibitory capability toward MMPs. In a seminal study, Arkadash et al. (2017) combined computational approaches for library design with YSD to identify high affinity and selectivity inhibitors of MT1-MMP. First, structural modeling and computational saturation mutagenesis of the TIMP-2 complex with multiple MMPs allowed prediction of binding energy changes that would result from mutations at each site (Sharabi et al., 2014). This systematic analysis identified specific positions at which mutations conferred high probability of improving affinity and specificity toward MT1-MMP. Among these, seven positions were selected for diversification within the focused library: positions 4, 35, 38, 68, 71, 97, and 99 (Arkadash et al., 2017). The top variant identified through YSD screening contained 5 mutations, improved MT1-MMP affinity by approximately 900-fold with a Ki value of 0.9 pM, and improved selectivity relative to a range of other MMPs by 2–4 orders of magnitude. Structural modeling of key mutations observed in the selected variants lent insight into affinity enhancements toward MT1-MMP. For example, the Ser4Arg mutation ameliorates the interaction of N-TIMP-2 and MT1-MMP by favoring intermolecular hydrogen bonds with Asn231 and Asn229 on MT1-MMP. Mutations of Ile35 to alternative hydrophobic amino acids Met, Pro, or Leu improved packing of MT1-MMP and N-TIMP-2 at the interface. Substitution of Asp at position Asn38 was predominant in selected variants and is predicted to enable formation of a salt bridge with MT1-MMP Lys41, thus stabilizing the interacting loop conformation in the context of the bound complex. The Ser68Asp substitution improves packing and stabilizes the N-TIMP-2 conformation by forming an intramolecular hydrogen bond with Ala70, while the nearby Ile71Val substitution likewise improves packing with MT1-MMP (Arkadash et al., 2017).
Arkadash et al. (2018) subsequently screened the same YSD library of N-TIMP-2 variants using a dual-target screening strategy to obtain selective inhibitors with high affinity for MMP-9 or MT1-MMP. Top variants to emerge from this sorting strategy showed strongly enhanced affinity and selectivity toward the intended targets. The best MT1-MMP-selective variants had ∼200-fold improved affinity toward MT1-MMP, with Ki values of 24–30 pM. Those selected for targeting MMP-9 showed little change in affinity toward the target, with Ki values of ∼1 nM, but nevertheless showed significant improvements in specificity. These biochemical results were consistent with the specificity and selectivity of the engineered N-TIMP-2 variants in cell-based assays. In assays with U87MG glioma cells, which depend upon MT1-MMP for activation of MMP-2, the MT1-MMP-selective N-TIMP-2 variant was uniquely able to reduce the activation of MMP-2 by up to 50% (Arkadash et al., 2018). In assays employing MCF-7 breast cancer cells overexpressing MMP-9, the MMP-9-selective N-TIMP-2 was uniquely capable of blocking gelatin degradation and inhibiting cellular migration. Conversely, in MCF-7 cells overexpressing MT1-MMP, the MT1-MMP-selective N-TIMP-2 variant uniquely inhibited cellular migration (Arkadash et al., 2018). These experiments demonstrate that the precision-engineering of TIMPs, guided by the YSD platform, can yield molecules with remarkable target specificity and potency in a cellular context. Moreover, they highlight the potential of tailored protein inhibitors to selectively modulate specific pathways in complex cellular environments, opening avenues for targeted therapeutic applications.
An even more narrowly targeted approach to N-TIMP-2 library design was reported by Shirian et al. (2018), toward a similar goal of developing selective inhibitors for MMP-9 and MT1-MMP. The library incorporated diversity at eight positions computationally predicted to enhance affinity and/or specificity toward the targets (positions 4, 6, 35, 38, 68, 71, 97, and 99). Furthermore, substitutions at each site were restricted to enrich for particular mutations that were shown in previous work to enhance binding specificity toward MT1-MMP as single mutations, resulting in a focused library of 68,480 unique sequences. The dual-target screening strategy employed was similar to that described above, involving simultaneous positive and negative selection with the two targets, each labeled with a different fluorophore. Library screening resulted in identification of distinct specificity signatures for N-TIMP-2 binding to MMP-9 and MT1-MMP. Hallmarks of MMP-9 specificity included Pro at position 4, Asn at position 38, and Ile at position 71, while hallmarks of MT1-MMP specificity included Gln at position 38 and Asn at position 71. The optimized MMP-9-selective N-TIMP-2 variant demonstrated a significant 1000-fold preference towards MMP-9 compared to MT1-MMP. As in the prior studies, computationally modeled structures contributed insights into the mechanisms responsible for conferring the observed shifts in selectivity (Shirian et al., 2018).
Beyond screening combinatorial libraries to extract a few top-performing variants for a particular target, YSD screening can also be combined with computational analyses to comprehensively map binding specificity landscapes of PPIs. For example, Aharon et al. (2020) developed an innovative approach for quantitatively mapping landscapes for interactions between N-TIMP-2 and three MMP targets (MMP-1, MMP-3, and MT1-MMP). The approach employed an N-TIMP-2 YSD library containing all possible single mutations at seven specificity-determining positions of the TIMP/MMP interface (N-TIMP-2 positions 4, 35, 38, 68, 71, 97, and 99). Multi-target library screening into high- and low-affinity populations was combined with advanced NGS analysis to build predictive models capable of comprehensively identifying hot-spots, cold-spots, and specificity-switch mutations that shape affinity and specificity (Aharon et al., 2020). The quantitative binding landscapes describing the interactions between N-TIMP-2 and these three MMPs showed surprising differences, despite the homologous structures of the targets. For example, the binding landscapes for N-TIMP-2/MMP-3 and N-TIMP-2/MMP-1 demonstrated the PPIs to be nearly fully optimized, with the majority of single mutations causing an affinity loss. For N-TIMP-2/MT1-MMP, by contrast, the native PPI was not tuned for optimal affinity, as evidenced by the wide range of binding affinities among mutants, where many mutations had the potential to increase binding affinity (Aharon et al., 2020). As illustrated by this study, integrating post-screening computational analyses with the YSD platform offers a compelling advance, enabling a more holistic understanding of PPIs that may in the future inform the design of TIMPs with yet greater potency and selectivity for specific MMP targets.
Recent efforts to engineer the TIMP scaffold have gone beyond simple residue substitutions, incorporating larger changes into the protein structure. For example, Bonadio et al. (2023) used computational design to incorporate a loop extension into N-TIMP-2, to facilitate new interactions with a nonconserved region of the MMP surface. This design, in which a 7-residue insertion extended the CD-loop, was then used as the basis for generating a focused combinatorial library of N-TIMP-2 variants. In the library, two residues predicted to interact directly with the nonconserved MMP surface patch were fully randomized, along with two loop-neighboring residues, while the five remaining residues of the engineered loop were partially diversified. YSD screening of this library for high affinity to the target MT1-MMP and low affinity to the anti-target MMP-3 was followed by NGS analysis to identify promising variants. The most optimized N-TIMP-2 variant showed Ki of 29 pM for MT1-MMP and 2.4 μM for MMP-3, a specificity improvement of 7500-fold compared to wild-type N-TIMP-2. By incorporating NGS data into an AlphaFold multiple sequence alignment, high-confidence structural models were generated and then validated via experimental mutagenesis, revealing how the loop extension interacts and discriminates between non-conserved residues on MT1-MMP and MMP-3. This study highlights how computational design and structural modeling can synergize with YSD screening of combinatorial libraries to reshape the basic structure of TIMPs into yet more selective inhibitors of individual MMPs.
Another recent direction in TIMP engineering has been toward the development of bivalent inhibitors capable of simultaneously targeting both an MMP and a second target for enhanced biological effects. In one such study, Yosef et al. (2018) centered their attention on simultaneously targeting synergistic roles of MT1-MMP and integrin αvβ3 in cancer progression. These molecules collaborate to foster invasion, metastasis, and angiogenesis, in part via a sequential activation pathway involving pro-MMP-2 activation by MT1-MMP and subsequent localization at the cell surface by αvβ3 integrin. This work built upon an engineered MT1-MMP-selective N-TIMP-2 variant previously developed by the group (Arkadash et al., 2017), and next used combinatorial library design and YSD screening to evolve a second epitope on N-TIMP2, distal to the MMP interaction site, to selectively bind to αvβ3 integrin. The apex of the N-TIMP-2 B-C loop is known to naturally bind to integrin α3β1 (Seo et al., 2011); here, this epitope was replaced in the combinatorial library with the diversified sequence X-X-X-R-G-D-X-X-X. Screening of the YSD library enabled identification of N-TIMP-2 variants with high affinity for integrin αvβ3 (Yosef et al., 2018). Monovalent and bivalent fusion constructs featuring one or two engineered N-TIMP-2 domains were evaluated in cell-based models of glioblastoma and in vitro angiogenesis, to select the top-performing agent, a bi-specific heterodimer targeting both MT1-MMP and αvβ3 integrin. This agent effectively curbed pro-MMP2 activation, invasiveness and capillary tube formation in glioblastoma and endothelial cell models, and reduced tumor growth in nude mice. The assays underscored the benefits of concurrent intervention, while the innovative engineering approach highlights the potential of the TIMP scaffold for bi-specific targeting, not limited to its potent native inhibition of MMP activity.
Another study from Yosef et al. (2021) aimed to produce a multi-specific agent capable of simultaneously intervening in oncogenic processes promoted by MMP-9 and its cell surface receptor CD44. These molecules have complex interactions crucial for cancer progression. MMP-9 not only dimerizes and colocalizes on the cell surface with CD44 through interactions of its hemopexin domain (PEX), but also cleaves CD44 to stimulate cell signaling and further induce MMP-9 expression. Together this interplay boosts cancer cell migration, invasiveness, and facilitates ECM degradation and metastatic potential. Yosef et al. (2021) developed a multi-specific inhibitor that targets both the PEX and catalytic domains of MMP-9 as well as the interaction with CD44 and its downstream oncogenic effects. First, they used the computationally designed N-TIMP-2 library and YSD screening strategy previously described by Arkadash et al. (2017) to identify a high-affinity inhibitor with Ki of 14 pM for inhibition of MMP-9. Then, they conjugated this inhibitor with PEX using a flexible linker to produce a multi-specific inhibitor that simultaneously targets the MMP9 catalytic domain, the PEX dimerization epitope, as well as CD44. In cell-based assays, the chimeric inhibitor significantly reduced MMP-9 secretion, suppressed MMP-9 gelatinolytic activity, and inhibited the ERK 1/2 phosphorylation that results from interaction of MMP-9 homodimers with CD44 (Yosef et al., 2021). This multi-specific inhibitor represents an advance in TIMP-based antagonists by simultaneously targeting multiple interactions to amplify the potential to disrupt oncogenic functions.
6 Conclusion
There has been considerable effort put into generating potent selective inhibitors of disease promoting MMPs, most of which target the catalytic sites of these enzymes (Fields, 2019). Progress has nonetheless been limited, however, due to the challenging constraint of very close structural homology among the 23 human MMP active sites. Three decades after YSD technology emerged, the method has transformed the investigation, engineering, and interpretation of a wide range of protein-protein interactions. Here, we have highlighted recent advancements in which protein engineering of TIMPs through methods combining YSD with structural and computational approaches have enabled creation of TIMP-based MMP inhibitors with enhanced affinity and specificity. The N-TIMP-2 domain has proven to possess a malleable scaffold, tolerant of mutations, for reengineering with altered protease specificity; it is also shown to be a versatile building block for assembly of bivalent and multivalent agents. While the C-terminal domain has been overlooked in most studies, this domain can contribute up to a third of the contact surface area at the TIMP/MMP interface (Batra et al., 2013). The work of Raeeszadeh-Sarmazdeh and colleagues has demonstrated how this interface can be synergistically optimized through structure-based library design and YSD screening to harness the potential of multidomain TIMPs for high-affinity, selective MMP inhibition (Raeeszadeh-Sarmazdeh et al., 2019; Raeeszadeh-Sarmazdeh et al., 2022). We anticipate that further developments in these approaches will continue to expand our understanding of TIMP/MMP specificity landscapes. One specific opportunity includes incorporating more sophisticated analyses of sequence and affinity enrichment data obtained from YSD screens, by harnessing state-of-the-art machine learning methods. Such methods may allow for accurate affinity predictions for complex combinations of mutations not directly observed within a library, thus expanding the size of the sequence space that can be effectively mapped. Another area for development includes the more detailed mapping of interactions of the four native TIMPs with effectors and binding partners beyond the MMP catalytic domains. A better structural understanding of these natural interactions will aid in the targeted ablation of unwanted off-target effects of engineered TIMP-based MMP inhibitors. These advances will ultimately enable creation of a toolbox of selective TIMP probes and inhibitors with therapeutic potential for a wide variety of diseases.
Author contributions
AS: Writing–original draft, Writing–review and editing. GK-T: Writing–original draft, Writing–review and editing. MC: Validation, Writing–review and editing. ER: Funding acquisition, Supervision, Writing–original draft, Writing–review and editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. We acknowledge support to ER from the U.S. National Institutes of Health grants R01GM132100, R01GM144393, R01CA258274, and R01HL157424 and U.S.-Israel Binational Science Foundation grant 2019303.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
References
Aharon, L., Aharoni, S. L., Radisky, E. S., and Papo, N. (2020). Quantitative mapping of binding specificity landscapes for homologous targets by using a high-throughput method. Biochem. J. 477 (9), 1701–1719. doi:10.1042/BCJ20200188
Amar, S., Smith, L., and Fields, G. B. (2017). Matrix metalloproteinase collagenolysis in health and disease. Biochim. Biophys. Acta Mol. Cell Res. 1864, 1940–1951. doi:10.1016/j.bbamcr.2017.04.015
Arkadash, V., Radisky, E. S., and Papo, N. (2018). Combinatorial engineering of N-TIMP2 variants that selectively inhibit MMP9 and MMP14 function in the cell. Oncotarget 9 (62), 32036–32053. doi:10.18632/oncotarget.25885
Arkadash, V., Yosef, G., Shirian, J., Cohen, I., Horev, Y., Grossman, M., et al. (2017). Development of high affinity and high specificity inhibitors of matrix metalloproteinase 14 through computational design and directed evolution. J. Biol. Chem. 292 (8), 3481–3495. doi:10.1074/jbc.M116.756718
Batra, J., and Radisky, E. S. (2011). Tissue inhibitors of metalloproteinases (TIMPs): inhibition of Zn-dependent metallopeptidases. Encycl. Inorg. Bioinorganic Chem. 2011, 1–10. doi:10.1002/9781119951438.eibc2168
Batra, J., Soares, A. S., Mehner, C., and Radisky, E. S. (2013). Matrix metalloproteinase-10/TIMP-2 structure and analyses define conserved core interactions and diverse exosite interactions in MMP/TIMP complexes. PLoS one 8 (9), e75836. doi:10.1371/journal.pone.0075836
Bissett, D., O'byrne, K. J., Von Pawel, J., Gatzemeier, U., Price, A., Nicolson, M., et al. (2005). Phase III study of matrix metalloproteinase inhibitor prinomastat in non-small-cell lung cancer. J. Clin. Oncol. 23 (4), 842–849. doi:10.1200/JCO.2005.03.170
Bode, W., Gomis-Rüth, F. X., and Stöckler, W. (1993). Astacins, serralysins, snake venom and matrix metalloproteinases exhibit identical zinc-binding environments (HEXXHXXGXXH and Met-turn) and topologies and should be grouped into a common family, the ‘metzincins. FEBS Lett. 331 (1-2), 134–140. doi:10.1016/0014-5793(93)80312-i
Boder, E. T., and Wittrup, K. D. (1997). Yeast surface display for screening combinatorial polypeptide libraries. Nat. Biotechnol. 15 (6), 553–557. doi:10.1038/nbt0697-553
Bonadio, A., and Shifman, J. M. (2021). Computational design and experimental optimization of protein binders with prospects for biomedical applications. Protein Eng. Des. Sel. 34, gzab020. doi:10.1093/protein/gzab020
Bonadio, A., Wenig, B. L., Hockla, A., Radisky, E. S., and Shifman, J. M. (2023). Designed loop extension followed by combinatorial screening confers high specificity to a broad matrix MetalloproteinaseInhibitor. J. Mol. Biol. 435 (13), 168095. doi:10.1016/j.jmb.2023.168095
Brew, K., and Nagase, H. (2010). The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity. Biochim. Biophys. Acta 1803 (1), 55–71. doi:10.1016/j.bbamcr.2010.01.003
Cerda-Costa, N., and Gomis-Ruth, F. X. (2014). Architecture and function of metallopeptidase catalytic domains. Protein Sci. 23 (2), 123–144. doi:10.1002/pro.2400
Cherf, G. M., and Cochran, J. R. (2015). Applications of yeast surface display for protein engineering. Methods Mol. Biol. 1319, 155–175. doi:10.1007/978-1-4939-2748-7_8
Clements, J. M., Cossins, J. A., Wells, G. M., Corkill, D. J., Helfrich, K., Wood, L. M., et al. (1997). Matrix metalloproteinase expression during experimental autoimmune encephalomyelitis and effects of a combined matrix metalloproteinase and tumour necrosis factor-alpha inhibitor. J. Neuroimmunol. 74, 85–94. doi:10.1016/s0165-5728(96)00210-x
Elkins, P. A., Ho, Y. S., Smith, W. W., Janson, C. A., D'alessio, K. J., Mcqueney, M. S., et al. (2002). Structure of the C-terminally truncated human ProMMP9, a gelatin-binding matrix metalloproteinase. Acta Crystallogr. D. Biol. Crystallogr. 58, 1182–1192. doi:10.1107/s0907444902007849
Fields, G. B. (2019). The rebirth of matrix metalloproteinase inhibitors: moving beyond the dogma. Cells 8 (9), 984. doi:10.3390/cells8090984
Fischer, T., Senn, N., and Riedl, R. (2019). Design and structural evolution of matrix metalloproteinase inhibitors. Chemistry 25 (34), 7960–7980. doi:10.1002/chem.201805361
Gai, S. A., and Wittrup, K. D. (2007). Yeast surface display for protein engineering and characterization. Curr. Opin. Struct. Biol. 17 (4), 467–473. doi:10.1016/j.sbi.2007.08.012
Gialeli, C., Theocharis, A. D., and Karamanos, N. K. (2011). Roles of matrix metalloproteinases in cancer progression and their pharmacological targeting. FEBS J. 278 (1), 16–27. doi:10.1111/j.1742-4658.2010.07919.x
Gomis-Ruth, F. X., Maskos, K., Betz, M., Bergner, A., Huber, R., Suzuki, K., et al. (1997). Mechanism of inhibition of the human matrix metalloproteinase stromelysin-1 by TIMP-1. Nature 389 (6646), 77–81. doi:10.1038/37995
Grillet, B., Pereira, R. V. S., Van Damme, J., Abu El-Asrar, A., Proost, P., and Opdenakker, G. (2023). Matrix metalloproteinases in arthritis: towards precision medicine. Nat. Rev. Rheumatol. 19 (6), 363–377. doi:10.1038/s41584-023-00966-w
Hu, J., Van Den Steen, P. E., Sang, Q. X., and Opdenakker, G. (2007). Matrix metalloproteinase inhibitors as therapy for inflammatory and vascular diseases. Nat. Rev. Drug Discov. 6 (6), 480–498. doi:10.1038/nrd2308
Itoh, Y. (2015). Membrane-type matrix metalloproteinases: their functions and regulations. Matrix Biol. 44, 207–223. doi:10.1016/j.matbio.2015.03.004
Klein, T., and Bischoff, R. (2011). Physiology and pathophysiology of matrix metalloproteases. Amino acids 41 (2), 271–290. doi:10.1007/s00726-010-0689-x
Krzeski, P., Buckland-Wright, C., Balint, G., Cline, G. A., Stoner, K., Lyon, R., et al. (2007). Development of musculoskeletal toxicity without clear benefit after administration of PG-116800, a matrix metalloproteinase inhibitor, to patients with knee osteoarthritis: a randomized, 12-month, double-blind, placebo-controlled study. Arthritis Res. Ther. 9 (5), R109. doi:10.1186/ar2315
Levin, M., Udi, Y., Solomonov, I., and Sagi, I. (2017). Next generation matrix metalloproteinase inhibitors - novel strategies bring new prospects. Biochim. Biophys. Acta Mol. Cell Res. 1864, 1927–1939. doi:10.1016/j.bbamcr.2017.06.009
Libson, A. M., Gittis, A. G., Collier, I. E., Marmer, B. L., Goldberg, G. I., and Lattman, E. E. (1995). Crystal structure of the haemopexin-like C-terminal domain of gelatinase A. Nat. Struct. Biol. 2 (11), 938–942. doi:10.1038/nsb1195-938
Lovejoy, B., Cleasby, A., Hassell, A. M., Longley, K., Luther, M. A., Weigl, D., et al. (1994). Structure of the catalytic domain of fibroblast collagenase complexed with an inhibitor. Science 263 (5145), 375–377. doi:10.1126/science.8278810
Malemud, C. J. (2006). Matrix metalloproteinases (MMPs) in health and disease: an overview. Front. Biosci. 11, 1696–1701. doi:10.2741/1915
Maskos, K. (2005). Crystal structures of MMPs in complex with physiological and pharmacological inhibitors. Biochimie 87 (3-4), 249–263. doi:10.1016/j.biochi.2004.11.019
Massova, I., Kotra, L. P., Fridman, R., and Mobashery, S. (1998). Matrix metalloproteinases: structures, evolution, and diversification. FASEB J. 12 (12), 1075–1095. doi:10.1096/fasebj.12.12.1075
Moore, M. J., Hamm, J., Dancey, J., Eisenberg, P. D., Dagenais, M., Fields, A., et al. (2003). Comparison of gemcitabine versus the matrix metalloproteinase inhibitor BAY 12-9566 in patients with advanced or metastatic adenocarcinoma of the pancreas: a phase III trial of the National Cancer Institute of Canada Clinical Trials Group. J. Clin. Oncol. 21 (17), 3296–3302. doi:10.1200/JCO.2003.02.098
Morgunova, E., Tuuttila, A., Bergmann, U., and Tryggvason, K. (2002). Structural insight into the complex formation of latent matrix metalloproteinase 2 with tissue inhibitor of metalloproteinase 2. Proc. Natl. Acad. Sci. U. S. A. 99 (11), 7414–7419. doi:10.1073/pnas.102185399
Murphy, G. (2011). Tissue inhibitors of metalloproteinases. Genome Biol. 12 (11), 233. doi:10.1186/gb-2011-12-11-233
Murphy, G., and Nagase, H. (2008). Progress in matrix metalloproteinase research. Mol. Asp. Med. 29 (5), 290–308. doi:10.1016/j.mam.2008.05.002
Nagase, H. (1996). Matrix metalloproteinases. Zinc metalloproteases health Dis. 1996, 173–224. doi:10.1201/9781482272765
Nagase, H., and Brew, K. (2002). Engineering of tissue inhibitor of metalloproteinases mutants as potential therapeutics. Arthritis Res. 4 (Suppl. 3), S51–S61. doi:10.1186/ar573
Nagase, H., Visse, R., and Murphy, G. (2006). Structure and function of matrix metalloproteinases and TIMPs. Cardiovasc Res. 69 (3), 562–573. doi:10.1016/j.cardiores.2005.12.002
Overall, C. M. (2002). Molecular determinants of metalloproteinase substrate specificity: matrix metalloproteinase substrate binding domains, modules, and exosites. Mol. Biotechnol. 22 (1), 51–86. doi:10.1385/MB:22:1:051
Overall, C. M., and Kleifeld, O. (2006). Tumour microenvironment - opinion: validating matrix metalloproteinases as drug targets and anti-targets for cancer therapy. Nat. Rev. Cancer 6 (3), 227–239. doi:10.1038/nrc1821
Park, S., Xu, Y., Stowell, X. F., Gai, F., Saven, J. G., and Boder, E. T. (2006). Limitations of yeast surface display in engineering proteins of high thermostability. Protein Eng. Des. Sel. 19 (5), 211–217. doi:10.1093/protein/gzl003
Ra, H. J., and Parks, W. C. (2007). Control of matrix metalloproteinase catalytic activity. Matrix Biol. 26 (8), 587–596. doi:10.1016/j.matbio.2007.07.001
Radisky, E. S., and Coban, M. (2021). “Enzymes matrix metalloproteinases,” in Encyclopedia of biological chemistry III. Editor J Jez. Third Edition (Oxford: Elsevier), 336–353.
Radisky, E. S., Raeeszadeh-Sarmazdeh, M., and Radisky, D. C. (2017). Therapeutic potential of matrix metalloproteinase inhibition in breast cancer. J. Cell Biochem. 118 (11), 3531–3548. doi:10.1002/jcb.26185
Raeeszadeh-Sarmazdeh, M., Coban, M., Mahajan, S., Hockla, A., and Sankaran, B., (2022). Engineering of tissue inhibitor of metalloproteinases TIMP-1 for fine discrimination between closely related stromelysins MMP-3 and MMP-10. J. Biol. Chem. 298 (3), 101654. doi:10.1016/j.jbc.2022.101654
Raeeszadeh-Sarmazdeh, M., Greene, K. A., Sankaran, B., Downey, G. P., Radisky, D. C., and Radisky, E. S. (2019). Directed evolution of the metalloproteinase inhibitor TIMP-1 reveals that its N-and C-terminal domains cooperate in matrix metalloproteinase recognition. J. Biol. Chem. 294 (24), 9476–9488. doi:10.1074/jbc.RA119.008321
Rosenfeld, L., Heyne, M., Shifman, J. M., and Papo, N. (2016). Protein engineering by combined computational and in vitro evolution approaches. Trends Biochem. Sci. 41 (5), 421–433. doi:10.1016/j.tibs.2016.03.002
Seo, D. W., Saxinger, W. C., Guedez, L., Cantelmo, A. R., Albini, A., and Stetler-Stevenson, W. G. (2011). An integrin-binding N-terminal peptide region of TIMP-2 retains potent angio-inhibitory and anti-tumorigenic activity in vivo. Peptides 32 (9), 1840–1848. doi:10.1016/j.peptides.2011.08.010
Sharabi, O., Shirian, J., Grossman, M., Lebendiker, M., Sagi, I., and Shifman, J. (2014). Affinity- and specificity-enhancing mutations are frequent in multispecific interactions between TIMP2 and MMPs. PLoS One 9 (4), e93712. doi:10.1371/journal.pone.0093712
Shirian, J., Arkadash, V., Cohen, I., Sapir, T., Radisky, E. S., Papo, N., et al. (2018). Converting a broad matrix metalloproteinase family inhibitor into a specific inhibitor of MMP-9 and MMP-14. FEBS Lett. 592 (7), 1122–1134. doi:10.1002/1873-3468.13016
Sparano, J. A., Bernardo, P., Stephenson, P., Gradishar, W. J., Ingle, J. N., Zucker, S., et al. (2004). Randomized phase III trial of marimastat versus placebo in patients with metastatic breast cancer who have responding or stable disease after first-line chemotherapy: eastern Cooperative Oncology Group trial E2196. J. Clin. Oncol. 22, 4683–4690. doi:10.1200/JCO.2004.08.054
Sternlicht, M. D., and Werb, Z. (2001). How matrix metalloproteinases regulate cell behavior. Annu. Rev. Cell Dev. Biol. 17, 463–516. doi:10.1146/annurev.cellbio.17.1.463
Tallant, C., Marrero, A., and Gomis-Ruth, F. X. (2010). Matrix metalloproteinases: fold and function of their catalytic domains. Biochim. Biophys. Acta 1803 (1), 20–28. doi:10.1016/j.bbamcr.2009.04.003
Visse, R., and Nagase, H. (2003). Matrix metalloproteinases and tissue inhibitors of metalloproteinases: structure, function, and biochemistry. Circ. Res. 92 (8), 827–839. doi:10.1161/01.RES.0000070112.80711.3D
Wilson, C. L., Ouellette, A. J., Satchell, D. P., Ayabe, T., Lopez-Boado, Y. S., Stratman, J. L., et al. (1999). Regulation of intestinal alpha-defensin activation by the metalloproteinase matrilysin in innate host defense. Science 286 (5437), 113–117. doi:10.1126/science.286.5437.113
Yosef, G., Arkadash, V., and Papo, N. (2018). Targeting the MMP-14/MMP-2/integrin αvβ3 axis with multispecific N-TIMP2-based antagonists for cancer therapy. J. Biol. Chem. 293 (34), 13310–13326. doi:10.1074/jbc.RA118.004406
Yosef, G., Hayun, H., and Papo, N. (2021). Simultaneous targeting of CD44 and MMP9 catalytic and hemopexin domains as a therapeutic strategy. Biochem. J. 478 (5), 1139–1157. doi:10.1042/BCJ20200628
Yu, W. H., Yu, S., Meng, Q., Brew, K., and Woessner, J. F. (2000). TIMP-3 binds to sulfated glycosaminoglycans of the extracellular matrix. J. Biol. Chem. 275 (40), 31226–31232. doi:10.1074/jbc.M000907200
Keywords: matrix metalloproteinase, tissue inhibitor of metalloproteinases, yeast surface display, protein engineering, directed evolution
Citation: Shoari A, Khalili-Tanha G, Coban MA and Radisky ES (2023) Structure and computation-guided yeast surface display for the evolution of TIMP-based matrix metalloproteinase inhibitors. Front. Mol. Biosci. 10:1321956. doi: 10.3389/fmolb.2023.1321956
Received: 15 October 2023; Accepted: 13 November 2023;
Published: 23 November 2023.
Edited by:
Alessandro Angelini, Ca’ Foscari University of Venice, ItalyReviewed by:
Athanasios Papakyriakou, National Centre of Scientific Research Demokritos, GreeceCopyright © 2023 Shoari, Khalili-Tanha, Coban and Radisky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Evette S. Radisky, cmFkaXNreS5ldmV0dGVAbWF5by5lZHU=