Trp42 rotamers report reduced flexibility when the inhibitor acetyl-pepstatin is bound to HIV-1 protease

The Q7K/L331/L631 HIV-1 protease mutant was expressed in Escherichia coli and the effect of binding a substrate-analog inhibitor, acetyl-pepstatin, was investigated by fluorescence spectroscopy and molecular dynamics. The dimeric enzyme has four intrinsic tryptophans, located at positions 6 and 42 in each monomer. Fluorescence spectra and acrylamide quenching experiments show two differently accessible Trp populations in the apoenzyme with k(q1) = 6.85 x 10(9) M(-1) s(-1) and k(q2) = 1.88 x 10(9) M(-1) s(-1), that merge into one in the complex with k(q) = 1.78 x 10(9) M(-1) s(-1). 500 ps trajectory analysis of Trp X1/X2 rotameric interconversions suggest a model to account for the observed Trp fluorescence. In the simulations, Trp6/Trp6B rotameric interconversions do not occur on this timescale for both HIV forms. In the apoenzyme simulations, however, both Trp42s and Trp42Bs are flipping between X1/X2 states; in the complexed form, no such interconverions occur. A detailed investigation of the local Trp environments sampled during the molecular dynamics simulation suggests that one of the apoenzyme Trp42B rotameric interconversions would allow indole-quencher contact, such as with nearby Tyr59. This could account for the short lifetime component. The model thus interprets the experimental data on the basis of the conformational fluctuations of Trp42s alone. It suggests that the rotameric interconversions of these Trps, located relatively far from the active site and at the very start of the flap region, becomes restrained when the apoenzyme binds the inhibitor. The model is thus consistent with associating components of the fluorescence decay in HIV-1 protease to ground state conformational heterogeneity.     

Tryptophan dynamics of the FK506 binding protein: time-resolved fluorescence and simulations.

The FK506-binding protein (FKBP12) is important in the immunosuppressant action of FK506 and rapamycin. We have investigated Trp side chain dynamics in FKBP12, with and without a bound immunosuppressant, by measuring the Trp time-resolved fluorescence anisotropy decay r(t). The r(t) for W59 in aqueous uncomplexed FKBP12 at 20 degrees C is well described by a single exponential with a recovered initial anisotropy, r(eff)o, of 0.192 and an overall rotational correlation time for the protein, phi p, of 4.7 ns; r(eff)o = 0.214 and phi p = 4.2 ns for the FKBP12/FK506 complex. Using an expression for the order parameter squared, namely S2 = r(eff)o/rTo, where rTo is the vitrified steady-state excitation anisotropy, we recovered an S2 of 0.75 for W59 fluorescence in uncomplexed FKBP12 and S2 approximately equal to 1 in the FKBP12/FK506 complex. Results obtained for the FKBP12/rapamycin complex are similar to those found for the FKBP12/FK506 complex. Minimum perturbation mapping simulations were performed on the free and complexed forms of FKBP12 and the results were generally in agreement with the experimental data.     

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.     

Protein dynamics. A time-resolved fluorescence, energetic and molecular dynamics study of ribonuclease T1

Studies using time-resolved fluorescence depolarization were performed on the internal motion of Trp 59 of ribonuclease T1 (EC 3.1.27.3) in the free enzyme, 2'-GMP-enzyme complex and 3'-GMP-enzyme complex. The Trp 59 motion was also studied in the free enzyme using molecular dynamics simulations. Energetic analysis of activation barriers to the Trp 59 motion was performed using both the transition state theory and Kramers' theory. The activation parameters showed a dependence on solvent viscosity indicating the transition state approach in aqueous solution to be inadequate. When taking solvent viscosity contributions into account agreement between the transition state and Kramers' theories was obtained. The results indicate the three enzyme forms to have different conformations with the free enzyme and 3'-GMP-enzyme complex being similar. Comparison of the experimental and theoretical results showed a good agreement on the Trp 59 motion in the free enzyme. Trp 59 appears to vibrate rapidly, with a relaxation time of the order of 1 ps, within free space in the protein matrix and to have a slower motion, with a relaxation time of the order of 100 ps, which is related to breathing of the surrounding protein matrix. Molecular dynamics results indicate high mobility in regions of the enzyme involved in the interaction with the guanine base of the inhibitor or substrate while much lower mobility occurred in residues involved in the catalytic mechanism of ribonuclease T1.     

Determination of rotational correlation times from deconvoluted fluorescence anisotropy decay curves. Demonstration with 6,7-dimethyl-8-ribityllumazine and lumazine protein from Photobacterium leiognathi as fluorescent indicators

The experimental and analytical protocols required for obtaining rotational correlation times of biological macromolecules from fluorescence anisotropy decay measurements are described. As an example, the lumazine protein from Photobacterium leiognathi was used. This stable protein (Mr 21 200) contains the noncovalently bound, natural fluorescent marker 6,7-dimethyl-8-ribityllumazine, which has in the bound state a long fluorescence lifetime (tau = 14 ns). Shortening of the fluorescence lifetime to 2.6 ns at room temperature was achieved by addition of the collisional fluorescence quencher potassium iodide. The shortening of tau had virtually no effect on the rotational correlation time of the lumazine protein (phi = 9.4 ns, 19 degrees C). The ability to measure biexponential anisotropy decay was tested by the addition of Photobacterium luciferase (Mr 80 000), which forms an equilibrium complex with lumazine protein. Under the experimental conditions used (2 degrees C) the biexponential anisotropy decay can best be described with correlation times of 20 and 60 ns, representing the uncomplexed and luciferase-associated lumazine proteins, respectively. The unbound 6,7-dimethyl-8-ribityllumazine itself (tau = 9 ns) was used as a model compound for determining correlation times in the picosecond time range. In the latter case rigorous deconvolution from the excitation profile was required to recover the correlation time, which was shorter (100-200 ps) than the measured laser excitation pulse width (500 ps).     

Picosecond resolution of tyrosine fluorescence and anisotropy decays by 2-GHz frequency-domain fluorometry

We extended the technique of frequency-domain fluorometry to an upper frequency limit of 2000 MHz. This was accomplished by using the harmonic content of a laser pulse train (3.76 MHz, 5 ps) from a synchronously pumped and cavity-dumped dye laser. We used a microchannel plate photomultiplier as the detector to obtain the 2-GHz bandwidth. This new instrument was used to examine tyrosine intensity and anisotropy decays from peptides and proteins. These initial data sets demonstrate that triply exponential tyrosine intensity decays are easily recoverable, even if the mean decay time is less than 1 ns. Importantly, the extended frequency range provides good resolution of rapid and/or multiexponential tyrosine anisotropy decays. Correlation times as short as 15 ps have been recovered for indole, with an uncertainty of +/- 3 ps. We recovered a doubly exponential anisotropy decay of oxytoxin (29 and 454 ps), which probably reflects torsional motions of the phenol ring and overall rotational diffusion, respectively. Also, a 40-ps component was found in the anisotropy decay of bovine pancreatic trypsin inhibitor, which may be due to rapid torsional motions of the tyrosine residues and/or energy transfer among these residues. The rapid component has an amplitude of 0.05, which is about 16% of the total anisotropy. The availability of 2-GHz frequency-domain data extends the measurable time scale for fluorescence to overlap with that of molecular dynamics calculations.     

5-Fluorotryptophan as dual probe for ground-state heterogeneity and excited-state dynamics in apoflavodoxin

The apoflavodoxin protein from Azotobacter vinelandii harboring three tryptophan (Trp) residues, was biosynthetically labeled with 5-fluorotryptophan (5-FTrp). 5-FTrp has the advantage that chemical differences in its microenvironment can be sensitively visualized via ¹⁹F NMR. Moreover, it shows simpler fluorescence decay kinetics. The occurrence of FRET was earlier observed via the fluorescence anisotropy decay of WT apoflavodoxin and the anisotropy decay parameters are in excellent agreement with distances between and relative orientations of all Trp residues. The anisotropy decay in 5-FTrp apoflavodoxin demonstrates that the distances and orientations are identical for this protein. This work demonstrates the added value of replacing Trp by 5-FTrp to study structural features of proteins via ¹⁹F NMR and fluorescence spectroscopy.     

Constrained analysis of fluorescence anisotropy decay:application to experimental protein dynamics

Hydrodynamic properties as well as structural dynamics of proteins can be investigated by the well-established experimental method of fluorescence anisotropy decay. Successful use of this method depends on determination of the correct kinetic model, the extent of cross-correlation between parameters in the fitting function, and differences between the timescales of the depolarizing motions and the fluorophore's fluorescence lifetime. We have tested the utility of an independently measured steady-state anisotropy value as a constraint during data analysis to reduce parameter cross correlation and to increase the timescales over which anisotropy decay parameters can be recovered accurately for two calcium-binding proteins. Mutant rat F102W parvalbumin was used as a model system because its single tryptophan residue exhibits monoexponential fluorescence intensity and anisotropy decay kinetics. Cod parvalbumin, a protein with a single tryptophan residue that exhibits multiexponential fluorescence decay kinetics, was also examined as a more complex model. Anisotropy decays were measured for both proteins as a function of solution viscosity to vary hydrodynamic parameters. The use of the steady-state anisotropy as a constraint significantly improved the precision and accuracy of recovered parameters for both proteins, particularly for viscosities at which the protein's rotational correlation time was much longer than the fluorescence lifetime. Thus, basic hydrodynamic properties of larger biomolecules can now be determined with more precision and accuracy by fluorescence anisotropy decay.     

Ligand-induced changes in the conformational dynamics of a bacterial cytotoxic endonuclease

Knowledge about the conformational dynamics of a protein is key to understanding its biochemical and biophysical properties. In the present work we investigated the dynamic properties of the enzymatic domain of DNase colicins via time-resolved fluorescence and anisotropy decay analysis in combination with steady-state acrylamide quenching experiments. The dynamic properties of the apoenzyme were compared to those of the E9 DNase ligated to the transition metal ion Zn(2+) and the natural inhibitor Im9. We further investigated the contributions of each of the two tryptophans within the E9 DNase (Trp22 and Trp58) using two single-tryptophan mutants (E9 W22F and E9 W58F). Wild-type E9 DNase, E9 W22F, and E9 W58F, as well as Im9, showed multiple lifetime decays. The time-resolved and steady-state fluorescence results indicated that complexation of E9 DNase with Zn(2+) induces compaction of the E9 DNase structure, accompanied by immobilization of Trp22 along with a reduced solvent accessibility for both tryptophans. Im9 binding resulted in immobilization of Trp22 along with a decrease in the longest lifetime component. In contrast, Trp58 experienced less restriction on complexation of E9 DNase with Im9 and showed an increase in the longest lifetime component. Furthermore, the results point out that the Im9-induced changes in the conformational dynamics of E9 DNase are predominant and occur independently of the Zn(2+)-induced conformational effects.     

Backbone dynamics of Tet repressor alpha8intersectionalpha9 loop

A set of single Trp mutants of class B Tet repressor (TetR), in which Trp residues are located from positions 159 to 167, has been engineered to investigate the dynamics of the loop joining the alpha-helices 8 and 9. The fluorescence anisotropy decay of most mutants can be described by the sum of three exponential components. The longest rotational correlation time, 30 ns at 10 degrees C, corresponds to the overall rotation of the protein. The shortest two components, on the subnanosecond and nanosecond time scale, are related to internal motions of the protein. The initial anisotropy, in the 0.16-0.22 range, indicates the existence of an additional ultrafast motion on the picosecond time scale. Examination of physical models for underlying motions indicates that librational motions of the Trp side chain within the rotameric chi(1) x chi(2) potential wells contribute to the picosecond depolarization process, whereas the subnanosecond and nanosecond depolarization processes are related to backbone dynamics. In the absence of inducer, the order parameters of these motions, about 0.90 and 0.80 for most positions, indicate limited flexibility of the loop backbone. Anhydrotetracycline binding to TetR induces an increased mobility of the loop on the nanosecond time scale. This suggests that entropic factors might play a role in the mechanism of allosteric transition.     

Experimentally verifying molecular dynamics simulations through fluorescence anisotropy measurements

The fluorescence anisotropy decay of the single tryptophan residue in phospholipase A2 was studied by use of differential polarized phase fluorometry and computer simulations of protein dynamics. The results enable the verification of a simulated dynamic event by direct experimental measurement on the same time scale. When all hydrogen atoms are modeled explicitly, the simulations agree well with the experimental measurements. However, the measurements contradict simulations in which nonpolar hydrogens are incorporated into "extended" or "united" atoms. These simulations predict an anisotropy decay in excess of measured values and appear to seriously underestimate the electrostatic interactions occurring between water and aromatic side chains. The results support the general validity of studying protein dynamics with the molecular-mechanics approach and illustrate a potentially serious deficiency of simulations which do not explicitly model all hydrogen atoms.     

Exploration of the correlation between solvation dynamics and internal dynamics of a protein

Protein function is intimately related to the dynamics of the protein as well as to the dynamics of the solvent shell around the protein. Although it has been argued extensively that protein dynamics is slaved to solvent dynamics, experimental support for this hypothesis is scanty. In this study, measurements of fluorescence anisotropy decay kinetics have been used to determine the motional dynamics of the fluorophore acrylodan linked to several locations in a small protein barstar in its various structural forms, including the native and unfolded states as well as the acid and protofibril forms. Fluorescence upconversion and streak camera measurements have been used to determine the solvation dynamics around the fluorophore. Both the motional dynamics and solvent dynamics were found to be dependent upon the location of the probe as well as on the structural form of the protein. While the (internal) motional dynamics of the fluorophore occur in the 0.1-3 ns time domain, the observed mean solvent relaxation times are in the range of 20-300 ps. A strong positive correlation between these two dynamical modes was found in spite of the significant difference in their time scales. This observed correlation is a strong indicator of the coupling between solvent dynamics and the dynamics in the protein.     

Structure and dynamics of the gammaM4 transmembrane domain of the acetylcholine receptor in lipid bilayers: insights into receptor assembly and function

A 28-mer peptide (gammaM4) corresponding to the fourth transmembrane segment of the nicotinic acetylcholine receptor (AChR) gamma-subunit, with a single tryptophan residue (Trp6), was reconstituted into lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), loaded with either high or low amounts of cholesterol, i.e., in the conjugated liquid-ordered and liquid-disordered phases, respectively, at room temperature. By making use of the Trp intrinsic fluorescence, both steady-state and time-resolved fluorescence techniques were employed, namely, red-edge excitation shift effect, decay-associated spectra (DAS), and time-resolved anisotropy. The results obtained here, together with previous studies on the same reconstituted peptide, indicate that: (i) Trp6 is strongly anchored in the bilayer with a defined transverse location; (ii) the modifications in the measured DAS are related to the complex result of a self-quenching process on the decay parameters; (iii) the wobbling movement of the indole moiety of Trp6 is fast but severely restricted in amplitude; and, (iv) in the liquid-ordered phase, the bilayer properties and the tilt angle of the peptide enhance peptide-peptide interactions, with the formation of peptide rich patches and possibly some anti-parallel helix-helix aggregates, showing different dynamics from that of the peptide in the liquid-disordered phase where the peptide is randomly distributed.     

Excitation dynamics in Photosystem I from Chlamydomonas reinhardtii. Comparative studies of isolated complexes and whole cells

Identical time-resolved fluorescence measurements with ~ 3.5-ps resolution were performed for three types of PSI preparations from the green alga, Chlamydomonas reinhardtii: isolated PSI cores, isolated PSI–LHCI complexes and PSI–LHCI complexes in whole living cells. Fluorescence decay in these types of PSI preparations has been previously investigated but never under the same experimental conditions. As a result we present consistent picture of excitation dynamics in algal PSI. Temporal evolution of fluorescence spectra can be generally described by three decay components with similar lifetimes in all samples (6–8 ps, 25–30 ps, 166–314 ps). In the PSI cores, the fluorescence decay is dominated by the two fastest components (~ 90%), which can be assigned to excitation energy trapping in the reaction center by reversible primary charge separation. Excitation dynamics in the PSI–LHCI preparations is more complex because of the energy transfer between the LHCI antenna system and the core. The average trapping time of excitations created in the well coupled LHCI antenna system is about 12–15 ps longer than excitations formed in the PSI core antenna. Excitation dynamics in PSI–LHCI complexes in whole living cells is very similar to that observed in isolated complexes. Our data support the view that chlorophylls responsible for the long-wavelength emission are located mostly in LHCI. We also compared in detail our results with the literature data obtained for plant PSI.     

Analysis of correlated motion in antibody combining sites from molecular dynamics simulations

The crystal structures of the NC6.8-antisweet taste ligand complex and the uncomplexed antibody structures display significant differences in the conformations of residues in the combining site. A molecular dynamics method was employed to understand the flexibility and correlated motion of key combining site residues in the uncomplexed antibody. The simulations reveal that residues that show conformational differences between the complex and uncomplexed structures display strong dynamical correlations. Extensive analysis of the dynamics trajectory using time correlation methods is presented.     

Insights into the functional role of protonation states in the HIV-1 protease-BEA369 complex: molecular dynamics simulations and free energy calculations

The molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method combined with molecular dynamics (MD) simulations were used to investigate the functional role of protonation in human immunodeficiency virus type 1 (HIV-1) protease complexed with the inhibitor BEA369. Our results demonstrate that protonation of two aspartic acids (Asp25/Asp25′) has a strong influence on the dynamics behavior of the complex, the binding free energy of BEA369, and inhibitor–residue interactions. Relative binding free energies calculated using the MM-PBSA method show that protonation of Asp25 results in the strongest binding of BEA369 to HIV-1 protease. Inhibitor–residue interactions computed by the theory of free energy decomposition also indicate that protonation of Asp25 has the most favorable effect on binding of BEA369. In addition, hydrogen-bond analysis based on the trajectories of the MD simulations shows that protonation of Asp25 strongly influences the water-mediated link of a conserved water molecule, Wat301. We expect that the results of this study will contribute significantly to binding calculations for BEA369, and to the design of high affinity inhibitors.     

Substrate-dependent conformational dynamics of the Escherichia coli membrane insertase YidC

The binding of Pf3 coat protein to the membrane insertase YidC from Escherichia coli induces a conformational change in the tertiary structure of the insertase, resulting in a quenching of the intrinsic tryptophan (Trp) fluorescence. Tryptophan mutants of YidC were generated to examine such conformational movements in detail with time-resolved and steady-state fluorescence spectroscopy. Ten of the 11 Trp residues within YidC were substituted to phenylalanines generating single Trp mutants either at position 354, 454, or 508. In addition, a double mutant with Trp residues at 332 and 334 was studied. Purified YidC mutants were reconstituted into DOPC/DOPG vesicles and titrated with a Trp-free mutant of Pf3 coat, enabling a detailed conformational study of the periplasmic P1, P2, and P3 domains of YidC before and after binding of substrate. Time-resolved fluorescence anisotropy revealed that the mobility of the residues W332/W334 and W508 was considerably increased after binding of Pf3 coat to the insertase. Furthermore, analysis of the fluorescence emission spectra and the decay times showed that all Trp residues are embedded in an equivalent environment that is a membrane/water interface.     

Site-specific tryptophan dynamics in class A amphipathic helical peptides at a phospholipid bilayer interface

The amphipathic helix plays a key role in many membrane-associating peptides and proteins. The dynamics of helices on membrane surfaces might be of importance to their function. The fluorescence anisotropy decay of tryptophan is a sensitive indicator of local, segmental, and global dynamics within a peptide or protein. We describe the use of frequency domain dynamic depolarization measurements to determine the site-specific tryptophan dynamics of single tryptophan amphipathic peptides bound to a phospholipid surface. The five 18-residue peptides studied are based on a class A amphipathic peptide that is known to associate at the interface of phospholipid bilayers. The peptides contain a single tryptophan located at positions 2, 3, 7, 12, or 14 in the sequence. Association of the peptides with egg phosphatidylcholine vesicles results in complex behavior of both the tryptophan intensity decay and the anisotropy decay. The anisotropy decays were biphasic and were fitted to an associated model where each lifetime component in the intensity decay is associated with a particular rotational correlation time from the anisotropy decay. In contrast, an unassociated model where all components of the intensity decay share common rotational modes was unable to provide an adequate fit to the data. Two correlation times were resolved from the associated analysis: one whose contribution to the anisotropy decay was dependent on the exposure of the tryptophan to the aqueous phase, and the other whose contribution reflected the position of the tryptophan in the sequence. The results are compared with existing x-ray structural data and molecular dynamics simulations of membrane-incorporated peptides.     

Analysis of internal motion of single tryptophan in Streptomyces subtilisin inhibitor from its picosecond time-resolved fluorescence.

A mode of internal motion of single tryptophan, Trp 86, of Streptomyces subtilisin inhibitor, was analyzed from its time-resolved fluorescence. The intensity and anisotropy decays of Trp 86 were measured in the picosecond range. These decays were analyzed with theoretical expressions derived assuming that the indole ring of tryptophan as an asymmetric rotor rotates around covalent bonds connecting indole with the peptide chain and an effective quencher of fluorescence of Trp 86 is the nearby SS bond of Cys 35-Cys 50. First, the intensity decays at 6 degrees, 20 degrees, and 40 degrees C were analyzed, and then the both decays of the intensity and anisotropy at 20 degrees C were simultaneously simulated with common parameters. Constants concerning geometrical structures of the protein used for the analysis were obtained from x-ray crystallographic data. Best fit between the observed and calculated decay curves was obtained by a nonlinear least squares method by adjusting a quenching constant averaged over the rotational angles, koq height of the potential energy, p, and three of six diffusion coefficients, Dxx, Dyy, Dzz, Dxy, Dyz, and Dzx, as variable parameters. The obtained results revealed that the internal motion of the indole ring became faster, the quenching rate of the fluorescence of Trp 86 was enhanced and the height of potential energy became lower at higher temperatures, and suggested that Trp 86 was wobbling around the long axis of the indole ring in the protein.     

Molecular dynamics of tryptophan in ribonuclease-T1. I. Simulation strategies and fluorescence anisotropy decay.

Molecular dynamics simulations of Ribonuclease-T1 (RNAse-T1) were performed using x-ray crystal coordinates for the enzyme and various simulation strategies. From each of the simulations, a predicted fluorescence anisotropy decay for the single-tryptophan residue was derived and compared with experimental values for the limiting anisotropy of this protein. Simulations conducted in vacuo demonstrated large displacements among some of the residues adjacent to the tryptophan side chain. As a consequence, the ring system rotates relatively unhindered through an angle far in excess of that implied by experimental data. In contrast, the explicit simulation of solvent within a stochastic boundary led to excellent agreement between simulation and experiment. In the case of RNAse-T1, the experimentally-determined limiting anisotropy is useful as a criterion of simulation accuracy in the vicinity of the tryptophan side chain.