Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site

Antifouling activities of proteinase K and α-amylase enzymes: Laboratory bioassays and in silico analysis

The application of enzymes as antifoulants is one of the environment-friendly strategies in biofouling management. In this study, antifouling activities of commercially available proteinase K and α-amylase enzymes were evaluated using barnacle larva and biofilm-forming bacteria as test organisms. The enzymes were also tested against barnacle cement protein through in silico analysis. The results showed that both enzymes inhibited the attachment of bacteria and settlement of barnacle larvae on the test surface. The lowest minimum inhibitory concentration of 0.312 mg ml-1 was exhibited by proteinase K against biofilm-forming bacteria. The calculated LC50 values for proteinase K and α-amylase against the barnacle nauplii were 91.8 and 230.96 mg ml-1 respectively. While α-amylase showed higher antibiofilm activity, proteinase K exhibited higher anti-larval settlement activity. Similarly, in silico analysis of the enzymes revealed promising anti-settlement activity, as the enzymes showed good binding scores with barnacle cement protein. Overall, the results suggested that the enzymes proteinase K and α-amylase could be used in antifouling coatings to reduce the settlement of biofouling on artificial materials in the marine environment.

Evaluating the Impact of the Tyr158 pKa on the Mechanism and Inhibition of InhA, the Enoyl-ACP Reductase from Mycobacterium tuberculosis

InhA, the Mycobacterium tuberculosis enoyl-ACP reductase, is a target for the tuberculosis (TB) drug isoniazid (INH). InhA inhibitors that do not require KatG activation avoid the most common mechanism of INH resistance, and there are continuing efforts to fully elucidate the enzyme mechanism to drive inhibitor discovery. InhA is a member of the short-chain dehydrogenase/reductase superfamily characterized by a conserved active site Tyr, Y158 in InhA. To explore the role of Y158 in the InhA mechanism, this residue has been replaced by fluoroTyr residues that increase the acidity of Y158 up to ∼3200-fold. Replacement of Y158 with 3-fluoroTyr (3-FY) and 3,5-difluoroTyr (3,5-F2Y) has no effect on kcatapp/KMapp nor on the binding of inhibitors to the open form of the enzyme (Kiapp), whereas both kcatapp/KMapp and Kiapp are altered by seven-fold for the 2,3,5-trifluoroTyr variant (2,3,5-F3Y158 InhA). 19F NMR spectroscopy suggests that 2,3,5-F3Y158 is ionized at neutral pH indicating that neither the acidity nor ionization state of residue 158 has a major impact on catalysis or on the binding of substrate-like inhibitors. In contrast, Ki*app is decreased 6- and 35-fold for the binding of the slow-onset inhibitor PT504 to 3,5-F2Y158 and 2,3,5-F3Y158 InhA, respectively, indicating that Y158 stabilizes the closed form of the enzyme adopted by EI*. The residence time of PT504 is reduced ∼four-fold for 2,3,5-F3Y158 InhA compared to wild-type, and thus, the hydrogen bonding interaction of the inhibitor with Y158 is an important factor in the design of InhA inhibitors with increased residence times on the enzyme.

The Role of H-Bonds in the Excited-State Properties of Multichromophoric Systems: Static and Dynamic Aspects

Given their importance, hydrogen bonds (H-bonds) have been the subject of intense investigation since their discovery. Indeed, H-bonds play a fundamental role in determining the structure, the electronic properties, and the dynamics of complex systems, including biologically relevant materials such as DNA and proteins. While H-bonds have been largely investigated for systems in their electronic ground state, fewer studies have focused on how the presence of H-bonds could affect the static and dynamic properties of electronic excited states. This review presents an overview of the more relevant progress in studying the role of H-bond interactions in modulating excited-state features in multichromophoric biomimetic complex systems. The most promising spectroscopic techniques that can be used for investigating the H-bond effects in excited states and for characterizing the ultrafast processes associated with their dynamics are briefly summarized. Then, experimental insights into the modulation of the electronic properties resulting from the presence of H-bond interactions are provided, and the role of the H-bond in tuning the excited-state dynamics and the related photophysical processes is discussed.

Number of Hydrogen Bonds per Unit Solvent Accessible Surface Area: A Descriptor of Functional States of Proteins

Proteins function close to native and near-native conformations. These states are evolutionarily selected to ensure the effect of mutations is minimized. The structural organization of a protein is hierarchical and modular, which reduces the dimensionality of the configurational space of the native states. Thus, finding appropriate descriptors that define the native state among all possible states of a protein is a problem of immense interest. The present study explores the correlation between solvent accessible surface areas (SASAs) and different intraprotein as well as protein-water hydrogen bonds of 55 single-chain globular proteins from four different structural classes (all α, all β, α+β, and α/β), 16 multichain proteins, and 4 macromolecular complexes. A systematic analysis of the solvent accessible surface area and intraprotein and protein-water hydrogen bonds suggests a linear relationship between SASAs and hydrogen bonds. The number of protein-water hydrogen bonds per unit SASA ranges from 3 to 4 for all the different structural protein classes. In contrast, the number of intramolecular hydrogen bonds per unit SASA, including the mainchain-mainchain, mainchain-sidechain, and sidechain-sidechain, varies between 0.75 to 2. The solvation free energy of a protein linearly decreases with SASA. Our study also shows that the solvation free energy/SASA varies from -75 to -105 kJ mol^-1 nm^-2 across all the native states studied here. The number conservancy of intraprotein hydrogen bonds per unit SASA possibly imparts structural stability to the native structure. On the other hand, 3-4 protein-water hydrogen bonds per unit SASA are possibly required to maintain a balance between the solubility and functionality of the native states. This study provides a basis for synthetic biologists to design new folds with improved functionality.

Metabolomics analysis of milk thistle lipids to identify drought-tolerant genes

Milk thistle is an oil and medicinal crop known as an alternative oil crop with a high level of unsaturated fatty acids, which makes it a favorable edible oil for use in food production. To evaluate the importance of Milk thistle lipids in drought tolerance, an experiment was performed in field conditions under three different water deficit levels (Field capacity (FC), 70% FC and 40% FC). After harvesting seeds of the plant, their oily and methanolic extracts were isolated, and subsequently, types and amounts of lipids were measured using GC-MS. Genes and enzymes engaged in biosynthesizing of these lipids were identified and their expression in Arabidopsis was investigated under similar conditions. The results showed that content of almost all measured lipids of milk thistle decreased under severe drought stress, but genes (belonged to Arabidopsis), which were involved in their biosynthetic pathway showed different expression patterns. Genes biosynthesizing lipids, which had significant amounts were selected and their gene and metabolic network were established. Two networks were correlated, and for each pathway, their lipids and respective biosynthesizing genes were grouped together. Four up-regulated genes including PXG3, LOX2, CYP710A1, PAL and 4 down-regulated genes including FATA2, CYP86A1, LACS3, PLA2-ALPHA were selected. The expression of these eight genes in milk thistle was similar to Arabidopsis under drought stress. Thus, PXG3, PAL, LOX2 and CYP86A1 genes that increased expression were selected for protein analysis. Due to the lack of protein structure of these genes in the milk thistle, modeling homology was performed for them. The results of molecular docking showed that the four proteins CYP86A1, LOX2, PAL and PXG3 bind to ligands HEM, 11O, ACT and LIG, respectively. HEM ligand was involved in production of secondary metabolites and dehydration tolerance, and HEM binding site remained conserved in various plants. CA ligands were involved in synthesis of cuticles and waxes. Overall, this study confirmed the importance of lipids in drought stress tolerance in milk thistle.

Mutation of the conserved Asp-Asp pair impairs the structure, function, and inhibition of CTX-M Class A β-lactamase

The Asp233-Asp246 pair is highly conserved in Class A β-lactamases, which hydrolyze β-lactam antibiotics. Here, we characterize its function using CTX-M-14 β-lactamase. The D233N mutant displayed decreased activity that is substrate-dependent, with reductions in k_cat /K_m ranging from 20% for nitrocefin to 6-fold for cefotaxime. In comparison, the mutation reduced the binding of a known reversible inhibitor by 10-fold. The mutant structures showed movement of the 213-219 loop and the loss of the Thr216-Thr235 hydrogen bond, which was restored by inhibitor binding. Mutagenesis of Thr216 further highlighted its contribution to CTX-M activity. These results demonstrate the importance of the aspartate pair to CTX-M hydrolysis of substrates with bulky side chains, while suggesting increased protein flexibility as a means to evolve drug resistance.

Differences in Gluco and Galacto Substrate-Binding Interactions in a Dual 6Pβ-Glucosidase/6Pβ-Galactosidase Glycoside Hydrolase 1 Enzyme from Bacillus licheniformis

Bacterial glycoside hydrolase 1 (GH1) enzymes with 6-phospho-β-galactosidase and 6-phospho-β-glucosidase activities have the important task of releasing phosphorylated and nonphosphorylated monosaccharides into the cytoplasm. Curiously, dual 6-phospho-β-galactosidase/6-phospho-β-glucosidase (dual-phospho) enzymes have broad specificity and are able to hydrolyze galacto- and gluco-derived substrates. This study investigates the structure and substrate specificity of a GH family 1 enzyme from Bacillus licheniformis, hereafter known as BlBglC. The enzyme structure has been solved, and sequence analysis, molecular dynamics simulations, and binding free energy calculations offered evidence of dual-phospho activity. Both test ligands p-nitrophenyl-β-d-galactoside-6-phosphate (PNP6Pgal) and p-nitrophenyl-β-d-glucoside-6-phosphate (PNP6Pglc) demonstrated strong binding to BlBglC although the pose and interactions of the PNP6Pglc triplicates were slightly more consistent. Interestingly, known specificity-inducing residues, Gln23 and Trp433, bind strongly to the ligand O3 hydroxyl group in the PNP6Pgal-BlBglC complex and to the ligand O4 hydroxyl group in the PNP6Pglc-BlBglC complex. Additionally, the BlBglC-His124 residue is a major contributor of hydrogen bonds to the PNP6Pgal O3 hydroxyl group but does not form any hydrogen bonds with PNP6Pglc. On the other hand, BlBglC residues Tyr173, Tyr301, Gln302, and Thr321 form hydrogen bonds with PNP6Pglc but not PNP6Pgal. These findings provide important details of the broad specificity of dual-phospho activity GH1 enzymes.

Combining flavin photocatalysis with parallel synthesis: a general platform to optimize peptides with non-proteinogenic amino acids

Most peptide drugs contain non-proteinogenic amino acids (NPAAs), born out through extensive structure-activity relationship (SAR) studies using solid-phase peptide synthesis (SPPS). Synthetically laborious and expensive to manufacture, NPAAs also can have poor coupling efficiencies allowing only a small fraction to be sampled by conventional SPPS. To gain general access to NPAA-containing peptides, we developed a first-generation platform that merges contemporary flavin photocatalysis with parallel synthesis to simultaneously make, purify, quantify, and even test up to 96 single-NPAA peptide variants via the unique combination of boronic acids and a dehydroalanine residue in a peptide. We showcase the power of our newly minted platform to introduce NPAAs of diverse chemotypes-aliphatic, aromatic, heteroaromatic-directly into peptides, including 15 entirely new residues, and to evolve a simple proteinogenic peptide into an unnatural inhibitor of thrombin by non-classical peptide SAR.

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.

How water-mediated hydrogen bonds affect chlorophyll a/b selectivity in Water-Soluble Chlorophyll Protein

The Water-Soluble Chlorophyll Protein (WSCP) of Brassicaceae is a remarkably stable tetrapyrrole-binding protein that, by virtue of its simple design, is an exceptional model to investigate the interactions taking place between pigments and their protein scaffold and how they affect the photophysical properties and the functionality of the complexes. We investigated variants of WSCP from Lepidium virginicum (Lv) and Brassica oleracea (Bo), reconstituted with Chlorophyll (Chl) b, to determine the mechanisms by which the different Chl binding sites control their Chl a/b specificities. A combined Raman and crystallographic investigation has been employed, aimed to characterize in detail the hydrogen-bond network involving the formyl group of Chl b. The study revealed a variable degree of conformational freedom of the hydrogen bond networks among the WSCP variants, and an unexpected mixed presence of hydrogen-bonded and not hydrogen-bonded Chls b in the case of the L91P mutant of Lv WSCP. These findings helped to refine the description of the mechanisms underlying the different Chl a/b specificities of WSCP versions, highlighting the importance of the structural rigidity of the Chl binding site in the vicinity of the Chl b formyl group in granting a strong selectivity to binding sites.

Hydrogen Bonds: Simple after All?

Hydrogen bonds play integral roles in biological structure, function, and conformational dynamics and are fundamental to life as it has evolved on Earth. However, our understanding of these fundamental and ubiquitous interactions has seemed fractured and incomplete, and it has been difficult to extract generalities and principles about hydrogen bonds despite thousands of papers published on this topic, perhaps in part because of the expanse of this subject and the density of studies. Fortunately, recent hydrogen bond proposals, discussions, and debates have stimulated new tests and models and have led to a remarkably simple picture of the structure of hydrogen bonds. This knowledge also provides clarity concerning hydrogen bond energetics, limiting and simplifying the factors that need be considered. Herein we recount the advances that have led to this simpler view of hydrogen bond structure, dynamics, and energetics. A quantitative predictive model for hydrogen bond length can now be broadly and deeply applied to evaluate current proposals and to uncover structural features of proteins, their conformational restraints, and their correlated motions. In contrast, a quantitative energetic description of molecular recognition and catalysis by proteins remains an important ongoing challenge, although our improved understanding of hydrogen bonds may aid in testing predictions from current and future models. We close by codifying our current state of understanding into five "Rules for Hydrogen Bonding" that may provide a foundation for understanding and teaching about these vital interactions and for building toward a deeper understanding of hydrogen bond energetics.

Electric Fields and Enzyme Catalysis

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

A Critical Test of the Electrostatic Contribution to Catalysis with Noncanonical Amino Acids in Ketosteroid Isomerase

The vibrational Stark effect (VSE) has been used to measure the electric field in the active site of ketosteroid isomerase (KSI). These measured fields correlate with ΔG(⧧) in a series of conventional mutants, yielding an estimate for the electrostatic contribution to catalysis (Fried et al. Science 2014, 346, 1510-1513). In this work we test this result with much more conservative variants in which individual Tyr residues in the active site are replaced by 3-chlorotyrosine via amber suppression. The electric fields sensed at the position of the carbonyl bond involved in charge displacement during catalysis were characterized using the VSE, where the field sensitivity has been calibrated by vibrational Stark spectroscopy, solvatochromism, and MD simulations. A linear relationship is observed between the electric field and ΔG(⧧) that interpolates between wild-type and more drastic conventional mutations, reinforcing the evaluation of the electrostatic contribution to catalysis in KSI. A simplified model and calculation are developed to estimate changes in the electric field accompanying changes in the extended hydrogen-bond network in the active site. The results are consistent with a model in which the O-H group of a key active site tyrosine functions by imposing a static electrostatic potential onto the carbonyl bond. The model suggests that the contribution to catalysis from the active site hydrogen bonds is of similar weight to the distal interactions from the rest of the protein. A similar linear correlation was also observed between the proton affinity of KSI's active site and the catalytic rate, suggesting a direct connection between the strength of the H-bond and the electric field it exerts.

Biophysical Characterization of Fluorotyrosine Probes Site-Specifically Incorporated into Enzymes: E. coli Ribonucleotide Reductase As an Example

Fluorinated tyrosines (F_nY’s, n = 2 and 3) have been site-specifically incorporated into E. coli class Ia ribonucleotide reductase (RNR) using the recently evolved M. jannaschii Y-tRNA synthetase/tRNA pair. Class Ia RNRs require four redox active Y’s, a stable Y radical (Y·) in the β subunit (position 122 in E. coli), and three transiently oxidized Y’s (356 in β and 731 and 730 in α) to initiate the radical-dependent nucleotide reduction process. F_nY (3,5; 2,3; 2,3,5; and 2,3,6) incorporation in place of Y₁₂₂-β and the X-ray structures of each resulting β with a diferric cluster are reported and compared with wt-β2 crystallized under the same conditions. The essential diferric-F_nY· cofactor is self-assembled from apo F_nY-β2, Fe²⁺, and O₂ to produce ∼1 Y·/β2 and ∼3 Fe³⁺/β2. The F_nY· are stable and active in nucleotide reduction with activities that vary from 5% to 85% that of wt-β2. Each F_nY·-β2 has been characterized by 9 and 130 GHz electron paramagnetic resonance and high-field electron nuclear double resonance spectroscopies. The hyperfine interactions associated with the ¹⁹F nucleus provide unique signatures of each F_nY· that are readily distinguishable from unlabeled Y·’s. The variability of the abiotic F_nY pK_a’s (6.4 to 7.8) and reduction potentials (−30 to +130 mV relative to Y at pH 7.5) provide probes of enzymatic reactions proposed to involve Y·’s in catalysis and to investigate the importance and identity of hopping Y·’s within redox active proteins proposed to protect them from uncoupled radical chemistry.

Regulation of protein-ligand binding affinity by hydrogen bond pairing

Hydrogen (H)-bonds potentiate diverse cellular functions by facilitating molecular interactions. The mechanism and the extent to which H-bonds regulate molecular interactions are a largely unresolved problem in biology because the H-bonding process continuously competes with bulk water. This interference may significantly alter our understanding of molecular function, for example, in the elucidation of the origin of enzymatic catalytic power. We advance this concept by showing that H-bonds regulate molecular interactions via a hitherto unappreciated donor-acceptor pairing mechanism that minimizes competition with water. On the basis of theoretical and experimental correlations between H-bond pairings and their effects on ligand binding affinity, we demonstrate that H-bonds enhance receptor-ligand interactions when both the donor and acceptor have either significantly stronger or significantly weaker H-bonding capabilities than the hydrogen and oxygen atoms in water. By contrast, mixed strong-weak H-bond pairings decrease ligand binding affinity due to interference with bulk water, offering mechanistic insight into why indiscriminate strengthening of receptor-ligand H-bonds correlates poorly with experimental binding affinity. Further support for the H-bond pairing principle is provided by the discovery and optimization of lead compounds targeting dietary melamine and Clostridium difficile toxins, which are not realized by traditional drug design methods. Synergistic H-bond pairings have therefore evolved in the natural design of high-affinity binding and provide a new conceptual framework to evaluate the H-bonding process in biological systems. Our findings may also guide wider applications of competing H-bond pairings in lead compound design and in determining the origin of enzymatic catalytic power.

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.

Structures of yeast peroxisomal Δ(3),Δ(2)-enoyl-CoA isomerase complexed with acyl-CoA substrate analogues: the importance of hydrogen-bond networks for the reactivity of the catalytic base and the oxyanion hole

Δ(3),Δ(2)-Enoyl-CoA isomerases (ECIs) catalyze the shift of a double bond from 3Z- or 3E-enoyl-CoA to 2E-enoyl-CoA. ECIs are members of the crotonase superfamily. The crotonase framework is used by many enzymes to catalyze a wide range of reactions on acyl-CoA thioesters. The thioester O atom is bound in a conserved oxyanion hole. Here, the mode of binding of acyl-CoA substrate analogues to peroxisomal Saccharomyces cerevisiae ECI (ScECI2) is described. The best defined part of the bound acyl-CoA molecules is the 3',5'-diphosphate-adenosine moiety, which interacts with residues of loop 1 and loop 2, whereas the pantetheine part is the least well defined. The catalytic base, Glu158, is hydrogen-bonded to the Asn101 side chain and is further hydrogen-bonded to the side chain of Arg100 in the apo structure. Arg100 is completely buried in the apo structure and a conformational change of the Arg100 side chain appears to be important for substrate binding and catalysis. The oxyanion hole is formed by the NH groups of Ala70 (loop 2) and Leu126 (helix 3). The O atoms of the corresponding peptide units, Gly69 O and Gly125 O, are both part of extensive hydrogen-bond networks. These hydrogen-bond networks are a conserved feature of the crotonase oxyanion hole and their importance for catalysis is discussed.

Ligand-Induced Proton Transfer and Low-Barrier Hydrogen Bond Revealed by X-ray Crystallography

Ligand binding can change the pKa of protein residues and influence enzyme catalysis. Herein, we report three ultrahigh resolution X-ray crystal structures of CTX-M β-lactamase, directly visualizing protonation state changes along the enzymatic pathway: apo protein at 0.79 Å, precovalent complex with nonelectrophilic ligand at 0.89 Å, and acylation transition state (TS) analogue at 0.84 Å. Binding of the noncovalent ligand induces a proton transfer from the catalytic Ser70 to the negatively charged Glu166, and the formation of a low-barrier hydrogen bond (LBHB) between Ser70 and Lys73, with a length of 2.53 Å and the shared hydrogen equidistant from the heteroatoms. QM/MM reaction path calculations determined the proton transfer barrier to be 1.53 kcal/mol. The LBHB is absent in the other two structures although Glu166 remains neutral in the covalent complex. Our data represents the first X-ray crystallographic example of a hydrogen engaged in an enzymatic LBHB, and demonstrates that desolvation of the active site by ligand binding can provide a protein microenvironment conducive to LBHB formation. It also suggests that LBHBs may contribute to stabilization of the TS in general acid/base catalysis together with other preorganized features of enzyme active sites. These structures reconcile previous experimental results suggesting alternatively Glu166 or Lys73 as the general base for acylation, and underline the importance of considering residue protonation state change when modeling protein-ligand interactions. Additionally, the observation of another LBHB (2.47 Å) between two conserved residues, Asp233 and Asp246, suggests that LBHBs may potentially play a special structural role in proteins.

Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability

Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.

Biosynthetic Incorporation of Fluorinated Amino Acids into Peptides and Proteins

Native and engineered protein biosynthetic machinery processes a wide range of fluorinated α-amino acids for incorporation into peptides and proteins, either as substitutes for structurally similar amino acids normally found in proteins, or as additional ones. In the former case, replacement occurs wherever the normal amino acid is encoded, while the latter method is site-specific. The fluorinated peptides have a diverse variety of interesting properties. The biochemical synthetic methods are straightforward, to the point that they should routinely be assessed as alternatives to traditional solid- and solution-phase peptide synthesis.