Evaluation of C-H cdots,three dots,centered O hydrogen bonds in native and misfolded proteins

Non-traditional C-H cdots, three dots, centered Y hydrogen bonds, in which a carbon atom acts as the hydrogen donor and an electronegative atom Y (Y=N, O or S) acts as the acceptor, have been reported in proteins, but their importance in protein structures is not well established. Here, we present the results of three computational tests that examine the significance of C-H cdots, three dots, centered Y bonds involving the C(alpha) in proteins. First, we compared the number of C(alpha)-H cdots, three dots, centered Y bonds in native structures with two sets of compact, energy-minimized decoy structures. The decoy structures contain about as many C(alpha)-H cdots, three dots, centered Y bonds as the native structures, indicating that the constraints of chain connectivity and compactness can lead to incidental formation of C(alpha)-H cdots, three dots, centered Y bonds. Secondly, we examined whether short C(alpha)-H cdots, three dots, centered Y bonds have a tendency to be linear, as is expected for a cohesive hydrogen-bonding interaction. The native structures do show this trend, but so does one of the decoy sets, suggesting that this criterion is also not sufficient to indicate a stabilizing interaction. Finally, we examined the preference for C(alpha)-H cdots, three dots, centered Y bond donors to be near to strong hydrogen bond acceptors. In the native proteins, the alpha protons attract strong acceptors like oxygen atoms more than weak acceptors. In contrast, hydrogen bond donors in the decoy structures do not distinguish between strong and weak acceptors. Thus, any individual C(alpha)-H cdots, three dots, centered Y bond may be fortuitous and occur due to the polypeptide connectivity and compactness. Taken collectively, however, C(alpha)-H cdots, three dots, centered Y bonds provide a weakly cohesive force that stabilizes proteins.     

A proton nuclear magnetic resonance study of the conformation of bovine anaphylatoxin C5a in solution

The solution conformation of bovine anaphylatoxin C5a has been investigated by nuclear magnetic resonance (NMR) spectroscopy. The 1H-NMR spectrum is assigned in a sequential manner using a variety of two-dimensional NMR techniques. A qualitative interpretation of the short range nuclear Overhauser enhancement data involving the NH, C alpha H and C beta H protons suggests that C5a has four helices comprising residues 5-11, 15-25, 33-39 and 46-61, and is composed of a globular head (residues 5-61) and a C-terminal tail. The polypeptide fold was determined by hybrid distance geometry-dynamical simulated annealing calculations on the basis of 203 approximate interproton distance restraints, 22 distance restraints for 11 intrahelical hydrogen bonds (identified on the basis of the pattern of short range NOEs and slowly exchanging backbone amide protons) and restraints for the 3 disulfide bridges. The overall polypeptide fold is similar to that of the sequence related human recombinant anaphylatoxin C5a [(1988) Proteins 3, 139-145].     

Topological nature of the genetic code

A model for topological coding of proteins is proposed. The model is based on the capacity of hydrogen bonds (property of connectivity) to fix conformations of protein molecules. The protein chain is modeled by an n -arc graph with the following elements: vertices (alpha -carbon atoms), structural edges (peptide bonds) and connectivity edges (virtual edges connecting non-adjacent atoms). It was shown that 64 conformations of the 4-arc graph can be described in the binary system by matrices of six variables which form a supermatrix containing four blocks. On the basis of correspondences between the pairs of variables in matrices and four letters of the genetic code matrices and supermatrix are converted, respectively, into the triplets and the table of the genetic code. An algorithm admitting computer programming is proposed for coding the n -arc graph and protein chain. Connectivity operators (polar amino acids) are assigned to blocks of triplets coding for cyclic conformations (G, A-in the second position), while anti-connectivity operators (non-polar amino acids) correspond to blocks of triplets coding for open conformations (C, U-in the second position). Amino acids coded by triplets differing by the first base have different structures. The third base for C, U and G, A is degenerated. Properties of the real genetic code are in full agreement with the model. The model provides an insight into the topological nature of the genetic code and can be used for development of algorithms for the prediction of the protein structure.     

The three-dimensional structure of guanine-specific ribonuclease F1 in solution determined by NMR spectroscopy and distance geometry.

Two-dimensional 1H-NMR studies have been performed on ribonuclease F1 (RNase F1), which contains 106 amino acid residues. Sequence-specific resonance assignments were accomplished for the backbone protons of 99 amino acid residues and for most of their side-chain protons. The three-dimensional structures were constructed on the basis of 820 interproton-distance restraints derived from NOE, 64 distance restraints for 32 hydrogen bonds and 33 phi torsion-angle restraints. A total of 40 structures were obtained by distance geometry and simulated-annealing calculations. The average root-mean-square deviation (residues 1-106) between the 40 converged structures and the mean structure obtained by averaging their coordinates was 0.116 +/- 0.018 nm for the backbone atoms and 0.182 +/- 0.015 nm for all atoms including the hydrogen atoms. RNase F1 was determined to be an alpha/beta-type protein. A well-defined structure constitutes the core region, which consists of a small N-terminal beta-sheet (beta 1, beta 2) and a central five-stranded beta-sheet (beta 3-beta 7) packed on a long helix. The structure of RNase F1 has been compared with that of RNase T1, which was determined by X-ray crystallography. Both belong to the same family of microbial ribonucleases. The polypeptide backbone fold of RNase F1 is basically identical to that of RNase T1. The conformation-dependent chemical shifts of the C alpha protons are well conserved between RNase F1 and RNase T1. The residues implicated in catalysis are all located on the central beta-sheet in a geometry similar to that of RNase T1.     

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.     

From determination of amino acid residue conformation to reconstruction of the spatial structure of proteins (from the data of nuclear Overhauser enhancement spectroscopy

A new approach to the calculation of the spatial protein structure based on the joint utilization of the theoretical conformational analysis method and nuclear Overhauser enhancement (NOE) spectroscopy data is proposed and verified. The quality in determining various molecule structural parameters is estimated in terms of the expected NOE spectral parameters derived from the X-ray analysis data of the avian pancreatic polypeptide. The proposed approach is shown to correctly determine such structural parameters of protein molecules as local amino acid residue conformations, reciprocal spatial orientation of the C alpha atoms neighbouring along amino acid sequence and reapproached segments of the polypeptide chain. Spatially remote molecule fragments are mainly responsible for the error in determining structural parameters.     

New structural insights into lectin-type proteins of the immune system.

New structural data have emerged for the ligand-binding sites of C-type lectin domains and C-type lectin-like domains of receptors of the immune system. These include binding sites for oligosaccharide or polypeptide ligands, or both oligosaccharide and polypeptide ligands. The structural basis for the binding of a lectin domain of the beta-trefoil family to different sulfooligosaccharide sequences has been revealed. Lectin activity has been documented for a beta/alpha TIM barrel fold that does not have the chitinase activity of the prototype enzyme with this fold.     

An automatic search for similar spatial arrangements of alpha-helices and beta-strands in globular proteins

A fast search algorithm to reveal similar polypeptide backbone structural motifs in proteins is proposed. It is based on the vector representation of a polypeptide chain fold in which the elements of regular secondary structures are approximated by linear segments (Abagyan and Maiorov, J. Biomol. Struct. Dyn. 5, 1267-1279 (1988)). The algorithm permits insertions and deletions in the polypeptide chain fragments to be compared. The fast search algorithm implemented in FASEAR program is used for collecting beta alpha beta supersecondary structure units in a number of alpha/beta proteins of Brookhaven Data Bank. Variation of geometrical parameters specifying backbone chain fold is estimated. It appears that the conformation of the majority of the fragments, although almost all of them are right-handed, is quite different from that of standard beta alpha beta units. Apart from searching for specific type of secondary structure motif, the algorithm allows automatically to identify new recurrent folding patterns in proteins. It may be of particular interest for the development of tertiary template approach for prediction of protein three-dimensional structure as well for constructing artificial polypeptides with goal-oriented conformation.     

Structures of intact glycoproteins from vibrational circular dichroism

Vibrational circular dichroism (VCD) spectra for the glycoproteins alpha1-acid glycoprotein (AGP) and bovine submaxillary mucin (BSM), have been measured in D2O solutions and for the films prepared from aqueous (H2O) buffer solutions in the 1800 to 900 cm(-1) region. The solution VCD results revealed that AGP has beta-sheet structure, along with a significant amount of alpha-helix as evidenced from a W pattern in the amide I region. The VCD of BSM solution suggested a polyproline II type structure, characterized by the appearance of strong negative couplet in the amide I region. The film VCD results on AGP and BSM suggested that the secondary structures of polypeptide fold in the film state are similar to those in the solution. The absence of any significant film VCD in the low frequency region (1200-900 cm(-1)), suggested that the dominant linkage for carbohydrate residues is likely to be a beta linkage. VCD spectroscopy gains importance in the secondary structural analysis of polypeptide fold in glycoproteins due to the absence of interfering VCD from the carbohydrate residues in the conformationally sensitive amide I region. Also, film VCD studies permit measurements in the low wavenumber region (1200-900 cm(-1)) that reveal the dominant type of linkage for carbohydrate residues. Such clear structural information is unlike that from ECD, where ECD bands of acylated amino sugar residues interfere with those of polypeptide backbone in the conformationally sensitive far-UV region.     

Reduced representation model of protein structure prediction: statistical potential and genetic algorithms.

A reduced representation model, which has been described in previous reports, was used to predict the folded structures of proteins from their primary sequences and random starting conformations. The molecular structure of each protein has been reduced to its backbone atoms (with ideal fixed bond lengths and valence angles) and each side chain approximated by a single virtual united-atom. The coordinate variables were the backbone dihedral angles phi and psi. A statistical potential function, which included local and nonlocal interactions and was computed from known protein structures, was used in the structure minimization. A novel approach, employing the concepts of genetic algorithms, has been developed to simultaneously optimize a population of conformations. With the information of primary sequence and the radius of gyration of the crystal structure only, and starting from randomly generated initial conformations, I have been able to fold melittin, a protein of 26 residues, with high computational convergence. The computed structures have a root mean square error of 1.66 A (distance matrix error = 0.99 A) on average to the crystal structure. Similar results for avian pancreatic polypeptide inhibitor, a protein of 36 residues, are obtained. Application of the method to apamin, an 18-residue polypeptide with two disulfide bonds, shows that it folds apamin to native-like conformations with the correct disulfide bonds formed.     

The structure of cyclin H: common mode of kinase activation and specific features.

The crystal structure of human cyclin H refined at 2.6 A resolution is compared with that of cyclin A. The core of the molecule consists of two repeats containing five helices each and forming the canonical cyclin fold also observed in TFIIB. One hundred and thirty-two out of the 217 C alpha atoms from the cyclin fold can be superposed with a root-mean-square difference of 1.8 A. The structural homology is even higher for the residues at the interface with the kinase, which is of functional significance, as shown by our observation that cyclin H binds to cyclin-dependent kinase 2 (cdk2) and that cyclin A is able to activate cdk7 in the presence of MAT1. Based on this superposition, a new signature sequence for cyclins was found. The specificity of the cyclin H molecule is provided mainly by two long helices which extend the cyclin fold at its N- and C-termini and pack together against the first repeat on the side opposite to the kinase. Deletion mutants show that the terminal helices are required for a functionally active cyclin H.     

When fold is not important: a common structural framework for adenine and AMP binding in 12 unrelated protein families

ATP is a ligand common to many proteins, yet it is unclear whether common recognition patterns do exist among the many different folds that bind ATP. Previously, it was shown that cAMP-dependent protein kinase, D-Ala:D-Ala ligase and the alpha-subunit of the alpha 2 beta 2 ribonucleotide reductase do share extensive common structural elements for ATP recognition although their folds are different. Here, we have made a survey of structures that bind ATP and compared them with the key features seen in these three proteins. Our survey shows that 12 different fold types share a specific recognition pattern for the adenine moiety, and 8 of these folds have a common structural framework for recognition of the AMP moiety of the ligand. The common framework consists of a tripeptide segment plus three additional residues, which provides similar polar and hydrophobic interactions between the protein and mononucleotide. Consensus interactions are represented by four key hydrogen bonds present in each fold type. Two of these four hydrogen bonds, together with three aliphatic residues, form a specific recognition pattern for the adenine moiety in all 12 folds. These similarities point to a structural-functional requirement shared by these different mononucleotide-binding proteins that represent at this time 28% of the adenine mononucleotide complexes found in the Brookhaven Protein Data Bank.     

Solution structure of a lipid transfer protein extracted from rice seeds. Comparison with homologous proteins.

Nuclear magnetic resonance (NMR) spectroscopy was used to determine the three dimensional structure of rice nonspecific lipid transfer protein (ns-LTP), a 91 amino acid residue protein belonging to the broad family of plant ns-LTP. Sequence specific assignment was obtained for all but three HN backbone 1H resonances and for more than 95% of the 1H side-chain resonances using a combination of 1H 2D NOESY; TOCSY and COSY experiments at 293 K. The structure was calculated on the basis of four disulfide bridge restraints, 1259 distance constraints derived from 1H-1H Overhauser effects, 72 phi angle restraints and 32 hydrogen-bond restraints. The final solution structure involves four helices (H1: Cys3-Arg18, H2: Ala25-Ala37, H3: Thr41-Ala54 and H4: Ala66-Cys73) followed by a long C-terminal tail (T) with no observable regular structure. N-capping residues (Thr2, Ser24, Thr40), whose side-chain oxygen atoms are involved in hydrogen bonds with i + 3 amide proton additionally stabilize the N termini of the first three helices. The fourth helix involving Pro residues display a mixture of alpha and 3(10) conformation. The rms deviation of 14 final structures with respect to the average structure is 1.14 +/- 0.16 A for all heavy atoms (C, N, O and S) and 0.72 +/- 0.01 A for the backbone atoms. The global fold of rice ns-LTP is close to the previously published structures of wheat, barley and maize ns-LTPs exhibiting nearly identical pattern of the numerous sequence specific interactions. As reported previously for different four-helix topology proteins, hydrophobic, hydrogen bonding and electrostatic mechanisms of fold stabilization were found for the rice ns-LTP. The sequential alignment of 36 ns-LTP primary structures strongly suggests that there is a uniform pattern of specific long-range interactions (in terms of sequence), which stabilize the fold of all plant ns-LTPs.     

Deterministic features of side-chain main-chain hydrogen bonds in globular protein structures

A total of 19 835 polar residues from a data set of 250 non-homologous and highly resolved protein crystal structures were used to identify side-chain main-chain (SC-MC) hydrogen bonds. The ratio of the number of SC-MC hydrogen bonds to the total number of polar residues is close to 1:2, indicating the ubiquitous nature of such hydrogen bonds. Close to 56% of the SC-MC hydrogen bonds are local involving side-chain acceptor/donor ('i') and a main-chain donor/acceptor within the window i-5 to i+5. These short-range hydrogen bonds form well defined conformational motifs characterized by specific combinations of backbone and side-chain torsion angles. (a) The Ser/Thr residues show the greatest preference in forming intra-helical hydrogen bonds between the atoms O(gamma)(i) and O(i-4). More than half the examples of such hydrogen bonds are found at the middle of alpha-helices rather than at their ends. The most favoured motif of these examples is alpha(R)alpha(R)alpha(R)alpha(R)(g(-)). (b) These residues also show great preference to form hydrogen bonds between O(gamma)(i) and O(i-3), which are closely related to the previous type and though intra-helical, these hydrogen bonds are more often found at the C-termini of helices than at the middle. The motif represented by alpha(R)alpha(R)alpha(R)alpha(R)(g(+)) is most preferred in these cases. (c) The Ser, Thr and Glu are the most frequently found residues participating in intra-residue hydrogen bonds (between the side-chain and main-chain of the same residue) which are characterized by specific motifs of the form beta(g(+)) for Ser/Thr residues and alpha(R)(g(-)g(+)t) for Glu/Gln. (d) The side-chain acceptor atoms of Asn/Asp and Ser/Thr residues show high preference to form hydrogen bonds with acceptors two residues ahead in the chain, which are characterized by the motifs beta (tt')alphaR and beta(t)alpha(R), respectively. These hydrogen bonded segments, referred to as Asx turns, are known to provide stability to type I and type I' beta-turns. (e) Ser/Thr residues often form a combination of SC-MC hydrogen bonds, with the side-chain donor hydrogen bonded to the carbonyl oxygen of its own peptide backbone and the side-chain acceptor hydrogen bonded to an amide hydrogen three residues ahead in the sequence. Such motifs are quite often seen at the beginning of alpha-helices, which are characterized by the beta(g(+))alpha(R)alpha(R) motif. A remarkable majority of all these hydrogen bonds are buried from the protein surface, away from the surrounding solvent. This strongly indicates the possibility of side-chains playing the role of the backbone, in the protein interiors, to satisfy the potential hydrogen bonding sites and maintaining the network of hydrogen bonds which is crucial to the structure of the protein.     

Novel protein folds and their nonsequential structural analogs

Newly determined protein structures are classified to belong to a new fold, if the structures are sufficiently dissimilar from all other so far known protein structures. To analyze structural similarities of proteins, structure alignment tools are used. We demonstrate that the usage of nonsequential structure alignment tools, which neglect the polypeptide chain connectivity, can yield structure alignments with significant similarities between proteins of known three-dimensional structure and newly determined protein structures that possess a new fold. The recently introduced protein structure alignment tool, GANGSTA, is specialized to perform nonsequential alignments with proper assignment of the secondary structure types by focusing on helices and strands only. In the new version, GANGSTA+, the underlying algorithms were completely redesigned, yielding enhanced quality of structure alignments, offering alignment against a larger database of protein structures, and being more efficient. We applied DaliLite, TM-align, and GANGSTA+ on three protein crystal structures considered to be novel folds. Applying GANGSTA+ to these novel folds, we find proteins in the ASTRAL40 database, which possess significant structural similarities, albeit the alignments are nonsequential and in some cases involve secondary structure elements aligned in reverse orientation. A web server is available at http://agknapp.chemie.fu-berlin.de/gplus for pairwise alignment, visualization, and database comparison.     

The contribution of C alpha-H...O hydrogen bonds to membrane protein stability depends on the position of the amide

Structural analyses of membrane proteins reveal a large number of C(alpha)-H...O contacts between transmembrane helices, presumed to be hydrogen bonds. Recent experiments produced conflicting results for the contribution of such hydrogen bonds to membrane protein stability. An FTIR study estimated an energy of -0.88 kcal/mol for the G79-C(alpha)-H...I76-O hydrogen bond in glycophorin A, whereas a mutagenesis study showed that the A51-C(alpha)-H...T24-O(gamma) hydrogen bond does not stabilize bacteriorhodopsin. Here, we reconcile these results using molecular mechanics calculations and an implicit membrane model (IMM1). With explicit hydrogen atoms, the potential energy of the G79-C(alpha)-H...I76-O interaction in GpA ranges from -0.54 to -0.9 kcal/mol and its contribution to stability (effective energy) from -0.49 to -0.83 kcal/mol, depending on the structural model used. The average values of these quantities in GpA-like motifs are similar. In bR, the corresponding numbers for the A51-C(alpha)-H...T24-O(gamma) interaction are +0.15 and +0.32 kcal/mol. The difference results from the different arrangement of the interacting groups and specifically the position of the acceptor with respect to the C(alpha) and N atoms. This conclusion likely applies to soluble proteins as well.     

Preferential binding of an unfolded protein to DsbA.

The oxidoreductase DsbA from the periplasm of escherichia coli introduces disulfide bonds into proteins at an extremely high rate. During oxidation, a mixed disulfide is formed between DsbA and the folding protein chain, and this covalent intermediate reacts very rapidly either to form the oxidized protein or to revert back to oxidized DsbA. To investigate its properties, a stable form of the intermediate was produced by reacting the C33A variant of DsbA with a variant of RNase T1. We find that in this stable mixed disulfide the conformational stability of the substrate protein is decreased by 5 kJ/mol, whereas the conformational stability of DsbA is increased by 5 kJ/mol. This reciprocal effect suggests strongly that DsbA interacts with the unfolded substrate protein not only by the covalent disulfide bond, but also by preferential non-covalent interactions. The existence of a polypeptide binding site explains why DsbA oxidizes protein substrates much more rapidly than small thiol compounds. Such a very fast reaction is probably important for protein folding in the periplasm, because the accessibility of the thiol groups for DsbA can decrease rapidly when newly exported polypeptide chains begin to fold.     

Structure of [4Fe-4S] ferredoxin from Bacillus thermoproteolyticus refined at 2.3 A resolution. Structural comparisons of bacterial ferredoxins

The structure of a low-potential ferredoxin isolated from Bacillus thermoproteolyticus has been refined by a restrained least-squares method. The final crystallographic R factor is 0.204 for 2906 reflections with F greater than 3 sigma F in the 6.0 to 2.3 A resolution range. The model contains 81 amino acid residues, one [4Fe-4S] cluster, and 59 water molecules. The root-mean-square deviation from ideal values for bond lengths is 0.018 A, and the mean coordinate error is estimated to be 0.25 A. The present ferredoxin is similar in the topology of the polypeptide backbone to the dicluster-type ferredoxins from Peptococcus aerogenes and Azotobacter vinelandii, but has considerable insertions and deletions of the peptide segments as well as different secondary structures. Although all but the C-terminal C zeta atoms of P. aerogenes ferredoxin superpose on the C alpha atoms of A. vinelandii ferredoxin, only 60% superpose on the C alpha atoms of B. thermoproteolyticus ferredoxin, with a root-mean-square distance of 0.82 A for each pair. The conformations of the peptide segments surrounding the [4Fe-4S] clusters in these three ferredoxins are all conserved. Moreover, the schemes for the NH...S hydrogen bonds in these ferredoxins are nearly identical. The site of the aromatic ring of Tyr27 in B. thermoproteolyticus ferredoxin is close spatially to that of Tyr28 in P. aerogenes ferredoxin with reference to the cluster, but these residues do not correspond in the spatial alignment of their polypeptide backbones. We infer that in monocluster-type ferredoxins, the side-chain at the 27th residue has a crucial effect on the stability of the cluster. Of the four cysteine residues that bind to the second Fe-S cluster in the dicluster-type ferredoxins, two are conserved in the monocluster-type ferredoxins from Desulfovibrio gigas. D. desulfuricans Norway, and Clostridium thermoaceticum. The tertiary structure of B. thermoproteolyticus ferredoxin suggests that in such monocluster-type ferredoxins these two cysteine residues, which in it correspond to Ala21 and Asp53, form a disulfide bridge.