Tryptophanyl fluorescence lifetime distribution of hyperthermophilic beta-glycosidase from molecular dynamics simulation: a comparison with the experimental data

A molecular dynamics simulation approach has been utilized to understand the unusual fluorescence emission decay observed for beta-glycosidase from the hyperthermophilic bacterium Solfolobus sulfotaricus (Sbeta gly), a tetrameric enzyme containing 17 tryptophanyl residues for each subunit. The tryptophanyl emission decay of Sbeta gly results from a bimodal distribution of fluorescence lifetimes with a short-lived component centered at 2.5 ns and a long-lived one at 7.4 ns (Bismuto E, Nucci R, Rossi M, Irace G, 1999, Proteins 27:71-79). From the examination of the trajectories of the side chains capable of causing intramolecular quenching for each tryptophan microenvironment and using a modified Stern-Volmer model for the emission quenching processes, we calculated the fluorescence lifetime for each tryptophanyl residue of Sbeta gly at two different temperatures, i.e., 300 and 365 K. The highest temperature was chosen because in this condition Sbeta gly evidences a maximum in its catalytic activity and is stable for a very long time. The calculated lifetime distributions overlap those experimentally determined. Moreover, the majority of trytptophanyl residues having longer lifetimes correspond to those originally identified by inspection of the crystallographic structure. The tryptophanyl lifetimes appear to be a complex function of several variables, such as microenvironment viscosity, solvent accessibility, the chemical structure of quencher side chains, and side-chain dynamics. The lifetime calculation by MD simulation can be used to validate a predicted structure by comparing the theoretical data with the experimental fluorescence decay results.     

The beta-glycosidase from the hyperthermophilic archaeon Sulfolobus solfataricus: enzyme activity and conformational dynamics at temperatures above 100 degrees C.

Enzymes from thermophilic organisms are stable and active at temperatures which rapidly denature mesophilic proteins. However, there is not yet a complete understanding of the structural basis of their thermostability and thermoactivity since for each protein there seems to exist special networks of interactions that make it stable under the desired conditions. Here we have investigated the activity and conformational dynamics above 100 degrees C of the beta-glycosidase isolated from the hyperthermophilic archaeon Sulfolobus solfataricus. This has been made possible using a special stainless steel optical pressure cell which allowed us to perform enzyme assays and fluorescence measurements up to 160 degrees C without boiling the sample. The beta-glycosidase from S. solfataricus showed maximal activity at 125 degrees C. The time-resolved fluorescence studies showed that the intrinsic tryptophanyl fluorescence emission of the protein was represented by a bimodal distribution with Lorential shape and that temperature strongly affected the protein conformational dynamics. Remarkably, the tryptophan emission reveals that the indolic residues remain shielded from the solvent even at 125 degrees C, as shown by shielding from quenching and restricted tryptophan solubility. The relationship between enzyme activity and protein structural dynamics is discussed.     

The thermophilic esterase from Archaeoglobus fulgidus: structure and conformational dynamics at high temperature

The esterase from the hyperthermophilic archaeon Archaeoglobus fulgidus is a monomeric protein with a molecular weight of about 35.5 kDa. The enzyme is barely active at room temperature, displaying the maximal enzyme activity at about 80 degrees C. We have investigated the effect of the temperature on the protein structure by Fourier-transform infrared spectroscopy. The data show that between 20 degrees C and 60 degrees C a small but significant decrease of the beta-sheet bands occurred, indicating a partial loss of beta-sheets. This finding may be surprising for a thermophilic protein and suggests the presence of a temperature-sensitive beta-sheet. The increase in temperature from 60 degrees C to 98 degrees C induced a decrease of alpha-helix and beta-sheet bands which, however, are still easily detected at 98 degrees C indicating that at this temperature some secondary structure elements of the protein remain intact. The conformational dynamics of the esterase were investigated by frequency-domain fluorometry and anisotropy decays. The fluorescence studies showed that the intrinsic tryptophanyl fluorescence of the protein was well represented by the three-exponential model, and that the temperature affected the protein conformational dynamics. Remarkably, the tryptophanyl fluorescence emission reveals that the indolic residues remained shielded from the solvent up to 80 degrees C, as shown from the emission spectra and by acrylamide quenching experiments. The relationship between enzyme activity and protein structure is discussed.     

Structural and dynamic aspects of beta-glycosidase from mesophilic and thermophilic bacteria by multitryptophanyl emission decay studies

The tryptophanyl emission decay of beta-glycosidase from the extremophilic archaeon Sulfolobus solfataricus (Sbetagly) has been investigated by frequency domain fluorometry. The data were analyzed in terms of sum of discrete lifetimes as well as in terms of quasi- continuous lifetime distributions of different shape. At neutral pH the emission decay is characterized by two components: a long-lived component, centered at 7.4 ns, and a short one at 2.7 ns, irrespective of the decay scheme used for the interpretation of the experimental results. The effects of an irreversible inhibitor, that is, cyclophellitol, and that of a powerful denaturant such as guanidinium hydrochloride on the dynamics of Sbetagly has been investigated by observing the changes induced in the two components of the tryptophanyl emission decay. The addition of cyclophellitol to native Sbetagly reduces the contribution of the short-lived component but does not affect the long-lived one. Increasing concentrations of guanidinium hydrochloride differently affect the contributions of the two emission components. Higher concentrations were required to unfold the molecular regions containing the long-lived indolic fluorophores. These results indicate that the long-lived contribution arises from tryptophanyl residues deeply clustered in the interior of the protein matrix, whereas the short-lived one includes residues located in less rigid and more solvent accessible regions, some of which might be located in functionally important parts of protein. The knowledge of the crystallographic structure of Sbetagly allowed us to evaluate some average parameters for each tryptophanyl microenvironment in the Sbetagly such as hydrophobicity, structural flexibility, and ability of side chains to act as fluorescence quenchers. These results permitted to divide the tryptophanyl fluorescence of Sbetagly in the contribution of two emitting groups: one consisting of eight closely clustered tryptophans, that is, Trp 33, 36, 60, 84, 151 174, 425, and 433, responsible for the long-lived emission component and the other one, composed of nine tryptophans nearer to the subunit surface, that is, Trp 12, 156, 192, 287, 288, 316, 361, 376, 455, associable to the short-lived emission component. Finally, the examination of the tryptophanyl emission decay of the mesophilic beta-galactosidase from Escherichia coli (Cbetagal) and the Arrhenius analysis of its dependence on temperature indicated that the tryptophanyl environments of the mesophilic enzyme are rather homogeneous in consequence of a larger protein dynamics.     

The effect of molecular confinement on the conformational dynamics of the native and partly folded state of apomyoglobin.

Inclusion in agarose gel significantly affects the conformational dynamics of native and acidic partly folded states of tuna apomyoglobin, a single tryptophan containing protein, as documented by frequency domain fluorometry investigations. The heterogeneity of the tryptophanyl emission decay increases on gel inclusion compared to that observed for free-in-solvent protein at both neutral and acidic pH, thus suggesting that the interconversion rate among conformational substates is somewhat reduced. The observation that this effect is much more pronounced for the partly folded state suggests that confined environments such as those existing in the living cells might favor the sequential folding process avoiding that structured intermediates rapidly convert into less structured ones.     

Perturbation of conformational dynamics,enzymatic activity,and thermostability of β-glycosidase from archaeon Sulfolobus solfataricus by pH and sodium dodecyl sulfate detergent

The conformational dynamics of β-glycosidase from Sulfolobus solfataricus was investigated by following the emission decay arising from the large number of tryptophanyl residues that are homogeneously dispersed in the primary structure. The fluorescence emission is characterized by a bimodal lifetime distribution, suggesting that the enzyme structure contains rigid and flexible regions, properly located in the macromolecule. The enzyme activity and thermostability appear to be related to the dynamic properties of these regions as evidenced by perturbation studies of the enzyme structure at alkaline pH and by addition of detergents such as SDS. The pH increase affects the protein dynamics with a remarkable loss of thermal stability and activity; these changes occur without any significant variation in the secondary structure as revealed by far-UV dichroic measurements. In the presence of 0.02% (w/v) SDS at alkaline pH, the enzymatic activity and thermostability are recovered. Under these conditions, the conformational dynamics appear to be similar to that evidenced at neutral pH. Further increases in SDS concentration, at alkaline pH, render the activity and thermostability of β-glycosidase similar to those observed in the absence of detergent. Proteins 27:71–79 © 1997 Wiley-Liss, Inc.     

A pulse fluorometry study of lipoamide dehydrogenase. Evidence for non-equivalent FAD centers.

The time dependence of the fluorescence of tryptophanyl and flavin residues in lipoamide dehydrogenase has been investigated with single-photon decay spectroscopy. When the two FAD molecules in the enzyme were directly excited the decay could only be analyzed in a sum of two exponentials with equal amplitudes. This phenomenon was observed at 4 degrees C (tau-1 = 0.8 ns, tau-2 = 4.7 ns) and at 20 degrees C (tau-1 = 0.8 ns, tau-2 = 3.4 ns) irrespective of the emission and excitation wavelengths. This result reveals a difference in the nature of the two FAD centers. By excitation at 290 nm the fluorescence decay curves of tryptophan and FAD were obtained. The decays are analyzed in terms of energy transfer from tryptophanyl to flavin residues. The results, which are in good agreement with those obtained previously with static fluorescence methods, show that one of the two tryptophanyl residues within the subunit transfers its excitation energy to the flavin located at a distance of 1.5 nm.     

A comparative molecular dynamics study of thermophilic and mesophilic ribonuclease HI enzymes

We studied a pair of homologous thermophilic and mesophilic ribonuclease HI enzymes by molecular dynamics simulations. Each protein was subjected to three 5 ns simulations in explicit water at both 310 K and 340 K. The thermophilic enzyme showed larger overall positional fluctuations at both temperatures, while only the mesophilic enzyme at the higher temperature showed significant instability. When the temperature is changed, the relative flexibility of different local segments on the two proteins changed differently. Principal component analysis showed that the simulations of the two proteins explored largely overlapping regions in the conformational space. However, at 340 K, the collective structure variations of the thermophilic protein are different from those of the mesophilic protein. Our results, although not in accordance with the view that hyperthermostability of proteins may originate from their conformational rigidity, are consistent with several recent experimental and simulation studies which showed that thermophilic proteins may be conformationally more flexible than their mesophilic counterparts. The decorrelation between conformational rigidity and hyperthermostability may be attributed to the temperature dependence and long range nature of electrostatic interactions that play more important roles in the structural stability of thermophilic proteins.     

Active-site residues move independently from the rest of the protein in a 200 ns molecular dynamics simulation of cytochrome P450 CYP119

The conformational dynamics of cytochrome P450 enzymes are critical to their catalytic activity. In this study, the correlated motion between residues in a 200 ns molecular dynamics trajectory of the thermophilic CYP119 was analyzed to parse out conformational relationships. Residues that are structurally related, for example residues within a helix, generally have highly correlated motion. In addition, clusters of non-adjacent residues that show correlated motion (“hot spots”) are seen in various regions, including at the base of the F and G helices that make up the most dynamic region of the enzyme. A modified k-means algorithm that clusters residues based on their correlated motion indicates that functionally related residues are in the same cluster (e.g., the catalytic threonines and the heme). Tightly coupled clusters form a solvent-exposed “shell” around the enzyme, whereas less coupling between clusters is seen in regions that are critical to ligand interactions, redox partner interactions, and catalysis. Most notably, we find that residues in the active site move independently from the rest of the enzyme, effectively insulating the catalytic machinery from other regions of the protein.     

Relation between stability, dynamics and enzyme activity in 3-phosphoglycerate kinases from yeast and Thermus thermophilus.

3-Phosphoglycerate kinases from yeast and the extreme thermophilic bacterium Thermus thermophilus HB8 have been used as models for investigating the relationship between stability, dynamics and activity. It was found that while at a given temperature the thermophilic protein is more stable, its conformational dynamics as measured by the ability of acrylamide to quench the fluorescence of a buried tryptophan as well as its specific activity, are both lower than for the mesophilic protein. As the temperature is increased, the thermodynamic stability of the thermophilic protein approaches that of the mesophilic protein at its working temperature. Its conformational dynamics and specific activity however were both shown to increase, until at the physiologically operational temperature, they become similar to those of the mesophilic enzyme at its operational temperature. These results confirm the proposal that a direct relationship and balance holds between thermodynamic stability, dynamics and specific activity in globular proteins. They demonstrate also the constraining effect of increased stability upon conformational dynamics and enzyme activity.     

Role of loop dynamics in thermal stability of mesophilic and thermophilic adenylosuccinate synthetase: a molecular dynamics and normal mode analysis study

Enzymes from thermophiles are poorly active at temperatures at which their mesophilic homologs exhibit high activity and attain corresponding active states at high temperatures. In this study, comparative molecular dynamics (MD) simulations, supplemented by normal mode analysis, have been performed on an enzyme Adenylosuccinate synthetase (AdSS) from E. coli (mesophilic) and P. horikoshii (thermophilic) systems to understand the effects of loop dynamics on thermal stability of AdSS. In mesophilic AdSS, both ligand binding and catalysis are facilitated through the coordinated movement of five loops on the protein. The simulation results suggest that thermophilic P. horikoshii preserves structure and catalytic function at high temperatures by using the movement of only a subset of loops (two out of five) for ligand binding and catalysis unlike its mesophilic counterpart in E. coli. The pre-arrangement of the catalytic residues in P. horikoshii is well-preserved and salt bridges remain stable at high temperature (363K). The simulations suggest a general mechanism (including pre-arrangement of catalytic residues, increased polar residue content, stable salt bridges, increased rigidity, and fewer loop movements) used by thermophilic enzymes to preserve structure and be catalytically active at elevated temperatures.     

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included α-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein. Figure Collective motions in C^α atoms of the active site of cold-active xylanase Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Frequency domain fluorescence studies of yeast phosphoglycerate kinase and its ternary complex

A frequency domain fluorescence study of yeast phosphoglycerate kinase has been performed to observe the effect of substrates on the structure and dynamics of the enzyme. At 20 degrees C and pH 7.2, a biexponential decay is observed for tryptophanyl emission. The short fluorescence lifetime (0.4 ns) component is associated with a spectrum having a 329-nm maximum and a 18.4-kJ/mol activation energy, Ea, for thermal quenching. The long-lifetime (3.5 ns) component has a 338-nm maximum and an Ea of only 7.9 kJ/mol. Tentatively we assign the short and long-lifetime components to Trp-333 and Trp-308. Binding of the substrates ATP and 3-phosphoglycerate leads to a significant increase in the fluorescence lifetime, the red shift of the emission spectrum and in the decrease in the Ea for both components. Acrylamide-quenching studies indicate that the two tryptophan residues have about the same degree of kinetic exposure to the quencher and that the binding of the substrates causes a very slight change in the quenching pattern. These fluorescence studies indicate that the binding of the substrates to phosphoglycerate kinase may influence the conformational dynamics around the two tryptophan residues located on one of the protein's domains.     

Linkage between dynamics and catalysis in a thermophilic-mesophilic enzyme pair

A fundamental question is how enzymes can accelerate chemical reactions. Catalysis is not only defined by actual chemical steps, but also by enzyme structure and dynamics. To investigate the role of protein dynamics in enzymatic turnover, we measured residue-specific protein dynamics in hyperthermophilic and mesophilic homologs of adenylate kinase during catalysis. A dynamic process, the opening of the nucleotide-binding lids, was found to be rate-limiting for both enzymes as measured by NMR relaxation. Moreover, we found that the reduced catalytic activity of the hyperthermophilic enzyme at ambient temperatures is caused solely by a slower lid-opening rate. This comparative and quantitative study of activity, structure and dynamics revealed a close link between protein dynamics and catalytic turnover.     

Multiple time scale backbone dynamics of homologous thermophilic and mesophilic ribonuclease HI enzymes

Backbone conformational fluctuations on multiple time scales in a cysteine-free Thermus thermophilus ribonuclease HI mutant (ttRNH(*)) are quantified using (15)N nuclear magnetic spin relaxation. Laboratory-frame relaxation data acquired at 310 K and at static magnetic field strengths of 11.7, 14.1 and 18.8 T are analysed using reduced spectral density mapping and model-free approaches. Chemical exchange line broadening is characterized using Hahn-echo transverse and multiple quantum relaxation data acquired over a temperature range of 290-320 K and at a static magnetic field strength of 14.1 T. Results for ttRNH(*) are compared to previously published data for a mesophilic homologue, Escherichia coli ribonuclease HI (ecRNH). Intramolecular conformational fluctuations on the picosecond-to-nanosecond time scale generally are similar for ttRNH(*) and ecRNH. beta-Strands 3 and 5 and the glycine-rich region are more rigid while the substrate-binding handle region and C-terminal tail are more flexible in ttRNH(*) than in ecRNH. Rigidity in the two beta-strands and the glycine-rich region, located along the periphery of the central beta-sheet, may be associated with the increased thermodynamic stability of the thermophilic enzyme. Chemical exchange line broadening, reflecting microsecond-to-millisecond time scale conformational changes, is more pronounced in ttRNH(*) than in ecRNH, particularly for residues in the handle and surrounding the catalytic site. The temperature dependence of chemical exchange show an increase of approximately 15 kJ/mol in the apparent activation energies for ttRNH(*) residues in the handle compared to ecRNH. Increased activation barriers, coupled with motion between alpha-helices B and C not present in ecRNH, may be associated with the reduced catalytic activity of the thermophilic enzyme at 310 K.     

Flexibility and enzymatic cold-adaptation: a comparative molecular dynamics investigation of the elastase family

Molecular dynamics simulations of representative mesophilic and psycrophilic elastases have been carried out at different temperatures to explore the molecular basis of cold adaptation inside a specific enzymatic family. The molecular dynamics trajectories have been compared and analyzed in terms of secondary structure, molecular flexibility, intramolecular and protein-solvent interactions, unravelling molecular features relevant to rationalize the efficient catalytic activity of psychrophilic elastases at low temperature. The comparative molecular dynamics investigation reveals that modulation of the number of protein-solvent interactions is not the evolutionary strategy followed by the psycrophilic elastase to enhance catalytic activity at low temperature. In addition, flexibility and solvent accessibility of the residues forming the catalytic triad and the specificity pocket are comparable in the cold- and warm-adapted enzymes. Instead, loop regions with different amino acid composition in the two enzymes, and clustered around the active site or the specificity pocket, are characterized by enhanced flexibility in the cold-adapted enzyme. Remarkably, the psycrophilic elastase is characterized by reduced flexibility, when compared to the mesophilic counterpart, in some scattered regions distant from the functional sites, in agreement with hypothesis suggesting that local rigidity in regions far from functional sites can be beneficial for the catalytic activity of psychrophilic enzymes.     

Fluorescence investigation on conformational state of rabbit muscle aldolase interacting with phosphatidylinositol liposomes

Evidence of conformational changes in rabbit muscle aldolase upon binding to phosphatidylinositol liposomes and the effect of the interaction on the thermal conformational transition are reported. Interaction with phosphatidylinositol liposomes significantly decreases the aldolase tryptophanyl fluorescence and shifts the maximum wavelength to higher values. The dynamic quenching constant for the aldolase fluorescence quenching by acrylamide in the presence of liposomes is much higher than that for unmodified enzyme; this signifies an increase in accessibility of some tryptophanyl residues to small polar molecules. Indirect interaction between single phospholipid molecules, small micelles or any soluble impurities able to penetrate into the protein molecule interior does not seem to be involved in the conformational rearrangement. Native and liposome-interaction-induced conformational states reveal different temperature dependences of the tryptophan residues exposure. The implications of the modification of the conformational state of the enzyme for its function in vivo are discussed.     

S-adenosylhomocysteine hydrolase from the archaeon Pyrococcus furiosus: biochemical characterization and analysis of protein structure by comparative molecular modeling

S-adenosylhomocysteine hydrolase (AdoHcyHD) is an ubiquitous enzyme that catalyzes the breakdown of S-adenosylhomocysteine, a powerful inhibitor of most transmethylation reactions, to adenosine and L-homocysteine. AdoHcyHD from the hyperthermophilic archaeon Pyrococcus furiosus (PfAdoHcyHD) was cloned, expressed in Escherichia coli, and purified. The enzyme is thermoactive with an optimum temperature of 95 degrees C, and thermostable retaining 100% residual activity after 1 h at 90 degrees C and showing an apparent melting temperature of 98 degrees C. The enzyme is a homotetramer of 190 kDa and contains four cysteine residues per subunit. Thiol groups are not involved in the catalytic process whereas disulfide bond(s) could be present since incubation with 0.8 M dithiothreitol reduces enzyme activity. Multiple sequence alignment of hyperthermophilic AdoHcyHD reveals the presence of two cysteine residues in the N-terminus of the enzyme conserved only in members of Pyrococcus species, and shows that hyperthermophilic AdoHcyHD lack eight C-terminal residues, thought to be important for structural and functional properties of the eukaryotic enzyme. The homology-modeled structure of PfAdoHcyHD shows that Trp220, Tyr181, Tyr184, and Leu185 of each subunit and Ile244 from a different subunit form a network of hydrophobic and aromatic interactions in the central channel formed at the subunits interface. These contacts partially replace the interactions of the C-terminal tail of the eukaryotic enzyme required for tetramer stability. Moreover, Cys221 and Lys245 substitute for Thr430 and Lys426, respectively, of the human enzyme in NAD-binding. Interestingly, all these residues are fairly well conserved in hyperthermophilic AdoHcyHDs but not in mesophilic ones, thus suggesting a common adaptation mechanism at high temperatures.     

Heterogeneity in the structural dynamics of Sulfolobus solfataricus beta-glycosidase revealed by electron paramagnetic resonance and frequency domain fluorometry

The local and global dynamics of the Sulfolobus solfataricus beta-glycosidase were studied by electron spin resonance and time-resolved fluorescence techniques. For electron paramagnetic resonance (EPR) investigations, the protein was covalently modified by the maleimido nitroxide spin label, which is specific for cysteine -SH groups, at position 344 and at position 101, where Ser-101 was changed into a cysteine by site-directed mutagenesis. The greater reactivity of exposed Cys-101 suggested the exclusive modification of this amino acid compared with Cys-344. The labeled proteins underwent temperature perturbation in the range 290-335 K and the values of the spin-label rotation correlation frequencies (nu(c)) ranged from 6 x 10(7) to 2 x 10(8) sec(-1) for the protein labeled at position C344 and from 5.62 x 10(7) to 1.10 x 10(8) sec(-1) for the protein labeled at C101. These rotation correlation values are related to the local dynamic characteristics of the protein matrix. The temperature dependence of rotation correlation frequencies expressed in terms of Arrhenius coordinates (log (nu(c)) vs. 1/T) for the protein labeled at C344 exhibited a linear dependence but with a change in the slope at 311 K. For the protein labeled at C101, no change in the slope was observed at the same temperature. General dynamic information was deduced from the analysis of the fluorescence emission decay of the tryptophanyl residues that are present in each region of the protein structure. Fluorescence data analysis highlighted a bimodal distribution of fluorescence lifetimes arising from the contribution of two emitting groups: one consisting of closely clustered tryptophans responsible for the long-lived emission component (7.1 nsec) and the other composed of tryptophans nearer to the protein surface, which can be associated to the short-lived component (2.5 nsec). The temperature dependence of lifetime distribution parameters linked to the long-lived and short-lived components, expressed in Arrhenius coordinates, showed two different points in which the change in the slope occurred (i.e., 328 K and 338 K, respectively). The Arrhenius analysis of data provided the activation energy relative to the conformational changes characterizing the local and global movements running through the protein matrix.     

HIV-1 integrase catalytic core: molecular dynamics and simulated fluorescence decays.

Two molecular dynamics simulations have been carried out on the HIV-1 integrase catalytic core starting from fully determined crystal structures. During the first one, performed in the absence of divalent cation (6-ns long), the catalytic core took on two main conformations. The conformational transition occurs at approximately 3.4 ns. In contrast, during the second one, in the presence of Mg(2+) (4-ns long), there were no such changes. The molecular dynamics simulations were used to compute the fluorescence intensity decays emitted by the four tryptophan residues considered as the only chromophores. The decay was computed by following, frame by frame, the amount of chromophores that remained excited at a certain time after light absorption. The simulation took into account the quenching through electron transfer to the peptide bond and the fluorescence resonance energy transfer between the chromophores. The fit to the experimental intensity decays obtained at 5 degrees C and at 30 degrees C is very good. The fluorescence anisotropy decays were also simulated. Interestingly, the fit to the experimental anisotropy decay was excellent at 5 degrees C and rather poor at 30 degrees C. Various hypotheses such as dimerization and abnormal increase of uncorrelated internal motions are discussed.