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.     

Classification of RNA structures based on hydrogen bond and base-base stacking patterns: application for NMR structures

A computational system, CSNA, for classifying RNA structures according to structural characters was developed. CSNA lists up all the hydrogen bonds and base-base stackings in the structures, and classifies the structures into sub-groups based on their patterns as the first step grouping. The frequency of each hydrogen bond or base-base stacking is calculated, the frequency score being defined as the sum of the frequency of existing hydrogen bonds or base-base stackings for each sub-group. Finally, the sub-groups are further classified into groups based on the frequency score defined in this study and the difference between the patterns. According to the frequency score, CSNA suggests a group that shares most frequently appearing hydrogen bonds and base-base stackings. CSNA was applied to the classification of the results of two individual simulated annealing calculations based on NMR information. It was found that CSNA could extract structures with lower energy without checking any energy term and could provide well converged groups as the lowest energy structures. Thus, CSNA could be a new tool for structural determination of nucleic acids.     

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.     

An orientation-dependent hydrogen bonding potential improves prediction of specificity and structure for proteins and protein-protein complexes

Hydrogen bonding is a key contributor to the specificity of intramolecular and intermolecular interactions in biological systems. Here, we develop an orientation-dependent hydrogen bonding potential based on the geometric characteristics of hydrogen bonds in high-resolution protein crystal structures, and evaluate it using four tests related to the prediction and design of protein structures and protein-protein complexes. The new potential is superior to the widely used Coulomb model of hydrogen bonding in prediction of the sequences of proteins and protein-protein interfaces from their structures, and improves discrimination of correctly docked protein-protein complexes from large sets of alternative structures.     

Hydrophobic tendencies of polar groups as a major force in molecular recognition

Proteins and nucleic acids are able to adopt their native conformation and perform their biological role only in the presence of water with which they actively interact in a mutually modifying way. Traditionally, hydrophobic effect has been considered to be the major factor stabilizing biopolymeric structures. However, solvent reorganization around polar groups is an event thermodynamically more unfavorable than solvent reorganization around nonpolar groups. Consequently, burial of polar groups with formation of complementary solute-solute hydrogen bonds out of contact with water is an energetically favorable process that also provides a major force driving macromolecular association and folding. In contrast to nonpolar groups, polar groups may form their complementary intra- or intersolute hydrogen bonds out of contact with water only provided that an appropriate solute structure has been formed with properly positioned hydrogen bond donors and acceptors. Formation of such structures is disfavored entropically and may not be possible due to steric reasons. However, the interior of a folded protein, alpha-helices and beta-sheets, double helical nucleic acid structures, and protein-ligand interfaces all provide rigid matrices where polar groups may form their complementary hydrogen bonds. For these structures, the inward drive of polar groups represents a considerable stabilizing factor.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

Mapping the lifetimes of local opening events in a native state protein.

The rate constants for the processes that lead to local opening and closing of the structures around hydrogen bonds in native proteins have been determined for most of the secondary structure hydrogen bonds in the four-helix protein acyl coenzyme A binding protein. In an analysis that combines these results with the energies of activation of the opening processes and the stability of the local structures, three groups of residues in the protein structure have been identified. In one group, the structures around the hydrogen bonds have frequent openings, every 600 to 1,500 s, and long lifetimes in the open state, around 1 s. In another group of local structures, the local opening is a very rare event that takes place only every 15 to 60 h. For these the lifetime in the open state is also around 1 s. The majority of local structures have lifetimes between 2,000 and 20,000 s and relatively short lifetimes of the open state in the range between 30 and 400 ms. Mapping of these groups of amides to the tertiary structure shows that the openings of the local structures are not cooperative at native conditions, and they rarely if ever lead to global unfolding. The results suggest a mechanism of hydrogen exchange by progressive local openings.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

Secondary structures and structural fluctuation in a dimeric protein, Streptomyces subtilisin inhibitor.

Based on the nuclear magnetic resonance assignments of a dimeric protein, Streptomyces subtilisin inhibitor (SSI), microscopic details of secondary structures in solution have been elucidated. The chemical shift index of C(alpha) signals, together with information on the hydrogen exchange rates of the backbone amide protons, were used to identify secondary structures. The locations of these secondary structures were found to be different in some critical points from those determined earlier by X-ray crystallography of the crystal. Notably, the beta3 strand is completely missing and the alpha2 helix is extended toward the C-terminus. Furthermore, hydrogen exchange experiments of individual peptide NH protons under strongly folding conditions revealed mechanisms of global and local structural fluctuation within the dimeric structure. It has been suggested that the global fluctuation of the monomeric unit occurs without affecting the accompanying monomer, in contrast to the equilibrium thermal unfolding, which is cooperative. Higher protection against hydrogen exchange for residues in part of the beta4 strand implies that this region might serve as a folding core.     

Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2

Protection factors obtained from equilibrium hydrogen exchange experiments are an important source of structural information on both native and nonnative states of proteins. We present a method for determining ensembles of protein structures by using hydrogen exchange data as restraints in molecular dynamics simulations in conjunction with an empirical force-field. The method is applied to determine the ensemble of structures representing the native state of chymotrypsin inhibitor 2 (CI2), including the rare, large fluctuations responsible for hydrogen exchange.     

The plasticity of a structural motif in RNA: Structural polymorphism of a kink turn as a function of its environment

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2′ at the −1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.     

Hydrogen bond connectivity patterns and hydrophobic interactions in crystal structures of small, acyclic peptides.

Crystal structures of all available unblocked linear peptides with two to five residues were retrieved from the Cambridge Structural Database and their intermolecular contacts and packing modes studied using molecular graphics. This survey reveals that interactions between hydrophobic portions of the molecules are critically important in determining the overall features of their crystal packing patterns. Distinct hydrophobic columns or layers are observed in almost all crystal structures. Analyses of the relationships between these interactions and crystal growth properties of small peptides are given. It is suggested that needle growth is promoted by hydrophobic packing, usually along a short crystallographic axis (4.6-6.0 angstroms). Also contributing to these morphologic characteristics are entropic factors associated with hydrophobic aggregation as well as tightly bound water molecules on hydrophobic faces. The paper also provides a comprehensive overview of hydrogen bond patterns in acyclic peptide crystals. It is demonstrated that one of their primary roles is to provide a scaffolding within which hydrophobic groups can aggregate. Even though there is a high density of hydrogen bonds in the crystals, often with complex patterns and networks, certain motifs are found to recur in a number of structures indicating specific hydrogen bond preferences. Water, for example, is an integral part of the hydrogen bond networks in these crystals, usually acting as the primary donor for main-chain carboxylate groups in peptide hydrates.     

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.     

HAAD: A Quick Algorithm for Accurate Prediction of Hydrogen Atoms in Protein Structures

Hydrogen constitutes nearly half of all atoms in proteins and their positions are essential for analyzing hydrogen-bonding interactions and refining atomic-level structures. However, most protein structures determined by experiments or computer prediction lack hydrogen coordinates. We present a new algorithm, HAAD, to predict the positions of hydrogen atoms based on the positions of heavy atoms. The algorithm is built on the basic rules of orbital hybridization followed by the optimization of steric repulsion and electrostatic interactions. We tested the algorithm using three independent data sets: ultra-high-resolution X-ray structures, structures determined by neutron diffraction, and NOE proton-proton distances. Compared with the widely used programs CHARMM and REDUCE, HAAD has a significantly higher accuracy, with the average RMSD of the predicted hydrogen atoms to the X-ray and neutron diffraction structures decreased by 26% and 11%, respectively. Furthermore, hydrogen atoms placed by HAAD have more matches with the NOE restraints and fewer clashes with heavy atoms. The average CPU cost by HAAD is 18 and 8 times lower than that of CHARMM and REDUCE, respectively. The significant advantage of HAAD in both the accuracy and the speed of the hydrogen additions should make HAAD a useful tool for the detailed study of protein structure and function. Both an executable and the source code of HAAD are freely available at http://zhang.bioinformatics.ku.edu/HAAD.     

The catalytic mechanism of dye-decolorizing peroxidase DyP may require the swinging movement of an aspartic acid residue

The dye-decolorizing peroxidase (DyP)-type peroxidase family is a unique heme peroxidase family. The primary and tertiary structures of this family are obviously different from those of other heme peroxidases. However, the details of the structure-function relationships of this family remain poorly understood. We show four high-resolution structures of DyP (EC1.11.1.19), which is representative of this family: the native DyP (1.40 ?), the D171N mutant DyP (1.42 ?), the native DyP complexed with cyanide (1.45 ?), and the D171N mutant DyP associated with cyanide (1.40 ?). These structures contain four amino acids forming the binding pocket for hydrogen peroxide, and they are remarkably conserved in this family. Moreover, these structures show that OD2 of Asp171 accepts a proton from hydrogen peroxide in compound I formation, and that OD2 can swing to the appropriate position in response to the ligand for heme iron. On the basis of these results, we propose a swing mechanism in compound I formation. When DyP reacts with hydrogen peroxide, OD2 swings towards an optimal position to accept the proton from hydrogen peroxide bound to the heme iron.     

A QM/MM analysis of the conformations of crystalline sucrose moieties

Both ab initio quantum mechanics (QM) and molecular mechanics (MM) were used to produce a hybrid energy surface for sucrose that simultaneously provides low energies for conformations that are observed in crystal structures and high energies for most unobserved structures. HF/6-31G* QM energies were calculated for an analogue based on tetrahydropyran (THP) and tetrahydrofuran (THF). Remaining contributions to the potential energy of sucrose were calculated with MM. To do this, the MM surface for the analogue was subtracted from the MM surface for the disaccharide, and the QM surface for the analogue was added. Prediction of the distribution of observable geometries was enhanced by reducing the strength of the hydrogen bonding. Reduced hydrogen-bonding strength is probably useful because many crystalline sucrose moieties do not have intramolecular hydrogen bonds between the fructose and glucose residues. Therefore, hydrogen bonding does not play a large role in determining the molecular conformation. On the hybrid energy surface that was constructed with a dielectric constant of 3.5, the average potential energy of 23 sucrose moieties from crystal structures is 1.16 kcal/mol, and the population of observed structures drops off exponentially as the energy increases.     

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.