Short range photoinduced electron transfer in proteins: QM-MM simulations of tryptophan and flavin fluorescence quenching in proteins

Factors deciding the assembly and thermostability of the DmrB cage

Dihydromethanopterin reductase (DmrB), is a naturally occurring cage protein found in various archaeal and a few bacterial species. It exists as 24mer with cubic geometry where 8 trimeric subunits are present at the corners of each cube. Each trimer is made up of three monomeric units and six FMN, where two molecules of FMN are present at the interface of each monomer. DmrB is involved in the conversion of dihydromethanopterin to tetrahydromethanopterin using FMN as a redox equivalent. In the present study, we have used spectroscopic and biochemical techniques along with complementary bio-informatic work to understand the assembly principles of the DmrB. Our results show a concentration dependant self-assembly of DmrB which is mediated by ionic interactions. The co-factor FMN stabilizes and preserves the secondary and quaternary structure of DmrB against thermal insult, indicating that the higher order assembly of DmrB is very thermostable. Our work provides an interesting piece of information regarding the role of the co-factors in the thermostability of these classes of cage proteins. The understanding of the assembly and disassembly of this thermostable cage would enable the downstream usage of this system in various nano-biotechnological applications.

Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding

Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.

Evidence of a Thermodynamic Ramp for Hole Hopping to Protect a Redox Enzyme from Oxidative Damage

Redox proteins and enzymes are at risk of irreversible oxidative damage from highly oxidizing intermediates generated in the active site in the case of unsuccessful functional reaction. Chains of tyrosine and/or tryptophan residues have been recently proposed to provide protection to the active site and the whole protein by delivering oxidizing equivalents (holes) out of the protein via a multistep hopping mechanism. In the present work we use a hybrid quantum/classical theoretical-computational methodology based on the perturbed matrix method and on molecular dynamics simulations to calculate the oxidation potential difference along a chain of tyrosine and tryptophan residues in a human redox enzyme of major importance, a superoxide dismutase, which acts as antioxidant defense. We show that the hole hopping is thermodynamically favored along such a chain and that the hopping propensity is strongly affected by the protein environment and in particular by the active site and its second coordination sphere.

Extending the perturbed matrix method beyond the dipolar approximation: comparison of different levels of theory

Some years ago we developed a theoretical-computational hybrid quantum/classical methodology, the Perturbed Matrix Method (PMM), to be used in conjunction with molecular dynamics simulations for the investigation of chemical processes in complex systems, that proved to be a valuable tool for the simulation of relevant experimental observables, e.g., spectroscopic signals, reduction potentials, kinetic constants. In typical PMM calculations the quantum sub-part of the system, the quantum centre, is embedded into an external perturbing field providing a perturbation operator explicitly calculated up to the dipolar terms. In this paper we further develop the PMM approach, beyond the dipolar terms in the perturbation operator expansion, by including explicitly the quadrupolar terms and/or by expanding the perturbation operator on each atom of the quantum centre. These different levels of the perturbation operator expansion, providing different levels of theory, have been tested by calculating three different spectroscopic observables: the spectral signal of liquid water and aqueous benzene due to the lowest energy electronic excitation and the infrared amide I band of aqueous trans-N-methylacetamide. All the systems tested show that, even though the previous PMM level of theory is already capable of reproducing the main features of the spectral signal, the higher levels of theory improve the quantitative reproduction of the spectral details.

Rapid and Programmable Protein Mutagenesis Using Plasmid Recombineering

Comprehensive and programmable protein mutagenesis is critical for understanding structure-function relationships and improving protein function. There is thus a need for robust and unbiased molecular biological approaches for the construction of the requisite comprehensive protein libraries. Here we demonstrate that plasmid recombineering is a simple and robust in vivo method for the generation of protein mutants for both comprehensive library generation as well as programmable targeting of sequence space. Using the fluorescent protein iLOV as a model target, we build a complete mutagenesis library and find it to be specific and comprehensive, detecting 99.8% of our intended mutations. We then develop a thermostability screen and utilize our comprehensive mutation data to rapidly construct a targeted and multiplexed library that identifies significantly improved variants, thus demonstrating rapid protein engineering in a simple protocol.

Alternative Electron-Transfer Channels Ensure Ultrafast Deactivation of Light-Induced Excited States in Riboflavin Binding Protein

Flavoproteins, containing flavin chromophores, are enzymes capable of transferring electrons at very high speeds. The ultrafast photoinduced electron-transfer (ET) kinetics of riboflavin binding protein to the excited riboflavin was studied by femtosecond spectroscopy and found to occur within a few hundred femtoseconds [ Zhong and Zewail, Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 11867-11872 ]. This ultrafast kinetics was attributed to the presence of two aromatic rings that could transfer the electron to riboflavin: the side chains of tryptophan 156 and tyrosine 75. However, the underlying ET mechanism remained unclear. Here, using a hybrid quantum mechanical-molecular dynamics approach, we perform ET dynamics simulations taking into account the motion of the protein and the solvent upon ET. This approach reveals that ET occurs via a major reaction channel involving tyrosine 75 (83%) and a minor one involving tryptophan 156 (17%). We also show that the protein environment is designed to ensure the fast quenching of the riboflavin excited state.

Excited state dynamics can be used to probe donor-acceptor distances for H-tunneling reactions catalyzed by flavoproteins

In enzyme systems where fast motions are thought to contribute to H-transfer efficiency, the distance between hydrogen donor and acceptor is a very important factor. Sub-ångstrom changes in donor-acceptor distance can have a large effect on the rate of reaction, so a sensitive probe of these changes is a vital tool in our understanding of enzyme function. In this study we use ultrafast transient absorption spectroscopy to investigate the photoinduced electron transfer rates, which are also very sensitive to small changes in distance, between coenzyme analog, NAD(P)H4, and the isoalloxazine center in the model flavoenzymes morphinone reductase (wild-type and selected variants) and pentaerythritol tetranitrate reductase (wild-type). It is shown that upon addition of coenzyme to the protein the rate of photoinduced electron transfer is increased. By comparing the magnitude of this increase with existing values for NAD(P)H4-FMN distances, based on charge-transfer complex absorbance and experimental kinetic isotope effect reaction data, we show that this method can be used as a sensitive probe of donor-acceptor distance in a range of enzyme systems.

MD + QM correlations with tryptophan fluorescence spectral shifts and lifetimes

Principles behind quenching of tryptophan (Trp) fluorescence are updated and extended in light of recent 100-ns and 1-μs molecular dynamics (MD) trajectories augmented with quantum mechanical (QM) calculations that consider electrostatic contributions to wavelength shifts and quenching. Four studies are summarized, including (1) new insight into the single exponential decay of NATA, (2) a study revealing how unsuspected rotamer transitions affect quenching of Trp when used as a probe of protein folding, (3) advances in understanding the origin of nonexponential decay from 100-ns simulations on 19 Trps in 16 proteins, and (4) the correlation of wavelength with lifetime for decay-associated spectra (DAS). Each study strongly reinforces the concept that-for Trp-electron transfer-based quenching is controlled much more by environment electrostatic factors affecting the charge transfer (CT) state energy than by distance dependence of electronic coupling. In each case, water plays a large role in unexpected ways.

Ultrafast real-time visualization of active site flexibility of flavoenzyme thymidylate synthase ThyX

In many bacteria the flavoenzyme thymidylate synthase ThyX produces the DNA nucleotide deoxythymidine monophosphate from dUMP, using methylenetetrahydrofolate as carbon donor and NADPH as hydride donor. Because all three substrates bind in close proximity to the catalytic flavin adenine dinucleotide group, substantial flexibility of the ThyX active site has been hypothesized. Using femtosecond time-resolved fluorescence spectroscopy, we have studied the conformational heterogeneity and the conformational interconversion dynamics in real time in ThyX from the hyperthermophilic bacterium Thermotoga maritima. The dynamics of electron transfer to excited flavin adenine dinucleotide from a neighboring tyrosine residue are used as a sensitive probe of the functional dynamics of the active site. The fluorescence decay spanned a full three orders of magnitude, demonstrating a very wide range of conformations. In particular, at physiological temperatures, multiple angstrom cofactor-residue displacements occur on the picoseconds timescale. These experimental findings are supported by molecular dynamics simulations. Binding of the dUMP substrate abolishes this flexibility and stabilizes the active site in a configuration where dUMP closely interacts with the flavin cofactor and very efficiently quenches fluorescence itself. Our results indicate a dynamic selected-fit mechanism where binding of the first substrate dUMP at high temperature stabilizes the enzyme in a configuration favorable for interaction with the second substrate NADPH, and more generally have important implications for the role of active site flexibility in enzymes interacting with multiple poly-atom substrates and products. Moreover, our data provide the basis for exploring the effect of inhibitor molecules on the active site dynamics of ThyX and other multisubstrate flavoenzymes.

Simulations of tryptophan fluorescence dynamics during folding of the villin headpiece

Protein folding kinetics is commonly monitored by changes in tryptophan (Trp) fluorescence intensity. Considerable recent discussion has centered on whether the fluorescence of the single Trp in the much-studied, fast-folding villin headpiece C-terminal domain (HP35) accurately reflects folding kinetics, given the general view that quenching is by a histidine cation (His(+)) one turn away in an α-helix (helix III) that forms early in the folding process, according to published MD simulations. To help answer this question, we ran 1.0 μs MD simulations on HP35 (N27H) and a faster-folding variant in its folded form at 300 K and used the coordinates and force field charges with quantum calculations to simulate fluorescence quenching caused by electron transfer to the local amide and to the His(+). The simulations demonstrate that quenching by His(+) in the fully formed helix III is possible only during certain Trp and His(+) rotamer and solvent conformations, the propensity of which is a variable that can allow Trp fluorescence to report the global folding rate, as recent experiments imply.

A temperature-dependent conformational change of NADH oxidase from Thermus thermophilus HB8

Using molecular dynamics simulations and steady-state fluorescence spectroscopy, we have identified a conformational change in the active site of a thermophilic flavoenzyme, NADH oxidase from Thermus thermophilus HB8 (NOX). The enzyme's far-UV circular dichroism spectrum, intrinsic tryptophan fluorescence, and apparent molecular weight measured by dynamic light scattering varied little between 25 and 75°C. However, the fluorescence of the tightly bound FAD cofactor increased approximately fourfold over this temperature range. This effect appears not to be due to aggregation, unfolding, cofactor dissociation, or changes in quaternary structure. We therefore attribute the change in flavin fluorescence to a temperature-dependent conformational change involving the NOX active site. Molecular dynamics simulations and the effects of mutating aromatic residues near the flavin suggest that the change in fluorescence results from a decrease in quenching by electron transfer from tyrosine 137 to the flavin.

Evidence for concerted electron proton transfer in charge recombination between FADH- and 306Trp• in Escherichia coli photolyase

Proton-coupled electron-transfer (PCET) is a mechanism of great importance in protein electron transfer and enzyme catalysis, and the involvement of aromatic amino acids in this process is of much interest. The DNA repair enzyme photolyase provides a natural system that allows for the study of PCET using a neutral radical tryptophan (Trp(•)). In Escherichia coli photolyase, photoreduction of the flavin adenine dinucleotide (FAD) cofactor in its neutral radical semiquinone form (FADH(•)) results in the formation of FADH(-) and (306)Trp(•). Charge recombination between these two intermediates requires the uptake of a proton by (306)Trp(•). The rate constant of charge recombination has been measured as a function of temperature in the pH range from 5.5 to 10.0, and the data are analyzed with both classical Marcus and semi-classical Hopfield electron transfer theory. The reorganization energy associated with the charge recombination process shows a pH dependence ranging from 2.3 eV at pH ≤ 7 and 1.2 eV at pH(D) 10.0. These findings indicate that at least two mechanisms are involved in the charge recombination reaction. Global analysis of the data supports the hypothesis that PCET during charge recombination can follow two different mechanisms with an apparent switch around pH 6.5. At lower pH, concerted electron proton transfer (CEPT) is the favorable mechanism with a reorganization energy of 2.1-2.3 eV. At higher pH, a sequential mechanism becomes dominant with rate-limiting electron-transfer followed by proton uptake which has a reorganization energy of 1.0-1.3 eV. The observed 'inverse' deuterium isotope effect at pH < 8 can be explained by a solvent isotope effect that affects the free energy change of the reaction and masks the normal, mass-related kinetic isotope effect that is expected for a CEPT mechanism. To the best of our knowledge, this is the first time that a switch in PCET mechanism has been observed in a protein.

Effects of the disappearance of one charge on ultrafast fluorescence dynamics of the FMN binding protein

Crystal structures of E13T (Glu13 was replaced by Thr13) and E13Q (Glu13 was replaced by Gln13) FMN binding proteins (FMN-bp) from Desulfovibrio vulgaris, strain Miyazaki F, were determined by the X-ray diffraction method. Geometrical factors related to photoinduced electron transfer from Trp32, Tyr35, and Trp106 to the excited isoalloxazine (Iso*) were compared among the three forms of FMN-bp. The rate of ET is considered to be fastest from Trp32 to Iso* in FMN-bp and then from Tyr35 and Trp106. The distances between Iso and Trp32 did not change appreciably (0.705-0.712 nm) among WT, E13T, and E13Q FMN-bps, though the distances between Iso and Tyr35 or Trp106 became a little shorter by ca. 0.01 nm in both mutated FMN-bps. The distances between the residue at 13 and the ET donors or acceptor in the mutated proteins, however, changed markedly, compared to WT. Hydrogen bonding pairs and distances between Iso and surrounding amino acids were not modified when Glu13 was replaced by Thr13 or Gln13. Effects of elimination of ionic charge at Glu13 on the ultrafast fluorescence dynamics in E13T and E13Q were investigated comparing to WT, by means of a fluorescence up-conversion method. Fluorescence lifetimes were tau(1) = 107 fs (alpha(1) = 0.86), tau(2) = 475 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.02) in E13T and tau(1) = 134 fs (alpha(1) = 0.85), alpha(2) = 746 fs (alpha(2) = 0.12), and tau(3) = 30 ps (alpha(3) = 0.03) in E13Q, which are compared to the reported lifetimes in WT, tau(1) = 168 fs (alpha(1) = 0.95) and alpha(2) = 1.4 ps (alpha(2) = 0.05). Average lifetimes (tau(AV) = Sigma(i=1)(2or3)alpha(i)tau(i)) were 0.75 ps in E13T, 1.10 ps in E13Q, and 0.23 ps in WT, which implies that tau(AV) was 3.3 times longer in E13T and 4.8 times longer in E13Q, compared to WT. The ultrafast fluorescence dynamics of WT did not change when solvent changed from H(2)O to D(2)O. Static ET rates (inverse of average lifetimes) were analyzed with static structures of the three systems of FMN-bp. Net electrostatic (ES) energies of Iso and Trp32, on which ET rates depend, were 0.0263 eV in WT, 0.322 eV in E13T, and 0.412 eV in E13Q. The calculated ET rates were in excellent agreement with the observed ones in all systems.

Femtosecond fluorescence spectra of tryptophan in human gamma-crystallin mutants: site-dependent ultrafast quenching

The eye lens Crystallin proteins are subject to UV irradiation throughout life, and the photochemistry of damage proceeds through the excited state; thus, their tryptophan (Trp) fluorescence lifetimes are physiologically important properties. The time-resolved fluorescence spectra of single Trps in human gammaD- and gammaS-Crystallins have been measured with both an upconversion spectrophotofluorometer on the 300 fs to 100 ps time scale, and a time correlated single photon counting apparatus on the 100 ps to 10 ns time scale, respectively. Three Trps in each wild type protein were replaced by phenylalanine, leading to single-Trp mutants: W68-only and W156-only of HgammaD- and W72-only and W162-only of HgammaS-Crystallin. These proteins exhibit similar ultrafast signatures: positive definite decay associated spectra (DAS) for 50-65 ps decay constants that indicate dominance of fast, heterogeneous quenching. The quenched population (judged by amplitude) of this DAS differs among mutants. Trps 68, 156 in human gammaD- and Trp72 in human gammaS-Crystallin are buried, but water can reach amide oxygen and ring HE1 atoms through narrow channels. QM-MM simulations of quenching by electron transfer predict heterogeneous decay times from 50-500 ps that agree with our experimental results. Further analysis of apparent radiative lifetimes allow us to deduce that substantial subpopulations of Trp are fully quenched in even faster (sub-300 fs) processes for several of the mutants. The quenching of Trp fluorescence of human gammaD- and gammaS-Crystallin may protect them from ambient light induced photo damage.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Simultaneous analysis of ultrafast fluorescence decays of FMN binding protein and its mutated proteins by molecular dynamic simulation and electron transfer theory

Ultrafast fluorescence decays of FMN binding proteins (FBP) from Desulfovibrio vulgaris (Miyazaki F) were analyzed with an electron transfer (ET) theory by Kakitani and Mataga (KM theory). Time-dependent distances among isoalloxazine (Iso) and Trp-32, Tyr-35, and Trp-106 in wild-type FBP (WT), among Iso and Tyr-32, Tyr-35, and Trp-106 in W32Y (Trp-32 was replaced by Tyr-32), and among Iso and Tyr-35 and Trp-106 in W32A (Trp-32 was replaced by Ala-32) were determined by molecular dynamic simulation (MD). Electrostatic energies between Iso anion and all other ionic groups, between Trp-32 cation and all other ionic groups, and between Tyr-32 cation and all other ionic groups were calculated in WT, W32Y, and W32A, from the MD coordinates. ET parameters contained in KM theory, such as frequency (nu 0), a coefficient of the ET process (beta), a critical distance of the ET process ( R 0), standard free energy related to the electron affinity of the excited Iso ( G Iso (0)), and the static dielectric constant in FBP species (epsilon 0), were determined with and without inclusion of the electrostatic energy, so as to fit the calculated fluorescence decays with the observed decays of all FBP species, by a nonlinear least-squares method according to the Marquardt algorithm. In the analyses the parameters, nu 0, beta, and R 0 were determined separately between Trp residues and Tyr residues among all FBP species. Calculated fluorescence intensities with the inclusion of the electrostatic energy fit quite well with the observed ones of all WT, W32Y, and W32A.

Unraveling the mechanisms of tryptophan fluorescence quenching in the triosephosphate isomerase from Giardia lamblia

In the native state several proteins exhibit a quenching of fluorescence of their tryptophans. We studied triosephosphate isomerase from Giardia lamblia (GlTIM) to dissect the mechanisms that account for the quenching of fluorescence of its Trp. GlTIM contains four Trp per monomer (Trp75, Trp162, Trp173, and Trp196) distributed throughout the 3D structure. The fluorescence of the denatured enzyme is 3-fold higher than that of native GlTIM. To ascertain the origin of this phenomenon, single and triple mutants of Trp per Phe were made. The intrinsic fluorescence was determined, and the data were interpreted on the basis of the crystal structure of the enzyme. Our data show that the fluorescence of all Trp residues is quenched through two different mechanisms. In one, fluorescence is quenched by aromatic-aromatic interactions due to the proximity and orientation of the indole groups of Trp196 and Trp162. The magnitude of the quenching of fluorescence in Trp162 is higher than in the other three Trp. Fluorescence quenching is also due to energy transfer to the charged residues that surround Trp 75, 173 and 196. Further analysis of the fluorescence of GlTIM showed that, among TIMs from other parasites, Trp at position 12 exhibits rather unique properties.

Comparison between Ultrafast Fluorescence Dynamics of FMN Binding Protein from Desulfovibrio vulgaris, Strain Miyazaki, in Solution vs Crystal Phases

Ultrafast fluorescence dynamics of FMN binding protein (FBP) from Desulfobivrio vulgaris, strain Miyaxaki F, were compared in solution and crystal phases. Fluorescence lifetimes of FBP were 167 fs (96%) and 1.5 ps (4%) in solution (tau(av) = 220 fs), and 730 fs (60%) and longer than 10 ps (40%) in crystals (tau(av) = 4.44 ps). The quenching of the fluorescence of flavin in the protein was considered to be due to photoinduced electron transfer (ET) from Trp or Tyr to the excited isoalloxazine (Iso) nearby. The average lifetime was 20 times longer in crystal vs in solution. Averaged distances between Iso and nearby Trp-32, Tyr-35, and Trp-106 were 8.42, 7.36, and 8.15 A in solution, respectively (obtained by NMR spectroscopy), and 7.05, 7.72, and 8.49 A in crystal, respectively (obtained by X-ray crystallography). The prolonged lifetime in crystal cannot be elucidated by the change in the distances between the states. It was suggested that the longer lifetime in crystal was ascribed to the absence of water molecules around FBP with rapid motional freedom, which may be the driving force for the ET in flavoproteins.

Fluorescence probe of Trp-cage protein conformation in solution and in gas phase

Measurements of protein unfolding in the absence of solvent, when combined with unfolding studies in solution, offer a unique opportunity to measure the effects of solvent on protein structure and dynamics. The experiments presented here rely on the fluorescence of an attached dye to probe the local conformational dynamics through interactions with a Trp residue and fields originating on charge sites. We present fluorescence measurements of thermal fluctuations accompanying conformational change of a miniprotein, Trp-cage, in solution and in gas phase. Molecular dynamics (MD) simulations are performed as a function of temperature, charge state, and charge location to elucidate the dye-protein conformational dynamics leading to the changes in measured fluorescence. The results indicate that the stability of the unsolvated protein is dominated by hydrogen bonds. Substituting asparagine for aspartic acid at position 9 results in a dramatic alteration of the solution unfolding curve, indicating that the salt bridge involving Lys8, Asp9, and Arg16 (+ - +) is essential for Trp-cage stability in solution. In contrast, this substitution results in minor changes in the unfolding curve of the unsolvated protein, showing that hydrogen bonds are the major contributor to the stability of Trp-cage in gas phase. Consistent with this hypothesis, the decrease in the number of hydrogen bonds with increasing temperature indicated by MD simulations agrees reasonably well with the experimentally derived enthalpies of conformational change. The simulation results display relatively compact conformations compared with NMR structures that are generally consistent with experimental results. The measured unfolding curves of unsolvated Trp-cage ions are invariant with the acetonitrile content of the solution from which they are formed, possibly as a result of conformational relaxation during or after desolvation. This work demonstrates the power of combined solution and gas-phase studies and of single-point mutations to identify specific noncovalent interactions which contribute to protein-fold stability. The combination of experiment and simulation is particularly useful because these approaches yield complementary information which can be used to deduce the details of structural changes of proteins in the gas phase.

Shedding Light on Biomolecule Conformational Dynamics Using Fluorescence Measurements of Trapped Ions

Biomolecule conformational change has been widely investigated in solution using several methods; however, much less experimental data about structural changes are available for completely isolated, gas-phase biomolecules. Studies of conformational change in unsolvated biomolecules are required to complement the interpretation of mass spectrometry measurements and in addition, can provide a means to directly test theoretical simulations of biomolecule structure and dynamics independent of a simulated solvent. In this Feature Article, we review our recent introduction of a fluorescence-based method for probing local conformational dynamics in unsolvated biomolecules through interactions of an attached dye with tryptophan (Trp) residues and fields originating on charge sites. Dye-derivatized biomolecule ions are formed by electrospray ionization and are trapped in a variable-temperature quadrupole ion trap in which they are irradiated with either continuous or short pulse lasers to excite fluorescence. Fluorescence is measured as a function of temperature for different charge states. Optical measurements of the dye fluorescence include average intensity changes, changes in the emission spectrum, and time-resolved measurements of the fluorescence decay. These measurements have been applied to the miniprotein, Trp-cage, polyproline peptides and to a beta-hairpin-forming peptide, and the results are presented as examples of the broad applicability and utility of these methods. Model fits to Trp-cage fluorescence data measured as a function of temperature provide quantitative information on the thermodynamics of conformational changes, which are reproduced well by molecular dynamics. Time-resolved measurements of the fluorescence decays of Trp-cage and small polyproline peptides definitively demonstrate the occurrence of fluorescence quenching by the amino acid Trp in unsolvated biomolecules.