Nuclear magnetic resonance of labile protons in the basic pancreatic trypsin inhibitor

NMR study of the active site of shark met-cyano myoglobins

The myoglobins from the sharks Galeorhinus japonicus and Musterus japonicus possess a distal Gln-E7 instead of the usually found His-E7. The met-cyano form of these shark myoglobins has been studied by 1H- and 15N-NMR in order to gain insight into the functional properties of the Gln-E7. The analysis of paramagnetic relaxation has provided the assignment of the resonance arising from one of the Gln-E7 N epsilon H labile protons, whilst the rate of its chemical exchange has been analyzed in detail by using a saturation transfer method. The hydrogen-bonding interaction between this proton and Fe-bound-CN(-) has been clearly manifested in the hyperfine shift of the Gln-E7 N epsilon H proton resonance as well as its chemical exchange behavior. The resonances of the Gln-E7 side-chain non-labile protons have been partly assigned on a basis of both scalar and dipolar connectivities. The analysis of the dipolar connectivities among the side-chain protons and the iron-proton distances determined from their paramagnetic relaxation rate has revealed that the side chain adopts a conformation with its carbonyl oxygen oriented away from the heme. Although (1)H-NMR spectra of these two myoglobins are essentially similar, a relatively large difference in the shift of Gln-E7 N(epsilon)H proton and Fe-bound C15 N- resonances between the two has been observed which is attributed to a differential hydrogen-bonding interaction between these proteins. The present study demonstrates the sensitivity of NMR parameters to the hydrogen-bonding interaction between coordinated ligand and a distal amino-acid side chain in paramagnetic hemoproteins.

Internal packing conditions and fluctuations of amino acid residues in globular proteins

In order to investigate the environmental conditions of amino acid residues in protein molecules, four kinds of packing studies (atomic, geometric, hydrophobic and hydration) were formulated and tested on two proteins; bovine pancreatic trypsin inhibitor (BPTI) and bovine pancreatic ribonuclease S (RNase S). The inter-relationship of these packings on the fluctuations of amino acid residues was analysed by comparing the packing results with the dynamical studies, such as the root-mean-square-deviation values of atomic displacements obtained from the trajectories of molecular dynamics simulation, temperature factor information from crystal structures and residue fluctuations in proteins from continuum model. These analyses yield information about the most fluctuating and most stabilizing residue sites. Comparison of the results obtained by these methods indicate a good agreement, specifying an inverse correlation between the residue packing and fluctuations. This kind of study is helpful in identifying the specific residue sites such as nucleation, receptor binding and antigenic determining sites which in a way indirectly correlates with the functional residues in protein molecules.

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.

Acid destabilization of a triple-helical peptide model of the macrophage scavenger receptor

Electrostatic interactions were studied in a triple-helical peptide, (POG)3PKGQKGEKG(POG)4, which contains a lysine-rich 9 residue sequence from the collagen-like domain of the macrophage scavenger receptor (MSR). This peptide adopts a stable triple-helical conformation only when the pH is higher than 4.5, corresponding to ionization of the Glu side chain. Modeling shows Glu forms ion pairs with one of the Lys residues, stabilizing the structure. Previously studied collagen-like peptides show relatively small contributions of electrostatic interactions to stability. The large magnitude of the pH mediated structural changes seen for this peptide suggests that specific placement of charged residues in the triple-helix conformation can generate strong electrostatic interactions.

Determination of the rotational dynamics and pH dependence of the hydrogen exchange rates of the arginine guanidino group using NMR spectroscopy

The dynamic behaviour of the guanidino group of arginine has been investigated quantitatively with the intention of providing a set of basis values for the interpretation of data acquired for arginine residues in proteins. At room temperature, a single broad resonance line is observed for the four η-NH(2) protons. Upon cooling the sample (≈10°C at 500 MHz), two η-NH(2) proton resonances are resolved which were shown by HMQC spectroscopy to be the result of slowed rotation about the N(ε)-C(ζ) partial double bond. The flip rate (k(NC)) about the N(ε)-C(ζ) bond was measured as a function of temperature using line-shape analysis of both (1)H and (15)N NMR spectra; at 25°C, k(NC) is between 900 and 1000 s(-1). The exchange broadening, due to N(ε)-C(ζ) bond flips, typically results in weak or missing signals for the η-NH(2) protons of arginine residues in HMQC or INEPT experiments recorded at room temperature, unless the motion is restricted in some way. In a related series of experiments, the pH dependence of the hydrogen exchange rates of the ε-NH and η-NH(2) protons of arginine was measured using saturation transfer (1)H NMR spectroscopy and compared with the equivalent NH(2) protons of the guanidinium ion. As expected, OH ion catalysis dominates over most of the pH range and proceeds at a rate close to the diffusion limit for both types of proton (k(OH)=2×10(9)-1×10(10) M(-1)s(-1), depending on conditions). At low pH values, however, catalysis by H(3)O(+) becomes important and leads to characteristic rate minima in the exchange versus pH profiles. Acid catalysis is significantly more effective for the η-NH(2) protons than for the ε-NH proton; at low ionic strength (50 mM KCl) the rate minima occurred at pH 3.6 and 2.3, respectively. Under these conditions, acid-catalysed rate constants (k(H)) of 706 M(-1)s(-1) (η-NH(2)) and 3 M(-1)s(-1) (ε-NH) were obtained at 25°C. At high ionic strength (1 M KCl) the rate of OH(-) ion catalysis is decreased slightly, whereas the H(3)O(+)-catalysed rate is unchanged. The k(OH) value of the free guanidinium ion is identical to that of the η-NH(2) protons but acid catalysis occurs less easily, leading to a rate minimum at pH 3.3.

Structure and dynamics of the neutrophil defensins NP-2, NP-5, and HNP-1: NMR studies of amide hydrogen exchange kinetics

The exchange kinetics for the slowly exchanging amide hydrogens in three defensins, rabbit NP-2, rabbit NP-5, and human HNP-1, have been measured over a range of pH at 25 degrees C using 1D and 2D NMR methods. These NHs have exchange rates 10(2) to 10(5) times slower than rates from unstructured model peptides. The observed distribution of exchange rates under these conditions can be rationalized by intramolecular hydrogen bonding of the individual NHs, solvent accessibility of the NHs, and local fluctuations in structure. The temperature dependencies of NH chemical shifts (NH temperature coefficients) were measured for the defensins and these values are consistent with the defensin structure. A comparison is made between NH exchange kinetics, NH solvent accessibility, and NH temperature coefficients of the defensins and other globular proteins. Titration of the histidine side chain in NP-2 was examined and the results are mapped to the three-dimensional structure.

Structural studies of two antiaggregant RGDW peptides by 1H and 13C NMR

The structural features of Arg-Gly-Asp-related sequences have been investigated by 1H and 13C NMR. Two linear peptides which inhibit platelet aggregation with a high efficiency have been studied: D-Arg-Gly-Asp-Trp and L-Arg-Gly-Asp-Trp. Analysis of pH titration effects, amide proton exchange rates and inter-proton distances obtained from ROESY spectra suggest that these small fragments predominantly adopt a type II' beta-turn structure in solution. Folding features of a non-active cyclic peptide based on the same sequence (cyclo-[Arg-Gly-Asp-Trp]2) have also been investigated. The biological relevance of these structures is discussed.

A conformational study of Lys-Arg-Asp-Ser and analogs, a series of potent antithrombotic peptides. An approach based on simulated annealing and 1H NMR

Simulated annealing techniques were used to explore the conformational space of the potent antithrombotic peptide L.Lys-L.Arg-L.Asp-L.Ser (KRDS) and of two analogs: D.Lys-L.Arg-L.Asp-L.Ser (KDRDS), which is inactive, and L.Lys-L.Arg-L.Glu-L.Glu (KREE), which exhibits a strong biological activity. For each peptide, a set of initial conformations was generated and submitted to simulated annealing, including a heating to 1000 K followed by a cooling to 300 K. 200 resulting conformations of each compound were analyzed and classified according to the network of electrostatic interactions involving charged side chains and charged C- and N-terminal groups. A reduced number of conformational classes was obtained and conformations corresponding to predominant classes were found to be in qualitative agreement with structural parameters deduced from 1H NMR spectra. A comparison between the classes of the active and non active peptide was achieved. Some conformations were found to be specific of active peptides.

Electronic effects on the fluorescence of tyrosine in small peptides

It is shown for a series of tyrosine-derivatives and tyrosine-containing peptides that the amide group in combination with electron-withdrawing substituents quenches the fluorescence of the phenol moiety. The ammonium group has the strongest electron-withdrawing effect and thus the largest influence on the quenching rate. The peptide group itself does not quench the fluorescence. In a series of peptides with an increasing number of alanines the decreasing quenching efficiency of the peptide group due to the greater distance of the ammonium group is demonstrated. In tyrosine-containing di- and tripeptides a linear correlation between the 13C-NMR chemical shift delta of the C alpha atom of various aliphatic amino acids and the fluorescence-quenching constant confirms the hypothesis that electron-withdrawing and -donating groups are modulating the fluorescence-quenching efficiency of the peptide group. In small peptides the fluorescence lifetime of tyrosine is characteristic for the neighboring amino acids. Using model substances the redox properties of a peptide group and the phenol ring were studied electrochemically. The highest occupied molecular orbital of the tyrosine (1.4 V vs saturated calomel electrode [SCE]) and the lowest unoccupied molecular orbital of the peptide group (-3.12 V vs SCE) have appropriate energies for a photoinduced electron transfer reaction. For solute-quenching experiments quencher molecules can be systematically selected.

On the calculation of pKas in proteins

This paper describes a general method to calculate the pKas of ionizable groups in proteins. Electrostatic calculations are carried out using the finite difference Poisson-Boltzmann (FDPB) method. A formal treatment of the calculation of pKas within the framework of the FDPB method is presented. The major change with respect to previous work is the specific incorporation of the complete charge distribution of both the neutral and charged forms of each ionizable group into the formalism. This is extremely important for the treatment of salt bridges. A hybrid statistical mechanical/Tanford-Roxby method, which is found to be significantly faster than previous treatments, is also introduced. This simplifies the problem of summing over the large number of possible ionization states for a complex polyion. Applications to BPTI and serine proteases suggest that the calculations can be quite reliable. However, the necessity of including bound waters in the treatment of the Asp-70... His-31 salt bridge in T4 lysozyme and experience with other proteins suggest that additional factors ultimately need to be considered in a comprehensive treatment of pKas in proteins.

Backbone dynamics of the glucocorticoid receptor DNA-binding domain

The extent of rapid (picosecond) backbone motions within the glucocorticoid receptor DNA-binding domain (GR DBD) has been investigated using proton-detected heteronuclear NMR spectroscopy on uniformly 15N-labeled protein fragments containing the GR DBD. Sequence-specific 15N resonance assignments, based on two- and three-dimensional heteronuclear NMR spectra, are reported for 65 of 69 backbone amides within the segment C440-A509 of the rat GR in a protein fragment containing a total of 82 residues (MW = 9200). Individual backbone 15N spin-lattice relaxation times (T1), rotating-frame spin-lattice relaxation times (T1 rho), and steady-state (1H)-15N nuclear Overhauser effects (NOEs) have been measured at 11.74 T for a majority of the backbone amide nitrogens within the segment C440-N506. T1 relaxation times and NOEs are interpreted in terms of a generalized order parameter (S2) and an effective correlation time (tau e) characterizing internal motions in each backbone amide using an optimized value for the correlation time for isotropic rotational motions of the protein (tau R = 6.3 ns). Average S2 order parameters are found to be similar (approximately 0.86 +/- 0.07) for various functional domains of the DBD. Qualitative inspection as well as quantitative analysis of the relaxation and NOE data suggests that the picosecond flexibility of the DBD backbone is limited and uniform over the entire protein, with the possible exception of residues S448-H451 of the first zinc domain and a few residues for which relaxation and NOE parameters were not obtained. in particular, we find no evidence for extensive rapid backbone motions within the second zinc domain. Our results therefore suggest that the second zinc domain is not disordered in the uncomplexed state of DBD, although the possibility of slowly exchanging (ordered) conformational states cannot be excluded in the present analysis.

Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.

Sequence specific 1H-NMR assignments and secondary structure of a carboxy-terminal functional fragment of apolipoprotein CII

The structural properties of a synthetic fragment of human apolipoprotein CII (apoCII) has been studied by circular dichroism and proton nuclear magnetic resonance. The fragment corresponds to the carboxy-terminal 30 amino acid residues and retains the ability of apoCII to activate lipoprotein lipase. Like native apoCII, the fragment has a tendency to self-associate in pure aqueous solution. Addition of 1,1,1,3,3,3-hexafluoro-2-isopropanol to aqueous solvent dissolves the aggregates and leads to an increase in the alpha-helical content of the peptide, probably by stabilizing transient helical structures. The resonances in the 1H-NMR spectrum of the fragment in 35% (CF3)2CHOH were assigned through standard procedures from nuclear Overhauser enhancement spectroscopy, correlated spectroscopy and total correlated spectroscopy experiments. The NMR data indicates the formation of a stable alpha helix spanning Ile66-Gly77. Another alpha helical turn may be formed between Lys55 and Ala59 and possibly span even further towards the carboxyl terminus. These structural elements are different from those previously predicted for this part of the sequence of apoCII.

Conformation and mobility of the RNA-binding N-terminal part of the intact coat protein of cowpea chlorotic mottle virus. A two-dimensional proton nuclear magnetic resonance study

Two-dimensional proton nuclear magnetic resonance (n.m.r.) experiments were performed on the coat protein of cowpea chlorotic mottle virus (molecular mass: 20.2 kDa) present as dimer (pH 7.5) or as capsid consisting of 180 protein monomers (pH 5.0). The spectra of both dimers and capsids showed resonances originating from the flexible N-terminal region of the protein. The complete resonance assignment of a synthetic pentacosapeptide representing this N terminus made it possible to interpret the spectra in detail. The capsid spectrum showed backbone amide proton resonances arising from the first eight residues having a flexible random coil conformation, and side-chain resonances arising from the first 25 N-terminal amino acids. The dimer spectrum showed also side-chain resonances of residues 26 to 33, which are flexible in the dimer but immobilized in the capsid. The n.m.r. experiments indicated that the conformation of the first 25 amino acids of the protein in dimers and capsids is comparable to the conformation of the synthetic peptide, which alternates among extended and helical conformations on the n.m.r. time-scale. It is suggested that the alpha-helical region, situated in the region between residues 10 and 20, binds to the RNA during assembly of the virus particle.

A two-dimensional 1H-NMR study on Megasphaera elsdenii flavodoxin in the oxidized state and some comparisons with the two-electron-reduced state

Assignments for the 137 amino acid residues of oxidized Megasphaera elsdenii flavodoxin have been made using the sequential resonance assignment procedure. Great benefit was experienced from assignments of the fully reduced protein. The secondary and tertiary structures of the typical alpha/beta protein remain virtually identical on going from the oxidized to the two-electron-reduced state as judged from two-dimensional NOE spectroscopy. However, functionally important conformation changes in the flavin-binding region do occur on reduction. Considerable reduction-state-dependent chemical shift variations of protons in the immediate vicinity of the isoalloxazine moiety take place. From analysis of these shifts, it can be concluded that ring current effects of the pyrazine part of the flavin diminish on two-electron reduction.

Hydrogen exchange kinetics in a membrane protein determined by 15N NMR spectroscopy: use of the INEPT experiment to follow individual amides in detergent-solubilized M13 coat protein

The coat protein of the filamentous coliphage M13 is a 50-residue polypeptide which spans the inner membrane of the Escherichia coli host upon infection. Amide hydrogen exchange kinetics have been used to probe the structure and dynamics of M13 coat protein which has been solubilized in sodium dodecyl sulfate (SDS) micelles. In a previous 1H nuclear magnetic resonance (NMR) study [O'Neil, J. D. J., & Sykes, B. D. (1988) Biochemistry 27, 2753-2762], multiple exponential analysis of the unresolved amide proton envelope revealed the existence of two slow "kinetic sets" containing a total of about 30 protons. The slower set (15-20 amides) originates from the hydrophobic membrane-spanning region and exchanges at least 10(5)-fold slower than the unstructured, non-H-bonded model polypeptide poly(DL-alanine). Herein we use 15N NMR spectroscopy of biosynthetically labeled coat protein to follow individual, assigned, slowly exchanging amides in or near the hydrophobic segment. The INEPT (insensitive nucleus enhancement by polarization transfer) experiment [Morris, G. A., & Freeman, R. (1979) J. Am. Chem. Soc. 101, 760-762] can be used to transfer magnetization to the 15N nucleus from a coupled proton; when 15N-labeled protonated protein is dissolved in 2H2O, the INEPT signal disappears with time as the amide protons are replaced by solvent deuterons. Amide hydrogen exchange is catalyzed by both H+ and OH- ions. Base catalysis is significantly more effective, resulting in a characteristic minimum rate in model peptides at pH approximately equal to 3. Rate versus pH profiles have been obtained by using the INEPT experiment for the amides of leucine-14, leucine-41, tyrosine-21, tyrosine-24, and valines-29, -30, -31, and -33 in M13 coat protein. The valine residues exchange most slowly and at very similar rates, showing an apparent 10(6)-fold retardation over poly(DL-alanine). A substantial basic shift in the pH of the minimum rate (up to 1.5 pH units) was also observed for some residues. Possible reasons for the shift include accumulation of catalytic H+ ions at the negatively charged micelle surface or destabilization of the negatively charged transition state of the base-catalyzed reaction by either charge or hydrophobic effects within the micelle. The time-dependent exchange-out experiment is suitable for slow exchange rates (kex), i.e., less than (1-2) x 10(-4) s-1.(ABSTRACT TRUNCATED AT 400 WORDS)

Structure and dynamics in proteins of pharmacological interest

This paper begins with a brief survey of the standard nuclear magnetic resonance (NMR) method for protein structure determination in solution, which has been applied successfully with numerous globular proteins with molecular weights in the range from 3,000 to 15,000. The results obtained show that for the core of globular proteins, the quality of the structures determined in solution can be comparable to that achieved with diffraction techniques using protein single crystals. In addition, they also indicate that a complete description of proteins in solution may include short-lived, transient structural features that could be of crucial importance for the functional properties. Several supplementary NMR techniques capable of characterizing diverse aspects of flexible polypeptide chains are discussed.

Conformational studies of nikkomycin X in aqueous solution

Nikkomycin X is a peptidyl-nucleoside antibiotic that inhibits chitin synthesis in fungi. Information on its conformation in aqueous solution was obtained from pH titration studies in which the nmr shifts of exchangeable and nonexchangeable protons were monitored. These studies and nuclear Overhauser effects support an unfolded or conformationally flexible structure for the antibiotic, and the syn and anti conformations of the nucleoside moiety were found to coexist. The conformation of the ribose ring was determined using a two-state model; a slight shift from type N to type S conformers occurred as the pH was raised from 1 to 6.

A two-dimensional 1H NMR study on Megasphaera elsdenii flavodoxin in the reduced state. Sequential assignments

Assignments for the 137 amino acid residues of Megasphaera elsdenii flavodoxin in the reduced state have been made using the sequential resonance assignment procedure. Several hydroxyl and sulfhydryl protons were observed at 41 degrees C at pH 8.3. Spin systems were sequentially assigned using phase-sensitive two-dimensional-correlated spectroscopy and phase-sensitive nuclear Overhauser enhancement spectroscopy. Spectra of the protein in H2O and of protein preparations either completely or partly exchanged against 2H2O were obtained. Use of the fast electron shuttle between the paramagnetic semiquinone and the diamagnetic hydroquinone state greatly simplified the NMR spectra, making it possible to assign easily the 1H resonances of amino acid residues located in the immediate neighbourhood of the isoalloxazine ring. The majority of the nuclear Overhauser effect contracts between the flavin and the apoprotein correspond to the crystal structure of the flavin domain of Clostridium MP flavodoxin, but differences are also observed. The assignments provide the basis for the structure determination of M. elsdenii flavodoxin in the reduced state as well as for assigning the resonances of the oxidized flavodoxin.

The concentration dependence of the diffusion coefficient for bovine pancreatic trypsin inhibitor: a dynamic light scattering study of a small protein

This paper reports the use of dynamic light scattering to investigate the concentration dependence of the diffusion coefficient for bovine pancreatic trypsin inhibitor (BPTI). BPTI is a small molecular weight protein (6511 Da) that has been the subject of numerous experimental studies. In addition to addressing questions that remain in the literature concerning the aggregation behavior of BPTI, we show that dynamic light scattering can be practically applied to proteins as small as BPTI, and that it can provide a useful means of parameterizing the solution behavior for proteins. We obtained values for the apparent diffusion coefficient of BPTI as a function of concentration over a range of pH values from 2.59 to 9.92 at an ionic strength of 0.3M, and over a range of ionic strength values from 0.1 to 0.5M at a pH of 7.0. The concentration dependence is linear for nearly all the conditions examined, even up to concentrations as high as 65 mg/mL. The average diffusion coefficient obtained at infinite dilution is 14.4 +/- 0.2 x 10(-7) cm2/s. This value agrees with that expected for a BPTI monomer hydrated with less than a monolayer of water. We used the theories of Felderhof, of Batchelor, and of Phillies, along with the DLVO theory to interpret the concentration dependence of the apparent diffusion coefficient. The variations observed with pH and ionic strength can be primarily attributed to screened coulombic interactions. In addition, there is an attractive interaction that is slightly stronger than the repulsive coulombic one, and that is essentially independent of pH and ionic strength. The attractive interactions appear to arise from nonspecific van der Waals interactions and do not lead to the formation of stable aggregates of BPTI.

A proton nuclear magnetic resonance assignment and secondary structure determination of recombinant human thioredoxin

Two-dimensional 1H NMR spectroscopy has been applied to a structural analysis of the reduced form of a recombinant human thioredoxin, a ubiquitous dithiol oxidoreductase recently isolated from an immunocompetent lymphoblastoid cell line. The sequential assignment of the spectrum, including all proline residues, has been accomplished by using experiments to demonstrate through-bond and through-space connectivities. The secondary structure has been determined by a qualitative interpretation of nuclear Overhauser effects, NH exchange data, and 3JHN alpha coupling constants. The secondary structure was found to be similar to that of the X-ray structure of Escherichia coli thioredoxin, consisting of a mixed five-stranded beta-sheet surrounded by four alpha-helices. The assignment and structural characterization of human thioredoxin was facilitated by the increased resolution and sensitivity afforded by a magnetic field strength of 600 MHz and required the use of two temperatures and two pH conditions to resolve ambiguities caused by a duplication of resonances. This duplication, extending from Phe-41 to Val-59, and including Lys-3-Ile-5, Val-24, Val-25, Asn-39, and Ile-101-Glu-103, appears to be due to heterogeneity arising from the presence or absence of the N-terminal methionine.

Ligand and proton exchange dynamics in recombinant human myoglobin mutants

Site-specific mutants of human myoglobin have been prepared in which lysine 45 is replaced by arginine (K45R) and aspartate 60 by glutamate (D60E), in order to examine the influence of these residues and their interaction on the dynamics of the protein. These proteins were studied by a variety of methods, including one and two-dimensional proton nuclear magnetic resonance spectroscopy, exchange kinetics for the distal and proximal histidine NH protons as a function of pH in the met cyano forms, flash photolysis of the CO forms, and ligand replacement kinetics. The electronic absorption and proton nuclear magnetic resonance spectra of the CO forms of these proteins are virtually identical, indicating that the structure of the heme pocket is unaltered by these mutations. There are, however, substantial changes in the dynamics of both CO binding and proton exchange for the mutant K45R, whereas the mutant D60E exhibits behavior indistinguishable from the reference human myoglobin. K45R has a faster CO bimolecular recombination rate and slower CO off-rate relative to the reference. The kinetics for CO binding are independent of pH (6.5 to 10) as well as ionic strength (0 to 1 M-NaCl). The exchange rate for the distal histidine NH is substantially lower for K45R than the reference, whereas the proximal histidine NH exchange rate is unaltered. The exchange behavior of the human proteins is similar to that reported for a comparison of the exchange rates for myoglobins having lysine at position 45 with sperm whale myoglobin, which has arginine at this position. This indicates that the differences in exchange rates reflects largely the Lys----Arg substitution. The lack of a simple correlation for the CO kinetics with this substitution means that these are sensitive to other factors as well. Specific kinetic models, whereby substitution of arginine for lysine at position 45 can affect ligand binding dynamics, are outlined. These experiments demonstrate that a relatively conservative change of a surface residue can substantially perturb ligand and proton exchange dynamics in a manner that is not readily predicted from the static structures.

Inelastic neutron scattering analysis of picosecond internal protein dynamics. Comparison of harmonic theory with experiment

The experimental inelastic neutron scattering spectrum of a protein, the bovine pancreatic trypsin inhibitor (BPTI), in a powder sample is presented together with the generalized density of states, G(omega), as a function of the frequency, omega, derived from the scattering data. The experimental results are compared with calculations from two different normal mode analyses of BPTI. One of these, based on an improved model, gives a calculated spectrum and density of states in general agreement with those obtained experimentally; the other, based on an earlier model, shows considerable disagreement. The important improvements in the newer normal mode analysis are the explicit treatment of all atoms (non-polar as well as polar hydrogens are included) and a modified truncation scheme for the long-range electrostatic interactions. The fact that the inelastic neutron scattering measurements can distinguish between the two theoretical models makes clear their utility for the analysis of protein dynamics.

Effects of thermal denaturation on the longitudinal relaxation time (T1) of water protons in protein solutions: study of the factors determining the T1 of water protons

The factors determining the longitudinal relaxation time (T1) of water protons in protein solutions were investigated by analyzing the effects of thermal denaturation on the T1 of the water protons. We treated the water protons and the protein protons "on a protein surface" as a dipole-dipole coupled two-spin system where relative translational diffusion is the dominant mechanism, and measured the change in the time development of the nuclear Overhauser effect (NOE) factors of the water protons. The T1 of the water protons was shortened markedly when the proteins were thermally denatured. Our analysis indicates that this relaxation enhancement is due to an increase in the value of the translational correlation time as well as the fraction of hydration water molecules, though the influence of "proton exchange" between the water protons and the labile protein protons cannot be completely neglected.

1H-NMR sequential assignments and cation-binding studies of spinach plastocyanin

The essentially complete assignment of the 1H-NMR spectrum of the Cu(i) form of spinach plastocyanin has been achieved using two-dimensional NMR techniques and sequence-specific resonance assignment procedures. A variety of pH and temperature conditions was utilised to overcome the problems of resonance overlap in the spectrum, degeneracy of C alpha H and solvent H2O chemical shifts, and cross-saturation of labile NH resonances. A qualitative analysis of the long-range nuclear Overhauser effects observed indicates that the backbone fold of spinach plastocyanin is very similar to that of poplar plastocyanin, whose structure has been solved by X-ray crystallography and differs in 22 of its 99 amino acid residues. The assignments provide a basis for further investigations into the structural and ion- and protein-binding properties of plastocyanin in solution.

Conformations of intermediates in the folding of the pancreatic trypsin inhibitor

Intermediates in the folding pathway of the bovine pancreatic trypsin inhibitor (PTI) have been examined by 1H nuclear magnetic resonance (n.m.r.). The intermediates were trapped during the reoxidation and consequent refolding of reduced PTI by alkylating free thiols; each intermediate contained different disulphide linkages. The n.m.r. spectra reveal that conformational features of the native protein are present in the intermediate containing just one of the three normal disulphide linkages (30-51). As additional normal disulphide bonds are formed, the conformation becomes more similar to that of the native protein. Introduction of additional but incorrect disulphide bonds does not lead to an increase in observable globular structure. A description of the folding process in terms of the conformations of the different intermediates is proposed. The significance of these results for the general mechanism of protein folding is outlined.

Nuclear magnetic resonance study of the isotope exchange of the proximal histidyl ring labile protons in hemoglobin A. The exchange rates and mechanisms of individual subunits in deoxy and oxy-hemoglobin

Proton nuclear magnetic resonance spectroscopy has been used to investigate the rates and mechanism of exchange with deuterium of the proximal histidyl imidazole labile ring proton in deoxy and oxy-hemoglobin A. The resolved signals for the two subunits indicate dynamic heterogeneity, with the exchange rate always faster in the alpha than the beta subunits, suggesting a lower dynamic stability for the alpha subunit. The activation energy for the exchange in both subunits (approximately 25 kcal; 1 cal = 4.184 J) indicates that exchange proceeds via an intermediate far from denaturation or global unfolding. The pH profiles for both hemoglobin states reflect the EX2 mechanism for both subunits. While the base catalysis expected for an iron-bound imidazole is observed in all cases, there are important differences in both rates and mechanisms between the subunits. In deoxy-hemoglobin, both base-catalyzed and water-assisted exchange contribute to the alpha subunit, but only the former to the beta subunit. For oxy-hemoglobin, the base-catalysis is retained for both subunits, but the slope is considerably less for the alpha relative to the beta subunit. Thus the two subunits in the two states of hemoglobin differ both in mechanisms and in the inherent dynamic stability reflected in any one mechanism. The relationships of the proximal histidyl ring NH exchange rates to previously characterized subsets of allosterically responsive protons in hemoglobin A is briefly discussed.

Protein volumes and hydration effects. The calculations of partial specific volumes, neutron scattering matchpoints and 280-nm absorption coefficients for proteins and glycoproteins from amino acid sequences

Amino acid sequences, carbohydrate compositions and residue volumes are used to compare critically calculations of partial specific volumes v, neutron scattering matchpoints and 280-nm absorption coefficients with experimental v values for proteins and glycoproteins. The v values that are obtained from amino acid densitometry underestimate experimental v values by 0.01-0.02 ml/g while the v values from crystallographic volumes overestimate the experimental v values by 0.04-0.05 ml/g. An intermediate consensus volume set of amino-acid-residue volumes is proposed in order to predict experimental v values using sequence information. The method is extended to carbohydrates and glycoproteins. Neutron scattering matchpoints can be calculated from crystallographic residue volumes on the basis of the non-exchange of 10% of the main-chain NH protons. Crystallographic results on protein-bound water are used to account for the experimental values of v and matchpoints. Finally, 280-nm absorption coefficients, A1%, 1 cm 280, of 5-27 are found to be well predicted by the Wetlaufer procedure based on the totals of Trp, Tyr and Cys residues. Average errors are +/- 0.7, and the experimental A(1%,1cm)280 values can be larger than the predicted values by 3%.

Protein dynamics investigated by neutron diffraction

The most correct model of the molecular structure of a protein molecule is one which describes it as having a number of variable conformational states. These states differ in degree over a large spectrum of structural variation ranging from individual atomic vibrational motion to significant tertiary denaturation. The application of the neutron diffraction techniques discussed above dealt with two classes of conformational fluctuation, "protein breathing" and "regional melting." The utility of the neutron technique stems from the ability to locate hydrogen atoms and to discriminate between hydrogens and deuteriums. This latter attribute allows for performing H/D exchange experiments by identifying individual sites of exchange. With this information it has been possible to discern which regions of the protein molecule undergo regional melting. Protein breathing was explored by analyzing the rotational properties of side chain methyl groups. Such information clearly suggested that most of these groups reside in "staggered" (low energy) conformation in the time averaged structure and are not greatly affected by local structural packing. Together these two classes of conformational fluctuation span nearly the full range of all motions that might play a role in biological activity. Information about such motions can be obtained from other types of physicochemical methods, but in most cases the interpretation of the data is considerably less definitive than that which can be obtained from a neutron analysis.

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.

Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor

The importance of including different energy contributions in calculations of electrostatic energies in proteins is examined by calculating the intrinsic pKa values of the acidic groups of bovine pancreatic trypsin inhibitor. It appears that such calculations provide a powerful and revealing test; the relevant solvation energies of the ionized acids are of the order of -70 kcal/mol (1 cal = 4.184 J), and microscopic calculations that do not attempt to simulate the complete protein dielectric effect (including the surrounding solvent) can underestimate the solvation energy by as much as 50 kcal/mol. Reproducing correctly, by the same set of parameters, the solvation energies of ionized acids in different sites of a protein cannot be accomplished by including only part of the key energy contributions. The problems associated with macroscopic calculations are also considered and illustrated by the specific case of bovine pancreatic trypsin inhibitor. A promising approach is shown to be provided by a refinement of the previously developed Protein Dipoles Langevin Dipoles model. This model seems to represent consistently the microscopic dielectric of the protein and the surrounding water molecules. The model overcomes the problems associated with the macroscopic models (by treating explicitly the solvent molecules) and avoids the convergence problems associated with all-atom solvent models (by treating the average solvent polarization rather than averaging the actual polarization energy). This paper describes in detail the actual implementation of the model and examines its performance in evaluating intrinsic pKa values. Preliminary microscopic considerations of charge-charge interactions are presented.

Mechanism of surface peptide proton exchange in bovine pancreatic trypsin inhibitor. Salt effects and O-protonation

The acid-catalyzed hydrogen exchange rate constants kH, and the base-catalyzed rate constants kOH, have been determined (in the preceding paper) for the 25 most rapidly exchanging NH groups of bovine pancreatic trypsin inhibitor. Most of these NH groups are at the protein-solvent interface. The correlation of kH, but not kOH, with the static accessibility and hydrogen bonding of the peptide carbonyl O atom indicates that the mechanism of acid catalysis in proteins involves O-protonation. Agreement between the ionic strength dependence observed for kH and kOH and the ionic strength dependence calculated for an O-protonation mechanism supports this conclusion. N-protonation for acid catalysis, as well as N-deprotonation for base catalysis, have traditionally been assumed in the mechanism of the chemical step in peptide amide proton exchange. A preference for the alternative O-protonation mechanism has far-reaching implications in the interpretation of protein hydrogen exchange kinetics. With an O-protonation mechanism, acid-catalyzed rates of surface NH groups are primarily a function of the average solvent accessibility of the carbonyl O atoms in the dynamic solution structure, while base-catalyzed rates of surface NH groups measure solvent accessibility of the peptide N. The relative dynamic accessibilities of peptide O atoms, as measured by relative values of kH (corrected for electrostatic effects), correlate with O static accessibilities in the crystal structure. A lower correlation of static accessibility of N atoms with kOH is observed for surface NH groups in peptide groups in which the carbonyl O is not hydrogen bonded. For some surface NH groups, the observed pH of minimum rate, pHmin, deviates widely from the pHmin of model compounds. This is explained as the combined result of electrostatic effects and of the differences in accessibility of the carbonyl O and N atoms that result in a change in the relative values of kH and kOH as compared to those of model peptides. A mechanism whereby exchange of interior sites is catalyzed by interactions of catalysis ions with protein surface atoms via charge transfer is suggested.

Calculations of electrostatic interactions in biological systems and in solutions

Correlating the structure and action of biological molecules requires knowledge of the corresponding relation between structure and energy. Probably the most important factors in such a structure– energy correlation are associated with electrostatic interactions. Thus the key requirement for quantative understanding of the action of biological molecules is the ability to correlate electrostatic interactions with structural information. To appreciate this point it is useful to compare the electrostatic energy of a charged amino acid in a polar solvent to the corresponding van der Waals energy. The electrostatic free energy, ΔGel, can be approximated (as will be shown in Section II) by the Born formula (ΔGel= –(166Q2/ā) (I – I/E)). Where ΔGelis given in kcal/mol,Qis the charge of the given group, in units of electron charge,āis the effective radius of the group, andEis the dielectric constant of the solvent. With an effective radius of charged amino acids of approximately 2 Å, Born's formula gives about – 80 kcal/mol for their energy in polar solvents whereEis larger than 10. This energy is two orders of magnitude larger than the van der Waals interaction of such groups and their surroundings.

Carbohydrate-containing derivatives of the trypsin-kallikrein inhibitor aprotinin from bovine organs. I. Modification with lactose, characterization and behaviour of the preparation in vivo

The trypsin-kallikrein inhibitor aprotinin was modified with lactose. The influence of reactant concentrations, temperature, reaction time and sodium borohydride on the carbohydrate residue content and the inhibiting activity of glycated aprotinin were studied. Glycation of aprotinin neither shifts the pH optimum of the inhibitor-trypsin association reaction nor does it alter the apparent dissociation constant Ki of the complex measured at pH optimum. Glycation by lactose stabilizes aprotinin against denaturation by increased temperature. The distribution of native and modified aprotinin in rat organs after endocardiac injection was studied. Fixation of glycated aprotinin increases 2.5- to 3-fold in liver and decreases 2-fold in kidneys during the observation time (5 min-2 h) compared to native aprotinin.

NMR studies of the variant-3 neurotoxin from Centruroides sculpturatus Ewing

We report a preliminary high-resolution proton nuclear magnetic resonance characterization of the variant-3 toxin from the scorpion Centruroides sculpturatus Ewing (range Southwestern USA). This toxin assumes a well defined folded conformation in aqueous solutions at room temperature and undergoes reversible thermal denaturation. A number of amide hydrogens exhibit exchange life times varying from several minutes to several hours. A few tentative assignments of the low field aromatic CH resonances has been made on the basis of 2D-COSY and NOE experiments. The upfield shifts exhibited by Trp-47 suggest a unique microenvironment for this residue. The NMR data suggest that there is some degree of correlation between the solution structure of the variant-3 toxin and its crystallographic structure. Our studies provide a basis for a detailed elucidation of the structure-function relationships of these interesting scorpion toxins which bind to the sodium channels of excitable membranes and delay sodium current inactivation.

Basic pancreatic trypsin inhibitor has unusual thermodynamic stability parameters

The stability parameters delta Gst, delta Hst and delta Sst of native basic pancreatic trypsin inhibitor (BPTI) have been characterized by microcalorimetric unfolding studies in various buffer solutions, at different pH values and in the presence of guanidine hydrochloride. The unfolding enthalpy of BPTI, in contrast ot other globular proteins, exhibits a very small dependence on temperature, which results in a characteristic different temperature dependence of the Gibbs energy of stabilization. BPTI has a very high specific Gibbs energy of stabilization, which renders the slow exchange rates of amide protons understandable. Comparison of the unfolding entropy of BPTI at 110 degrees C with corresponding values of other proteins, revealed that the delta S values of BPTI are lower by 2.9 J/(K X residue). This lower value of the unfolding entropy is in good agreement with predictions of a theoretical study by Poland & Scheraga (1965) where the influence of crosslinks on the configurational entropy has been studied. Additionally, we were able to calculate an interaction enthalpy per site of -5.6 kJ/mol based on the measurements of unfolding of BPTI in 6 M-guanidine hydrochloride.

Hydrogen exchange and structural dynamics of proteins and nucleic acids

Though the structures presented in crystallographic models of macromolecules appear to possess rock-like solidity, real proteins and nucleic acids are not particularly rigid. Most structural work to date has centred upon the native state of macromolecules, the most probable macromolecular form. But the native state of a molecule is merely its most abundant form, certainly not its only form. Thermodynamics requires that all other possible structural forms, however improbable, must also exist, albeit with representation corresponding to the factor exp( —Gi/RT) for each state of free energyGi(see Moelwyn-Hughes, 1961), and one appreciates that each molecule within a population of molecules will in time explore the vast ensemble ofpossiblestructural states.