DEER distance measurement between a spin label and a native FAD semiquinone in electron transfer flavoprotein

Structural insights into the semiquinone form of human Cytochrome P450 reductase by DEER distance measurements between a native flavin and a spin labelled non-canonical amino acid

The flavoprotein Cytochrome P450 reductase (CPR) is the unique electron pathway from NADPH to Cytochrome P450 (CYPs). The conformational dynamics of human CPR in solution, which involves transitions from a "locked/closed" to an "unlocked/open" state, is crucial for electron transfer. To date, however, the factors guiding these changes remain unknown. By Site-Directed Spin Labelling coupled to Electron Paramagnetic Resonance spectroscopy, we have incorporated a non-canonical amino acid onto the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) domains of soluble human CPR, and labelled it with a specific nitroxide spin probe. Taking advantage of the endogenous FMN cofactor, we successfully measured for the first time, the distance distribution by DEER between the semiquinone state FMNH• and the nitroxide. The DEER data revealed a salt concentration-dependent distance distribution, evidence of an "open" CPR conformation at high salt concentrations exceeding previous reports. We also conducted molecular dynamics simulations which unveiled a diverse ensemble of conformations for the "open" semiquinone state of the CPR at high salt concentration. This study unravels the conformational landscape of the one electron reduced state of CPR, which had never been studied before.

Molecular Dynamics Simulations Based on Newly Developed Force Field Parameters for Cu2+ Spin Labels Provide Insights into Double-Histidine-Based Double Electron-Electron Resonance

Electron paramagnetic resonance (EPR) in combination with the recently developed double-histidine (dHis)-based Cu²⁺ spin labeling has provided valuable insights into protein structure and conformational dynamics. To relate sparse distance constraints measured by EPR to protein fluctuations in solution, modeling techniques are needed. In this work, we have developed force field parameters for Cu²⁺-nitrilotriacetic and Cu²⁺-iminodiacetic acid spin labels. We employed molecular dynamics (MD) simulations to capture the atomic-level details of dHis-labeled protein fluctuations. The interspin distances extracted from 200 ns MD trajectories show good agreement with the experimental results. The MD simulations also illustrate the dramatic rigidity of the Cu²⁺ labels compared to the standard nitroxide spin label. Further, the relative orientations between spin-labeled sites were measured to provide insight into the use of double electron-electron resonance (DEER) methods for such labels. The relative mean angles, as well as the standard deviations of the relative angles, agree well in general with the spectral simulations published previously. The fluctuations of relative orientations help rationalize why orientation selectivity effects are minimal at X-band frequencies, but observable at the Q-band for such labels. In summary, the results show that by combining the experimental results with MD simulations precise information about protein conformations as well as flexibility can be obtained.

Pulsed EPR Dipolar Spectroscopy on Spin Pairs with one Highly Anisotropic Spin Center: The Low-Spin FeIII Case

Pulsed electron paramagnetic resonance (EPR) dipolar spectroscopy (PDS) offers several methods for measuring dipolar coupling constants and thus the distance between electron spin centers. Up to now, PDS measurements have been mostly applied to spin centers whose g-anisotropies are moderate and therefore have a negligible effect on the dipolar coupling constants. In contrast, spin centers with large g-anisotropy yield dipolar coupling constants that depend on the g-values. In this case, the usual methods of extracting distances from the raw PDS data cannot be applied. Here, the effect of the g-anisotropy on PDS data is studied in detail on the example of the low-spin Fe3+ ion. First, this effect is described theoretically, using the work of Bedilo and Maryasov (Appl. Magn. Reson. 2006, 30, 683-702) as a basis. Then, two known Fe3+ /nitroxide compounds and one new Fe3+ /trityl compound were synthesized and PDS measurements were carried out on them using a method called relaxation induced dipolar modulation enhancement (RIDME). Based on the theoretical results, a RIDME data analysis procedure was developed, which facilitated the extraction of the inter-spin distance and the orientation of the inter-spin vector relative to the Fe3+ g-tensor frame from the RIDME data. The accuracy of the determined distances and orientations was confirmed by comparison with MD simulations. This method can thus be applied to the highly relevant class of metalloproteins with, for example, low-spin Fe3+ ions.

EPR spectroscopy on flavin radicals in flavoproteins

Flavin semiquinone redox states are important intermediates in a broad variety of reactions catalyzed by flavoproteins. As paramagnetic states they can be favorably probed by EPR spectroscopy in all its flavors. This review summarizes recent results in the characterization of flavin radicals. On the one hand, flavin radical states, e.g., trapped as reaction intermediates, can be characterized using modern pulsed EPR methods to unravel their electronic structure and to gain information about the surrounding environment and its changes on protein action. On the other hand, short-lived intermediate flavin radical states generated, e.g., photochemically, can be followed by time-resolved EPR, which allows a direct tracking of flavin-dependent reactions with a temporal resolution reaching nanoseconds.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

EPR on Flavoproteins

Flavoproteins often employ radical mechanisms in their enzymatic reactions. This involves paramagnetic species, which can ideally be investigated with electron paramagnetic resonance (EPR) spectroscopy. In this chapter we focus on the example of flavin-based photoreceptors and discuss, how different EPR methods have been used to extract information about the flavin radical's electronic state, its binding pocket, electron-transfer pathways, and about the protein's tertiary and quaternary structure.

Elucidating a relationship between conformational sampling and drug resistance in HIV-1 protease

Enzyme targets in rapidly replicating systems, such as retroviruses, commonly respond to drug-selective pressure with mutations arising in the active site pocket that limit inhibitor effectiveness by introducing steric hindrance or by eliminating essential molecular interactions. However, these primary mutations are disposed to compromising pathogenic fitness. Emerging secondary mutations, which are often found outside of the binding cavity, may or can restore fitness while maintaining drug resistance. The accumulated drug pressure selected mutations could have an indirect effect in the development of resistance, such as altering protein flexibility or the dynamics of protein-ligand interactions. Here, we show that accumulation of mutations in a drug-resistant HIV-1 protease (HIV-1 PR) variant, D30N/M36I/A71V, changes the fractional occupancy of the equilibrium conformational sampling ensemble. Correlations are made among populations of the conformational states, namely, closed-like, semiopen, and open-like, with inhibition constants, as well as kinetic parameters. Mutations that stabilize a closed-like conformation correlate with enzymes of lowered activity and with higher affinity for inhibitors, which is corroborated by a further increase in the fractional occupancy of the closed state upon addition of inhibitor or substrate-mimic. Cross-resistance is found to correlate with combinations of mutations that increase the population of the open-like conformations at the expense of the closed-like state while retaining native-like occupancy of the semiopen population. These correlations suggest that at least three states are required in the conformational sampling model to establish the emergence of drug resistance in HIV-1 PR. More importantly, these results shed light on a possible mechanism whereby mutations combine to impart drug resistance while maintaining catalytic activity.

Inhibitor-induced conformational shifts and ligand-exchange dynamics for HIV-1 protease measured by pulsed EPR and NMR spectroscopy

Double electron-electron resonance (DEER) spectroscopy was utilized to investigate shifts in conformational sampling induced by nine FDA-approved protease inhibitors (PIs) and a nonhydrolyzable substrate mimic for human immunodeficiency virus type 1 protease (HIV-1 PR) subtype B, subtype C, and CRF_01 A/E. The ligand-bound subtype C protease has broader DEER distance profiles, but trends for inhibitor-induced conformational shifts are comparable to those previously reported for subtype B. Ritonavir, one of the strong-binding inhibitors for subtypes B and C, induces less of the closed conformation in CRF_01 A/E. (1)H-(15)N heteronuclear single-quantum coherence (HSQC) spectra were acquired for each protease construct titrated with the same set of inhibitors. NMR (1)H-(15)N HSQC titration data show that inhibitor residence time in the protein binding pocket, inferred from resonance exchange broadening, shifting or splitting correlates with the degree of ligand-induced flap closure measured by DEER spectroscopy. These parallel results show that the ligand-induced conformational shifts resulting from protein-ligand interactions characterized by DEER spectroscopy of HIV-1 PR obtained at the cryogenic temperature are consistent with more physiological solution protein-ligand interactions observed by solution NMR spectroscopy.

DEER-Stitch: combining three- and four-pulse DEER measurements for high sensitivity, deadtime free data

Over approximately the last 15 years the electron paramagnetic resonance (EPR) technique of double electron electron resonance (DEER) has attracted considerable attention since it allows for the precise measurement of the dipole-dipole coupling between radicals and thus can lead to distance information between pairs of radicals separated by up to ca. 8 nm. The "deadtime free" 4-pulse DEER sequence is widely used but can suffer from poor sensitivity if the electron spin-echo decays too quickly to allow collection of a sufficiently long time trace. In this paper we present a method which takes advantage of the much greater sensitivity that the 3-pulse sequence offers over the 4-pulse sequence since the measured electron spin-echo intensity (for equal sequence lengths) is greater. By combining 3- and 4-pulse DEER time traces using a method coined DEER-Stitch (DEERS) accurate dipole-dipole coupling measurements can be made which combine the sensitivity of the 3-pulse DEER sequence with the deadtime free advantage of the 4-pulse DEER sequence. To develop the DEER-Stitch method three systems were measured: a semi-rigid bis-nitroxide labeled nanowire, the bis-nitroxide labeled protein CD55 with a distance between labels of almost 8 nm and a dimeric copper amine oxidase from Arthrobacter globiformis (AGAO).

Sodium-dependent movement of covalently bound FMN residue(s) in Na(+)-translocating NADH:quinone oxidoreductase

Na(+)-translocating NADH:quinone oxidoreductase (Na(+)-NQR) is a component of respiratory electron-transport chain of various bacteria generating redox-driven transmembrane electrochemical Na(+) potential. We found that the change in Na(+) concentration in the reaction medium has no effect on the thermodynamic properties of prosthetic groups of Na(+)-NQR from Vibrio harveyi, as was revealed by the anaerobic equilibrium redox titration of the enzyme's EPR spectra. On the other hand, the change in Na(+) concentration strongly alters the EPR spectral properties of the radical pair formed by the two anionic semiquinones of FMN residues bound to the NqrB and NqrC subunits (FMN(NqrB) and FMN(NqrC)). Using data obtained by pulse X- and Q-band EPR as well as by pulse ENDOR and ELDOR spectroscopy, the interspin distance between FMN(NqrB) and FMN(NqrC) was found to be 15.3 Å in the absence and 20.4 Å in the presence of Na(+), respectively. Thus, the distance between the covalently bound FMN residues can vary by about 5 Å upon changes in Na(+) concentration. Using these results, we propose a scheme of the sodium potential generation by Na(+)-NQR based on the redox- and sodium-dependent conformational changes in the enzyme.

1H dynamic nuclear polarization based on an endogenous radical

We demonstrate a 15-fold enhancement of solid-state NMR signals via dynamic nuclear polarization (DNP) based on a stable, naturally occurring radical in a protein: the flavin mononucleotide (FMN) semiquinone of flavodoxin. The line width of flavodoxin's EPR signal suggests that the dominant DNP mechanism is the solid effect, consistent with the field-dependent DNP enhancement profile. The magnitude of the enhancement as well as the bulk-polarization build-up time constant (τ(B)) with which it develops are dependent on the isotopic composition of the protein. Deuteration of the protein to 85% increased the nuclear longitudinal relaxation time T(1n) and τ(B) by factors of five and seven, respectively. Slowed dissipation of polarization can explain the 2-fold higher maximal enhancement than that obtained in proteated protein, based on the endogenous semiquinone. In contrast, the long τ(B) of TOTAPOL-based DNP in nonglassy samples was not accompanied by a similarly important long T(1n), and in this case the enhancement was greatly reduced. The low concentrations of radicals occurring naturally in biological systems limit the magnitude of DNP enhancement that is attainable by this means. However, our enhancement factors of up to 15 can nonetheless make an important difference to the feasibility of applying solid-state NMR to biochemical systems. We speculate that DNP based on endogenous radicals may facilitate MAS NMR characterization of biochemical complexes and even organelles, and could also serve as a source of additional structural and physiological information.

Determination of the distance between the Mo(V) and Fe(III) heme centers of wild type human sulfite oxidase by pulsed EPR spectroscopy

Intramolecular electron transfer (IET) between the molybdenum and heme centers of vertebrate sulfite oxidase (SO) is proposed to be a key step in the catalytic cycle of the enzyme. However, the X-ray crystallographic distance between these centers, R(MoFe) = 32.3 Å, appears to be too long for the rapid IET rates observed in liquid solution. The Mo and heme domains are linked by a flexible tether, and it has been proposed that dynamic interdomain motion brings the two metal centers closer together and thereby facilitates rapid IET. To date, there have been no direct distance measurements for SO in solution that would support or contradict this model. In this work, pulsed electron-electron double resonance (ELDOR) and relaxation induced dipolar modulation enhancement (RIDME) techniques were used to obtain information about R(MoFe) in the Mo(V)Fe(III) state of wild type recombinant human SO in frozen glassy solution. Surprisingly, the data obtained suggest a fixed structure with R(MoFe) = 32 Å, similar to that determined by X-ray crystallography for chicken SO, although the orientation of the R(MoFe) radius-vector with respect to the heme center was found to be somewhat different. The implications of these findings for the flexible tether model are discussed.

Electron transfer flavoprotein domain II orientation monitored using double electron-electron resonance between an enzymatically reduced, native FAD cofactor, and spin labels

Human electron transfer flavoprotein (ETF) is a soluble mitochondrial heterodimeric flavoprotein that links fatty acid β-oxidation to the main respiratory chain. The crystal structure of human ETF bound to medium chain acyl-CoA dehydrogenase indicates that the flavin adenine dinucleotide (FAD) domain (αII) is mobile, which permits more rapid electron transfer with donors and acceptors by providing closer access to the flavin and allows ETF to accept electrons from at least 10 different flavoprotein dehydrogenases. Sequence homology is high and low-angle X-ray scattering is identical for Paracoccus denitrificans (P. denitrificans) and human ETF. To characterize the orientations of the αII domain of P. denitrificans ETF, distances between enzymatically reduced FAD and spin labels in the three structural domains were measured by double electron-electron resonance (DEER) at X- and Q-bands. An FAD to spin label distance of 2.8 ± 0.15 nm for the label in the FAD-containing αII domain (A210C) agreed with estimates from the crystal structure (3.0 nm), molecular dynamics simulations (2.7 nm), and rotamer library analysis (2.8 nm). Distances between the reduced FAD and labels in αI (A43C) were between 4.0 and 4.5 ± 0.35 nm and for βIII (A111C) the distance was 4.3 ± 0.15 nm. These values were intermediate between estimates from the crystal structure of P. denitrificans ETF and a homology model based on substrate-bound human ETF. These distances suggest that the αII domain adopts orientations in solution that are intermediate between those which are observed in the crystal structures of free ETF (closed) and ETF bound to a dehydrogenase (open).

Studying biomolecular complexes with pulsed electron–electron double resonance spectroscopy

The function of biomolecules is intrinsically linked to their structure and the complexes they form during function. Techniques for the determination of structures and dynamics of these nanometre assemblies are therefore important for an understanding on the molecular level. PELDOR (pulsed electron–electron double resonance) is a pulsed EPR method that can be used to reliably and precisely measure distances in the range 1.5–8 nm, to unravel orientations and to determine the number of monomers in complexes. In conjunction with site-directed spin labelling, it can be applied to biomolecules of all sizes in aqueous solutions or membranes. PELDOR is therefore complementary to the methods of X-ray crystallography, NMR and FRET (fluorescence resonance energy transfer) and is becoming a powerful method for structural determination of biomolecules. In the present review, the methods of PELDOR are discussed and examples where PELDOR has been used to obtain structural information on biomolecules are summarized.

Rotamer libraries of spin labelled cysteines for protein studies

Studies of structure and dynamics of proteins using site-directed spin labelling rely on explicit modelling of spin label conformations. The large computational effort associated with such modelling with molecular dynamics (MD) simulations can be avoided by a rotamer library approach based on a coarse-grained representation of the conformational space of the spin label. We show here that libraries of about 200 rotamers, obtained by iterative projection of a long MD trajectory of the free spin label onto a set of canonical dihedral angles, provide a representation of the underlying trajectory adequate for EPR distance measurements. Rotamer analysis was performed on selected X-ray structures of spin labelled T4 lysozyme mutants to characterize the spin label rotamer ensemble on a single protein site. Furthermore, predictions based on the rotamer library approach are shown to be in nearly quantitative agreement with electron paramagnetic resonance (EPR) distance data on the Na(+)/H(+) antiporter NhaA and on the light-harvesting complex LHCII whose structures are known from independent cryo electron microscopy and X-ray studies, respectively. Suggestions for the selection of labelling sites in proteins are given, limitations of the approach discussed, and requirements for further development are outlined.

Nature of the energy landscape for gated electron transfer in a dynamic redox protein

Conformational control limits most electron transfer (ET) reactions in biology, but we lack general insight into the extent of conformational space explored, and specifically the properties of the associated energy landscape. Here we unite electron-electron double resonance (ELDOR) studies of the diradical (disemiquinoid) form of human cytochrome P450 reductase (CPR), a nicotinamide adenine phosphate dinucleotide (NADPH)-linked diflavin oxidoreductase required for P450 enzyme reduction, with functional studies of internal ET to gain new insight into the extent and properties of the energy landscape for conformationally controlled ET. We have identified multiple conformations of disemiquinoid CPR, which point to a rugged energy landscape for conformational sampling consistent with functional analysis of ET using high-pressure stopped-flow, solvent, and temperature perturbation studies. Crystal structures of CPR have identified discrete "closed" and "open" states, but we emphasize the importance of a continuum of conformational states across the energy landscape. Within the landscape more closed states that favor internal ET are formed by nucleotide binding. Open states that enable P450 enzymes to gain access to electrons located in the FMN-domain are favored in the absence of bound coenzyme. The extent and nature of energy landscapes are therefore accessible through the integration of ELDOR spectroscopy with functional studies. We suggest this is a general approach that can be used to gain new insight into energy landscapes for biological ET mediated by conformational sampling mechanisms.