Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Mössbauer and electron paramagnetic resonance studies of chloroperoxidase following mechanism-based inactivation with allylbenzene

We have used Mössbauer and electron paramagnetic resonance (EPR) spectroscopy to study a heme-N-alkylated derivative of chloroperoxidase (CPO) prepared by mechanism-based inactivation with allylbenzene and hydrogen peroxide. The freshly prepared inactivated enzyme (“green CPO”) displayed a nearly pure low-spin ferric EPR signal with g = 1.94, 2.15, 2.31. The Mössbauer spectrum of the same species recorded at 4.2 K showed magnetic hyperfine splittings, which could be simulated in terms of a spin Hamiltonian with a complete set of hyperfine parameters in the slow spin fluctuation limit. The EPR spectrum of green CPO was simulated using a three-term crystal field model including g-strain. The best-fit parameters implied a very strong octahedral field in which the three 2T2 levels of the (3d)5 configuration in green CPO were lowest in energy, followed by a quartet. In native CPO, the 6A1 states follow the 2T2 ground state doublet. The alkene-mediated inactivation of CPO is spontaneously reversible. Warming of a sample of green CPO to 22°C for increasing times before freezing revealed slow conversion of the novel EPR species to two further spin S = ½ ferric species. One of these species displayed g = 1.82, 2.25, 2.60 indistinguishable from native CPO. By subtracting spectral components due to native and green CPO, a third species with g = 1.86, 2.24, 2.50 could be generated. The EPR spectrum of this “quasi-native CPO,” which appears at intermediate times during the reactivation, was simulated using best-fit parameters similar to those used for native CPO.

Large kinetic isotope effects in enzymatic proton transfer and the role of substrate oscillations

We propose an interpretation of the experimental findings of Klinman and coworkers [Cha, Y., Murray, C. J. & Klinman, J. P. (1989) Science 243, 1325–1330; Grant, K. L. & Klinman, J. P. (1989) Biochemistry 28, 6597–6605; and Bahnson, B. J. & Klinman, J. P. (1995) Methods Enzymol. 249, 373–397], who showed that proton transfer reactions that are catalyzed by bovine serum amine oxidase proceed through tunneling. We show that two different tunneling models are consistent with the experiments. In the first model, the proton tunnels from the ground state. The temperature dependence of the kinetic isotope effect is caused by a thermally excited substrate mode that modulates the barrier, as has been suggested by Borgis and Hynes [Borgis, D. & Hynes, J. T. (1991) J. Chem. Phys. 94, 3619–3628]. In the second model, there is both over-the-barrier transfer and tunneling from excited states. Finally, we propose two experiments that can distinguish between the possible mechanisms.

The change in hydrogen bond strength accompanying charge rearrangement: Implications for enzymatic catalysis

The equilibrium for formation of the intramolecular hydrogen bond (KHB) in a series of substituted salicylate monoanions was investigated as a function of ΔpKa, the difference between the pKa values of the hydrogen bond donor and acceptor, in both water and dimethyl sulfoxide. The dependence of log KHB upon ΔpKa is linear in both solvents, but is steeper in dimethyl sulfoxide (slope = 0.73) than in water (slope = 0.05). Thus, hydrogen bond strength can undergo substantially larger increases in nonaqueous media than aqueous solutions as the charge density on the donor or acceptor atom increases. These results support a general mechanism for enzymatic catalysis, in which hydrogen bonding to a substrate is strengthened as charge rearranges in going from the ground state to the transition state; the strengthening of the hydrogen bond would be greater in a nonaqueous enzymatic active site than in water, thus providing a rate enhancement for an enzymatic reaction relative to the solution reaction. We suggest that binding energy of an enzyme is used to fix the substrate in the low-dielectric active site, where the strengthening of the hydrogen bond in the course of a reaction is increased.

Unraveling the mechanism of the lactose permease of Escherichia coli

We studied the effect of pH on ligand binding in wild-type lactose permease or mutants in the four residues—Glu-269, Arg-302, His-322, and Glu-325—that are the key participants in H+ translocation and coupling between sugar and H+ translocation. Although wild-type permease or mutants in Glu-325 and Arg-302 exhibit marked decreases in affinity at alkaline pH, mutants in either His-322 or Glu-269 do not titrate. The results offer a mechanistic model for lactose/H+ symport. In the ground state, the permease is protonated, the H+ is shared between His-322 and Glu-269, Glu-325 is charge-paired with Arg-302, and substrate is bound with high affinity at the outside surface. Substrate binding induces a conformational change that leads to transfer of the H+ from His-322/Glu-269 to Glu-325 and reorientation of the binding site to the inner surface with a decrease in affinity. Glu-325 then is deprotonated on the inside because of rejuxtaposition with Arg-302. The His-322/Glu-269 complex then is reprotonated from the outside surface to reinitiate the cycle.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.

On the acidity and reactivity of HNO in aqueous solution and biological systems

The gas phase and aqueous thermochemistry and reactivity of nitroxyl (nitrosyl hydride, HNO) were elucidated with multiconfigurational self-consistent field and hybrid density functional theory calculations and continuum solvation methods. The pKa of HNO is predicted to be 7.2 ± 1.0, considerably different from the value of 4.7 reported from pulse radiolysis experiments. The ground-state triplet nature of NO− affects the rates of acid-base chemistry of the HNO/NO− couple. HNO is highly reactive toward dimerization and addition of soft nucleophiles but is predicted to undergo negligible hydration (Keq = 6.9 × 10−5). HNO is predicted to exist as a discrete species in solution and is a viable participant in the chemical biology of nitric oxide and derivatives.

An unusual pathway of excitation energy deactivation in carotenoids: Singlet-to-triplet conversion on an ultrafast timescale in a photosynthetic antenna

Carotenoids are important biomolecules that are ubiquitous in nature and find widespread application in medicine. In photosynthesis, they have a large role in light harvesting (LH) and photoprotection. They exert their LH function by donating their excited singlet state to nearby (bacterio)chlorophyll molecules. In photosynthetic bacteria, the efficiency of this energy transfer process can be as low as 30%. Here, we present evidence that an unusual pathway of excited state relaxation in carotenoids underlies this poor LH function, by which carotenoid triplet states are generated directly from carotenoid singlet states. This pathway, operative on a femtosecond and picosecond timescale, involves an intermediate state, which we identify as a new, hitherto uncharacterized carotenoid singlet excited state. In LH complex-bound carotenoids, this state is the precursor on the reaction pathway to the triplet state, whereas in extracted carotenoids in solution, this state returns to the singlet ground state without forming any triplets. We discuss the possible identity of this excited state and argue that fission of the singlet state into a pair of triplet states on individual carotenoid molecules constitutes the mechanism by which the triplets are generated. This is, to our knowledge, the first ever direct observation of a singlet-to-triplet conversion process on an ultrafast timescale in a photosynthetic antenna.