Backbone dynamics of oxidized and reduced D. vulgaris flavodoxin in solution

Structure, dynamics, and redox reactivity of an all-purpose flavodoxin

The flavodoxin of Rhodopseudomonas palustris CGA009 (Rp9Fld) supplies highly reducing equivalents to crucial enzymes such as hydrogenase, especially when the organism is iron-restricted. By acquiring those electrons from photodriven electron flow via the bifurcating electron transfer flavoprotein, Rp9Fld provides solar power to vital metabolic processes. To understand Rp9Fld's ability to work with diverse partners, we solved its crystal structure. We observed the canonical flavodoxin (Fld) fold and features common to other long-chain Flds but not all the surface loops thought to recognize partner proteins. Moreover, some of the loops display alternative structures and dynamics. To advance studies of protein-protein associations and conformational consequences, we assigned the 19F NMR signals of all five tyrosines (Tyrs). Our electrochemical measurements show that incorporation of 3-19F-Tyr in place of Tyr has only a modest effect on Rp9Fld's redox properties even though Tyrs flank the flavin on both sides. Meanwhile, the 19F probes demonstrate the expected paramagnetic effect, with signals from nearby Tyrs becoming broadened beyond detection when the flavin semiquinone is formed. However, the temperature dependencies of chemical shifts and linewidths reveal dynamics affecting loops close to the flavin and regions that bind to partners in a variety of systems. These coincide with patterns of amino acid type conservation but not retention of specific residues, arguing against detailed specificity with respect to partners. We propose that the loops surrounding the flavin adopt altered conformations upon binding to partners and may even participate actively in electron transfer.

Assignments of 19F NMR resonances and exploration of dynamics in a long-chain flavodoxin

Flavodoxin is a small protein that employs a non-covalently bound flavin to mediate single-electron transfer at low potentials. The long-chain flavodoxins possess a long surface loop that is proposed to interact with partner proteins. We have incorporated ¹⁹F-labeled tyrosine in long-chain flavodoxin from Rhodopseudomonas palustris to gain a probe of possible loop dynamics, exploiting the presence of a Tyr in the long loop in addition to Tyr residues near the flavin. We report ¹⁹F resonance assignments for all four Tyrs, and demonstration of a pair of resonances in slow exchange, both corresponding to a Tyr adjacent to the flavin. We also provide evidence for dynamics affecting the Tyr in the long loop. Thus, we show that ¹⁹F NMR of ¹⁹F-Tyr labeled flavodoxin holds promise for monitoring possible changes in conformation upon binding to partner proteins.

Conformational dynamics of Escherichia coli flavodoxins in apo- and holo-states by solution NMR spectroscopy

Flavodoxins are a family of small FMN-binding proteins that commonly exist in prokaryotes. They utilize a non-covalently bound FMN molecule to act as the redox center during the electron transfer processes in various important biological pathways. Although extensive investigations were performed, detailed molecular mechanisms of cofactor binding and electron transfer remain elusive. Herein we report the solution NMR studies on Escherichia coli flavodoxins FldA and YqcA, belonging to the long-chain and short-chain flavodoxin subfamilies respectively. Our structural studies demonstrate that both proteins show the typical flavodoxin fold, with extensive conformational exchanges observed near the FMN binding pocket in their apo-forms. Cofactor binding significantly stabilizes both proteins as revealed by the extension of secondary structures in the holo-forms, and the overall rigidity shown by the backbone dynamics data. However, the 50 s loops of both proteins in the holo-form still show conformational exchanges on the µs-ms timescales, which appears to be a common feature in the flavodoxin family, and might play an important role in structural fine-tuning during the electron transfer reactions.

NMR spectroscopy on flavins and flavoproteins

(1)H-, (11)B-, (13)C-, (15)N-, (17)O-, (19)F-, and (31)P-NMR chemical shifts of flavocoenzymes and derivatives of it, as well as of alloxazines and isoalloxazinium salts, from NMR experiments performed under various experimental conditions (e.g., dependence of the chemical shifts on temperature, concentration, solvent polarity, and pH) are reported. Also solid-state (13)C- and (15)N-NMR experiments are described revealing the anisotropic values of corresponding chemical shifts. These data, in combination with a number of coupling constants, led to a detailed description of the electronic structure of oxidized and reduced flavins. The data also demonstrate that the structure of oxidized flavin can assume a configuration deviating from coplanarity, depending on substitutions in the isoalloxazine ring, while that of reduced flavin exhibits several configurations, from almost planar to quite bended. The complexes formed between oxidized flavin and metal ions or organic molecules revealed three coordination sites with metal ions (depending on the chemical nature of the ion), and specific interactions between the pyrimidine moiety of flavin and organic molecules, mimicking specific interactions between apoflavoproteins and their coenzymes. Most NMR studies on flavoproteins were performed using (13)C- and (15)N-substituted coenzymes, either specifically enriched in the pterin moiety of flavin or uniformly labeled flavins. The chemical shifts of free flavins are used as a guide in the interpretation of the chemical shifts observed in flavoproteins. Although the hydrogen-bonding pattern in oxidized and reduced flavoproteins varies considerably, no correlation is obvious between these patterns and the corresponding redox potentials. In all reduced flavoproteins the N(1)H group of the flavocoenzyme is deprotonated, an exception is thioredoxin reductase. Three-dimensional structures of only a few flavoproteins, mostly belonging to the family of flavodoxins, have been solved. Also the kinetics of unfolding and refolding of flavodoxins has been investigated by NMR techniques. In addition, (31)P-NMR data of all so far studied flavoproteins and some (19)F-NMR spectra are discussed.

Molecular dynamics simulation of energy migration between tryptophan residues in apoflavodoxin

By performing molecular dynamics (MD) simulations of apoflavodoxin over the same timescale as fluorescence anisotropy decay measurements, the anisotropy model of two unidirectional FRET steps from two tryptophan residues to a third one can be reproduced from the tryptophan atomic coordinates in the MD trajectory.

Distant residues mediate picomolar binding affinity of a protein cofactor

Numerous proteins require cofactors to be active. Computer simulations suggest that cooperative interaction networks achieve optimal cofactor binding. There is a need for the experimental identification of the residues crucial for stabilizing these networks and thus for cofactor binding. Here we investigate the electron transporter flavodoxin, which contains flavin mononucleotide as non-covalently bound cofactor. We show that after binding flavin mononucleotide with nanomolar affinity, the protein relaxes extremely slowly (time constant ~5 days) to an energetically more favourable state with picomolar-binding affinity. Rare small-scale openings of this state are revealed through H/D exchange of N(3)H of flavin. We find that H/D exchange can pinpoint amino acids that cause tight cofactor binding. These hitherto unknown residues are dispersed throughout the structure, and many are located distantly from the flavin and seem irrelevant to flavodoxin's function. Quantification of the thermodynamics of ligand binding is important for understanding, engineering, designing and evolving ligand-binding proteins.

Flavodoxin requires tight binding of its FMN cofactor to be active, but the residues involved are unknown. In this biophysical study, FMN binding is shown to change from nanomolar to picomolar affinity on extremely slow protein relaxation and the residues responsible for cofactor binding are identified.

Mapping solvation dynamics at the function site of flavodoxin in three redox states

Flavoproteins are unique redox coenzymes, and the dynamic solvation at their function sites is critical to the understanding of their electron-transfer properties. Here, we report our complete characterization of the function-site solvation of holoflavodoxin in three redox states and of the binding-site solvation of apoflavodoxin. Using intrinsic flavin cofactor and tryptophan residue as the local optical probes with two site-specific mutations, we observed distinct ultrafast solvation dynamics at the function site in the three states and at the related recognition site of the cofactor, ranging from a few to hundreds of picoseconds. The initial ultrafast motion in 1-2.6 ps reflects the local water-network relaxation around the shallow, solvent-exposed function site. The second relaxation in 20-40 ps results from the coupled local water-protein fluctuation. The third dynamics in hundreds of picoseconds is from the intrinsic fluctuation of the loose loops flanking the cofactor at the function site. These solvation dynamics with different amplitudes well correlate with the redox states from the oxidized form, to the more rigid semiquinone and to the much looser hydroquinone. This observation of the redox control of local protein conformation plasticity and water network flexibility is significant, and such an intimate relationship is essential to the biological function of interprotein electron transfer.

Tryptophan-tryptophan energy migration as a tool to follow apoflavodoxin folding

Submolecular details of Azotobacter vinelandii apoflavodoxin (apoFD) (un)folding are revealed by time-resolved fluorescence anisotropy using wild-type protein and variants lacking one or two of apoFD's three tryptophans. ApoFD equilibrium (un)folding by guanidine hydrochloride follows a three-state model: native <--> unfolded <--> intermediate. In native protein, W128 is a sink for Förster resonance energy transfer (FRET). Consequently, unidirectional FRET with a 50-ps transfer correlation time occurs from W167 to W128. FRET from W74 to W167 is much slower (6.9 ns). In the intermediate, W128 and W167 have native-like geometry because the 50-ps transfer time is observed. However, non-native structure exists between W74 and W167 because instead of 6.9 ns the transfer correlation time is 2.0 ns. In unfolded apoFD this 2.0-ns transfer correlation time is also detected. This decrease in transfer correlation time is a result of W74 and W167 becoming solvent accessible and randomly oriented toward one another. Apparently W74 and W167 are near-natively separated in the folding intermediate and in unfolded apoFD. Both tryptophans may actually be slightly closer in space than in the native state, even though apoFD's radius increases substantially upon unfolding. In unfolded apoFD the 50-ps transfer time observed for native and intermediate folding states becomes 200 ps as W128 and W167 are marginally further separated than in the native state. Apparently, apoFD's unfolded state is not a featureless statistical coil but contains well-defined substructures. The approach presented is a powerful tool to study protein folding.

Cofactor-Apoprotein Hydrogen Bonding in Oxidized and Fully Reduced Flavodoxin Monitored by Trans-Hydrogen-Bond Scalar Couplings

Hydrogen bonding plays a key role in the tight binding of the FMN cofactor and the regulation of its redox properties in flavodoxins. Hydrogen bonding interactions can be directly observed in solution by multidimensional heteronuclear NMR spectroscopy through the scalar couplings between donor and acceptor nuclei. Here we report on the detection of intermolecular trans-hydrogen-bond couplings ((h)J) between the flavin ring system and the backbone of Desulfovibrio vulgaris flavodoxin in the oxidized and the two-electron reduced states. For this purpose, experiments are adapted from pulse sequences previously applied to determining (h)J coupling constants in nucleic acid-base pairs and proteins. The resulting (h2)J(N,N), (h4)J(N,N), (h3)J(C,N), and (h1)J(H,N) couplings involve the (15)N(1), (13)C(2), and (15)N(3) nuclei of the pyrimidine moiety of FMN, whereas no such interactions are detectable for (13)C(4) and (15)N(5). Several long-range (15)N-(15)N, (13)C-(15)N, and (1)H-(15)N J-coupling constants within the flavin are obtained as "by-products". The magnitudes of both (h)J and regular J couplings are found to be dependent on the redox state. In general, good correlations between (h)J coupling constants and donor-group (1)H chemical shifts and also crystallographic donor-acceptor distances are observed.

Improved measurement of (3)J(H(alpha)(i),N(i+1)) coupling constants in H(2)O dissolved proteins

A modification to the recently proposed alpha/beta-HN(CO)CA-J TROSY pulse sequence (P. Permi et al., J. Magn. Reson. 146, 255-259 (2000)) makes it possible to determine (3)J(H(alpha)(i), N(i+1)) coupling constants from a single E.COSY-type cross-peak pattern rather than from two (1)H(alpha) spin-state-edited subspectra. Advantages are increased (15)N resolution, critical to extracting accurate (1)H(alpha)-(15)N coupling constants, and minimized differential relaxation due to nested (13)C(alpha) and (15)N evolution periods. Application of the improved pulse sequence to Desulfovibrio vulgaris flavodoxin results in (3)J(H(alpha)(i), N(i+1)) values being systematically larger than those obtained with the original scheme. Parametrization of the coupling dependence on the protein backbone torsion angle psi yields the Karplus relation (3)J(H(alpha)(i), N(i+1))=-1.00 cos(2)(psi-120 degrees )+0.65 cos(psi-120 degrees )-0.15 Hz, with a residual root-mean-square difference of 0.13 Hz between measured and back-calculated coupling constants. The curve compares with data derived from ubiquitin (A. C. Wang and A. Bax, J. Am. Chem. Soc. 117, 1810-1813 (1995)), although spanning a slightly larger range of J values in flavodoxin. The orientation of the Ala39/Ser40 peptide link, forming a type-II beta-turn in flavodoxin, is twisted against X-ray-derived torsions by approximately 10 degrees in the NMR structure as evident from the analysis of straight phi- and psi-related (3)J coupling constants. The remaining deviation of some experimental values from the prediction is likely to be due to strong hydrogen bonding, substituent effects, or the additional dependence on the adjacent torsions straight phi.

Efficient measurement of (3)J(N,Cgamma) and (3)J(C',Cgamma) coupling constants of aromatic residues in (13)C, (15)N-labeled proteins

An NMR pulse sequence is proposed for the simultaneous determination of side chain chi1 torsion-angle related (3)J(N,Cgamma) and (3)J(C', Cgamma) couplings in aromatic amino acid spin systems. The method is of the quantitative J correlation type and takes advantage of attenuated (15)N and (1)H transverse relaxation by means of the TROSY principle. Unlike previously developed schemes for the measurement of either of the two coupling types, spectra contain internal reference peaks that are usually recorded in separate experiments. Therefore, the desired information is extracted from a single rather than four data sets. The new method is demonstrated with uniformly (13)C/(15)N labeled Desulfovibrio vulgaris flavodoxin, which contains 14 aromatic out of 147 total amino acid residues.

Comparison of backbone dynamics of oxidized and reduced putidaredoxin by 15N NMR relaxation measurements

The backbone dynamics of uniformly 15N-labeled reduced and oxidized putidaredoxin (Pdx) have been studied by 2D 15N NMR relaxation measurements. 15N T1 and T2 values and 1H-15N NOEs have been measured for the diamagnetic region of the protein. These data were analyzed by using a model-free dynamics formalism to determine the generalized order parameters (S2), the effective correlation time for internal motions (tau e), and the 15N exchange broadening contributions (Rex) for each residue, as well as the overall correlation time (tau(m)). Order parameters for the reduced Pdx are generally higher than for the oxidized Pdx, and there is increased mobility on the microsecond to millisecond time scale for the oxidized Pdx, in comparison with the reduced Pdx. These results clearly indicate that the oxidized protein exhibits higher mobility than the reduced one, which is in agreement with the recently published redox-dependent dynamics studied by amide proton exchange. In addition, we observed very high T1/T2 ratios for residues 33 and 34, giving rise to a large Rex contribution. Residue 34 is believed to be involved in the binding of Pdx to cytochrome P450cam (CYP101). The differences in the backbone dynamics are discussed in relation to the oxidation states of Pdx, and their impact on electron transfer. The entropy change occurring on oxidation of reduced Pdx has been calculated from the order parameters of the two forms.

pH-dependent spectroscopic changes associated with the hydroquinone of FMN in flavodoxins

Photoreduction with a 5-deazaflavin as the catalyst was used to convert flavodoxins from Desulfovibrio vulgaris, Megasphaera elsdenii, Anabaena PCC 7119, and Azotobacter vinelandii to their hydroquinone forms. The optical spectra of the fully reduced flavodoxins were found to vary with pH in the pH range of 5.0-8.5. The changes correspond to apparent pKa values of 6.5 and 5.8 for flavodoxins from D. vulgaris and M. elsdenii, respectively, values that are similar to the apparent pKa values reported earlier from the effects of pH on the redox potential for the semiquinone-hydroquinone couples of these two proteins (7 and 5.8, respectively). The changes in the spectra resemble those occurring with the free two-electron-reduced flavin for which the pKa is 6.7, but they are red-shifted compared with those of the free flavin. The optical changes occurring with flavodoxins from D. vulgaris and A. vinelandii flavodoxins are larger than those of free reduced FMN. The absorbance of the free and bound flavin increases in the region of 370-390 nm (Delta epsilon = 1-1.8 mM-1 cm-1) with increases of pH. Qualitatively similar pH-dependent changes occur when FMN in D. vulgaris flavodoxin is replaced by iso-FMN, and in the following mutants of D. vulgaris flavodoxin in which the residues mutated are close to the isoalloxazine of the bound flavin: D95A, D95E, D95A/D127A, W60A, Y98S, W60M/Y98W, S96R, and G61A. The 13C NMR spectrum of reduced D. vulgaris [2,4a-13C2]FMN flavodoxin shows two peaks. The peak due to C(4a) is unaffected by pH, but the peak due to C(2) broadens with decreasing pH; the apparent pKa for the change is 6.2. It is concluded that a decrease in pH induces a change in the electronic structure of the reduced flavin due to a change in the ionization state of the flavin, a change in the polarization of the flavin environment, a change in the hydrogen-bonding network around the flavin, and/or possibly a change in the bend along the N(5)-N(10) axis of the flavin. A change in the ionization state of the flavin is the simplest explanation, with the site of protonation differing from that of free FMNH-. The pH effect is unlikely to result from protonation of D95 or D127, the negatively charged amino acids closest to the flavin of D. vulgaris flavodoxin, because the optical changes observed with alanine mutants at these positions are similar to those occurring with the wild-type protein.

Probing the backbone dynamics of oxidized and reduced rat microsomal cytochrome b5 via 15N rotating frame NMR relaxation measurements: biological implications

Rotating frame 15N relaxation NMR experiments have been performed to study the local mobility of the oxidized and reduced forms of rat microsomal cytochrome b5, in the microsecond to millisecond time range. Measurements of rotating frame relaxation rates (R1rho) were performed as a function of the effective magnetic field amplitude by using off-resonance radio frequency irradiation. Detailed analysis of the two data sets resulted in the identification of slow motions along the backbone nitrogens for both oxidation states of the protein. The local mobility of reduced and oxidized cytochrome b5 turned out to be significantly different; 28 backbone nitrogens of the oxidized form were shown to participate in a conformational exchange process, while this number dropped to 12 in the reduced form. The correlation time, tauex, for the exchange processes could be determined for 21 and 9 backbone nitrogens for oxidized and reduced cytochrome b5, respectively, with their values ranging between 70 and 280 microseconds. The direct experimental evidence provided in this study for the larger mobility of the oxidized form of the protein is consistent with the different backbone NH solvent exchangeability recently documented for the two oxidation states [Arnesano, F., et al. (1998) Biochemistry 37, 173-184]. Our experimental observations may have significant biological implications. The differential local mobility between the two oxidation states is proposed to be an important factor controlling the molecular recognition processes in which cytochrome b5 is involved.

Modulation of the redox potentials of FMN in Desulfovibrio vulgaris flavodoxin: thermodynamic properties and crystal structures of glycine-61 mutants

Mutants of the electron-transfer protein flavodoxin from Desulfovibrio vulgaris were made by site-directed mutagenesis to investigate the role of glycine-61 in stabilizing the semiquinone of FMN by the protein and in controlling the flavin redox potentials. The spectroscopic properties, oxidation-reduction potentials, and flavin-binding properties of the mutant proteins, G61A/N/V and L, were compared with those of wild-type flavodoxin. The affinities of all of the mutant apoproteins for FMN and riboflavin were less than that of the wild-type apoprotein, and the redox potentials of the two 1-electron steps in the reduction of the complex with FMN were also affected by the mutations. Values for the dissociation constants of the complexes of the apoprotein with the semiquinone and hydroquinone forms of FMN were calculated from the redox potentials and the dissociation constant of the oxidized complex and used to derive the free energies of binding of the FMN in its three oxidation states. These showed that the semiquinone is destabilized in all of the mutants, and that the extent of destabilization tends to increase with increasing bulkiness of the side chain at residue 61. It is concluded that the hydrogen bond between the carbonyl of glycine-61 and N(5)H of FMN semiquinone in wild-type flavodoxin is either absent or severely impaired in the mutants. X-ray crystal structure analysis of the oxidized forms of the four mutant proteins shows that the protein loop that contains residue 61 is moved away from the flavin by 5-6 A. The hydrogen bond formed between the backbone nitrogen of aspartate-62 and O(4) of the dimethylisoalloxazine of the flavin in wild-type flavodoxin is absent in the mutants. Reliable structural information was not obtained for the reduced forms of the mutant proteins, but if the mutants change conformation when the flavin is reduced to the semiquinone, to facilitate hydrogen bonding between N(5)H and the carbonyl of residue 61, then the change must be different from that known to occur in wild-type flavodoxin.