The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G·A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K_D = 10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G·A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G·A pairs to the stability of the K-turn.     

Observation of a unique pattern of bifurcated hydrogen bonds in the crystal structures of the N-glycoprotein linkage region models

Elucidation of the intra- and intermolecular carbohydrate-protein interactions would greatly contribute toward obtaining a better understanding of the structure-function correlations of the protein-linked glycans. The weak interactions involving C-H...O have recently been attracting immense attention in the domain of biomolecular recognition. However, there has been no report so far on the occurrence of C-H...O hydrogen bonds in the crystal structures of models and analogs of N-glycoproteins. We present herein an analysis of C-H...O interactions in the crystal structures of all N-glycoprotein linkage region models and analogs. The study reveals a cooperative network of bifurcated hydrogen bonds consisting of N-H...O and C-H...O interactions seen uniquely for the models. The cooperative network consists of two antiparallel chains of bifurcated hydrogen bonds, one involving N1-H, C2'-H and O1' of the aglycon moiety and the other involving N2-H, C1-H and O1' of the sugar. Such bifurcated hydrogen bonds between the core glycan and protein are likely to play an important role in the folding and stabilization of proteins.     

The Role of Backbone Hydrogen Bonds in the Transition State for Protein Folding of a PDZ Domain

Backbone hydrogen bonds are important for the structure and stability of proteins. However, since conventional site-directed mutagenesis cannot be applied to perturb the backbone, the contribution of these hydrogen bonds in protein folding and stability has been assessed only for a very limited set of small proteins. We have here investigated effects of five amide-to-ester mutations in the backbone of a PDZ domain, a 90-residue globular protein domain, to probe the influence of hydrogen bonds in a β-sheet for folding and stability. The amide-to-ester mutation removes NH-mediated hydrogen bonds and destabilizes hydrogen bonds formed by the carbonyl oxygen. The overall stability of the PDZ domain generally decreased for all amide-to-ester mutants due to an increase in the unfolding rate constant. For this particular region of the PDZ domain, it is therefore clear that native hydrogen bonds are formed after crossing of the rate-limiting barrier for folding. Moreover, three of the five amide-to-ester mutants displayed an increase in the folding rate constant suggesting that the hydrogen bonds are involved in non-native interactions in the transition state for folding.     

Energetics and cooperativity of tertiary hydrogen bonds in RNA structure.

Tertiary interactions that allow RNA to fold into intricate three-dimensional structures are being identified, but little is known about the thermodynamics of individual interactions. Here we quantify the tertiary structure contributions of individual hydrogen bonds in a "ribose zipper" motif of the recently crystallized Tetrahymena group I intron P4-P6 domain. The 2'-hydroxyls of P4-P6 nucleotides C109/A184 and A183/G110 participate in forming the "teeth" of the zipper. These four nucleotides were substituted in all combinations with their 2'-deoxy and (separately) 2'-methoxy analogues, and thermodynamic effects on the tertiary folding DeltaG degrees ' were assayed by the Mg2+ dependence of electrophoretic mobility in nondenaturing gels. The 2'-deoxy series showed a consistent trend with an average contribution to the tertiary folding DeltaG degrees' of -0.4 to -0.5 kcal/mol per hydrogen bond. Contributions were approximately additive, reflecting no cooperativity among the hydrogen bonds. Each "tooth" of the ribose zipper (comprising two hydrogen bonds) thus contributes about -1.0 kcal/mol to the tertiary folding DeltaG degrees'. Single 2'-methoxy substitutions destabilized folding by approximately 1 kcal/mol, but the trend reversed with multiple 2'-methoxy substitutions; the folding DeltaG degrees' for the quadruple 2'-methoxy derivative was approximately unchanged relative to wild-type. On the basis of these data and on temperature-gradient gel results, we conclude that entropically favorable hydrophobic interactions balance enthalpically unfavorable hydrogen bond deletions and steric clashes for multiple 2'-methoxy substitutions. Because many of the 2'-deoxy derivatives no longer have the characteristic hydrogen-bond patterns of the ribose zipper motif but simply have individual long-range ribose-base or ribose-ribose hydrogen bonds, we speculate that the energetic value of -0.4 to -0.5 kcal/mol per tertiary hydrogen bond may be more generally applicable to RNA folding.     

YbxF and YlxQ are bacterial homologs of L7Ae and bind K-turns but not K-loops

The archaeal protein L7Ae and eukaryotic homologs such as L30e and 15.5kD comprise the best characterized family of K-turn-binding proteins. K-turns are an RNA motif comprised of a bulge flanked by canonical and noncanonical helices. They are widespread in cellular RNAs, including bacterial gene-regulatory RNAs such as the c-di-GMP-II, lysine, and SAM-I riboswitches, and the T-box. The existence in bacteria of K-turn-binding proteins of the L7Ae family has not been proven, although two hypothetical proteins, YbxF and YlxQ, have been proposed to be L7Ae homologs based on sequence conservation. Using purified, recombinant proteins, we show that Bacillus subtilis YbxF and YlxQ bind K-turns (K(d) ~270 nM and ~2300 nM, respectively). Crystallographic structure determination demonstrates that both YbxF and YlxQ adopt the same overall fold as L7Ae. Unlike the latter, neither bacterial protein recognizes K-loops, a structural motif that lacks the canonical helix of the K-turn. This property is shared between the bacterial and eukaryal family members. Comparison of our structure of YbxF in complex with the K-turn of the SAM-I riboswitch and previously determined structures of archaeal and eukaryal homologs bound to RNA indicates that L7Ae approaches the K-turn at a unique angle, which results in a considerably larger RNA-protein interface dominated by interactions with the noncanonical helix of the K-turn. Thus, the inability of the bacterial and eukaryal L7Ae homologs to bind K-loops probably results from their reliance on interactions with the canonical helix. The biological functions of YbxF and YlxQ remain to be determined.     

Elbow Flexibility of the kt38 RNA Kink-Turn Motif Investigated by Free-Energy Molecular Dynamics Simulations

Kink-turns (K-turns) are common structural motifs that can introduce sharp kinks into double-stranded RNA, and have been proposed to mediate large-scale motions in the ribosome. K-turns consist of a bulge loop region flanked by trans sugar-Hoogsteen G:A pairs, and the sharp kink conformation is stabilized by A-minor interactions (adenine contacting a G:C basepair in the minor groove). Umbrella-sampling molecular dynamics simulations were used to disrupt an A-minor interaction in the ribosomal kt38 turn and to calculate the associated free-energy change. Coupling of umbrella sampling with replica exchanges between neighboring umbrella-sampling intervals could further improve the convergence of the free-energy calculations. The simulations revealed a coupled A-minor disruption and global opening of the K-turn motif, and allowed us to characterize several intermediate A-minor conformations. The calculated free-energy profile indicated a meta-stable, semi-open structure of slightly higher free energy (∼1 kcal mol⁻¹), separated by a small free-energy barrier (∼1.5 kcal mol⁻¹) from the closed (highly kinked) form. Both K-turn states are stabilized by distinct variants of the A-minor interaction. Further opening of the K-turn structure required significantly larger free-energy changes. The semi-open form had a reduced kink angle compatible with experimental data on K-turn solution structures, and opening was coupled to a continuous global unwinding of the K-turn motif. The range of free-energy changes associated with kt38 opening and unwinding are compatible with the idea that K-turns may facilitate biologically relevant motions during large-scale ribosome dynamics.     

Plasticity of the RNA kink turn structural motif

The kink turn (K-turn) is an RNA structural motif found in many biologically significant RNAs. While most examples of the K-turn have a similar fold, the crystal structure of the Azoarcus group I intron revealed a novel RNA conformation, a reverse kink turn bent in the direction opposite that of a consensus K-turn. The reverse K-turn is bent toward the major grooves rather than the minor grooves of the flanking helices, yet the sequence differs from the K-turn consensus by only a single nucleotide. Here we demonstrate that the reverse bend direction is not solely defined by internal sequence elements, but is instead affected by structural elements external to the K-turn. It bends toward the major groove under the direction of a tetraloop–tetraloop receptor. The ability of one sequence to form two distinct structures demonstrates the inherent plasticity of the K-turn sequence. Such plasticity suggests that the K-turn is not a primary element in RNA folding, but instead is shaped by other structural elements within the RNA or ribonucleoprotein assembly.     

Folding of the SAM aptamer is determined by the formation of a K-turn-dependent pseudoknot

The S-adenosylmethionine (SAM) riboswitch is one of the most recurrent riboswitches found in bacteria and has three known different natural aptamers. The Bacillus subtilis yitJ SAM riboswitch aptamer is organized around a four-way junction which is characterized by the presence of a pseudoknot and a K-turn motif. By replacing the adenine involved in a Watson-Crick base pair at position 138 in the core region of the aptamer with the fluorescent analogue 2-aminopurine (2AP), we show that the ligand-induced reorganization of the aptamer strongly attenuates 2AP fluorescence. The fluorescence quenching process is specific to SAM on the basis of the observation that the structural analogue S-adenosylhomocysteine does not promote a similar effect. We find that the pseudoknot is important for the reorganization of the core domain and that the K-turn motif also has a marked influence on the core domain reorganization, most probably through its important role in pseudoknot formation. Finally, we show that SAM riboswitch ligand binding is facilitated by the L7Ae K-turn binding protein, which suggests that K-turn motifs may be protein anchor sites used by riboswitches to promote RNA folding.     

Properties of antigenomic hepatitis delta virus ribozyme cis- and trans- analogs

A series of permuted variants of antigenomic HDV ribozyme and trans-acting variants were constructed. The catalytic activity study of the ribozymes has shown that all the variants were capable of self-cleaving with equally biphasic kinetics. Ribonuclease and Fe(II)-EDTA cleavage have provided evidence that all designed ribozymes fold according to the pseudoknot model and the conformations of the initial and cleaved ribozyme are different. A scheme of HDV ribozyme self-cleavage reaction was suggested. The role of hydrogen bonds in the reaction was evaluated by substitution of ribose in the ribozyme for deoxyribose. It was found that the 2'-OH group of U23 and C27 is critical for the reaction to occur; the 2'-OH group of U32 and U39 is important, while 2'-OH groups of other nucleotides of loop 3, stem 4 and stem 1 are unimportant for the cleavage activity.     

Refined structure of alpha-lytic protease at 1.7 A resolution. Analysis of hydrogen bonding and solvent structure

The structure of alpha-lytic protease, a serine protease produced by the bacterium Lysobacter enzymogenes, has been refined at 1.7 A resolution. The conventional R-factor is 0.131 for the 14,996 reflections between 8 and 1.7 A resolution with I greater than or equal to 2 sigma (I). The model consists of 1391 protein atoms, two sulfate ions and 156 water molecules. The overall root-meansquare error is estimated to be about 0.14 A. The refined structure was compared with homologous enzymes alpha-chymotrypsin and Streptomyces griseus protease A and B. A new sequence numbering was derived based on the alignment of these structures. The comparison showed that the greatest structural homology is around the active site residues Asp102, His57 and Ser195, and that basic folding pathways are maintained despite chemical changes in the hydrophobic cores. The hydrogen bonds in the structure were tabulated and the distances and angles of interaction are similar to those found in small molecules. The analysis also revealed the presence of close intraresidue interactions. There are only a few direct intermolecular hydrogen bonds. Most intermolecular interactions involve bridging solvent molecules. The structural importance of hydrogen bonds involving the side-chain of Asx residues is discussed. All the negatively charged groups have a counterion nearby, while the excess positively charged groups are exposed to the solvent. One of the sulfate ions is located near the active site, whereas the other is close to the N terminus. Of the 156 water molecules, only seven are not involved in a hydrogen bond. Six of these have polar groups nearby, while the remaining one is in very weak density. There are nine internal water molecules, consisting of two monomers, two dimers and one trimer. No significant second shell of solvent is observed.     

Structural bases of lectin-carbohydrate affinities: comparison with protein-folding energetics.

We have made a comparative structure based analysis of the thermodynamics of lectin-carbohydrate (L-C) binding and protein folding. Examination of the total change in accessible surface area in those processes revealed a much larger decrease in free energy per unit of area buried in the case of L-C associations. According to our analysis, this larger stabilization of L-C interactions arises from a more favorable enthalpy of burying a unit of polar surface area, and from higher proportions of polar areas. Hydrogen bonds present at 14 L-C interfaces were identified, and their overall characteristics were compared to those reported before for hydrogen bonds in protein structures. Three major factors might explain why polar-polar interactions are stronger in L-C binding than in protein folding: (1) higher surface density of hydrogen bonds; (2) better hydrogen-bonding geometry; (3) larger proportion of hydrogen bonds involving charged groups. Theoretically, the binding entropy can be partitioned into three main contributions: entropy changes due to surface desolvation, entropy losses arising from freezing rotatable bonds, and entropic effects that result from restricting translation and overall rotation motions. These contributions were estimated from structural information and added up to give calculated binding entropies. Good correlation between experimental and calculated values was observed when solvation effects were treated according to a parametrization developed by other authors from protein folding studies. Finally, our structural parametrization gave calculated free energies that deviate from experimental values by 1.1 kcal/mol on the average; this amounts to an uncertainty of one order of magnitude in the binding constant.     

The 2'-OH group of the peptidyl-tRNA stabilizes an active conformation of the ribosomal PTC

The ribosome accelerates the rate of peptidyl transfer by >10(6)-fold relative to the background rate. A widely accepted model for this rate enhancement invokes entropic effects whereby the ribosome and the 2'-OH of the peptidyl-tRNA substrate precisely position the reactive moieties through an extensive network of hydrogen bonds that allows proton movement through them. Some studies, however, have called this model into question because they find the 2'-OH of the peptidyl-tRNA to be dispensable for catalysis. Here, we use an in vitro reconstituted translation system to resolve these discrepancies. We find that catalysis is at least 100-fold slower with the dA76-substituted peptidyl-tRNA substrate and that the peptidyl transferase centre undergoes a slow inactivation when the peptidyl-tRNA lacks the 2'-OH group. Additionally, the 2'-OH group was found to be critical for EFTu binding and peptide release. These findings reconcile the conflict in the literature, and support a model where interactions between active site residues and the 2'-OH of A76 of the peptidyl-tRNA are pivotal in orienting substrates in this active site for optimal catalysis.     

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.     

Mechanism of protein folding II. Lysozyme and phospholipase

Refolding of hen egg-white lysozyme assuming the formation of secondary structures (-helices and β-sheets) is carried out by the method presented in the previous paper (N. Saitô et al., Proteins; Struct. Funct. Genet. 3 (1988) 199–208). To do this, the hydrophobic interactions between the hydrophobic residues which are located at the key positions for folding and can be identified without te knowledge of the native structure, and the nonbonded interactions between every pair of atoms (except hydrogen) or groups are introduced successively from short-to medium-distance pairs. The search for the energy minimum by these interactions can afford a conformation of especially the mutual arrangements between neighboring secondary structures. When these local structures are accomplished, some of the long-distance amino-acid pairs come close together and then the possible interactions (hydrophobic, nonbonded) are introduced. The three-dimensional structure of lysozyme thus obtained is shown to have locally correct arrangements of the secondary structures, but mutual relations between long-distance parts of the chain are not similar to the native structure. The introduction of disulfide bonds between appropriate cysteine residues is necessary to reach the native structure. The choice of cysteine pairs for disulfide bonding is made by the criterion given in the paper to follow (K. Watanabe, A. Nakamura, Y. Fukuda and N. Saitô, Biophys. Chem. 40 (1991) 293). The same treatment is applied to bovine pancreatic phospholipase with 7 disulfide bonds. The formation of the antiparallel β-structures from neighboring β-strands and the problem of the folding order are also discussed.     

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH₂− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.     

Molecular and crystal structures of N-arylglycopyranosylamines formed by reaction between sulfanilamide and D-ribose, D-arabinose and D-mannose

The X-ray crystal structures of three monosaccharide derivatives prepared by the reaction of sulfanilamide with D-ribose, D-arabinose, and D-mannose have been determined. The derivatives are N-(p-sulfamoylphenyl)-alpha-D-ribopyranosylamine (1), N-(p-sulfamoylphenyl)-alpha-D-arabinopyranosylamine (2), and N-(p-sulfamoylphenyl)-beta-D-mannopyranosylamine monohydrate (3). The monosaccharide ring of 1 and 2 has the 1C4 conformation, stabilized in 1 by an intramolecular hydrogen bond from 0-2 to 0-4. Compound 3 has the 4C1 conformation at the monosaccharide ring and the gt conformation at the C-6-O-6 side chain. Occupancy of the water molecule in the crystal of 3 actually examined was 22%. The degree of interaction between sulfamoyl groups and monosaccharide moieties varies from structure to structure. The packing arrangement of 2 involves hydrogen bonding between sulfamoyl groups and monosaccharide hydroxyl groups, but interactions of this type are fewer in 1, and in 3 the hydrogen bonds are either strictly between monosaccharide hydroxyl groups or strictly between sulfamoyl groups. Pairs of hydrogen bonds (two-point contacts) link neighboring molecules in all three structures, between screw-axially related molecules in 1 and 2 and between translationally related molecules in 3. The contact in 3 defined by the O-3-H...O-5 and O-6-H...O-4 hydrogen bonds is found in several other N-aryl-beta-D-mannopyranosylamine crystal structures and is apparently an especially favorable mode of intermolecular interaction in these compounds.     

Site-specific modification of functional groups in genomic hepatitis delta virus (HDV) ribozyme.

Human hepatitis delta (HDV) ribozyme is one of small ribozymes, such as hammerhead and hairpin ribozymes, etc. Its secondary structure shows pseudoknot structure composed of four stems (I to IV) and three single-stranded regions (SSrA, -B and -C). The 3D structure of 3'-cleaved product of genomic HDV ribozyme provided extensive information about tertiary hydrogen bonding interactions between nucleotide bases, phosphate oxygens and 2'OHs including new stem structure P1.1. To analyze the role of these hydrogen bond networks in the catalytic reaction, site-specific atomic-level modifications (such as deoxynucleotides, deoxyribosyl-2-aminopurine, deoxyribosylpurine, 7-deaza-ribonucleotide and inosine) were incorporated in the smallest trans-acting HDV ribozyme (47-mer). Kinetic analysis of these ribozyme variants demonstrated the importance of the two W-C base pairs of P1.1 for cleavage; in addition, the results suggest that all hydrogen bond interactions detected in the crystal structure involving 2'-OH and N7 atoms are present in the active ribozyme structure. In most of the variants, the relative reduction in kobs caused by substitution of the 2'-OH group correlated with the number of hydrogen bonds affected by the substitution. However G74 and C75 may have more than one hydrogen bond involving the 2'-OH in both the trans- and cis-acting HDV ribozyme. Moreover, in variants in which N7 was deleted, kobs was reduced 5- to 15-fold, it may suggest that N7 assists in coordinating Mg2+ ions or water molecules which bind with weak affinity in the active structure.     

Interaction geometry involving planar groups in protein-protein interfaces

The geometry of interactions of planar residues is nonrandom in protein tertiary structures and gives rise to conventional, as well as nonconventional (X--H...pi, X--H...O, where X = C, N, or O) hydrogen bonds. Whether a similar geometry is maintained when the interaction is across the protein-protein interface is addressed here. The relative geometries of interactions involving planar residues, and the percentage of contacts giving rise to different types of hydrogen bonds are quite similar in protein structures and the biological interfaces formed by protein chains in homodimers and protein-protein heterocomplexes--thus pointing to the similarity of chemical interactions that occurs during protein folding and binding. However, the percentage is considerably smaller in the nonspecific and nonphysiological interfaces that are formed in crystal lattices of monomeric proteins. The C--H...O interaction linking the aromatic and the peptide groups is quite common in protein structures as well as the three types of interfaces. However, as the interfaces formed by crystal contacts are depleted in aromatic residues, the weaker hydrogen bond interactions would contribute less toward their stability.