Prediction of C alpha-H...O and C alpha-H...pi interactions in proteins using recurrent neural network

In this study, an attempt has been made to develop a method for predicting weak hydrogen bonding interactions, namely, C alpha-H...O and C alpha-H...pi interactions in proteins using artificial neural network. Both standard feed-forward neural network (FNN) and recurrent neural networks (RNN) have been trained and tested using five-fold cross-validation on a non-homologous dataset of 2298 protein chains where no pair of sequences has more than 25% sequence identity. It has been found that the prediction accuracy varies with the separation distance between donor and acceptor residues. The maximum sensitivity achieved with RNN for C alpha-H...O is 51.2% when donor and acceptor residues are four residues apart (i.e. at delta D-A = 4) and for C alpha-H...pi is 82.1% at delta D-A = 3. The performance of RNN is increased by 1-3% for both types of interactions when PSIPRED predicted protein secondary structure is used. Overall, RNN performs better than feed-forward networks at all separation distances between donor-acceptor pair for both types of interactions. Based on the observations, a web server CHpredict (available at http://www.imtech.res.in/raghava/chpredict/) has been developed for predicting donor and acceptor residues in C alpha-H...O and C alpha-H...pi interactions in proteins.     

Tyrosine hydrogen bonds make a large contribution to protein stability

The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.     

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.     

Investigations on unconventional hydrogen bonds in RNA binding proteins: the role of CH...O=C interactions

We have investigated the roles played by C-H...O=C interactions in RNA binding proteins. There was an average of 78 CH...O=C interactions per protein and also there was an average of one significant CH...O=C interactions for every 6 residues in the 59 RNA binding proteins studied. Main chain-Main chain (MM) CH...O=C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH...O=C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH...O=C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH...O=C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH...O=C interacting polar amino acid residues were solvent exposed while, majority of the CH...O=C interacting non polar residues were excluded from the solvent. Long and medium-range CH...O=C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH...O=C interacting residues had a higher conservation score. Significant percentage of CH...O=C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH...O=C interactions and hence these residues may also contribute additional stability to RNA binding proteins.     

Non-classical hydrogen bonds in interleukins: the role of C-H...O interactions

The role of classical hydrogen bonds in the structural stability of biological macro-molecules is well understood. In the present study, we explore the influence of C-H...O interactions in relation to other environmental preferences in interleukins. Main chain-main chain interactions are predominant. Pro residues might stabilize helices and strands by C-H...O H-bonds in interleukins. Majority of the C-H...O interacting residues were solvent exposed. 62% of C-H...O interactions was long-range interactions. The results presented in this study might be useful for structural stability studies in interleukins.     

Functional role of C(alpha)-H...O hydrogen bonds between transmembrane alpha-helices in photosystem I

The crystal structure at 2.5A resolution of the membrane-intrinsic, homotrimeric photosystem I (PSI) isolated from the thermophilic cyanobacterium Synechococcus elongatus shows that each monomer is composed of 12 protein subunits of which nine are embedded in the membrane and feature a total of 34 transmembrane alpha-helices (TMH). Hence, PSI provides an ideal case to study "conventional" and C(alpha)-H...O hydrogen bonds between TMH engaged in intra- and intersubunit interactions. Of the total of 75 C(alpha)-H...O hydrogen bonds between TMHs, 72 are intrasubunit and only three are intersubunit. The two largest subunits PsaA and PsaB are each folded into 11 TMHs showing 29 and 24 intrasubunit C(alpha)-H...O hydrogen bonds, respectively, that are not distributed randomly but many of them flank chlorophyll a (Chl a) co-ordinating amino acids, suggesting stabilisation of these structural segments. As major constituent of the trimerisation domain, subunit PsaL is located next to the 3-fold axis relating the three monomers of PSI. PsaL features a unique number of 19 intrasubunit C(alpha)-H...O hydrogen bonds that connect two of its three TMHs but there are no intersubunit C(alpha)-H...O hydrogen bonds between the three PsaL. Of the three intersubunit C(alpha)-H...O hydrogen bonds, two are formed between PsaA and PsaB and one between PsaB and PsaM. The large number of 75 C(alpha)-H...O hydrogen bonds contrasts the 49 conventional hydrogen bonds, indicating that the former and van der Waals contacts determine association and orientation of TMHs in PSI.     

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.     

Kinase inhibitors and the case for CH...O hydrogen bonds in protein-ligand binding

Although the hydrogen bond is known to be an important mediator of intermolecular interactions, there has yet to be an analysis of the role of CH...O hydrogen bonds in protein-ligand complexes. In this work, we present evidence for such nonstandard hydrogen bonds from a survey of aromatic ligands in 184 kinase crystal structures and 358 high-resolution structures from the Protein Data Bank. CH groups adjacent to the positively charged nitrogen of nicotinamide exhibit geometric preferences strongly suggestive of hydrogen bonding interactions, as do heterocyclic CH groups in kinase ligands, while other aromatic CH groups do not exhibit these characteristics. Ab initio calculations reveal a considerable range of CH...O hydrogen bonding potentials among different aromatic ring systems, with nicotinamide and heterocycles preferred in kinase inhibitors showing particularly favorable interactions. These results provide compelling evidence for the existence of CH...O hydrogen bonds in protein-ligand interactions, as well as information on the relative strength of various aromatic CH donors. Such knowledge will be of considerable value in protein modeling, ligand design, and structure-activity analysis.     

The occurrence of C--H...O hydrogen bonds in alpha-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.     

Is the hydrophobic effect stabilizing or destabilizing in proteins? The contribution of disulphide bonds to protein stability

It has been recently concluded that the hydrophobic effect, hitherto regarded as a major driving force in the folding of proteins, destabilizes the folded state relative to the unfolded state. We summarize the properties of the hydrophobic effect obtained from solvent transfer experiments and show that the recent conclusion is an artifact of crosslinking in the unfolded state, caused by disulphide bonds, metals or cofactors. We show that, for the proteins in the data set, crosslinks surprisingly destabilize folded structures entropically, but stabilize them enthalpically to a greater extent. We also calculate non-polar surface areas of these unfolded proteins. These surface areas are decreased by crosslinks. The unfolded state of proteins lacking constraints, such as myoglobin, is well approximated by a mixture of residues containing alpha-helical and beta-sheet dihedral angles. Surface areas of unfolded proteins cannot be obtained by summing the surface areas of individual residues, since this ignores any unavoidable side-chain-side-chain interactions.     

C-H...O hydrogen bonds in the nuclear receptor RARgamma--a potential tool for drug selectivity

Hydrogen bonds between polarized atoms play a crucial role in protein interactions and are often used in drug design, which usually neglects the potential of C-H...O hydrogen bonds. The 1.4 A resolution crystal structure of the ligand binding domain of the retinoic acid receptor RARgamma complexed with the retinoid SR11254 reveals several types of C-H...O hydrogen bonds. A striking example is the hydroxyl group of the ligand that acts as an H bond donor and acceptor, leading to a synergy between classical and C-H...O hydrogen bonds. This interaction introduces both specificity and affinity within the hydrophobic ligand pocket. The similarity of intraprotein and protein-ligand C-H...O interactions suggests that such bonds should be considered in rational drug design approaches.     

Understanding the role of hydrogen bonds in water dynamics and protein stability

The mechanisms of cold and pressure denaturation of proteins are a matter of debate, but it is commonly accepted that water plays a fundamental role in the process. It has been proposed that the denaturation process is related to an increase of hydrogen bonds among hydration water molecules. Other theories suggest that the causes of denaturation are the density fluctuations of surface water, or the destabilization of hydrophobic contacts as a consequence of water molecule inclusions inside the protein, especially at high pressures. We review some theories that have been proposed to give insight into this problem, and we describe a coarse-grained model of water that compares well with experiments for proteins’ hydration water. We introduce its extension for a homopolymer in contact with the water monolayer and study it by Monte Carlo simulations in an attempt to understand how the interplay of water cooperativity and interfacial hydrogen bonds affects protein stability.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

The contribution of electrostatics to hydrogen exchange in the unfolded protein state

Although electrostatics have long been recognized to play an important role in hydrogen exchange (HX) with solvent, the quantitative assessment of its magnitude in the unfolded state has hitherto been lacking. This limits the utility of HX as a quantitative method to study protein stability, folding, and dynamics. Using the intrinsically disordered human protein α-synuclein as a proxy for the unfolded state, we show that a hybrid mean-field approach can effectively compute the electrostatic potential at all backbone amide positions along the chain. From the electrochemical potential, a fourfold reduction in hydroxide concentration near the protein backbone is predicted for the C-terminal domain, a prognosis that is in direct agreement with experimentally derived protection factors from NMR spectroscopy. Thus, impeded HX for the C-terminal region of α-synuclein is not the result of intramolecular hydrogen bonding and/or structure formation.     

The effect of 17O on the relaxation of an amide proton within a hydrogen bond

Summary The relaxation rates of the multiple-quantum coherence for the amide hydrogen of Gly¹³ in ras p21·GDP were determined in the presence and absence of ¹⁷O labeling in the -phosphate of GDP. No significant difference could be observed between labeled and unlabeled samples, in spite of the fact that the hydrogen bond from the amide group of Gly¹³ to the -phosphate is shorter than is typical, based on its chemical shift. For macromolecules in which an oxygen atom is the acceptor of a hydrogen bond, dipolar coupling between ¹⁷O and hydrogen is unlikely to be observable, except for extremely short H-bonds.Abbreviations p21 21-kD protein product of the human H-ras gene - GMPPCP guanylyl [,-methylene]diphosphate - HPLC high-performance liquid chromatography - CSA chemical shift anisotropy - HMQC heteronuclear multiple-quantum coherence     

The contribution of polar group burial to protein stability is strongly context-dependent

We previously suggested that proteins gain more stability from the burial and hydrogen bonding of polar groups than from the burial of nonpolar groups (Pace, C. N. (2001) Biochemistry 40, 310-313). To study this further, we prepared eight Thr-to-Val mutants of RNase Sa, four in which the Thr side chain is hydrogen-bonded and four in which it is not. We measured the stability of these mutants by analyzing their thermal denaturation curves. The four hydrogen-bonded Thr side chains contribute 1.3 +/- 0.9 kcal/mol to the stability; those that are not still contribute 0.4 +/- 0.9 kcal/mol to the stability. For 40 Thr-to-Val mutants of 11 proteins, the average decrease in stability is 1.0 +/- 1.0 kcal/mol when the Thr side chain is hydrogen-bonded and 0.0 +/- 0.5 kcal/mol when it is not. This is clear evidence that hydrogen bonds contribute favorably to protein stability. In addition, we prepared four Val-to-Thr mutants of RNase Sa, measured their stability, and determined their crystal structures. In all cases, the mutants are less stable than the wild-type protein, with the decreases in stability ranging from 0.5 to 4.4 kcal/mol. For 41 Val-to-Thr mutants of 11 proteins, the average decrease in stability is 1.8 +/- 1.3 kcal/mol and is unfavorable for 40 of 41 mutants. This shows that placing an [bond]OH group at a site designed for a [bond]CH3 group is very unfavorable. So, [bond]OH groups can contribute favorably to protein stability, even if they are not hydrogen-bonded, if the site was selected for an [bond]OH group, but they will make an unfavorable contribution to stability, even if they are hydrogen-bonded, when they are placed at a site selected for a [bond]CH3 group. The contribution that polar groups make to protein stability depends strongly on their environment.     

Correlation of inhibitor effects on enzyme activity and thermal stability for the integral membrane protein fatty acid amide hydrolase

The melting curves of fatty acid amide hydrolase (FAAH) in the presence of 29 reversible inhibitors were measured using a thiol-reactive fluorophore. The thermal stability (T(m)) of the FAAH/inhibitor complex varied significantly depending on the chemical characteristics of the inhibitors, notably variations in the head group. Two separate distributions were observed when T(m) was plotted against K(i). The majority of the inhibitors showed a positive correlation between binding affinity and T(m), however inhibitors with a pyridine carboxylic acid moiety in the head group fell in a distinct and uncorrelated distribution when tail groups were varied.     

Direct detection of N-H[...]O=C hydrogen bonds in biomolecules by NMR spectroscopy

A nuclear magnetic resonance (NMR) experiment is described for the direct detection of N-H[...]O=C hydrogen bonds (H-bonds) in 15N and 13C isotope-labeled biomolecules. This quantitative 'long-range' HNCO-COSY (correlation spectroscopy) experiment detects and quantifies electron-mediated scalar couplings across the H-bond (H-bond scalar couplings), which connect the magnetically active (15)N and (13)C nuclei on both sides of the H-bond. Detectable H-bonds comprise the canonical backbone H-bonds in proteins as well as other H-bonds in proteins and nucleic acids with N-H donors and O=C (carbonylic or carboxylic) acceptors. Unlike other NMR observables, which provide only indirect evidence of the presence of H-bonds, the H-bond scalar couplings identify all partners of the H-bond, the donor, the donor proton and the acceptor, in a single experiment. The size of the scalar couplings can be related to H-bond geometries. The time required to detect the N-H[...]O=C H-bonds in small proteins (< or = approximately 10 kDa) is typically on the order of 1 d at millimolar concentrations, whereas H-bond detection for larger proteins (< or = approximately 30 kDa) may be possible within several days depending on concentration, isotope composition, magnetic field strength and molecular weight. The proteins ubiquitin (8.6 kDa), dimeric RANTES (2 x 8.5 kDa) and MAP30 (30 kDa) are used as examples to illustrate this procedure.