skip to main content

Title: When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins
Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural ( i.e. , X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups ( e.g. , hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic more » benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Chemical Science
Page Range or eLocation-ID:
1147 to 1162
Sponsoring Org:
National Science Foundation
More Like this
  1. Hydrogen bonds (HB)s are the most abundant motifs in biological systems. They play a key role in determining protein–ligand binding affinity and selectivity. We designed two pharmaceutically beneficial HB databases, database A including ca. 12,000 protein–ligand complexes with ca. 22,000 HBs and their geometries, and database B including ca. 400 protein–ligand complexes with ca. 2200 HBs, their geometries, and bond strengths determined via our local vibrational mode analysis. We identified seven major HB patterns, which can be utilized as a de novo QSAR model to predict the binding affinity for a specific protein–ligand complex. Glycine was reported as the most abundant amino acid residue in both donor and acceptor profiles, and N–H⋯O was the most frequent HB type found in database A. HBs were preferred to be in the linear range, and linear HBs were identified as the strongest. HBs with HB angles in the range of 100–110°, typically forming intramolecular five-membered ring structures, showed good hydrophobic properties and membrane permeability. Utilizing database B, we found a generalized Badger’s relationship for more than 2200 protein–ligand HBs. In addition, the strength and occurrence maps between each amino acid residue and ligand functional groups open an attractive possibility for a novel drug-designmore »approach and for determining drug selectivity and affinity, and they can also serve as an important tool for the hit-to-lead process.« less
  2. Abstract

    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,more »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.

    « less
  3. The sodium/potassium-ATPase (NKA) is the enzyme that establishes gradients of sodium and potassium across the plasma membrane. NKA activity is tightly regulated for different physiological contexts through interactions with single-span transmembrane peptides, the FXYD proteins. This diverse family of regulators has in common a domain containing a Phe-X-Tyr-Asp (FXYD) motif, two conserved glycines, and one serine residue. In humans, there are seven tissue-specific FXYD proteins that differentially modulate NKA kinetics as appropriate for each system, providing dynamic responsiveness to changing physiological conditions. Our understanding of how FXYD proteins contribute to homeostasis has benefitted from recent advances described in this review: biochemical and biophysical studies have provided insight into regulatory mechanisms, genetic models have uncovered remarkable complexity of FXYD function in integrated physiological systems, new posttranslational modifications have been identified, high-resolution structural studies have revealed new details of the regulatory interaction with NKA, and new clinical correlations have been uncovered. In this review, we address the structural determinants of diverse FXYD functions and the special roles of FXYDs in various physiological systems. We also discuss the possible roles of FXYDs in protein trafficking and regulation of non-NKA targets.

  4. Hemoglobins (Hbs) of crocodilians are reportedly characterized by unique mechanisms of allosteric regulatory control, but there are conflicting reports regarding the importance of different effectors, such as chloride ions, organic phosphates, and CO 2 . Progress in understanding the unusual properties of crocodilian Hbs has also been hindered by a dearth of structural information. Here, we present the first comparative analysis of blood properties and Hb structure and function in a phylogenetically diverse set of crocodilian species. We examine mechanisms of allosteric regulation in the Hbs of 13 crocodilian species belonging to the families Crocodylidae and Alligatoridae. We also report new amino acid sequences for the α- and β-globins of these taxa, which, in combination with structural analyses, provide insights into molecular mechanisms of allosteric regulation. All crocodilian Hbs exhibited a remarkably strong sensitivity to CO 2 , which would permit effective O 2 unloading to tissues in response to an increase in metabolism during intense activity and diving. Although the Hbs of all crocodilians exhibit similar intrinsic O 2 -affinities, there is considerable variation in sensitivity to Cl − ions and ATP, which appears to be at least partly attributable to variation in the extent of NH 2 -terminalmore »acetylation. Whereas chloride appears to be a potent allosteric effector of all crocodile Hbs, ATP has a strong, chloride-independent effect on Hb-O 2 affinity only in caimans. Modeling suggests that allosteric ATP binding has a somewhat different structural basis in crocodilian and mammalian Hbs.« less
  5. Abstract

    The ground‐state energy and density of 4 low‐energy conformations of the formic acid dimer were calculated via partition density functional theory (PDFT). The differences between isolated and PDFT monomer densities display similar deformation patterns for primary and secondary hydrogen bonds (HBs) among all 4 dimers. In contrast, the partition potential shows no transferable features in the bonding regions. These observations highlight the global character of the partition potential and the cooperative effect that occurs when a dimer is bound via more than 1 HB. We also provide numerical confirmation of the intuitive (but unproven) observation that fragment deformation energies are larger for systems with larger binding energies.