Increased backbone mobility in beta-barrel enhances entropy gain driving binding of N-TIMP-1 to MMP-3

Designed loop extension followed by combinatorial screening confers high specificity to a broad matrix metalloproteinaseinhibitor

Matrix metalloproteinases (MMPs) are key drivers of various diseases, including cancer. Development of selective probes and drugs capable of selectively inhibiting the individual members of the large MMP family remains a persistent challenge. The inhibitory N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP2), a natural broad MMP inhibitor, can provide a scaffold for protein engineering to create more selective MMP inhibitors. Here, we pursued a unique approach harnessing both computational design and combinatorial screening to confer high binding specificity toward a target MMP in preference to an anti-target MMP. We designed a loop extension of N-TIMP2 to allow new interactions with the non-conserved MMP surface and generated an efficient focused library for yeast surface display, which was then screened for high binding to the target MMP-14 and low binding to anti-target MMP-3. Deep sequencing analysis identified the most promising variants, which were expressed, purified, and tested for selectivity of inhibition. Our best N-TIMP2 variant exhibited 29 pM binding affinity to MMP-14 and 2.4 µM affinity to MMP-3, revealing 7500-fold greater specificity than WT N-TIMP2. High-confidence structural models were obtained by including NGS data in the AlphaFold multiple sequence alignment. The modeling together with experimental mutagenesis validated our design predictions, demonstrating that the loop extension packs tightly against non-conserved residues on MMP-14 and clashes with MMP-3. This study demonstrates how introduction of loop extensions in a manner guided by target protein conservation data and loop design can offer an attractive strategy to achieve specificity in design of protein ligands.

Transient association between proteins elicits alteration of dynamics at sites far away from interfaces

Proteins are known to undergo structural changes upon binding to partner proteins. However, the prevalence, extent, location, and function of change in protein dynamics due to transient protein-protein interactions is not well documented. Here, we have analyzed a dataset of 58 protein-protein complexes of known three-dimensional structure and structures of their corresponding unbound forms to evaluate dynamics changes induced by binding. Fifty-five percent of cases showed significant dynamics change away from the interfaces. This change is not always accompanied by an observed structural change. Binding of protein partner is found to alter inter-residue communication within the tertiary structure in about 90% of cases. Also, residue motions accessible to proteins in unbound form were not always maintained in the bound form. Further analyses revealed functional roles for the distant site where dynamics change was observed. Overall, the results presented here strongly suggest that alteration of protein dynamics due to binding of a partner protein commonly occurs.

Thermodynamic profiles of the interactions of suramin, chondroitin sulfate, and pentosan polysulfate with the inhibitory domain of TIMP-3

Extracellular levels of soluble TIMP-3 are low, reflecting its binding by extracellular matrix (ECM) components including sulfated glycosaminoglycans (SGAGs) and endocytosis via low density lipoprotein receptor-related protein 1. Since TIMP-3 inhibits ECM degradation, the ability of SGAGs to elevate extracellular TIMP-3 is significant for osteoarthritis treatment. Previous studies of such interactions have utilized immobilized TIMP-3 or ligands. Here, we report the thermodynamics of the interactions of the sGAG-binding N-domain of TIMP-3 with chondroitin sulfate, pentosan polysulfate, and suramin in solution using isothermal titration calorimetry. All three interactions are driven by a favorable negative enthalpy change combined with an unfavorable decrease in entropy. The heat capacity changes (ΔC_p ) for all of the interactions are zero, indicating an insignificant contribution from hydrophobic interactions.

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.

Therapeutic Potential of Matrix Metalloproteinase Inhibition in Breast Cancer

Matrix metalloproteinases (MMPs) are a family of zinc endopeptidases that cleave nearly all components of the extracellular matrix as well as many other soluble and cell-associated proteins. MMPs have been implicated in normal physiological processes, including development, and in the acquisition and progression of the malignant phenotype. Disappointing results from a series of clinical trials testing small molecule, broad spectrum MMP inhibitors as cancer therapeutics led to a re-evaluation of how MMPs function in the tumor microenvironment, and ongoing research continues to reveal that these proteins play complex roles in cancer development and progression. It is now clear that effective targeting of MMPs for therapeutic benefit will require selective inhibition of specific MMPs. Here, we provide an overview of the MMP family and its biological regulators, the tissue inhibitors of metalloproteinases (TIMPs). We then summarize recent research from model systems that elucidate how specific MMPs drive the malignant phenotype of breast cancer cells, including acquisition of cancer stem cell features and induction of the epithelial-mesenchymal transition, and we also outline clinical studies that implicate specific MMPs in breast cancer outcomes. We conclude by discussing ongoing strategies for development of inhibitors with therapeutic potential that are capable of selectively targeting the MMPs most responsible for tumor promotion, with special consideration of the potential of biologics including antibodies and engineered proteins based on the TIMP scaffold. J. Cell. Biochem. 118: 3531-3548, 2017. © 2017 The Authors. Journal of Cellular Biochemistry Published by Wiley Periodicals, Inc.

Advances in NMR Methods To Map Allosteric Sites: From Models to Translation

The last five years have witnessed major developments in the understanding of the allosteric phenomenon, broadly defined as coupling between remote molecular sites. Such advances have been driven not only by new theoretical models and pharmacological applications of allostery, but also by progress in the experimental approaches designed to map allosteric sites and transitions. Among these techniques, NMR spectroscopy has played a major role given its unique near-atomic resolution and sensitivity to the dynamics that underlie allosteric couplings. Here, we highlight recent progress in the NMR methods tailored to investigate allostery with the goal of offering an overview of which NMR approaches are best suited for which allosterically relevant questions. The picture of the allosteric "NMR toolbox" is provided starting from one of the simplest models of allostery (i.e., the four-state thermodynamic cycle) and continuing to more complex multistate mechanisms. We also review how such an "NMR toolbox" has assisted the elucidation of the allosteric molecular basis for disease-related mutations and the discovery of novel leads for allosteric drugs. From this overview, it is clear that NMR plays a central role not only in experimentally validating transformative theories of allostery, but also in tapping the full translational potential of allosteric systems.

Thermodynamic Basis of Selectivity in the Interactions of Tissue Inhibitors of Metalloproteinases N-domains with Matrix Metalloproteinases-1, -3, and -14

The four tissue inhibitors of metalloproteinases (TIMPs) are potent inhibitors of the many matrixins (MMPs), except that TIMP1 weakly inhibits some MMPs, including MMP14. The broad-spectrum inhibition of MMPs by TIMPs and their N-domains (NTIMPs) is consistent with the previous isothermal titration calorimetric finding that their interactions are entropy-driven but differ in contributions from solvent and conformational entropy (ΔSsolv, ΔSconf), estimated using heat capacity changes (ΔCp). Selective engineered NTIMPs have potential applications for treating MMP-related diseases, including cancer and cardiomyopathy. Here we report isothermal titration calorimetric studies of the effects of selectivity-modifying mutations in NTIMP1 and NTIMP2 on the thermodynamics of their interactions with MMP1, MMP3, and MMP14. The weak inhibition of MMP14 by NTIMP1 reflects a large conformational entropy penalty for binding. The T98L mutation, peripheral to the NTIMP1 reactive site, enhances binding by increasing ΔSsolv but also reduces ΔSconf However, the same mutation increases NTIMP1 binding to MMP3 in an interaction that has an unusual positive ΔCp This indicates a decrease in solvent entropy compensated by increased conformational entropy, possibly reflecting interactions involving alternative conformers. The NTIMP2 mutant, S2D/S4A is a selective MMP1 inhibitor through electrostatic effects of a unique MMP-1 arginine. Asp-2 increases reactive site polarity, reducing ΔCp, but increases conformational entropy to maintain strong binding to MMP1. There is a strong negative correlation between ΔSsolv and ΔSconf for all characterized interactions, but the data for each MMP have characteristic ranges, reflecting intrinsic differences in the structures and dynamics of their free and inhibitor-bound forms.

Energetic and flexibility properties captured by long molecular dynamics simulations of a membrane-embedded pMHCII-TCR complex

Although crystallographic data have provided important molecular insight into the interactions in the pMHC-TCR complex, the inherent features of this structural approach cause it to only provide a static picture of the interactions. While unbiased molecular dynamics simulations (UMDSs) have provided important information about the dynamic structural behavior of the pMHC-TCR complex, most of them have modeled the pMHC-TCR complex as soluble, when in physiological conditions, this complex is membrane bound; therefore, following this latter UMDS protocol might hamper important dynamic results. In this contribution, we performed three independent 300 ns-long UMDSs of the pMHCII-TCR complex anchored in two opposing membranes to explore the structural and energetic properties of the recognition of pMHCII by the TCR. The conformational ensemble generated through UMDSs was subjected to clustering and Cartesian principal component analyses (cPCA) to explore the dynamical behavior of the pMHCII-TCR association. Furthermore, based on the conformational population sampled through UMDSs, the effective binding free energy, per-residue free energy decomposition, and alanine scanning mutations were explored for the native pMHCII-TCR complex, as well as for 12 mutations (p1-p12MHCII-TCR) introduced in the native peptide. Clustering analyses and cPCA provide insight into the rocking motion of the TCR onto pMHCII, together with the presence of new electrostatic interactions not observed through crystallographic methods. Energetic results provide evidence of the main contributors to the pMHC-TCR complex formation as well as the key residues involved in this molecular recognition process.

Fluctuations within folded proteins: implications for thermodynamic and allosteric regulation

Folded protein structures are both stable and dynamic. Historically, our clearest window into these structures came from X-ray crystallography, which generally provided a static image of each protein's singular "folded state", highlighting its stability. Deviations away from that crystallographic structure were difficult to quantify, and as a result, their potential functional consequences were often neglected. However, several dynamical and statistical studies now highlight the structural variability that is present within the protein's folded state. Here we review mounting evidence of the importance of these structural rearrangements; both experiment and computation indicate that folded proteins undergo substantial fluctuations that can greatly influence their function. Crucially, recent studies have shown that structural elements of proteins, especially their side-chain degrees of freedom, fluctuate in ways that generate significant conformational heterogeneity. The entropy associated with these motions contributes to the folded structure's thermodynamic stability. In addition, since these fluctuations can shift in response to perturbations such as ligand binding, they may play an important role in the protein's capacity to respond to environmental cues. In one compelling example, the entropy associated with side-chain fluctuations contributes significantly to regulating the binding of calmodulin to a set of peptide ligands. The neglect of fluctuations within proteins' native states was often justified by the dense packing within folded proteins, which has inspired comparisons with crystalline solids. Many liquids, however, can achieve similarly dense packing yet fluidity is maintained through correlated molecular motions. Indeed, the studies we discuss favor comparison of folded proteins not with solids but instead with dense liquids, where the internal side chain fluidity is facilitated by collective motions that are correlated over long distances. These correlated rearrangements can enable allosteric communication between different parts of a protein, through subtle and varied channels. Such long-range correlations appear to be an innate feature of proteins in general, manifest even in molecules lacking known allosteric regulators and arising robustly from the physical nature of their internal environment. Given their ubiquity, it is only to be expected that, over time, nature has refined some subset of these correlated motions and put them to use. Native state fluctuations increasingly appear to be vital for proteins' natural functions. Understanding the diversity, origin, and range of these rearrangements may provide novel routes for rationally manipulating biomolecular activity.

TNF-α stimulates MMP-3 production via PGE2 signalling through the NF-kB and p38 MAPK pathway in a murine cementoblast cell line

BACKGROUND

Cementoblasts are considered to play an important role in the homeostasis of periodontal tissues under both physiologic and pathologic conditions. Matrix metalloproteinases (MMPs) is the key family of enzymes participating in extracellular matrix remodelling. In the present study, the effects and regulatory mechanisms of tumour necrosis factor (TNF)-α on the expression of MMPs and their inhibitors (tissue inhibitor of metalloproteinases; TIMPs) were investigated.

MATERIALS AND METHODS

OCCM-30, an immortalised murine cementoblast cell line, was stimulated with TNF-α at 1 and 10ng/ml for 24h. The expression of Mmp-2, Mmp-3, Mmp-13, Mmp-14, Timp-1, and Timp-2 as well as PGE2 was determined. Inhibitors of MAPKs, PI3K/Akt, NF-kB and Cox-2 were employed to reveal possible TNF-α induced regulatory signalling pathway(s). The mRNA and protein expression were analysed by (semi)quantitative real-time PCR and enzyme-linked immunosorbent assay (ELISA), respectively.

RESULTS

TNF-α dose-dependently stimulated MMP-3 expression by cementoblasts. This was found for mRNA as well as protein expression. No significant differences were found in the mRNA expression of Mmp-2, Mmp-13, Mmp-14, Timp-1, and Timp-2 upon TNF-α stimulation. The level of PGE2, however, was significantly increased along with MMP-3. Treatment with a selective Cox-2 inhibitor resulted in partial suppression of TNF-α-induced Mmp-3 mRNA expression. Addition of PGE2 enhanced Mmp-3 mRNA in a dose dependent manner, suggesting an inductive effect of TNF-α partly via PGE2. The up-regulation of Mmp-3 by TNF-α was completely suppressed by a combination of NF-kB and p38 MAPK inhibitors, while partial suppression was found with each inhibitor. The effect of PGE2 on Mmp-3 expression was abolished by treating cells with an NF-kB inhibitor; a p38 MAPK inhibitor had only a small effect.

CONCLUSIONS

The present study indicates that cementoblasts respond to TNF-α by increasing MMP-3 production partially via PGE2 and signalling through the NF-kB and p38 MAPK pathway. MMP-3 may participate in periodontal tissue degradation/remodelling.

NMR structure of the HWE kinase associated response regulator Sma0114 in its activated state

Bacterial receiver domains modulate intracellular responses to external stimuli in two-component systems. Sma0114 is the first structurally characterized representative from the family of receiver domains that are substrates for histidine-tryptophan-glutamate (HWE) kinases. We report the NMR structure of Sma0114 bound by Ca(2+) and BeF3(-), a phosphate analogue that stabilizes the activated state. Differences between the NMR structures of the inactive and activated states occur in helix α1, the active site loop that connects strand β3 and helix α3, and in the segment from strand β5 to helix α5 of the 455 (α4-β5-α5) face. Structural rearrangements of the 455 face typically make receiver domains competent for binding downstream target molecules. In Sma0114 the structural changes accompanying activation result in a more negatively charged surface for the 455 face. Coupling between the 455 face and active site phosphorylation is usually mediated through the rearrangement of a threonine and tyrosine residue, in a mechanism called Y-T coupling. The NMR structure indicates that Sma0114 lacks Y-T coupling and that communication between the active site and the 455 face is achieved through a conserved lysine residue that stabilizes the acyl phosphate in receiver domains. (15)N-NMR relaxation experiments were used to investigate the backbone dynamics of the Sma0114 apoprotein, the binary Sma0114·Ca(2+) complex, and the ternary Sma0114·Ca(2+)·BeF3(-) complex. The loss of entropy due to ligand binding at the active site is compensated by increased flexibility in the 455 face. The dynamic character of the 455 face in Sma0114, which results in part from the replacement of helix α4 by a flexible loop, may facilitate induced-fit recognition of target molecules.

Structural insights into the role of the cyclic backbone in a squash trypsin inhibitor

MCoTI-II is a head-to-tail cyclic peptide with potent trypsin inhibitory activity and, on the basis of its exceptional proteolytic stability, is a valuable template for the design of novel drug leads. Insights into inhibitor dynamics and interactions with biological targets are critical for drug design studies, particularly for protease targets. Here, we show that the cyclization and active site loops of MCoTI-II are flexible in solution, but when bound to trypsin, the active site loop converges to a single well defined conformation. This finding of reduced flexibility on binding is in contrast to a recent study on the homologous peptide MCoTI-I, which suggested that regions of the peptide are more flexible upon binding to trypsin. We provide a possible explanation for this discrepancy based on degradation of the complex over time. Our study also unexpectedly shows that the cyclization loop, not present in acyclic homologues, facilitates potent trypsin inhibitory activity by engaging in direct binding interactions with trypsin.

Intrinsic dynamics of an extended hydrophobic core in the S. cerevisiae RNase III dsRBD contributes to recognition of specific RNA binding sites

The Saccharomyces cerevisiae RNase III enzyme Rnt1p preferentially binds to double-stranded RNA hairpin substrates with a conserved (A/u)GNN tetraloop fold, via shape-specific interactions by its double-stranded RNA-binding domain (dsRBD) helix α1 to the tetraloop minor groove. To investigate whether conformational flexibility in the dsRBD regulates the binding specificity, we determined the backbone dynamics of the Rnt1p dsRBD in the free and AGAA hairpin-bound states using NMR spin-relaxation experiments. The intrinsic microsecond-to-millisecond timescale dynamics of the dsRBD suggests that helix α1 undergoes conformational sampling in the free state, with large dynamics at some residues in the α1-β1 loop (α1-β1 hinge). To correlate free dsRBD dynamics with structural changes upon binding, we determined the solution structure of the free dsRBD used in the previously determined RNA-bound structures. The Rnt1p dsRBD has an extended hydrophobic core comprising helix α1, the α1-β1 loop, and helix α3. Analysis of the backbone dynamics and structures of the free and bound dsRBD reveals that slow-timescale dynamics in the α1-β1 hinge are associated with concerted structural changes in the extended hydrophobic core that govern binding of helix α1 to AGAA tetraloops. The dynamic behavior of the dsRBD bound to a longer AGAA hairpin reveals that dynamics within the hydrophobic core differentiate between specific and nonspecific sites. Mutations of residues in the α1-β1 hinge result in changes to the dsRBD stability and RNA-binding affinity and cause defects in small nucleolar RNA processing invivo. These results reveal that dynamics in the extended hydrophobic core are important for binding site selection by the Rnt1p dsRBD.

Entropy-enthalpy transduction caused by conformational shifts can obscure the forces driving protein-ligand binding

Molecular dynamics simulations of unprecedented duration now can provide new insights into biomolecular mechanisms. Analysis of a 1-ms molecular dynamics simulation of the small protein bovine pancreatic trypsin inhibitor reveals that its main conformations have different thermodynamic profiles and that perturbation of a single geometric variable, such as a torsion angle or interresidue distance, can select for occupancy of one or another conformational state. These results establish the basis for a mechanism that we term entropy-enthalpy transduction (EET), in which the thermodynamic character of a local perturbation, such as enthalpic binding of a small molecule, is camouflaged by the thermodynamics of a global conformational change induced by the perturbation, such as a switch into a high-entropy conformational state. It is noted that EET could occur in many systems, making measured entropies and enthalpies of folding and binding unreliable indicators of actual thermodynamic driving forces. The same mechanism might also account for the high experimental variance of measured enthalpies and entropies relative to free energies in some calorimetric studies. Finally, EET may be the physical mechanism underlying many cases of entropy-enthalpy compensation.

Pentosan polysulfate increases affinity between ADAMTS-5 and TIMP-3 through formation of an electrostatically driven trimolecular complex

The semi-synthetic sulfated polysaccharide PPS (pentosan polysulfate) increases affinity between the aggrecan-degrading ADAMTSs (adamalysins with thrombospondin motifs) and their endogenous inhibitor, TIMP (tissue inhibitor of metalloproteinases)-3. In the present study we demonstrate that PPS mediates the formation of a high-affinity trimolecular complex with ADAMTS-5 and TIMP-3. A TIMP-3 mutant that lacks extracellular-matrix-binding ability was insensitive to this affinity increase, and truncated forms of ADAMTS-5 that lack the Sp (spacer) domain had reduced PPS-binding ability and sensitivity to the affinity increase. PPS molecules composed of 11 or more saccharide units were 100-fold more effective than those of eight saccharide units, indicating the involvement of extended or multiple protein-interaction sites. The formation of a high-affinity trimolecular complex was completely abolished in the presence of 0.4 M NaCl. These results suggest that PPS enhances the affinity between ADAMTS-5 and TIMP-3 by forming electrostatically driven trimolecular complexes under physiological conditions.

Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins

BACKGROUND

Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question "What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?"

RESULTS

A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex.Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors.

CONCLUSIONS

The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.

Conformational dynamics and thermodynamics of protein–ligand binding studied by NMR relaxation

Protein conformational dynamics can be critical for ligand binding in two ways that relate to kinetics and thermodynamics respectively. First, conformational transitions between different substates can control access to the binding site (kinetics). Secondly, differences between free and ligand-bound states in their conformational fluctuations contribute to the entropy of ligand binding (thermodynamics). In the present paper, I focus on the second topic, summarizing our recent results on the role of conformational entropy in ligand binding to Gal3C (the carbohydrate-recognition domain of galectin-3). NMR relaxation experiments provide a unique probe of conformational entropy by characterizing bond-vector fluctuations at atomic resolution. By monitoring differences between the free and ligand-bound states in their backbone and side chain order parameters, we have estimated the contributions from conformational entropy to the free energy of binding. Overall, the conformational entropy of Gal3C increases upon ligand binding, thereby contributing favourably to the binding affinity. Comparisons with the results from isothermal titration calorimetry indicate that the conformational entropy is comparable in magnitude to the enthalpy of binding. Furthermore, there are significant differences in the dynamic response to binding of different ligands, despite the fact that the protein structure is virtually identical in the different protein–ligand complexes. Thus both affinity and specificity of ligand binding to Gal3C appear to depend in part on subtle differences in the conformational fluctuations that reflect the complex interplay between structure, dynamics and ligand interactions.

Matrix metalloproteinase-10 (MMP-10) interaction with tissue inhibitors of metalloproteinases TIMP-1 and TIMP-2: binding studies and crystal structure

Matrix metalloproteinase 10 (MMP-10, stromelysin-2) is a secreted metalloproteinase with functions in skeletal development, wound healing, and vascular remodeling; its overexpression is also implicated in lung tumorigenesis and tumor progression. To understand the regulation of MMP-10 by tissue inhibitors of metalloproteinases (TIMPs), we have assessed equilibrium inhibition constants (K(i)) of putative physiological inhibitors TIMP-1 and TIMP-2 for the active catalytic domain of human MMP-10 (MMP-10cd) using multiple kinetic approaches. We find that TIMP-1 inhibits the MMP-10cd with a K(i) of 1.1 × 10(-9) M; this interaction is 10-fold weaker than the inhibition of the similar MMP-3 (stromelysin-1) catalytic domain (MMP-3cd) by TIMP-1. TIMP-2 inhibits the MMP-10cd with a K(i) of 5.8 × 10(-9) M, which is again 10-fold weaker than the inhibition of MMP-3cd by this inhibitor (K(i) = 5.5 × 10(-10) M). We solved the x-ray crystal structure of TIMP-1 bound to the MMP-10cd at 1.9 Å resolution; the structure was solved by molecular replacement and refined with an R-factor of 0.215 (R(free) = 0.266). Comparing our structure of MMP-10cd·TIMP-1 with the previously solved structure of MMP-3cd·TIMP-1 (Protein Data Bank entry 1UEA), we see substantial differences at the binding interface that provide insight into the differential binding of stromelysin family members to TIMP-1. This structural information may ultimately assist in the design of more selective TIMP-based inhibitors tailored for specificity toward individual members of the stromelysin family, with potential therapeutic applications.

Long-range intra-protein communication can be transmitted by correlated side-chain fluctuations alone

Allosteric regulation is a key component of cellular communication, but the way in which information is passed from one site to another within a folded protein is not often clear. While backbone motions have long been considered essential for long-range information conveyance, side-chain motions have rarely been considered. In this work, we demonstrate their potential utility using Monte Carlo sampling of side-chain torsional angles on a fixed backbone to quantify correlations amongst side-chain inter-rotameric motions. Results indicate that long-range correlations of side-chain fluctuations can arise independently from several different types of interactions: steric repulsions, implicit solvent interactions, or hydrogen bonding and salt-bridge interactions. These robust correlations persist across the entire protein (up to 60 Å in the case of calmodulin) and can propagate long-range changes in side-chain variability in response to single residue perturbations.

Allosteric regulation occurs when the function of one part of a protein changes in response to a signal recognized by another part of the protein. Such intra-protein communication is essential for many biochemical processes, allowing the cell to adapt its behavior to a dynamic environment. Most studies of the information conveyance underlying allostery have to date focused on the role of backbone motions in mediating large structural changes. Here we focus instead on more subtle contributions, arising from fluctuations of side-chain torsions. Using a model for side-chain bond rotations in the tightly packed environment imposed by native backbone conformations, we observed significant sensitivity of side-chain organization to small, localized perturbations. This susceptibility arises from correlations among side-chain motions that can propagate information within a protein in complex, heterogeneous ways. Specifically, we found appreciable correlations even between side-chains distant from one another, so that the effect of a minor perturbation at one site on the protein could be observed in the altered fluctuations of side-chains throughout the protein. In conclusion, we have demonstrated that the statistical mechanics of correlated side-chain fluctuations within a model of the folded protein provides the basis for an unconventional but potentially important means of allostery.

Expanding the proteome: disordered and alternatively folded proteins

Proteins provide much of the scaffolding for life, as well as undertaking a variety of essential catalytic reactions. These characteristic functions have led us to presuppose that proteins are in general functional only when well structured and correctly folded. As we begin to explore the repertoire of possible protein sequences inherent in the human and other genomes, two stark facts that belie this supposition become clear: firstly, the number of apparent open reading frames in the human genome is significantly smaller than appears to be necessary to code for all of the diverse proteins in higher organisms, and secondly that a significant proportion of the protein sequences that would be coded by the genome would not be expected to form stable three-dimensional (3D) structures. Clearly the genome must include coding for a multitude of alternative forms of proteins, some of which may be partly or fully disordered or incompletely structured in their functional states. At the same time as this likelihood was recognized, experimental studies also began to uncover examples of important protein molecules and domains that were incompletely structured or completely disordered in solution, yet remained perfectly functional. In the ensuing years, we have seen an explosion of experimental and genome-annotation studies that have mapped the extent of the intrinsic disorder phenomenon and explored the possible biological rationales for its widespread occurrence. Answers to the question 'why would a particular domain need to be unstructured?' are as varied as the systems where such domains are found. This review provides a survey of recent new directions in this field, and includes an evaluation of the role not only of intrinsically disordered proteins but also of partially structured and highly dynamic members of the disorder-order continuum.

Phage display of tissue inhibitor of metalloproteinases-2 (TIMP-2): identification of selective inhibitors of collagenase-1 (metalloproteinase 1 (MMP-1))

Tissue inhibitor of metalloproteinases-2 (TIMP-2) is a broad spectrum inhibitor of the matrix metalloproteinases (MMPs), which function in extracellular matrix catabolism. Here, phage display was used to identify variants of human TIMP-2 that are selective inhibitors of human MMP-1, a collagenase whose unregulated action is linked to cancer, arthritis, and fibrosis. Using hard randomization of residues 2, 4, 5, and 6 (L1) and soft randomization of residues 34-40 (L2) and 67-70 (L3), a library was generated containing 2 × 10(10) variants of TIMP-2. Five clones were isolated after five rounds of selection with MMP-1, using MMP-3 as a competitor. The enriched phages selectively bound MMP-1 relative to MMP-3 and contained mutations only in L1. The most selective variant (TM8) was used to generate a second library in which residues Cys(1)-Gln(9) were soft-randomized. Four additional clones, selected from this library, showed a similar affinity for MMP-1 as wild-type TIMP-2 but reduced affinity for MMP-3. Variants of the N-terminal domain of TIMP-2 (N-TIMP-2) with the sequences of the most selective clones were expressed and characterized for inhibitory activity against eight MMPs. All were effective inhibitors of MMP-1 with nanomolar K(i) values, but TM8, containing Ser(2) to Asp and Ser(4) to Ala substitutions, was the most selective having a nanomolar K(i) value for MMP-1 but no detectable inhibitory activity toward MMP-3 and MMP-14 up to 10 μM. This study suggests that phage display and selection with other MMPs may be an effective method for discovering tissue inhibitor of metalloproteinase variants that discriminate between specified MMPs as targets.

Entropy increases from different sources support the high-affinity binding of the N-terminal inhibitory domains of tissue inhibitors of metalloproteinases to the catalytic domains of matrix metalloproteinases-1 and -3

The avid binding of tissue inhibitors of metalloproteinases (TIMPs) to matrix metalloproteinases (MMPs) is crucial for the regulation of pericellular and extracellular proteolysis. The interactions of the catalytic domain (cd) of MMP-1 with the inhibitory domains of TIMP-1 and TIMP-2 (N-TIMPs) and MMP-3cd with N-TIMP-2 have been characterized by isothermal titration calorimetry and compared with published data for the N-TIMP-1/MMP-3cd interaction. All interactions are largely driven by increases in entropy but there are significant differences in the profiles for the interactions of both N-TIMPs with MMP-1cd as compared with MMP-3cd; the enthalpy change ranges from small for MMP-1cd to highly unfavorable for MMP-3cd (-0.1 ± 0.7 versus 6.0 ± 0.5 kcal mol(-1)). The heat capacity change (ΔC(p)) of binding to MMP-1cd (temperature dependence of ΔH) is large and negative (-210 ± 20 cal K(-1) mol(-1)), indicating a large hydrophobic contribution, whereas the ΔC(p) values for the binding to MMP-3cd are much smaller (-53 ± 3 cal K(-1) mol(-1)), and some of the entropy increase may arise from increased conformational entropy. Apart from differences in ionization effects, it appears that the properties of the MMP may have a predominant influence in the thermodynamic profiles for these N-TIMP/MMP interactions.

Protein flexibility and conformational entropy in ligand design targeting the carbohydrate recognition domain of galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

Comparison of entropic contributions to binding in a "hydrophilic" versus "hydrophobic" ligand-protein interaction

In the present study we characterize the thermodynamics of binding of histamine to recombinant histamine-binding protein (rRaHBP2), a member of the lipocalin family isolated from the brown-ear tick Rhipicephalus appendiculatus. The binding pocket of this protein contains a number of charged residues, consistent with histamine binding, and is thus a typical example of a "hydrophilic" binder. In contrast, a second member of the lipocalin family, the recombinant major urinary protein (rMUP), binds small hydrophobic ligands, with a similar overall entropy of binding in comparison with rRaHBP2. Having extensively studied ligand binding thermodynamics for rMUP previously, the data we obtained in the present study for HBP enables a comparison of the driving forces for binding between these classically distinct binding processes in terms of entropic contributions from ligand, protein, and solvent. In the case of rRaHBP2, we find favorable entropic contributions to binding from desolvation of the ligand; however, the overall entropy of binding is unfavorable due to a dominant unfavorable contribution arising from the loss of ligand degrees of freedom, together with the sequestration of solvent water molecules into the binding pocket in the complex. This contrasts with binding in rMUP where desolvation of the protein binding pocket makes a minor contribution to the overall entropy of binding given that the pocket is substantially desolvated prior to binding.

The intrinsic protein flexibility of endogenous protease inhibitor TIMP-1 controls its binding interface and affects its function

Protein flexibility is thought to play key roles in numerous biological processes, including antibody affinity maturation, signal transduction, and enzyme catalysis, yet only limited information is available regarding the molecular details linking protein dynamics with function. A single point mutation at the distal site of the endogenous tissue inhibitor of metalloproteinase 1 (TIMP-1) enables this clinical target protein to tightly bind and inhibit membrane type 1 matrix metalloproteinase (MT1-MMP) by increasing only the association constant. The high-resolution X-ray structure of this complex determined at 2 A could not explain the mechanism of enhanced binding and pointed to a role for protein conformational dynamics. Molecular dynamics (MD) simulations reveal that the high-affinity TIMP-1 mutants exhibit significantly reduced binding interface flexibility and more stable hydrogen bond networks. This was accompanied by a redistribution of the ensemble of substrates to favorable binding conformations that fit the enzyme catalytic site. Apparently, the decrease in backbone flexibility led to a lower entropy cost upon formation of the complex. This work quantifies the effect of a single point mutation on the protein conformational dynamics and function of TIMP-1. Here we argue that controlling the intrinsic protein dynamics of MMP endogenous inhibitors may be utilized for rationalizing the design of selective novel protein inhibitors for this class of enzymes.

NMR and bioinformatics discovery of exosites that tune metalloelastase specificity for solubilized elastin and collagen triple helices

The catalytic domain of metalloelastase (matrix metalloproteinase-12 or MMP-12) is unique among MMPs in exerting high proteolytic activity upon fibrils that resist hydrolysis, especially elastin from lungs afflicted with chronic obstructive pulmonary disease or arteries with aneurysms. How does the MMP-12 catalytic domain achieve this specificity? NMR interface mapping suggests that α-elastin species cover the primed subsites, a strip across the β-sheet from β-strand IV to the II-III loop, and a broad bowl from helix A to helix C. The many contacts may account for the comparatively high affinity, as well as embedding of MMP-12 in damaged elastin fibrils in vivo. We developed a strategy called BINDSIght, for bioinformatics and NMR discovery of specificity of interactions, to evaluate MMP-12 specificity without a structure of a complex. BINDSIght integration of the interface mapping with other ambiguous information from sequences guided choice mutations in binding regions nearer the active site. Single substitutions at each of ten locations impair specific activity toward solubilized elastin. Five of them impair release of peptides from intact elastin fibrils. Eight lesions also impair specific activity toward triple helices from collagen IV or V. Eight sites map to the "primed" side in the III-IV, V-B, and S1' specificity loops. Two map to the "unprimed" side in the IV-V and B-C loops. The ten key residues circumscribe the catalytic cleft, form an exosite, and are distinctive features available for targeting by new diagnostics or therapeutics.

Interactions between PTB RRMs induce slow motions and increase RNA binding affinity

Polypyrimidine tract binding protein (PTB) participates in a variety of functions in eukaryotic cells, including alternative splicing, mRNA stabilization, and internal ribosomal entry site-mediated translation initiation. Its mechanism of RNA recognition is determined in part by the novel geometry of its two C-terminal RNA recognition motifs (RRM3 and RRM4), which interact with each other to form a stable complex (PTB1:34). This complex itself is unusual among RRMs, suggesting that it performs a specific function for the protein. In order to understand the advantage it provides to PTB, the fundamental properties of PTB1:34 are examined here as a comparative study of the complex and its two constituent RRMs. Both RRM3 and RRM4 adopt folded structures that NMR data show to be similar to their structure in PRB1:34. The RNA binding properties of the domains differ dramatically. The affinity of each separate RRM for polypyrimidine tracts is far weaker than that of PTB1:34, and simply mixing the two RRMs does not create an equivalent binding platform. (15)N NMR relaxation experiments show that PTB1:34 has slow, microsecond motions throughout both RRMs including the interdomain linker. This is in contrast to the individual domains, RRM3 and RRM4, where only a few backbone amides are flexible on this time scale. The slow backbone dynamics of PTB1:34, induced by packing of RRM3 and RRM4, could be essential for high-affinity binding to a flexible polypyrimidine tract RNA and also provide entropic compensation for its own formation.

The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity

Tissue inhibitors of metalloproteinases (TIMPs) are widely distributed in the animal kingdom and the human genome contains four paralogous genes encoding TIMPs 1 to 4. TIMPs were originally characterized as inhibitors of matrix metalloproteinases (MMPs), but their range of activities has now been found to be broader as it includes the inhibition of several of the disintegrin-metalloproteinases, ADAMs and ADAMTSs. TIMPs are therefore key regulators of the metalloproteinases that degrade the extracellular matrix and shed cell surface molecules. Structural studies of TIMP-MMP complexes have elucidated the inhibition mechanism of TIMPs and the multiple sites through which they interact with target enzymes, allowing the generation of TIMP variants that selectively inhibit different groups of metalloproteinases. Engineering such variants is complicated by the fact that TIMPs can undergo changes in molecular dynamics induced by their interactions with proteases. TIMPs also have biological activities that are independent of metalloproteinases; these include effects on cell growth and differentiation, cell migration, anti-angiogenesis, anti- and pro-apoptosis, and synaptic plasticity. Receptors responsible for some of these activities have been identified and their signaling pathways have been investigated. A series of studies using mice with specific TIMP gene deletions has illuminated the importance of these molecules in biology and pathology.

Resonance assignments for stromelysin complexed with a beta-sulfonyl hydroxamate inhibitor

The sequence specific (1)H, (13)C, and (15)N resonance assignments for stromelysin, a Matrix metalloproteinase, are reported in this article. Almost 70% of assignable backbone and side-chain atoms were assigned in this highly dynamic protein.

Ligand-induced changes in the structure and dynamics of Escherichia coli peptide deformylase

Peptide deformylase (PDF) is an enzyme that is responsible for removing the formyl group from nascently synthesized polypeptides in bacteria, attracting much attention as a potential target for novel antibacterial agents. Efforts to develop potent inhibitors of the enzyme have progressed on the basis of classical medicinal chemistry, combinatorial chemistry, and structural approaches, yet the validity of PDF as an antibacterial target hangs, in part, on the ability of inhibitors to selectively target this enzyme in favor of structurally related metallohydrolases. We have used (15)N NMR spectroscopy and isothermal titration calorimetry to investigate the high-affinity interaction of EcPDF with actinonin, a naturally occurring potent EcPDF inhibitor. Backbone amide chemical shifts, residual dipolar couplings, hydrogen-deuterium exchange, and (15)N relaxation reveal structural and dynamic effects of ligand binding in the immediate vicinity of the ligand-binding site as well as at remote sites. A comparison of the crystal structures of free and actinonin-bound EcPDF with the solution data suggests that most of the consequences of the ligand binding to the protein are lost or obscured during crystallization. The results of these studies improve our understanding of the thermodynamic global minimum and have important implications for structure-based drug design.

Calculation of proteins' total side-chain torsional entropy and its influence on protein-ligand interactions

Despite the high density within a typical protein fold, the ensemble of sterically permissible side-chain repackings is vast. Here, we examine the extent of this variability that survives energetic biases due to van der Waals interactions, hydrogen bonding, salt bridges, and solvation. Monte Carlo simulations of an atomistic model exhibit thermal fluctuations among a diverse set of side-chain arrangements, even with the peptide backbone fixed in its crystallographic conformation. We have quantified the torsional entropy of this native-state ensemble, relative to that of a noninteracting reference system, for 12 small proteins. The reduction in entropy per rotatable bond due to each kind of interaction is remarkably consistent across this set of molecules. To assess the biophysical importance of these fluctuations, we have estimated side-chain entropy contributions to the binding affinity of several peptide ligands with calmodulin. Calculations for our fixed-backbone model correlate very well with experimentally determined binding entropies over a range spanning more than 80 kJ/(mol x 308 K).

Configurational entropy in protein-peptide binding: computational study of Tsg101 ubiquitin E2 variant domain with an HIV-derived PTAP nonapeptide

Configurational entropy is thought to influence biomolecular processes, but there are still many open questions about this quantity, including its magnitude, its relationship to molecular structure, and the importance of correlation. The mutual information expansion (MIE) provides a novel and systematic approach to extracting configurational entropy changes due to correlated motions from molecular simulations. We present the first application of the MIE method to protein-ligand binding using multiple molecular dynamics simulations to study the association of the ubiquitin E2 variant domain of the protein Tsg101 and an HIV-derived nonapeptide. This investigation utilizes the second-order MIE approximation, which accounts for correlations between all pairs of degrees of freedom. The computed change in configurational entropy is large and has a major contribution from changes in pairwise correlation. The results also reveal intricate structure-entropy relationships. Thus, the present analysis suggests that in order for a model of binding to be accurate, it must include a careful accounting of configurational entropy changes.

Inactivation of N-TIMP-1 by N-terminal acetylation when expressed in bacteria

The high-affinity binding of tissue inhibitors of metalloproteinases (TIMPs) to matrix metalloproteinases (MMPs) is essential for regulation of the turnover of the extracellular matrix during development, wound healing, and progression of inflammatory diseases, such as cancer, atherosclerosis, and arthritis. Bacterially expressed N-terminal inhibitory domains of TIMPs (N-TIMPs) have been used extensively for biochemical and biophysical study of interactions with MMPs. Titration of N-TIMP-1 expressed in E. coli indicates, however, that only about 42% of the protein is active as an MMP inhibitor. The separation of inactive from fully active N-TIMP-1 has been achieved both by MMP affinity and by high-resolution cation exchange chromatography at an appropriate pH, based on a slight difference of charge. Purification by cation exchange chromatography with a Mono S column enriches the active portion of N-TIMP-1 to >95%, with K(i) of 1.5 nM for MMP-12. Mass spectra reveal that the inactive form differs from active N-TIMP-1 in being N-terminally acetylated, underscoring the importance of the free alpha-NH(2) of Cys1 for MMP inhibition. N(alpha)-acetylation of the CTCVPP sequence broadens the N-terminal sequence motifs reported to be susceptible to alpha-amino acetylation by E. coli N-acetyl transferases.

Residual ligand entropy in the binding of p-substituted benzenesulfonamide ligands to bovine carbonic anhydrase II

In studies on the thermodynamics of ligand-protein interactions, it is often assumed that the configurational and conformational entropy of the ligand is zero in the bound state (i.e., the ligand is rigidly fixed in the binding pocket). However, there is little direct experimental evidence for this assumption, and in the case of binding of p-substituted benzenesulfonamide inhibitors to bovine carbonic anhydrase II (BCA II), the observed thermodynamic binding signature derived from isothermal titration calorimetry experiments leads indirectly to the conclusion that a considerable degree of residual entropy remains in the bound ligand. Specifically, the entropy of binding increases with glycine chain length n, and strong evidence exists that this thermodynamic signature is not driven by solvent reorganization. By use of heteronuclear (15)N NMR relaxation measurements in a series (n = 1-6) of (15)N-glycine-enriched ligands, we find that the observed thermodynamic binding signature cannot be explained by residual ligand dynamics in the bound state, but rather results from the indirect influence of ligand chain length on protein dynamics.

Dynamic characterisation of the netrin-like domain of human type 1 procollagen C-proteinase enhancer and comparison to the N-terminal domain of tissue inhibitor of metalloproteinases (TIMP)

The backbone mobility of the C-terminal domain of procollagen C-proteinase enhancer (NTR PCOLCE1), part of a connective tissue glycoprotein, was determined using 15N NMR spectroscopy. NTR PCOLCE1 has been shown to be a netrin-like domain and adopts an OB-fold such as that found in the N-terminal domain of tissue inhibitors of metalloproteinases-1 (N-TIMP-1), N-TIMP-2, the laminin-binding domain of agrin and the C-terminal domain of complement protein C5. NMR relaxation dynamics of NTR PCOLCE1 highlight conformational flexibility in the N-terminus, strand A and the proximal CD loop. This region in N-TIMP is known to be essential for inhibitory activity against the matrix metalloproteinases and suggests that this region is of equal importance for NTR PCOLCE1, although the specific functional activity of the NTR PCOLCE1 domain is still unknown. Dynamics observed within the structural core of NTR PCOLCE1 that are not observed in N-TIMP molecules suggest that although the two domains have a similar architecture, the NTR PCOLCE1 domain will show different thermodynamic properties on binding and hence the target molecule could be somewhat different from that observed for the TIMPs. ModelFree order parameters show that NTR PCOLCE1 has more flexibility than both N-TIMP-1 and N-TIMP-2.

Compensatory and long-range changes in picosecond-nanosecond main-chain dynamics upon complex formation: 15N relaxation analysis of the free and bound states of the ubiquitin-like domain of human plexin-B1 and the small GTPase Rac1

The formation of a complex between Rac1 and the cytoplasmic domain of plexin-B1 is one of the first documented cases of a direct interaction between a small guanosine 5'-triphosphatase (GTPase) and a transmembrane receptor. Structural studies have begun to elucidate the role of this interaction for the signal transduction mechanism of plexins. Mapping of the Rac1 GTPase surface that contacts the Rho GTPase binding domain of plexin-B1 by solution NMR spectroscopy confirms the plexin domain as a GTPase effector protein. Regions neighboring the GTPase switch I and II regions are also involved in the interaction and there is considerable interest to examine the changes in protein dynamics that take place upon complex formation. Here we present main-chain nitrogen-15 relaxation measurements for the unbound proteins as well as for the Rho GTPase binding domain and Rac1 proteins in their complexed state. Derived order parameters, S2, show that considerable motions are maintained in the bound state of plexin. In fact, some of the changes in S2 on binding appear compensatory, exhibiting decreased as well as increased dynamics. Fluctuations in Rac1, already a largely rigid protein on the picosecond-nanosecond timescale, are overall diminished, but isomerization dynamics in the switch I and II regions of the GTPase are retained in the complex and appear to be propagated to the bound plexin domain. Remarkably, fluctuations in the GTPase are attenuated at sites, including helices alpha6 (the Rho-specific insert helix), alpha7 and alpha8, that are spatially distant from the interaction region with plexin. This effect of binding on long-range dynamics appears to be communicated by hinge sites and by subtle conformational changes in the protein. Similar to recent studies on other systems, we suggest that dynamical protein features are affected by allosteric mechanisms. Altered protein fluctuations are likely to prime the Rho GTPase-plexin complex for interactions with additional binding partners.

Restricted mobility of conserved residues in protein-protein interfaces in molecular simulations

Conserved residues in protein-protein interfaces correlate with residue hot-spots. To obtain insight into their roles, we have studied their mobility. We have performed 39 explicit solvent simulations of 15 complexes and their monomers, with the interfaces varying in size, shape, and function. The dynamic behavior of conserved residues in unbound monomers illustrates significantly lower flexibility as compared to their environment, suggesting that already before binding they are constrained in a boundlike configuration. To understand this behavior, we have analyzed the inter- and intrachain hydrogen-bond residence-time in the interfaces. We find that conserved residues are not involved significantly in hydrogen bonds across the interface as compared to nonconserved. However, the monomer simulations reveal that conserved residues contribute dominantly to hydrogen-bond formation before binding. Packing of conserved residues across the trajectories is significantly higher before and after the binding, rationalizing their lower mobility. Backbone torsional angle distributions show that conserved residues assume restricted regions of space and the most visited conformations in the bound and unbound trajectories are similar, suggesting that conserved residues are preorganized. Combined with previous studies, we conclude that conserved residues, hot spots, anchor, and interface-buried residues may be similar residues, fulfilling similar roles.

Extraction of configurational entropy from molecular simulations via an expansion approximation

A method is presented for extracting the configurational entropy of solute molecules from molecular dynamics simulations, in which the entropy is computed as an expansion of multidimensional mutual information terms, which account for correlated motions among the various internal degrees of freedom of the molecule. The mutual information expansion is demonstrated to be equivalent to estimating the full-dimensional configurational probability density function (PDF) using the generalized Kirkwood superposition approximation (GKSA). While the mutual information expansion is derived to the full dimensionality of the molecule, the current application uses a truncated form of the expansion in which all fourth- and higher-order mutual information terms are neglected. Truncation of the mutual information expansion at the nth order is shown to be equivalent to approximating the full-dimensional PDF using joint PDFs with dimensionality of n or smaller by successive application of the GKSA. The expansion method is used to compute the absolute (classical) configurational entropy in a basis of bond-angle-torsion internal coordinates for several small molecules as well as the change in entropy upon binding for a small host-guest system. Convergence properties of the computed entropy values as a function of simulation time are investigated and comparisons are made with entropy values from the second generation Mining Minima software. These comparisons demonstrate a deviation in -TS of no more than about 2 kcal/mol for all cases in which convergence has been obtained.

Constraining specificity in the N-domain of tissue inhibitor of metalloproteinases-1; gelatinase-selective inhibitors

The tissue inhibitors of metalloproteinases (TIMPs) are endogenous inhibitors of the matrix metalloproteinases (MMPs). Since unregulated MMP activities are linked to arthritis, cancer, and atherosclerosis, TIMP variants that are selective inhibitors of disease-related MMPs have potential therapeutic value. The structures of TIMP/MMP complexes reveal that most interactions with the MMP involve the N-terminal pentapeptide of TIMP and the C-D beta-strand connector which occupy the primed and unprimed regions of the active site. The loop between beta-strands A and B forms a secondary interaction site for some MMPs, ranging from multiple contacts in the TIMP-2/membrane type-1 (MT1)-MMP complex to none in the TIMP-1/MMP-1 complex. TIMP-1 and its inhibitory domain, N-TIMP-1, are weak inhibitors of MT1-MMP; inhibition is not improved by grafting the longer AB loop from TIMP-2 into N-TIMP-1, but this change impairs binding to MMP-3 and MMP-7. Mutational studies with N-TIMP-1 suggest that its weak inhibition of MT1-MMP, as compared to other N-TIMPs, arises from multiple (>3) sequence differences in the interaction site. Substitutions for Thr2 of N-TIMP-1 strongly influence MMP selectivity; Arg and Gly, that generally reduce MMP affinity, have less effect on binding to MMP-9. When the Arg mutation is added to the N-TIMP-1(AB2) mutant, it produces a gelatinase-specific inhibitor with Ki values of 2.8 and 0.4 nM for MMP-2 and -9, respectively. Interestingly, the Gly mutant has a Ki of 2.1 nM for MMP-9 and >40 muM for MMP-2, indicating that engineered TIMPs can discriminate between MMPs in the same subfamily.

Structures of MART-126/27-35 Peptide/HLA-A2 complexes reveal a remarkable disconnect between antigen structural homology and T cell recognition

Small structural changes in peptides presented by major histocompatibility complex (MHC) molecules often result in large changes in immunogenicity, supporting the notion that T cell receptors are exquisitely sensitive to antigen structure. Yet there are striking examples of TCR recognition of structurally dissimilar ligands. The resulting unpredictability of how T cells will respond to different or modified antigens impacts both our understanding of the physical bases for TCR specificity as well as efforts to engineer peptides for immunomodulation. In cancer immunotherapy, epitopes and variants derived from the MART-1/Melan-A protein are widely used as clinical vaccines. Two overlapping epitopes spanning amino acid residues 26 through 35 are of particular interest: numerous clinical studies have been performed using variants of the MART-1 26-35 decamer, although only the 27-35 nonamer has been found on the surface of targeted melanoma cells. Here, we show that the 26-35 and 27-35 peptides adopt strikingly different conformations when bound to HLA-A2. Nevertheless, clonally distinct MART-1(26/27-35)-reactive T cells show broad cross-reactivity towards these ligands. Simultaneously, however, many of the cross-reactive T cells remain unable to recognize anchor-modified variants with very subtle structural differences. These dichotomous observations challenge our thinking about how structural information on unligated peptide/MHC complexes should be best used when addressing questions of TCR specificity. Our findings also indicate that caution is warranted in the design of immunotherapeutics based on the MART-1 26/27-35 epitopes, as neither cross-reactivity nor selectivity is predictable based on the analysis of the structures alone.