The Role of Conformational Dynamics and Allostery in the Disease Development of Human Ferritin

Determining the three-dimensional structure of myoglobin, the first solved structure of a protein, fundamentally changed the way protein function was understood. Even more revolutionary was the information that came afterward: protein dynamics play a critical role in biological functions. Therefore, understanding conformational dynamics is crucial to obtaining a more complete picture of protein evolution. We recently analyzed the evolution of different protein families including green fluorescent proteins (GFPs), β-lactamase inhibitors, and nuclear receptors, and we observed that the alteration of conformational dynamics through allosteric regulation leads to functional changes. Moreover, proteome-wide conformational dynamics analysis of more than 100 human proteins showed that mutations occurring at rigid residue positions are more susceptible to disease than flexible residue positions. These studies suggest that disease-associated mutations may impair dynamic allosteric regulations, leading to loss of function. Thus, in this study, we analyzed the conformational dynamics of the wild-type light chain subunit of human ferritin protein along with the neutral and disease forms. We first performed replica exchange molecular dynamics simulations of wild-type and mutants to obtain equilibrated dynamics and then used perturbation response scanning (PRS), where we introduced a random Brownian kick to a position and computed the fluctuation response of the chain using linear response theory. Using this approach, we computed the dynamic flexibility index (DFI) for each position in the chain for the wild-type and the mutants. DFI quantifies the resilience of a position to a perturbation and provides a flexibility/rigidity measurement for a given position in the chain. The DFI analysis reveals that neutral variants and the wild-type exhibit similar flexibility profiles in which experimentally determined functionally critical sites act as hinges in controlling the overall motion. However, disease mutations alter the conformational dynamic profile, making hinges more loose (i.e., softening the hinges), thus impairing the allosterically regulated dynamics.     

Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP

α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.     

Two-dimensional NMR as a probe of structural similarity applied to mutants of cytochrome c

Using site-directed mutagenesis, it is possible to prepare many mutants of a protein in a short time, and to uncover differences in function. To understand the changes in function, it is essential to understand the effect(s) of the mutation in terms of structural and dynamic changes. It is particularly important to establish a rapid method for comparing the structure of the mutants with that of the wild-type protein. We propose that a combination of overlayed and difference two-dimensional NOE spectra between the wild-type and mutant protein provide a rapid method for determination of structural similarity. The observation of differences other than those due directly to the field effects of the exchanged side chain allow both local and distant conformational changes to be assessed. Here we compare NOESY spectra from a mutant of yeast iso-1-ferrocytochrome c in which the invariant residue Phe-82 has been changed to a Tyr. We conclude that NMR can show subtle changes in protein structure. Specifically, we show the change must involve the reorientation of the side chain of Leu-85 which is proximal to the mutation. The dynamics of the aromatic side chain at position 82 are shown not to give rise to measurable differences between the wild-type and mutant protein. Structural changes are not propagated to a measurable degree in other parts of the protein.     

Ca2+ binding and conformational change in two series of point mutations to the individual Ca(2+)-binding sites of calmodulin.

Two series of site-directed mutations to the individual Ca(2+)-binding sites of Drosophila melanogaster calmodulin have been generated and studied. In each mutant, a conserved glutamic acid residue at position 12 in all of the Ca(2+)-binding loops has been mutated in one site. In one series the residue is changed to glutamine; in the second series the change is to lysine. The Ca(2+)-binding properties of these mutants and the wild-type protein under pseudo-physiological conditions are presented. In addition, Ca(2+)-induced changes to the environment of the single tyrosine residue (Tyr-138) have been studied for some of the mutants. Ca2+ binding to the wild-type protein is best modeled as two pairs of sites with a higher affinity pair that shows strong cooperativity. For all but one of these eight mutant proteins, only three Ca(2+)-binding events can be detected. In three of the amino-terminal mutants, the three residual sites are (i) a pair of relatively high affinity sites and (ii) a weakened low affinity site. For all four carboxyl-terminal mutations, the residual sites are three relatively low affinity sites. In general, mutations to sites 2 and 4 prove more deleterious than mutations to sites 1 and 3. The Ca(2+)-induced conformational changes in the vicinity of Tyr-138 are relatively undisturbed by mutations of site 1. However, the changes to Tyr-138 in the carboxyl-terminal site mutants indicate that upon disruption of the cooperative binding at the high affinity sites, conformational change in the carboxyl terminus occurs in two phases. It appears that binding of Ca2+ to either carboxyl-terminal site can elicit the first phase of the response but the second phase is almost abolished when site 4 is the mutated site. The final conformations of site 3 and 4 mutants are thus significantly different.     

Structure and dynamics of the anti-AMCV scFv(F8): effects of selected mutations on the antigen combining site

The recombinant antibody fragment scFv(F8), which recognizes the coat protein of the plant virus AMCV, is characterized by peculiar high in vitro stability and functional folding even in reducing environments, making it fit for designing stable antibodies with desired properties. Mutagenesis and functional analysis evidenced two residues, at positions 47 and 58 of the V(H) chain, playing a crucial role in the antigen binding recognition. Here, we used a computational procedure to assess the effects of these mutations on the stability, structure and dynamics of the antigen-binding site. Structural models of the wild type scFv(F8) and of its H47 and H58 mutants were built by homology modelling and assessed by multiple 15.5ns of molecular dynamics simulations. Computational results indicate that the 47H substitution strongly affects the CDR-H(2) conformation, destabilizes the V(H)/V(L) interface and confers high conformational flexibility to the antigen-binding site, leading the mutant to functional loss. The mutation at position H58 strenghtens the binding site, bestowing a high antigen specificity on the mutant. The essential dynamics and the analysis of the protein-solvent interface further corroborate the correspondence between the extent of the structurally-determined flexibility of the binding site with the different functional behaviours proved by the wild-type and its mutants. These results may have useful implications for structure-based design of antibody combining site.     

Backbone and side chain dynamics of mutant calmodulin-peptide complexes

The mechanism of long-range coupling of allosteric sites in calcium-saturated calmodulin (CaM) has been explored by characterizing structural and dynamics effects of mutants of calmodulin in complex with a peptide corresponding to the smooth muscle myosin light chain kinase calmodulin-binding domain (smMLCKp). Four CaM mutants were examined: D95N and D58N, located in Ca2+-binding loops; and M124L and E84K, located in the target domain-binding site of CaM. Three of these mutants have altered allosteric coupling either between Ca2+-binding sites (D58N and D95N) or between the target- and Ca2+-binding sites (E84K). The structure and dynamics of the mutant calmodulins in complex with smMLCKp were characterized using solution NMR. Analysis of chemical shift perturbations was employed to detect largely structural perturbations. 15N and 2H relaxation was employed to detect perturbations of the dynamics of the backbone and methyl-bearing side chains of calmodulin. The least median squares method was found to be robust in the detection of perturbed sites. The main chain dynamics of calmodulin are found to be largely unresponsive to the mutations. Three mutants show significantly perturbed dynamics of methyl-bearing side chains. Despite the pseudosymmetric location of Ca2+-binding loop mutations D58N and D95N, the dynamic response of CaM is asymmetric, producing long-range perturbation in D58N and almost none in D95N. The mutations located at the target domain-binding site have quite different effects. For M124L, a local perturbation of the methyl dynamics is observed, while the E84K mutation produces a long-range propagation of dynamic perturbations along the target domain-binding site.     

Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication

In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.     

A Potential Yeast Actin Allosteric Conduit Dependent on Hydrophobic Core Residues Val-76 and Trp-79

Intramolecular allosteric interactions responsible for actin conformational regulation are largely unknown. Previous work demonstrated that replacing yeast actin Val-76 with muscle actin Ile caused decreased nucleotide exchange. Residue 76 abuts Trp-79 in a six-residue linear array beginning with Lys-118 on the surface and ending with His-73 in the nucleotide cleft. To test if altering the degree of packing of these two residues would affect actin dynamics, we constructed V76I, W79F, and W79Y single mutants as well as the Ile-76/Phe-79 and Ile-76/Tyr-79 double mutants. Tyr or Phe should decrease crowding and increase protein flexibility. Subsequent introduction of Ile should restore packing and dampen changes. All mutants showed decreased growth in liquid medium. W79Y alone was severely osmosensitive and exhibited vacuole abnormalities. Both properties were rescued by Ile-76. Phe-79 or Tyr decreased the thermostability of actin and increased its nucleotide exchange rate. These effects, generally greater for Tyr than for Phe, were reversed by introduction of Ile-76. HD exchange showed that the mutations caused propagated conformational changes to all four subdomains. Based on results from phosphate release and light-scattering assays, single mutations affected polymerization in the order of Ile, Phe, and Tyr from least to most. Introduction of Ile-76 partially rescued the polymerization defects caused by either Tyr-79 or Phe-79. Thus, alterations in crowding of the 76–79 residue pair can strongly affect actin conformation and behavior, and these results support the theory that the amino acid array in which they are located may play a central role in actin regulation.     

Insight into the conformational dynamics of specific regions of porcine pancreatic phospholipase A2 from a time-resolved fluorescence study of a genetically inserted single tryptophan residue

The effects of Ca2+ and substrate analogue binding on the conformational dynamics of porcine pancreas phospholipase A2 (PLA2) in different regions was explored by combining site-directed mutagenesis and time-resolved fluorescence measurements. The single tryptophan residue (Trp-3) of the wild-type protein (W3), in the alpha-helix A, was replaced by a phenylalanine residue (W3F), whereafter Trp was substituted either for leucine-31 (W31), located in the calcium binding loop, or for phenylalanine-94 (W94), located at the "back side" of the enzyme. Furthermore, mutants lacking the 62-66 sequence were constructed with the Trp at position 3 (delta W3) or 31 (delta W31). The total fluorescence intensity decays of Trp in each protein, in the protein-calcium and the protein-calcium-substrate analogue complexes, analyzed by the maximum entropy method (MEM) can be interpreted as distributions of separated lifetime classes. In the case of the W94 mutant, a major short-lived excited-state population (tau approximately 50 ps) is observed, probably deactivated by the interaction with two proximate disulfide bridges via a radiationless process. For the four other mutants, the respective barycenters of the four lifetime classes display comparable values, but the amplitude distributions are different for Trp-3 and Trp-31. The rotational mobility of the Trp residue varies along the peptide chain. Trp-3 experiences only a fast hindered motion. Trp-31 is sensitive to an additional local flexibility that is absent in the N-terminal part of the protein. The largest wobbling angle is observed at position 94. No effect of calcium binding occurs on the lifetime distribution of the Trp-3 and Trp-94 residues. Their mobilities are not affected. In contrast, calcium binding displays a strong influence on the excited-state population distribution of Trp-31. A major population decaying with the longest lifetime is selected in the W31 protein and contributes to approximately 50% of the decay. The local flexibility and the amplitude of motion of Trp-31 is wider in the protein-calcium complex than in the unliganded protein. Binding of the monomeric substrate analogue n-dodecylphosphocholine (C12PN) in the presence of calcium slightly affects the Trp-3 excited-state population distribution and its mobility. Trp-31 is more sensitive to this binding. In particular, a more restricted rotation of the Trp-31 residue and a decrease of the peptide local flexibility as protein-calcium complexes are observed in both the W31 and delta W31 mutants.(ABSTRACT TRUNCATED AT 400 WORDS)     

Hydrophobic core mutations in CI2 globally perturb fast side-chain dynamics similarly without regard to position

Protein dynamics is currently an area of intense research because of its importance as complementary information to the huge quantity of available data relating protein structure and function. Because it is usually the influence of dynamics on function that is studied, the physical determinants of the distribution of flexibility in proteins have not been explored as thoroughly. In the present NMR study, an expanded suite of five (2)H relaxation experiments was used to characterize the picosecond-to-nanosecond side-chain dynamics of chymotrypsin inhibitor 2 (CI2) and five hydrophobic core mutants, some of which are members of the folding nucleus. Because CI2 is a homologue of the serine protease inhibitor eglin c, which has already been extensively characterized in terms of its dynamics, it was possible to compare not only side-chain dynamics but also the responses of these dynamics to analogous mutations. Remarkably, each of the five core mutations in CI2 led to similar and reproducible increases in side-chain flexibility throughout the entire structure. Although the expanded suite of (2)H relaxation experiments does not affect model selection for the vast majority of residues, it did enable the detection of increasing levels of nanosecond-scale motions in CI2's reactive site binding loop as the L68 side chain was progressively shortened by mutation. Collectively, we observed that the CI2 mutants are more dynamically similar to each other than to the more rigid wild-type CI2, from which we propose that wild-type CI2 has been optimized to a specific level of rigidity which may aid in its function as a serine protease inhibitor. We also observed that the pattern of side-chain dynamics of CI2 is quantitatively similar to eglin c, but that this similarity is lost upon mutating both proteins at an equivalent position. Finally, (15)N relaxation was used to characterize the backbone dynamics of wild-type and mutant CI2. Interestingly, mutation at folding nucleus positions led to widespread increases in backbone flexibility, whereas non-folding-nucleus positions led to increases in flexibility in the C-terminal half of the protein only.     

Disease-associated mutations impacting BC-loop flexibility trigger long-range transthyretin tetramer destabilization and aggregation

Hereditary transthyretin amyloidosis (ATTR) is an autosomal dominant disease characterized by the extracellular deposition of the transport protein transthyretin (TTR) as amyloid fibrils. Despite the progress achieved in recent years, understanding why different TTR residue substitutions lead to different clinical manifestations remains elusive. Here, we studied the molecular basis of disease-causing missense mutations affecting residues R34 and K35. R34G and K35T variants cause vitreous amyloidosis, whereas R34T and K35N mutations result in amyloid polyneuropathy and restrictive cardiomyopathy. All variants are more sensitive to pH-induced dissociation and amyloid formation than the wild-type (WT)-TTR counterpart, specifically in the variants deposited in the eyes amyloid formation occurs close to physiological pHs. Chemical denaturation experiments indicate that all the mutants are less stable than WT-TTR, with the vitreous amyloidosis variants, R34G and K35T, being highly destabilized. Sequence-induced stabilization of the dimer–dimer interface with T119M rendered tetramers containing R34G or K35T mutations resistant to pH-induced aggregation. Because R34 and K35 are among the residues more distant to the TTR interface, their impact in this region is therefore theorized to occur at long range. The crystal structures of double mutants, R34G/T119M and K35T/T119M, together with molecular dynamics simulations indicate that their strong destabilizing effect is initiated locally at the BC loop, increasing its flexibility in a mutation-dependent manner. Overall, the present findings help us to understand the sequence-dynamic-structural mechanistic details of TTR amyloid aggregation triggered by R34 and K35 variants and to link the degree of mutation-induced conformational flexibility to protein aggregation propensity.     

Exploring the role of structure and dynamics in the function of chymotrypsin inhibitor 2

Increasing awareness of the possible role of internal dynamics in protein function has led to the development of new methods for experimentally characterizing protein dynamics across multiple time scales, especially using NMR spectroscopy. A few analyses of the conformational dynamics of proteins ranging from nonallosteric single domains to multidomain allosteric enzymes are now available; however, demonstrating a connection between dynamics and function remains difficult on account of the comparative lack of studies examining both changes in dynamics and changes in function in response to the same perturbations. In previous work, we characterized changes in structure and dynamics on the ps–ns time scale resulting from hydrophobic core mutations in chymotrypsin inhibitor 2 and found that there are moderate, persistent global changes in dynamics in the absence of gross structural changes (Whitley et al., Biochemistry 2008;47:8566–8576). Here, we assay those and additional mutants for inhibitory ability toward the serine proteases elastase and chymotrypsin to determine the effects of mutation on function. Results indicate that core mutation has only a subtle effect on CI2 function. Using chemical shifts, we also studied the effect of complex formation on CI2 structure and found that perturbations are greatest at the complex interface but also propagate toward CI2's hydrophobic core. The structure–dynamics–function data set completed here suggests that dynamics plays a limited role in the function of this small model system, although we do observe a correlation between nanosecond-scale reactive loop motions and inhibitory ability for mutations at one key position in the hydrophobic core.     

Analysis of side chain mobility among protein G B1 domain mutants with widely varying stabilities

"Host-guest" studies of the B1 domain from Streptococcal protein G have been used previously to establish a thermodynamic scale for the beta-sheet-forming propensities of the 20 common amino acids. To investigate the contribution of side chain conformational entropy to the relative stabilities of B1 domain mutants, we have determined the dynamics of side chain methyl groups in 10 of the 20 mutants used in a previous study. Deuterium relaxation rates were measured using two-dimensional NMR techniques for 13CH2D groups. Analysis of the relaxation data using the Lipari-Szabo model-free formalism showed that mutations introduced at the guest position caused small but statistically significant changes in the methyl group dynamics. In addition, there was a low level of covariation of the Lipari-Szabo order parameters among the 10 mutants. The variations in conformational free energy estimated from the order parameters were comparable in magnitude to the variations in global stability of the 10 mutants but did not correlate with the global stability of the domain or with the structural properties of the guest amino acids. The data support the view that conformational entropy in the folded state is one of many factors that can influence the folding thermodynamics of proteins.     

Ensemble-based computational approach discriminates functional activity of p53 cancer and rescue mutants

The tumor suppressor protein p53 can lose its function upon single-point missense mutations in the core DNA-binding domain (“cancer mutants”). Activity can be restored by second-site suppressor mutations (“rescue mutants”). This paper relates the functional activity of p53 cancer and rescue mutants to their overall molecular dynamics (MD), without focusing on local structural details. A novel global measure of protein flexibility for the p53 core DNA-binding domain, the number of clusters at a certain RMSD cutoff, was computed by clustering over 0.7 µs of explicitly solvated all-atom MD simulations. For wild-type p53 and a sample of p53 cancer or rescue mutants, the number of clusters was a good predictor of in vivo p53 functional activity in cell-based assays. This number-of-clusters (NOC) metric was strongly correlated (r² = 0.77) with reported values of experimentally measured ΔΔG protein thermodynamic stability. Interpreting the number of clusters as a measure of protein flexibility: (i) p53 cancer mutants were more flexible than wild-type protein, (ii) second-site rescue mutations decreased the flexibility of cancer mutants, and (iii) negative controls of non-rescue second-site mutants did not. This new method reflects the overall stability of the p53 core domain and can discriminate which second-site mutations restore activity to p53 cancer mutants.     

A new concept to reveal protein dynamics based on energy dissipation

Protein dynamics is essential for its function, especially for intramolecular signal transduction. In this work we propose a new concept, energy dissipation model, to systematically reveal protein dynamics upon effector binding and energy perturbation. The concept is applied to better understand the intramolecular signal transduction during allostery of enzymes. The E. coli allosteric enzyme, aspartokinase III, is used as a model system and special molecular dynamics simulations are designed and carried out. Computational results indicate that the number of residues affected by external energy perturbation (i.e. caused by a ligand binding) during the energy dissipation process shows a sigmoid pattern. Using the two-state Boltzmann equation, we define two parameters, the half response time and the dissipation rate constant, which can be used to well characterize the energy dissipation process. For the allostery of aspartokinase III, the residue response time indicates that besides the ACT2 signal transduction pathway, there is another pathway between the regulatory site and the catalytic site, which is suggested to be the β15-αK loop of ACT1. We further introduce the term "protein dynamical modules" based on the residue response time. Different from the protein structural modules which merely provide information about the structural stability of proteins, protein dynamical modules could reveal protein characteristics from the perspective of dynamics. Finally, the energy dissipation model is applied to investigate E. coli aspartokinase III mutations to better understand the desensitization of product feedback inhibition via allostery. In conclusion, the new concept proposed in this paper gives a novel holistic view of protein dynamics, a key question in biology with high impacts for both biotechnology and biomedicine.     

Identification of regions of potential flexibility in protein structures: folding units and correlations with intron positions

An algorithm has been developed to estimate flexibility for potential hinge motion at specified residues, that is, the mutual movement of two domains by rotation around a set of main-chain dihedral angles with torsion angles of neighboring side chains as variables. Such conformational changes must occur without severe atomic collisions. Flexible hinges have been found that satisfy such criteria. Sequence flexibility charts were obtained by plotting the flexibility of each residue against the residue number. Such charts were calculated for 10 proteins (ovomucoid third domain, cytochrome c, lysozyme, hemoglobin β-chain, α-chymotrypsin, elastase, carboxypeptidase A, dihydrofolate reductase, triosephosphate isomerase, and alcohol dehydrogenase) taken from the Protein Data Bank. The first step of unfolding is likely to occur at the hinge point with the largest flexibility. Following this idea, the polypeptide chain can be dissected into several folding units according to the sequence flexibility chart. When two domains are separated by conformational changes at such a hinge, the sequence flexibility chart for each domain changes, and it is recalculated and used to indicate subsequent unfolding steps. In this process of iterative estimation of flexibile hinges, some well-isolated hinges, or the border line between flexible and inflexible regions, were found to be directly at or close to the positions of splice junctions in the eukaryotic genes. Of a total of 45 splice junctions in the 10 proteins examined in this paper, 38 junctions can be identified as flexible hinges between folding units. We suggest that the iterative estimation of flexible hinges may define an array of possible folding/unfolding paths, and that the exon–intron arrangement in the gene may be closely correlated with the folding process of the protein.     

Methionine and alanine substitutions show that the formation of wild-type-like structure in the carboxy-terminal domain of T4 lysozyme is a rate-limiting step in folding.

In an attempt to identify a systematic relation between the structure of a protein and its folding kinetics, the rate of folding was determined for 20 mutants of T4 lysozyme in which a bulky, buried, nonpolar wild-type residue (Leu, Ile, Phe, Val, or Met) was substituted with alanine. Methionine, which approximated the size of the original side chain but which is of different shape and flexibility, was also substituted at most of the same sites. Mutations that substantially destabilize the protein and are located in the carboxy-terminal domain generally slow the rate of folding. Destabilizing mutations in the amino-terminal domain, however, have little effect on the rate of folding. Mutations that have little effect on stability tend to have little effect on the rate, no matter where they are located. These results suggest that, at the rate-limiting step, elements of structure in the C-terminal domain are formed and have a structure similar to that of the fully folded protein. Consistent with this, two variants that somewhat increase the rate of folding (Phe104 --> Met and Val149 --> Met) are located within the carboxy-terminal domain and maintain or improve packing with very little perturbation of the wild-type structure.     

The conformational flexibility of the carboxy terminal residues 105-114 is a key modulator of the catalytic activity and stability of macrophage migration inhibitory factor

Macrophage migration inhibitory factor (MIF) is a multifunctional protein and a major mediator of innate immunity. Although X-ray crystallography revealed that MIF exists as a homotrimer, its oligomerization state in vivo and the factors governing its oligomerization and stability remain poorly understood. The C-terminal region of MIF is highly conserved and participates in several intramolecular interactions that suggest a role in modulating the stability and biochemical activity of MIF. To determine the importance of these interactions, point mutations (A48P, L46A), insertions (P107) at the monomer-monomer interfaces, and C-terminal deletion (Delta 110-114NSTFA and Delta 105-114NVGWNNSTFA) variants were designed and their structural properties, thermodynamic stability, oligomerization state, catalytic activity and receptor binding were characterized using a battery of biophysical methods. The C-terminal deletion mutants DeltaC5 huMIF 1-109 and DeltaC10 huMIF 1-104 were enzymatically inactive and thermodynamically less stable than wild type MIF. Analytical ultracentrifugation studies demonstrate that both C-terminal mutants sediment as trimers and exhibit similar binding to CD74 as the wild type protein. Disrupting the conformation of the C-terminal region 105-114 and increasing its conformational flexibility through the insertion of a proline residue at position 107 was sufficient to reproduce the structural, biochemical and thermodynamic properties of the deletion mutants. P107 MIF forms an enzymatically inactive trimer and exhibits reduced thermodynamic stability relative to the wild type protein. To provide a rationale for the changes induced by these mutations at the molecular level, we also performed molecular dynamics simulations on these mutants in comparison to the wild type MIF. Together, our studies demonstrate that intersubunit interactions involving the C-terminal region 105-114, including a salt-bridge interaction between Arg73 of one monomer and the carboxy terminus of a neighboring monomer, play critical roles in modulating tertiary structure stabilization, enzymatic activity, and thermodynamic stability of MIF, but not its oligomerization state and receptor binding properties. Our results suggest that targeting the C-terminal region could provide new strategies for allosteric modulation of MIF enzymatic activity and the development of novel inhibitors of MIF tautomerase activity.     

The mobility of an HIV-1 integrase active site loop is correlated with catalytic activity.

Replication of HIV-1 requires the covalent integration of the viral cDNA into the host chromosomal DNA directed by the virus-encoded integrase protein. Here we explore the importance of a protein surface loop near the integrase active site using protein engineering and X-ray crystallography. We have redetermined the structure of the integrase catalytic domain (residues 50-212) using an independent phase set at 1.7 A resolution. The structure extends helix alpha4 on its N-terminal side (residues 149-154), thus defining the position of the three conserved active site residues. Evident in this and in previous structures is a conformationally flexible loop composed of residues 141-148. To probe the role of flexibility in this loop, we replaced Gly 140 and Gly 149, residues that appear to act as conformational hinges, with Ala residues. X-ray structures of the catalytic domain mutants G149A and G140A/G149A show further rigidity of alpha4 and the adjoining loop. Activity assays in vitro revealed that these mutants are impaired in catalysis. The DNA binding affinity, however, is minimally affected by these mutants as assayed by UV cross-linking. We propose that the conformational flexibility of this active site loop is important for a postbinding catalytic step.