This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.106.093625v1 92/12/4401    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Noronha, M. Articles by Maçanita, A. L. Search for Related Content PubMed PubMed Citation Articles by Noronha, M. Articles by Maçanita, A. L.

Tracking Local Conformational Changes of Ribonuclease A Using Picosecond Time-Resolved Fluorescence of the Six Tyrosine Residues

Melinda Noronha ^* , João C. Lima , Emanuele Paci , Helena Santos  and António L. Maçanita ^*

^* Centro de Química Estrutural, Instituto Superior Técnico, Universidade Técnica de Lisboa, Lisbon, Portugal; Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa, Oeiras, Portugal; REQUIMTE/CQFB, Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, Caparica, Portugal; Department of Physics and Astronomy, University of Leeds, Leeds, United Kingdom

Correspondence: Address reprint requests to António L. Maçanita, Instituto Superior Técnico, Av. Rovisco Pais, 1049-001 Lisboa, Portugal. Tel.: 351-21-8419606; E-mail: macanita{at}ist.utl.pt.

	ABSTRACT

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
The six tyrosine residues of ribonuclease A (RNase A) are used as individual intrinsic probes for tracking local conformational changes during unfolding. The fluorescence decays of RNase A are well described by sums of three exponentials with decay times (₁ = 1.7 ns, ₂ = 180 ps, and ₃ = 30 ps) and preexponential coefficients (A₁ = 1, A₂ = 1, and A₃ = 4) at pH 7, 25°C. The decay times are controlled by photo-induced electron transfer from individual tyrosine residues to the nearest disulphide (–SS–), bridge, which is distance (R) dependent. We assign ₁ to Tyr-76 (R = 12.8 Å), ₂ to Tyr-115 (R = 6.9 Å), and ₃ to Tyr-25, Tyr-73, Tyr-92, and Tyr-97 (all four at R = 5.5 ± 0.3 Å) at 23°C. On the basis of this assignment, the results show that, upon thermal or chemical unfolding only Tyr-25, Tyr-92, and Tyr-76 undergo significant displacement from their nearest –SS– bridge. Despite reporting on different regions of the protein, the concordance between the transition temperatures, T_m, obtained from Tyr-76 (T_m = 59.2°C) and Tyr-25 and Tyr-92 (T_m = 58.2°C) suggests a single unfolding event in this temperature range that affects all these regions similarly.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Fluorescence spectroscopy is widely used as a tool for protein folding studies (1). When multiple tryptophan (Trp) or tyrosine (Tyr) residues are present, the steady-state fluorescence is the sum of the signals from all the fluorescent residues, located in different parts of the protein. Time-resolved fluorescence spectroscopy (TRFS) has the potential to track several Trp/Tyr residues given that the decay times can be attributed to individual residues.

The use of TRFS has been hindered, in part, because several studies have shown the appearance of more than one decay time for proteins with a single Trp (2) or Tyr (3) residue. However, we have recently shown that the fluorescence of Staphylococcus aureus nuclease A (SNase), with a single Trp (4), and that of bovine ubiquitin (UBQ), with a single tyrosine (Tyr) residue (5), decay as a single exponential below the unfolding temperature region. Above this temperature region, the decays become multiexponential due to the fluorescence contribution of the unfolded protein, and from the preexponential coefficients of these decays, the mol fractions of the folded and unfolded proteins could be accurately evaluated and used to calculate equilibrium constants, enthalpies, entropies, and specific heats of unfolding.

In this work, we explored the application of TRFS to discriminate the fluorescence signal of the multichromophoric protein, ribonuclease A (RNase A), and to track local conformational changes resulting from thermal and chemical unfolding. RNase A was selected because it is a benchmark protein for biophysical studies, being one of the most well-characterized proteins (6,7). Furthermore, some studies have employed site-directed mutagenesis to introduce Trp residues in different parts of the protein (8), thus providing information in the site-specific region around the Trp residue in the mutant. However, these substitutions can result in structural changes with respect to the wild-type protein, as indirectly observed on comparison of changes in the thermodynamic properties associated with the unfolding process such as the transition temperature (T_m), enthalpy (H), and entropy (S) (8). Consequently, the possibility of obtaining information on the sites around the six Tyr residues in the unaltered protein would be beneficial.

The results of this work allowed us to assign the fluorescence decay times of the multiexponential decays of RNase A to its different Tyr residues. We also concluded that electron transfer is the mechanism responsible for the observed short decay times. Finally, on the basis of the assignment, we were able to use TRFS to probe the conformational changes, resulting from thermal and chemical unfolding that occur in the different regions where the Tyr residues are located and couple the experimental data to molecular dynamic (MD) simulations.

	EXPERIMENTAL

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Sample preparation
Bovine pancreatic ribonuclease A (M_w = 13,700 Da) was purchased from Boehringer Mannheim (Basel, Switzerland) and from Sigma (ribonuclease A type XII-A; St. Louis, MO) and dialyzed against 10 mM sodium phosphate, pH 7.3. Protein concentration was determined by optical spectroscopy, with an extinction coefficient of 9800 Lmol^–1cm^–1 at 278 nm. Guanidine hydrochloride (99% pure) was purchased from Fluka (Milwaukee, WI). Sodium dihydrogen phosphate and disodium hydrogen phosphate were of analytical grade and were purchased from Merck (Rahway, NJ). N-acetyl-tyrosinamide (NAYA) was purchased from Sigma.

Absorption and fluorescence spectroscopy
Ultraviolet (UV)-absorption spectra were measured on an Olis (Bogart, GA) UV-Vis DW2 double-beam spectrophotometer with 1.0 nm resolution equipped with a circulating water bath and a thermocouple placed inside the cell for temperature control. Steady-state fluorescence excitation and emission spectra were measured using a SPEX (Edison, NJ) Fluorog 212I spectrofluorimeter. All spectra were collected in the sample-to-reference (S/R) mode and corrected for optics, detector wavelength dependence (emission spectra), and lamp intensity wavelength dependence (excitation spectra). Fluorescence was collected in right angle geometry.

Fluorescence quantum yield of RNase A, at pH 7 and at 25°C, was determined by comparison to parent compound N-acetyl-tyrosinamide (NAYA), whose quantum yield was in turn determined using a standard solution of p-terphenyl in cyclohexane (0.77) (9).

Time-resolved fluorescence measurements were carried out using the single-photon counting technique as previously described (10). Excitation of samples was carried out with the frequency-tripled output of an actively mode-locked picosecond Ti-Sapphire laser (Spectra Physics Tsunami; Mountain View, CA), pumped by a solid-state laser (Spectra Physics Millennia Xs). The repetition rate was set to 4 MHz by passage through an optical-acoustic modulator (Pulse Selector 3980, Spectra Physics). Light pulses were monitored with a fast photodiode, filtered with a constant fraction discriminator (Canberra 2126), and used as stop signals in a time-to-amplitude converter (Canberra 2145 TAC). Excitation was vertically polarized and emission was collected at 90° geometry, passed through a polarizer at 54.7° (Spindler & Hoyer (Göttingen, Germany) Glan-Thompson laser prism polarizer, Franklin, MA) and a monochromator (Jobin-Yvon (Longjumeau, France) H20 Vis;), and detected with a microchannel plate photomultiplier (MCP-PT Hamamatsu R3809u-50; Hamamatsu, Japan). Inverted START-STOP configuration was used in the acquisition. The experimental instrumental response function for all excitation wavelengths was in the 38–42-ps range. Alternate collection of pulse profile and sample decays was performed (10³ counts at the maximum per cycle) until 5 x 10³ (typical) to 3 x 10⁴ total detected counts had been accumulated at the maximum of the fluorescence signal. The fluorescence decays were deconvoluted on a PC, using Striker's SAND program (11) (LINUX version), which allows individual and global analysis of the decays with individual shift optimizations.

Molecular orbital calculations
Gas phase ionization potentials (IP) and electron affinities (EA) were calculated from the formation enthalpies of the molecule and of its oxidized and reduced forms, respectively (12). Molecular geometries were first calculated using the MM⁺ molecular mechanics method and then optimized with the semiempirical AM1 method at the UHF level, using HyperChem version 5.0 software by Hypercube (Gainesville, FL) (13). Enthalpies of formation were obtained by AM1 from single point calculations on the optimized geometries.

Molecular dynamics
Simulations were performed with the program CHARMM (14) using the all-atom force-field CHARMM22 (15) and a generalized Born continuum model to include the effects of the solvent (16). Starting from the x-ray crystal structure, 1FS3 (17), we first performed a short steepest descent minimization to remove possible steric clashes, followed by a slow heating to the desired temperatures, and then performed canonical simulations lasting between 2 and 10 ns. The temperatures chosen were 23, 87, and 177°C. At 23°C, the protein remains stable and the root mean-square deviation (RMSD) reaches a plateau at 2.5 Å after 0.5 ns and never exceeds 3 Å during the whole duration of the simulation (4 ns). At the higher temperatures, the RMSD increases, but the protein never completely unfolds due to the presence of the four disulphide bonds. At 177°C, the RMSD from the native structure is, on average, 5.5 Å and never exceeds 7 Å. Langevin dynamics were used with a low friction coefficient to accelerate the sampling relative to a water-like friction. The integration step was 2 fs for all of the simulations performed.

	RESULTS

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
Absorption and steady-state fluorescence of RNase A
The absorption and emission spectra of RNase A at pH 7.2 and 23°C (Fig. 1 a) show maximum absorption at _max = 278 nm and maximum emission at _em = 306 nm. The absorption spectrum is 1 nm blue-shifted with respect to NAYA in the nonpolar solvent dioxane (_max = 279 nm) and 3 nm red-shifted with respect to NAYA in water (_max = 275 nm) (9).

View larger version (15K): [in this window] [in a new window]   FIGURE 1  (a) Absorption () and emission () spectra of RNase A in 10 mM sodium phosphate, pH 7.2, at 23°C. The emission spectra of RNase A at pH 2.2 (darker line) at 23°C show no change in the fluorescence intensity from pH 7 to pH 2.2. (Inset) Fluorescence quantum yield of RNase A in phosphate buffer at 23°C from pH 7.2 to pH 1. (b) Absorption spectra of RNase A at 25 and 80°C show a 2.5-nm blue shift and reduction in molar absorption coefficient () at 80°C due to exposure to water.

Studies carried out by Wills and Szabo (19) demonstrated that the fluorescence of Tyr is quenched in the presence of hydrogen-bond acceptors such as acetate. We have made a similar observation for Tyr in UBQ, whose fluorescence is quenched as a result of proton transfer to a nearby carboxylate group of Glu-51 (5). Since previous studies have suggested that, at pH 5, Tyr-25 of RNase A is hydrogen-bonded to the carboxylate of Asp-14 (7,20), we have analyzed the fluorescence response to the pH change. Upon acidification to pH 2.2, there is no visible change in the fluorescence spectra and intensity relative to pH 7.2 (Fig. 1 a).

Temperature effects on the absorption and fluorescence spectra
At 80°C, the absorption spectrum of RNase A shows a decrease in the molar absorption coefficient and a blue shift of 1 nm, indicating the exposure of buried tyrosine residues (Fig. 1 b). Fig. 2 a shows a plot of the differential absorbance (A) of RNase A measured at = 287 nm and pH 7.2 as a function of temperature (T). From the decrease in the absorbance in the 50–70°C temperature range, a melting temperature of T_m = 62°C was estimated.

View larger version (9K): [in this window] [in a new window]   FIGURE 2  (a) Variation of differential absorbance (A) of RNase A measured at = 287 nm with increasing temperature. (b) Fluorescence quantum yield (F) of RNase A at pH 7 in 10 mM phosphate buffer as a function of increasing temperature.

Denaturant effects on the absorption and fluorescence spectra
Similar to the effect of increasing temperature, the absorption spectra of RNase A show a 7% decrease in the molar extinction coefficient () with increasing concentrations of added GuHCl and a 2-nm blue shift from = 278 nm in the absence of GuHCl to = 276 nm at 4.2 M GuHCl, indicating the exposure of buried tyrosine residues. At the same time, the fluorescence quantum yield increases from _F = 0.013 at 0 M GuHCl to _F = 0.039 at 4.2 M GuHCl, indicating a decrease in the efficiency of the fluorescence quenching.

Time-resolved fluorescence of RNase A
The fluorescence decays of RNase A in 10 mM phosphate buffer at pH 7 (23°C), measured with excitation at _exc = 275 nm and emission at _em = 295 nm, can be fitted with a sum of three exponentials, providing decay times of ₁ = 1.72 ± 0.02 ns, ₂ = 180 ± 15 ps, and ₃ = 30 ± 8 ps (Fig. 3). The errors represent the maximum variation in the values of the decay times obtained from five fluorescence decays measured with different samples and two emission wavelengths.

View larger version (26K): [in this window] [in a new window]   FIGURE 3  Fluorescence decay of RNase A in 10 mM phosphate buffer at pH 7 and at 23°C with a timescale of 24.26 ps/channel is best fit with a sum of three exponentials (a) with decay times 1 = 1.72 ns, 2 = 180 ps, and 3 = 30 ps as seen from the autocorrelation (A.C.), weighted residuals (W.R.), and 2 of 1.01. On normalization of the ratios of preexponential coefficients, a distribution of 1:1:4 is observed. (b) Fluorescence decay of RNase A in 10 mM phosphate buffer at pH 7 and at 23°C with a timescale of 24.26 ps/channel fit with a sum of two exponentials, 2 = 1.49 with decay times 1 = 1.69 ns and 2 = 90 ps and preexponential coefficients of A1 = 1.47 and A2 = 4.53.

(1)

with

Since the fluorescence decays of NAYA are single exponential in low dielectric constant solvents (e.g., mixtures of dioxane:water below 70% v/v water concentration) (9), the fluorescence decay of a protein with six tyrosine residues in different environments could potentially be expected to display up to six different decay times (_i). The respective preexponential coefficients (A_i) would then be proportional to the concentration of each tyrosine residue, given that each decay time can be properly resolved (1,2) (see Discussion). In this case, we observe only three lifetimes and normalization of the preexponential coefficients for the six tyrosine residues present () provides a relative ratio of 0.9:1.0:4.1 for the three decay times. The result indicates that four tyrosines are strongly quenched (₃ = 30 ps), one partially quenched (₂ = 180 ps), and the sixth only slightly quenched (₁ = 1.72 ns). Excitation of RNase A at different wavelengths gave similar results, except when the emission was collected at wavelengths longer than 305 nm. At these wavelengths, a residual 3.5-ns component, clearly due to an impurity, has been detected.

In summary, the foregoing results show that the fluorescence decays of RNase A can be described as triple exponentials and that the preexponential coefficients indicate that four tyrosines have a very short decay time (30–40 ps), a fifth tyrosine a relatively short decay time (165–195 ps), and the last a relatively long lifetime (1.70–1.74 ns).

Temperature and denaturant effects on the time-resolved fluorescence of RNase A
The fluorescence decays of RNase A remain triple exponential up to 80°C. The decay times (_i) and preexponential coefficients (A_i) are plotted in Fig. 4, a and b, as a function of temperature.

View larger version (20K): [in this window] [in a new window]   FIGURE 4  (a) The longest decay time of RNase A 1 (squares) decreases from 1.72 to 0.82 ns, 2 (open circles) increases from 180 to 400 ps, whereas 3 (solid circles) increases from 30 to 80 ps from 23 to 80°C. The values of 1 fall in between the decay times of NAYA in mixtures of 20:80 dioxane/water and neat water (solid circles and solid lines). (b) The normalized preexponential coefficient A1 (squares) is constant within experimental error, A2 (open circles) increases from 1 to 3 and A3 (solid circles) decreases from 4 to 2 from 23 to 80°C. Solid lines shown are theoretical fits assuming the six-Tyr model (Eq. 13), from which values of and were calculated.

Significantly, the preexponential coefficients show an alteration in the relative ratios from 1:1:4 at 23°C to 1:3:2 at 70°C. This indicates that two of the four strongly quenched tyrosines become only partially quenched when the protein unfolds (Fig. 4 b).

The chemical unfolding of RNase A employing GuHCl was also followed using time-resolved fluorescence. From 0 to 4.2 M GuHCl, the long decay time ₁ increases slightly from 1.71 to 1.86 ns, ₂ increases from 165 to 754 ps, and ₃ increases from 34 to 104 ps (Fig. 5 a). The preexponential coefficient A₁ increases from 1 to 1.3, A₂ from 1 to 2.7, and A₃ remains approximately constant at 4 up to 2.6 M GuHCl, and sharply decreases to 2 in the 2.6–4.2 M GuHCl concentration range (Fig. 5 b). As in the case of thermal unfolding, chemical unfolding of RNase A induces a change in the environment of two highly quenched tyrosine residues in the native state, which become partially quenched in the unfolded state.

View larger version (11K): [in this window] [in a new window]   FIGURE 5  (a) Upon increasing the concentration of GuHCl from 0 to 4.2 M, 1 (squares) increases from 1.71 to 1.85 ns, 2 (open circles) increases from 160 to 755 ps, and 3 (solid circles) increases from 33 to 104 ps. (b) The normalized preexponential coefficient A1 (squares) increases from 0.9 to 1.3, A2 (open circles) increases from 1 to 2.7, and A3 (solid circles) decreases from 4 to 2 from 0 to 4.2 M GuHCl. The continuous line for A1 is a fit using a sigmoidal function; the solid lines shown for A2 and A3 are theoretical fits assuming the six-Tyr model from which values of and m = 12.9 kJ/mol/M were determined.

	DISCUSSION

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

(2)

The preexponential coefficients represent the fluorescence intensity of each tyrosine residue at t = 0, and can be used to evaluate the mol fraction of the residue in the excited and ground states after appropriate corrections (4). Briefly, the fluorescence intensity of each tyrosine residue at t = 0, I_i(0), is equal to the concentration of the corresponding excited species [Y_i*](0) multiplied by its radiative rate constant For a decay measured at a given emission wavelength, _em, the respective preexponential coefficient A_i is equal to the fluorescence intensity of the residue () times the fraction of the emission, f_i(_em), that is detected at _em, modulated by an instrumental constant, c_inst (Eq. 3).

(3)

with:

(4)

where is the total fluorescence emission.

Thus, Eq. 3 permits determination of the excited-state mol fractions of the residues when the k_fi and f_i(_em) are known. The radiative rate constant of NAYA does not change from dioxane (3.5 x 10⁷ s^–1) to water (3.6 x 10⁷ s^–1); hence, k_fi of each tyrosine residue is expected to have similar values. Also, the emission spectra of NAYA in dioxane and water show only a very small bathochromic shift (2 nm) and, as such, the normalized spectra are almost coincident. Consequently, the fi of each tyrosine residue can be considered to be the same. Thus, the preexponential coefficients (A_i) are proportional to the excited state concentrations at t = 0 without further correction.

The ground-state concentrations of the tyrosine residues are proportional to the excited-state concentrations at t = 0 divided by their respective molar extinction coefficients (_i) at the excitation wavelength. In the case of NAYA, an 20% change in _max is observed upon going from dioxane to water. However, an isosbestic wavelength exists at = 265 nm (9). As a result, despite the fact that the six tyrosine residues are located in different regions of the protein, with selective excitation at _exc = 265 nm, all the tyrosine residues have the same molar extinction coefficient and absorb the same fraction of the incident excitation intensity. In conclusion, the preexponential coefficients obtained from the fluorescence decays of RNase A with excitation at = 265 nm, are directly proportional to the ground-state concentration of each tyrosine residue with no correction required for the fraction of light emitted (f_i(_em)), radiative rate constant (k_f), or molar extinction coefficient.

In the case of native RNase A, the fluorescence decay is best fit with a sum of three exponential functions. Therefore, the 1:1:4 ratio of the normalized preexponential coefficients () implies that four tyrosine residues are strongly quenched and one partially quenched, whereas the third exhibits a decay time close to that of NAYA in water (slightly quenched). However, the foregoing interpretation relies on the assumption that all tyrosine residues possess single exponential decays.

The fluorescence decay of NAYA in dioxane-water mixtures exhibits single exponential decays over a wide temperature range below 70% water content. Above this water content, a second shorter decay time (400 ps) starts to develop and its normalized preexponential ultimately reaches the value of 0.12 in water (9). Very similar observations were made with the single-tyrosine protein UBQ, but the preexponential value was greater (5). The observation of this double-exponential decay, which results from electron transfer from one of the three cresol rotamers to the protein backbone, depends on three factors: i), the relative mol fractions of the ground-state rotamer population, the rotamer conformations being determined by rotation about the C-C_ß bond (); ii), the rotational rate constants for interconversion between the rotamers, k_r and k_r'; and iii), the photo-induced electron transfer rate constant, k_ET, which depends strongly on solvent polarity (9). In dioxane/water mixtures below 70% water, the decays are single exponentials because k_ET is too small.

The structure of RNase A obtained by x-ray crystallography or by NMR (22) indicates that three tyrosines (Tyr-25, Tyr-73, and Tyr-97) are buried in the interior of the protein, two (Tyr-92 and Tyr-115) are 50% exposed and one (Tyr-76) is 70–80% exposed. Therefore, it is quite reasonable to assume that at least five of the tyrosines will exhibit single-exponential decay at 23°C. Although the sixth or most-exposed tyrosine might potentially show double exponential decay, the fact that the preexponential coefficient of the longest decay time, assigned to this most exposed tyrosine (see below), is approximately equal to unity over the entire temperature range studied demonstrates that all of the tyrosines undergo single-exponential decay. This fact also implies that the mol fraction of rotamers undergoing electron transfer to the backbone is low.

Mechanisms of fluorescence quenching
When electron transfer to the protein backbone would be the sole mechanism of fluorescence quenching, the decay of RNase A at 23°C would be expected to be described by a sum of three exponentials with decay times of ₁ 5 ns, corresponding to three tyrosine residues buried in the interior of the protein, a second decay time ₂ 4 ns corresponding to the 50% exposed tyrosines, and a third decay time (₃) corresponding to the decay time of NAYA in 70–80% water/dioxane (₃ 2 ns). The much shorter decay times observed in the fluorescence decay of RNase A at 23°C clearly show that at least five of the tyrosine residues are highly quenched. Specifically, the shortest decay time of RNase A (₃) indicates that four tyrosines suffer an 150-fold reduction in fluorescence intensity () as a result of quenching.

In the presence of quenching, the reciprocal fluorescence decay time of each individual tyrosine residue (_i) is equal to the sum of the radiative ( k_fi), the radiationless (k_nri), and the quenching (k_qi) rate constants, which leads to Eq. 5, where _0i represents the unquenched lifetime of each Tyr residue, assumed equal to the lifetime of NAYA in the appropriate dioxane/water mixture (9)

(5)

From this equation, values of k_q1 = 7 x 10⁷ s^–1, k_q2 = 5.3 x 10⁹ s^–1, and k_q3 = 3.0 x 10¹⁰ s^–1 are obtained for the pseudounimolecular quenching rate constants of the three groups (1:1:4) of tyrosines, respectively.

Quenching by proton transfer
Among the fluorescence quenching mechanisms known for tyrosine within proteins, quenching by proton transfer to a carboxylate of an aspartic or glutamic amino acid residue has been observed in bovine ubiquitin (5). Studies have suggested that Tyr-25 is hydrogen bonded to Asp-14. As a result, quenching of this tyrosine could occur as a result of excited-state proton transfer (ESPT). To check if ESPT might be responsible for tyrosine fluorescence quenching in RNase A, the fluorescence quantum yield was measured as a function of pH from pH 7 to pH 1 at 23°C (Fig. 1 a, inset). Because ESPT does not occur from tyrosine to protonated carboxylate groups, a substantial increase in the fluorescence yield of RNase A would be expected at sufficiently low pH values if ESPT were the mechanism responsible for fluorescence quenching. The fluorescence quantum yield showed no change between pH 7 and pH 2 and only a 20% increase in fluorescence intensity below pH 2. Previous studies have shown that, at pH 2, RNase A has a T_m = 30°C (23). The lack of a change in fluorescence intensity in the pH range from 7 to 2 (Fig. 1 a, inset) indicates that, even though Tyr-25 may be hydrogen bonded to Asp-14, the quenching of Tyr-25 does not result from proton transfer to nearby carboxylate groups. The modest increase in fluorescence intensity below pH 2 at 25°C is attributed to the pH-induced denaturation of RNase A, as previously stated by Cowgill (24).

Quenching by disulfide bridges
Three decades ago, Cowgill classified the fluorescence of the tyrosine residues in RNase A in the following manner: Tyr-25, Tyr-97, and Tyr-92 as being buried in the interior, hydrogen bonded to carbonyl groups, and completely quenched, whereas Tyr-73, Tyr-76, and Tyr-115 are partially quenched by different mechanisms varying from quenching by the hydrated peptide carbonyl group, by disulphide bonds, and by resonance energy transfer to other quenched tyrosines (25). However, the x-ray crystal structure of the protein shows that Tyr-92 is neither completely buried nor hydrogen bonded to a carbonyl group. Hence, the quenching of this tyrosine residue cannot be described by the mechanism proposed by Cowgill. In an era when time-resolved fluorescence spectroscopy and crystal structure analyses were in their initial stages, his work on this protein remains noteworthy. However, the discrepancy observed for Tyr-92 suggests that the classification adopted for tyrosine residues in RNase A requires revision. In this context, Swadesh et al. subsequently proposed that the quenching of tyrosine in RNase A is a result of complex formation between tyrosine residues and disulphide bridges (26).

Mechanism of quenching by disulfide bridges
Different mechanisms of fluorescence quenching of tyrosine or tryptophan residues by disulphide bridges have been proposed, including excited-state electron transfer (27), energy transfer (28), and short-range interactions that permit vibrational dissipation of the excitation energy (29,30).

The standard Gibbs energy of electron transfer from tyrosine to a disulphide bridge can be estimated from the difference between the oxidation potential of cresol ( versus saturated calomel electrode in acetonitrile) (31) and the reduction potential of cystine ( versus saturated calomel electrode in water) (32) given by Eq. 6, assuming that cresol and cystine are representative models for the donor (tyrosine ring) and acceptor (disulphide bridge), respectively. On the basis of these values it is clear that disulfides are excellent electron acceptors and that the excited-state electron transfer from tyrosine to disulfide bridges is strongly exergonic ( from Eq. 6).

(6)

In Eq. 6, the excitation energy () was taken from the mid-energy value of the absorption and emission spectra of NAYA in dioxane-water mixtures. G⁰() (Eq. 7) accounts for the Coulombic interaction of the ion pair (G_ionic) and corrects the values for transfer from acetonitrile to water (G_solv), in Joule, (33).

(7)

However, the experimental reduction potential of cystine refers to a two-electron, two-proton irreversible reaction and not a single electron transfer as is in the case of the excited-state electron transfer from tyrosine to disulphide bridges in RNase A. Consequently, we also evaluated by using the experimental ionization potential of cresol (IP = 8.33 eV) and experimental electron affinities of model compounds (dimethyl disulfide, EA = 1.87 eV (34) or diethyl disulfide, EA = 1.95 eV (35) (Eq. 8)).

(8)

Once again, the values of (e.g., for dimethyl disulfide) indicate that disulfides are excellent electron acceptors and that the excited-state electron transfer from tyrosine to disulfide bridges is strongly exergonic ( from Eq. 8). In the case of Eq. 8, G_solv corrects the IP-EA values from gas phase to water (i.e., _ref = 1 in Eq. 7).

The values of suggest that excited-state electron transfer in the tyrosine-disulfide bridge should be practically activationless. In fact, although the value of may lie in the Marcus inverted region ( in Eq. 9) (36,37) when a value 1.6 eV is assumed for the reorganization energy (38), it is well known that Eq. 9 (Marcus' classic theory) overestimates the increase in the Gibbs activation energy of electron transfer with decreasing in the inverted energy region (39).

(9)

The absence of a temperature dependence of the two shortest decay times of RNase A in the pretransition temperature range (23–50°C) also indicates a practically activationless quenching of the corresponding five tyrosines. Hence, the electron transfer rate constant (k_ET) from tyrosine to a nearby disulfide bridge should depend essentially on the effective donor-acceptor distance R (Eq. 10).

(10)

In Eq. 10, R is the center-to-center distance, R₀ is the sum of the van der Waals radii of the acceptor and donor, L is the orbital overlap parameter, and k₀ is the preexponential coefficient, given by Eq. 11, where J₀(R₀) is the coupling matrix element at the contact distance R₀.

(11)

The center-to-center distances from the six tyrosines of RNase A to the closest disulfide bridge, obtained from the x-ray crystal structure (Table 1) and MD simulations at 23°C (Fig. 6), show that four tyrosines (25, 73, 92, and 97) are within 5–6 Å, Tyr-115 is at 7 Å and Tyr-76 at 12 Å. Clearly, then, electron transfer will be the dominant deactivation process for the first five tyrosines.

View larger version (22K): [in this window] [in a new window]   FIGURE 6  Comparison of the distances from the center of the tyrosine ring to the center of the disulfide bridge (R) for (a) Tyr-25 from disulfide bridge C26–C84 at 23°C and (b) at 87°C, (c) Tyr-73 from C58–C110 at 23°C and (d) at 87°C, (e) Tyr-92 from C40–C95 at 23°C and (f) at 87°C, (g) Tyr-97 from C40–C95 at 23°C and (h) at 87°C, and (i) Tyr-115 from C58–C110 at 23°C and (j) at 87°C.

h

Under the assumption that is close to zero or simply identical for all tyrosines, the electron-transfer rate constant will depend on three parameters k₀, L, and R (Eq. 10). A plot of ln k_ET vs. R should therefore be linear. Strictly, the experimental kET represents an average kET and not its value at the average R. However, the maximum variation in R for all tyrosines is around 0.3 Å except for Tyr-92, which shows variations up to 0.8 Å. Studies on electron transfer in proteins between redox centers have shown that kET is significantly modulated by the intervening medium (42) where kET is faster between covalently bonded redox centers than through vacuum. On a closer look at the crystal structure, all tyrosines seem equivalent with respect to their medium except for Tyr-115, which has the –CH2– group of one of the cysteines in between the tyrosine ring and the center of the disulphide bridge Cys-58–Cys-110. However, no covalent bond is apparent directly between the tyrosine rings and disulphide bridges from the crystal structure. A logarithmic plot of the k_ET values for the three groups of tyrosine in RNase A versus the tyrosine-disulfide distances provides estimates of the parameters of Eq. 10 of k₀ = 2.4 x 10¹² s^–1 and L = 2.45 Å. These values are consistent with values of these parameters reported for other electron-transfer pairs (42,43).

Reconstruction of fluorescence decays
Employing the values of k₀ and L, we calculated values of the fluorescence decay times _i(R_i) of the six tyrosines and used these to construct a theoretical hexa-exponential decay of RNase A with 38,000 counts at the peak channel. The theoretical fluorescence decay was convoluted with the experimental pulse and then analyzed using sums of two and three exponentials. A sum of three exponentials was found to be sufficient to fit the theoretical hexa-exponential decay (Fig. 7). Thus, the decay times of four of the tyrosines are close enough to each other to be experimentally indistinguishable and appear as a single decay time (₃) with a fourfold larger preexponential coefficient A₃. The decay times of Tyr-115 (₂) and Tyr-76 (₁) are sufficiently different to require another two exponential terms in the fit. The recovered decay times and respective preexponentials of ₁ =1.75 ns (A₁ = 1.05) ₂ =127 ps (A₂ = 0.93) ₃ = 40 ps (A₃ = 4.02) are similar to those obtained from experimental decays.

View larger version (22K): [in this window] [in a new window]   FIGURE 7  Analysis of a theoretical hexa-exponential fluorescence decay of RNase A with 38,000 counts at the maximum. The six decay times were calculated with Eqs. 5 and 10 using values of k0 = 2.4 x 1012 s–1 and L = 2.45 Å and the distances of each tyrosine residue presented in Table 1 at 23°C. The decay can be fitted with a triple exponential with decay times and preexponentials of 1 = 1.75 ns (A1 = 1.05) 2 = 127 ps (A2 = 0.93) 3 = 40 ps (A3 = 4.02).

Thermal unfolding
The study of the parent compound NAYA in water as a function of temperature showed a gradual decrease in the decay times with increasing temperature (Fig. 4 a). In the case of RNase A, ₁ decreased gradually from 20 to 50°C, and sharply from 50 to 60°C, whereas its preexponential coefficient remains constant. The sharper decrease in ₁ reflects the exposure of the slightly buried Tyr-76 to water upon denaturation. This decrease results from the fact that the decay time of the unfolded state (corresponding to Tyr-76 in 95% water). Using the pre- and posttransitional baselines of ₁, we calculated a T_m = 59.2°C. This value essentially reflects the response of the protein to heating in the region where Tyr-76 is located, which is on a loop between ß-strand III and ß-strand IV (Fig. 8 a). Additionally, in the temperature range from 20 to 50°C, we expect k_ET to remain constant or increase slightly with temperature. However, a plot of k_ET, calculated from 1/₁ – 1/₀, where ₀ is the decay time of NAYA in a mixture of 20:80 dioxane/water from 20 to 50°C, shows a reduction in k_ET with increasing temperature, which would only be possible if the distance (R) between the tyrosine residue and the nearest disulfide bridge increases with temperature. No increase in R is apparent from the MD simulations, which predict a constant value of R 12.8 Å for Tyr-76 from 23 to 87°C (Table 1). The decrease in k_ET from 20 to 50°C suggests a pretransition, observed from the decay time ₁ that is not adequately predicted by the MD simulations.

View larger version (77K): [in this window] [in a new window]   FIGURE 8  Stereoscopic views of one of the structures of RNase A from the simulations carried out at (a) 23°C and (b) 177°C.

Summarizing, Tyr-73 and Tyr-97 remain 5.3 Å from the closest disulfide bridge and continue to be strongly quenched, whereas Tyr-25 and Tyr-92 increase their distance from the closest disulfide bridge from 5.4 Å to 7.5 and 6.2 Å, respectively, and are only partially quenched upon denaturation. This suggests that the regions associated with Tyr-25, namely -helix II, and with Tyr-92, in the loop in between ß-strand IV and ß-strand V, lose their structure partially or completely, whereas ß-strand III (Tyr-73) and ß-strand V (Tyr-97) maintain their structure even at 80°C.

Chemical unfolding
In the case of chemical unfolding, the temperature effect on the decay times is absent, leaving only the displacement of the tyrosine residues away from their disulphide bridges and exposure to water as the causes of changes in the decay times. Below 2 M GuHCl, there is no appreciable change in the decay times, whereas above this concentration there is an increase in all the decay times (₁ increases from 1.70 to 1.87 ns, ₂ from 0.16 to 0.76 ns, and ₃ from 30 to 100 ps), indicating a significant increase in some of the tyrosine-disulphide bridge distances (Fig. 5 a).

Previous studies have suggested that the denaturant-induced unfolded state of RNase A exhibits less conformational order than the thermally unfolded state. Thus, MD simulations were carried out at 177°C in an attempt to evaluate the distances from the tyrosine residues to their nearest disulphide bridge (R) in a more expanded unfolded state. Fig. 9 shows an increase in R for Tyr-76 during part of the simulation to 13.5 Å (compared to 12.8 Å at 23 and 87°C), which would correspond to ₁ = 1.86 ns.

View larger version (18K): [in this window] [in a new window]   FIGURE 9  Distances from the center of the tyrosine ring to the center of the disulfide bridge (R) for (a) Tyr-92 from C40–C95 at 23°C, (b) 87°C, and (c) 177°C from molecular dynamics simulation started from the native state.

Thermodynamics of RNase A unfolding from preexponential coefficients
In a previous study on thermal unfolding of ubiquitin, a protein with a single tyrosine, we showed that it is possible to determine the exact mol fractions (x_i) from the preexponential coefficients (A_i) associated with the decay times (_i) of the native (N) and unfolded (U) ubiquitin (5). In the case of ubiquitin, the native state was described by a single decay time, whereas two decay times were observed for the unfolded state. For RNase A, our experimental data show good fits with triple exponentials over the entire temperature range. This implies that the observed decay times are distributions of decay times of the native and unfolded states that are experimentally unresolved. However, the preexponential coefficients contain information on the mol fractions of the native folded (x_N) and unfolded (x_U) states. Thus, instead of the appearance of additional decay times, a change in the ratios of the preexponential coefficients is observed during the transition. Thermal and chemical unfolding both show a change in A₃ from 4 to 2. This means that, when the protein is in the native state (x_N = 1) A₃ = 4, and when it is in the unfolded state (x_N = 0) A₃ = 2. Therefore, A₃ can be expressed in terms of the native molar fraction x_N (Eq. 12), from which the equilibrium constant for unfolding, K_U, can be determined and, in turn, the enthalpy (H) and entropy (S) of unfolding (Eq. 13).

(12)

(13)

On fitting Eq. 13 to the temperature dependence of A₃, values of T_m = 58.2°C, H = 461.4 kJ/mol, and S = 1.39 kJ/mol/K were obtained (see Fig. 4 b). These values are in good agreement with previously published differential scanning calorimetry (DSC) results for thermal unfolding of RNase A at pH 5.5 (T_m = 60.4°C and H = 453.4 kJ/mol) (44).

Likewise, the GuHCl-induced unfolding provided values of C_m = 2.93 M, and m = 12.9 kJ/mol/M, in reasonable agreement with values cited in the literature at pH 6.98: C_m = 2.99 M, and m = 12.9 kJ/mol/M (45).

The interesting aspect of the agreement of the TRFS and DSC thermal unfolding results is that our TRFS values refer to thermal unfolding occurring in the regions of RNase A around Tyr-25 and Tyr-92, whereas the latter were derived from DSC measurements, thus reflecting the change in C_p for the whole protein.

The picture emerging from this observation is that the unfolding of RNase A occurs via a "domino effect", in which all of the pieces undergo the transformation in a sequence but yet in one single process. Consistent with this picture is the fact that the T_m values derived from ₁ vs. T (59.2°C), reflecting Tyr-76 (located on a loop between ß-strands III and IV), from A₃ vs. T (58.2°C) for two tyrosine residues (Tyr-25, located on helix II, and Tyr-92, located on a loop between ß-strands IV and V), from A vs. T (62.2°C) reflecting the exposure to water of the three buried tyrosine residues (Tyr-25, Tyr-73, ß-strand III, and Tyr-97, ß-strand), and from DSC (60.4°C) (the whole protein) are all within the experimental DSC error (46).

Information on the unfolded conformation
There are several studies on the unfolded state of RNase A (for a review see Neira and Rico (6) and Scheraga et al. (7)) that indicate the presence of residual structure. This is confirmed by our results and expected from the presence of the four disulphide bridges, which prevent the complete expansion of the protein.

An interesting result of this work is that, of the four highly quenched Tyr residues, Tyr-73 and Tyr-97 remain in close proximity to their nearest disulphide bridges whereas Tyr-25 and Tyr-92 move away. Because the disulfide bridges nearest to Tyr-73 and Tyr-97 are located in secondary structures in regions of the protein different from the Tyr residues, this would tend to indicate that, at 80°C (above the T_m 60°C), tertiary contacts persist for Tyr-73 and Tyr-97.

On comparison of RNase A structures from the MD simulations at 23°C and 177°C, it is clear that the distance R between Tyr-25 and its nearest disulphide bridge Cys-26–Cys-84 results from the disruption of -helix I, which causes the loop between -helix I and -helix II to shift toward -helix I. In the case of Tyr-92 the partial loss in structure of ß-strand V apparently causes an increase in the flexibility of the loop that contains Tyr-92, which results in a shift away from Cys-40–Cys-95.

Even more curious is the fact that, despite the partial loss in secondary structure of ß-strand V, which contains Tyr-97, this strand still remains in close proximity to Cys-40–Cys-95. The C-terminal part of RNase A, which consists of Tyr-73 and Tyr-115, keeps its secondary structure, which maintains the two tyrosine residues in close proximity to their disulphide bridges. The last observation is in agreement with previous studies, which have indicated that Tyr-73 is located in a hydrophobic core of exceptional stability against unfolding.

	CONCLUSIONS

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
The fluorescence decay of Ribonuclease A at pH 7 in the 20–80°C temperature range can be fit with triple-exponential decay functions despite the presence of six tyrosines. Below the unfolding temperature, four tyrosines are equally strongly quenched, one is partially quenched, and the sixth only slightly quenched. Upon unfolding two out of the four strongly quenched tyrosines become less quenched.

The quenching mechanism of these tyrosines in RNase A is ascribed to photo-induced electron transfer from the tyrosine residues to nearby disulphide bridges. The electron-transfer rate constant k_ET is exponentially dependent on the distance R from each tyrosine to the nearest disulphide bridge, using R-values obtained from the x-ray crystal structure and from the dynamic picture of the native and unfolded states provided by all-atom molecular dynamic simulations at 23, 87, and 177°C.

The three decay times were assigned as follows: ₁ is the decay time of Tyr-76 (at R = 12.8 Å), ₂ is the decay time of Tyr-115 (at R = 6.9 Å), and ₃ is the mean decay time of Tyr-25, Tyr-73, Tyr-92, and Tyr-97 (all at R = 5.5 ± 0.3 Å). Upon unfolding, the distances from two of the tyrosines (Tyr-25 and Tyr-92) to the nearest disulfide bridge increase from 5.5 Å to 7 Å and ₂ becomes the average decay time of Tyr-115, Tyr-25, and Tyr-92, whereas ₃ is the average decay time of Tyr-73 and Tyr-97.

From the preexponential coefficients of the fluorescence decays, the unfolding equilibrium constants (K_U) of RNase A were calculated, as a function of temperature or denaturant (GuHCl) concentration, providing values of T_m = 58.2°C and C_m = 2.93 M. The successful application of this method to the unfolding of RNase A rests on the fortunate case that, of the six tyrosine residues present, only two change significantly during unfolding.

This work establishes a new methodology of using TRFS as a tool for protein studies. The discrimination between and assignment of the six tyrosines allowed us to follow the unfolding in specific and different regions of the protein, evidencing an emerging picture of a single unfolding event that affects all the protein regions evenly.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
In memory of his outstanding scientific and human qualities, and for being a source of inspiration for us, the authors dedicate this work to Prof. António V. Xavier who passed away on May 7, 2006.

The work was supported by the Fundação para a Ciência e a Tecnologia (FCT), Portugal, Projects POCI/QUI/56585/04 and POCI/BIA-PRO/57263/04 and by the European Commission, 5^th Framework Program contract QLK3-CT-2000-00640. M. Noronha acknowledges a PhD grant from FCT, SRFH/BD/9096/2002.

Submitted on July 19, 2006; accepted for publication February 1, 2007.

	REFERENCES

TOP ABSTRACT INTRODUCTION EXPERIMENTAL RESULTS DISCUSSION CONCLUSIONS ACKNOWLEDGEMENTS REFERENCES

 
1. Royer, C. A. 2006. Probing protein folding and conformational changes with fluorescence. Chem. Rev. 106:1769–1784.[CrossRef][Medline]

2. Beecham, J. M., and L. Brand. 1985. Time-resolved fluorescence of proteins. Annu. Rev. Biochem. 54:43–71.[CrossRef][Medline]

3. Ferreira, S. T., L. Stella, and E. Gratton. 1994. Conformational dynamics of bovine Cu, Zn superoxide-dismutase revealed by time-resolved fluorescence spectroscopy of the single tyrosine residue. Biophys. J. 66:1185–1196.[Abstract/Free Full Text]

4. Faria, T. Q., J. C. Lima, M. Bastos, A. L. Maçanita, and H. Santos. 2004. Protein stabilization by osmolytes from hyperthermophiles: effect of mannosylglycerate on the thermal unfolding of recombinant nuclease a from Staphylococcus aureus studied by picosecond time-resolved fluorescence and calorimetry. J. Biol. Chem. 279:48680–48691.[Abstract/Free Full Text]

5. Noronha, M., J. C. Lima, M. Bastos, H. Santos, and A. L. Maçanita. 2004. Unfolding of ubiquitin studied by picosecond time-resolved fluorescence of the tyrosine residue. Biophys. J. 87:2609–2620.[Abstract/Free Full Text]

6. Neira, J. L., and M. Rico. 1997. Folding studies on ribonuclease A, a model protein. Fold. Des. 2:R1–R11.[CrossRef][Medline]

7. Scheraga, H. A., W. J. Wedemeyer, and E. Welker. 2001. Bovine pancreatic ribonuclease A: oxidative and conformational folding studies. Methods Enzymol. 341:189–221.[Medline]

8. Navon, A., V. Ittah, J. H Laity, H. A. Scheraga, E. Haas, and E. E. Gussakovsky. 2001. Local and long-range interactions in the thermal unfolding