Analysis of hydride transfer and cofactor fluorescence decay in mutants of dihydrofolate reductase: possible evidence for participation of enzyme molecular motions in catalysis

Modulating Enzyme Function via Dynamic Allostery within Biliverdin Reductase B

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic "dials" to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.

Biliverdin Reductase B Dynamics Are Coupled to Coenzyme Binding

Biliverdin reductase B (BLVRB) is a newly identified cellular redox regulator that catalyzes the NADPH-dependent reduction of multiple substrates. Through mass spectrometry analysis, we identified high levels of BLVRB in mature red blood cells, highlighting the importance of BLVRB in redox regulation. The BLVRB conformational changes that occur during conezyme/substrate binding and the role of dynamics in BLVRB function, however, remain unknown. Through a combination of NMR, kinetics, and isothermal titration calorimetry studies, we determined that BLVRB binds its coenzyme 500-fold more tightly than its substrate. While the active site of apo BLVRB is highly dynamic on multiple timescales, active site dynamics are largely quenched within holo BLVRB, in which dynamics are redistributed to other regions of the enzyme. We show that a single point mutation of Arg78➔Ala leads to both an increase in active site micro-millisecond motions and an increase in the microscopic rate constants of coenzyme binding. This demonstrates that altering BLVRB active site dynamics can directly cause a change in functional characteristics. Our studies thus address the solution behavior of apo and holo BLVRB and identify a role of enzyme dynamics in coenzyme binding.

Relationship of femtosecond-picosecond dynamics to enzyme-catalyzed H-transfer

At physiological temperatures, enzymes exhibit a broad spectrum of conformations, which interchange via thermally activated dynamics. These conformations are sampled differently in different complexes of the protein and its ligands, and the dynamics of exchange between these conformers depends on the mass of the group that is moving and the length scale of the motion, as well as restrictions imposed by the globular fold of the enzymatic complex. Many of these motions have been examined and their role in the enzyme function illuminated, yet most experimental tools applied so far have identified dynamics at time scales of seconds to nanoseconds, which are much slower than the time scale for H-transfer between two heavy atoms. This chemical conversion and other processes involving cleavage of covalent bonds occur on picosecond to femtosecond time scales, where slower processes mask both the kinetics and dynamics. Here we present a combination of kinetic and spectroscopic methods that may enable closer examination of the relationship between enzymatic C-H → C transfer and the dynamics of the active site environment at the chemically relevant time scale. These methods include kinetic isotope effects and their temperature dependence, which are used to study the kinetic nature of the H-transfer, and 2D IR spectroscopy, which is used to study the dynamics of transition-state- and ground-state-analog complexes. The combination of these tools is likely to provide a new approach to examine the protein dynamics that directly influence the chemical conversion catalyzed by enzymes.

Structure, dynamics, and catalytic function of dihydrofolate reductase

Molecular motions are widely regarded as contributing factors in many aspects of protein function. The enzyme dihydrofolate reductase (DHFR), and particularly that from Escherichia coli, has become an important system for investigating the linkage between protein dynamics and catalytic function, both because of the location and timescales of the motions observed and because of the availability of a large amount of structural and mechanistic data that provides a detailed context within which the motions can be interpreted. Changes in protein dynamics in response to ligand binding, conformational change, and mutagenesis have been probed using numerous experimental and theoretical approaches, including X-ray crystallography, fluorescence, nuclear magnetic resonance (NMR), molecular dynamics simulations, and hybrid quantum/classical dynamics methods. These studies provide a detailed map of changes in conformation and dynamics throughout the catalytic cycle of DHFR and give new insights into the role of protein motions in the catalytic activity of this enzyme.

Conformational dynamics of the active site loop of S-adenosylmethionine synthetase illuminated by site-directed spin labeling

S-adenosylmethionine synthetase (ATP: L-methionine S-adenosyltransferase, methionine adenosyltransferase, a.k.a. MAT) is one of numerous enzymes that have a flexible polypeptide loop that moves to gate access to the active site in a motion that is closely coupled to catalysis. Crystallographic studies of this tetrameric enzyme have shown that the loop is closed in the absence of bound substrates. However, the loop must open to allow substrate binding and a variety of data indicate that the loop is closed during the catalytic steps. Previous kinetic studies indicate that during turnover loop motion occurs on a time scale of 10(-2)s, ca. 10-fold faster than chemical transformations and turnover. Site-directed spin labeling has been used to introduce nitroxide groups at two positions in the loop to illuminate how the motion of the loop is affected by substrate binding. The two loop mutants constructed, G105C and D107C, retain wild type levels of MAT activity; attachment of a methanethiosulfonate spin label to convert the cysteine to the "R1" residue reduced the k(cat) only for the labeled D107R1 form (7-fold). The K(m) value for methionine increased 2- to 4-fold for the cysteine mutants and 2- to 7-fold for the labeled proteins, whereas the K(m) for ATP was changed by at most 2-fold. EPR spectra for both labeled proteins are nearly identical and show the presence of two major spin label environments with rotational diffusion rates differing by approximately 10-fold; the slower rate is ca. 4-fold faster than the estimated protein rotational rate. The spectra are not altered by addition of substrates or products. At both positions the less mobile conformation constitutes ca. 65% of the total species, indicating an equilibrium that only slightly favors one form, that in which the label is more immobilized. The equilibrium constant that relates the two forms is comparable to the equilibrium constant of 1.5 for a conformational change that was previously deduced from the viscosity dependence of the rate of AdoMet formation. The results suggest that the motion of the loop may be an intrinsic property of the protein and not be strictly ligand modulated.

Flexibility analysis of enzyme active sites by crystallographic temperature factors

Protein flexibility is inherent to protein structural behavior. Experimental evidence for protein flexibility is extensive both in solution and in the solid state. A major question is whether the flexibility observed in enzymes is simply an inherent property of proteins that must always be borne in mind or is essential for catalysis or substrate binding. The temperature factors or B-values, as determined crystallographically, are linearly related to the mean square displacement of an atom and give an indication of atomic flexibility in the crystalline state. In this paper, we describe the frequency distributions of the normalized B-factor (B'-factor) for the active site and non-active site residues in the selected 69 apo-enzymes. This analysis was performed over the entire sequences and for different structural subsets defined by the three-dimensional structure of proteins, as alpha-helices, beta-structures and coil conformation and buried and non-buried residues. The results show that in all cases, the active site residues predominantly occur in region of low B'-factor and the non-active site residues have a tendency to exist in the high B'-factor region. This observation suggests that the active site residues, in general, are less flexible than the non-active site residues and therefore the vibrational and the fast collective motions of the C(alpha) atoms of proteins appear not to have clear biological significance.

One site fits both: a model for the ternary complex of folate + NADPH in R67 dihydrofolate reductase, a D2 symmetric enzyme

R67 dihydrofolate reductase (DHFR) is a novel enzyme that confers resistance to the antibiotic trimethoprim. The crystal structure of R67 DHFR displays a toroidal structure with a central active-site pore. This homotetrameric protein exhibits 222 symmetry, with only a few residues from each chain contributing to the active site, so related sites must be used to bind both substrate (dihydrofolate) and cofactor (NADPH) in the productive R67 DHFR.NADPH.dihydrofolate complex. Whereas the site of folate binding has been partially resolved crystallographically, an interesting question remains: how can the highly symmetrical active site also bind and orient NADPH for catalysis? To model this ternary complex, we employed DOCK and SLIDE, two methods for docking flexible ligands into proteins using quite different algorithms. The bound pteridine ring of folate (Fol I) from the crystal structure of R67 DHFR was used as the basis for docking the nicotinamide-ribose-Pi (NMN) moiety of NADPH. NMN was positioned by both DOCK and SLIDE on the opposite side of the pore from Fol I, where it interacts with Fol I at the pore's center. Numerous residues serve dual roles in binding. For example, Gln 67 from both the B and D subunits has several contacts with the pteridine ring, while the same residue from the A and C subunits has several contacts with the nicotinamide ring. The residues involved in dual roles are generally amphipathic, allowing them to make both hydrophobic and hydrophilic contacts with the ligands. The result is a 'hot spot' binding surface allowing the same residues to co-optimize the binding of two ligands, and orient them for catalysis.

X-Ray crystal structures of Candida albicans dihydrofolate reductase: high resolution ternary complexes in which the dihydronicotinamide moiety of NADPH is displaced by an inhibitor

X-ray crystallographic analysis of 5-(4'-substituted phenyl)sulfanyl-2,4-diaminoquinazoline inhibitors in ternary complex with Candida albicans dihydrofolate reductase (DHFR) and NADPH revealed two distinct modes of binding. The two compounds with small 4'-substituents (H and CH3) were found to bind with the phenyl group oriented in the plane of the quinazoline ring system and positioned adjacent to the C-helix. In contrast, the more selective inhibitors with larger 4'-substituents (tert-butyl and N-morpholino) were bound to the enzyme with the phenyl group perpendicular to the quinazoline ring and positioned in the region of the active site that typically binds the dihydronicotinamide moiety of NADPH. The cofactor appeared bound to DHFR but with the disordered dihydronicotinamide swung away from the protein surface and into solution. This unusual inhibitor binding mode may play an important role in the high DHFR selectivity of these compounds and also may provide new ideas for inhibitor design.

Mechanisms for auto-inhibition and forced product release in glycine N-methyltransferase: crystal structures of wild-type, mutant R175K and S-adenosylhomocysteine-bound R175K enzymes

Glycine N-methyltransferase (S-adenosyl-l-methionine: glycine methyltransferase, EC 2.1.1.20; GNMT) catalyzes the AdoMet-dependent methylation of glycine to form sarcosine (N-methylglycine). Unlike most methyltransferases, GNMT is a tetrameric protein showing a positive cooperativity in AdoMet binding and weak inhibition by S-adenosylhomocysteine (AdoHcy). The first crystal structure of GNMT complexed with AdoMet showed a unique "closed" molecular basket structure, in which the N-terminal section penetrates and corks the entrance of the adjacent subunit. Thus, the apparent entrance or exit of the active site is not recognizable in the subunit structure, suggesting that the enzyme must possess a second, enzymatically active, "open" structural conformation. A new crystalline form of the R175K enzyme has been grown in the presence of an excess of AdoHcy, and its crystal structure has been determined at 3.0 A resolution. In this structure, the N-terminal domain (40 amino acid residues) of each subunit has moved out of the active site of the adjacent subunit, and the entrances of the active sites are now opened widely. An AdoHcy molecule has entered the site occupied in the "closed" structure by Glu15 and Gly16 of the N-terminal domain of the adjacent subunit. An AdoHcy binds to the consensus AdoMet binding site observed in the other methyltransferase. This AdoHcy binding site supports the glycine binding site (Arg175) deduced from a chemical modification study and site-directed mutagenesis (R175K). The crystal structures of WT and R175K enzymes were also determined at 2.5 A resolution. These enzyme structures have a closed molecular basket structure and are isomorphous to the previously determined AdoMet-GNMT structure. By comparing the open structure to the closed structure, mechanisms for auto-inhibition and for the forced release of the product AdoHcy have been revealed in the GNMT structure. The N-terminal section of the adjacent subunit occupies the AdoMet binding site and thus inhibits the methyltransfer reaction, whereas the same N-terminal section forces the departure of the potentially potent inhibitor AdoHcy from the active site and thus facilitates the methyltransfer reaction. Consequently GNMT is less active at a low level of AdoMet concentration, and is only weakly inhibited by AdoHcy. These properties of GNMT are particularly suited for regulation of the cellular AdoMet/AdoHcy ratio.

Optical spectroscopy of nicotinoprotein alcohol dehydrogenase from Amycolatopsis methanolica: a comparison with horse liver alcohol dehydrogenase and UDP-galactose epimerase

The NADH absorbance spectrum of nicotinoprotein (NADH-containing) alcohol dehydrogenase from Amycolatopsis methanolica has a maximum at 326 nm. Reduced enzyme-bound pyridine dinucleotide could be reversibly oxidized by acetaldehyde. The fluorescence excitation spectrum for NADH bound to the enzyme has a maximum at 325 nm. Upon excitation at 290 nm, energy transfer from tryptophan to enzyme-bound NADH was negligible. The fluorescence emission spectrum (excitation at 325 nm) for NADH bound to the enzyme has a maximum at 422 nm. The fluorescence intensity is enhanced by a factor of 3 upon binding of isobutyramide (Kd = 59 microM). Isobutyramide acts as competitive inhibitor (Ki = 46 microM) with respect to the electron acceptor NDMA (N,N-dimethyl-p-nitrosoaniline), which binds to the enzyme containing the reduced cofactor. The nonreactive substrate analogue trifluoroethanol acts as a competitive inhibitor with respect to the substrate ethanol (Ki = 1.6 microM), which binds to the enzyme containing the oxidized cofactor. Far-UV circular dichroism spectra of the enzyme containing NADH and the enzyme containing NAD+ were identical, indicating that no major conformational changes occur upon oxidation or reduction of the cofactor. Near-UV circular dichroism spectra of NADH bound to the enzyme have a minimum at 323 nm (Deltaepsilon = -8.6 M-1 cm-1). The fluorescence anisotropy decay of enzyme-bound NADH showed no rotational freedom of the NADH cofactor. This implies a rigid environment as well as lack of motion of the fluorophore. The average fluorescence lifetime of NADH bound to the enzyme is 0.29 ns at 20 degreesC and could be resolved into at least three components (in the range 0.13-0.96 ns). Upon binding of isobutyramide to the enzyme-containing NADH, the average excited-state lifetime increased to 1.02 ns and could be resolved into two components (0.37 and 1.11 ns). The optical spectra of NADH bound to nicotinoprotein alcohol dehydrogenase have blue-shifted maxima compared to other NADH-dehydrogenase complexes, but comparable to that observed for NADH bound to horse liver alcohol dehydrogenase. The fluorescence lifetime of NADH bound to the nicotinoprotein is very short compared to enzyme-bound NADH complexes, also compared to NADH bound to horse liver alcohol dehydrogenase. The cofactor-protein interaction in the nicotinoprotein alcohol dehydrogenase active site is more rigid and apolar than that in horse liver alcohol dehydrogenase. The optical properties of NADH bound to nicotinoprotein alcohol dehydrogenase differ considerably from NADH (tightly) bound to UDP-galactose epimerase from Escherichia coli. This indicates that although both enzymes have NAD(H) as nonexchangeable cofactor, the NADH binding sites are quite different.

Effects of fluorine substitution on the structure and dynamics of complexes of dihydrofolate reductase (Escherichia coli)

Fluorine NMR experiments with a protein containing fluorinated amino acid analogs can often be used to probe structure and dynamics of the protein as well as conformational changes produced by binding of small molecules. The relevance of NMR experiments with fluorine-containing materials to characteristics of the corresponding native (nonfluorinated) proteins depends upon the extent to which these characteristics are altered by the presence of fluorine. The present work uses molecular dynamics simulations to explore the effects of replacement of tryptophan by 6-fluorotryptophan in folate and methotrexate complexes of the enzyme dihydrofolate reductase (DHFR) (Escherichia coli). Simulations of the folate-native enzyme complex produce local correlation times and order parameters that are generally in good agreement with experimental values. Simulations of the corresponding fluorotryptophan-containing system indicate that the structure and dynamics of this complex are scarcely changed by the presence of fluorinated amino acids. Calculations with the pharmacologically important methotrexate-enzyme complex predict dynamical behavior of the protein similar to that of the folate complex for both the fluorinated and native enzyme. It thus appears that, on the time scale sampled by these computer simulations, substitution of 6-fluorotryptophan for tryptophan has little effect on either the structures or dynamics of DHFR in these complexes.

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.

Dynamics of the dihydrofolate reductase-folate complex: catalytic sites and regions known to undergo conformational change exhibit diverse dynamical features

Backbone and tryptophan side-chain dynamics of uniformly 15N-labeled Escherichia coli dihydrofolate reductase were determined for the binary folate complex. The 15N T1 and T2 relaxation times and [1H]-15N heteronuclear NOEs were measured for 118 protonated backbone nitrogen atoms. The generalized order parameter (S2), the effective correlation time for internal motions (tau e), and the contribution to spin-spin relaxation through 15N exchange broadening (Rex) were determined for each residue by model-free analysis. Back-calculation of the relaxation rates for each resonance showed that the calculated dynamical parameters accurately predict the experimental data. Diverse dynamical features were evident in the DHFR backbone. Six sites exhibited order parameters significantly below the weighted mean S2 value (for the complex) of 0.81 +/- 0.002: residues G67 and D69 of the adenosine binding domain, and "hinge" residues K38 and V88, exhibited low S2 (0.29 < or = S2 < or = 0.6) and high tau e values (700 ps < or = tau e < or = 2 ns), as did sites within the beta A-alpha B loop and the beta F-beta G loop. Thus, large amplitude backbone motions, on the picosecond and nanosecond time scales, occurred at regions implicated in transition-state stabilization and in ligand-dependent conformational change. Significant Rex values (> or = 1 s-1) were determined for 45% of assigned resonances, many of which arise from residues surrounding the folate binding site. The mean S2 values of the occupied folate binding site and the unoccupied NADPH binding site were similar, indicating the backbone of the latter is at least as conformationally restricted as that of the occupied folate site. We conclude that the observed time-dependent structural fluctuations of the binary complex are in fact associated with catalytic properties of the molecule.

Protein fluorescence, dynamics and function: exploration of analogy between electronically excited and biocatalytic transition states

With the advent and development of time-resolved spectroscopic techniques and substantial progress in understanding of photophysical and photochemical phenomena, a new goal may be achieved: modeling of biochemical reaction or its elementary step by a photochemical event occurring within the probe, bound to a protein molecule. The probe may be located in a well-determined site of the protein matrix and report on the modulation of the reaction rate by the matrix and by the surrounding solvent, or by interactions in multiprotein complexes and in biomembranes. The advantages of this approach are obvious: in contrast to ordinary biochemical reaction, the excited-state reaction may be started by a short light pulse, and its kinetics may be observed directly with high resolution in time. In addition, if the reaction rate is influenced by the dynamics of the protein matrix, these dynamics may be studied simultaneously with the reaction, by using the same or a similar probe and within the same time range. In this review, the prospects for application of probes exhibiting electron transfer, proton transfer, molecular rotations and isomerizations are presented and discussed. The general problem of photochemical modeling of biochemical reactions is discussed. This modeling may result in deeper understanding of enzyme catalyzed reaction mechanisms.

Accuracy of protein flexibility predictions

Protein structural flexibility is important for catalysis, binding, and allostery. Flexibility has been predicted from amino acid sequence with a sliding window averaging technique and applied primarily to epitope search. New prediction parameters were derived from 92 refined protein structures in an unbiased selection of the Protein Data Bank by developing further the method of Karplus and Schulz (Naturwissenschaften 72:212-213, 1985). The accuracy of four flexibility prediction techniques was studied by comparing atomic temperature factors of known three-dimensional protein structures to predictions by using correlation coefficients. The size of the prediction window was optimized for each method. Predictions made with our new parameters, using an optimized window size of 9 residues in the prediction window, were giving the best results. The difference from another previously used technique was small, whereas two other methods were much poorer. Applicability of the predictions was also tested by searching for known epitopes from amino acid sequences. The best techniques predicted correctly 20 of 31 continuous epitopes in seven proteins. Flexibility parameters have previously been used for calculating protein average flexibility indices which are inversely correlated to protein stability. Indices with the new parameters showed better correlation to protein stability than those used previously; furthermore they had relationship even when the old parameters failed.

Dynamics of a flexible loop in dihydrofolate reductase from Escherichia coli and its implication for catalysis

Apo-dihydrofolate reductase from Escherichia coli samples two distinct environments slowly on the NMR time scale at room temperature. Several assigned resonances belong to residues in, or proximal to, a loop (loop I) which is comprised of residues 9-24. This exchange process was altered (either removed or made fast on the NMR time scale) by deleting three hairpin turn forming residues from the loop and filling the gap with a single glycine [Li, L., Falzone, C. J., Wright, P. E., & Benkovic, S. J. (1992) Biochemistry 31, 7826-7833]. An approximate value of 35 s-1 for the exchange rate associated with loop I in apo-DHFR was obtained in two-dimensional nuclear Overhauser spectra by analyzing the time dependence of the cross-peak volume for N epsilon H of Trp-22, a residue which is located in this loop and which has resolved cross-peaks. Owing to the critical role that this loop plays in catalysis, the correspondence between this rate of conformational exchange and off-rates for tetrahydrofolate and the reduced nicotinamide cofactor from product and substrate complexes suggests that loop movement may be a limiting factor in substrate turnover.

Methotrexate resistance in an in vivo mouse tumor due to a non-active-site dihydrofolate reductase mutation

A series of methotrexate (MTX)-resistant L1210 leukemia murine ascites tumors were developed in vivo and analyzed for drug resistance. Three of 20 tumors studied expressed an altered dihydrofolate reductase (DHFR) and each was identical, having a C to T base transition at nucleotide 46 in the DHFR gene as demonstrated by PCR and direct sequencing. This transition results in a Gly to Trp substitution at amino acid 15 of the enzyme. Purified altered enzyme displays significantly lower binding affinity for the antifolates MTX, trimetrexate, edatrexate, and trimethoprim with respective Ki values 165-, 76-, 30-, and 28-fold higher than values obtained for enzyme isolated from parental tumor (wild-type enzyme). Substrate (dihydrofolate) and cofactor (NADPH) binding is also diminished for the mutant enzyme, although to a lesser extent (17.3- and 3.6-fold higher Km, respectively). Gly-15 is highly conserved for all vertebrate species of DHFR but has no known interaction(s), either directly or indirectly, with bound cofactor, substrate, or inhibitor. Protein molecular modeling reveals that the affected residue is 9-12 A away from the enzyme active site and located in a region analogous to the mobile Met-20 loop domain characterized for Escherichia coli DHFR.

Flavin dynamics in oxidized Clostridium beijerinckii flavodoxin as assessed by time-resolved polarized fluorescence

The time-resolved fluorescence characteristics of flavin in oxidized flavodoxin isolated from the anaerobic bacterium Clostridium beijerinckii have been examined. The fluorescence intensity decays were analyzed using the maximum-entropy method. It is demonstrated that there exist large differences in fluorescence behaviour between free and protein-bound FMN. Three fluorescence lifetime components are found in oxidized flavodoxin, two of which are not present in the fluorescence-intensity decay of free FMN. The main component is distributed at 30 ps, with relative contribution of 90%. Another minor component (4% contribution) is distributed at 0.5 ns. The third component is distributed at 4.8 ns (6%), coinciding with the main distribution present in the fluorescence decay of free FMN. The results allowed us to determine the dissociation constant, Kd = 2.61 x 10(-10) M (at 20 degrees C). Collisional fluorescence-quenching experiments revealed that the flavin moiety responsible for the longest fluorescence lifetime is, at least partially, exposed to the solvent. The shortest lifetime is not affected significantly, indicating that it possibly originates from an active-site conformation in which the flavin is more or less buried in the protein and not accessible to iodide. The fluorescence anisotropy behaviour of free and protein-bound FMN was examined and analyzed with the maximum-entropy method. It was found that an excess of apoflavodoxin is required to detect differences between free and protein-bound FMN. In free FMN one single distribution of rotational correlation times is detected, whereas in flavodoxin the anisotropy decay is composed of more than one distribution. Associative analysis of fluorescence anisotropy decays shows that part of the 4.8 ns fluorescence lifetime present in the flavodoxin fluorescence decay, is coupled to a rotational correlation time similar to that of free FMN in solution, while another part of this lifetime is coupled to a longer correlation time of about 1 ns. This finding is in accordance with earlier studies [Barman, B. G. & Tollin, G. (1972) Biochemistry 11, 4746-4754] in which it was proposed that the first binding step of the flavin to the protein involves the phosphate group rather than another part of the FMN. The two shortest fluorescence lifetimes, which do not carry information on the long-term rotational behaviour of the protein, seem nonetheless to be associated with a longer rotational correlation time which is comparable to overall protein tumbling. These lifetime components probably originate from a complex in which the flavin-ring system is more or less immobilized within the protein matrix.

Characterization of the backbone dynamics of an anti-digoxin antibody VL domain by inverse detected 1H-15N NMR: comparisons with X-ray data for the Fab

The dynamic behavior of the polypeptide backbone of a recombinant antidigoxin antibody VL domain has been characterized by measurements of 15NT1 and T2 relaxation times, 1H-15N NOE values, and 1H-2H exchange rates. These data were acquired with 2D inverse detected heteronuclear 1H-15N NMR methods. The relaxation data are interpreted in terms of model free spectral density functions and exchange contributions to transverse relaxation rates R2 (= 1/T2). All characterized residues display low-amplitude picosecond time-scale librational motions. Fifteen residues undergo conformational changes on the nanosecond timescale, and 24 residues have significant R2 exchange contributions, which reflect motions on the microsecond to millisecond time-scale. For several residues, microsecond to millisecond motions of nearby aromatic rings are postulated to account for some or all of their observed R2 exchange contributions. The measured 1H-2H exchange rates are correlated with hydrogen bonding patterns and distances from the solvent accessible surface. The degree of local flexibility indicated by the NMR measurements is compared to crystallographic B-factors derived from X-ray analyses of the native Fab and the Fab/digoxin complex. In general, both the NMR and X-ray data indicate enhanced flexibility in the turns, hypervariable loops, and portions of beta-strands A, B, and G. However, on a residue-specific level, correlations among the various NMR data, and between the NMR and X-ray data, are often absent. This is attributed to the different dynamic processes and environments that influence the various observables. The combined data indicate that certain regions of the VL domain, including the three hypervariable loops, undergo dynamic changes upon VL:VH association and/or complexation with digoxin. Overall, the 26-10 VL domain exhibits relatively low flexibility on the ps-ns timescale. The possible functional consequences of this result are considered.

Flavin dynamics in reduced flavodoxins. A time-resolved polarized fluorescence study

The time-resolved fluorescence and fluorescence anisotropy characteristics of reduced flavin mononucleotide in solution as well as bound in flavodoxins isolated from the bacteria Desulfovibrio gigas, Desulfovibrio vulgaris, Clostridium beijerinckii MP and Megasphaera elsdenii have been examined. All fluorescence and fluorescence anisotropy decays were analyzed by two different methods: (a) least-squares fitting with a sum of exponentials and (b) the maximum entropy method to yield distributed lifetimes and correlation times. The results of both approaches are in excellent agreement. The fluorescence decay of the free as well as protein-bound reduced flavin chromophore is made up of three components. The shortest component proves to be relatively sensitive to the environment and can therefore be used as a diagnostic tool to probe the microenvironment of the reduced isoalloxazine ring system. The other two longer fluorescence lifetime components are insensitive to the chromophore environment and seem therefore to be related to intrinsic, photophysical properties of the reduced chromophore. Fluorescence anisotropy decays show that the flavin mononucleotide in all four reduced flavodoxins is immobilized within the protein matrix, as indicated by the recovery of a single rotational correlation time, reflecting the rotation of the whole protein. No indications are found that rapid structural fluctuations occur in reduced flavodoxins, and the mechanism of electron transfer from flavodoxin to other redox proteins seems to involve immobilized reduced flavin.