A low-barrier hydrogen bond in the catalytic triad of serine proteases? Theory versus experiment

Significance of Histidine Hydrogen-Deuterium Exchange Mass Spectrometry in Protein Structural Biology

Histidine residues play crucial roles in shaping the function and structure of proteins due to their unique ability to act as both acids and bases. In other words, they can serve as proton donors and acceptors at physiological pH. This exceptional property is attributed to the side-chain imidazole ring of histidine residues. Consequently, determining the acid-base dissociation constant (Ka) of histidine imidazole rings in proteins often yields valuable insights into protein functions. Significant efforts have been dedicated to measuring the pKa values of histidine residues in various proteins, with nuclear magnetic resonance (NMR) spectroscopy being the most commonly used technique. However, NMR-based methods encounter challenges in assigning signals to individual imidazole rings and require a substantial amount of proteins. To address these issues associated with NMR-based approaches, a mass-spectrometry-based method known as histidine hydrogen-deuterium exchange mass spectrometry (His-HDX-MS) has been developed. This technique not only determines the pKa values of histidine imidazole groups but also quantifies their solvent accessibility. His-HDX-MS has proven effective across diverse proteins, showcasing its utility. This review aims to clarify the fundamental principles of His-HDX-MS, detail the experimental workflow, explain data analysis procedures and provide guidance for interpreting the obtained results.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

The peptide bond rupture mechanism in the serine proteases: an in silico study based on sequential scale models

Given the importance of serine proteases for biochemical processes, we have studied the peptide bond rupture mechanism using three sequential scale models as representations of the KLK5 enzyme (a protein overexpressed in ovarian cancer). The first model contains the basic functional groups of the residues that conform to the catalytic triad present in serine proteases; the second model contains some additional residues and, finally, the last representation includes all atoms of the KLK5 protein together with 10.000 explicit water molecules. This separation into three scale models allows us to separate the intrinsic reactivity of the catalytic triad from the process taking place in the enzyme. The methodologies employed in this work include full DFT calculations with a dielectric continuum in the first two models and a multi-level setup with a Quantum Mechanics/Molecular Mechanics (QM/MM) partition in the whole protein system. Our results show that the peptide-bond rupture mechanism is a stepwise process involving two proton transfer reactions. The rate-determining step is the second proton transfer from the imidazole group to the amidic nitrogen of the substrate. In addition, we find that the simplest model does not provide accurate results compared to the full protein system. This can be attributed to the electronic stabilization conferred by the residues around the reaction site. Interestingly, the energy profile obtained with the second scale model with additional residues shows the same trends as the full system and could therefore be considered an appropriate model system. It could be used for studying the peptide bond rupture mechanism in case full QM/MM calculations cannot be performed, or as a rapid tool for screening purposes.

Effective prediction of short hydrogen bonds in proteins via machine learning method

Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.

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.

Protonation State of a Histidine Residue in Human Oligopeptide Transporter 1 (hPEPT1) Regulates hPEPT1-Mediated Efflux Activity

To clarify the role of an amino acid residue in the pH-dependent efflux process in Chinese hamster ovary (CHO) cells expressing the human oligopeptide transporter hPEPT1 (CHO/hPEPT1), we determined the effect of extracellular pH on the hPEPT1-mediated efflux process. The efflux of glycylsarcosine (Gly-Sar), a typical substrate for hPEPT1, was determined using an infinite dilution method after cells were preloaded with [³H]-Gly-Sar. The efflux of [³H]-Gly-Sar was stimulated by 5 mM unlabeled hPEPT1 substrates in the medium. This trans-stimulation phenomenon showed that hPEPT1 mediated the efflux of [³H]-Gly-Sar from CHO/hPEPT1 and that hPEPT1 is a bi-directional transporter. We then determined the effect of extracellular pH (varying from 8.0 to 3.5) on the efflux activity. The efflux activity by hPEPT1 decreased with the decrease in extracellular pH. The Henderson-Hasselbälch-type equation, which fitted well to the pH-profile of efflux activity, indicated that a single amino acid residue with a pK_a value of approximately 5.7 regulates the efflux activity. The pH-profile of the efflux activity remained almost unchanged irrespective of the proton gradient across the plasma membrane. In addition, the chemical modification of the histidine residue with diethylpyrocarbonate completely abolished the efflux activity from cells, which could be prevented by the presence of 10 mM Gly-Sar. These data indicate that the efflux process of hPEPT1 is also regulated in a pH-dependent manner by the protonation state of a histidine residue located at or near the substrate recognition site facing the extracellular space.

Proton Bridging in Catalysis by and Inhibition of Serine Proteases of the Blood Cascade System

Inquiries into the participation of short hydrogen bonds in stabilizing transition states and intermediate states in the thrombin, factor Xa, plasmin and activated protein C-catalyzed reactions revealed that specific binding of effectors at Sn, n = 1-4 and S'n, n = 1-3 and at remote exosites elicit complex patterns of hydrogen bonding and involve water networks. The methods employed that yielded these discoveries include; (1) kinetics, especially partial or full kinetic deuterium solvent isotope effects with short cognate substrates and also with the natural substrates, (2) kinetic and structural probes, particularly low-field high-resolution nuclear magnetic resonance (1H NMR), of mechanism-based inhibitors and substrate-mimic peptide inhibitors. Short hydrogen bonds form at the transition states of the catalytic reactions at the active site of the enzymes as they do with mechanism-based covalent inhibitors of thrombin. The emergence of short hydrogen bonds at the binding interface of effectors and thrombin at remote exosites has recently gained recognition. Herein, I describe our contribution, a confirmation of this discovery, by low-field 1H NMR. The principal conclusion of this review is that proton sharing at distances below the sum of van der Waals radii of the hydrogen and both donor and acceptor atoms contribute to the remarkable catalytic prowess of serine proteases of the blood clotting system and other enzymes that employ acid-base catalysis. Proton bridges also play a role in tight binding in proteins and at exosites, i.e., allosteric sites, of enzymes.

Combination of "Buttressing" and "Clothespin" Effects for Reaching the Shortest NHN Hydrogen Bond in Proton Sponge Cations

A series of previously unknown 2,4,5-tri- and 2,4,5,7-tetrasubstituted 1,8-bis(dimethylamino)naphthalenes and their salts with HBF₄ containing bulky spherically shaped substituents (Me, Br, and SiMe₃) in the naphthalene ring has been synthesized. Using XRD analysis of 11 samples, the influence of the so-called "buttressing" and "clothespin" effects on their molecular structure and the NHN hydrogen bond geometry in the solid cations were investigated. The combined action of both effects has been shown to significantly increase the compression of the hydrogen bond. As a result, the previous record of the hydrogen bond shortness (N···N = 2.524 Å) has been surpassed in favor of 2.502 Å found for the tetrafluoroborate of 2,4,5,7-tetramethyl-1,8-bis(dimethylamino)naphthalene. The molecular structure of the latter differs by perfect symmetry and practically barrier-free proton transfer in the [NHN]⁺ bond. On the basis of the results of quantum-chemical calculations, it is suggested that the value of 2.502 Å likely represents or lies very close to the theoretical limit for the short hydrogen bonds between the amine-type nitrogen atoms.

Degradation of a Main Plastic Pollutant Polyethylene Terephthalate by Two Distinct Proteases (Neprilysin and Cutinase-like Enzyme)

In this DFT study, hydrolysis of polyethylene terephthalate (PET), a major cause of plastic pollution, by two distinct enzymes, neprilysin (NEP, a mononuclear metalloprotease) and cutinase-like enzyme (CLE, a serine protease), has been investigated. These enzymes utilize different mechanisms for the degradation of PET. NEP uses either the metal-bound hydroxide attack (MH) mechanism or reverse protonation (RP) mechanism, while CLE utilizes a general acid/base mechanism that includes acylation and deacylation processes. Additionally, the RP mechanism of NEP can proceed through three pathways, RP0, RP1, and RP2. The DFT calculations predict that, among all these mechanisms, the MH mechanism is the energetically most favorable one for the NEP enzyme. In comparison, CLE catalyzes this reaction with a significantly higher barrier. These results suggest that the Lewis acid and nucleophile activations provided by the Zn metal center of NEP are more effective than the hydrogen bonding interactions afforded by the catalytic Ser85-His180-Asp165 triad of CLE. They have provided intrinsic details regarding PET degradation and will pave the way for the design of efficient metal-based catalysts for this critical reaction.

Gut bacteria are essential for normal cuticle development in herbivorous turtle ants

Across the evolutionary history of insects, the shift from nitrogen-rich carnivore/omnivore diets to nitrogen-poor herbivorous diets was made possible through symbiosis with microbes. The herbivorous turtle ants Cephalotes possess a conserved gut microbiome which enriches the nutrient composition by recycling nitrogen-rich metabolic waste to increase the production of amino acids. This enrichment is assumed to benefit the host, but we do not know to what extent. To gain insights into nitrogen assimilation in the ant cuticle we use gut bacterial manipulation, ¹⁵N isotopic enrichment, isotope-ratio mass spectrometry, and ¹⁵N nuclear magnetic resonance spectroscopy to demonstrate that gut bacteria contribute to the formation of proteins, catecholamine cross-linkers, and chitin in the cuticle. This study identifies the cuticular components which are nitrogen-enriched by gut bacteria, highlighting the role of symbionts in insect evolution, and provides a framework for understanding the nitrogen flow from nutrients through bacteria into the insect cuticle.

The hydrogen-bond configuration modulates the energy transfer efficiency in helical protein nanotubes

Energy transport in proteins is critical to a variety of physical, chemical, and biological processes in living organisms. While strenuous efforts have been made to study vibrational energy transport in proteins, thermal transport processes across the most fundamental building blocks of proteins, i.e. helices, are not well understood. This work studies energy transport in a group of "isomer" helices. The π-helix is shown to have the highest thermal conductivity, 110% higher than that of the α-helix and 207% higher than that of the 310-helix. The H-bond connectivity is found to govern thermal transport mechanisms including the phonon spectral energy density, dispersion, mode-specific transport, group velocity, and relaxation time. The energy transport is strongly correlated with the H-bond strength which is also modulated by the H-bond connectivity. These fundamental insights provide a novel perspective for understanding energy transfer in proteins and guiding a rational molecule-level design of novel materials with configurable H-bonds.

Single-Atom Catalysts across the Periodic Table

Isolated atoms featuring unique reactivity are at the heart of enzymatic and homogeneous catalysts. In contrast, although the concept has long existed, single-atom heterogeneous catalysts (SACs) have only recently gained prominence. Host materials have similar functions to ligands in homogeneous catalysts, determining the stability, local environment, and electronic properties of isolated atoms and thus providing a platform for tailoring heterogeneous catalysts for targeted applications. Within just a decade, we have witnessed many examples of SACs both disrupting diverse fields of heterogeneous catalysis with their distinctive reactivity and substantially enriching our understanding of molecular processes on surfaces. To date, the term SAC mostly refers to late transition metal-based systems, but numerous examples exist in which isolated atoms of other elements play key catalytic roles. This review provides a compositional encyclopedia of SACs, celebrating the 10th anniversary of the introduction of this term. By defining single-atom catalysis in the broadest sense, we explore the full elemental diversity, joining different areas across the whole periodic table, and discussing historical milestones and recent developments. In particular, we examine the coordination structures and associated properties accessed through distinct single-atom-host combinations and relate them to their main applications in thermo-, electro-, and photocatalysis, revealing trends in element-specific evolution, host design, and uses. Finally, we highlight frontiers in the field, including multimetallic SACs, atom proximity control, and possible applications for multistep and cascade reactions, identifying challenges, and propose directions for future development in this flourishing field.

Noninvasive Epidermal Metabolite Profiling

A buffer placed in brief contact in the skin was assayed by ¹H NMR spectroscopy. We found that this passive extraction of the skin surface yields abundant metabolites. Metabolites of the skin surface originate from a variety of sources, including the sweat gland, which produces lactate from the glucose received from its capillary bed. Little is known about how metabolites resident on and within the skin surface respond to a metabolic or hemodynamic perturbation. As a possible application of epidermal metabolite profiling, we asked whether metabolites extracted from the skin surface are indicative of heart failure. The levels of lactate and other molecules were significantly lower in patients in heart failure than in individuals who reported healthy heart function, possibly due to reduced blood flow to the sweat gland resulting in a lack of tissue perfusion. Most amino acids were unchanged in levels, except for glycine and serine that increased as a percentage of all amino acids. These results have the potential in the long term to help decide the extent to which a patient has heart failure for which objective measures are lacking. Moreover, the results suggest that epidermal metabolite profiling may be useful for other assessments of human health.

Improving production of Streptomyces griseus trypsin for enzymatic processing of insulin precursor

Background

Trypsin has many applications in food and pharmaceutical manufacturing. Although commercial trypsin is usually extracted from porcine pancreas, this source carries the risks of infectivity and immunogenicity. Microbial Streptomyces griseus trypsin (SGT) is a prime alternative because it possesses efficient hydrolysis activity without such risks. However, the remarkable hydrolysis efficiency of SGT causes autolysis, and five autolysis sites, R21, R32, K122, R153, and R201, were identified from its autolysate.

Results

The tbcf (K101A, R201V) mutant was screened by a directed selection approach for improved activity in flask culture (60.85 ± 3.42 U mL⁻¹, increased 1.5-fold). From the molecular dynamics simulation, in the K101A/R201V mutant the distance between the catalytical residues D102 and H57 was shortened to 6.5 Å vs 7.0 Å in the wild type, which afforded the improved specific activity of 1527.96 ± 62.81 U mg⁻¹. Furthermore, the production of trypsin was increased by 302.8% (689.47 ± 6.78 U mL⁻¹) in a 3-L bioreactor, with co-overexpression of chaperones SSO2 and UBC1 in Pichia pastoris.

Conclusions

SGT protein could be a good source of trypsin for insulin production. As a result of the hydrolysates analysis and direct selection, the activity of the tbcf (K101A, R201V) mutant increased 1.5-fold. Furthermore, the production of trypsin was improved threefold by overexpressing chaperone protein in Pichia pastoris. Future studies should investigate the application of SGT to insulin and pharmaceutical manufacturing.

Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects

Electronic and nuclear quantum effects determine the symmetry and highly downfield ¹H NMR chemical shifts of short hydrogen bonds.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

The rhomboid protease GlpG has weak interaction energies in its active site hydrogen bond network

Rhomboid proteases are membrane-integrated enzymes that hydrolyze peptide bonds in the transmembrane domains of protein substrates. Gaffney and Hong experimentally determine interaction energies between active site residues to reveal weak coupling, which may explain the slow proteolysis mediated by GlpG.

Intramembrane rhomboid proteases are of particular interest because of their function to hydrolyze a peptide bond of a substrate buried in the membrane. Crystal structures of the bacterial rhomboid protease GlpG have revealed a catalytic dyad (Ser201-His254) and oxyanion hole (His150/Asn154/the backbone amide of Ser201) surrounded by the protein matrix and contacting a narrow water channel. Although multiple crystal structures have been solved, the catalytic mechanism of GlpG is not completely understood. Because it is a serine protease, hydrogen bonding interactions between the active site residues are thought to play a critical role in the catalytic cycle. Here, we dissect the interaction energies among the active site residues His254, Ser201, and Asn154 of Escherichia coli GlpG, which form a hydrogen bonding network. We combine double mutant cycle analysis with stability measurements using steric trapping. In mild detergent, the active site residues are weakly coupled with interaction energies (ΔΔG_Inter) of ‒1.4 kcal/mol between His254 and Ser201 and ‒0.2 kcal/mol between Ser201 and Asn154. Further, by analyzing the propagation of single mutations of the active site residues, we find that these residues are important not only for function but also for the folding cooperativity of GlpG. The weak interaction between Ser and His in the catalytic dyad may partly explain the unusually slow proteolysis by GlpG compared with other canonical serine proteases. Our result suggests that the weak hydrogen bonds in the active site are sufficient to carry out the proteolytic function of rhomboid proteases.

Direct evidence of a low barrier hydrogen bond in the catalytic triad of a Serine protease

Serine proteases are one of the largest groups of enzymes, found in both eukaryotes and prokaryotes, and are responsible for many different functions. The detailed information about the hydrogen-bonds in the catalytic triad (Asp…His…Ser) of these enzymes is of importance in order to fully understand the mechanism of action. The aspartate of the triad is hydrogen bonded to the histidine but the exact nature of this bond has been under discussion for some time. It is either a common short ionic hydrogen bond (SIHB) or a delocalized low barrier hydrogen bond (LBHB) were the hydrogen bond is shorter. So far, the evidence for LBHB in proteins have not been conclusive. Here we show clear NMR evidence that LBHB does exist in NS3, a serine protease from Dengue. The one bond coupling constant between the hydrogen and nitrogen was shown to be only 52 Hz instead of the usual 90 Hz. This together with a 1H chemical shift of 19.93 ppm is evidence that the hydrogen bond distance between His and Asp is shorter than for SIHB. Our result clearly shows the existence of LBHB and will help in understanding the mechanism of the catalytic triad in the important group of serine proteases.

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.

Solvent effect on hydrogen bonded Tyr⋯Asp⋯Arg triads: Enzymatic catalyzed model system

The hydrogen bond plays a vital role in structural arrangement, intermediate state stabilization, materials function, and biological activity of certain enzymatic reactions. The solvent and electronic effects on hydrogen bonds are illustrated employing the polarizable contimuum model at B3LYP/6-311++G(d,p) level. Geometry optimizations reflect the significant solvent and electronic effect. The proton departs spontaneously upon oxidation from the hydroxyl group of tyrosyl in hydrogen bonded Tyr⋯Asp⋯Arg triads in both gas phase and solvents. The electron transfer isomers are observed for anionic triads, no matter what the solvent is. The difference of distance between two hydrogen bonds is enlarged in solvent as compared to that in gas phase. The electronic effect on IR spectra is distinctive. The tyrosyl fragment tends to be oxidized and the arginine moiety is easier to bind an excess electron. The variations of chemical shift and spin-spin coupling constant are more significant upon electron transfer than upon solvent dielectric constant. The augmentation of solvent dielectric constant stabilizes the system, enhances the difference of isomers, and increases the vertical ionization potential and vertical electron affinity values.

Hydration effect on proton transfer in melamine-cyanuric acid complex

Self-assembly of melamine-cyanuric acid (MC) leads to urinary tract calculi and renal failure. The hydration effects on molecular geometry, the IR spectra, the frontier molecular orbital, the energy barrier of proton transfer (PT), as well as the stability of MC were explored by density functional theory (DFT) calculations. The intramolecular PT breaks the big π-conjugated ring of melamine or converts the p-π conjugation (:N-C'=O) to π-π conjugation (O=C-N=C') of cyanuric acid. The intermolecular PT varies the coupling between melamine and cyanuric acid from pure hydrogen bonds (Na…HNd and NH…O) to the cooperation of cation…anion electrostatic interaction (NaH(+)…Nd (-)) and two NH…O hydrogen bonds. Distinct IR spectra shifts occur for Na…HNd stretching mode upon PT, i.e., blue-shift upon intramolecular PT and red-shift upon intermolecular PT. It is expected that the PT would inhibit the generation of rosette-like structure or one-dimensional tape conformer for the MC complexes. Hydration obviously effects the local geometric structure around the water binding site, as well as the IR spectra of NH…O and N…HN hydrogen bonds. Hydration decreases the intramolecular PT barrier from ~45 kcal mol(-1) in anhydrous complex to ~11.5 kcal mol(-1) in trihydrated clusters. While, the hydration effects on intermolecular PT barrier is slight. The relative stability of MC varies slightly by hydration due to the strong hydrogen bond interaction between melamine and cyanuric acid fragments. Graphical Abstract Hydration effect on proton transfer in melamine-cyanuric acid complex.

Assessing the Indirect Photochemical Transformation of Dissolved Combined Amino Acids through the Use of Systematically Designed Histidine-Containing Oligopeptides

Photooxidation is an important abiotic transformation pathway for amino acids (AAs) in sunlit waters. Although dissolved free AAs are well studied, the photooxidation of dissolved combined AAs (DCAAs) remains poorly investigated. This study is a systematic investigation of the effect of neighboring photostable AA residues (i.e., aliphatic, cationic, anionic, or aromatic residues) on the environmental indirect photochemical transformation of histidine (His) in His-containing oligopeptides. The pKa values of His residues in the studied oligopeptides were found to be between 4.3 and 8.1. Accordingly, the phototransformation rate constants of the His-containing oligopeptides were highly pH-dependent in an environmentally relevant pH range with higher reactivity for neutral His than for the protonated species. The photostable AA residues significantly modulated the photoreactivity of oligopeptides either through altering the accessibility of His to photochemically produced oxidants or through shifting the pKa values of His residues. In addition, the influence of neighboring photostable AA residues on the sorption-enhanced phototransformation of oligopeptides in solutions containing chromophoric dissolved organic matter (CDOM) was assessed. The constituent photostable AA residues promoted sorption of His-containing oligopeptides to CDOM macromolecules through electrostatic attraction, hydrophobic effects, and/or low-barrier hydrogen bonds, and subsequently increased the apparent phototransformation rate constants by up to 2 orders of magnitude.

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.

Contribution of a low-barrier hydrogen bond to catalysis is not significant in ketosteroid isomerase

Low-barrier hydrogen bonds (LBHBs) have been proposed to have important influences on the enormous reaction rate increases achieved by many enzymes. Δ(5)-3-ketosteroid isomerase (KSI) catalyzes the allylic isomerization of Δ(5)-3-ketosteroid to its conjugated Δ(4)-isomers at a rate that approaches the diffusion limit. Tyr14, a catalytic residue of KSI, has been hypothesized to form an LBHB with the oxyanion of a dienolate steroid intermediate generated during the catalysis. The unusual chemical shift of a proton at 16.8 ppm in the nuclear magnetic resonance spectrum has been attributed to an LBHB between Tyr14 Oη and C3-O of equilenin, an intermediate analogue, in the active site of D38N KSI. This shift in the spectrum was not observed in Y30F/Y55F/D38N and Y30F/Y55F/Y115F/D38N mutant KSIs when each mutant was complexed with equilenin, suggesting that Tyr14 could not form LBHB with the intermediate analogue in these mutant KSIs. The crystal structure of Y30F/Y55F/Y115F/D38N-equilenin complex revealed that the distance between Tyr14 Oη and C3-O of the bound steroid was within a direct hydrogen bond. The conversion of LBHB to an ordinary hydrogen bond in the mutant KSI reduced the binding affinity for the steroid inhibitors by a factor of 8.1-11. In addition, the absence of LBHB reduced the catalytic activity by only a factor of 1.7-2. These results suggest that the amount of stabilization energy of the reaction intermediate provided by LBHB is small compared with that provided by an ordinary hydrogen bond in KSI.

Atomic resolution structure of a lysine-specific endoproteinase fromLysobacter enzymogenessuggests a hydroxyl group bound to the oxyanion hole

Lysobacter enzymogeneslysyl endoproteinase (LysC) is a trypsin-type serine protease with a high pH optimum that hydrolyses all Lys-Xaa peptide bonds. The high specificity of LysC renders it useful for biotechnological purposes. The K30R variant of a related lysyl endoproteinase fromAchromobacter lyticushas favourable enzymatic properties that might be transferrable to LysC. To visualize structural differences in the substrate-binding sites, the crystal structures of wild-type and the K30R variant of LysC were determined. The mutation is located at a distance of 12 Å from the catalytic triad and subtly changes the surface properties of the substrate-binding site. The high pH optimum of LysC can be attributed to electrostatic effects of an aromatic Tyr/His stack on the catalytic aspartate and is a general feature of this enzyme subfamily. LysC crystals in complex with the covalent inhibitorNα-p-tosyl-lysyl chloromethylketone yielded data to 1.1 and 0.9 Å resolution, resulting in unprecedented precision of the active and substrate-binding sites for this enzyme subfamily. Error estimates on bond lengths and difference electron density indicate that instead of the expected oxyanion a hydroxyl group binds to the partially solvent-exposed oxyanion hole. Protonation of the alkoxide catalytic intermediate might be a recurring feature during serine protease catalysis.

Hydrogen bonding involving side chain exchangeable groups stabilizes amyloid quarternary structure

The amyloid β-peptide (Aβ) is the major structural component of amyloid fibrils in the plaques of brains of Alzheimer's disease patients. Numerous studies have addressed important aspects of secondary and tertiary structure of fibrils. In electron microscopic images, fibrils often bundle together. The mechanisms which drive the association of protofilaments into bundles of fibrils are not known. We show here that amino acid side chain exchangeable groups like e.g. histidines can provide useful restraints to determine the quarternary assembly of an amyloid fibril. Exchangeable protons are only observable if a side chain hydrogen bond is formed and the respective protons are protected from exchange. The method relies on deuteration of the Aβ peptide. Exchangeable deuterons are substituted with protons, before fibril formation is initiated.

A QM/MM MD study of the pH-dependent ring-opening catalysis and lid motif flexibility in glucosamine 6-phosphate deaminase

QM/MM MD and MM MD simulations reveal pH-dependent proton-shuttle ring-opening mechanisms of GlcN6P and dynamical behavior of the lid motif inSmuNagB.

Proton bridging in the interactions of thrombin with hirudin and its mimics

Thrombin is the pivotal serine protease enzyme in the blood cascade system and thus a target of drug design for control of its activity. The most efficient nonphysiologic inhibitor of thrombin is hirudin, a naturally occurring small protein. Hirudin and its synthetic mimics employ a range of hydrogen bonding, salt bridging, and hydrophobic interactions with thrombin to achieve tight binding with K(i) values in the nano- to femtomolar range. The one-dimensional (1)H nuclear magnetic resonance spectrum recorded at 600 MHz reveals a resonance 15.33 ppm downfield from silanes in complexes between human α-thrombin and r-hirudin in pH 5.6-8.8 buffers and between 5 and 35 °C. There is also a resonance between 15.17 and 15.54 ppm seen in complexes of human α-thrombin with hirunorm IV, hirunorm V, an Nα(Me)Arg peptide, RGD-hirudin, and Nα-2-naphthylsulfonyl-glycyl-DL-4-amidinophenylalanyl-piperidide acetate salt (NAPAP), while there is no such low-field resonance observed in a complex of porcine trypsin and NAPAP. The chemical shifts suggest that these resonances represent H-bonded environments. H-Donor-acceptor distances in the corresponding H-bonds are estimated to be <2.7 Å. Addition of Phe-Pro-Arg-chloromethylketone (PPACK) to a complex of human α-thrombin with r-hirudin results in an additional signal at 18.03 ppm, which is 0.10 ppm upfield from the observed signal [Kovach, I. M., et al. (2009) Biochemistry 48, 7296-7304] for thrombin covalently modified with PPACK. In contrast, the peak at 15.33 ppm remains unchanged. The fractionation factors for the thrombin-hirudin complexes are near 1.0 within 20% error. The most likely site of the short H-bond in complexes of thrombin with the hirudin family of inhibitors is in the hydrophobic patch of the C-terminus of hirudin where Glu(57') and Glu(58') are embedded and interact with Arg(75) and Arg(77) and their solvate water (on thrombin). Glu(57') and Glu(58') present in the hirudin family of inhibitors make up a key binding epitope of fibrinogen, thrombin's prime substrate, which lends substantial interest to the short hydrogen bond as a binding element at the fibrinogen recognition site.

The catalytic aspartate is protonated in the Michaelis complex formed between trypsin and an in vitro evolved substrate-like inhibitor: a refined mechanism of serine protease action

The mechanism of serine proteases prominently illustrates how charged amino acid residues and proton transfer events facilitate enzyme catalysis. Here we present an ultrahigh resolution (0.93 Å) x-ray structure of a complex formed between trypsin and a canonical inhibitor acting through a substrate-like mechanism. The electron density indicates the protonation state of all catalytic residues where the catalytic histidine is, as expected, in its neutral state prior to the acylation step by the catalytic serine. The carboxyl group of the catalytic aspartate displays an asymmetric electron density so that the O_δ2–C_γ bond appears to be a double bond, with O_δ2 involved in a hydrogen bond to His-57 and Ser-214. Only when Asp-102 is protonated on O_δ1 atom could a density functional theory simulation reproduce the observed electron density. The presence of a putative hydrogen atom is also confirmed by a residual mF_obs − DF_calc density above 2.5 σ next to O_δ1. As a possible functional role for the neutral aspartate in the active site, we propose that in the substrate-bound form, the neutral aspartate residue helps to keep the pK_a of the histidine sufficiently low, in the active neutral form. When the histidine receives a proton during the catalytic cycle, the aspartate becomes simultaneously negatively charged, providing additional stabilization for the protonated histidine and indirectly to the tetrahedral intermediate. This novel proposal unifies the seemingly conflicting experimental observations, which were previously seen as either supporting the charge relay mechanism or the neutral pK_a histidine theory.

Structural and kinetic analysis of Bacillus subtilis N-acetylglucosaminidase reveals a unique Asp-His dyad mechanism

Three-dimensional structures of NagZ of Bacillus subtilis, the first structures of a two-domain β-N-acetylglucosaminidase of family 3 of glycosidases, were determined with and without the transition state mimicking inhibitor PUGNAc bound to the active site, at 1.84- and 1.40-Å resolution, respectively. The structures together with kinetic analyses of mutants revealed an Asp-His dyad involved in catalysis: His(234) of BsNagZ acts as general acid/base catalyst and is hydrogen bonded by Asp(232) for proper function. Replacement of both His(234) and Asp(232) with glycine reduced the rate of hydrolysis of the fluorogenic substrate 4'-methylumbelliferyl N-acetyl-β-D-glucosaminide 1900- and 4500-fold, respectively, and rendered activity pH-independent in the alkaline range consistent with a role of these residues in acid/base catalysis. N-Acetylglucosaminyl enzyme intermediate accumulated in the H234G mutant and β-azide product was formed in the presence of sodium azide in both mutants. The Asp-His dyad is conserved within β-N-acetylglucosaminidases but otherwise absent in β-glycosidases of family 3, which instead carry a "classical" glutamate acid/base catalyst. The acid/base glutamate of Hordeum vulgare exoglucanase (Exo1) superimposes with His(234) of the dyad of BsNagZ and, in contrast to the latter, protrudes from a second domain of the enzyme into the active site. This is the first report of an Asp-His catalytic dyad involved in hydrolysis of glycosides resembling in function the Asp-His-Ser triad of serine proteases. Our findings will facilitate the development of mechanism-based inhibitors that selectively target family 3 β-N-acetylglucosaminidases, which are involved in bacterial cell wall turnover, spore germination, and induction of β-lactamase.

Characterization of Bafinivirus main protease autoprocessing activities

The production of functional nidovirus replication-transcription complexes involves extensive proteolytic processing by virus-encoded proteases. In this study, we characterized the viral main protease (M(pro)) of the type species, White bream virus (WBV), of the newly established genus Bafinivirus (order Nidovirales, family Coronaviridae, subfamily Torovirinae). Comparative sequence analysis and mutagenesis data confirmed that the WBV M(pro) is a picornavirus 3C-like serine protease that uses a Ser-His-Asp catalytic triad embedded in a predicted two-β-barrel fold, which is extended by a third domain at its C terminus. Bacterially expressed WBV M(pro) autocatalytically released itself from flanking sequences and was able to mediate proteolytic processing in trans. Using N-terminal sequencing of autoproteolytic processing products we tentatively identified Gln↓(Ala, Thr) as a substrate consensus sequence. Mutagenesis data provided evidence to suggest that two conserved His and Thr residues are part of the S1 subsite of the enzyme's substrate-binding pocket. Interestingly, we observed two N-proximal and two C-proximal autoprocessing sites in the bacterial expression system. The detection of two major forms of M(pro), resulting from processing at two different N-proximal and one C-proximal site, in WBV-infected epithelioma papulosum cyprini cells confirmed the biological relevance of the biochemical data obtained in heterologous expression systems. To our knowledge, the use of alternative M(pro) autoprocessing sites has not been described previously for other nidovirus M(pro) domains. The data presented in this study lend further support to our previous conclusion that bafiniviruses represent a distinct group of viruses that significantly diverged from other phylogenetic clusters of the order Nidovirales.

Catalytic reaction mechanism of acetylcholinesterase determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations

Acetylcholinesterase (AChE) is a remarkably efficient serine hydrolase responsible for the termination of impulse signaling at cholinergic synapses. By employing Born-Oppenheimer molecular dynamics simulations with a B3LYP/6-31G(d) QM/MM potential and the umbrella sampling method, we have characterized its complete catalytic reaction mechanism for hydrolyzing neurotransmitter acetylcholine (ACh) and determined its multistep free-energy reaction profiles for the first time. In both acylation and deacylation reaction stages, the first step involves the nucleophilic attack on the carbonyl carbon, with the triad His447 serving as the general base, and leads to a tetrahedral covalent intermediate stabilized by the oxyanion hole. From the intermediate to the product, the orientation of the His447 ring needs to be adjusted very slightly, and then, the proton transfers from His447 to the product, and the break of the scissile bond happens spontaneously. For the three-pronged oxyanion hole, it only makes two hydrogen bonds with the carbonyl oxygen at either the initial reactant or the final product state, but the third hydrogen bond is formed and stable at all transition and intermediate states during the catalytic process. At the intermediate state of the acylation reaction, a short and low-barrier hydrogen bond (LBHB) is found to be formed between two catalytic triad residues His447 and Glu334, and the spontaneous proton transfer between two residues has been observed. However, it is only about 1-2 kcal/mol stronger than the normal hydrogen bond. In comparison with previous theoretical investigations of the AChE catalytic mechanism, our current study clearly demonstrates the power and advantages of employing Born-Oppenheimer ab initio QM/MM MD simulations in characterizing enzyme reaction mechanisms.

Combined high-resolution neutron and X-ray analysis of inhibited elastase confirms the active-site oxyanion hole but rules against a low-barrier hydrogen bond

To help resolve long-standing questions regarding the catalytic activity of the serine proteases, the structure of porcine pancreatic elastase has been analyzed by high-resolution neutron and X-ray crystallography. To mimic the tetrahedral transition intermediate, a peptidic inhibitor was used. A single large crystal was used to collect room-temperature neutron data to 1.65 A resolution and X-ray data to 1.20 A resolution. Another crystal provided a low-temperature X-ray data set to 0.94 A resolution. The neutron data are to higher resolution than previously reported for a serine protease and the X-ray data are comparable with other studies. The neutron and X-ray data show that the hydrogen bond between His57 and Asp102 (chymotrypsin numbering) is 2.60 A in length and that the hydrogen-bonding hydrogen is 0.80-0.96 A from the histidine nitrogen. This is not consistent with a low-barrier hydrogen which is predicted to have the hydrogen midway between the donor and acceptor atom. The observed interaction between His57 and Asp102 is essentially a short but conventional hydrogen bond, sometimes described as a short ionic hydrogen bond. The neutron analysis also shows that the oxygen of the oxopropyl group of the inhibitor is present as an oxygen anion rather than a hydroxyl group, supporting the role of the "oxyanion hole" in stabilizing the tetrahedral intermediate in catalysis.

Proton bridging in the interactions of thrombin with small inhibitors

Thrombin is the pivotal serine protease enzyme in the blood cascade system. Phe-Pro-Arg-chloromethylketone (PPACK), phosphate, and phosphonate ester inhibitors form a covalent bond with the active-site Ser of thrombin. PPACK, a mechanism-based inhibitor, and the phosphate/phosphonate esters form adducts that mimic intermediates formed in reactions catalyzed by thrombin. Therefore, the dependence of the inhibition of human alpha-thrombin on the concentration of these inhibitors, pH, and temperature was investigated. The second-order rate constant (ki/Ki) and the inhibition constant (Ki) for inhibition of human alpha-thrombin by PPACK are (1.1 +/- 0.2) x 10(7) M(-1) s(-1) and (2.4 +/- 1.3) x 10(-8) M, respectively, at pH 7.00 in 0.05 M phosphate buffer and 0.15 M NaCl at 25.0 +/- 0.1 degrees C, in good agreement with previous reports. The activation parameters at pH 7.00 in 0.05 M phosphate buffer and 0.15 M NaCl are as follows: DeltaH = 10.6 +/- 0.7 kcal/mol, and DeltaS = 9 +/- 2 cal mol(-1) degrees C(-1). The pH dependence of the second-order rate constants of inhibition is bell-shaped. Values of pKa1 and pKa2 are 7.3 +/- 0.2 and 8.8 +/- 0.3, respectively, at 25.0 +/- 0.1 degrees C. A phosphate and a phosphonate ester inhibitor gave higher values, 7.8 and 8.0 for pKa1 and 9.3 and 8.6 for pKa2, respectively. They inhibit thrombin more than 6 orders of magnitude less efficiently than PPACK does. The deuterium solvent isotope effect for the second-order rate constant at pH 7.0 and 8.3 at 25.0 +/- 0.1 degrees C is unity within experimental error in all three cases, indicating the absence of proton transfer in the rate-determining step for the association of thrombin with the inhibitors, but in a 600 MHz 1H NMR spectrum of the inhibition adduct at pH 6.7 and 30 degrees C, a peak at 18.10 ppm with respect to TSP appears with PPACK, which is absent in the 1H NMR spectrum of a solution of the enzyme between pH 5.3 and 8.5. The peak at low field is an indication of the presence of a short-strong hydrogen bond (SSHB) at the active site in the adduct. The deuterium isotope effect on this hydrogen bridge is 2.2 +/- 0.2 (phi = 0.45). The presence of an SSHB is also established with a signal at 17.34 ppm for a dealkylated phosphate adduct of thrombin.

Low-barrier hydrogen bond in photoactive yellow protein

Low-barrier hydrogen bonds (LBHBs) have been proposed to play roles in protein functions, including enzymatic catalysis and proton transfer. Transient formation of LBHBs is expected to stabilize specific reaction intermediates. However, based on experimental results and theoretical considerations, arguments against the importance of LBHB in proteins have been raised. The discrepancy is caused by the absence of direct identification of the hydrogen atom position. Here, we show by high-resolution neutron crystallography of photoactive yellow protein (PYP) that a LBHB exists in a protein, even in the ground state. We identified approximately 87% (819/942) of the hydrogen positions in PYP and demonstrated that the hydrogen bond between the chromophore and E46 is a LBHB. This LBHB stabilizes an isolated electric charge buried in the hydrophobic environment of the protein interior. We propose that in the excited state the fast relaxation of the LBHB into a normal hydrogen bond is the trigger for photo-signal propagation to the protein moiety. These results give insights into the novel roles of LBHBs and the mechanism of the formation of LBHBs.