Unusual 1H NMR chemical shifts support (His) C(epsilon) 1...O==C H-bond: proposal for reaction-driven ring flip mechanism in serine protease catalysis

Exploring the pH dependence of an improved PETase

Enzymatic recycling of plastic and especially of polyethylene terephthalate (PET) has shown great potential to reduce its negative impact on our society. PET hydrolases (PETases) have been optimized using rational design and machine learning, but the mechanistic details of the PET depolymerization process remain unclear. Belonging to the carboxylic-ester hydrolase family with a canonical Ser-His-Asp catalytic triad, their observed alkaline pH optimum is generally thought to be related to the protonation state of the catalytic His. Here, we explore this aspect in the context of LCCICCG, an optimized PETase, derived from the leaf-branch compost cutinase enzyme. We use NMR to identify the dominant tautomeric structure of the six histidines. Five show surprisingly low pKa values below 4.0, whereas the catalytic H242 in the active enzyme displays a pKa value that varies from 4.9 to 4.7 when temperatures increase from 30°C to 50°C. Whereas the hydrolytic activity of the enzyme toward a soluble substrate can be modeled by the corresponding protonation/deprotonation curve, an important discrepancy is found when the substrate is the solid plastic. This opens the way to further mechanistic understanding of the PETase activity and underscores the importance of studying the enzyme at the liquid-solid interface.

C-H Groups as Donors in Hydrogen Bonds: A Historical Overview and Occurrence in Proteins and Nucleic Acids

Hydrogen bonds constitute a unique type of non-covalent interaction, with a critical role in biology. Until fairly recently, the canonical view held that these bonds occur between electronegative atoms, typically O and N, and that they are mostly electrostatic in nature. However, it is now understood that polarized C-H groups may also act as hydrogen bond donors in many systems, including biological macromolecules. First recognized from physical chemistry studies, C-H…X bonds were visualized with X-ray crystallography sixty years ago, although their true significance has only been recognized in the last few decades. This review traces the origins of the field and describes the occurrence and significance of the most important C-H…O bonds in proteins and nucleic acids.

Competition between Binding to Various Sites of Substituted Imidazoliums

The imidazolium cation has a number of different sites that can interact with a nucleophile. Adding a halogen atom (X) or a chalcogen (YH) group introduces the possibility of an NX···nuc halogen or NY···nuc chalcogen bond, which competes against the various H-bonds (NH and CH donors) as well as the lone pair···π interaction wherein the nucleophile lies above the plane of the cation. Substituted imidazoliums are paired with the NH₃ base, and the various different complexes are evaluated by density functional theory (DFT) calculations. The strength of XB and YB increases quickly along with the size and polarizability of the X/Y atom, and this sort of bond is the strongest for the heavier Br, I, Se, and Te atoms, followed by the NH···N H-bond, but this order reverses for Cl and S. The various CH···N H-bonds are comparable to one another and to the lone pair···π bond, all with interaction energies of 10-13 kcal/mol, values which show very little dependence upon the substituent placed on the imidazolium.

Crystal structure of the feruloyl esterase from Lentilactobacillus buchneri reveals a novel homodimeric state

Ferulic acid is a common constituent of the plant cell-wall matrix where it decorates and can crosslink mainly arabinoxylans to provide structural reinforcement. Microbial feruloyl esterases (FAEs) specialize in catalyzing hydrolysis of the ester bonds between phenolic acids and sugar residues in plant cell-wall polysaccharides such as arabinoxylan to release cinnamoyl compounds. Feruloyl esterases from lactic acid bacteria (LAB) have been highlighted as interesting enzymes for their potential applications in the food and pharmaceutical industries; however, there are few studies on the activity and structure of FAEs of LAB origin. Here, we report the crystal structure and biochemical characterization of a feruloyl esterase (LbFAE) from Lentilactobacillus buchneri, a LAB strain that has been used as a silage additive. The LbFAE structure was determined in the absence and presence of product (FA) and reveals a new type of homodimer association not previously observed for fungal or bacterial FAEs. The two subunits associate to restrict access to the active site such that only single FA chains attached to arabinoxylan can be accommodated, an arrangement that excludes access to FA cross-links between arabinoxylan chains. This narrow specificity is further corroborated by the observation that no FA dimers are produced, only FA, when feruloylated arabinoxylan is used as substrate. Docking of arabinofuranosyl-ferulate in the LbFAE structure highlights the restricted active site and lends further support to our hypothesis that LbFAE is specific for single FA side chains in arabinoxylan.

To what extent do bond length and angle govern the 13C and 1H NMR response to weak CH⋯O hydrogen bonds? A case study of caffeine and theophylline cocrystals

Weak hydrogen bonds are important structure-directing elements in supramolecular chemistry and biochemistry. We consider here weak CH⋯O hydrogen bonds in a series of cocrystals of theophylline and caffeine and assess to what extent the CH⋯O distance and angle govern the observed ¹³C and ¹H isotropic chemical shifts. Gauge-including projector-augmented wave density functional theory (GIPAW DFT) calculations consistently predict a decrease in the ¹³C and ¹H magnetic shielding constants upon hydrogen bond formation on the order of 2-5 ppm (¹³C) and 1-2 ppm (¹H). These trends are reproduced using the machine-learning approach implemented in ShiftML. Experimental ¹³C and ¹H chemical shifts obtained for powdered samples using one-dimensional NMR spectroscopy as well as heteronuclear correlation (HETCOR) spectroscopy correlate well with the GIPAW DFT results. However, the experimental ¹³C NMR response only correlates moderately well with the hydrogen bond length and angle, while the experimental ¹H chemical shifts only show very weak correlations to these local structural elements. DFT computations on isolated imidazole-formaldehyde models show that the ¹³C and ¹H chemical shifts generally decrease with the C⋯O distance but show no clear dependence on the CH⋯O angle. These results demonstrate that the ¹³C and ¹H response to weak CH⋯O hydrogen bonding is influenced significantly by additional weak contacts within cocrystal heterodimeric units.

Weakly Bound Dimer of a Diaryloxygermylene Derived from a tBuPh2Si-Substituted 2,2′-Methylenediphenol

Novel diaryloxygermylenes have been prepared by the reaction of Lappert’s germylene, Ge[N(SiMe3)2]2, with 2,2′-methylenediphenols bearing different substituents. The bulkiness of the substituents on the ortho positions of the phenolic oxygen (6 and 6′ positions) affects the structure of the products both in the solid-state and in solution. When the ortho substituents are SitBuPh2, the diaryloxygemylene crystalizes as a weakly bound dimer with intermolecular Ge…O distances of ca. 3.0 Å and exists as a monomer in solution. In contrast, the germylene with SiMePh2 groups as the ortho substituents form a tightly bound dimer featuring a Ge2O2 rhombus with cis-oriented terminal aryloxy groups in the crystalline state, which is confirmed to be maintained in solution through the VT (variable-temperature)-1H NMR studies. To the best of our knowledge, the former dimeric structure is unprecedented in the family of dioxytetrylenes.

Characterization and Photofragmentation Studies of the Benzimidazole Homodimer: Evidence for Excited-State Charge-Coupled Proton Transfer

The spectroscopic characterization of the benzimidazole (BIM) homodimer was carried out in a molecular beam in the ground state as well as in the cationic state using the R2PI and RIDIR methods. Primarily, interest in the dimer was due to the observation of a proton-transferred BIM fragment at energies well below its thermodynamic threshold (i.e., barely above the ionization energy of the dimer where fragmentation was not expected). The detailed photofragmentation studies of the homodimer combined with spectroscopic observations and quantum chemical computations of the excited states established that the proton transfer from one subunit to the other occurs via conical intersections connecting the locally excited state, the charge-transfer state, and the ground state. In this study, we have also determined the N-H···N hydrogen bond dissociation energy in the ground state and in the cationic state to be 10.36 ± 0.14 and 27.55 ± 0.20 kcal mol^-1, respectively. Incidentally, this happens to be the first such report on the dissociation energy of the N-H···N hydrogen bond.

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s^-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

Proving and Probing the Presence of the Elusive C-H⋅⋅⋅O Hydrogen Bond in Liquid Solutions at Room Temperature

Hydrogen bonds (H bonds) play a major role in defining the structure and properties of many substances, as well as phenomena and processes. Traditional H bonds are ubiquitous in nature, yet the demonstration of weak H bonds that occur between a highly polarized C-H group and an electron-rich oxygen atom, has proven elusive. Detailed here are linear and nonlinear IR spectroscopy experiments that reveal the presence of H bonds between the chloroform C-H group and an amide carbonyl oxygen atom in solution at room temperature. Evidence is provided for an amide solvation shell featuring two clearly distinguishable chloroform arrangements that undergo chemical exchange with a time scale of about 2 ps. Furthermore, the enthalpy of breaking the hydrogen bond is found to be 6-20 kJ mol^-1 . Ab-initio computations support the findings of two distinct solvation shells formed by three chloroform molecules, where one thermally undergoes hydrogen-bond making and breaking.

Structural and Functional Characterization of Sulfonium Carbon-Oxygen Hydrogen Bonding in the Deoxyamino Sugar Methyltransferase TylM1

The N-methyltransferase TylM1 from Streptomyces fradiae catalyzes the final step in the biosynthesis of the deoxyamino sugar mycaminose, a substituent of the antibiotic tylosin. The high-resolution crystal structure of TylM1 bound to the methyl donor S-adenosylmethionine (AdoMet) illustrates a network of carbon-oxygen (CH···O) hydrogen bonds between the substrate's sulfonium cation and residues within the active site. These interactions include hydrogen bonds between the methyl and methylene groups of the AdoMet sulfonium cation and the hydroxyl groups of Tyr14 and Ser120 in the enzyme. To examine the functions of these interactions, we generated Tyr14 to phenylalanine (Y14F) and Ser120 to alanine (S120A) mutations to selectively ablate the CH···O hydrogen bonding to AdoMet. The TylM1 S120A mutant exhibited a modest decrease in its catalytic efficiency relative to that of the wild type (WT) enzyme, whereas the Y14F mutation resulted in an approximately 30-fold decrease in catalytic efficiency. In contrast, site-specific substitution of Tyr14 by the noncanonical amino acid p-aminophenylalanine partially restored activity comparable to that of the WT enzyme. Correlatively, quantum mechanical calculations of the activation barrier energies of WT TylM1 and the Tyr14 mutants suggest that substitutions that abrogate hydrogen bonding with the AdoMet methyl group impair methyl transfer. Together, these results offer insights into roles of CH···O hydrogen bonding in modulating the catalytic efficiency of TylM1.

Understanding the structural and functional properties of carbohydrate esterases with a special focus on hemicellulose deacetylating acetyl xylan esterases

Acetyl and methyl esterifications are two major naturally found substitutions in the plant cell-wall polysaccharides. The non-cellulosic plant cell-wall polysaccharides such as pectin and hemicellulose are differentially esterified by the O-acetyl and methyl groups to cease the action of various hydrolytic enzymes secreted by different fungi and bacterial species. Thus, microorganisms have emerged with a special class of enzymes known as carbohydrate esterases (CE). The CE catalyse O-de, N-deacetylation of acetylated saccharide residues (esters or amides, where sugars play the role of alcohol/amine/acid). Carbohydrate active enzyme (CAZy) database has classified CE into 16 classes, of which hemicellulose deacetylating CE were grouped into eight classes (CE-1 to CE-7 and CE-16). Various plant biomass degrading fungi and bacteria secretes acetyl xylan esterases (AcXE); however, these enzymes exhibit varied substrate specificities. AcXE and xylanases-coupled pretreatment methods exhibit significant applications, such as enhancing animal feedstock, baking industry, production of food additives, paper and pulp, xylitol production and biorefinery-based industries, respectively. Thus, understanding the structural and functional properties of acetyl xylan esterase will significantly aid in developing the efficient AcXE with wide range of industrial applications.

Model molecules to classify CHO hydrogen-bonds

We developed a set of conformationally locked molecules each of which makes a single CHO H-bond/short contact and has different electron density at the acceptor oxygen atom. The downfield shift of the 1H NMR signals due to the hydrogen involved in the CHO H-bond varied from 0.93-1.6 ppm, and the magnitude of Δδ is in correlation with the hybridization state of the acceptor oxygen and with the CHO H-bond strengths quantified using a computational method.

An Intramolecular CAr-H•••O=C Hydrogen Bond and the Configuration of Rotenoids

Over the past half a century, the structure and configuration of the rotenoids, a group of natural products showing multiple promising bioactivities, have been established by interpretation of their NMR and electronic circular dichroism spectra and confirmed by analysis of single-crystal X-ray diffraction data. The chemical shift of the H-6' ¹H NMR resonance has been found to be an indicator of either a cis or trans C/D ring system. In the present study, four structures representing the central rings of a cis-, a trans-, a dehydro-, and an oxadehydro-rotenoid have been plotted using the Mercury program based on X-ray crystal structures reported previously, with the conformations of the C/D ring system, the local bond lengths or interatomic distances, hydrogen bond angles, and the H-6' chemical shift of these compounds presented. It is shown for the first time that a trans-fused C/D ring system of rotenoids is preferred for the formation of a potential intramolecular C_6'-H_6'•••O=C₄ H-bond, and that such H-bonding results in the ¹H NMR resonance for H-6' being shifted downfield.

Conformational Heterogeneity and the Role of the C(2)-H Donor in Mono- and Dihydrated Clusters of Benzoxazole

The significance of the heteroatom in influencing the acidity and binding affinity of the C(2)-H donor in five-membered heterocyclic rings is explored. The water clusters of benzoxazole (BOX) are studied in a supersonic jet by IR-UV double resonance spectroscopy and compared with those of benzimidazole (BIM) and its N-methyl derivative (MBIM). Two conformers are identified for the monohydrated cluster, both of which are O-H···N bound and differ in their orientation with respect to the intermolecular hydrogen bond. IR spectroscopy of the doubly hydrated cluster shows the presence of an intensity enhanced C(2)-H (carbon atom between the heteroatoms in the five-membered ring) stretching mode in addition to two red-shifted bound OH stretches, indicating that the water molecules form a hydrogen-bonded bridge encompassing the N acceptor and the weakly activated C(2)-H bond in oxazole. Comparison of the topological parameters of electron density, natural bond orbital analyses, and computed binding energies of the analogous hydrated structures of BOX, BIM, and MBIM indicates that the C(2)-H bond in the former is a more potent H-bond donor.

Unraveling the stereochemical and dynamic aspects of the catalytic site of bacterial peptidyl-tRNA hydrolase

Bacterial peptidyl-tRNA hydrolase (Pth; EC 3.1.1.29) hydrolyzes the peptidyl-tRNAs accumulated in the cytoplasm and thereby prevents cell death by alleviating tRNA starvation. X-ray and NMR studies of Vibrio cholerae Pth (VcPth) and mutants of its key residues involved in catalysis show that the activity and selectivity of the protein depends on the stereochemistry and dynamics of residues H24, D97, N118, and N14. D97-H24 interaction is critical for activity because it increases the nucleophilicity of H24. The N118 and N14 have orthogonally competing interactions with H24, both of which reduce the nucleophilicity of H24 and are likely to be offset by positioning of a peptidyl-tRNA substrate. The region proximal to H24 and the lid region exhibit slow motions that may assist in accommodating the substrate. Helix α3 exhibits a slow wobble with intermediate time scale motions of its N-cap residue N118, which may work as a flypaper to position the scissile ester bond of the substrate. Overall, the dynamics of interactions between the side chains of N14, H24, D97, and N118, control the catalysis of substrate by this enzyme.

Torsional and Electronic Factors Control the C-H⋅⋅⋅O Interaction

The precise role of non‐conventional hydrogen bonds such as the C−H⋅⋅⋅O interaction in influencing the conformation of small molecules remains unresolved. Here we survey a series of β‐turn mimetics using X‐ray crystallography and NMR spectroscopy in conjunction with quantum calculation, and conclude that favourable torsional and electronic effects are important for the population of states with conformationally influential C−H⋅⋅⋅O interactions. Our results also highlight the challenge in attempting to deconvolute a myriad of interdependent noncovalent interactions in order to focus on the contribution of a single one. Within a small molecule that is designed to resemble the complexity of the environment within peptides and proteins, the interplay of different steric burdens, hydrogen‐acceptor/‐donor properties and rotational profiles illustrate why unambiguous conclusions based solely on NMR chemical shift data are extremely challenging to rationalize.

Water bridges anchored by a C–H⋯O hydrogen bond: the role of weak interactions in molecular solvation

Hydrogen-bonded water bridges are re-directed from a polar NH bond to a weakly activated C(2)–H bond upon N-methylation. Infrared spectra, supported by ab initio calculations, provide direct evidence of the role of the C(2)–H donor in the solvation of the imidazole ring.

Conformational preferences of monohydrated clusters of imidazole derivatives revisited

IR-UV double resonance spectroscopy was used to identify the conformers of monohydrated benzimidazole andN-methylbenzimidazole in a supersonic jet. A new OH–N bound conformer relevant to histidine containing proteins was discovered. The long standing differences in the literature about the relative energies and abundance of the monohydrated imidazole derivatives have also been resolved.

ω-Turn: a novel β-turn mimic in globular proteins stabilized by main-chain to side-chain C−H···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10-membered hydrogen-bonded ring involving the main-chain C − O in a β-turn can be formed using a side-chain carbonyl group leading to Asx-turn. We show that the N − H component of hydrogen bond can be replaced by a C(γ) -H group in the side chain, culminating in a nonconventional C − H···O interaction. Because of its shape this β-turn mimic is designated as ω-turn, which is found to occur ∼ three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C − H···O interaction occurring between the terminal residues, constraining the torsion angles ϕi + 1, ψi + 1, ϕi + 2 and χ'1(i + 2) (using the interacting C(γ) atom). Based on these angles there are two types of ω-turns, each of which can be further divided into two groups. C(β) -branched side-chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal-binding sites. N-linked glycosylation occurs at the consensus pattern Asn-Xaa-Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω-turn, which may be the recognition site for protein modification. Location between two β-strands is the most common occurrence in protein tertiary structure, and being generally exposed ω-turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design.

Mechanism of Orlistat Hydrolysis by the Thioesterase of Human Fatty Acid Synthase

Fatty acid synthase (FASN), the sole protein capable of de novo synthesis of free fatty acids, is overexpressed in a wide variety of human cancers and is associated with poor prognosis and aggressiveness of these cancers. Orlistat, an FDA-approved drug for obesity treatment that inhibits pancreatic lipases in the GI tract, also inhibits the thioesterase (TE) of human FASN. The cocrystal structure of TE with orlistat shows a pseudo TE dimer containing two different forms of orlistat in the active site, an intermediate that is covalently bound to a serine residue (Ser²³⁰⁸) and a hydrolyzed and inactivated product. In this study, we attempted to understand the mechanism of TE-catalyzed orlistat hydrolysis by examining the role of the hexyl tail of the covalently bound orlistat in water activation for hydrolysis using molecular dynamics simulations. We found that the hexyl tail of the covalently bound orlistat undergoes a conformational transition, which is accompanied by destabilization of a hydrogen bond between a hydroxyl moiety of orlistat and the catalytic His²⁴⁸¹ of TE that in turn leads to an increased hydrogen bonding between water molecules and His²⁴⁸¹ and increased chance for water activation to hydrolyze the covalent bond between orlistat and Ser²³⁰⁸. Thus, the conformation of the hexyl tail of orlistat plays an important role in orlistat hydrolysis. Strategies that stabilize the hexyl tail may lead to the design of more potent irreversible inhibitors that target FASN and block TE activity with greater endurance.

Cysteine-mediated dynamic hydrogen-bonding network in the active site of Pin1

Enzymes catalyze a plethora of chemical reactions that are tightly regulated and intricately coupled in biology. Catalysis of phosphorylation-dependent cis-trans isomerization of peptidyl-prolyl bonds, which act as conformational switches in regulating many post-phosphorylation processes, is considered to be one of the most critical. Pin1 is a cis-trans isomerase of peptidyl-prolyl(ω-) bonds of phosphorylated-Ser/Thr-Pro motifs and has been implicated in many diseases. Structural and experimental studies are still unable to resolve the mechanistic role and protonation states of two adjacent histidines (His59 and His157) and a cysteine (Cys113) in the active site of Pin1. Here, we show that the protonation state of Cys113 mediates a dynamic hydrogen-bonding network in the active site of Pin1, involving the two adjacent histidines and several other residues that are highly conserved and necessary for catalysis. We have used detailed free energy calculations and molecular dynamics simulations, complementing previous experiments, to resolve the ambiguities in the orientations of the histidines and protonation states of these key active site residues, details that are critical for fully understanding the mechanism of Pin1 and necessary for developing potent inhibitors. Importantly, Cys113 is shown to alternate between the unprotonated and neutral states, unprotonated in free Pin1 and neutral in substrate-bound Pin1. Our results are consistent with experiments and provide an explanation for the chemical reactivity of free Pin1 that is suggested to be necessary for the regulation of the enzyme.

Complete backbone and DENQ side chain NMR assignments in proteins from a single experiment: implications to structure-function studies

Resonance assignment is the first and the most crucial step in all nuclear magnetic resonance (NMR) investigations on structure-function relationships in biological macromolecules. Often, the assignment exercise has to be repeated several times when specific interactions with ligands, substrates etc., have to be elucidated for understanding the functional mechanisms. While the protein backbone serves to provide a scaffold, the side chains interact directly with the ligands. Such investigations will be greatly facilitated, if there are rapid methods for obtaining exhaustive information with minimum of NMR experimentation. In this context, we present here a pulse sequence which exploits the recently introduced technique of parallel detection of multiple nuclei, e.g. (1)H and (13)C, and results in two 3D-data sets simultaneously. These yield complete backbone resonance assignment ((1)H(N), (15)N, (13)CO, (1)Hα/(13)Cα, and (1)Hβ/(13)Cβ chemical shifts) and side chain assignment of D, E, N and Q residues. Such an exhaustive assignment has the potential of yielding accurate 3D structures using one or more of several algorithms which calculate structures of the molecules very reliably on the basis of NMR chemical shifts alone. The side chain assignments of D, E, N, and Q will be extremely valuable for interaction studies with different ligands; D and E side chains are known to be involved in majority of catalytic activities. Utility of this experiment has been demonstrated with Ca(2+) bound M-crystallin, which contains largely D, E, N and Q residues at the metal binding sites.

Structure of a proteolytically resistant analogue of (NLys)5SFTI-1 in complex with trypsin: evidence for the direct participation of the Ser214 carbonyl group in serine protease-mediated proteolysis

Peptide-peptoid hybrids are found to be potent inhibitors of serine proteases. These engineered peptidomimetics benefit from both types of units of the biopolymeric structure: the natural inhibitor part serves as a good binding template, while the P1-positioned peptoid component provides complete resistance towards proteolysis. In this report, the mechanism of proteolytic resistance of a P1 peptoid-containing analogue is postulated based on the crystal structure of the (NLys)(5)-modified sunflower trypsin inhibitor SFTI-1 in complex with bovine trypsin solved at 1.29 Å resolution. The structural differences between the (NLys)(5)SFTI-1-trypsin complex and the native SFTI-1-trypsin complex are surprisingly small and reveal the key role of the carbonyl group of the Ser214 residue of the enzyme, which is crucial for binding of the inhibitor and plays a crucial role in proteolysis mediated by serine proteases. The incorporated NLys5 peptoid residue prevents Ser214 from forming a hydrogen bond to the P1 residue, and in turn Gln192 does not form a hydrogen bond to the carbonyl group of the P2 residue. It also increases the distance between the Ser214 carbonyl group and the Ser195 residue, thus preventing proteolysis. The hybrid inhibitor structure reported here provides insight into protein-protein interaction, which can be efficiently and selectively probed with the use of peptoids incorporated within endogenous peptide ligands.

C–H/O interactions of nucleic bases with a water molecule: a crystallographic and quantum chemical study

The C–H/O interactions of nucleic bases are substantially stronger than the C–H/O interactions of benzene and pyridine. These results can be very important for molecular recognition of DNA and RNA.

Conservation and functional importance of carbon-oxygen hydrogen bonding in AdoMet-dependent methyltransferases

S-adenosylmethionine (AdoMet)-based methylation is integral to metabolism and signaling. AdoMet-dependent methyltransferases belong to multiple distinct classes and share a catalytic mechanism that arose through convergent evolution; however, fundamental determinants underlying this shared methyl transfer mechanism remain undefined. A survey of high-resolution crystal structures reveals that unconventional carbon-oxygen (CH···O) hydrogen bonds coordinate the AdoMet methyl group in different methyltransferases irrespective of their class, active site structure, or cofactor binding conformation. Corroborating these observations, quantum chemistry calculations demonstrate that these charged interactions formed by the AdoMet sulfonium cation are stronger than typical CH···O hydrogen bonds. Biochemical and structural studies using a model lysine methyltransferase and an active site mutant that abolishes CH···O hydrogen bonding to AdoMet illustrate that these interactions are important for high-affinity AdoMet binding and transition-state stabilization. Further, crystallographic and NMR dynamics experiments of the wild-type enzyme demonstrate that the CH···O hydrogen bonds constrain the motion of the AdoMet methyl group, potentially facilitating its alignment during catalysis. Collectively, the experimental findings with the model methyltransferase and structural survey imply that methyl CH···O hydrogen bonding represents a convergent evolutionary feature of AdoMet-dependent methyltransferases, mediating a universal mechanism for methyl transfer.

First steps in growth of a polypeptide toward β-sheet structure

The full conformational energy surface is examined for a molecule in which a dipeptide is attached to the same spacer group as another peptide chain, so as to model the seminal steps of β-sheet formation. This surface is compared with the geometrical preferences of the isolated dipeptide to extract the perturbations induced by interactions with the second peptide strand. These interpeptide interactions remove any tendency of the dipeptide to form a C5 ring structure, one of its two normally stable geometries. A C7 structure, the preferred conformation of the isolated dipeptide, remains as the global minimum in the full molecule. However, the stability of this structure is highly dependent upon interpeptide H-bonds with the second chain. The latter forces include not only the usual NH···O interaction, but also a pair of CH···O H-bonds. The secondary minimum is also of C7 type and likewise depends in part upon CH···O H-bonds for its stability. The latter interactions also play a part in the tertiary minimum. A two-strand β-sheet structure is not yet in evidence for this small model system, requiring additional peptide units to be added to each chain.

The structure of the karrikin-insensitive protein (KAI2) in Arabidopsis thaliana

KARRIKIN INSENSITIVE 2 (KAI2) is an α/β hydrolase involved in seed germination and seedling development. It is essential for plant responses to karrikins, a class of butenolide compounds derived from burnt plant material that are structurally similar to strigolactone plant hormones. The mechanistic basis for the function of KAI2 in plant development remains unclear. We have determined the crystal structure of Arabidopsis thaliana KAI2 in space groups P2₁ 2₁ 2₁ (a  = 63.57 Å, b  = 66.26 Å, c  = 78.25 Å) and P2₁ (a  = 50.20 Å, b  = 56.04 Å, c  = 52.43 Å, β  = 116.12°) to 1.55 and 2.11 Å respectively. The catalytic residues are positioned within a large hydrophobic pocket similar to that of DAD2, a protein required for strigolactone response in Petunia hybrida. KAI2 possesses a second solvent-accessible pocket, adjacent to the active site cavity, which offers the possibility of allosteric regulation. The structure of KAI2 is consistent with its designation as a serine hydrolase, as well as previous data implicating the protein in karrikin and strigolactone signalling.

Carbon-oxygen hydrogen bonding in biological structure and function

Carbon-oxygen (CH···O) hydrogen bonding represents an unusual category of molecular interactions first documented in biological structures over 4 decades ago. Although CH···O hydrogen bonding has remained generally underappreciated in the biochemical literature, studies over the last 15 years have begun to yield direct evidence of these interactions in biological systems. In this minireview, we provide a historical context of biological CH···O hydrogen bonding and summarize some major advancements from experimental studies over the past several years that have elucidated the importance, prevalence, and functions of these interactions. In particular, we examine the impact of CH···O bonds on protein and nucleic acid structure, molecular recognition, and enzyme catalysis and conclude by exploring overarching themes and unresolved questions regarding unconventional interactions in biomolecular structure.

Functional and structural characterization of a thermostable acetyl esterase from Thermotoga maritima

TM0077 from Thermotoga maritima is a member of the carbohydrate esterase family 7 and is active on a variety of acetylated compounds, including cephalosporin C. TM0077 esterase activity is confined to short-chain acyl esters (C2-C3), and is optimal around 100°C and pH 7.5. The positional specificity of TM0077 was investigated using 4-nitrophenyl-β-D-xylopyranoside monoacetates as substrates in a β-xylosidase-coupled assay. TM0077 hydrolyzes acetate at positions 2, 3, and 4 with equal efficiency. No activity was detected on xylan or acetylated xylan, which implies that TM0077 is an acetyl esterase and not an acetyl xylan esterase as currently annotated. Selenomethionine-substituted and native structures of TM0077 were determined at 2.1 and 2.5 Å resolution, respectively, revealing a classic α/β-hydrolase fold. TM0077 assembles into a doughnut-shaped hexamer with small tunnels on either side leading to an inner cavity, which contains the six catalytic centers. Structures of TM0077 with covalently bound phenylmethylsulfonyl fluoride and paraoxon were determined to 2.4 and 2.1 Å, respectively, and confirmed that both inhibitors bind covalently to the catalytic serine (Ser188). Upon binding of inhibitor, the catalytic serine adopts an altered conformation, as observed in other esterase and lipases, and supports a previously proposed catalytic mechanism in which Ser hydroxyl rotation prevents reversal of the reaction and allows access of a water molecule for completion of the reaction.

Crystal structures of enterovirus 71 3C protease complexed with rupintrivir reveal the roles of catalytically important residues

EV71 is the primary pathogenic cause of hand-foot-mouth disease (HFMD), but an effective antiviral drug currently is unavailable. Rupintrivir, an inhibitor against human rhinovirus (HRV), has potent antiviral activities against EV71. We determined the high-resolution crystal structures of the EV71 3C(pro)/rupintrivir complex, showing that although rupintrivir interacts with EV71 3C(pro) similarly to HRV 3C(pro), the C terminus of the inhibitor cannot accommodate the leaving-group pockets of EV71 3C(pro). Our structures reveal that EV71 3C(pro) possesses a surface-recessive S2' pocket that is not present in HRV 3C(pro) that contributes to the additional substrate binding affinity. Combined with mutagenic studies, we demonstrated that catalytic Glu71 is irreplaceable for maintaining the overall architecture of the active site and, most importantly, the productive conformation of catalytic His40. We discovered the role of a previously uncharacterized residue, Arg39 of EV71 3C(pro), that can neutralize the negative charge of Glu71, which may subsequently assist deprotonation of His40 during proteolysis.

Direct evidence for methyl group coordination by carbon-oxygen hydrogen bonds in the lysine methyltransferase SET7/9

SET domain lysine methyltransferases (KMTs) are S-adenosylmethionine (AdoMet)-dependent enzymes that catalyze the site-specific methylation of lysyl residues in histone and non-histone proteins. Based on crystallographic and cofactor binding studies, carbon-oxygen (CH · · · O) hydrogen bonds have been proposed to coordinate the methyl groups of AdoMet and methyllysine within the SET domain active site. However, the presence of these hydrogen bonds has only been inferred due to the uncertainty of hydrogen atom positions in x-ray crystal structures. To experimentally resolve the positions of the methyl hydrogen atoms, we used NMR (1)H chemical shift coupled with quantum mechanics calculations to examine the interactions of the AdoMet methyl group in the active site of the human KMT SET7/9. Our results indicated that at least two of the three hydrogens in the AdoMet methyl group engage in CH · · · O hydrogen bonding. These findings represent direct, quantitative evidence of CH · · · O hydrogen bond formation in the SET domain active site and suggest a role for these interactions in catalysis. Furthermore, thermodynamic analysis of AdoMet binding indicated that these interactions are important for cofactor binding across SET domain enzymes.