HBAT: a complete package for analysing strong and weak hydrogen bonds in macromolecular crystal structures

The program HBAT is a tool to automate the analysis of potential hydrogen bonds and similar type of weak interactions like halogen bonds and non-canonical interactions in macromolecular structures, available in Brookhaven Protein Database (PDB) file format. HBAT is written using PERL and TK languages. The program generates an MSOFFICE Excel compatible output file for statistical analysis. HBAT identify potential interactions based on geometrical criteria. A series of analysis reports like frequency tables, geometry distribution tables, furcations list are generated. A user friendly GUI offers freedom to select several parameters and options. Graphviz based visualization of hydrogen bond networks in 2D helps to study the cooperativity and anticooperativity geometry in hydrogen bond. HBAT supports post docking interaction analysis between PDB files for any target/receptor (in PDB files) and docked ligands/poses (in SDF). This tool can be implemented in active site interaction analysis, structure based drug design and molecular dynamics simulations.     

Geometric cooperativity and anticooperativity of three-body interactions in native proteins

Characterizing multibody interactions of hydrophobic, polar, and ionizable residues in protein is important for understanding the stability of protein structures. We introduce a geometric model for quantifying 3-body interactions in native proteins. With this model, empirical propensity values for many types of 3-body interactions can be reliably estimated from a database of native protein structures, despite the overwhelming presence of pairwise contacts. In addition, we define a nonadditive coefficient that characterizes cooperativity and anticooperativity of residue interactions in native proteins by measuring the deviation of 3-body interactions from 3 independent pairwise interactions. It compares the 3-body propensity value from what would be expected if only pairwise interactions were considered, and highlights the distinction of propensity and cooperativity of 3-body interaction. Based on the geometric model, and what can be inferred from statistical analysis of such a model, we find that hydrophobic interactions and hydrogen-bonding interactions make nonadditive contributions to protein stability, but the nonadditive nature depends on whether such interactions are located in the protein interior or on the protein surface. When located in the interior, many hydrophobic interactions such as those involving alkyl residues are anticooperative. Salt-bridge and regular hydrogen-bonding interactions, such as those involving ionizable residues and polar residues, are cooperative. When located on the protein surface, these salt-bridge and regular hydrogen-bonding interactions are anticooperative, and hydrophobic interactions involving alkyl residues become cooperative. We show with examples that incorporating 3-body interactions improves discrimination of protein native structures against decoy conformations. In addition, analysis of cooperative 3-body interaction may reveal spatial motifs that can suggest specific protein functions.     

Detection of hydrogen-bond signature patterns in protein families

SUMMARY: We have developed a WWW server, HBPRINT, for the identification of hydrogen-bond signature patterns in protein families from their structures. The server calculates (a) common hydrogen bonds between two structures (b) a hydrogen-bond fingerprint in a set of structural neighbours and (c) details of conserved hydrogen bonds. The server also enables the visualization of the hydrogen bond network comprising the signature pattern. AVAILABILITY: HBPRINT and a tutorial are available from http://144.16.93.115/hb_page/index.html.     

Hydrogen-bond cooperativity in protein secondary structure

Hydrogen bonding in the α-helix and β-sheet has been studied by ab initio molecular orbital calculations carried out on complexes of formamide. Hydrogen-bond geometries were taken from x-ray crystallography of polypeptides. Positive cooperativity is found in all cases. The limiting value for infinite chains is obtained by use of a double-reciprocal plot and indicates an increase in the effective bond strength of 25% over that of a single isolated bond. Parallel calculations based on a classical electrostatic model yield qualitatively similar trends but underestimate the cooperativity by half. Charge redistribution accompanying cooperativity is characterized by a new type of charge-density difference plot, the cooperativity map. The magnitude and distance over which cooperativity acts suggest several significant biological consequences. Thus the average of α-helices and the number of β-sheet strands found in protein may be influenced by cooperativity. Cooperativity in the interpeptide hydrogen bond may also be partly responsible for the rapid formation of secondary structure in renaturing proteins and help stabilize secondary structure relative to the random-coil conformation.     

Protein Structure Validation and Refinement Using Amide Proton Chemical Shifts Derived from Quantum Mechanics

We present the ProCS method for the rapid and accurate prediction of protein backbone amide proton chemical shifts - sensitive probes of the geometry of key hydrogen bonds that determine protein structure. ProCS is parameterized against quantum mechanical (QM) calculations and reproduces high level QM results obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is interfaced with the PHAISTOS protein simulation program and is used to infer statistical protein ensembles that reflect experimentally measured amide proton chemical shift values. Such chemical shift-based structural refinements, starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and trans-hydrogen bond (^h3 J_NC'') spin-spin coupling constants that are in excellent agreement with experiment. We show that the structural sensitivity of the QM-based amide proton chemical shift predictions is needed to obtain this agreement. The ProCS method thus offers a powerful new tool for refining the structures of hydrogen bonding networks to high accuracy with many potential applications such as protein flexibility in ligand binding.     

Theoretical investigation into the cooperativity effect between the intermolecular π∙π and H-bonding interactions in the curcumin∙cytosine∙H<Subscript>2</Subscript>O system

In order to reveal the mechanism of drug action and design of DNA/RNA-targeted drugs containing aromatic rings, the cooperativity effects between the intermolecular π???π and H-bonding interactions in curcumin(drug)???cytosine(DNA/RNA base)???H₂O were investigated by the B3LYP-D3 and MP2(full) methods with the 6–311++G(2d,p) basis set. The π???π interaction plays an important role in stabilizing the linear ternary complexes with the cooperativity effects, and the cyclic structures suffer the anticooperativity effects. The cooperativity or anticooperativity effects are notable, which could lead to a possible significant change in drug activity. The hydration is essentially the cooperativity or anticooperativity effect. These results were confirmed by the atoms in molecules (AIM), reduced density gradient (RDG), and surface electrostatic potentials analyses. The cyclic complexes are more stable, from which it can be deduced that the drug always links with the DNA/RNA base and H₂O by the π???π or H-bonding interactions, and only in this way can the drug activity be shown. Therefore, the designed DNA/RNA-targeted drugs should possess a certain number of hydrophilic groups in contact with the DNA/RNA base and H₂O to reconcile drug activity by the cooperativity effect between the π???π and H-bonding interactions, as is in agreement with many of the drugs in use.

Open image in new window

Graphical abstract RDG isosurface of ternary complex

     

ProSAT: functional annotation of protein 3D structures

ProSAT (for Protein Structure Annotation Tool) is a tool to facilitate interactive visualization of non-structure-based functional annotations in protein 3D structures. It performs automated mapping of the functional annotations onto the protein structure and allows functional sites to be readily identified upon visualization. The current version of ProSAT can be applied to large datasets of protein structures for fast visual identification of active and other functional sites derived from the SwissProt and Prosite databases.     

NMR of hydrogen bonding in cold-shock protein A and an analysis of the influence of crystallographic resolution on comparisons of hydrogen bond lengths

Hydrogen bonding in cold-shock protein A of Escherichia coli has been investigated using long-range HNCO spectroscopy. Nearly half of the amide protons involved in hydrogen bonds in solution show no measurable protection from exchange in water, cautioning against a direct correspondence between hydrogen bonding and hydrogen exchange protection. The N to O atom distance across a hydrogen bond, R(NO), is related to the size of the (3h)J(NC') trans hydrogen bond coupling constant and the amide proton chemical shift. Both NMR parameters show poorer agreement with the 2.0-A resolution X-ray structure of the cold-shock protein studied by NMR than with a 1.2-A resolution X-ray structure of a homologous cold-shock protein from the thermophile B. caldolyticus. The influence of crystallographic resolution on comparisons of hydrogen bond lengths was further investigated using a database of 33 X-ray structures of ribonuclease A. For highly similar structures, both hydrogen bond R(NO) distance and Calpha coordinate root mean square deviations (RMSD) show systematic increases as the resolution of the X-ray structure used for comparison decreases. As structures diverge, the effects of coordinate errors on R(NO) distance and Calpha coordinate root mean square deviations become progressively smaller. The results of this study are discussed with regard to the influence of data precision on establishing structure similarity relationships between proteins.     

Positioning hydrogen atoms by optimizing hydrogen-bond networks in protein structures

A method is presented that positions polar hydrogen atoms in protein structures by optimizing the total hydrogen bond energy. For this goal, an empirical hydrogen bond force field was derived from small molecule crystal structures. Bifurcated hydrogen bonds are taken into account. The procedure also predicts ionization states of His, Asp, and Glu residues. During optimization, sidechain conformations of His, Gln, and Asn residues are allowed to change their last χ angle by 180° to compensate for crystallographic misassignments. Crystal structure symmetry is taken into account where appropriate. The results can have significant implications for molecular dynamics simulations, protein engineering, and docking studies. The largest impact, however, is in protein structure verification: over 85% of protein structures tested can be improved by using our procedure. Proteins 26:363–376 © 1996 Wiley-Liss, Inc.     

Steroid ring hydroxylation patterns govern cooperativity in human bile acid binding protein

Human ileal bile acid binding protein (I-BABP) is a member of the intracellular lipid binding protein family. This protein is thought to function in the transcellular transport and enterohepatic circulation of bile salts. Human I-BABP binds two molecules of glycocholate, the physiologically most abundant bile salt, with modest intrinsic affinity but a remarkably high degree of positive cooperativity. Here we report a calorimetric analysis for the binding of a broad panel of bile salts to human I-BABP. The interaction of I-BABP with nine physiologically relevant derivatives of cholic acid, chenodeoxycholic acid, and deoxycholic acid in their conjugated (glycine and taurine) and unconjugated forms was monitored by isothermal titration calorimetry. All bile salts bound to I-BABP with a 2:1 stoichiometry and similar overall affinity, but the derivatives of cholic acid displayed much higher Hill coefficients, a measure of macroscopic positive cooperativity. To test whether the cooperativity was dependent on individual structural features of the bile salt side chain, a series of side-chain-extended bile salts that lacked a hydrogen bond donor or acceptor at C-24 were chemically synthesized. These synthetic variants exhibited the same energetic and cooperativity profile as the naturally occurring bile salts. Our findings indicate that cooperativity in bile salt-I-BABP recognition is governed by the pattern of steroid B- and C-ring hydroxylation and not the presence or type of side-chain conjugation.     

On the evolutionary conservation of hydrogen bonds made by buried polar amino acids: the hidden joists,braces and trusses of protein architecture

Background  

The hydrogen bond patterns between mainchain atoms in protein structures not only give rise to regular secondary structures but also satisfy mainchain hydrogen bond potential. However, not all mainchain atoms can be satisfied through hydrogen bond interactions that arise in regular secondary structures; in some locations sidechain-to-mainchain hydrogen bonds are required to provide polar group satisfaction. Buried polar residues that are hydrogen-bonded to mainchain amide atoms tend to be highly conserved within protein families, confirming that mainchain architecture is a critical restraint on the evolution of proteins. We have investigated the stabilizing roles of buried polar sidechains on the backbones of protein structures by performing an analysis of solvent inaccessible residues that are entirely conserved within protein families and superfamilies and hydrogen bonded to an equivalent mainchain atom in each family member.     

An improved hydrogen bond potential: impact on medium resolution protein structures

A new semi-empirical force field has been developed to describe hydrogen-bonding interactions with a directional component. The hydrogen bond potential supports two alternative target angles, motivated by the observation that carbonyl hydrogen bond acceptor angles have a bimodal distribution. It has been implemented as a module for a macromolecular refinement package to be combined with other force field terms in the stereochemically restrained refinement of macromolecules. The parameters for the hydrogen bond potential were optimized to best fit crystallographic data from a number of protein structures. Refinement of medium-resolution structures with this additional restraint leads to improved structure, reducing both the free R-factor and over-fitting. However, the improvement is seen only when stringent hydrogen bond selection criteria are used. These findings highlight common misconceptions about hydrogen bonding in proteins, and provide explanations for why the explicit hydrogen bonding terms of some popular force field sets are often best switched off.     

Learning about protein hydrogen bonding by minimizing contrastive divergence

Defining the strength and geometry of hydrogen bonds in protein structures has been a challenging task since early days of structural biology. In this article, we apply a novel statistical machine learning technique, known as contrastive divergence, to efficiently estimate both the hydrogen bond strength and the geometric characteristics of strong interpeptide backbone hydrogen bonds, from a dataset of structures representing a variety of different protein folds. Despite the simplifying assumptions of the interatomic energy terms used, we determine the strength of these hydrogen bonds to be between 1.1 and 1.5 kcal/mol, in good agreement with earlier experimental estimates. The geometry of these strong backbone hydrogen bonds features an almost linear arrangement of all four atoms involved in hydrogen bond formation. We estimate that about a quarter of all hydrogen bond donors and acceptors participate in these strong interpeptide hydrogen bonds.     

Strong and weak hydrogen bonds in drug-DNA complexes: A statistical analysis

A statistical analysis of strong and weak hydrogen bonds in the minor groove of DNA was carried out for a set of 70 drug-DNA complexes. The terms ‘strong’ and ‘weak’ pertain to the inherent strengths and weakness of the donor and acceptor fragments rather than to any energy considerations. The dataset was extracted from the protein data bank (PDB). The analysis was performed with an in-house software, hydrogen bond analysis tool (HBAT). In addition to strong hydrogen bonds such as O—H⋯O and N—H⋯O, the ubiquitous presence of weak hydrogen bonds such as C—H⋯O is implicated in molecular recognition. On an average, there are 1.4 weak hydrogen bonds for every strong hydrogen bond. For both categories of interaction, the N(3) of purine and the O(2) of pyrimidine are favoured acceptors. Donor multifurcation is common with the donors generally present in the drug molecules, and shared by hydrogen bond acceptors in the minor groove. Bifurcation and trifurcation are most commonly observed. The metrics for strong hydrogen bonds are consistent with established trends. The geometries are variable for weak hydrogen bonds. A database of recognition geometries for 26 literature amidinium-based inhibitors of Human African Trypanosomes (HAT) was generated with a docking study using seven inhibitors which occur in published crystal structures included in the list of 70 complexes mentioned above, and 19 inhibitors for which the drug-DNA complex crystal structures are unknown. The virtual geometries so generated correlate well with published activities for these 26 inhibitors, justifying our assumption that strong and weak hydrogen bonds are optimized in the active site.     

PDBToSDF: Create ligand structure files from PDB file

Protein Data Bank (PDB) file contains atomic data for protein and ligand in protein-ligand complexes. Structure data file (SDF) contains data for atoms, bonds, connectivity and coordinates of molecule for ligands. We describe PDBToSDF as a tool to separate the ligand data from pdb file for the calculation of ligand properties like molecular weight, number of hydrogen bond acceptors, hydrogen bond receptors easily.     

Strong Ionic Hydrogen Bonding Causes a Spectral Isotope Effect in Photoactive Yellow Protein

Standard hydrogen bonds are of great importance for protein structure and function. Ionic hydrogen bonds often are significantly stronger than standard hydrogen bonds and exhibit unique properties, but their role in proteins is not well understood. We report that hydrogen/deuterium exchange causes a redshift in the visible absorbance spectrum of photoactive yellow protein (PYP). We expand the range of interpretable isotope effects by assigning this spectral isotope effect (SIE) to a functionally important hydrogen bond at the active site of PYP. The inverted sign and extent of this SIE is explained by the ionic nature and strength of this hydrogen bond. These results show the relevance of ionic hydrogen bonding for protein active sites, and reveal that the inverted SIE is a novel, to our knowledge, tool to probe ionic hydrogen bonds. Our results support a classification of hydrogen bonds that distinguishes the properties of ionic hydrogen bonds from those of both standard and low barrier hydrogen bonds, and show how this classification helps resolve a recent debate regarding active site hydrogen bonding in PYP.     

Strong Ionic Hydrogen Bonding Causes a Spectral Isotope Effect in Photoactive Yellow Protein

Standard hydrogen bonds are of great importance for protein structure and function. Ionic hydrogen bonds often are significantly stronger than standard hydrogen bonds and exhibit unique properties, but their role in proteins is not well understood. We report that hydrogen/deuterium exchange causes a redshift in the visible absorbance spectrum of photoactive yellow protein (PYP). We expand the range of interpretable isotope effects by assigning this spectral isotope effect (SIE) to a functionally important hydrogen bond at the active site of PYP. The inverted sign and extent of this SIE is explained by the ionic nature and strength of this hydrogen bond. These results show the relevance of ionic hydrogen bonding for protein active sites, and reveal that the inverted SIE is a novel, to our knowledge, tool to probe ionic hydrogen bonds. Our results support a classification of hydrogen bonds that distinguishes the properties of ionic hydrogen bonds from those of both standard and low barrier hydrogen bonds, and show how this classification helps resolve a recent debate regarding active site hydrogen bonding in PYP.     

New plot showing the existence of cooperativity, anticooperativity, and an arbitrary value of the neighbor-exclusion parameter in systems of ligand binding to linear biopolymers

A new way of plotting the isotherm for ligand binding to linear biopolymers is presented. In this plot the isotherm for noncooperative binding of ligands of length n (n-mers) becomes a straight line and the existence of cooperativity and anticooperativity between bound ligands is detected by appearance of opposite convexity of the curved isotherm. It is also usable for cases of fractional n values. Usefulness of the new plot in determining precise mechanisms of binding is shown using experimental data.     

Contribution of hydrogen bonds to the conformational stability of human lysozyme: calorimetry and X-ray analysis of six Ser --> Ala mutants.

To further examine the contribution of hydrogen bonds to the conformational stability of the human lysozyme, six Ser to Ala mutants were constructed. The thermodynamic parameters for denaturation of these six Ser mutant proteins were investigated by differential scanning calorimetry (DSC), and the crystal structures were determined by X-ray analysis. The denaturation Gibbs energy (DeltaG) of the Ser mutant proteins was changed from 2.0 to -5.7 kJ/mol, compared to that of the wild-type protein. With an analysis in which some factors that affected the stability due to mutation were considered, the contribution of hydrogen bonds to the stability (Delta DeltaGHB) was extracted on the basis of the structures of the mutant proteins. The results showed that hydrogen bonds between protein atoms and between a protein atom and a water bound with the protein molecule favorably contribute to the protein stability. The net contribution of one intramolecular hydrogen bond to protein stability (DeltaGHB) was 8.9 +/- 2.6 kJ/mol on average. However, the contribution to the protein stability of hydrogen bonds between a protein atom and a bound water molecule was smaller than that for a bond between protein atoms.     

Analysis of the dynamic equilibrium of histone oligomers in a solution. The nature of forces stabilizing the (H2A-H2B-H3-H4)2 octamer structure

Dynamic equilibrium analysis of the (H2A-H2B-H3-H4)2 histone octamer with lower oligomers was performed in 2 M NaCl. Calculated data on the relative content of histone oligomers upon changing protein concentration in solution are given. The red shift of lambda max for histone tyrosine fluorescence spectra is shown to be due to hydrogen bond formation by tyrosyl OH-groups. Analysis of free energy changes of histone oligomers upon association (delta G = -17,37 +/- 0,14 kcal/mole) as well as the effect of urea on histone octamer dissociation made it possible to conclude that virtually all tyrosyls in octamer form hydrogen bonds. Intermolecular hydrogen bonds formed by tyrosyls contribute substantially to octamer stabilization. The (H2A-H2B) dimer positive cooperativity in association with the (H3-H4)2 tetramer was found. This cooperativity is caused by interaction between association sites with a two order increase in an apparent constant of dimers with tetramer association. The histone octamer was determined to be of asymmetric structure due to unequivolency of the two binding sites for the (H2A-H2B) dimers.