Tissue inhibitor of metalloproteinases-1 undergoes microsecond to millisecond motions at sites of matrix metalloproteinase-induced fit

Conformational dynamics in interleukin 17A and 17F functional complexes is a key determinant of receptor A affinity and specificity

The proinflammatory cytokines IL-17A and IL-17F have been identified as key drivers of a range of human inflammatory diseases, such as psoriasis, which has led to several therapeutic antibodies targeted at IL-17A. The two cytokines have been shown to tightly associate as functional homo and hetero dimers, which induce signalling via the formation of a cell surface signalling complex with a single copy of both IL-17RA and IL-17RC. Striking differences in affinity have been observed for IL-17RA binding to IL-17AA, IL-17AF and IL-17FF, however, the functional significance and molecular basis for this has remained unclear. We have obtained comprehensive backbone NMR assignments for full length IL-17AA (79%), IL-17AF (93%) and IL-17FF (89%), which show that the dimers adopt almost identical backbone topologies in solution to those observed in reported crystal structures. Analysis of the line widths and intensities of assigned backbone amide NMR signals has revealed striking differences in the conformational plasticity and dynamics of IL-17AA compared to both IL-17AF and IL-17FF. Our NMR data indicate that a number of regions of IL-17AA are interconverting between at least two distinct conformations on a relatively slow timescale. Such conformational heterogeneity has previously been shown to play an important role in the formation of many high affinity protein-protein complexes. The locations of the affected IL-17AA residues essentially coincides with the regions of both IL-17A and IL-17F previously shown to undergo significant structural changes on binding to IL-17RA. Substantially less conformational exchange was revealed by the NMR data for IL-17FF and IL-17AF. We propose that the markedly different conformational dynamic properties of the distinct functional IL-17 dimers plays a key role in determining their affinities for IL-17RA, with the more dynamic and plastic nature of IL-17AA contributing to the significantly tighter affinity observed for binding to IL-17RA. In contrast, the dynamic properties are expected to have little influence on the affinity of IL-17 dimers for IL-17RC, which has recently been shown to induce only small structural changes in IL-17FF upon binding.

Universality and Specificity in Protein Fluctuation Dynamics

We investigate the universal scaling of protein fluctuation dynamics with a site-specific diffusive model of protein motion, which predicts an initial subdiffusive regime in the configurational relaxation. The long-time dynamics of proteins is controlled by an activated regime. We argue that the hierarchical free energy barriers set the time scales of biological processes and establish an upper limit to the size of single protein domains. We find it compelling that the scaling behavior for the protein dynamics is in close agreement with the Kardar-Parisi-Zhang scaling exponents.

Predicting protein dynamics from structural ensembles

The biological properties of proteins are uniquely determined by their structure and dynamics. A protein in solution populates a structural ensemble of metastable configurations around the global fold. From overall rotation to local fluctuations, the dynamics of proteins can cover several orders of magnitude in time scales. We propose a simulation-free coarse-grained approach which utilizes knowledge of the important metastable folded states of the protein to predict the protein dynamics. This approach is based upon the Langevin Equation for Protein Dynamics (LE4PD), a Langevin formalism in the coordinates of the protein backbone. The linear modes of this Langevin formalism organize the fluctuations of the protein, so that more extended dynamical cooperativity relates to increasing energy barriers to mode diffusion. The accuracy of the LE4PD is verified by analyzing the predicted dynamics across a set of seven different proteins for which both relaxation data and NMR solution structures are available. Using experimental NMR conformers as the input structural ensembles, LE4PD predicts quantitatively accurate results, with correlation coefficient ρ = 0.93 to NMR backbone relaxation measurements for the seven proteins. The NMR solution structure derived ensemble and predicted dynamical relaxation is compared with molecular dynamics simulation-derived structural ensembles and LE4PD predictions and is consistent in the time scale of the simulations. The use of the experimental NMR conformers frees the approach from computationally demanding simulations.

Using NMR spectroscopy to elucidate the role of molecular motions in enzyme function

Conformational motions play an essential role in enzyme function, often facilitating the formation of enzyme-substrate complexes and/or product release. Although considerable debate remains regarding the role of molecular motions in the conversion of enzymatic substrates to products, numerous examples have found motions to be crucial for optimization of enzyme scaffolds, effective substrate binding, and product dissociation. Conformational fluctuations are often rate-limiting to enzyme catalysis, primarily through product release, with the chemical reaction occurring much more quickly. As a result, the direct involvement of motions at various stages along the enzyme reaction coordinate remains largely unknown and untested. In the following review, we describe the use of solution NMR techniques designed to probe various timescales of molecular motions and detail examples in which motions play a role in propagating catalytic effects from the active site and directly participate in essential aspects of enzyme function.

Conformational heterogeneity in antibody-protein antigen recognition: implications for high affinity protein complex formation

Specific, high affinity protein-protein interactions lie at the heart of many essential biological processes, including the recognition of an apparently limitless range of foreign proteins by natural antibodies, which has been exploited to develop therapeutic antibodies. To mediate biological processes, high affinity protein complexes need to form on appropriate, relatively rapid timescales, which presents a challenge for the productive engagement of complexes with large and complex contact surfaces (∼600-1800 Å(2)). We have obtained comprehensive backbone NMR assignments for two distinct, high affinity antibody fragments (single chain variable and antigen-binding (Fab) fragments), which recognize the structurally diverse cytokines interleukin-1β (IL-1β, β-sheet) and interleukin-6 (IL-6, α-helical). NMR studies have revealed that the hearts of the antigen binding sites in both free anti-IL-1β Fab and anti-IL-6 single chain variable exist in multiple conformations, which interconvert on a timescale comparable with the rates of antibody-antigen complex formation. In addition, we have identified a conserved antigen binding-induced change in the orientation of the two variable domains. The observed conformational heterogeneity and slow dynamics at protein antigen binding sites appears to be a conserved feature of many high affinity protein-protein interfaces structurally characterized by NMR, suggesting an essential role in protein complex formation. We propose that this behavior may reflect a soft capture, protein-protein docking mechanism, facilitating formation of high affinity protein complexes on a timescale consistent with biological processes.

Mycobacterium tuberculosis RNA polymerase-binding protein A (RbpA) and its interactions with sigma factors

RNA polymerase-binding protein A (RbpA), encoded by Rv2050, is specific to the actinomycetes, where it is highly conserved. In the pathogen Mycobacterium tuberculosis, RbpA is essential for growth and survival. RbpA binds to the β subunit of the RNA polymerase where it activates transcription by unknown mechanisms, and it may also influence the response of M. tuberculosis to the current frontline anti-tuberculosis drug rifampicin. Here we report the solution structure of RbpA and identify the principle sigma factor σ(A) and the stress-induced σ(B) as interaction partners. The protein has a central ordered domain with a conserved hydrophobic surface that may be a potential protein interaction site. The N and C termini are highly dynamic and are involved in the interaction with the sigma factors. RbpA forms a tight complex with the N-terminal domain of σ(B) via its N- and C-terminal regions. The interaction with sigma factors may explain how RbpA stabilizes sigma subunit binding to the core RNA polymerase and thereby promotes initiation complex formation. RbpA could therefore influence the competition between principal and alternative sigma factors and hence the transcription profile of the cell.

Relative stabilities of conserved and non-conserved structures in the OB-fold superfamily

The OB-fold is a diverse structure superfamily based on a beta-barrel motif that is often supplemented with additional non-conserved secondary structures. Previous deletion mutagenesis and NMR hydrogen exchange studies of three OB-fold proteins showed that the structural stabilities of sites within the conserved beta-barrels were larger than sites in non-conserved segments. In this work we examined a database of 80 representative domain structures currently classified as OB-folds, to establish the basis of this effect. Residue-specific values were obtained for the number of Calpha-Calpha distance contacts, sequence hydrophobicities, crystallographic B-factors, and theoretical B-factors calculated from a Gaussian Network Model. All four parameters point to a larger average flexibility for the non-conserved structures compared to the conserved beta-barrels. The theoretical B-factors and contact densities show the highest sensitivity. Our results suggest a model of protein structure evolution in which novel structural features develop at the periphery of conserved motifs. Core residues are more resistant to structural changes during evolution since their substitution would disrupt a larger number of interactions. Similar factors are likely to account for the differences in stability to unfolding between conserved and non-conserved structures.

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.

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.

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.

Threonine 98, the pivotal residue of tissue inhibitor of metalloproteinases (TIMP)-1 in metalloproteinase recognition

Tissue inhibitors of metalloproteinases (TIMPs) are the endogenous modulators of the zinc-dependent mammalian matrix metalloproteinases (MMPs) and their close associates, proteinases of the ADAM (a disintegrin and metalloproteinase) and ADAM with thrombospondin repeats families. There are four variants of TIMPs, and each has its defined set of metalloproteinase (MP) targets. TIMP-1, in particular, is inactive against several of the membrane-type MMPs (MT-MMPs), MMP-19, and the ADAM proteinase TACE (tumor necrosis factor-alpha-converting enzyme, ADAM-17). The molecular basis for such inactivity is unknown. Previously, we showed that TIMP-1 could be transformed into an active inhibitor against MT1-MMP by the replacement of threonine 98 residue with leucine (T98L). Here, we reveal that the T98L mutation has in fact transformed TIMP-1 into a versatile inhibitor against an array of MPs otherwise insensitive to wild-type TIMP-1; examples include TACE, MMP-19, and MT5-MMP. Using T98L as the scaffold, we created a TIMP-1 variant that is fully active against TACE. The binding affinity of the mutant (V4S/TIMP-3-AB-loop/V69L/T98L) (K (app)(i) 0.14 nm) surpassed that of TIMP-3 (K (app)(i) 0.22 nm), the only natural TIMP inhibitor of the enzyme. The requirement for leucine is absolute for the transformation in inhibitory pattern. On the other hand, the mutation has minimal impact on the MPs already well inhibited by wild-type TIMP-1, such as gelatinase-A and stromelysin-1. Not only have we unlocked the molecular basis for the inactivity of TIMP-1 against several of the MPs, but also our findings fundamentally modify the current beliefs on the molecular mechanism of TIMP-MP recognition and selectivity.

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.

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

The high-affinity inhibition of stromelysin 1 (MMP-3) by tissue inhibitor of metalloproteinases 1 (TIMP-1) helps control tissue remodeling and tumor development. The interaction of N-TIMP-1 with the catalytic domain of MMP-3 has been investigated by titration calorimetry and 15N NMR. Their unfavorable enthalpy of binding of +6.5 kcal mol(-1) is unusual among protein-protein associations, deviates from structure-based prediction, and is compensated by a net entropy increase providing at least 18 kcal mol(-1) of favorable free energy of binding at a 1M reference state. The small heat capacity of binding agrees well with the heat capacity predicted from 65% of the surface buried on binding being polar, and suggests that the hydrophobic effect can account for only part of the entropy of binding. Using NMR, binding-induced changes in the backbone of N-TIMP-1 were checked as one possible source of conformational entropy changes. MMP binding slightly increases rigidity in some contact sites in TIMP-1 but increases mobility remotely in the otherwise rigid beta-barrel core of N-TIMP-1, increasing 15N relaxation evidence of pico- to nanosecond and micro- to millisecond fluctuations of beta-strands A-F. Residual dipolar couplings suggest dynamic deviations from X-ray coordinates of the complex. These suggest that the beta-barrel has small backbone conformational fluctuations, while segments of strands betaB, betaE and betaF might experience fluctuations only in their backbone environment. This is a distinctive example of affinity between two well-structured proteins being enhanced by increased conformational entropy in the reservoir of a folding core.

Backbone dynamics and thermodynamics of Borrelia outer surface protein A

Nuclear spin relaxation experiments performed at 298K, 308K and 318K are used to characterize the intramolecular dynamics and thermodynamics of outer surface protein A (OspA), a key protein in the life-cycle of Borrelia burgdorferi, the causative agent of Lyme disease. It has recently been demonstrated that OspA specifically binds to the gut of the intermediate tick host (Ixodes scapularis), and that this interaction is mediated, at least in part, by residues in the C-terminal domain of OspA that are largely inaccessible to solvent in all X-ray structures of this protein. Our analysis of 15N relaxation parameters in OspA shows that the putative-binding region contains and is surrounded by flexible residues, which could facilitate accessibility to solvent and ligands. In addition, residues with similar activation energies are clustered in a manner that suggests locally collective motions. We have used molecular modeling to show that these collective motions are consistent with a hinge-bending mechanism that exposes residues implicated in binding. Characteristic temperatures describing the energy landscape of the OspA backbone are derived from the temperature dependence of the N-H bond vector order parameters, and a comparison is made between the N and C-terminal globular domains and the unusual single-layer beta-sheet connecting them. The average characteristic temperatures in the three regions indicate that, with an increase in temperature, a larger increase in accessible conformational states occurs for N-H bond vectors in the single-layer central beta-sheet than for bond vectors in the globular N and C-terminal domains. These conformational states are accessible without disruption of hydrogen bonds, providing a conformational entropic gain, upon increase in temperature, without a significant enthalpic penalty. This increase in heat capacity may help to explain the unexpected thermal stability of the unusual single-layer beta-sheet.

Hsc70-interacting HPD loop of the J domain of polyomavirus T antigens fluctuates in ps to ns and micros to ms

The backbone dynamics of the J domain from polyomavirus T antigens have been investigated using 15N NMR relaxation and molecular dynamics simulation. Model-free relaxation analysis revealed picosecond to nanosecond motions in the N terminus, the I-II loop, the C-terminal end of helix II through the HPD loop to the beginning of helix III, and the C-terminal end of helix III to the C terminus. The backbone dynamics of the HPD loop and termini are dominated by motions with moderately large amplitudes and correlation times of the order of a nanosecond or longer. Conformational exchange on the microsecond to millisecond timescale was identified in the HPD loop, the N and C termini, and the I-II loop. A 9.7ns MD trajectory manifested concerted swings of the HPD loop. Transitions between major and minor conformations of the HPD loop featured distinct patterns of change in backbone dihedral angles and hydrogen bonds. Fraying of the C-terminal end of helix II and the N-terminal end of helix III correlated with displacements of the HPD loop. Correlation of crankshaft motions of Gly46 and Gly47 with the collective motions of the HPD loop suggested an important role of the two glycine residues in the mobility of the loop. Fluctuations of the HPD loop correlated with relative reorientation of side-chains of Lys35 and Asp44 that interact with Hsc70.

The role of backbone motions in ligand binding to the c-Src SH3 domain

The Src homology 3 (SH3) domain of pp60(c-src) (Src) plays dual roles in signal transduction, through stabilizing the repressed form of the Src kinase and through mediating the formation of activated signaling complexes. Transition of the Src SH3 domain between a variety of binding partners during progression through the cell cycle requires adjustment of a delicate free energy balance. Although numerous structural and functional studies of SH3 have provided an in-depth understanding of structural determinants for binding, the origins of binding energy in SH3-ligand interactions are not fully understood. Considering only the protein-ligand interface, the observed favorable change in standard enthalpy (DeltaH=-9.1 kcal/mol) and unfavorable change in standard entropy (TDeltaS=-2.7 kcal/mol) upon binding the proline-rich ligand RLP2 (RALPPLPRY) are inconsistent with the predominantly hydrophobic interaction surface. To investigate possible origins of ligand binding energy, backbone dynamics of free and RLP2-bound SH3 were performed via (15)N NMR relaxation and hydrogen-deuterium (H/(2)H) exchange measurements. On the ps-ns time scale, assuming uncorrelated motions, ligand binding results in a significant reduction in backbone entropy (-1.5(+/-0.6) kcal/mol). Binding also suppresses motions on the micros-ms time scale, which may additionally contribute to an unfavorable change in entropy. A large increase in protection from H/(2)H exchange is observed upon ligand binding, providing evidence for entropy loss due to motions on longer time scales, and supporting the notion that stabilization of pre-existing conformations within a native state ensemble is a fundamental paradigm for ligand binding. Observed changes in motion on all three time scales occur at locations both near and remote from the protein-ligand interface. The propagation of ligand binding interactions across the SH3 domain has potential consequences in target selection through altering both free energy and geometry in intact Src, and suggests that looking beyond the protein-ligand interface is essential in understanding ligand binding energetics.