Quantitative comparison of the hydrogen bond network of A-state and native ubiquitin by hydrogen bond scalar couplings

The backbone hydrogen bond (H-bond) network of the partially folded A-state of ubiquitin (60% methanol, 40% water, pH 2) has been characterized quantitatively by (h3)J(NC)(') H-bond scalar couplings between the (15)N nuclei of amino acid H-bond donors and the (13)C carbonyl nuclei of the acceptors. Results on (h3)J(NC)(') couplings and the amide proton ((1)H(N)) chemical shifts for the A-state are compared quantitatively to the native state. The (h3)J(NC)(') correlations of the A-state show intact, nativelike H-bonds of the first beta-hairpin beta1/beta2 and the alpha-helix, albeit at lower strength, whereas the H-bonds in the C-terminal part change from a pure beta-structure to an all alpha-helical H(N)(i)-->O(i-4) connectivity pattern. A residue-specific analysis reveals that the conformations within the conserved secondary structure segments are much more homogeneous in the A-state than in the native state. Thus, the strong asymmetry of (h3)J(NC)(') couplings and (1)H(N) chemical shifts between the interior and exterior sides of the native state alpha-helix vanishes in the A-state. This indicates that the bend of this helix around the native state hydrophobic core is released in the homogeneous solvent environment of the A-state. Similarly, an irregularity in the behavior of H-bond I3-->L15 in hairpin beta1/beta2, which results from strong contacts to strand beta5 in the native state, is absent in the A-state. These findings rationalize the behavior of the (1)H(N) chemical shifts in both states and indicate that the A-state is in many aspects similar to the onset of thermal denaturation of the native state.     

Solvent-induced differentiation of protein backbone hydrogen bonds in calmodulin

In apo and holoCaM, almost half of the hydrogen bonds (H-bonds) at the protein backbone expected from the corresponding NMR or X-ray structures were not detected by h3JNC' couplings. The paucity of the h3JNC' couplings was considered in terms of dynamic features of these structures. We examined a set of seven proteins and found that protein-backbone H-bonds form two groups according to the h3JNC' couplings measured in solution. H-bonds that have h3JNC' couplings above the threshold of 0.2 Hz show distance/angle correlation among the H-bond geometrical parameters, and appear to be supported by the backbone dynamics in solution. The other H-bonds have no such correlation; they populate the water-exposed and flexible regions of proteins, including many of the CaM helices. The observed differentiation in a dynamical behavior of backbone H-bonds in apo and holoCaM appears to be related to protein functions.     

Quantification of H/D Isotope Effects on Protein Hydrogen-bonds by h3JNC′ and 1JNC′ Couplings and Peptide Group 15N and 13C′ Chemical Shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.     

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.     

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

A nuclear magnetic resonance (NMR) experiment is described for the direct detection of N-H[...]N hydrogen bonds (H-bonds) in 15N isotope-labeled biomolecules. This quantitative HNN-COSY (correlation spectroscopy) experiment detects and quantifies electron-mediated scalar couplings across the H-bond (H-bond scalar couplings), which connect magnetically active (15)N nuclei of the H-bond donor and acceptor. Detectable H-bonds comprise the imino H-bonds in canonical Watson-Crick base pairs, many H-bonds in unusual nucleic acid base pairs and H-bonds between protein backbone or side-chain N-H donor and N acceptor moieties. 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 and as a time average to H-bond dynamics. The time required to detect the H-bonds is typically less than 1 d at millimolar concentrations for samples of molecular weight < or = approximately 25 kDa. A C15N/13C-labeled potato spindle tuber viroid T1 RNA domain is used as an example to illustrate this procedure.     

Hydrogen bonds in rubredoxins from mesophilic and hyperthermophilic organisms

The extent and strength of the hydrogen bond networks in rubredoxins from the hyperthermophile Pyrococcus furiosus (PfRd), and its mesophilic analogue Clostridium pasteurianum (CpRd), are examined and compared using NMR spectroscopy. NMR parameters examined in this study include through-hydrogen bond (h3)J(NC)(') scalar couplings and (1)H, (13)C, and (15)N chemical shifts, as well as covalent (1)J(NH) and (1)J(NC)(') scalar couplings. These parameters have allowed the characterization in solution of 12 hydrogen bonds in each protein. Despite a 83% sequence homology and a low RMSD for the backbone heavy atoms (0.648 A) in the crystalline state, subtle, but definite, changes have been identified in the detailed hydrogen-bonding patterns. CpRd shows an increased number of hydrogen bonds in the triple-stranded beta-sheet and an additional hydrogen bond in the multiple-turn segment including residues 14-32. On the other hand, PfRd exhibits an overall strengthening of N-H...O=C hydrogen bonds in the loops involved at the metal binding site as well as evidence for an additional NH...S(Cys) hydrogen bond involving the alanine residue 44. These data, as well as temperature dependence of the NMR parameters, suggest that the particular NMR hydrogen bond pattern found in the hyperthermophile rubredoxin leads to an increased stabilization at the metal binding pocket. It seems to result from a subtle redistribution of hydrogen-bonding interactions between the triple-stranded beta-sheet and the actual metal binding site.     

D/H amide kinetic isotope effects reveal when hydrogen bonds form during protein folding

We have exploited a procedure to identify when hydrogen bonds (H-bonds) form under two-state folding conditions using equilibrium and kinetic deuterium/hydrogen amide isotope effects. Deuteration decreases the stability of equine cytochrome c and the dimeric and crosslinked versions of the GCN4-p1 coiled coil by approximately 0. 5 kcal mol-1. For all three systems, the decrease in equilibrium stability is reflected by a decrease in refolding rates and a near equivalent increase in unfolding rates. This apportionment indicates that approximately 50% of the native H-bonds are formed in the transition state of these helical proteins. In contrast, an alpha/beta protein, mammalian ubiquitin, exhibits a small isotope effect only on unfolding rates, suggesting its folding pathway may be different. These four proteins recapitulate the general trend that approximately 50% of the surface buried in the native state is buried in the transition state, leading to the hypothesis that H-bond formation in the transition state is cooperative, with alpha-helical proteins forming a number of H-bonds proportional to the amount of surface buried in the transition state.     

Correlation of the sweetness of variants of the protein brazzein with patterns of hydrogen bonds detected by NMR spectroscopy

In sequence-function investigations, approaches are needed for rapidly screening protein variants for possible changes in conformation. Recent NMR methods permit direct detection of hydrogen bonds through measurements of scalar couplings that traverse hydrogen bonds (trans-hydrogen bond couplings). We have applied this approach to screen a series of five single site mutants of the sweet protein brazzein with altered sweetness for possible changes in backbone hydrogen bonding with respect to wild-type. Long range, three-dimensional data correlating connectivities among backbone 1HN, 15N, and 13C' atoms were collected from the six brazzein proteins labeled uniformly with carbon-13 and nitrogen-15. In wild-type brazzein, this approach identified 17 backbone hydrogen bonds. In the mutants, altered magnitudes of the couplings identified hydrogen bonds that were strengthened or weakened; missing couplings identified hydrogen bonds that were broken, and new couplings indicated the presence of new hydrogen bonds. Within the series of brazzein mutants investigated, a pattern was observed between sweetness and the integrity of particular hydrogen bonds. All three "sweet" variants exhibited the same pattern of hydrogen bonds, whereas all three "non-sweet" variants lacked one hydrogen bond at the middle of the alpha-helix, where it is kinked, and one hydrogen bond in the middle of beta-strands II and III, where they are twisted. Two of the non-sweet variants lack the hydrogen bond connecting the N and C termini. These variants showed greater mobility in the N- and C-terminal regions than wild-type brazzein.     

Ligand-induced strain in hydrogen bonds of the c-Src SH3 domain detected by NMR

Changes in the molecular conformation of proteins can result from a variety of perturbations, and can play crucial roles in the regulation of biological activity. A new solution NMR method has been applied to monitor ligand-induced changes in hydrogen bond geometry in the chicken c-Src SH3 domain. The structural response of this domain to ligand binding has been investigated by measuring trans-hydrogen bond (15)N-(13)C' scalar couplings in the free state and when bound to the high affinity class I ligand RLP2, containing residues RALPPLPRY. A comparison between hydrogen bonds in high resolution X-ray structures of this domain and those observed via (h3)J(NC') couplings in solution shows remarkable agreement. Two backbone-to-side-chain hydrogen bonds are observed in solution, and each appears to play a role in stabilization of loop structure. Reproducible ligand-induced changes in trans-hydrogen bond scalar couplings are observed across the domain that translate into changes in hydrogen bond length ranging between 0.02 to 0.12 A. The observed changes can be rationalized by an induced fit mechanism in which hydrogen bonds across the protein participate in a compensatory response to forces imparted at the protein-ligand interface. Upon ligand binding, mutual intercalation of the two Leu-Pro segments of the ligand between three aromatic side-chains protruding from the SH3 surface wedges apart secondary structural elements within the SH3 domain. This disruption is transmitted in a domino-like effect across the domain through networks of hydrogen bonded peptide planes. The unprecedented resolution obtained demonstrates the ability to characterize subtle structural rearrangements within a protein upon perturbation, and represents a new step in the endeavor to understand how hydrogen bonds contribute to the stabilization and function of biological macromolecules.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Solution structure of microcin J25, the single macrocyclic antimicrobial peptide from Escherichia coli.

The three-dimensional solution structure of microcin J25, the single cyclic representative of the microcin antimicrobial peptide class produced by enteric bacteria, was determined using two-dimensional 1H NMR spectroscopy and molecular modeling. This hydrophobic 21-residue peptide exhibits potent activity directed to Gram-negative bacteria. Its primary structure, cyclo(-V1GIGTPISFY10GGGAGHVPEY20F-), has been determined previously [Blond, A., Péduzzi, J., Goulard, C., Chiuchiolo, M. J., Barthélémy, M., Prigent, Y., Salomón, R.A., Farías, R.N., Moreno, F. & Rebuffat, S. (1999) Eur. J. Biochem., 259, 747-755]. Conformational parameters (3JNHCalphaH coupling constants, quantitative nuclear Overhauser enhancement data, chemical shift deviations, temperature coefficients of amide protons, NH-ND exchange rates) were obtained in methanol solution. Structural restraints consisting of 190 interproton distances inferred from NOE data, 11 phi backbone dihedral angle and 9 chi1 angle restraints derived from the coupling constants and three hydrogen bonds in agreement with the amide exchange rates were used as input for simulated annealing calculations and energy minimization in the program XPLOR. Microcin J25 adopts a well-defined compact structure consisting of a distorted antiparallel beta sheet, which is twisted and folded back on itself, thus resulting in three loops. Residues 7-10 and 17-20 form the more regular part of the beta sheet. The region encompassing residues Gly11-His16 consists of a distorted beta hairpin, which divides into two small loops and is stabilized by an inverse gamma turn and a type I' beta turn. The reversal of the chain leading to the Phe21-Pro6 loop results from a mixed beta/gamma turn. A cavity, in which the hydrophilic Ser8 side-chain is confined, is delimited by two crab pincer-like regions that comprise residues 6-8 and 18-1.     

Contribution of hydrogen bonding to protein stability estimated from isotope effects

An unresolved issue in structural biology concerns the relative contribution of H bonds to protein stability. We use the small molecules 4-acetamidobenzoic acid and N-acetylanthranilic acid as model compounds to relate the energetic contribution from hydrogen bonds (H bonds) to the deuterium/hydrogen amide isotope effect. N-Acetylanthranilic acid models carbonyl-amide H bonds formed during protein folding; 4-acetamidobenzoic acid models the unfolded state in which the amide H bonds to water. NMR is used to measure shifts in the pK(a) of the ionizable carboxyl group when the amides of the compounds are either protonated or deuterated. From the pK(a) shift, we obtain a quantitative scale factor: SF = partial partial differential(DeltaG(HB))/partial partial differential(RT ln Phi), where DeltaG(HB) is the change in free energy of an H bond upon isotope substitution and Phi is the fractionation factor. Isotope effect data also are reported for a small globular protein, lambda repressor, using the "C(m) experiment". The protein's isotope effect, which reports on the shape of the energy well, is converted to H-bonding free energy by applying the scale factor. We estimate that amide-related H bonds (amide-carbonyl and amide-water) contribute favorably to protein stability by approximately 30-50 kcal/mol in lambda repressor, GCN4 coiled coil, and cytochrome c but unfavorably by approximately 6 kcal/mol in ubiquitin. The results indicate that H-bond strength varies from one protein to another and presumably at different sites within the same protein.     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.     

Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

Hydrogen exchange of primary amide protons in basic pancreatic trypsin inhibitor: evidence for NH2 group rotation in buried asparagine side chains

Hydrogen-deuterium exchange is measured for the buried primary amide groups of Asn-43 and Asn-44 in bovine pancreatic trypsin inhibitor. Amide protons trans and cis to the amide carbonyl oxygen (HE and HZ, respectively) exchange at indistinguishable rates. Uncorrelated exchange of HE and HZ is established for both residues by following the nuclear Overhauser enhancement from HE to HZ during the deuterium exchange. The exchange of Asn-43 and Asn-44 side-chain protons differs qualitatively from exchange of primary amide groups in fully solvated model compounds, for which HE generally exchanges faster than HZ. The equal rates for the buried primary amide HE and HZ in BPTI are not a consequence of coupled exchange. The data indicate rapid rotation around the CO-NH2 bond for both Asn-43 and Asn-44 and suggest considerable lability of intramolecular hydrogen bonds. The side chain of Asn-43 has all of its polar atoms integrated into the very stable hydrogen-bonded structure of the protein. Asn-44 is hydrogen-bonded to side chains and to a buried water molecule. Solvent isotope exchange is several orders of magnitude more restricted by protein secondary and tertiary structure than the CO-NH2 rotation, indicating that N delta H2 groups flip many times before hydrogen isotope exchange occurs.     

Porcine pancreatic phospholipase A2: sequence-specific 1H and 15N NMR assignments and secondary structure

The solution structure of porcine pancreatic phospholipase A2 (124 residues, 14 kDa) has been studied by two-dimensional homonuclear 1H and two- and three-dimensional heteronuclear 15N-1H nuclear magnetic resonance spectroscopy. Backbone assignments were made for 117 of the 124 amino acids. Short-range nuclear Overhauser effect (NOE) data show three alpha-helices from residues 1-13, 40-58, and 90-109, an antiparallel beta-sheet for residues 74-85, and a small antiparallel beta-sheet between residues 25-26 and 115-116. A 15N-1H heteronuclear multiple-quantum correlation experiment was used to monitor amide proton exchange over a period of 22 h. In total, 61 amide protons showed slow or intermediate exchange, 46 of which are located in the three large helices. Helix 90-109 was found to be considerably more stable than the other helices. For the beta-sheets, four hydrogen bonds could be identified. The secondary structure of porcine PLA in solution, as deduced from NMR, is basically the same as the structure of porcine PLA in the crystalline state. Differences were found in the following regions, however. Residues 1-6 in the first alpha-helix are less structured in solution than in the crystal structure. Whereas in the crystal structure residues 24-29 are involved both in a beta-sheet with residues 115-117 and in a hairpin turn, the expected hydrogen bonds between residues 24-117 and 25-29 do not show slow exchange behavior. This and the absence of several expected NOEs imply that this region has a less well defined structure in solution. Finally, the hydrogen bond between residues 78-81, which is part of a beta-sheet, does not show slow exchange behavior.     

The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

Observation of the closing of individual hydrogen bonds during TFE-induced helix formation in a peptide

Helix formation of an S-peptide analog, comprising the first 20 residues of Ribonuclease A and two additional N-terminal residues, was studied by measuring hydrogen bond (H-bond) (h3)J(NC') scalar couplings as a function of 2,2,2-trifluoroethanol (TFE) concentration. The (h3)J(NC') couplings give direct evidence for the closing of individual backbone N-H***O = C H-bonds during the TFE-induced formation of secondary structure. Whereas no (h3)J(NC') correlations could be detected without TFE, alpha-helical (i,i +4) H-bond correlations were observed for the amides of residues A5 to M15 in the presence of TFE. The analysis of individual coupling constants indicates that alpha-helix formation starts at the center of the S-peptide around residue E11 and proceeds gradually from there to both peptide ends as the TFE concentration is increased. At 60% to 90% TFE, well-formed alpha-helical H-bonds were observed for the amides hydrogens of residues K9 to Q13, whereas H-bonds of residues T5 to A8, H14, and M15 are affected by fraying. No intramolecular backbone H-bonds are present at and beyond the putative helix stop signal D16. As the (h3)J(NC') constants represent ensemble averages and the dependence of (h3)J(NC') on H-bond lengths is very steep, the size of the individual (h3)J(NC') coupling constants can be used as a measure for the population of a closed H-bond. These individual populations are in agreement with results derived from the Lifson-Roig theory for coil-to-helix transitions. The present work shows that the closing of individual H-bonds during TFE-induced helix formation can be monitored by changes in the size of H-bond scalar couplings.     

Effects of amino acid substitutions at beta 131 on the structure and properties of hemoglobin: evidence for communication between alpha 1 beta 1- and alpha 1 beta 2-subunit interfaces

Substitutions of Asn, Glu, and Leu for Gln at the beta131 position of the hemoglobin molecule result in recombinant hemoglobins (rHbs) with moderately lowered oxygen affinity and high cooperativity compared to human normal adult hemoglobin (Hb A). The mutation site affects the hydrogen bonds present at the alpha(1)beta(1)-subunit interface between alpha103His and beta131Gln as well as that between alpha122His and beta35Tyr. NMR spectroscopy shows that the hydrogen bonds are indeed perturbed; in the case of rHb (beta131Gln --> Asn) and rHb (beta131Gln --> Leu), the perturbations are propagated to the other alpha(1)beta(1)-interface H-bond involving alpha122His and beta35Tyr. Proton exchange measurements also detect faster exchange rates for both alpha(1)beta(1)-interface histidine side chains of the mutant rHbs in 0.1 M sodium phosphate buffer at pH 7.0 than for those of Hb A under the same conditions. In addition, the same measurements in 0.1 M Tris buffer at pH 7.0 show a much slower exchange rate for mutant rHbs and Hb A. One of the mutants, rHb (beta131Gln --> Asn), shows the conformational exchange of its interface histidines, and exchange rate measurements have been attempted. We have also conducted studies on the reactivity of the SH group of beta93Cys (a residue located in the region of the alpha(1)beta(2)-subunit interface) toward p-mercuribenzoate, and our results show that low-oxygen-affinity rHbs have a more reactive beta93Cys than Hb A in the CO form. Our results indicate that there is communication between the alpha(1)beta(1)- and alpha(1)beta(2)-subunit interfaces, and a possible communication pathway for the cooperative oxygenation of Hb A that allows the alpha(1)beta(1)-subunit interface to modulate the functional properties in conjunction with the alpha(1)beta(2) interface is proposed.     

Assessing the native state conformational distribution of ubiquitin by peptide acidity

At equilibrium, every energetically feasible conformation of a protein occurs with a non-zero probability. Quantitative analysis of protein flexibility is thus synonymous with determining the proper Boltzmann-weighting of this conformational distribution. The exchange reactivity of solvent-exposed amide hydrogens greatly varies with conformation, while the short-lived peptide anion intermediate implies an insensitivity to the dynamics of conformational motion. Amides that are well-exposed in model conformational ensembles of ubiquitin vary a million-fold in exchange rates which continuum dielectric methods can predict with an rmsd of 3. However, the exchange rates for many of the more rarely exposed amides are markedly overestimated in the PDB-deposited 2K39 and 2KN5 ubiquitin ensembles, while the 2NR2 ensemble predictions are largely consistent with those of the Boltzmann-weighted conformational distribution sampled at the level of 1%. The correlation between the fraction of solvent-accessible conformations for a given amide hydrogen and the exchange rate constant for that residue provides a useful monitor of the degree of completeness with which a given ensemble has sampled the energetically accessible conformational space. These exchange predictions correlate with the degree to which each ensemble deviates from a set of 46 ubiquitin X-ray structures. Kolmogorov-Smirnov analysis for the distribution of intra- and inter-ensemble pairwise structural rmsd values assisted the identification of a subensemble of 2K39 that eliminates the overestimations of hydrogen exchange rates observed for the full ensemble. The relative merits of incorporating experimental restraints into the conformational sampling process are compared to using these restraints as filters to select subpopulations consistent with the experimental data.