help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hellings, M.
Right arrow Articles by Engelborghs, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hellings, M.
Right arrow Articles by Engelborghs, Y.
Biophysical Journal 85:1894-1902 (2003)
© 2003 The Biophysical Society

The Dead-End Elimination Method, Tryptophan Rotamers, and Fluorescence Lifetimes

Mario Hellings *, Marc De Maeyer {dagger}, Stefan Verheyden *, Qiang Hao {ddagger}, Els J. M. Van Damme §, Willy J. Peumans ¶ and Yves Engelborghs *

* Laboratory of Biomolecular Dynamics, Catholic University of Leuven, Heverlee, Belgium; {dagger} Laboratory of Biomolecular Modelling, Catholic University of Leuven, Heverlee, Belgium; {ddagger} Institute of Biosciences and Technology, Texas A&M University, System Health Sciences Center, Houston, Texas; § Department of Molecular Biotechnology, University of Ghent, Ghent, Belgium; and Laboratory of Phytopathology and Plant Protection, Catholic University of Leuven, Heverlee, Belgium

Correspondence: Address reprint requests to Dr. Yves Engelborghs, Katholieke Universiteit Leuven, Laboratory of Biomolecular Dynamics, Celestijnenlaan 200D, 3001 Heverlee, Belgium. Tel.: 321-632-7160; Fax: 321-632-7974; E-mail: yves.engelborghs{at}fys.kuleuven.ac.be; or Dr. Marc De Maeyer, Katholieke Universiteit Leuven, Laboratory of Biomolecular Modelling, Celestijnenlaan 200D, 3001 Heverlee, Belgium. Tel.: 321-632-7521; Fax: 321-632-7974; E-mail: Marc.DeMaeyer{at}fys.kuleuven.ac.be.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The Dead-End Elimination method was used to identify 40 low energy microconformations of 16 tryptophan residues in eight proteins. Single Trp-mutants of these proteins all show a double- or triple-exponential fluorescence decay. For ten of these lifetimes the corresponding rotameric state could be identified by comparing the bimolecular acrylamide quenching constant (kq) and the relative solvent exposure of the side chain in that microstate. In the absence of any identifiable quencher, the origin of the lifetime heterogeneity is interpreted in terms of the electron transfer process from the indole C{varepsilon}3 atom to the carbonyl carbon of the peptide bond. Therefore it is expected that a shorter [C{varepsilon}3-C=O] distance leads to a shorter lifetime as observed for these ten rotamers. Applying the same rule to the other 30 lifetimes, a link with their corresponding rotameric state could also be made. In agreement with the theory of Marcus and Sutin, the nonradiative rate constant shows an exponential relationship with the [C{varepsilon}3-C=O] distance for the 40 datapoints.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Fluorescence studies of proteins often provide useful information regarding the environment and the dynamics of the region surrounding naturally occurring fluorescent amino acids like tryptophan. However, this necessitates an exact knowledge of the relation between the fluorescent properties and the protein environment. This is often problematical because the fluorescence decay of a single tryptophan residue in a protein is usually multiexponential. This multiexponential decay is either the consequence of excited state reactions or of the existence of different tryptophan microconformations (Szabo and Rayner, 1980Go; Bajzer and Prendergast, 1993Go; Brown and Royer, 1997Go; Dahms et al., 1995Go; Ababou and Bombarda, 2001Go). Hudson proposed a model that involves reversible ionization of the excited tryptophan residue due to collisional transfer of an electron to a neighboring residue (Hudson, 1999Go). Other groups have also suggested the involvement of an electron transfer process as the principal quenching mechanism in proteins and peptides (Antonini et al., 1997Go; Chen and Barkley, 1998Go; Ababou and Bombarda, 2001Go). In the absence of quenching side chains, the carbonyl group of the peptide itself is thought to be the responsible quencher (Ricci and Nesta, 1976Go; Chen et al., 1996Go; Adams et al., 2002Go). The carbonyl group of the peptide itself will act as an acceptor with the indole C{varepsilon}3 atom of tryptophan as donor (Sillen et al., 2000Go). According to Milner-White (1997)Go, the carbon of the peptide carbonyl contains a net-positive charge, and is therefore available to accept an electron. The dependency of electron transfer on the distance between the acceptor and the donor, as well as on the medium in between, are supposed to be responsible for different lifetimes for different rotamer positions of tryptophan (Chang et al., 1983Go; Goldman et al., 1995Go; Sillen et al., 2000Go; Ababou and Bombarda, 2001Go; Adams et al., 2002Go).

Sillen et al. (2000)Go presented a method that allowed the calculation of the lifetimes of tryptophan residues on the basis of spectral and structural data. In their article the authors determined an exponential relation between the nonradiative rate constant knr and the distance between the indole C{varepsilon}3 atom and the carbonyl carbon of the peptide bond. This suggests a mechanism of electron transfer as the main determinant for the value of the nonradiative rate constant. This relation was based on the equation derived for electron transfer by Marcus and Sutin (1985)Go. The authors used a molecular dynamics simulation map to identify the different rotamers of tryptophan. Further exploration of this method has to be done.

Here we report a new method to identify the different microconformations of tryptophan together with a method to link the rotamers with their respective lifetime on the basis of the correlation between the values of the rate constants for quenching with an external quencher (kq) and the rotamer surface exposed to the solvent. In this study we use the dead-end elimination (DEE) method to determine the possible rotamers of tryptophan. This method is able to find the global optimal arrangement of a collection of side chains attached to a fixed main-chain structure (De Maeyer et al., 2000Go). The DEE method detects and eliminates iteratively those rotamers that cannot be members of the Global Minimum Energy Conformation (GMEC; Desmet et al., 1992Go). The rotamer elimination occurs on the basis of energy criteria balancing rotamers against each other by comparing their main-chain interaction energy and a lower or higher limit for their interactions with the other side chains of the protein (De Maeyer et al., 2000Go). In this study the DEE method is applied on eight different proteins, whose single tryptophan mutants are available. The fact that the tryptophans of the eight proteins are buried inside the molecule strengthens the use of the DEE method to determine the possible rotamers of the Trp-residue. Indeed, it has been proven that the DEE method is quite accurate for buried residues (De Maeyer et al., 1997Go). The time-resolved fluorescence properties of the single tryptophan mutants of seven of these proteins have already been extensively investigated. These proteins, their tryptophan residues and the reference that describes them are listed in Table 1.


View this table:
[in this window]
[in a new window]
 
TABLE 1  The used proteins, together with the position of the tryptophan residues, the PDB:IDs of the starting structures and mutations that have been made in the starting structures

 
The eighth protein we use is trichosanthin (TCS). Trichosanthin from Trichosanthes kirilowii, a type I ribosome-inactivating protein is a monomeric 247-residue, 27,128 Da globular protein (Collins et al., 1990Go; Xiong et al., 1994Go). It contains a single tryptophan residue at position 192 that is buried inside the molecule (Zhang et al., 1993Go). The steady-state fluorescence properties of the Trp192 have been reported previously (Hao et al., 1995Go). Here we report the time-resolved fluorescence properties of TCS and a detailed quenching study with acrylamide to investigate the validity, in this case, of the proposed Hudson model. The values of kq, which also reflect the accessible surface area (ASA) of tryptophan, provided us with a procedure to link the different rotamer clusters with their respective lifetime.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Materials
Root tuber of Trichosanthes kirilowii (TianHuaFen) were obtained from a local drug store in China. All the chromatography materials were purchased from Pharmacia LKB. Acrylamide and n-acetyltryptophanamide (NATA) were purchased from Sigma Chemical (St Louis, MO). All other chemicals were of analytical grade. Distilled and deionized water was used throughout the experiments.

Protein purification
Trichosanthin was essentially purified as described previously (Zhang et al., 1993Go). The purity of the protein was examined by SDS-PAGE. The activity of the protein was determined by the method of Endo et al. (1987)Go by running the targeted RNA on a 1,2% agarose gel.

Multifrequency phase fluorometry
Fluorescence lifetimes were measured using automatic multifrequency phase fluorometry between 1.6 MHz and 1 GHz as described previously (Sillen et al., 2000Go). NATA in ddH2O with a fluorescence lifetime of 3.059 ns (at 22°C) was used as reference. Data analysis was performed using a nonlinear least-squares algorithm (Bevington, 1969Go). Measurements performed at different emission wavelengths (320–380 nm) were analyzed simultaneously with global analysis to improve the recovery of lifetimes and amplitude fractions (Beechem et al., 1983Go). The data were fitted using the modified Levenberg-Marquardt algorithm (More and Sorensen, 1983Go) assuming that the fluorescence lifetimes are independent of the wavelength. This procedure allows the recovery of the fractional intensities with good accuracy.

Fluorescence quenching
The quenching of trichosanthin with increasing amounts of quencher was performed by adding aliquots of a freshly prepared 1-M acrylamide stock solution to the cuvette after which the changes in the fluorescence lifetimes and corresponding amplitude fractions were monitored. Fluorescence lifetime as function of quencher concentration [Q] was fitted by the classical Stern-Volmer equation,

(1)
where {tau}0 ({tau}) is the lifetime in absence (presence) of quencher.

Starting structures
For each protein a full hydrogen model and energy minimized x-ray structure was used as the starting structure. For each structure only the monomer is taken into account. Some residues like Cys, HisH+, Tyr, and Cystine are able to strongly quench fluorescence (Vos and Engelborghs, 1994Go; Chen and Barkley, 1998Go). Fluorescence data are available for single Trp-mutants of the eight proteins in which the neighboring quenching residues are mutated in nonquenching residues. Because we only take electron transfer to the peptide carbonyl into account, it is necessary to mutate in silico these neighboring quenching residues in the starting structures, so that it is possible to link the modeled structures with the experimental data. Table 1 shows the PDB:IDs of the starting structures and the mutations that have been made.

Dead-end elimination (DEE) method
To determine the possible conformations of tryptophan we used the dead-end elimination method. This method allows the determination of the GMEC of a side-chain collection in a fixed backbone configuration, using a side-chain rotamer library (Desmet et al., 1992Go). Tryptophan was placed in 1296 different starting values (for each {chi}-angle every 10°). Throughout this article, {chi}1 is defined by the internal rotation around bond C{alpha}-Cß and {chi}2 by the bond Cß-C{gamma}. Within a shell of 8 Å around the fixed tryptophan, all side-chain residues were allowed to evolve. For each of the starting structures the GMEC was calculated with DEE. The backbone is kept fixed during the calculation. The rotamer library we used was an enhanced version of the libraries of De Maeyer et al. (1997)Go. The energy parameters are based on the CHARMM force field (Brooks et al., 1983Go) but with reduced Van der Waals radii to account for small backbone movement (Gordon et al., 1999Go). For each GMEC, the nonbonded energy of tryptophan (ETrpnb) in relation to the remaining structure of the protein and the total nonbonded energy (Etotnb) of the whole modeled protein was computed. The distance between the indole C{varepsilon}3 atom and the carbonyl carbon of the peptide bond is also determined. The DEE algorithm has been implemented in the Brugel package (Delhaise et al., 1984Go). All calculations were performed on a single processor of a 4-processor Silicon Graphics Origin 2100 (R10000, 250 MHz).

Solvent-accessible area calculation
The solvent-accessible surfaces were calculated with the Survol option in the Brugel software package (Delhaise et al., 1984Go), using a probe sphere with a standard radius of 1.4 Å (Alard et and Wodak, 1991Go).

Determination of the radiative rate constant
Quantum yields were determined relative to tryptophan in water according to the method of Parker and Rees (1960)Go:

(2)
where {int}I is the integrated intensity over the wavelength region 300–450 nm. A is the absorbance at 295 nm, and the quantum yield QTrp for tryptophan in water is taken as 0.14 (Kirby and Steiner, 1970Go).

The average radiative rate constant is calculated by dividing the quantum yield by a wavelength-independent amplitude average lifetime (Sillen and Engelborghs, 1998Go),

(3)
with {sum}{tau}i{alpha}i the average lifetime, {alpha}i is a wavelength-independent amplitude fraction and is defined as (Sillen and Engelborghs, 1998Go),

(4)
The fluorescence intensity at time zero I0 ({lambda}) of each lifetime is integrated over the wavelength region 300–450 nm and then normalized.

Calculation of the nonradiative rate constant
The nonradiative rate constant (knr) is calculated as

(5)
where kr is the radiative rate constant, which is calculated as shown before, and {tau} is the fluorescence lifetime.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Identification of the tryptophan rotamers
The total nonbonded energies of the 1296 structures, as determined earlier, are sorted and only those with a negative Etotnb value are kept. This removes from the list of candidates all those with bad contacts. For the remaining structures a rough sorting is performed allowing only internal clashes below 36 kcal/mol. These structures might contain small, short contacts due to the discretization of the rotamer library. To alleviate these repulsions, the structures are further subjected to 100 steps of steepest descent energy minimization with fixed backbone. Because we are only interested in energetically low tryptophan rotamers, the remaining structures are sorted on basis of their ETrpnb value as shown in Fig. 1 for the case of Trp130 of Colicin A. In agreement with Street and Mayo, (1999)Go we eliminate those structures where the Van der Waals energy of any atom pair exceeds a threshold of 3 kcal/mol, this is around the 125th structure. The remaining structures represent approximately the best 25% of the energetically possible structures and provide us with clusters of Trp rotamer positions (see further). This procedure is repeated for each Trp of the eight modeled proteins.



View larger version (9K):
[in this window]
[in a new window]
 
FIGURE 1  ETrpnb profile of Trp130 of Colicin A.

 
Dead-end elimination
The rotamer clusters found for Trp130 of Colicin A are shown in Fig. 2. For each of the three observed clusters the mean distance between the indole C{varepsilon}3 atom and the C atom of the peptide carbonyl is calculated. Table 2 shows the number and range of the rotamer clusters found for each modeled structure together with their mean distance between the indole C{varepsilon}3 atom and the carbonyl carbon of the peptide bond. In most cases we find for tryptophan an equal number of rotamer clusters and fluorescence lifetimes.



View larger version (12K):
[in this window]
[in a new window]
 
FIGURE 2  Rotamer clusters found for Trp130 of Colicin A, using the Van der Waals cutoff value of 3 kcal/mol.

 

View this table:
[in this window]
[in a new window]
 
TABLE 2  Rotamer clusters for the tryptophans of the eight used proteins

 
Fluorescence lifetimes
The fluorescence decay of Trp192 in native TCS was measured at emission wavelengths ranging from 320 nm to 380 nm in 10-nm intervals. The decay curves were analyzed by single-curve analysis of the data at each wavelength and global analysis of data for seven emission wavelengths. The sum of two exponential terms was adequate to fit the data and yielded lifetime components of 1755 ns and a longer one of 4016 ns. Table 3 shows the results of the global fit for TCS.


View this table:
[in this window]
[in a new window]
 
TABLE 3  Fluorescence lifetimes and amplitude fractions of the Trp192 in trichosanthin; the measurements were performed at 20°C

 
Fluorescence quenching of TCS with acrylamide
The time-resolved fluorescence quenching by acrylamide was monitored to investigate the nature of the quenching mechanism. The effect of addition of acrylamide was not the same on the different lifetime components. Table 3 shows the results of quenching TCS with acrylamide. The two lifetimes show a linear dependence for the acrylamide concentration. Stern-Volmer plots for acrylamide quenching of the TCS lifetime components are shown in Fig. 3. The values for the bimolecular quenching constant kq are 1.3 M-1 ± 0.1 M-1 ns-1 and 0.137 ± 0.06 M-1 ns-1 for {tau}1 and {tau}2, respectively. The value of kq for the short lifetime component is 10-fold larger than that of the longer one, indicating that the short lifetime component is more solvent-accessible. For a fully exposed tryptophan residue (e.g., free Trp in water), kq can be as high as 6.4 M-1 ns-1 (Vos et al., 1995Go). The ratio between {alpha}1 and {alpha}2 is shown in Table 4 and remains practically constant during the quenching experiment.



View larger version (9K):
[in this window]
[in a new window]
 
FIGURE 3  Acrylamide quenching of tryptophan fluorescence of TCS. (A) Stern-Volmer plot for the shortest lifetime component. (B) Stern-Volmer plot for the longest lifetime component.

 

View this table:
[in this window]
[in a new window]
 
TABLE 4  Ratio between the amplitude fractions

 
Accessible surface area
The ASA for a fully exposed Trp is 264 Å. Table 5 shows the surface-accessible areas for every Trp in the x-ray structure. The ASA for every rotamer cluster of Trp192 of TCS and the three tryptophans of colicin A are shown in Table 6.


View this table:
[in this window]
[in a new window]
 
TABLE 5  Accessible surface areas of each tryptophan residue in the PDB structure

 

View this table:
[in this window]
[in a new window]
 
TABLE 6  ASA for W192 of TCS, W86, W130, and W140 of Colicin A

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The fluorescence decay of the single tryptophan residue in trichosanthin is double-exponential. For T4-lysozyme, a single Trp protein which also shows a double-exponential decay, Hudson (1999)Go proposed a new model. In this model it is assumed that reversible electron transfer can occur between tryptophan and a neighboring amino acid side chain. This mechanism gives rise to two exponentials: Eq. 6, where the amplitudes are determined in a complex way by the different rate constants involved; and Eq. 7.

(6)

with

W* is the excited tryptophan residue which decays to the ground state W by both radiative and nonradiative processes (kA0 = kr + knr). kAB and kBA are respectively the rate constants of the forward and reverse electron transfer reaction. Q is a collisional electron acceptor and H is either the ground state of tryptophan or another low-energy species (e.g., the triplet state; Hudson, 1999Go):

(7)

When we fill in the results of the fluorescence lifetime measurements of trichosanthin with no acrylamide present, we find a value of 0.2926 for (kA0 + kAB). If Hudson's model would be applicable on TCS, then adding an external quencher Qext, e.g., acrylamide would change the ratio of the amplitudes, as in

(8)
with kqA the bimolecular quenching constant of acrylamide on the excited state of tryptophan.

The quenching experiments we performed provide us with values for k1, k2, and kqA. Together with the value of 0.2926 for (kA0 + kAB), we can calculate the ratio of the amplitudes for each concentration of acrylamide, in which case Hudson's model would be valid for TCS. These values are shown in Table 4. We clearly see that the ratio of the amplitudes should raise upon addition of acrylamide. However, in our experiments the ratio of the amplitudes remains practically constant (see Table 4). We therefore can exclude reversible electron transfer in this case. We can explain the results by the existence of different rotamers of Trp, where the amplitudes reflect populations of different rotamers. These ground-state rotamers are the result from rotations about the C{alpha}-Cß and Cß-C{gamma} bonds. Each rotamer has a specific local environment with different distances and/or orientations to quencher groups, yielding distinct nonradiative decay rates. Since the Trp fluorescence is not damaged upon radiation, electron transfer has to be reversible. The fact that the Hudson model is not applicable for this protein indicates that the reverse reaction is much slower than the nanosecond timescale. Usually reverse electron transfer reactions are orders-of-magnitude smaller than forward reactions (Mataga et al., 1988aGo,bGo; Paddon-Row et al.,1988Go). For the other proteins we did not eliminate the possibility of fast reversible electron transfer, but on the basis of the reasoning above we think it is very unlikely.

Except for the case of the two tryptophans in NSCP, the number of rotamer clusters found for a tryptophan residue is equal to the number of Trp fluorescence lifetimes. The challenge is to find a procedure to link the rotamer clusters with their respective lifetimes. In this case we are able to link 10 lifetimes with 10 rotamers. To link to the rotamer clusters with their fluorescence lifetimes, we calculated the accessible surface area of the clusters and compared them with the values of the bimolecular quenching constants. In the case of trichosanthin, we find two rotamer clusters, one with an ASA of 13.5 Å2 and another with an ASA of 3.0 Å2, as shown in Table 6. The quenching experiments with acrylamide yield two bimolecular quenching constants and indicate that the short lifetime is more solvent-accessible. According to these results the cluster with the highest ASA should be linked with the short fluorescence lifetime of tryptophan. If we compare the accessible areas of the tryptophans of Colicin A (Table 6) with the bimolecular quenching constants found in quenching experiments with acrylamide of these tryptophans (Table 3, Vos et al., 1995Go) we see the same similarity as in the case of trichosanthin. In this set of 10 rotamers we also observed that the short lifetimes correlate with a short distance between C{varepsilon}3 and the peptide carbonyl C atom, in agreement with the mechanism of electron transfer proposed by Marcus and Sutin. Because of this observed relation, we linked the 40 found rotamer clusters for Trp of the eight proteins with their respectable lifetimes in the same way. In the future, additional quenching experiments of those proteins where up to now no acrylamide quenching has been performed, can further confirm the match between both methods.

Recently, the quantum yield was correlated to the environment of the tryptophan in the protein, i.e., the electric field, by using quantum mechanical calculations (Callis and Vivian, 2003Go). Although these calculations provide the most general answer, we try here to find out if the exponential relation for electron transfer, as proposed by Marcus and Sutin (1985)Go, and as attempted before (Sillen et al., 2000Go), can be applied. It should be noted that the tryptophan residues we study in this sample are selected for the absence of special features like quenching groups, charge, stacked, or perpendicular aromatic groups.

The radiative rate constant is experimentally determined for all the proteins from the quantum yield and the average fluorescence lifetime. The knowledge of both the fluorescence lifetime and radiative rate constant for each rotamer cluster makes it possible to calculate the nonradiative rate constant for each rotamer cluster with Eq. 5.

These nonradiative rate constants are fitted to their respective distance between C{varepsilon}3 of Trp and the carbonyl carbon of the peptide bond by the following equation, derived for electron transfer by Marcus and Sutin:

(9)
where ket (R) is the rate for electron transfer, k0 is the rate of electron transfer at the Van der Waals contact distance (R0 = 3 Å), ß is the range parameter that is depended on the medium. Fluorescence lifetime measurements show an experimental error, and because knr is calculated on the basis of these lifetime measurements, this error will stay. Therefore, we use a weighted-fit with weight factors given as follows:

(10)
with Sknr the standard deviation of knr, as shown in Eq. 10,

(11)
And with S{tau} the standard deviation of the fluorescence lifetime {tau} and Skr the standard deviation of the radiative rate constant kr. This fit yields a value of 5.046 ± 0.335 ns-1 for k0 and a value of 1.497 ± 0.074 Å-1 for ß and is shown in Fig. 4.



View larger version (9K):
[in this window]
[in a new window]
 
FIGURE 4  The nonradiative rate constant as function of the distance between the C{varepsilon}3 atom of Trp and the C atom of the peptide carbonyl. {blacksquare}, 10 datapoints for which the relation between kq and ASA could be determined. •, 30 datapoints assigned on the basis of the relation of short-lifetime, short C{varepsilon}3-C=O distance.

 
If we compare these values with the value of 25 ± 3 ns-1 for k0 and 1.9 ± 0.1 Å-1, determined by Sillen et al. (2000)Go (with only four data points), we see that both parameters are smaller. Recently Laboulais et al. (2001)Go simulated the fluorescence decay of the tryptophans in the HIV-1 integrase catalytic core. Among the many possible quenching processes, they only took the quenching by the peptide bond into account using the formula of Sillen et al. (2000)Go. In their work, the best fits were obtained with a value for ß ~ 1.9 Å-1 in agreement with Sillen et al. (2000)Go. In contrast the electron transfer rate constant k0 had to be lowered from 25 ns-1 (Sillen et al., 2000Go) to a value ~3.7 ns-1. Adams et al. (2002)Go recently published a value of 330 ns-1 for the pre-exponential factor, B of NATA, using a ß-value of 1.7 Å-1 (Beratan et al., 1991Go), and assuming that there is no orientation dependency. If they take the orientation dependence into account, this value raises to 1100 ns-1. Translation of this pre-exponential factor B (with R0 = 0 Å) to a k0 (with R0 = 3 Å) value is what we use in formula 4, which gives a value of 2.03 ns-1 for k0 with no orientation dependency and a value of 6.1 ns-1 for k0 with orientation dependency. The k0 value of 5 ns-1 that we determined approximates the k0 values determined by Laboulais et al. (2001)Go and Adams et al. (2002)Go. The experimental correlation that we observe in this way, linking the nonradiative rate constant and the rotamer conformation, can be used to predict the lifetime in the absence of any special feature. It is clear that only a quantum mechanical calculation can give a more general answer.

For Ribonuclease T1, using molecular mechanics (Ababou and Bombarda, 2001Go) only identified one tryptophan rotamer and therefore these authors rejected the rotamer model. This in contrast with the two rotamers revealed here by the DEE method. Contrary to other methods like energy minimization, molecular mechanics, or molecular dynamics, the DEE method explores the complete conformational space and detects the GMEC. The classic method of an energy map cannot explore all possible rotamer combinations by brute force. As an example, a Trp in contact with six other residues (having only two rotatable {chi}-angles) and allowing {chi}-angle deviations in steps of 10° explodes into 6.1021 possible conformations. If one rotamer calculation takes 1 ms, the total exploration of the conformational space would require 2.1011 years. In addition, molecular dynamics runs are in general too short to be sure that all possible energy minima are visited. If we compare the energy minima we find with DEE for Trp57 in NSCP with those determined by a molecular dynamics simulation performed by Sillen et al. (2000)Go, we can see that the two most populated conformations are found with both techniques.


    CONCLUSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
In conclusion it can be stated that for 10 lifetimes we were able to obtain a fully objective identification of the corresponding rotamer, on the basis of the correlation between the quenching constant (kq) and the accessible surface. For these points of data the relation holds that a short lifetime corresponds to short C{varepsilon}3-C=O distance, suggesting the electron transfer mechanism. Therefore the same correlation was used to identify the rotamers of 30 additional lifetimes. The whole set of 40 datapoints allows us to determine a good average of k0 and ß. Only for the tryptophans of NSCP is an additional lifetime found, which might be an indication for additional sources of heterogeneity.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
This research is supported by the National Fund for Scientific Research of the Flanders (project G-0092-01).

Submitted on November 17, 2002; accepted for publication May 30, 2003.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 CONCLUSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Ababou, A., and E. Bombarda. 2001. On the involvement of electron transfer reactions in the fluorescence decay kinetics heterogeneity of proteins. Protein Sci. 10:2102–2113.[Abstract/Free Full Text]

Adams, P. D., Y. Chen, K. Ma, M. G. Zagorski, F. D. Sönnichsen, M. L. McLaughin, and M. D. Barkley. 2002. Intramolecular quenching of tryptophan fluorescence by the peptide bond in cyclic hexapeptides. J. Am. Chem. Soc. 124:9278–9286.[Medline]

Alard, P., and S. J. Wodak. 1991. Detection of cavities in a set of interpenetrating spheres. J. Comput. Chem. 12:918–922.

Antonini, P. S., W. Hillen, N. Ettner, W. Hinrichs, P. Fantucci, S. M. Doglia, J.-A. Bousquet, and M. Chabbert. 1997. Molecular mechanics analysis of Tet repressor Trp43 fluorescence. Biophys. J. 72:1800–1811.[Abstract/Free Full Text]

Bajzer, Z., and F. G. Prendergast. 1993. A model for multi-exponential tryptophan fluorescence intensity decay in proteins. Biophys. J. 65:2313–2323.[Abstract/Free Full Text]

Beechem, J. M., J. R. Knutson, J. B. A. Ross, B. W. Turner, and L. Brand. 1983. Global resolution of heterogeneous decay by phase modulation fluorometry: mixtures and proteins. Biochemistry. 22:6054–6058.

Beratan, D. N., J. N. Betts, and J. N. Onuchic. 1991. Protein electron transfer rates set by the bridging secondary and tertiary structure. Science. 252:1285–1288.[Abstract/Free Full Text]

Bevington, P. R. 1969. Gradient-expansion algorithm. In Data Reduction and Error Analyses for the Physical Sciences. McGraw-Hill, New York. pp.235–240.

Brooks, B. R., R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus. 1983. CHARMM: a program for molecular energy, minimization, and dynamics calculations. J. Comput. Chem. 4:187–217.

Brown, M. P., and C. Royer. 1997. Fluorescence spectroscopy as a tool to investigate protein interactions. Curr. Opin. Biotechnol. 8:45–49.[Medline]

Callis, P. R., and J. T. Vivian. 2003. Understanding the variable fluorescence quantum yield of tryptophan in proteins using QM-MM simulations. Quenching by charge transfer to the peptide backbone. Chem. Phys. Lett. 369:409–414.

Chang, M. W., J. W. Petrich, D. B. McDonald, and G. R. Fleming. 1983. Nonexponential fluorescence decay of tryptophan and N-terminal tryptophyl peptides. J. Am. Chem. Soc. 105:3819–3824.

Chen, Y., H. T. Yu, and M. D. Barkley. 1996. The peptide bond quenches indole fluorescence. J. Am. Chem. Soc. 118:9271–9278.

Chen, Y., and M. D. Barkley. 1998. Toward understanding tryptophan fluorescence in proteins. Biochemistry. 37:9976–9982.[Medline]

Collins, E. J., J. D. Robertus, M. LoPresti, K. L. Stone, K. R. Williams, P. Wu, K. Hwang, and M. Piatak. 1990. Primary amino acid sequence of alpha-trichosanthin and molecular models for abrin A-chain and alpha-trichosanthin. J. Biol. Chem. 265:8665–8669.[Abstract/Free Full Text]

Dahms, T. E. S., K. J. Willis, and A. G. Szabo. 1995. Conformational heterogeneity of tryptophan in a protein crystal. J. Am. Chem. Soc. 117:2321–2326.

De Beuckeleer, K., G. Volckaert, and Y. Engelborghs. 1999. Time-resolved fluorescence and phosphorescence properties of the individual tryptophan residues of barnase: evidence for protein-protein interactions. Protein Struct. Funct. Gen. 36:42–53.

Delhaise, P., M. Bardiaux, and S. Wodak. 1984. Interactive computer animation of macromolecules. J. Mol. Graph. 2:103–106.[Medline]

Desmet, J., M. De Maeyer, B. Hazes, and I. Lasters. 1992. The dead-end elimination theorem and its use in protein side-chain positioning. Nature. 356:539–542.

De Maeyer, M., J. Desmet, and I. Lasters. 1997. All in one: a highly detailed rotamer library improves both the accuracy and speed in the modelling of sidechains by dead-end elimination. Fold. Des. 2:53–66.[Medline]

De Maeyer, M., J. Desmet, and I. Lasters. 2000. The dead-end elimination theorem: mathematical aspects, implementation, optimizations, evaluation, and performance. Methods Mol. Biol. 143:265–304.[Medline]

Endo, Y., K. Mitsui, M. Motizuki, and K. Tsurugi. 1987. The mechanism of action of ricin and related toxic lectins on eukaryotic ribosomes. The site and the characteristics of modification in 28 S ribosomal RNA caused by the toxins. J. Biol. Chem. 262:5908–5912.[Abstract/Free Full Text]

Goldman, C., P. G. Pascutti, P. Piquini, and A. S. Ito. 1995. On the contribution of electron transfer reactions to quenching of tryptophan fluorescence. J. Chem. Phys. 103:10614–10620.

Gordon, D. B., S. A. Marshall, and S. L. Mayo. 1999. Energy functions for protein design. Curr. Opin. Struct. Biol. 9:509–513.[Medline]

Hao, Q., Y. Zhang, H. Yang, G. Liu, Z. Huang, B. Liu, Q. Yao, and Q. Li. 1995. Fluorescence spectroscopic study of the interaction of adenine and nucleotide with trichosanthin. Biochem. Mol. Biol. Int. 36:889–895.[Medline]

Hennecke, J., A. Sillen, M. Huber-Wunderlich, Y. Engelborghs, and R. Glockshuber. 1997. Quenching of tryptophan fluorescence by the active-site disulfide bridge in the DsbA protein from Escherichia coli. Biochemistry. 36:6391–6400.[Medline]

Hudson, B. S. 1999. An ionization/recombination mechanism for complexity of the fluorescence of tryptophan in proteins. Acc. Chem. Res. 32:297–300.

Kirby, E. P., and R. F. S. Steiner. 1970. The influence of solvent and temperature upon the fluorescence of indole derivatives. J. Phys. Chem. 74:4480–4490.

Laboulais, C., E. Deprez, H. Leh, J. F. Mouscadet, J. C. Brochon, and M. Le Bret. 2001. HIV-1 integrase catalytic core: molecular dynamics and simulated fluorescence decays. Biophys. J. 81:473–489.[Abstract/Free Full Text]

Marcus, R. A., and N. Sutin. 1985. Electron transfer in chemistry and biology. Biochim. Biophys. Acta. 811:265–322.

Mataga, N., Y. Kanda, T. Asahi, H. Miyasaka, T. Okada, and T. Kakitani. 1988a. Mechanisms of the strongly exothermic charge separation reaction in the excited singlet-state-picosecond laser photolysis studies on aromatic hydrocarbon tetracyanoethylene and aromatic hydrocarbon pyromellitic dianhydride systems in polar solutions. Chem. Phys. 127:239–248.

Mataga, N., T. Asahi, Y. Kanda, T. Okada, and T. Kakitani. 1988b. The bell-shaped energy-gap dependence of the charge recombination reaction of geminate radical ion-pairs produced by fluorescence quenching reaction in acetonitrile solution. J. Phys. Chem. 88:5138–5145.

Milner-White, E. J. 1997. The partial charge of the nitrogen atom in peptide bonds. Protein Sci. 6:2477–2482.[Abstract]

More, J. J., and D. C. Sorensen. 1983. Computing a trust region step. SIAM J. Sci. Stat. Comput. 4:553–572.

Paddon-Row, M. N., A. M. Oliver, J. M. Warman, K. J. Smit, M. P. Dehaas, H. Oevering, and J. W. Verhoeven. 1988. Factors affecting charge separation and recombination in photoexcited rigid donor insulator acceptor compounds. J. Phys. Chem. 92:9658–9662.

Ricci, R. W., and J. M. Nesta. 1976. Inter- and intramolecular quenching of indole fluorescence by carbonyl compounds. J. Phys. Chem. 80:974–980.

Parker, C. A., and W. T. Rees. 1960. Corrections of fluorescence spectra and the measurement of fluorescence quantum efficiency. Analyst. 85:587–600.

Sillen, A., and Y. Engelborghs. 1998. The correct use of "average" fluorescence parameters. Photochem. Photobiol. 67:475–486.

Sillen, A., J. Hennecke, D. Roethlisberger, R. Glockshuber, and Y. Engelborghs. 1999. Fluorescence quenching in the DsbA protein from Escherichia coli: complete picture of the excited-state energy pathway and evidence for the reshuffling dynamics of the microstates of tryptophan. Protein Struct. Funct. Gen. 37:253–263.

Sillen, A., J. F. Díaz, and Y. Engelborghs. 2000. A step toward the prediction of the fluorescence lifetimes of tryptophan residues in proteins based on structural and spectral data. Protein Sci. 9:158–169.[Abstract]

Street, A. G., and S. L. Mayo. 1999. Intrinsic ß-sheet propensities result from van der Waals interactions between side chains and the local backbone. Proc. Natl. Acad. Sci. USA. 96:9074–9076.[Abstract/Free Full Text]

Szabo, A. G., and D. M. Rayner. 1980. Fluorescence decay of tryptophan conformers in aqueous solution. J. Am. Chem. Soc. 102:554–563.

Vos, R., and Y. Engelborghs. 1994. A fluorescence study of tryptophan-histidine interaction in the peptide anantin and in solution. Photochem. Photobiol. 60:24–32.[Medline]

Vos, R., Y. Engelborghs, J. Izard, and D. Baty. 1995. Fluorescence study of the three tryptophan residues of the pore-forming domain of Colicin A using multifrequency phase fluorometry. Biochemistry. 35:1734–1743.

Willaert, K., R. Loewenthal, J. Sancho, M. Froeyen, A. R. Fersht, and Y. Engelborghs. 1992. Determination of the excited-state lifetimes of the tryptophan residues in barnase, via multifrequency phase fluorometry of tryptophan mutants. Biochemistry. 31:711–716.[Medline]

Xiong, J. P., Z. X. Xia, and Y. Wang. 1994. Crystal structure of trichosanthin-NADPH complex at 1.7 Å resolution reveals active-site architecture. Nat. Struct. Biol. 1:695–700.[Medline]

Zhang, Y., W. Li, Q. Hao, G. Yu, Q. Li, and Q. Yao. 1993. Tb(III) and Eu(III) as fluorescent probes to investigate the metal-binding sites of trichosanthin. Biochem. Biophys. Res. Commun. 197:407–414.[Medline]




This article has been cited by other articles:


Home page
Biophys. JHome page
R. W. Alston, M. Lasagna, G. R. Grimsley, J. M. Scholtz, G. D. Reinhart, and C. N. Pace
Tryptophan Fluorescence Reveals the Presence of Long-Range Interactions in the Denatured State of Ribonuclease Sa
Biophys. J., March 15, 2008; 94(6): 2288 - 2296.
[Abstract] [Full Text] [PDF]


Home page
Protein Sci.Home page
G. Maglia, A. Jonckheer, M. De Maeyer, J.-M. Frere, and Y. Engelborghs
An unusual red-edge excitation and time-dependent Stokes shift in the single tryptophan mutant protein DD-carboxypeptidase from Streptomyces: The role of dynamics and tryptophan rotamers
Protein Sci., February 1, 2008; 17(2): 352 - 361.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
S. L. C. Moors, M. Hellings, M. De Maeyer, Y. Engelborghs, and A. Ceulemans
Tryptophan Rotamers as Evidenced by X-Ray, Fluorescence Lifetimes, and Molecular Dynamics Modeling
Biophys. J., August 1, 2006; 91(3): 816 - 823.
[Abstract] [Full Text] [PDF]


Home page
Protein Sci.Home page
K. R.F. Somers, P. Kruger, S. Bucikiewicz, M. De Maeyer, Y. Engelborghs, and A. Ceulemans
Protein simulations: The absorption spectrum of barnase point mutants
Protein Sci., July 1, 2004; 13(7): 1823 - 1831.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
C.-P. Pan and M. D. Barkley
Conformational Effects on Tryptophan Fluorescence in Cyclic Hexapeptides
Biophys. J., June 1, 2004; 86(6): 3828 - 3835.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hellings, M.
Right arrow Articles by Engelborghs, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hellings, M.
Right arrow Articles by Engelborghs, Y.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2003 by the Biophysical Society.