Ab Initio Prediction of Tryptophan Fluorescence Quenching by Protein Electric Field Enabled Electron Transfer

Multifarious analytical capabilities of the UV/Vis protein fluorescence in blood plasma

Autofluorescence of blood plasma has been broadly considered as a prospective disease screening method. However, the assessment of such intrinsic fluorescence is mostly phenomenological, and its origin is still not fully understood, complicating its use in the clinical practice. Here we present the detailed evaluation of analytical capabilities, variability, and formation of blood plasma protein fluorescence based on the open dataset of excitation-emission matrices measured for ∼300 patients with suspected colorectal cancer, and our supporting model experiments. Using high-resolution size-exclusion chromatography coupled with comprehensive spectral analysis, we demonstrate, for the first time, the dominant role of HSA in the formation of blood plasma fluorescence in the visible spectral range (excitation wavelength >350 nm), presumably caused by its oxidative modifications. Furthermore, the diagnostic value of the tryptophan emission, as well as of the tyrosine fluorescence and visible fluorescence of proteins is shown by building a tree-based classification model that uses a small subset of physically interpretable fluorescence features for distinguishing between the control group and cancer patients with >80% accuracy. The obtained results extend current understanding and approaches used for the analysis of blood plasma fluorescence and pave the way for novel autofluorescence-based disease screening methods.

Fluorescence and quantum mechanical approach on the interaction of amides and their role on the stability and coexistence of the rotamer conformations of L-tryptophan in aqueous solution

Photophysical investigation on the fluorescence decay characteristics of L-tryptophan and a derivative N-acetyl-L-tryptophanamide (NATA) with alkyl amides were carried out in water. L-tryptophan exists in the zwitterionic form and exhibits a biexponential lifetime which is correlated to the existence of rotamer structures. Addition of formamide (F) and dimethylformamide (DMF) results in a decrease in the fluorescence lifetime and its proportion of the most stable structure of L-tryptophan wherein acetamide (ACM) results in an increase of the same. Interestingly, all the amides result in the formation of the lifetime of the rotamer whose lifetime doesn't exist initially and the lifetime and its distribution increases irrespective of the nature of amide. The interaction between L-tryptophan and amide is attributed to hydrogen-bonding such that these interactions influence the relative proportion of the existence of individual rotamers in the presence of amides.Strikingly, in the case of NATA that does not exhibit rotamer structures; the fluorescence lifetime is quenched in the presence of F, whereas ACM and DMF result in a larger fold of enhancement resulting in two different lifetimes. The variation in the fluorescence lifetime and amplitude of the various conformers of L-tryptophan and of NATA is completely governed by the concentration of the amides in solution such that the microenvironment surrounding the fluorophores are completely reorganised. The hydrogen-bonding functional groups in amides that are responsible for the coexistence of rotamers are elucidated and well supported by quantum mechanical (QM) studies. Time-correlated single-photon counting(TCSPC) technique is used as a probe as well as marker in establishing the variation in the lifetime properties of L-tryptophan and NATA with non-fluorescent hydrogen-bonding solutes in water which promotes this as fascinating field of research in the context of fluorescence properties of a complicated amino acid-like tryptophan.

Ionic Hydrogen Bonds and Lipid Packing Defects Determine the Binding Orientation and Insertion Depth of RecA on Multicomponent Lipid Bilayers

We describe a computational and experimental approach for probing the binding properties of the RecA protein at the surface of anionic membranes. Fluorescence measurements indicate that RecA behaves differently when bound to phosphatidylglycerol (PG)- and cardiolipin (CL)-containing liposomes. We use a multistage computational protocol that integrates an implicit membrane/solvent model, the highly mobile mimetic membrane model, and the full atomistic membrane model to study how different anionic lipids perturb RecA binding to the membrane. With anionic lipids studied here, the binding interface involves three key regions: the N-terminal helix, the DNA binding loop L2, and the M-M7 region. The nature of binding involves both electrostatic interactions between cationic protein residues and lipid polar/charged groups and insertion of hydrophobic residues. The L2 loop contributes more to membrane insertion than the N-terminal helix. More subtle aspects of RecA-membrane interaction are influenced by specific properties of anionic lipids. Ionic hydrogen bonds between the carboxylate group in phosphatidylserine and several lysine residues in the C-terminal region of RecA stabilize the parallel (∥) binding orientation, which is not locally stable on PG- and CL-containing membranes despite similarity in the overall charge density. Lipid packing defects, which are more prevalent in the presence of conical lipids, are observed to enhance the insertion depth of hydrophobic motifs. The computational finding that RecA binds in a similar orientation to PG- and CL-containing membranes is consistent with the fact that PG alone is sufficient to induce RecA polar localization, although CL might be more effective because of its tighter binding to RecA. The different fluorescence behaviors of RecA upon binding to PG- and CL-containing liposomes is likely due to the different structures and flexibility of the C-terminal region of RecA when it binds to different anionic phospholipids.

Combined use of optical spectroscopy and computational methods to study the binding and the photoinduced conformational modification of proteins when NMR and X-ray structural determinations are not an option

The functions of proteins depend on their interactions with various ligands and these interactions are controlled by the structure of the polypeptides. If one can manipulate the structure of proteins, their functions can in principle be modulated. The issue of protein structure-function relationship is not only a central problem in biophysics, but is becoming clear that the ability to "artificially" modify the structure of proteins could be relevant in fields beyond the biomedical area to provide, for instance, light responses in proteins which would not possess such properties in their native state. This chapter presents an overview of the combination of optical electronic and vibrational spectroscopy with various computational methods to investigate the binding between photoactive ligands and proteins.

Fluorescence and phosphorescence of tryptophan in peptides of different length and sequence

To interpret accurately protein fluorescence and phosphorescence, it is essential to achieve a better understanding of the luminescence properties of tryptophan (Trp, or W) in peptides. In published literature data on luminescence of peptides of varied length are scarce. This article describes studies of fluorescence and phosphorescence properties of the eight Trp-containing synthetic peptides: WAK, AWK, SWA, KYLWE, AVSWK, WVSWAK, WAKLAWE, and AVSWAKLARE. The aim was to investigate which factors influence the fluorescence yield and phosphorescence-spectra and lifetimes. Absorption spectra, room temperature fluorescence emission and corresponding excitation spectra and time-resolved phosphorescence spectra (77K) have been recorded; the dependence of the fluorescence quantum yield on the specific peptide and its variation with the wavelength of excitation has been studied. The changes in fluorescence yield and shape of phosphorescence spectra are explained in terms of internal electron and proton transfer. The structured phosphorescence spectrum originates from proton transfer occurring upon excitation of Trp, while electron transfer gives rise to a non-structured luminescence spectrum. There is also electron transfer from higher vibronic S1 states. In the peptides there is higher probability of electron transfer than in Trp alone. The obtained data are interpreted in light of the peptides' sequence, length and conformation.

Fluorescence Characteristics of [(Benzoyloxy)methyl]anthracene Donor-Acceptor Systems

A theoretical investigation on fluorescence properties of [(benzoyloxy)methyl]anthracene derivatives containing different groups (OCH3, CH3, H, CF3, F, CN, and NO2) as substituent on the phenyl ring is reported. Electron transfer rate constants for the molecules were calculated using Marcus theory. Theoretical electron transfer rate constants agreed with experimentally observed trend of fluorescence quenching. Large electron transfer rate constants were obtained for molecules containing strongly electron withdrawing groups as the substituent on the phenyl ring. Calculations were conducted at Hartree-Fock and density functional (HF/6-31G(d) and B3LYP/6-31G(d)) levels of theory. Density functional theory predicted spurious charge transfer excited states for molecules containing NO2, CN, and CF3 as the substituent on the phenyl ring.

Exploring protein solution structure: Second moments of fluorescent spectra report heterogeneity of tryptophan rotamers

Trp fluorescent spectra appear as a log-normal function but are usually analyzed with λmax, full width at half maximum, and the first moment of incomplete spectra. Log-normal analyses have successfully separated fluorescence contributions from some multi-Trp proteins but deviations were observed in single Trp proteins. The possibility that disparate rotamer environments might account for these deviations was explored by moment spectral analysis of single Trp mutants spanning the sequence of tear lipocalin as a model. The analysis required full width Trp spectra. Composite spectra were constructed using log-normal analysis to derive the inaccessible blue edge, and the experimentally obtained spectra for the remainder. First moments of the composite spectra reflected the site-resolved secondary structure. Second moments were most sensitive for spectral deviations. A novel parameter, derived from the difference of the second moments of composite and simulated log-normal spectra correlated with known multiple heterogeneous rotamer conformations. Buried and restricted side chains showed the most heterogeneity. Analyses applied to other proteins further validated the method. The rotamer heterogeneity values could be rationalized by known conformational properties of Trp residues and the distribution of nearby charged groups according to the internal Stark effect. Spectral heterogeneity fits the rotamer model but does not preclude other contributing factors. Spectral moment analysis of full width Trp emission spectra is accessible to most laboratories. The calculations are informative of protein structure and can be adapted to study dynamic processes.

The broken ring: reduced aromaticity in Lys-Trp cations and high pH tautomer correlates with lower quantum yield and shorter lifetimes

Several nonradiative processes compete with tryptophan fluorescence emission. The difficulty in spectral interpretation lies in associating specific molecular environmental features with these processes and thereby utilizing the fluorescence spectral data to identify the local environment of tryptophan. Here, spectroscopic and molecular modeling study of Lys-Trp dipeptide charged species shows that backbone-ring interactions are undistinguished. Instead, quantum mechanical ground state isosurfaces reveal variations in indole π electron distribution and density that parallel charge (as a function of pK(1), pK(2), and pK(R)) on the backbone and residues. A pattern of aromaticity-associated quantum yield and fluorescence lifetime changes emerges. Where quantum yield is high, isosurfaces have a charge distribution similar to the highest occupied molecular orbital (HOMO) of indole, which is the dominant fluorescent ground state of the (1)L(a) transition dipole moment. Where quantum yield is low, isosurface charge distribution over the ring is uneven, diminished, and even found off ring. At pH 13, the indole amine is deprotonated, and Lys-Trp quantum yield is extremely low due to tautomer structure that concentrates charge on the indole amine; the isosurface charge distribution bears scant resemblance to the indole HOMO. Such greatly diminished fluorescence has been observed for proteins where the indole nitrogen is hydrogen bonded, lending credence to the association of aromaticity changes with diminished quantum yield in proteins as well. Thus tryptophan ground state isosurfaces are an indicator of indole aromaticity, signaling the partition of excitation energy between radiative and nonradiative processes.

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.

Implementing electrostatic polarization cannot fill the gap between experimental and theoretical measurements for the ultrafast fluorescence decay of myoglobin

Over the past few years, time-dependent ultrafast fluorescence spectroscopy method has been applied to the study of protein dynamics. However, observations from these experiments are in a controversy with other experimental studies. Participating of theoretical methods in this debate has not reconciled the contradiction, because the predicted initial relaxation from computer simulations is one-order faster than the ultrafast fluorescence spectroscopy experiment. In those simulations, pairwise force fields are employed, which have been shown to underestimate the roughness of the free energy landscape. Therefore, the relaxation rate of protein and water molecules under pairwise force fields is falsely exaggerated. In this work, we compared the relaxations of tryptophan/environment interaction under linear response approximation employing pairwise, polarized, and polarizable force fields. Results show that although the relaxation can be slowed down to a certain extent, the large gap between experiment and theory still cannot be filled.

Effect of short- and long-range interactions on trp rotamer populations determined by site-directed tryptophan fluorescence of tear lipocalin

In the lipocalin family, the conserved interaction between the main α-helix and the β-strand H is an ideal model to study protein side chain dynamics. Site-directed tryptophan fluorescence (SDTF) has successfully elucidated tryptophan rotamers at positions along the main alpha helical segment of tear lipocalin (TL). The rotamers assigned by fluorescent lifetimes of Trp residues corroborate the restriction expected based on secondary structure. Steric conflict constrains Trp residues to two (t, g⁻) of three possible χ₁ (t, g⁻, g⁺) canonical rotamers. In this study, investigation focused on the interplay between rotamers for a single amino acid position, Trp 130 on the α-helix and amino acids Val 113 and Leu 115 on the H strand, i.e. long range interactions. Trp130 was substituted for Phe by point mutation (F130W). Mutations at positions 113 and 115 with combinations of Gly, Ala, Phe residues alter the rotamer distribution of Trp130. Mutations, which do not distort local structure, retain two rotamers (two lifetimes) populated in varying proportions. Replacement of either long range partner with a small amino acid, V113A or L115A, eliminates the dominance of the t rotamer. However, a mutation that distorts local structure around Trp130 adds a third fluorescence lifetime component. The results indicate that the energetics of long-range interactions with Trp 130 further tune rotamer populations. Diminished interactions, evident in W130G113A115, result in about a 22% increase of α-helix content. The data support a hierarchic model of protein folding. Initially the secondary structure is formed by short-range interactions. TL has non-native α-helix intermediates at this stage. Then, the long-range interactions produce the native fold, in which TL shows α-helix to β-sheet transitions. The SDTF method is a valuable tool to assess long-range interaction energies through rotamer distribution as well as the characterization of low-populated rotameric states of functionally important excited protein states.

Pan1 is an intrinsically disordered protein with homotypic interactions

The yeast scaffold protein Pan1 contains two EH domains at its N-terminus, a predicted coiled-coil central region, and a C-terminal proline-rich domain. Pan1 is also predicted to contain regions of intrinsic disorder, characteristic of proteins that have many binding partners. In vitro biochemical data suggest that Pan1 exists as a dimer, and we have identified amino acids 705 to 848 as critical for this homotypic interaction. Tryptophan fluorescence was used to further characterize Pan1 conformational states. Pan1 contains four endogenous tryptophans, each in a distinct region of the protein: Trp(312) and Trp(642) are each in an EH domain, Trp(957) is in the central region, and Trp(1280) is a critical residue in the Arp2/3 activation domain. To examine the local environment of each of these tryptophans, three of the four tryptophans were mutagenized to phenylalanine to create four proteins, each with only one tryptophan residue. When quenched with acrylamide, these single tryptophan mutants appeared to undergo collisional quenching exclusively and were moderately accessible to the acrylamide molecule. Quenching with iodide or cesium, however, revealed different Stern-Volmer constants due to unique electrostatic environments of the tryptophan residues. Time-resolved fluorescence anisotropy data confirmed structural and disorder predictions of Pan1. Further experimentation to fully develop a model of Pan1 conformational dynamics will assist in a deeper understanding of the mechanisms of endocytosis.

Fluorescence of tryptophan in designed hairpin and Trp-cage miniproteins: measurements of fluorescence yields and calculations by quantum mechanical molecular dynamics simulations

The quantum yield of tryptophan (Trp) fluorescence was measured in 30 designed miniproteins (17 β-hairpins and 13 Trp-cage peptides), each containing a single Trp residue. Measurements were made in D(2)O and H(2)O to distinguish between fluorescence quenching mechanisms involving electron and proton transfer in the hairpin peptides, and at two temperatures to check for effects of partial unfolding of the Trp-cage peptides. The extent of folding of all the peptides also was measured by NMR. The fluorescence yields ranged from 0.01 in some of the Trp-cage peptides to 0.27 in some hairpins. Fluorescence quenching was found to occur by electron transfer from the excited indole ring of the Trp to a backbone amide group or the protonated side chain of a nearby histidine, glutamate, aspartate, tyrosine, or cysteine residue. Ionized tyrosine side chains quenched strongly by resonance energy transfer or electron transfer to the excited indole ring. Hybrid classical/quantum mechanical molecular dynamics simulations were performed by a method that optimized induced electric dipoles separately for the ground and excited states in multiple π-π* and charge-transfer (CT) excitations. Twenty 0.5 ns trajectories in the tryptophan's lowest excited singlet π-π* state were run for each peptide, beginning by projections from trajectories in the ground state. Fluorescence quenching was correlated with the availability of a CT or exciton state that was strongly coupled to the π-π* state and that matched or fell below the π-π* state in energy. The fluorescence yields predicted by summing the calculated rates of charge and energy transfer are in good accord with the measured yields.

Tryptophan rotamer distribution revealed for the α-helix in tear lipocalin by site-directed tryptophan fluorescence

Rotamer libraries are a valuable tool for protein structure determination, modeling, and design. Site-directed tryptophan fluorescence (SDTF) was used in combination with the rotamer model for the fluorescence intensity decays to solve α-helical conformations of proteins in solution. Single Trp mutations located in an α-helical segment of human tear lipocalin were explored for structure assignment. Along with fluorescence λ(max) values, the rotamer model assignment of fluorescence lifetimes fits the backbone conformation. Typically, Trp fluorescence in proteins shows three lifetimes. However, for the α-helix, two lifetimes assigned to t and g(-) rotamers were satisfactory to describe Trp fluorescence intensity decays. The g(+) rotamer is not feasible in the α-helix due to steric restriction. Trp rotamer distributions obtained by fluorescence were compared with the rotamer library derived from X-ray crystallography data of proteins. The Trp rotamer distributions vary for solvent exposed and buried (tertiary interaction) sites. A new strategy using the rotamer distribution with SDTF (RD-SDTF) removes the limitation of regular SDTF and other labeling techniques, in which site-specific differences, e.g., accessibility, are presumed. The RD-SDTF technique does not rely on environmental differences of side chains and is able to detect α-helical structure where all side chains are exposed to solvent. Potentially, this technique is applicable to various proteins including membrane proteins, which are rich in α-helix motif.

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.

Atomistic molecular simulations of protein folding

Theory and experiment have provided answers to many of the fundamental questions of protein folding; a remaining challenge is an accurate, high-resolution picture of folding mechanism. Atomistic molecular simulations with explicit solvent are the most promising method for providing this information, by accounting more directly for the physical interactions that stabilize proteins. Although simulations of folding with such force fields are extremely challenging, they have become feasible as a result of recent advances in computational power, accuracy of the energy functions or 'force fields', and methods for improving sampling of folding events. I review the recent progress in these areas, and highlight future challenges and questions that we may hope to address with these methods. I also attempt to place atomistic models into the context of the energy landscape view of protein folding, and coarse-grained simulations.

Making connections between ultrafast protein folding kinetics and molecular dynamics simulations

Determining the rate of forming the truly folded conformation of ultrafast folding proteins is an important issue for both experiments and simulations. The double-norleucine mutant of the 35-residue villin subdomain is the focus of recent computer simulations with atomistic molecular dynamics because it is currently the fastest folding protein. The folding kinetics of this protein have been measured in laser temperature-jump experiments using tryptophan fluorescence as a probe of overall folding. The conclusion from the simulations, however, is that the rate determined by fluorescence is significantly larger than the rate of overall folding. We have therefore employed an independent experimental method to determine the folding rate. The decay of the tryptophan triplet-state in photoselection experiments was used to monitor the change in the unfolded population for a sequence of the villin subdomain with one amino acid difference from that of the laser temperature-jump experiments, but with almost identical equilibrium properties. Folding times obtained in a two-state analysis of the results from the two methods at denaturant concentrations varying from 1.5-6.0 M guanidinium chloride are in excellent agreement, with an average difference of only 20%. Polynomial extrapolation of all the data to zero denaturant yields a folding time of 220 (+100,-70) ns at 283 K, suggesting that under these conditions the barrier between folded and unfolded states has effectively disappeared--the so-called "downhill scenario."

Efavirenz binding to HIV-1 reverse transcriptase monomers and dimers

Efavirenz (EFV) is a nonnucleoside reverse transcriptase inhibitor (NNRTI) of HIV-1 reverse transcriptase (RT) used for the treatment of AIDS. RT is a heterodimer composed of p66 and p51 subunits; p51 is produced from p66 by C-terminal truncation by HIV protease. The monomers can form p66/p66 and p51/p51 homodimers as well as the p66/p51 heterodimer. Dimerization and efavirenz binding are coupled processes. In the crystal structure of the p66/p51-EFV complex, the drug is bound to the p66 subunit. The binding of efavirenz to wild-type and dimerization-defective RT proteins was studied by equilibrium dialysis, tryptophan fluorescence, and native gel electrophoresis. A 1:1 binding stoichiometry was determined for both monomers and homodimers. Equilibrium dissociation constants are approximately 2.5 microM for both p66- and p51-EFV complexes, 250 nM for the p66/p66-EFV complex, and 7 nM for the p51/p51-EFV complex. An equilibrium dissociation constant of 92 nM for the p66/p51-EFV complex was calculated from the thermodynamic linkage between dimerization and inhibitor binding. Binding and unbinding kinetics monitored by fluorescence were slow. Progress curve analyses revealed a one-step, direct binding mechanism with association rate constants k(1) of approximately 13.5 M(-1) s(-1) for monomers and heterodimer and dissociation rate constants k(-1) of approximately 9 x 10(-5) s(-1) for monomers. A conformational selection mechanism is proposed to account for the slow association rate. These results show that efavirenz is a slow, tight-binding inhibitor capable of binding all forms of RT and suggest that the NNRTI binding site in monomers and dimers is similar.

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.

Crystalline indole at high pressure: chemical stability, electronic, and vibrational properties

Vibrational and electronic spectra of crystalline indole were measured up to 25.5 GPa at room temperature in a diamond anvil cell. In particular, Fourier transform infrared (FTIR) spectra in the mid-infrared region and two-photon excitation profiles and fluorescence spectra in the region of the HOMO-LUMO transitions were obtained. The analysis of the FTIR spectra revealed a large red-shift of the N-H stretching mode with increasing pressure, indicating the strengthening of the H-bond between the NH group and the pi electron density of nearest neighbor molecules. The frequencies of four vibronic bands belonging to the (1)L(a) and (1)L(b) systems were obtained as a function of pressure. Comparison with literature data shows that the crystal acts as a highly polar environment with regard to the position of the (1)L(b) origin and of the fluorescence maximum, which are largely red-shifted with respect to the gas phase or to solutions in apolar solvents. A large, and increasing with pressure, frequency difference between the (1)L(b) origin and the blue edge of the fluorescence spectrum suggests that the emitting state is (1)L(a), that is known to be more stabilized than (1)L(b) by dipolar relaxation. Crystalline indole was found to be very stable with respect to pressure-induced reactivity. Only traces of a reaction product, containing saturated C-H bonds, are detected after a full compression-decompression cycle. In addition, differently from many unsaturated compounds at high pressure, irradiation with light matching a two-photon absorption for a HOMO-LUMO transition has no enhancing effect on reactivity. The chemical stability of indole at high pressure is ascribed to the crystal structure, where nearest neighbor molecules, formig H-bonds, are not in a favorable position to react, while reaction between equivalent molecules, for which a superposition of the pi electron clouds would be possible, is hindered by H-bonded molecules. Consistently, no excimer emission was observed except at the cell opening at the end of the compression-decompression run. Extremely limited chemical reactivity and excimer formation likely occur at crystal defects, evidencing the strict connection between the two phenomena.

Mechanism of the very efficient quenching of tryptophan fluorescence in human gamma D- and gamma S-crystallins: the gamma-crystallin fold may have evolved to protect tryptophan residues from ultraviolet photodamage

Proteins exposed to UV radiation are subject to irreversible photodamage through covalent modification of tryptophans (Trps) and other UV-absorbing amino acids. Crystallins, the major protein components of the vertebrate eye lens that maintain lens transparency, are exposed to ambient UV radiation throughout life. The duplicated beta-sheet Greek key domains of beta- and gamma-crystallins in humans and all other vertebrates each have two conserved buried Trps. Experiments and computation showed that the fluorescence of these Trps in human gammaD-crystallin is very efficiently quenched in the native state by electrostatically enabled electron transfer to a backbone amide [Chen et al. (2006) Biochemistry 45, 11552-11563]. This dispersal of the excited state energy would be expected to minimize protein damage from covalent scission of the excited Trp ring. We report here both experiments and computation showing that the same fast electron transfer mechanism is operating in a different crystallin, human gammaS-crystallin. Examination of solved structures of other crystallins reveals that the Trp conformation, as well as favorably oriented bound waters, and the proximity of the backbone carbonyl oxygen of the n - 3 residues before the quenched Trps (residue n), are conserved in most crystallins. These results indicate that fast charge transfer quenching is an evolved property of this protein fold, probably protecting it from UV-induced photodamage. This UV resistance may have contributed to the selection of the Greek key fold as the major lens protein in all vertebrates.

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.

Mechanism of the efficient tryptophan fluorescence quenching in human gammaD-crystallin studied by time-resolved fluorescence

Human γD-crystallin (HγD-Crys) is a two-domain, β-sheet eye lens protein found in the lens nucleus. Its long-term solubility and stability are important to maintain lens transparency throughout life. HγD-Crys has four highly conserved buried tryptophans (Trps), with two in each of the homologous β-sheet domains. In situ, these Trps will be absorbing ambient UV radiation that reaches the lens. The dispersal of the excited-state energy to avoid covalent damage is likely to be physiologically relevant for the lens crystallins. Trp fluorescence is efficiently quenched in native HγD-Crys. Previous steady-state fluorescence measurements provide strong evidence for energy transfer from Trp42 to Trp68 in the N-terminal domain and from Trp130 to Trp156 in the C-terminal domain [Chen, J., et al. (2006) Biochemistry 45, 11552−11563]. Hybrid quantum mechanical−molecular mechanical (QM-MM) simulations indicated that the fluorescence of Trp68 and Trp156 is quenched by fast electron transfer to the amide backbone. Here we report additional information obtained using time-resolved fluorescence spectroscopy. In the single-Trp-containing proteins (Trp42-only, Trp68-only, Trp130-only, and Trp156-only), the highly quenched Trp68 and Trp156 have very short lifetimes, τ ∼0.1 ns, whereas the moderately fluorescent Trp42 and Trp130 have longer lifetimes, τ ∼3 ns. In the presence of the energy acceptor (Trp68 or Trp156), the lifetime of the energy donor (Trp42 or Trp130) decreased from ∼3 to ∼1 ns. The intradomain energy transfer efficiency is 56% in the N-terminal domain and is 71% in the C-terminal domain. The experimental values of energy transfer efficiency are in good agreement with those calculated theoretically. The absence of a time-dependent red shift in the time-resolved emission spectra of Trp130 proves that its local environment is very rigid. Time-resolved fluorescence anisotropy measurements with the single-Trp-containing proteins, Trp42-only and Trp130-only, indicate that the protein rotates as a rigid body and no segmental motion is detected. A combination of energy transfer with electron transfer results in short excited-state lifetimes of all Trps, which, together with the high rigidity of the protein matrix around Trps, could protect HγD-Crys from excited-state reactions causing permanent covalent damage.

Irradiation of the porphyrin causes unfolding of the protein in the protoporphyrin IX/beta-lactoglobulin noncovalent complex

Porphyrins such as protoporphyrin IX (PPIX) are known to occasionally cause conformational changes in proteins for which they are specific ligands. It has also been established that irradiation of porphyrins noncovalently intercalated between bases or bound to one of the grooves can cause conformational effects on DNA. Conversely, there is no evidence reported in the literature of conformational changes caused by noncovalently bound PPIX to globular proteins for which the porphyrin is not a specific ligand. This study shows that the irradiation of the porphyrin in the PPIX/lactoglobulin noncovalent complex indeed causes a local and limited (approximately 7%) unfolding of the protein near the location of Trp19. This event causes the intrinsic fluorescence spectrum of the protein to shift to the red by 2 nm and the average decay lifetime to lengthen by approximately 0.5 ns. The unfolding of lactoglobulin occurs only at pH >7 because of the increased instability of the protein at alkaline pH. The photoinduced unfolding does not depend on the presence of O2 in solution; therefore, it is not mediated by formation of singlet oxygen and is likely the result of electron transfer between the porphyrin and amino acid residues.