Hydrogen exchange kinetics of surface peptide amides in bovine pancreatic trypsin inhibitor

The acid and base catalytic rate constants, kH, obs and kOH, obs and the pH at the minimum rate, pHmin, of 25 rapidly exchanging protons in bovine pancreatic trypsin inhibitor have been determined. Here we report the labeling procedure giving 1H nuclear magnetic resonance spectral resolution of seven additional rapidly exchanging NH protons and the pH dependence of their chemical shifts. Values of kH,obs kOH,obs and pHmin are given for Ala16, Gly28 and Arg53 NH groups, the only backbone amide protons with static accessibility of more than zero in the crystal structure not previously reports, and for Gly56 NH, buried at the C terminus of an alpha-helix. All four protons reported here have pH min greater than or equal to 3. Conclusions of the previous study predict that peptide protons with pHmin higher than those of model compounds have greater static accessibility of the peptide O than of the peptide N atom. The locations in the crystal structure of the four NH groups whose exchange rates are reported here are in qualitative agreement with these predictions. The ionic strength dependence of Ala16 at pH 5.5 shows a sharp increase in the exchange rate with decreasing salt concentration, as expected for base-catalyzed exchange in a positive electrostatic field.     

On the pH dependence of amide proton exchange rates in proteins.

We have analyzed the pH dependencies of published amide proton exchange rates (kex) in three proteins: bovine pancreatic trypsin inhibitor (BPTI), bull seminal plasma proteinase inhibitor IIA (BUSI IIA), and calbindin D9K. The base-catalyzed exchange rate constants (kOH) of solvent exposed amides in BPTI are lower for residues with low peptide carbonyl exposure, showing that the environment around the carbonyl oxygen influences kOH. We also examined the possible importance of an exchange mechanism that involves formations of imidic acid intermediates along chains of hydrogen-bonded peptides in the three proteins. By invoking this "relayed imidic acid exchange mechanism," which should be essentially acid-catalyzed, we can explain the surprisingly high pHmin (the pH value at which kex reaches a minimum) found for the non-hydrogen-bonded amide protons in the beta-sheet in BPTI. The successive increase of pHmin along a chain of hydrogen-bonded peptides from the free amide to the free carbonyl, observed in BPTI, can be explained as an increasing contribution of the proposed mechanism in this direction of the chain. For BUSI IIA (pH 4-5) and calbindin D9K (pH 6-7) the majority of amide protons with negative pH dependence of kex are located in chains of hydrogen-bonded peptides; this situation is shown to be consistent with the proposed mechanism.     

pH dependence of hydrogen exchange from backbone peptide amides in apamin

The kinetics of hydrogen exchange of the 11 most protected backbone amides of bee venom apamin have been measured between pH 1 and pH 8.5 by using time-resolved and saturation-transfer NMR spectroscopy. The five amides most protected from base-catalyzed exchange, those of residues 5 and 12-15, show highly correlated exchange behavior in the base-catalyzed regime. It is proposed that the intramolecular hydrogen bonds stabilizing these amides define a stable cooperative unit of secondary structure in apamin (a C-terminal helix and an N-terminal beta-turn). This conformational unit is further stabilized (by 5-6 kJ mol-1) on titration of the Glu-7 side-chain carboxyl group. The relative contributions of specific intramolecular interactions to this conformational stabilization are estimated. The pHminima in the pH-dependent single amide exchange curves are compared with values predicted by correcting for sequence-dependent contributions to amide exchange rates [Molday, R. S., Englander, S. W., & Kallen, R. G. (1972) Biochemistry 11, 150-158]. The lack of correlation suggests that the "open" conformers from which amide exchange occurs are nonrandom. This conclusion is dependent on the assumption that acid-catalyzed exchange occurs via N-protonation so that residual conformational effects on exchange rates in the open conformers will affect acid- and base-catalyzed rates in approximately equal and opposite ways. A strong correlation between the measured pHminima and the amide proton chemical shifts is observed, however, and this may be most easily accommodated if acid-catalyzed exchange occurs by the imidic acid mechanism (via amide O-protonation).     

Rapid amide proton exchange rates in peptides and proteins measured by solvent quenching and two-dimensional NMR.

In an effort to develop a more versatile quenched hydrogen exchange method for studies of peptide conformation and protein-ligand interactions, the mechanism of amide proton exchange for model peptides in DMSO-D2O mixtures was investigated by NMR methods. As in water, H-D exchange rates in the presence of 90% or 95% DMSO exhibit characteristic acid- and base-catalyzed processes and negligible water catalysis. However, the base-catalyzed rate is suppressed by as much as four orders of magnitude in 95% DMSO. As a result, the pH at which the exchange rate goes through a minimum is shifted up by about two pH units and the minimum exchange rate is approximately 100-fold reduced relative to that in D2O. The solvent-dependent decrease in base-catalyzed exchange rates can be attributed primarily to a large increase in pKa values for the NH group, whereas solvent effects on pKW seem less important. Addition of toluene and cyclohexane resulted in improved proton NMR chemical shift dispersion. The dramatic reduction in exchange rates observed in the solvent mixture at optimal pH makes it possible to apply 2D NMR for NH exchange measurements on peptides under conditions where rates are too rapid for direct NMR analysis. To test this solvent-quenching method, melittin was exchanged in D2O (pH 3.2, 12 degrees C), aliquots were quenched by rapid freezing, lyophilized, and dissolved in quenching buffer (70% DMSO, 25% toluene, 4% D2O, 1% cyclohexane, 75 mM dichloroacetic acid) for NMR analysis. Exchange rates for 21 amide protons were measured by recording 2D NMR spectra on a series of samples quenched at different times. The results are consistent with a monomeric unfolded conformation of melittin at acidic pH. The ability to trap labile protons by solvent quenching makes it possible to extend amide protection studies to peptide ligands or labile protons on the surface of a protein involved in macromolecular interactions.     

Solvent exchange of buried water and hydrogen exchange of peptide NH groups hydrogen bonded to buried waters in bovine pancreatic trypsin inhibitor

Solvent exchange of 18O-labeled buried water in bovine pancreatic trypsin inhibitor (BPTI), trypsin, and trypsin-BPTI complex is measured by high-precision isotope ratio mass spectrometry. Buried water is labeled by equilibration of the protein in 18O-enriched water. Protein samples are then rapidly dialyzed against water of normal isotope composition by gel filtration and stored. The exchangeable 18O label eluting with the protein in 10-300 s is determined by an H2O-CO2 equilibration technique. Exchange of buried waters with solvent water is complete before 10-15 s in BPTI, trypsin, and BPTI-trypsin, as well as in lysozyme and carboxypeptidase measured as controls. When in-exchange dialysis and storage are carried out at pH greater than or equal to 2.5, trypsin-BPTI and trypsin, but not free BPTI, have the equivalent of one 18O atom that exchanges slowly (after 300 s and before several days). This oxygen is probably covalently bound to a specific site in trypsin. When in-exchange dialysis and storage are carried out at pH 1.1, the equivalent of three to seven 18O atoms per molecule is associated with the trypsin-BPTI complex, apparently due to nonspecific covalent 18O labeling of carboxyl groups at low pH. In addition to 18O exchange of buried waters, the hydrogen isotope exchange of buried NH groups H bonded to buried waters was also measured. Their base-catalyzed exchange rate constants are on the order of NH groups that in the crystal are exposed to solvent (static accessibility greater than 0) and hydrogen-bonded main chain O, and their pH min is similar to that for model compounds. The pH dependence of their exchange rate constants suggests that direct exchange with water may significantly contribute to their observed exchange rate.     

Hydrogen kinetics of peptide amide protons at the bovine pancreatic trypsin inhibitor protein-solvent interface

Hydrogen exchange rate constants of the 25 most rapidly exchanging peptide amide protons in bovine pancreatic trypsin inhibitor have been determined over a range of pH that spans pH min, the pH of minimum rate. Most of these are on the protein surface, exposed to solvent and not hydrogen bonded in the crystal structure. Contrary to commonly held assumptions, the exchange kinetics of surface NH groups are not equivalent to the kinetics of NH groups in peptides in the extended configuration. All surface NH groups exchange more slowly than NH groups in model peptides, with rate constants distributed over a range of more than two orders of magnitude. In addition, their pH min values vary widely. For most of the surface NH groups, pH min is lower than in model compounds and, for several, pH min is less than 1. These results indicate that the local environment of the surface peptide groups when the exchange event occurs is very different from that of extended peptides. Analysis based on consideration of an O-protonation mechanism for acid catalysis and of electrostatic effects on exchange kinetics further indicates (see the accompanying paper) that, in general, exchange of surface NH groups occurs from a conformation of the protein approximated by the crystal structure. The 1H-2H exchange rate constants were measured from 300 MHz nuclear magnetic resonance spectra in which assigned surface N1H resonances are resolved by the use of partially deuterated protein samples. A marked pH dependence of the chemical shifts observed in the pH range 1 to 4.5 for several surface NH groups reflects the titration of nearby carboxyl groups.     

Hydrogen-deuterium exchange in nucleosides and nucleotides. A mechanism for exchange of the exocyclic amino hydrogens of adenosine.

The pH dependence of the apparent first-order rate constant for the exchange of the exocyclic amino hydrogens of adenosine with deuterium from the solvent was measured by stopped-flow ultraviolet spectroscopy. This dependence shows acid catalysis, base catalysis, and spontaneous exchange at neutral pH values. A study of the effect of several buffers on the rates of exchange reveals both general acid and general base catalytic behavior for the exchange process. We propose a general mechanism for the exchange which requires N-1 protonated adenosine as an intermediate for the acid-catalyzed exchange and amidine anion for the base-catalyzed exchange. In both cases the rate-limiting step is the base-catalyzed abstraction of a proton from the exocyclic amino moiety. Evaluation of the rate constants predicts the equilibrium for the exocyclic amino/imino tautomers to be 6.3 times 10(3):1.     

The prohormone proenkephalin possesses differential conformational features of subdomains revealed by rapid H-D exchange mass spectrometry

Active enkephalin and related peptide hormones or neurotransmitters are generated by proteolytic processing of inactive prohormone precursors. Little is known about the relative accessibilities of prohormone cleavage sites and conformations of subdomains that undergo proteolytic processing. Therefore, this study investigated the conformational features of the prohormone proenkephalin (PE) by rapid hydrogen-deuterium exchange mass spectrometry (DXMS). DXMS analyzes rates of hydrogen exchange of the polypeptide backbone of PE with deuterium from D(2) O (heavy water) by mass spectrometry, accomplished at sub-second and multisecond time periods. Results showed differential accessibilities of cleavage sites and adjacent subdomains of PE to the aqueous environment. Importantly, protease cleavage sites of PE with greater relative accessibilities correspond to sites most readily cleaved by processing proteases to generate active peptide neurotransmitters. For comparison, peptides derived from PE (by pepsin digestion) displayed greater accessibility to the solvent environment, illustrated by their higher rates of H-D exchange compared to that of intact PE protein. The more limited H-D exchange accessibilities of PE protein, compared to peptides derived from PE, indicate that PE possesses tertiary conformation. These results demonstrate that differential tertiary conformations of PE subdomains undergo ordered proteolytic processing to generate active enkephalin peptides for cell-cell communication in the nervous and endocrine systems.     

Hydrogen--deuterium exchange kinetics of the amide protons of oxytocin studied by nuclear magnetic resonance

The contribution of intramolecular hydrogen bonding to the solution structure of oxytocin was evaluated by study of amide hydrogen exchange rates in D2O by Fourier transform 1H NMR spectroscopy. Resolution enhancement filtering was employed in the determination of individual pseudo-first-order rate constants. Apparent barriers to exchange of 0.5 and 0.6 kcal mol-1 were measured for Asn5 and Cys6 peptide NH, respectively. The slowing is best explained by steric hindrance to solvent access in the case of Asn5, while for the Cys6 participation in a weak intramolecular hydrogen bond is possible. Fourfold acceleration of base-catalyzed exchange was observed for Tyr2 NH; it is proposed that this is the result of electronic effects induced by hydrogen bonding of Cys1 C=0, either to Cys6 NH or to the N-terminal amino group. Exchange proceeds near the random coil limit for each of the remaining residues. Comparison with exchange data for the model tripeptide N-acetyl-L-prolyl-L-leucylglycinamide demonstrates no evidence of noncovalent association of the tocin ring with the tripeptide tail of the hormone.     

Hydratation of the rubredoxin polypeptide chain from data on the spatial structure

The hydration water distribution around the main chain protein rubredoxin has been analysed using the crystal data at high resolution obtained earlier. The analysis was based on the consideration of all nearest neighbour atoms around the N and O atoms of peptide groups. The atoms which can form hydrogen binds were the subject of final analysis. The nitrogen atom of a peptide NH group has only one vacancy for neighbours. The oxygen atom of a peptide CO group has one, two or more neighbours, some of them are oxygen-water atoms. About 27% of NH and 53% of CO peptide groups are hydrated, that corresponds to 0.12 H2O per gram of protein. A detailed analysis shows that NH and CO groups of the main chain are hydrated according to the principle of maximum possible in situ saturation of hydrogen bonds. Thus the peptide groups incorporated in the peptide hydrogen bond network were not hydrated as a rule. Consequently, for rubredoxin a pleated sheet region, some regions for the large and small main chain loops, and Fe-containing pocket are not hydrated. A method for evaluation of the main chain hydration is proposed when the coordinates of protein atoms are available.     

Hydrogen bonding monitored by deuterium isotope effects on carbonyl 13C chemical shift in BPTI: intra-residue hydrogen bonds in antiparallel beta-sheet.

Deuterium isotope effects on carbonyl 13C magnetic shielding were measured for the backbone carbonyl groups in BPTI (basic pancreatic trypsin inhibitor), and interpreted as a measure of hydrogen bond energies. The effects originate from peptide amide proton deuterium substitution and were observed on carbonyl carbons separated by two or three covalent bonds from the amide H/D. Two-bond isotope effects depend on the energy of the hydrogen bond donated by NH/D. Calibration of the effect with model compound data leads to hydrogen bond enthalpies less than 4.7 kcal/mol. Isotope effects over three bonds from the amide H/D to the carbonyl carbon of the same amino acid residue are observed for seven carbonyl resonances in BPTI. The three-bond isotope effects are highly related to the various backbone conformations. The largest effects are observed for residues with an approximate syn- periplanar conformation of the H-N-C alpha-C = O atoms, as realized for many residues in the BPTI antiparallel beta-sheet. The residues that show measurable three-bond effects have unusually short distances between H and O. The size of this effect decreases rapidly with increased O..H distance in the open five-membered ring. This observation suggests appreciable interactions in these rings.     

The environment of amide groups in protein–ligand complexes: H-bonds and beyond

A comprehensive structural analysis of interactions involving amide NH and C=O groups in protein-ligand complexes has been performed based on 3,275 published crystal structures (resolution < or =2.5 A). Most of the amide C=O and NH groups at the protein-ligand interface are highly buried within the binding site and involved in H-bonds with corresponding counter-groups. Small percentages of C=O and NH groups are solvated or embedded in hydrophobic environments. In particular, C=O groups show a higher propensity to be solvated or embedded in a hydrophobic environment than NH groups do. A small percentage of carbonyl groups is involved in weak hydrogen bonds with CH. Cases of dipolar interactions, involving carbonyl oxygen and electrophilic carbon atoms, such as amide, amidinium, guanidium groups, are also identified. A higher percentage of NH are in contact with aromatic carbons, interacting either through hydrogen bonds (preferably with the NH group pointing towards a ring carbon atom) or through stacking between amide plane and ring plane. Comprehensive studies such as the present one are thought to be important for future improvements in the molecular design area, in particular for the development of new scoring functions. [Figure: see text].     

Thiazolium C(2)-proton exchange: structure-reactivity correlations and the pKa of thiamin C(2)-H revisited

Rate constants for C(2)-proton exchange from thiamin, N(1')-methylthiamin, and several 3-substituted-4-methylthiazolium ions catalyzed by D2O and deuterioxide ion were determined by 1H NMR at 30 degrees C and ionic strength 2.0 M. Values of pKa for the thiazolium ions, including thiamin itself, were found to be in the range pKa = 17-19; the pKa values for N(1')-protonated thiamin and free thiamin C(2)-H in H2O are 17.7 and 18.0, respectively. The pKa value for N(1')-protonated thiamin was calculated from the observed rate constant for the pD-independent reaction with D2O after correction for a secondary solvent deuterium isotope effect of kH2O/kD2O = 2.6. The pKa value for free thiamin was calculated from the rate constant for catalysis by OD- after correction by a factor of 3.3 = 8/2.4 for an 8-fold negative deviation of kOD from the Br?nsted plot of slope 1.0 for general base catalysis and a secondary solvent isotope effect of kOD/kOH = 2.4. Values of k-a = 2 X 10(10) and 3 X 10(9) M-1 s-1 were assumed for diffusion-controlled protonation of the C(2) ylide in the reverse direction by H3O+ and H2O, respectively. The Hammett rho I value for the exchange reaction catalyzed by deuterioxide ion or D2O is 8.4 +/- 0.2. There is no positive deviation of the rate constants for free or N(1')-substituted thiamin analogues in either Hammett correlation. This shows that the aminopyrimidinyl group does not provide significant intramolecular catalysis of nonenzymic C(2)-proton removal in the coenzyme.     

Peptide--water association in peptide crystals.

The structures of 37 peptide crystals, containing 78 water-peptide hydrogen bonds and 77 other hydrogen bonds involving water, were surveyed to identify the geometry of peptide backbone hydration. In the sample, hydration of peptide carbonyl occurred more frequently than hydration of peptide N--H. The most probable value of the C'=O ... O water angle was near 138 degrees, considerably greater than the 120 degrees to the axis of a lone electron pair on the carbonyl oxygen. Associated water oxygens tended to be in the plane of the peptide bond, bui--H and Ci+1=O atoms, was common in glycine-containing cyclic hexapeptides. The distribution of angles between two hydrogen bonds at a single water molecule, as defined by the three nonhydrogen atoms involved, was centered near the tetrahedral angle.     

The interpretation of protein structures: estimation of static accessibility

A program is described for drawing the van der Waal's surface of a protein molecule. An extension of the program permits the accessibility of atoms, or groups of atoms, to solvent or solute molecules of specified size to be quantitatively assessed. As defined in this study, the accessibility is proportional to surface area. The accessibility of all atoms in the twenty common amino acids in model tripeptides of the type Ala-X-Ala are given for defined conformation. The accessibilities are also given for all atoms in ribonuclease-S, lysozyme and myogoblin. Internal cavities are defined and discussed. Various summaries of these data are provided. Forty to fifty per cent of the surface area of each protein is occupied by non-polar atoms. The actual numerical results are sensitive to the values chosen for the van der Waal's radii of the various groups. Since there is uncertainty over the correct values for these radii, the derived numbers should only be used as a qualitative guide at this stage.     

Acid proteases from species of Mucor. IV. Hydrogen-tritium exchange of Mucor miehei protease.

The kinetics of hydrogen-tritium exchange were studied in the range pH-3 for both the fully and partially tritiated protein. Exchange constants for an intermediate class and slow class of hydrogens were determined and found to give a parabolic curve characteristic of acid and base catalysis about the observed pHmin of 4.03. The anomalous rate retardation on the acid portion of the curve was attributed to electrostatic interactions which could be evaluated quantitatively from the titration data. Partial tritation and pH cross-over experiments indicated that the rank order was pH-independent thus eliminating the possiblitity of a major conformational change. Consequently, the data are most likely explicable in terms of restricted solvent accessibility.     

Atom-wise statistics and prediction of solvent accessibility in proteins

In this work, we explore a novel method to broaden the scope of sequence-based predictions of solvent accessibility or accessible surface area (ASA) to the atomic level. All 167 heavy atoms from the 20 types of amino acid residues in proteins have been studied. An analysis of ASA distribution of these atomic groups in different proteins has been performed and rotamer-style libraries have been developed. We observe that the ASA of some atomic groups (e.g., backbone C and N atoms) can be estimated from the sequence environment within a mean absolute error of 2-3 angstroms(2). However, some side chain atoms such as CG in Pro, NH1 in Arg and NE2 in Gln show a strong variability making it more difficult to estimate their ASA from sequence environment. In general, the prediction of ASA becomes more difficult for atomic positions at the side chain extremities of long amino acid residues (aromatic side chain terminals being the exception). Several atomic groups are frequently exposed to solvent. Some of them have a bimodal distribution, suggesting two stable conformations in terms of their solvent exposure. More detailed understanding and prediction of solvent accessibility, i.e., at an atomic level is expected to help in bioinformatics approaches to structure prediction, functional relevance of atomic solvent accessibilities and other interaction analyses.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

General Acid-Base Catalysis in Nucleobase Amino Proton Exchange: Cytidine

Abstract

A useful property of DMSO solvent has been exploited to reveal a new catalytic route for cytidine amino proton exchange, relevant to exchange in the macromolecular state, but hidden in aqueous solution. Additional exchange mechanisms in aqueous monomeric cytidine (and adenosine) are obscured by the formation of a fast-exchanging endocyclic-protonated intermediate, which dominates the kinetics. Endocyclic nucleobase protonation could be circumvented in the presence of buffer conjugate acid by the use of DMSO/water solvent, permitting the first unequivocal observation buffer acid-catalyzed exchange from the neutral, unprotonated nucleobase, i.e., general acid catalysis. Because buffer ionization is greatly reduced in DMSO through anion desolvation, nucleobase protonation is supressed m the presence of buffer acid. Evidence is presented to describe this catalytic route as one involving hydrogen bond formation between the buffer acid and the endocyclic protonation site, C(N-3). Since this same configuration is found in Watson-Crick hydrogen bonding, experiments are presented to demonstrate faster cytidine amino proton exchange with the formation of the G-C base pair in DMSO. The importance of this mechanism in past aqueous monomer studies and in the interpretation of macromolecular (DNA) hydrogen exchange is discussed.     

Probing protein structure by solvent perturbation of NMR spectra: the surface accessibility of bovine pancreatic trypsin inhibitor.

In the absence of specific interactions, the relative attenuation of protein NMR signals due to added stable free radicals such as TEMPOL should reflect the solvent accessibility of the molecular surface. The quantitative correlation between observed attenuation and surface accessibility was investigated with a model system, i.e., the small protein bovine pancreatic trypsin inhibitor. A detailed discussion is presented on the reliability and limits of the approach, and guidelines are provided for data acquisition, treatment, and interpretation. The NMR-derived accessibilities are compared with those obtained from x-ray diffraction and molecular dynamics data. Although the time-averaged accessibilities from molecular dynamics are ideally suited to fit the NMR data, better agreement was observed between the paramagnetic attenuations of the fingerprint cross-peaks of homonuclear proton spectra and the total NH and H alpha accessibilities calculated from x-ray coordinates, than from time-averaged molecular dynamics simulations. In addition, the solvent perturbation response appears to be a promising approach for detecting the thermal conformational evolution of secondary structure elements in proteins.