Conformational substates in azurin

Heterogeneity of protein substates visualized by spin-label EPR

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na,K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both α-helical and β-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to ∼ 250 K or higher. Above ∼ 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na,K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.

Ligand dynamics in heme proteins observed by Fourier transform infrared-temperature derivative spectroscopy

Fourier transform infrared (FTIR) spectroscopy is a powerful tool for the investigation of protein-ligand interactions in heme proteins. Nitric oxide and carbon monoxide are attractive physiologically relevant ligands because their bond stretching vibrations give rise to strong mid-infrared absorption bands that can be measured with exquisite sensitivity and precision using photolysis difference spectroscopy at cryogenic temperatures. These stretching bands are fine-tuned by electrostatic interactions with the environment and, therefore, ligands can be utilized as local probes of structure and dynamics. Bound to the heme iron, the ligand stretching bands are susceptible to changes in the iron-ligand bond and the electric field at the active site. Upon photolysis, the vibrational bands display changes due to ligand relocation to docking sites within the protein, rotational motions of the ligand in these sites and protein conformational changes. Photolysis difference spectra taken over a wide temperature range (3-300K) using specific temperature protocols for sample photodissociation can provide detailed insights into both protein and ligand dynamics. Moreover, temperature-derivative spectroscopy (TDS) has proven to be a particularly powerful technique to study protein-ligand interactions. The FTIR-TDS technique has been extensively applied to studies of carbon monoxide binding to heme proteins, whereas measurements with nitric oxide are still scarce. Here we describe infrared cryo-spectroscopy and present a variety of applications to the study of protein-ligand interactions in heme proteins. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Reaction of the Co(II)-substrate radical pair catalytic intermediate in coenzyme B12-dependent ethanolamine ammonia-lyase in frozen aqueous solution from 190 to 217 K

The decay kinetics of the aminoethanol-generated Co(II)-substrate radical pair catalytic intermediate in ethanolamine ammonia-lyase from Salmonella typhimurium have been measured on timescales of <10(5) s in frozen aqueous solution from 190 to 217 K. X-band continuous-wave electron paramagnetic resonance (EPR) spectroscopy of the disordered samples has been used to continuously monitor the full radical pair EPR spectrum during progress of the decay after temperature step reaction initiation. The decay to a diamagnetic state is complete and no paramagnetic intermediate states are detected. The decay exhibits three kinetic regimes in the measured temperature range, as follows. i), Low temperature range, 190 < or = T < or = 207 K: the decay is biexponential with constant fast (0.57 +/- 0.04) and slow (0.43 +/- 0.04) phase amplitudes. ii), Transition temperature range, 207 < T < 214 K: the amplitude of the slow phase decreases to zero with a compensatory rise in the fast phase amplitude, with increasing temperature. iii), High temperature range, T > or = 214 K: the decay is monoexponential. The observed first-order rate constants for the monoexponential (k(obs,m)) and the fast phase of the biexponential decay (k(obs,f)) adhere to the same linear relation on an lnk versus T(-1) (Arrhenius) plot. Thus, k(obs,m) and k(obs,f) correspond to the same apparent Arrhenius prefactor and activation energy (logA(app,f) (s(-1)) = 13.0, E(a,app,f) = 15.0 kcal/mol), and therefore, a common decay mechanism. We propose that k(obs,m) and k(obs,f) represent the native, forward reaction of the substrate through the radical rearrangement step. The slow phase rate constant (k(obs,s)) for 190 < or = T < or = 207 K obeys a different linear Arrhenius relation (logA(app,s) (s(-1)) = 13.9, E(a,app,s) = 16.6 kcal/mol). In the transition temperature range, k(obs,s) displays a super-Arrhenius increase with increasing temperature. The change in E(a,app,s) with temperature and the narrow range over which it occurs suggest an origin in a liquid/glass or dynamical transition. A discontinuity in the activation barrier for the chemical reaction is not expected in the transition temperature range. Therefore, the transition arises from a change in the properties of the protein. We propose that a protein dynamical contribution to the reaction, which is present above the transition temperature, is lost below the transition temperature, owing to an increase in the activation energy barrier for protein motions that are coupled to the reaction. For both the fast and slow phases of the low temperature decay, the dynamical transition in protein motions that are obligatorily coupled to the reaction of the Co(II)-substrate radical pair lies below 190 K.

Influence of distal residue B10 on CO dynamics in myoglobin and neuroglobin

For many years, myoglobin has served as a paradigm for structure-function studies in proteins. Ligand binding and migration within myoglobin has been studied in great detail by crystallography and spectroscopy, showing that gaseous ligands such as O(2), CO, and NO not only bind to the heme iron but may also reside transiently in three internal ligand docking sites, the primary docking site B and secondary sites C and D. These sites affect ligand association and dissociation in specific ways. Neuroglobin is another vertebrate heme protein that also binds small ligands. Ligand migration pathways in neuroglobin have not yet been elucidated. Here, we have used Fourier transform infrared temperature derivative spectroscopy at cryogenic temperatures to compare the influence of the side chain volume of amino acid residue B10 on ligand migration to and rebinding from docking sites in myoglobin and neuroglobin.

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.

Dominant features of protein reaction dynamics: conformational relaxation and ligand migration

Here, we review the dominant aspects of protein dynamics as revealed by studying hemoproteins using the combination of laser flash photolysis, kinetic spectroscopy and low temperature. The first breakthrough was the finding that geminate ligand rebinding with myoglobin is highly non-exponential at temperature T<200 K, providing evidence for the trapping of a large number of protein statistical substates. Another major advance was the introduction of a "model free" approach to analyze polychromatic kinetics in terms of their rate spectrum rather than to fit the data to some arbitrarily predefined kinetic scheme. Kinetic processes are identified and quantified directly from the rate spectrum without a priori assumptions. In recent years, further progresses were achieved by using xenon gas as a soft external perturbing agent that competes with ligand rebinding pathways by occupying hydrophobic protein cavities. The first part of this paper introduces several basic principles that are spread throughout a vast literature. The second part describes the main conclusions regarding conformational relaxation and ligand migration in hemoproteins obtained by combining these approaches.

The Solution Structure of [Cu(aq)]2+ and Its Implications for Rack-Induced Bonding in Blue Copper Protein Active Sites

The structure of [Cu(aq)]2+ has been investigated by using full multiple-scattering theoretical (MXAN) analysis of the copper K-edge X-ray absorption (XAS) spectrum and density functional theory (DFT) to test both ideal Td and square-planar four-coordinate, five-coordinate square-pyramidal, and six-coordinate octahedral [Cu(aq)]2+ models. The best fit was an elongated five-coordinate square pyramid with four Cu-O(eq) bonds (2 x 1.98 +/- 0.03 A and 2 x 1.95 +/- 0.03 A) and a long Cu-O(ax) bond (2.35 +/- 0.05 A). The four equatorial ligands were D2d-distorted from the mean equatorial plane by +/-(17 +/- 4) degrees, so that the overall symmetry of [Cu(H2O)5]2+ is C2v. The four-coordinate MXAN fit was nearly as good, but the water ligands (4 x 1.96 +/- 0.02 A) migrated +/-(13 +/- 4) degrees from the mean equatorial plane, making the [Cu(H2O)4]2+ model again D2d-distorted. Spectroscopically calibrated DFT calculations were carried out on the C2v elongate square-pyramidal and D2d-distorted four-coordinate MXAN copper models, providing comparative electronic structures of the experimentally observed geometries. These calculations showed 0.85e spin on Cu(II) and 0.03e electron spin on each of the four equatorial water oxygens. All covalent bonding was restricted to the equatorial plane. In the square-pyramidal model, the electrostatic Cu-O(ax) bond was worth only 96.8 kJ mol(-1), compared to 304.6 kJ mol(-1) for each Cu-O(eq) bond. Both MXAN and DFT showed the potential well of the axial bond to be broad and flat, allowing large low-energy excursions. The irregular geometry and D2d-distorted equatorial ligand set sustained by unconstrained [Cu(H2O)5]2+ warrants caution in drawing conclusions regarding structural preferences from small molecule crystal structures and raises questions about the site-structural basis of the rack-induced bonding hypothesis of blue copper proteins. Further, previously neglected protein folding thermodynamic consequences of the rack-bonding hypothesis indicate an experimental disconfirmation.

Disentangling ligand migration and heme pocket relaxation in cytochrome P450cam

In this work we show that ligand migration and active site conformational relaxation can occur independently of each other in hemoproteins. The complicated kinetics of carbon monoxide rebinding with cytochrome P450cam display up to five distinct processes between 77 K and 300 K. They were disentangled by using a combination of three approaches: 1), the competition of the ligand with xenon for the occupation of internal protein cavities; 2), the modulation of the amount of distal steric hindrance within the heme pocket by varying the nature of the substrate; and 3), molecular mechanics calculations to support the proposed heme-substrate relaxation mechanism and to seek internal cavities. In cytochrome P450cam, active site conformational relaxation results from the displacement of the substrate toward the heme center upon photodissociation of the ligand. It is responsible for the long, puzzling bimodal nature of the rebinding kinetics observed down to 77 K. The relaxation rate is strongly substrate-dependent. Ligand migration is slower and is observed only above 135 K. Migration and return rates are independent of the substrate.

Photodynamics of red fluorescent proteins studied by fluorescence correlation spectroscopy

Red fluorescent proteins are important tools in fluorescence-based life science research. Recently, we have introduced eqFP611, a red fluorescent protein with advantageous properties from the sea anemone Entacmaea quadricolor. Here, we have studied the submillisecond light-driven intramolecular dynamics between bright and dark states of eqFP611 and, for comparison, drFP583 (DsRed) by using fluorescence correlation spectroscopy on protein solutions. A three-state model with one dark and two fluorescent states describes the power-dependence of the flickering dynamics of both proteins at different excitation wavelengths. It involves two light-driven conformational transitions. We have also studied the photodynamics of individual (monomeric) eqFP611 molecules immobilized on surfaces. The flickering rates and dark state fractions of eqFP611 bound to polyethylene glycol-covered glass surfaces were identical to those measured in solution, showing that the bound FPs behaved identically. A second, much slower flickering process was observed on the 10-ms timescale. Deposition of eqFP611 molecules on bare glass surfaces yielded bright fluorescence without any detectable flickering and a >10-fold decreased photobleaching yield. These observations underscore the intimate connection between protein motions and photophysical processes in fluorescent proteins.

Slaving: solvent fluctuations dominate protein dynamics and functions

Protein motions are essential for function. Comparing protein processes with the dielectric fluctuations of the surrounding solvent shows that they fall into two classes: nonslaved and slaved. Nonslaved processes are independent of the solvent motions; their rates are determined by the protein conformation and vibrational dynamics. Slaved processes are tightly coupled to the solvent; their rates have approximately the same temperature dependence as the rate of the solvent fluctuations, but they are smaller. Because the temperature dependence is determined by the activation enthalpy, we propose that the solvent is responsible for the activation enthalpy, whereas the protein and the hydration shell control the activation entropy through the energy landscape. Bond formation is the prototype of nonslaved processes; opening and closing of channels are quintessential slaved motions. The prevalence of slaved motions highlights the importance of the environment in cells and membranes for the function of proteins.

The reactions of copper proteins with nitric oxide

Nitric oxide (NO) can act as a ligand for copper atoms and may also engage in redox chemistry with the metal once bound. Furthermore NO posses an unpaired electron which can couple with the unpaired electron on Cu2+. These properties have been exploited to probe the active sites of copper-containing enzymes and proteins. We review these studies. In addition to the use as a spectroscopic probe for the active site we draw attention to the rapid reactions of NO at the copper sites in Cytochrome c oxidase (CcO) and laccase. These reactions in CcO occur in the ms time range, at low NO concentrations and in the presence of oxygen and may therefore be of physiological relevance to the control of respiration. Finally we speculate on the wider role that NO may play in regulation of an important group of Type 2 copper containing enzymes.

Two-dimensional distributions of activation enthalpy and entropy from kinetics by the maximum entropy method

The maximum entropy method (MEM) is used to numerically invert the kinetics of ligand rebinding at low temperatures to obtain the underlying two-dimensional distribution of activation enthalpies and entropies, g(H,S). A global analysis of the rebinding of carbon monoxide (CO) to myoglobin (Mb), monitored in the Soret band at temperatures from 60 to 150 K, is performed using a Newton-Raphson optimization algorithm. The MEM approach describes the data much better than traditional least-squares analyses, reducing chi 2 by an order of magnitude. The MEM resolves two barrier distributions suggestive of rebinding to different bound conformations of MbCO, the so-called A1 and A3 substates, whose activation barriers have been independently estimated from kinetics monitored in the infrared. The distribution corresponding to A3 possesses higher activation entropies, also consistent with infrared measurements. Within an A substate, correlations of S and H are recovered qualitatively from simulated data but can be difficult to obtain from experimental data. When the rebinding measured at 60 K is excluded from the inversion, two peaks are no longer clearly resolved. Thus, data of very high quality are required to unambiguously determine the kinetic resolvability of subpopulations and the shape of the barrier distribution for a single A substate.

Ligand binding and protein dynamics in cupredoxins

Type 1 copper sites bind nitric oxide (NO) in a photolabile complex. We have studied the NO binding properties of the type 1 copper sites in two cupredoxins, azurin and halocyanin, by measuring the temperature dependence of the ligand binding equilibria and the kinetics of the association reaction after photodissociation over a wide range of temperature (80-280 K) and time (10(-6)-10(2) s). In both proteins, we find nonexponential kinetics below 200 K that do not depend on the NO concentration. Consequently, this process is interpreted as geminate recombination. In azurin, the rebinding can be modeled with the Arrhenius law using a single pre-exponential factor of 10(8.3) s-1 and a Gaussian distribution of enthalpy barriers centered at 22 kJ/mol with a width [full width at half-maximum (FWHM)] of 11 kJ/mol. In halocyanin, a more complex behavior is observed. About 97% of the rebinding population can also be characterized by a Gaussian distribution of enthalpy barriers at 12 kJ/mol with a width of 6.0 kJ/mol (FWHM). The pre-exponential of this population is 1.6 x 10(12) s-1 at 100 K. After the majority population has rebound, a power-law phase that can be modeled with a gamma-distribution of enthalpy barriers is observed. Between 120 and 180 K, an additional feature that can be interpreted as a relaxation of the barrier distribution toward higher barriers shows up in the kinetics. Above 200 K, a slower, exponential rebinding appears in both cupredoxins. Since the kinetics depend on the NO concentration, this process is identified as bimolecular rebinding.(ABSTRACT TRUNCATED AT 250 WORDS)

Mössbauer spectroscopy on nonequilibrium states of myoglobin: a study of r-t relaxation

A frozen solution of 57Fe-enriched metmyoglobin was irradiated by x rays at 77 K. Mössbauer spectra showed a reduction of Fe(III) high spin by thermalized electrons and a production of a metastable Fe(II) low spin myoglobin complex with H2O at its sixth coordination site. The relaxation of the intermediate was investigated by Mössbauer spectroscopy as a function of temperature and time. The relaxation process starts above 140 K and is fully completed at approximately 200 K. At temperatures between 140 and 200 K, the relaxation lasts for hours and is nonexponential in time. Up to 180 K, the process can be described satisfactorily by a gamma distribution of activation enthalpies with an Arrhenius relation for the rate coefficient. The temperature and time dependence of the Mössbauer parameters indicates structural changes in the active center of the protein as early as 109 K that continue for several hours at higher temperatures. Above 180 K, structural rearrangements involving the whole protein molecule lead to a shift and narrowing of the barrier height distribution.

Ligand binding to anti-fluorescyl antibodies: stability of the antigen binding site

The problem of protein stability is addressed with spectroscopic studies of equilibrium and kinetic properties of the binding of fluorescein to high-affinity monoclonal anti-fluorescyl antibodies (Mab 4-4-20), Fab fragments, and single-chain antibodies (SCA). SCA molecules contain only the variable domains of the antibody and an amino acid linker. The influence of glycerol on the antigen binding reaction is studied by circular dichroism, fluorescence, and absorption spectroscopy. The presence of glycerol in the solvent lowers the affinity of SCA for the ligand drastically, and the affinity even decreases toward lower temperatures. This effect is not observed in Fab and Mab. Analysis of the temperature jump kinetics shows that the dissociation reaction can be modeled as a two-state transition. The CD spectra indicate that the domain structure of the SCA remains unaltered in the presence of glycerol. Therefore, it is concluded that glycerol promotes the dissociation of the two variable domains of SCA. In Fab and Mab, the constant domains provide additional stabilization of the molecular structure at the antigen binding site.

Ligand binding to heme proteins: II. Transitions in the heme pocket of myoglobin

Phenomena occurring in the heme pocket after photolysis of carbonmonoxymyoglobin (MbCO) below about 100 K are investigated using temperature-derivative spectroscopy of the infrared absorption bands of CO. MbCO exists in three conformations (A substrates) that are distinguished by the stretch bands of the bound CO. We establish connections among the A substates and the substates of the photoproduct (B substates) using Fourier-transform infrared spectroscopy together with kinetic experiments on MbCO solution samples at different pH and on orthorhombic crystals. There is no one-to-one mapping between the A and B substates; in some cases, more than one B substate corresponds to a particular A substate. Rebinding is not simply a reversal of dissociation; transitions between B substates occur before rebinding. We measure the nonequilibrium populations of the B substates after photolysis below 25 K and determine the kinetics of B substate transitions leading to equilibrium. Transitions between B substates occur even at 4 K, whereas those between A substates have only been observed above about 160 K. The transitions between the B substates are nonexponential in time, providing evidence for a distribution of substates. The temperature dependence of the B substate transitions implies that they occur mainly by quantum-mechanical tunneling below 10 K. Taken together, the observations suggest that the transitions between the B substates within the same A substate reflect motions of the CO in the heme pocket and not conformational changes. Geminate rebinding of CO to Mb, monitored in the Soret band, depends on pH. Observation of geminate rebinding to the A substates in the infrared indicates that the pH dependence results from a population shift among the substates and not from a change of the rebinding to an individual A substate.

Structural relaxation and nonexponential kinetics of CO-binding to horse myoglobin. Multiple flash photolysis experiments

The geminate recombination kinetics of CO-myoglobin strongly deviates from single exponential behavior in contrast to what is expected for unimolecular reactions (1). At low temperatures, this result was attributed to slowly exchanging conformational states which differ substantially in barrier height for ligand binding. Above 160 K the kinetics apparently slow down with temperature increase. Agmon and Hopfield (2) explain this result in terms of structural relaxation perpendicular to the reaction coordinate, which enhances the activation energy. In their model, structural relaxation homogenizes the kinetic response. Recently, Steinbach et al. (3) proposed a relaxation model which conserves the kinetic inhomogeneity. Below we test these conjectures by single and multiple excitation experiments. This method allows for discrimination between parallel (inhomogeneous) and sequential (homogeneous) kinetic schemes. The kinetic anomaly above 160 K is shown to result from a homogeneous, structurally relaxed intermediate. However a second anomaly is found above 210 K concerning the inhomogeneous phase which may indicate either a shift in activation energy or entropy.