Effects of point mutation in a flexible loop on the stability and enzymatic function of Escherichia coli dihydrofolate reductase

Identification of molecular basis that underlie enzymatic specificity of AzoRo from Rhodococcus opacus 1CP: A potential NADH:quinone oxidoreductase

Azo dyes are important to various industries such as textile industries. However, these dyes are known to comprise toxic, mutagenic, and carcinogenic representatives. Several approaches have already been employed to mitigate the problem such as the use of enzymes. Azoreductases have been well-studied in its capability to reduce azo dyes. AzoRo from Rhodococcus opacus 1CP has been found to be accepting only methyl red as a substrate, surmising that the enzyme may have a narrow active site. To determine the active site configuration of AzoRo at atomic level and identify the key residues involved in substrate binding and enzyme specificity, we have determined the crystal structure of holo-AzoRo and employed a rational design approach to generate AzoRo variants. The results reported here show that AzoRo has a different configuration of the active site when compared with other bacterial NAD(P)H azoreductases, having other key residues playing a role in the substrate binding and restricting the enzyme activity towards different azo dyes. Moreover, it was observed that AzoRo has only about 50% coupling yield to methyl red and p-benzoquinone - giving rise to the possibility that NADH oxidation still occurs even during catalysis. Results also showed that AzoRo is more active and more efficient towards quinones (about four times higher than methyl red).

Site-Specific Tryptophan Labels Reveal Local Microsecond-Millisecond Motions of Dihydrofolate Reductase

Many enzymes are known to change conformations during their catalytic cycle, but the role of these protein motions is not well understood. Escherichia coli dihydrofolate reductase (DHFR) is a small, flexible enzyme that is often used as a model system for understanding enzyme dynamics. Recently, native tryptophan fluorescence was used as a probe to study micro- to millisecond dynamics of DHFR. Yet, because DHFR has five native tryptophans, the origin of the observed conformational changes could not be assigned to a specific region within the enzyme. Here, we use DHFR mutants, each with a single tryptophan as a probe for temperature jump fluorescence spectroscopy, to further inform our understanding of DHFR dynamics. The equilibrium tryptophan fluorescence of the mutants shows that each tryptophan is in a different environment and that wild-type DHFR fluorescence is not a simple summation of all the individual tryptophan fluorescence signatures due to tryptophan–tryptophan interactions. Additionally, each mutant exhibits a two-phase relaxation profile corresponding to ligand association/dissociation convolved with associated conformational changes and a slow conformational change that is independent of ligand association and dissociation, similar to the wild-type enzyme. However, the relaxation rate of the slow phase depends on the location of the tryptophan within the enzyme, supporting the conclusion that the individual tryptophan fluorescence dynamics do not originate from a single collective motion, but instead report on local motions throughout the enzyme.

Crowders Steal Dihydrofolate Reductase Ligands through Quinary Interactions

Dihydrofolate reductase (DHFR) reduces dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. Due to its role in one carbon metabolism, chromosomal DHFR is the target of the antibacterial drug, trimethoprim. Resistance to trimethoprim has resulted in a type II DHFR that is not structurally related to the chromosomal enzyme target. Because of its metabolic significance, understanding DHFR kinetics and ligand binding behavior in more cell-like conditions, where the total macromolecule concentration can be as great as 300 mg/mL, is important. The progress-curve kinetics and ligand binding properties of the drug target (chromosomal E. coli DHFR) and the drug resistant (R67 DHFR) enzymes were studied in the presence of macromolecular cosolutes. There were varied effects on NADPH oxidation and binding to the two DHFRs, with some cosolutes increasing affinity and others weakening binding. However, DHF binding and reduction in both DHFRs decreased in the presence of all cosolutes. The decreased binding of ligands is mostly attributed to weak associations with the macromolecules, as opposed to crowder effects on the DHFRs. Computer simulations found weak, transient interactions for both ligands with several proteins. The net charge of protein cosolutes correlated with effects on NADP⁺ binding, with near neutral and positively charged proteins having more detrimental effects on binding. For DHF binding, effects correlated more with the size of binding pockets on the protein crowders. These nonspecific interactions between DHFR ligands and proteins predict that the in vivo efficiency of DHFRs may be much lower than expected from their in vitro rates.

Effect of Asp122 Mutation on the Hydride Transfer in E. coli DHFR Demonstrates the Goldilocks of Enzyme Flexibility

Dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate (THF) in the presence of NADPH. The key hydride transfer step in the reaction is facilitated by a combination of enzyme active site preorganization and correlated protein motions in the Michaelis-Menten (E:NADPH:DHF) complex. The present theoretical study employs mutagenesis to examine the relation between structural and functional properties of the enzyme. We mutate Asp122 in Escherichia coli DHFR, which is a conserved amino acid in the DHFR family. The consequent effect of the mutation on enzyme catalysis is examined from an energetic, structural and short-time dynamic perspective. Our investigations suggest that the structural and short-time dynamic perturbations caused by Asp122X mutations (X = Asn, Ser, Ala) are along the reaction coordinate and lower the rate of hydride transfer. Importantly, analysis of the correlated and principle component motions in the enzyme suggest that the mutation alters the coupled motions that are present in the wild-type enzyme. In the case of D122N and D122S, the mutations inhibit coupled motion, whereas in the case of D122A, the mutation enhances coupled motion, although all mutations result in similar rate reduction. These results emphasize a Goldilocks principle of enzyme flexibility, that is, enzymes should neither be too rigid nor too flexible.

Reduction of Folate by Dihydrofolate Reductase from Thermotoga maritima

Mammalian dihydrofolate reductases (DHFRs) catalyze the reduction of folate more efficiently than the equivalent bacterial enzymes do, despite typically having similar efficiencies for the reduction of their natural substrate, dihydrofolate. In contrast, we show here that DHFR from the hyperthermophilic bacterium Thermotoga maritima can catalyze reduction of folate to tetrahydrofolate with an efficiency similar to that of reduction of dihydrofolate under saturating conditions. Nuclear magnetic resonance and mass spectrometry experiments showed no evidence of the production of free dihydrofolate during either the EcDHFR- or TmDHFR-catalyzed reductions of folate, suggesting that both enzymes perform the two reduction steps without release of the partially reduced substrate. Our results imply that the reaction proceeds more efficiently in TmDHFR than in EcDHFR because the more open active site of TmDHFR facilitates protonation of folate. Because T. maritima lives under extreme conditions where tetrahydrofolate is particularly prone to oxidation, this ability to salvage folate may impart an advantage to the bacterium by minimizing the squandering of a valuable cofactor.

Applications of Normal Mode Analysis Methods in Computational Protein Design

Recent advances in coarse-grained normal mode analysis methods make possible the large-scale prediction of the effect of mutations on protein stability and dynamics as well as the generation of biologically relevant conformational ensembles. Given the interplay between flexibility and enzymatic activity, the combined analysis of stability and dynamics using the Elastic Network Contact Model (ENCoM) method has ample applications in protein engineering in industrial and medical applications such as in computational antibody design. Here, we present a detailed tutorial on how to perform such calculations using ENCoM.

Protein motions and dynamic effects in enzyme catalysis

The role of protein motions in promoting the chemical step of enzyme catalysed reactions remains a subject of considerable debate. Here, a unified view of the role of protein dynamics in dihydrofolate reductase catalysis is described. Recently the role of such motions has been investigated by characterising the biophysical properties of isotopically substituted enzymes through a combination of experimental and computational analyses. Together with previous work, these results suggest that dynamic coupling to the chemical coordinate is detrimental to catalysis and may have been selected against during DHFR evolution. The full catalytic power of Nature's catalysts appears to depend on finely tuning protein motions in each step of the catalytic cycle.

Ligand-Dependent Conformational Dynamics of Dihydrofolate Reductase

Enzymes are known to change among several conformational states during turnover. The role of such dynamic structural changes in catalysis is not fully understood. The influence of dynamics in catalysis can be inferred, but not proven, by comparison of equilibrium structures of protein variants and protein-ligand complexes. A more direct way to establish connections between protein dynamics and the catalytic cycle is to probe the kinetics of specific protein motions in comparison to progress along the reaction coordinate. We have examined the enzyme model system dihydrofolate reductase (DHFR) from Escherichia coli with tryptophan fluorescence-probed temperature-jump spectroscopy. We aimed to observe the kinetics of the ligand binding and ligand-induced conformational changes of three DHFR complexes to establish the relationship among these catalytic steps. Surprisingly, in all three complexes, the observed kinetics do not match a simple sequential two-step process. Through analysis of the relationship between ligand concentration and observed rate, we conclude that the observed kinetics correspond to the ligand binding step of the reaction and a noncoupled enzyme conformational change. The kinetics of the conformational change vary with the ligand's identity and presence but do not appear to be directly related to progress along the reaction coordinate. These results emphasize the need for kinetic studies of DHFR with highly specific spectroscopic probes to determine which dynamic events are coupled to the catalytic cycle and which are not.

Structure of a thermostable methionine adenosyltransferase from Thermus thermophilus HB27 reveals a novel fold of the flexible loop

Methionine adenosyltransferases (MATs) are the family of enzymes which synthesize S-adenosylmethionine (AdoMet), the major biological methyl donor.

Keep on moving: discovering and perturbing the conformational dynamics of enzymes

CONSPECTUS: Because living organisms are in constant motion, the word "dynamics" can hold many meanings to biologists. Here we focus specifically on the conformational changes that occur in proteins and how studying these protein dynamics may provide insights into enzymatic catalysis. Advances in integrating techniques such as X-ray crystallography, nuclear magnetic resonance, and electron cryomicroscopy (cryo EM) allow us to model the dominant structures and exchange rates for many proteins and protein complexes. For proteins amenable to atomic resolution techniques, the major questions shift from simply describing the motions to discovering their role in function. Concurrently, there is an increasing need for using perturbations to test predictive models of dynamics-function relationships. Examples are the catalytic cycles of dihydrofolate reductase (DHFR) and cyclophilin A (CypA). In DHFR, mutations that alter the ability of the active site to sample productive higher energy states on the millisecond time scale reduce the rate of hydride transfer significantly. Recently identified rescue mutations restore function, but the mechanism by which they do so remains unclear. The exact role of any changes in the dynamics remains an open question. For CypA, a network of side chains that exchange between two conformations is critical for catalysis. Mutations that lock the network in one state also reduce catalytic activity. A further understanding of enzyme dynamics of well-studied enzymes such as dihydrofolate reductase and cyclophilin A will lead to improvement in ability to modulate the functions of computationally designed enzymes and large macromolecular machines. In designed enzymes, directed evolution experiments increase catalytic rates. Detailed X-ray studies suggest that these mutations likely limit the conformational space explored by residues in the active site. For proteins where atomic resolution information is currently inaccessible, other techniques such as cryo-EM and high-resolution single molecule microscopy continue to advance. Understanding the conformational dynamics of larger systems such as protein machines will likely become more accessible and provide new opportunities to rationally modulate protein function.

Loop interactions during catalysis by dihydrofolate reductase from Moritella profunda

Dihydrofolate reductase (DHFR) is often used as a model system to study the relation between protein dynamics and catalysis. We have studied a number of variants of the cold-adapted DHFR from Moritella profunda (MpDHFR), in which the catalytically important M20 and FG loops have been altered, and present a comparison with the corresponding variants of the well-studied DHFR from Escherichia coli (EcDHFR). Mutations in the M20 loop do not affect the actual chemical step of transfer of hydride from reduced nicotinamide adenine dinucleotide phosphate to the substrate 7,8-dihydrofolate in the catalytic cycle in either enzyme; they affect the steady state turnover rate in EcDHFR but not in MpDHFR. Mutations in the FG loop also have different effects on catalysis by the two DHFRs. Despite the two enzymes most likely sharing a common catalytic cycle at pH 7, motions of these loops, known to be important for progression through the catalytic cycle in EcDHFR, appear not to play a significant role in MpDHFR.

A coarse-grained elastic network atom contact model and its use in the simulation of protein dynamics and the prediction of the effect of mutations

Normal mode analysis (NMA) methods are widely used to study dynamic aspects of protein structures. Two critical components of NMA methods are coarse-graining in the level of simplification used to represent protein structures and the choice of potential energy functional form. There is a trade-off between speed and accuracy in different choices. In one extreme one finds accurate but slow molecular-dynamics based methods with all-atom representations and detailed atom potentials. On the other extreme, fast elastic network model (ENM) methods with Cα-only representations and simplified potentials that based on geometry alone, thus oblivious to protein sequence. Here we present ENCoM, an Elastic Network Contact Model that employs a potential energy function that includes a pairwise atom-type non-bonded interaction term and thus makes it possible to consider the effect of the specific nature of amino-acids on dynamics within the context of NMA. ENCoM is as fast as existing ENM methods and outperforms such methods in the generation of conformational ensembles. Here we introduce a new application for NMA methods with the use of ENCoM in the prediction of the effect of mutations on protein stability. While existing methods are based on machine learning or enthalpic considerations, the use of ENCoM, based on vibrational normal modes, is based on entropic considerations. This represents a novel area of application for NMA methods and a novel approach for the prediction of the effect of mutations. We compare ENCoM to a large number of methods in terms of accuracy and self-consistency. We show that the accuracy of ENCoM is comparable to that of the best existing methods. We show that existing methods are biased towards the prediction of destabilizing mutations and that ENCoM is less biased at predicting stabilizing mutations.

Solvent environments significantly affect the enzymatic function of Escherichia coli dihydrofolate reductase: comparison of wild-type protein and active-site mutant D27E

To investigate the contribution of solvent environments to the enzymatic function of Escherichia coli dihydrofolate reductase (DHFR), the salt-, pH-, and pressure-dependence of the enzymatic function of the wild-type protein were compared with those of the active-site mutant D27E in relation to their structure and stability. The salt concentration-dependence of enzymatic activity indicated that inorganic cations bound to and inhibited the activity of wild-type DHFR at neutral pH. The BaCl2 concentration-dependence of the (1)H-(15)N HSQC spectra of the wild-type DHFR-folate binary complex showed that the cation-binding site was located adjacent to the Met20 loop. The insensitivity of the D27E mutant to univalent cations, the decreased optimal pH for its enzymatic activity, and the increased Km and Kd values for its substrate dihydrofolate suggested that the substrate-binding cleft of the mutant was slightly opened to expose the active-site side chain to the solvent. The marginally increased fluorescence intensity and decreased volume change due to unfolding of the mutant also supported this structural change or the modified cavity and hydration. Surprisingly, the enzymatic activity of the mutant increased with pressurization up to 250MPa together with negative activation volumes of -4.0 or -4.8mL/mol, depending on the solvent system, while that of the wild-type was decreased and had positive activation volumes of 6.1 or 7.7mL/mol. These results clearly indicate that the insertion of a single methylene at the active site could substantially change the enzymatic reaction mechanism of DHFR, and solvent environments play important roles in the function of this enzyme.

A distal mutation perturbs dynamic amino acid networks in dihydrofolate reductase

Correlated networks of amino acids have been proposed to play a fundamental role in allostery and enzyme catalysis. These networks of amino acids can be traced from surface-exposed residues all the way into the active site, and disruption of these networks can decrease enzyme activity. Substitution of the distal Gly121 residue in Escherichia coli dihydrofolate reductase results in an up to 200-fold decrease in the hydride transfer rate despite the fact that the residue is located 15 Å from the active-site center. In this study, nuclear magnetic resonance relaxation experiments are used to demonstrate that dynamics on the picosecond to nanosecond and microsecond to millisecond time scales are changed significantly in the G121V mutant of dihydrofolate reductase. In particular, picosecond to nanosecond time scale dynamics are decreased in the FG loop (containing the mutated residue at position 121) and the neighboring active-site loop (the Met20 loop) in the mutant compared to those of the wild-type enzyme, suggesting that these loops are dynamically coupled. Changes in methyl order parameters reveal a pathway by which dynamic perturbations can be propagated more than 25 Å across the protein from the site of mutation. All of the enzyme complexes, including the model Michaelis complex with folate and nicotinamide adenine dinucleotide phosphate bound, assume an occluded ground-state conformation, and we do not observe sampling of a higher-energy closed conformation by (15)N R2 relaxation dispersion experiments. This is highly significant, because it is only in the closed conformation that the cofactor and substrate reactive centers are positioned for reaction. The mutation also impairs microsecond to millisecond time scale fluctuations that have been implicated in the release of product from the wild-type enzyme. Our results are consistent with an important role for Gly121 in controlling protein dynamics critical for enzyme function and further validate the dynamic energy landscape hypothesis of enzyme catalysis.

Coupling effects of distal loops on structural stability and enzymatic activity of Escherichia coli dihydrofolate reductase revealed by deletion mutants

Residues distal from the active site in dihydrofolate reductase (DHFR) have regulatory roles in catalytic reaction and also folding stability. The couplings of the distal residues to the ones in the active site have been analyzed using site-directed mutants. To expand our understanding of the structural and functional influences of distal residue mutation, we explored the structural stability and enzymatic activity of deletion mutants. Deletion has greater structural and dynamical impacts on the corresponding part than site-directed mutation does. Thus, deletion amplifies the effects caused by distal mutations, which should make the mutual couplings among the distant residues more apparent. We focused on residues 52, 67, 121, and 145 in the four distinct loops of DHFR. All the single-residue deletion mutants showed marked reduction in stability, except for Delta52 in an alphaC-betaC loop. Double deletion mutants showed that the loop alphaC-betaC has nonadditive couplings with the betaF-betaG and betaG-betaH loops regarding stability. Single deletion to the loops alphaC-betaC or betaC-betaD resulted in considerable activity reduction, demonstrating that the loops couple to the residues near the active site. The four loops were shown to be functionally interdependent from the double deletion experiments.

Ensemble-based signatures of energy propagation in proteins: a new view of an old phenomenon

The ability of a protein to transmit the energetic effects of binding from one site to another constitutes the underlying basis for allosterism and signal transduction. Despite clear experimental evidence indicating the ability of proteins to transmit the effects of binding, the means by which this propagation is facilitated is not well understood. Using our previously developed ensemble-based description of the equilibrium, we investigated the physical basis of energy propagation and identified several fundamental and general aspects of energetic coupling between residues in a protein. First, partitioning of a conformational ensemble into four distinct sub-ensembles allows for explanation of the range of experimentally observed coupling behaviors (i.e., positive, neutral, and negative coupling between various regions of the protein structure). Second, the relative thermodynamic properties of these four sub-ensembles define the energetic coupling between residues as either positive, neutral, or negative. Third, analysis of the structural and thermodynamic features of the states within each sub-ensemble reveals significant variability. This third result suggests that a quantitative description of energy propagation in proteins requires an understanding of the structural and energetic features of more than just one or a few low-energy states, but also of many high-energy states. Such findings illuminate the difficulty in interpreting energy propagation in proteins in terms of a structural pathway that physically links coupled sites.

Effects of Mutation at Methionine-42 of Escherichia coli Dihydrofolate Reductase on Stability and Function: Implication of Hydrophobic Interactions

Methionine-42, distal to the active site of Escherichia coli dihydrofolate reductase, was substituted by site-directed mutagenesis with 14 amino acids (Ala, Cys, Glu, Gln, Gly, His, Ile, Leu, Pro, Ser, Thr, Trp, Tyr, and Val) to elucidate its role in the stability and function of this enzyme. Far-ultraviolet circular dichroism spectra of these mutants showed a distinctive negative peak at around 230 nm beside 220 nm, depending on the hydrophobicity of the amino acids introduced. The fluorescence intensity also increased in an order similar to that of the amino acids. These spectroscopic data suggest that the mutations do not affect the secondary structure, but strongly perturb the exciton coupling between Trp47 and Trp74. The free energy of urea unfolding, deltaG(o)u, increased with increases in the side-chain hydrophobicity in the range 2.96-6.40 kcal x mol(-1), which includes the value for the wild-type enzyme (6.08 kcal x mol(-1)). The steady-state kinetic parameters, Km and kcat, also increased with increases in the side-chain hydrophobicity, with the M42W mutant showing the largest increases in Km (35-fold) and kcat (4.3-fold) compared with the wild-type enzyme. These results demonstrate that site 42 distal to the active site plays an important role in the stability and function of this enzyme, and that the main effect of the mutations is to modify of hydrophobic interactions with the residues surrounding this position.

Peroxynitrite-mediated oxidation of the C85S/C152E mutant of dihydrofolate reductase from Escherichia coli: functional and structural effects

Peroxynitrite is a potent reactive oxygen species that is believed to mediate deleterious protein modifications in a wide variety of neurodegenerative disorders. In this study, we have analysed the effects of oxidative damage induced by peroxynitrite on a cysteine-free mutant of dihydrofolate reductase (SE-DHFR), from a functional and a structural point of view. The peroxynitrite-mediated oxidation results in the inhibition, concentration-dependent, of the catalytic activity. This effect is strongly influenced by the HCO(3)(-)/CO(2) buffering system, that we observed to significantly affect the yield of protein oxidation by modulating the peroxynitrite-induced modification of aromatic residues. Because of this effect, in presence of bicarbonate system, we have observed a protection of enzymatic activity of SE-DHFR with regard to peroxynitrite. The thermodynamic stability of the oxidized protein has been studied in comparison with the non-oxidized protein by differential scanning calorimetry. The thermodynamic parameters obtained showed a decrease of stability of SE-DHFR upon oxidation, evaluated in terms of Gibbs free energy of about 1.25 kcal/mol at 25 degrees C, with respect to the non-oxidized protein. Together, these data indicate that structural and functional alterations induced by peroxynitrite may play a direct role in compromising DHFR function in multiple pathological conditions.

Pivotal Role of Gly 121 in Dihydrofolate Reductase from Escherichia coli: The Altered Structure of a Mutant Enzyme May Form the Basis of Its Diminished Catalytic Performance

The structure and folding of dihydrofolate reductase (DHFR) from Escherichia coli and the mutant G121V-DHFR, in which glycine 121 in the exterior FG loop was replaced with valine, were studied by molecular dynamics simulations and CD and fluorescence spectroscopy. The importance of residue 121 for the chemical step during DHFR catalysis had been demonstrated previously. High-temperature MD simulations indicated that while DHFR and G121V-DHFR followed similar unfolding pathways, the strong contacts between the M20 loop and the FG loop in DHFR were less stable in the mutant. These contacts have been proposed to be involved in a coupled network of interactions that influence the protein dynamics and promote catalysis [Benkovic, S. J., and Hammes-Schiffer, S. (2003) Science 301, 1196-1202]. CD spectroscopy of DHFR and G121V-DHFR indicated that the two proteins existed in different conformations at room temperature. While the thermally induced unfolding of DHFR was highly cooperative with a midpoint at 51.6 +/- 0.7 degrees C, G121V-DHFR exhibited a gradual decrease in its level of secondary structure without a clear melting temperature. Temperature-induced unfolding and renaturation from the urea-denatured state revealed that both proteins folded via highly fluorescent intermediates. The formation of these intermediates occurred with relaxation times of 149 +/- 4.5 and 256 +/- 13 ms for DHFR and G121V-DHFR, respectively. The fluorescence intensity for the intermediates formed during refolding of G121V-DHFR was approximately twice that of the wild-type. While the fluorescence intensity then slowly decayed for DHFR toward a state representing the native protein, G121V-DHFR appeared to be trapped in a highly fluorescent state. These results suggest that the reduced catalytic activity of G121V-DHFR is the consequence of nonlocal structural effects that may result in a perturbation of the network of promoting motions.

Biochemical characteristics of guanine nucleotide binding protein alpha-subunit recombinant protein and three mutants: investigation of a domain motion involved in GDP-GTP exchange

Our previous studies using molecular dynamics have shown a hinge bending motion between the helical and the GTPase domains of GalphaT (Mello et al., 1998). The hypothesis that this motion is allowed by residues Gly56 and Gly179 and that this motion may affect the ligand exchange was tested in this work. Mutations of Gly 56 were carried out and the mutant proteins were expressed in Sf9 cells using the Baculovirus expression system. The recombinant proteins were purified using Ni-NTA affinity chromatography. The results for the (GDP/GTP) exchange assays showed that G56S and double mutants (D55G/G56S) proteins differ significantly from the wild type and D55G mutant forms. The Kd values for GTPgammaS binding of those mutants have decreased by approximately 10-fold. No difference in the GTPase activity was detected for the mutants. Thus, the biochemical results obtained support the conclusions of the computational studies.

The pepsin residue glycine-76 contributes to active-site loop flexibility and participates in catalysis

Glycine residues are known to contribute to conformational flexibility of polypeptide chains, and have been found to contribute to flexibility of some loops associated with enzymic catalysis. A comparison of porcine pepsin in zymogen, mature and inhibited forms revealed that a loop (a flap), consisting of residues 71--80, located near the active site changed its position upon substrate binding. The loop residue, glycine-76, has been implicated in the catalytic process and thought to participate in a hydrogen-bond network aligning the substrate. This study investigated the role of glycine-76 using site-directed mutagenesis. Three mutants, G76A, G76V and G76S, were constructed to increase conformational restriction of a polypeptide chain. In addition, the serine mutant introduced a hydrogen-bonding potential at position 76 similar to that observed in human renin. All the mutants, regardless of amino acid size and polarity, had lower catalytic efficiency and activated more slowly than the wild-type enzyme. The slower activation process was associated directly with altered proteolytic activity. Consequently, it was proposed that a proteolytic cleavage represents a limiting step of the activation process. Lower catalytic efficiency of the mutants was explained as a decrease in the flap flexibility and, therefore, a different pattern of hydrogen bonds responsible for substrate alignment and flap conformation. The results demonstrated that flap flexibility is essential for efficient catalytic and activation processes.

Effect of ligand binding on the flexibility of dihydrofolate reductase as revealed by compressibility

The partial specific volume, v, and adiabatic compressibility, beta(s), of Escherichia coli dihydrofolate reductase were measured at 30 degrees C in the presence of various ligands (folate, dihydrofolate, tetrahydrofolate, NADPH, NADP, methotrexate, and KCl). Binding of these ligands (binary and ternary complexes) brought about large changes of v (0.734-0.754 cm(3) g(-1)) and beta(s) (6. 6x10(-6)-9.8x10(-6) bar(-1)), keeping a linear relationship between the two parameters. The values of v and beta(s) increased with an increase in internal cavity, V(cav), and a decrease in accessible surface area, ASA, which were calculated from the X-ray crystal structures of the complexes. A large variation of V(cav) relative to ASA by ligand binding suggested that the cavity is a dominant factor and the effect of hydration might be small for the ligand-induced changes of v and beta(s). The beta(s) values of the binary and ternary complexes suggested a characteristic conformational flexibility of the kinetic intermediates in the enzyme reaction coordinate. Comparison of beta(s) with the cavity distribution in the crystal structures revealed that the flexibility of the intermediates was mainly determined by the total cavity volume with minor contributions of the number, position, and size of cavities. These results demonstrate that the compressibility is a useful measure of the conformational flexibility of the intermediates in the enzyme reaction and that the combined study of compressibility and X-ray crystallography gives new insight into the protein dynamics through the behavior of the cavities.

A test of the relationship between sequence and structure in proteins: excision of the heme binding site in apocytochrome b5

The water-soluble domain of rat hepatic holocytochrome b5 is an alphabeta protein containing elements of secondary structure in the sequence beta1-alpha1-beta4-beta3-alpha2-alpha3-beta5- alpha4-alpha5-beta2-alpha6. The heme group is enclosed by four helices, a2, a3, a4, and a5. To test the hypothesis that a small b hemoprotein can be constructed in two parts, one forming the heme site, the other an organizing scaffold, a protein fragment corresponding to beta1-alpha1-beta4-beta3-lambda-beta2-alpha6 was prepared, where lambda is a seven-residue linker bypassing the heme binding site. The fragment ("abridged b5") was found to contain alpha and beta secondary structure by circular dichroism spectroscopy and tertiary structure by Trp fluorescence emission spectroscopy. NMR data revealed a species with spectral properties similar to those of the full-length apoprotein. This folded form is in slow equilibrium on the chemical shift time scale with other less folded species. Thermal denaturation, as monitored by circular dichroism, absorption, and fluorescence spectroscopy, as well as size-exclusion chromatography-fast protein liquid chromatography (SEC-FPLC), confirmed the coexistence of at least two distinct conformational ensembles. It was concluded that the protein fragment is capable of adopting a specific fold likely related to that of cytochrome b5, but does not achieve high thermodynamic stability and cooperativity. Abridged b5 demonstrates that the spliced sequence contains the information necessary to fold the protein. It suggests that the dominating influence to restrict the conformational space searched by the chain is structural propensities at a local level rather than internal packing. The sequence also holds the properties necessary to generate a barrier to unfolding.

Deletion of a highly motional residue affects formation of the Michaelis complex for Escherichia coli dihydrofolate reductase

Analysis of the dihydrofolate reductase (DHFR) complex with folate by two-dimensional heteronuclear (1H-15N) nuclear magnetic relaxation revealed that isolated residues exhibit diverse backbone fluctuations on the nanosecond to picosecond time scale [Epstein, D. M., Benkovic, S. J., and Wright, P. E. (1995) Biochemistry 34, 11037-11048]. These dynamical features may be significant in forming the Michaelis complex. Of these residues, glycine 121 displays large-amplitude backbone motions on the nanosecond time scale. This amino acid, strictly conserved for prokaryotic DHFRs, is located at the center of the betaF-betaG loop. To investigate the catalytic importance of this residue, we report the effects of Gly121 deletion and glycine insertion into the modified betaF-betaG loop. Relative to wild type, deletion of Gly121 dramatically decreases the rate of hydride transfer 550-fold and the strength of cofactor binding 20-fold for NADPH and 7-fold for NADP+. Furthermore, DeltaG121 DHFR requires conformational changes dependent on the initial binary complex to attain the Michaelis complex poised for hydride transfer. Surprisingly, the insertion mutants displayed a significant decrease in both substrate and cofactor binding. The introduction of glycine into the modified betaF-betaG loop, however, generally eliminated conformational changes required by DeltaG121 DHFR to attain the Michaelis complex. Taken together, these results suggest that the catalytic role for the betaF-betaG loop includes formation of liganded complexes and proper orientation of substrate and cofactor. Through a transient interaction with the Met20 loop, alterations of the betaF-betaG loop can orchestrate proximal and distal effects on binding and catalysis that implicate a variety of enzyme conformations participating in the catalytic cycle.

Mutagenesis of histidine 26 demonstrates the importance of loop-loop and loop-protein interactions for the function of iso-1-cytochrome c

In yeast iso-1-cytochrome c, the side chain of histidine 26 (His26) attaches omega loop A to the main body of the protein by forming a hydrogen bond to the backbone atom carbonyl of glutamic acid 44. The His26 side chain also forms a stabilizing intra-loop interaction through a hydrogen bond to the backbone amide of asparagine 31. To investigate the importance of loop-protein attachment and intra-loop interactions to the structure and function of this protein, a series of site-directed and random-directed mutations were produced at His26. Yeast strains expressing these variant proteins were analyzed for their ability to grow on non-fermentable carbon sources and for their intracellular production of cytochrome c. While the data show that mutations at His26 lead to slightly decreased intracellular amounts of cytochrome c, the level of cytochrome c function is decreased more. The data suggest that cytochrome c reductase binding is affected more than cytochrome c oxidase or lactate dehydrogenase binding. We propose that mutations at this residue increase loop mobility, which, in turn, decreases the protein's ability to bind redox partners.

Effect of limited proteolysis on the stability and enzymatic activity of human placental S-adenosylhomocysteine hydrolase

Human placental S-adenosylhomocysteine (AdoHcy) hydrolase was subjected to limited papain digestion. The multiple cleavage sites in the enzyme were identified to be Lys94-Ala95, Tyr100-Ala101, Glu243-Ile244, Met367-Ala368, Gln369-Ile370, and Gly382-Val383. Despite multiple cleavage sites in the backbone of the protein, the digested enzyme was able to maintain its quaternary structure and retain its full catalytic activity. The enzyme activity of the partially digested AdoHcy hydrolase was essentially identical to that of the native enzyme at several pH values. The thermal stabilities of the native and partially digested enzymes were only slightly different at all temperatures tested. The stability of both native and partially digested enzymes were examined in guanidine hydrochloride and equilibrium unfolding transitions were monitored by CD spectroscopy and tryptophan fluorescence spectroscopy. The results of these experiments can be summarized as follows: (1) CD spectroscopic analysis showed that the overall secondary and tertiary structures of the partially digested enzyme are essentially identical with those of the native enzyme; and (2) tryptophan fluorescence spectroscopic analysis indicated that there are small differences in the environments of surface-exposed tryptophan residues between the partially digested enzyme and the native enzyme under unfolding conditions. The differences in the free energy of unfolding, delta(delta Gu) [delta Gu(native)-delta Gu(digested)], is approximately 1.3 kcal/mol. When NAD+ was removed from the partially digested enzyme, the secondary and tertiary structures of the apo form of the digested AdoHcy hydrolase were completely lost and the enzymatic activity could not be recovered by incubation with excess NAD+. These results suggest that AdoHcy hydrolase exists as a very compact enzyme with extensive intramolecular bonding, which contributes significantly to the overall global protein stabilization. Identification of the surface-exposed peptide bonds, which are susceptible to papain digestion, has provided some constraints on the spatial orientations of subunits of the enzyme. This information, in turn, has provided supplemental data for X-ray crystallographic studies currently ongoing in our laboratories.

Analysis of the structure and stability of omega loop A replacements in yeast iso-1-cytochrome c

Omega (omega)-loop A, residues 18-32 in wild-type yeast iso-1-cytochrome c, has been deleted and replaced with loop sequences from three other cytochromes c and one from esterase. Yeast expressing a partial loop deletion do not contain perceptible amounts of holoprotein as measured by low-temperature spectroscopy and cannot grow on nonfermentable media. Strains expressing loop replacement mutations accumulate holoprotein in vivo, but the protein function varies depending on the sequence and length of the replacement loop; in vivo expression levels do not correlate with their thermal denaturation temperatures. In vitro spectroscopic studies of the loop replacement proteins indicate that all fold into a native-like cytochrome c conformation, but are less stable than the wild-type protein. Decreases in thermal stability are caused by perturbation of loop C backbone in one case and a slight reorganization of the protein hydrophobic core in another case, rather than rearrangement of the loop A backbone. A single-site mutation in one of the replacement mutants designed to relieve inefficient hydrophobic core packing caused by the new loop recovers some, but not all, of the lost stability.

Flexible loop in the structure of S-adenosylmethionine synthetase crystallized in the tetragonal modification

S-Adenosylmethionine synthetase (MAT, ATP:L-methionine S-adenosyltransferase, E.C.2.5.1.6.) plays a central metabolic role in all organisms. MAT catalyzes the two-step reaction which synthesizes S-adenosylmethionine (AdoMet), pyrophosphate (PPi) and orthophosphate (Pi) from ATP and L-methionine. AdoMet is the primary methyl group donor in biological systems. MAT from Escherichia coli was crystallized in the tetragonal modification with space group P4(3)2(1)2 using the same conditions as previously yielded crystals of the hexagonal system [Takusagawa, et al., (1996), J. Biol. Chem. 171, 136-147], except for the crystallization temperature. The structure has been determined by molecular replacement at 3.2 A resolution. The overall structure of the tetrameric MAT in the tetragonal modification is essentially the same as the structure found in the hexagonal modification. However there are two remarkable differences between the structures of two modifications. One is the contents in the active sites (holoform vs. apo-form), and the other is the conformation of the flexible loop over the active site (open vs. closed). These differences in the crystal structures are caused solely by the difference in crystallization temperatures (26 degrees C vs. 4 degrees C). We have interpreted the structural data obtained from the X-ray analyses in conjunction with the results of the mechanistic and sequencing studies in terms of possible dynamic motion of the flexible loop. When a substrate/product binds in the active site (hexagonal modification), the loop becomes disordered, apparently due to flexibility at the entrance of the active site as if it acts as a "mobile loop" during the catalytic reaction. On the other hand, when the temperature is decreased, the dynamic motion of the flexible loop may be reduced, and the loop residues enter the active site and close its entrance (tetragonal modification). Thus, the active site of the tetragonal modification is empty despite the crystals being grown in mother liquor containing a large concentration of phosphate (100 mM). There is no significant displacement of amino acid residues in the active site between the holo and apo forms, suggesting that the flexible loop plays an important role in determination of the contents in the active site. Since the functionally important amino acid residues in the active site are all conserved throughout various species, the structures of the active sites and the mechanism of the catalysis are probably essentially identical in the enzymes from a wide range of organisms. However, the substrate KM and Vmax values of MATs from various species are distributed over a wide range. The amino acid residues in the flexible loop regions are poorly conserved throughout various species. Therefore, the wide differences in catalysis rates of MATs from various speeches may be due to the differences in the composition of the flexible loop.

Thermal denaturation of bacterial and bovine dihydrofolate reductases and their complexes with NADPH, trimethoprim and methotrexate

Scanning microcalorimetry was used for the study of thermal denaturation of E.coli and bovine liver dihydrofolate reductases (cDHFR and bDHFR, respectively) and their complexes with NADPH, trimethoprim (TMP) and methotrexate (MTX) at pH 6.8. It was shown that the denaturation temperature of bDHFR is 7.2 degrees C less than that of cDHFR and that ionic strength is equally important for the thermostability and cooperativity of the denaturation process of the two proteins. Binding of antifolate compounds significantly stabilizes DHFR against heat denaturation. The stabilizing effect and the transition cooperativity depend on the nature of the inhibitor, the presence of NADPH and the origin of the enzyme. The dependence of calorimetric denaturation enthalpy (calculated per gram of protein) on denaturation temperature for DHFRs, their complexes with NADPH and binary/ternary complexes with TMP/MTX fits to the same straight line with the slope of 0.66 J/Kg. This relatively high value indicates an essential role of hydrophobic contacts in the stabilization of DHFR structure. The change of denaturation temperatures in binary complexes with MTX/TMP (in comparison with the free enzymes) is as much as 14.2 degrees C/8.5 degrees C and 13.3 degrees C/3.2 degrees C for cDHFR and bDHFR, respectively. The same change in ternary complexes with MTX/TMP is much more pronounced and equals to 21.9 degrees C/16.8 degrees C and 29.0 degrees C/16.4 degrees C. The vast difference of binary and ternary complexes thermostability demonstrates the important role of cofactor in the stabilization of enzyme. Moving from binary to ternary systems causes a significant increase in denaturation temperatures, even when corresponding association constants do not change (cDHFR binary/ternary complexes with MTX) or increases very slightly (bDHFR binary/ternary complexes with TMP). In all other cases the increase of denaturation temperature for each protein in complex with ligands correlates with the association constant for the corresponding complex.

Site-directed mutagenesis studies of the high-affinity streptavidin-biotin complex: contributions of tryptophan residues 79, 108, and 120

We report the functional characterization of site-directed biotin binding-site mutants of recombinant core streptavidin. The mutagenesis studies were aimed at characterizing the contributions of Trp residues known to contact biotin that have been postulated to control the exceptional binding affinity observed in this system. The functional properties of single site-directed mutants replacing Trp residues with Phe or Ala at positions 79, 108, and 120 were investigated by quantitating the EC50 binding parameters of these mutants to biotin and 2-iminobiotin in an ELISA format. The biotin EC50 for all mutants was the same as wild-type streptavidin, demonstrating that their delta Ka values relative to wild type were < 10(6). The conservative W79F and W108F mutants displayed only a 2- to 3-fold increase in EC50 for 2-iminobiotin, corresponding to an estimated delta Ka < 10, while the W120F mutant displayed a much greater alteration in 2-iminobiotin EC50, corresponding to an estimated delta Ka of 10(2). These delta Ka values are likely to reflect similar changes for biotin. The 2-iminobiotin EC50 values for the Ala mutants fell outside the accessible concentration range of the ELISA assay, demonstrating that these mutations lowered the Ka by a factor of 10(4) to 10(6). Direct estimation of biotin Ka values for W79A, W120A, and W120F in an ultrafiltration binding assay yielded Ka values of 4.3 x 10(7) M-1, 8.6 x 10(6) M-1, and > 5 x 10(9) M-1, respectively, in excellent agreement with the ELISA estimates of delta Ka with 2-iminobiotin as a reporter ligand. The results of these preliminary functional studies suggest that these aromatic side chains contribute significantly to the streptavidin-biotin binding free energy.