Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3β

Glycogen synthase kinase 3β (GSK-3β) targets specific signaling pathways in response to distinct upstream signals. We used structural and functional studies to dissect how an upstream phosphorylation step primes the Wnt signaling component β-catenin for phosphorylation by GSK-3β and how scaffolding interactions contribute to this reaction. Our crystal structure of GSK-3β bound to a phosphoprimed β-catenin peptide confirmed the expected binding mode of the phosphoprimed residue adjacent to the catalytic site. An aspartate phosphomimic in the priming site of β-catenin adopted an indistinguishable structure but reacted approximately 1000-fold slower than the native phosphoprimed substrate. This result suggests that substrate positioning alone is not sufficient for catalysis and that native phosphopriming interactions are necessary. We also obtained a structure of GSK-3β with an extended peptide from the scaffold protein Axin that bound with greater affinity than that of previously crystallized Axin fragments. This structure neither revealed additional contacts that produce the higher affinity nor explained how substrate interactions in the GSK-3β active site are modulated by remote Axin binding. Together, our findings suggest that phosphopriming and scaffolding produce small conformational changes or allosteric effects, not captured in the crystal structures, that activate GSK-3β and facilitate β-catenin phosphorylation. These results highlight limitations in our ability to predict catalytic activity from structure and have potential implications for the role of natural phosphomimic mutations in kinase regulation and phosphosite evolution.

Good Vibrations Report on the DNA Quadruplex Binding of an Excited State Amplified Ruthenium Polypyridyl IR Probe

The nitrile containing Ru(II)polypyridyl complex [Ru(phen)2(11,12-dCN-dppz)]2+ (1) is shown to act as a sensitive infrared probe of G-quadruplex (G4) structures. UV-visible absorption spectroscopy reveals enantiomer sensitive binding for the hybrid htel(K) and antiparallel htel(Na) G4s formed by the human telomer sequence d[AG3(TTAG3)3]. Time-resolved infrared (TRIR) of 1 upon 400 nm excitation indicates dominant interactions with the guanine bases in the case of Λ-1/htel(K), Δ-1/htel(K), and Λ-1/htel(Na) binding, whereas Δ-1/htel(Na) binding is associated with interactions with thymine and adenine bases in the loop. The intense nitrile transient at 2232 cm-1 undergoes a linear shift to lower frequency as the solution hydrogen bonding environment decreases in DMSO/water mixtures. This shift is used as a sensitive reporter of the nitrile environment within the binding pocket. The lifetime of 1 in D2O (ca. 100 ps) is found to increase upon DNA binding, and monitoring of the nitrile and ligand transients as well as the diagnostic DNA bleach bands shows that this increase is related to greater protection from the solvent environment. Molecular dynamics simulations together with binding energy calculations identify the most favorable binding site for each system, which are in excellent agreement with the observed TRIR solution study. This study shows the power of combining the environmental sensitivity of an infrared (IR) probe in its excited state with the TRIR DNA "site effect" to gain important information about the binding site of photoactive agents and points to the potential of such amplified IR probes as sensitive reporters of biological environments.

Bond Bundle Analysis of Ketosteroid Isomerase

Bond bundle analysis is used to investigate enzymatic catalysis in the ketosteroid isomerase (KSI) active site. We identify the unique bonding regions in five KSI systems, including those exposed to applied oriented electric fields and those with amino acid mutations, and calculate the precise redistribution of electron density and other regional properties that accompanies either enhancement or inhibition of KSI catalytic activity. We find that catalytic enhancement results from promoting both inter- and intra-molecular electron density redistribution, between bond bundles and bond wedges within the KSI-docked substrate molecule, in the forward direction of the catalyzed reaction. Though the redistribution applies to both types of perturbed systems and is thus suggestive of a general catalytic role, we observe that bond properties (e.g., volume vs energy vs electron count) can respond independently and disproportionately depending on the type of perturbation. We conclude that the resulting catalytic enhancement/inhibition proceeds via different mechanisms, where some bond properties are utilized more by one type of perturbation than the other. Additionally, we find that the correlations between bond wedge properties and catalyzed reaction barrier energies are additive to predict those of bond bundles and atomic basins, providing a rigorous grounding for connecting changes in local charge density to resulting shifts in reaction barrier energy.

Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra

The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH₂O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.

High throughput and quantitative enzymology in the genomic era

Accurate predictions from models based on physical principles are the ultimate metric of our biophysical understanding. Although there has been stunning progress toward structure prediction, quantitative prediction of enzyme function has remained challenging. Realizing this goal will require large numbers of quantitative measurements of rate and binding constants and the use of these ground-truth data sets to guide the development and testing of these quantitative models. Ground truth data more closely linked to the underlying physical forces are also desired. Here, we describe technological advances that enable both types of ground truth measurements. These advances allow classic models to be tested, provide novel mechanistic insights, and place us on the path toward a predictive understanding of enzyme structure and function.

Novel protein science enabled by total chemical synthesis

Chemical synthesis is a well-established method for the preparation in the research laboratory of multiple-tens-of-milligram amounts of correctly folded, high purity protein molecules. Chemically synthesized proteins enable a broad spectrum of novel protein science. Racemic mixtures consisting of d-protein and l-protein enantiomers facilitate crystallization and determination of protein structures by X-ray diffraction. d-Proteins enable the systematic development of unnatural mirror image protein molecules that bind with high affinity to natural protein targets. The d-protein form of a therapeutic target can also be used to screen natural product libraries to identify novel small molecule leads for drug development. Proteins with novel polypeptide chain topologies including branched, circular, linear-loop, and interpenetrating polypeptide chains can be constructed by chemical synthesis. Medicinal chemistry can be applied to optimize the properties of therapeutic protein molecules. Chemical synthesis has been used to redesign glycoproteins and for the a priori design and construction of covalently constrained novel protein scaffolds not found in nature. Versatile and precise labeling of protein molecules by chemical synthesis facilitates effective application of advanced physical methods including multidimensional nuclear magnetic resonance and time-resolved FTIR for the elucidation of protein structure-activity relationships. The chemistries used for total synthesis of proteins have been adapted to making artificial molecular devices and protein-inspired nanomolecular constructs. Research to develop mirror image life in the laboratory is in its very earliest stages, based on the total chemical synthesis of d-protein forms of polymerase enzymes.

Proton Transfer at the Interaction Interface of Graphene Oxide

Proton transfer plays a crucial role in a variety of biological phenomena. The transformation of nanomaterials in the environment and biology makes probing the potential proton transfer between nanomaterials and biomolecules a crucial issue, but it still remains a significant challenge. Here, we report proton transfer at the interface of graphene oxide (GO) by studying the GO-induced vibrational changes of interfacial water and carboxyl-terminated self-assembled monolayer (SAM) with surface-enhanced infrared absorption spectroscopy. In addition to simply acting as a macromolecular buffer in solution, the GO sheet behaves as a two-dimensional hydrogen-bonded exchangeable proton pool to dissociate and transfer protons at the interface with a suitable Brønsted base pair, which may bear a significant potential toxic origin for biological systems with proton-coupled reactions.

Rapid Sampling of Hydrogen Bond Networks for Computational Protein Design

Hydrogen bond networks play a critical role in determining the stability and specificity of biomolecular complexes, and the ability to design such networks is important for engineering novel structures, interactions, and enzymes. One key feature of hydrogen bond networks that makes them difficult to rationally engineer is that they are highly cooperative and are not energetically favorable until the hydrogen bonding potential has been satisfied for all buried polar groups in the network. Existing computational methods for protein design are ill-equipped for creating these highly cooperative networks because they rely on energy functions and sampling strategies that are focused on pairwise interactions. To enable the design of complex hydrogen bond networks, we have developed a new sampling protocol in the molecular modeling program Rosetta that explicitly searches for sets of amino acid mutations that can form self-contained hydrogen bond networks. For a given set of designable residues, the protocol often identifies many alternative sets of mutations/networks, and we show that it can readily be applied to large sets of residues at protein-protein interfaces or in the interior of proteins. The protocol builds on a recently developed method in Rosetta for designing hydrogen bond networks that has been experimentally validated for small symmetric systems but was not extensible to many larger protein structures and complexes. The sampling protocol we describe here not only recapitulates previously validated designs with performance improvements but also yields viable hydrogen bond networks for cases where the previous method fails, such as the design of large, asymmetric interfaces relevant to engineering protein-based therapeutics.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

Proton transfer aiding phase transitions in oxalic acid dihydrate under pressure

Oxalic acid dihydrate, an important molecular solid in crystal chemistry, ecology and physiology, has been studied for nearly 100 years now. The most debated issues regarding its proton dynamics have arisen due to an unusually short hydrogen bond between the acid and water molecules. Using combined in situ spectroscopic studies and first-principles simulations at high pressures, we show that the structural modification associated with this hydrogen bond is much more significant than ever assumed. Initially, under pressure, proton migration takes place along this strong hydrogen bond at a very low pressure of 2 GPa. This results in the protonation of water with systematic formation of dianionic oxalate and hydronium ion motifs, thus reversing the hydrogen bond hierarchy in the high pressure phase II. The resulting hydrogen bond between a hydronium ion and a carboxylic group shows remarkable strengthening under pressure, even in the pure ionic phase III. The loss of cooperativity of hydrogen bonds leads to another phase transition at ∼ 9 GPa through reorientation of other hydrogen bonds. The high pressure phase IV is stabilized by a strong hydrogen bond between the dominant CO2 and H2O groups of oxalate and hydronium ions, respectively. These findings suggest that oxalate systems may provide useful insights into proton transfer reactions and assembly of simple molecules under extreme conditions.

Evaluation of the Catalytic Contribution from a Positioned General Base in Ketosteroid Isomerase

Proton transfer reactions are ubiquitous in enzymes and utilize active site residues as general acids and bases. Crystal structures and site-directed mutagenesis are routinely used to identify these residues, but assessment of their catalytic contribution remains a major challenge. In principle, effective molarity measurements, in which exogenous acids/bases rescue the reaction in mutants lacking these residues, can estimate these catalytic contributions. However, these exogenous moieties can be restricted in reactivity by steric hindrance or enhanced by binding interactions with nearby residues, thereby resulting in over- or underestimation of the catalytic contribution, respectively. With these challenges in mind, we investigated the catalytic contribution of an aspartate general base in ketosteroid isomerase (KSI) by exogenous rescue. In addition to removing the general base, we systematically mutated nearby residues and probed each mutant with a series of carboxylate bases of similar pKa but varying size. Our results underscore the need for extensive and multifaceted variation to assess and minimize steric and positioning effects and determine effective molarities that estimate catalytic contributions. We obtained consensus effective molarities of ∼5 × 10(4) M for KSI from Comamonas testosteroni (tKSI) and ∼10(3) M for KSI from Pseudomonas putida (pKSI). An X-ray crystal structure of a tKSI general base mutant showed no additional structural rearrangements, and double mutant cycles revealed similar contributions from an oxyanion hole mutation in the wild-type and base-rescued reactions, providing no indication of mutational effects extending beyond the general base site. Thus, the high effective molarities suggest a large catalytic contribution associated with the general base. A significant portion of this effect presumably arises from positioning of the base, but its large magnitude suggests the involvement of additional catalytic mechanisms as well.

Role of conserved Met112 residue in the catalytic activity and stability of ketosteroid isomerase

Ketosteroid isomerase (3-oxosteroid Δ(5)-Δ(4)-isomerase, KSI) from Pseudomonas putida catalyzes allylic rearrangement of the 5,6-double bond of Δ(5)-3-ketosteroid to 4,5-position by stereospecific intramolecular transfer of a proton. The active site of KSI is formed by several hydrophobic residues and three catalytic residues (Tyr14, Asp38, and Asp99). In this study, we investigated the role of a hydrophobic Met112 residue near the active site in the catalysis, steroid binding, and stability of KSI. Replacing Met112 with alanine (yields M112A) or leucine (M112L) decreased the kcat by 20- and 4-fold, respectively. Compared with the wild type (WT), M112A and M112L KSIs showed increased KD values for equilenin, an intermediate analogue; these changes suggest that loss of packing at position 112 might lead to unfavorable steroid binding, thereby resulting in decreased catalytic activity. Furthermore, M112A and M112L mutations reduced melting temperature (Tm) by 6.4°C and 2.5°C, respectively. These changes suggest that favorable packing in the core is important for the maintenance of stability in KSI. The M112K mutation decreased kcat by 2000-fold, compared with the WT. In M112K KSI structure, a new salt bridge was formed between Asp38 and Lys112. This bridge could change the electrostatic potential of Asp38, and thereby contribute to the decreased catalytic activity. The M112K mutation also decreased the stability by reducing Tm by 4.1°C. Our data suggest that the Met112 residue may contribute to the catalytic activity and stability of KSI by providing favorable hydrophobic environments and compact packing in the catalytic core.

Effect of the Hydrogen Bond in Photoinduced Water Dissociation: A Double-Edged Sword

Photoinduced water dissociation on rutile-TiO2 was investigated using various methods. Experimental results reveal that the water dissociation occurs via transferring an H atom to a bridge bonded oxygen site and ejecting an OH radical to the gas phase during irradiation. The reaction is strongly suppressed as the water coverage increases. Further scanning tunneling microscopy study demonstrates that hydrogen bonds between water molecules have a dramatic effect on the reaction. Interestingly, a single hydrogen bond in water dimer enhances the water dissociation reaction, while one-dimensional hydrogen bonds in water chains inhibit the reaction. Density functional theory calculations indicate that the effect of hydrogen bonds on the OH dissociation energy is likely the origin of this remarkable behavior. The results suggest that avoiding a strong hydrogen bond network between water molecules is crucial for water splitting.

Dissecting Proton Delocalization in an Enzyme's Hydrogen Bond Network with Unnatural Amino Acids

Extended hydrogen bond networks are a common structural motif of enzymes. A recent analysis proposed quantum delocalization of protons as a feature present in the hydrogen bond network spanning a triad of tyrosines (Y(16), Y(32), and Y(57)) in the active site of ketosteroid isomerase (KSI), contributing to its unusual acidity and large isotope shift. In this study, we utilized amber suppression to substitute each tyrosine residue with 3-chlorotyrosine to test the delocalization model and the proton affinity balance in the triad. X-ray crystal structures of each variant demonstrated that the structure, notably the O-O distances within the triad, was unaffected by 3-chlorotyrosine substitutions. The changes in the cluster's acidity and the acidity's isotope dependence in these variants were assessed via UV-vis spectroscopy and the proton sharing pattern among individual residues with (13)C nuclear magnetic resonance. Our data show pKa detuning at each triad residue alters the proton delocalization behavior in the H-bond network. The extra stabilization energy necessary for the unusual acidity mainly comes from the strong interactions between Y(57) and Y(16). This is further enabled by Y(32), which maintains the right geometry and matched proton affinity in the triad. This study provides a rich picture of the energetics of the hydrogen bond network in enzymes for further model refinement.

Molecular quantum mechanical gradients within the polarizable embedding approach--application to the internal vibrational Stark shift of acetophenone

We present an implementation of analytical quantum mechanical molecular gradients within the polarizable embedding (PE) model to allow for efficient geometry optimizations and vibrational analysis of molecules embedded in large, geometrically frozen environments. We consider a variational ansatz for the quantum region, covering (multiconfigurational) self-consistent-field and Kohn-Sham density functional theory. As the first application of the implementation, we consider the internal vibrational Stark effect of the C=O group of acetophenone in different solvents and derive its vibrational linear Stark tuning rate using harmonic frequencies calculated from analytical gradients and computed local electric fields. Comparisons to PE calculations employing an enlarged quantum region as well as to a non-polarizable embedding scheme show that the inclusion of mutual polarization between acetophenone and water is essential in order to capture the structural modifications and the associated frequency shifts observed in water. For more apolar solvents, a proper description of dispersion and exchange-repulsion becomes increasingly important, and the quality of the optimized structures relies to a larger extent on the quality of the Lennard-Jones parameters.

Accurate ab initio calculations of O-HO and O-H(-)O proton chemical shifts: towards elucidation of the nature of the hydrogen bond and prediction of hydrogen bond distances

The inability to determine precisely the location of labile protons in X-ray molecular structures has been a key barrier to progress in many areas of molecular sciences. We report an approach for predicting hydrogen bond distances beyond the limits of X-ray crystallography based on accurate ab initio calculations of O-HO proton chemical shifts, using a combination of DFT and contactor-like polarizable continuum model (PCM). Very good linear correlation between experimental and computed (at the GIAO/B3LYP/6-311++G(2d,p) level of theory) chemical shifts were obtained with a large set of 43 compounds in CHCl3 exhibiting intramolecular O-HO and intermolecular and intramolecular ionic O-H(-)O hydrogen bonds. The calculated OH chemical shifts exhibit a strong linear dependence on the computed (O)HO hydrogen bond length, in the region of 1.24 to 1.85 Å, of -19.8 ppm Å(-1) and -20.49 ppm Å(-1) with optimization of the structures at the M06-2X/6-31+G(d) and B3LYP/6-31+G(d) level of theory, respectively. A Natural Bond Orbitals (NBO) analysis demonstrates a very good linear correlation between the calculated (1)H chemical shifts and (i) the second-order perturbation stabilization energies, corresponding to charge transfer between the oxygen lone pairs and σ antibonding orbital and (ii) Wiberg bond order of the O-HO and O-H(-)O hydrogen bond. Accurate ab initio calculations of O-HO and O-H(-)O (1)H chemical shifts can provide improved structural and electronic description of hydrogen bonding and a highly accurate measure of distances of short and strong hydrogen bonds.

Live-cell vibrational imaging of choline metabolites by stimulated Raman scattering coupled with isotope-based metabolic labeling

Choline is a small molecule that occupies a key position in the biochemistry of all living organisms. Recent studies have strongly implicated choline metabolites in cancer, atherosclerosis and nervous system development. To detect choline and its metabolites, existing physical methods such as magnetic resonance spectroscopy and positron emission tomography are often limited by the poor spatial resolution and substantial radiation dose. Fluorescence imaging, although with submicrometer resolution, requires introduction of bulky fluorophores and thus is difficult in labeling the small choline molecule. By combining the emerging bond-selective stimulated Raman scattering microscopy with metabolic incorporation of deuterated choline, herein we have achieved high resolution imaging of choline-containing metabolites in living mammalian cell lines, primary hippocampal neurons and the multicellular organism C. elegans. Different subcellular distributions of choline metabolites are observed between cancer cells and non-cancer cells, which may reveal a functional difference in the choline metabolism and lipid-mediated signaling events. In neurons, choline incorporation is visualized within both soma and neurites, where choline metabolites are more evenly distributed compared to proteins. Furthermore, choline localization is also observed in the pharynx region of C. elegans larvae, consistent with its organogenesis mechanism. These applications demonstrate the potential of isotope-based stimulated Raman scattering microscopy for future choline-related disease detection and development monitoring in vivo.

Joint neutron crystallographic and NMR solution studies of Tyr residue ionization and hydrogen bonding: Implications for enzyme-mediated proton transfer

Human carbonic anhydrase II (HCA II) uses a Zn-bound OH(-)/H2O mechanism to catalyze the reversible hydration of CO2. This catalysis also involves a separate proton transfer step, mediated by an ordered solvent network coordinated by hydrophilic residues. One of these residues, Tyr7, was previously shown to be deprotonated in the neutron crystal structure at pH 10. This observation indicated that Tyr7 has a perturbed pKa compared with free tyrosine. To further probe the pKa of this residue, NMR spectroscopic measurements of [(13)C]Tyr-labeled holo HCA II (with active-site Zn present) were preformed to titrate all Tyr residues between pH 5.4-11.0. In addition, neutron studies of apo HCA II (with Zn removed from the active site) at pH 7.5 and holo HCA II at pH 6 were conducted. This detailed interrogation of tyrosines in HCA II by NMR and neutron crystallography revealed a significantly lowered pKa of Tyr7 and how pH and Tyr proximity to Zn affect hydrogen-bonding interactions.

Site-specific protonation kinetics of acidic side chains in proteins determined by pH-dependent carboxyl (13)C NMR relaxation

Proton-transfer dynamics plays a critical role in many biochemical processes, such as proton pumping across membranes and enzyme catalysis. The large majority of enzymes utilize acid-base catalysis and proton-transfer mechanisms, where the rates of proton transfer can be rate limiting for the overall reaction. However, measurement of proton-exchange kinetics for individual side-chain carboxyl groups in proteins has been achieved in only a handful of cases, which typically have involved comparative analysis of mutant proteins in the context of reaction network modeling. Here we describe an approach to determine site-specific protonation and deprotonation rate constants (kon and koff, respectively) of carboxyl side chains, based on (13)C NMR relaxation measurements as a function of pH. We validated the method using an extensively studied model system, the B1 domain of protein G, for which we measured rate constants koff in the range (0.1-3) × 10(6) s(-1) and kon in the range (0.6-300) × 10(9) M(-1) s(-1), which correspond to acid-base equilibrium dissociation constants (Ka) in excellent agreement with previous results determined by chemical shift titrations. Our results further reveal a linear free-energy relationship between log kon and pKa, which provides information on the free-energy landscape of the protonation reaction, showing that the variability among residues in these parameters arises primarily from the extent of charge stabilization of the deprotonated state by the protein environment. We find that side-chain carboxyls with extreme values of koff or kon are involved in hydrogen bonding, thus providing a mechanistic explanation for the observed stabilization of the protonated or deprotonated state.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

Revealing the functional states in the active site of BLUF photoreceptors from electrochromic shift calculations

Photoexcitation with blue light of the flavin chromophore in BLUF photoreceptors induces a switch into a metastable signaling state that is characterized by a red-shifted absorption maximum. The red shift is due to a rearrangement in the hydrogen bond pattern around Gln63 located in the immediate proximity of the isoalloxazine ring system of the chromophore. There is a long-lasting controversy between two structural models, named Q63A and Q63J in the literature, on the local conformation of the residues Gln63 and Tyr21 in the dark state of the photoreceptor. As regards the mechanistic details of the light-activation mechanism, rotation of Gln63 is opposed by tautomerism in the Q63A and Q63J models, respectively. We provide a structure-based simulation of electrochromic shifts of the flavin chromophore in the wild type and in various site-directed mutants. The excellent overall agreement between experimental and computed data allows us to evaluate the two structural models. Compelling evidence is obtained that the Q63A model is incorrect, whereas the Q63J is fully consistent with the present computations. Finally, we confirm independently that a keto-enol tautomerization of the glutamine at position 63, which was proposed as molecular mechanism for the transition between the dark and the light-adapted state, explains the measured 10 to 15 nm red shift in flavin absorption between these two states of the protein. We believe that the accurateness of our results provides evidence that the BLUF photoreceptors absorption is fine-tuned through electrostatic interactions between the chromophore and the protein matrix, and finally that the simplicity of our theoretical model is advantageous as regards easy reproducibility and further extensions.

Ground-state proton transfer kinetics in green fluorescent protein

Proton transfer plays an important role in the optical properties of green fluorescent protein (GFP). While much is known about excited-state proton transfer reactions (ESPT) in GFP occurring on ultrafast time scales, comparatively little is understood about the factors governing the rates and pathways of ground-state proton transfer. We have utilized a specific isotopic labeling strategy in combination with one-dimensional ¹³C nuclear magnetic resonance (NMR) spectroscopy to install and monitor a ¹³C directly adjacent to the GFP chromophore ionization site. The chemical shift of this probe is highly sensitive to the protonation state of the chromophore, and the resulting spectra reflect the thermodynamics and kinetics of the proton transfer in the NMR line shapes. This information is complemented by time-resolved NMR, fluorescence correlation spectroscopy, and steady-state absorbance and fluorescence measurements to provide a picture of chromophore ionization reactions spanning a wide time domain. Our findings indicate that proton transfer in GFP is described well by a two-site model in which the chromophore is energetically coupled to a secondary site, likely the terminal proton acceptor of ESPT, Glu222. Additionally, experiments on a selection of GFP circular permutants suggest an important role played by the structural dynamics of the seventh β-strand in gating proton transfer from bulk solution to the buried chromophore.

Putative hydrogen bond to tyrosine M208 in photosynthetic reaction centers from Rhodobacter capsulatus significantly slows primary charge separation

Slow, ∼50 ps, P* → P⁺H_A^– electron transfer is observed in Rhodobacter capsulatus reaction centers (RCs) bearing the native Tyr residue at M208 and the single amino acid change of isoleucine at M204 to glutamic acid. The P* decay kinetics are unusually homogeneous (single exponential) at room temperature. Comparative solid-state NMR of [4′-¹³C]Tyr labeled wild-type and M204E RCs show that the chemical shift of Tyr M208 is significantly altered in the M204E mutant and in a manner consistent with formation of a hydrogen bond to the Tyr M208 hydroxyl group. Models based on RC crystal structure coordinates indicate that if such a hydrogen bond is formed between the Glu at M204 and the M208 Tyr hydroxyl group, the −OH would be oriented in a fashion expected (based on the calculations by Alden et al., J. Phys. Chem.1996, 100, 16761–16770) to destabilize P⁺B_A^– in free energy. Alteration of the environment of Tyr M208 and B_A by Glu M204 via this putative hydrogen bond has a powerful influence on primary charge separation.

Uncovering the determinants of a highly perturbed tyrosine pKa in the active site of ketosteroid isomerase

Within the idiosyncratic enzyme active-site environment, side chain and ligand pKa values can be profoundly perturbed relative to their values in aqueous solution. Whereas structural inspection of systems has often attributed perturbed pKa values to dominant contributions from placement near charged groups or within hydrophobic pockets, Tyr57 of a Pseudomonas putida ketosteroid isomerase (KSI) mutant, suggested to have a pKa perturbed by nearly 4 units to 6.3, is situated within a solvent-exposed active site devoid of cationic side chains, metal ions, or cofactors. Extensive comparisons among 45 variants with mutations in and around the KSI active site, along with protein semisynthesis, (13)C NMR spectroscopy, absorbance spectroscopy, and X-ray crystallography, was used to unravel the basis for this perturbed Tyr pKa. The results suggest that the origin of large energetic perturbations are more complex than suggested by visual inspection. For example, the introduction of positively charged residues near Tyr57 raises its pKa rather than lowers it; this effect, and part of the increase in the Tyr pKa from the introduction of nearby anionic groups, arises from accompanying active-site structural rearrangements. Other mutations with large effects also cause structural perturbations or appear to displace a structured water molecule that is part of a stabilizing hydrogen-bond network. Our results lead to a model in which three hydrogen bonds are donated to the stabilized ionized Tyr, with these hydrogen-bond donors, two Tyr side chains, and a water molecule positioned by other side chains and by a water-mediated hydrogen-bond network. These results support the notion that large energetic effects are often the consequence of multiple stabilizing interactions rather than a single dominant interaction. Most generally, this work provides a case study for how extensive and comprehensive comparisons via site-directed mutagenesis in a tight feedback loop with structural analysis can greatly facilitate our understanding of enzyme active-site energetics. The extensive data set provided may also be a valuable resource for those wishing to extensively test computational approaches for determining enzymatic pKa values and energetic effects.

Thermodynamic framework for identifying free energy inventories of enzyme catalytic cycles

Pauling's suggestion that enzymes are complementary in structure to the activated complexes of the reactions they catalyze has provided the conceptual basis to explain how enzymes obtain their fantastic catalytic prowess, and has served as a guiding principle in drug design for over 50 y. However, this model by itself fails to predict the magnitude of enzymes' rate accelerations. We construct a thermodynamic framework that begins with the classic concept of differential binding but invokes additional terms that are needed to account for subtle effects in the catalytic cycle's proton inventory. Although the model presented can be applied generally, this analysis focuses on ketosteroid isomerase (KSI) as an example, where recent experiments along with a large body of kinetic and thermodynamic data have provided strong support for the noncanonical thermodynamic contribution described. The resulting analysis precisely predicts the free energy barrier of KSI's reaction as determined from transition-state theory using only empirical thermodynamic data. This agreement is suggestive that a complete free energy inventory of the KSI catalytic cycle has been identified.