Library

Notes: The equilibrium between the cis and trans forms of X-Pro peptide bonds can readily be measured in the 13C nmr spectra. In the present paper we investigate how observation of this equilibrium could be used as an nmr probe for conformational studies of flexible polypeptide chains. The experiments include studies by 13C nmr of a series of linear oligopeptides containing different X-L-Pro peptide bonds, with X = Gly, L-Ala, L-Leu, L-Phe, D-Ala, D-Leu, and D-Phe. Overall the study confirms that X-Pro peptide bonds can generally be useful as 13C nmr probes reporting the formation of nonrandom conformation in flexible polypeptide chains. It was found that the cis-trans equilibrium of X-Pro is greatly affected by the side chain of X and the configuration of the α-carbon atom of X. On the basis of these observations some general rules are suggested for a practical applications of the X-Pro nmr probes in conformational studies of polypeptide chains.

Notes: The molecular conformations of the linear oligopeptides H-(L-Ala)n-L-Pro-OH, with n = 1,2 and 3, have been investigated. 13C nmr observation of the equilibrium between the cis and trans forms of the Ala-Pro peptide bond indicated the occurrence of nonrandom conformations in solutions of these flexible peptides. The formation of the nonrandom species containing the cis form of the Ala-Pro bond was found to depend on the deprotonation of the carboxylic acid group of proline, the solvent, and the ionic strength in aqueous solution. The influence of intramolecular hydrogen bonding on the relative conformational energies of the species containing the cis and trans Ala-Pro peptide bond was studied by comparison of the peptides H-(Ala)n-Pro-OH with analogous molecules where hydrogen bond formation was excluded by the covalent structure. In earlier work a hydrogen bond between the protonated terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue had been suggested to stabilize conformations including trans proline. For the systems described here this hypothesis can be ruled out, since the cis:trans ratio is identical for molecules with methyl ester protected and free protonated terminal carboxylic acid groups of proline. Direct evidence for hydrogen bond formation between the deprotonated terminal carboxylic acid group and the amide proton of the penultimate amino acid residue in the molecular species containing cis proline was obtained from 1H nmr studies. However, the cis:trans ratio of the Ala-Pro bond was not affected by N-methylation of the penultimate amino acid residue, which prevents formation of this hydrogen bond. Overall the experimental observations lead to the conclusion that the relative energies of the peptide conformations including cis or trans proline are mainly determined by intramolecular electrostatic interactions, whereas in the molecules considered, intramolecular hydrogen bonding is a consequence of specific peptide backbone conformations rather than a cause for the occurrence of energetically favored species. Independent support for this conclusion was obtained from model consideration which indicated that electrostatic interactions between the terminal carboxylic acid group and the carbonyl oxygen of the penultimate amino acid residue could indeed account for the observed relative conformational energies of the species containing cis and trans proline, respectively.

Notes: The 1H-nmr chemical shifts and the spin-spin coupling constants of the common amino acid residues were measured in solutions of the linear tetrapeptides H-Gly-Gly-X-L-Ala-OH in D2O and H2O, the influence of X on the nmr parameters of the neighboring residues Gly 2 and Ala 4 was investigated. The titration parameters for the side chains of Asp, Glu, Lys, Tyr, and His were determined. The pKa values obtained in D2O, with the use of pH-meter readings with a combination glass electrode uncorrected for istope effects, were 0.06 pH units higher in the acidic range and 0.10 pH units higher in the basic range than the corresponding pKa values in H2O. This suggests that the present data are suitable “random-coil” 1H-nmr parameters for conformational studies of polypeptide chains in D2O and H2O solutions.

Notes: This paper shows that backbone amide proton titration shifts in polypeptide chains are a very sensitive manifestation of intramolecular hydrogen bonding between carboxylate groups and backbone amide protons. The population of specific hydrogen-bonded structures in the ensemble of species that constitutes the conformation of a flexible nonglobular linear peptide can be determined from the extent of the titration shifts. As an illustration, an investigation of the molecular conformation of the linear peptide H-Gly-Gly-L-Glu-L-Ala-OH is described. The proposed use of amide proton titration shifts for investigating polypeptide conformation is based on 360-MHz 1H-nmr studies of selected linear oligopeptides in H2O solutions. It was found that only a very limited number of amide protons in a polypeptide chain show sizable intrinsic intration shifts arising from through-bond interactions with ionizable groups. These are the amide proton of the C-terminal amino acid residue, the amide protons of Asp and the residues following Asp, and possibly the amide proton of the residue next to the N-terminus. Since the intrinsic titration shifts are upfield, the downfield titration shifts arising from conformation-dependent through-space interactions, in particular hydrogen bonding between the amide protons and carboxylate groups, can readily be identified.

Notes: Empirical conformational energy calculations with the use of ECEPP energy functions have been carried out for linear dipeptides H-X-L-Pro-OH, with X = Gly, L-Ala, D-Ala, L-Leu, D-Leu, L-Phe, and D-Phe, in different states of protonation of the end groups. The results of these calculations are compared with the previously reported experimental equilibrium populations for the cis and trans isomers of the X-Pro bond in the different species. For all the protonation states of the seven dipeptides, the calculated nonbonded interactions and the conformational entropy term lead to a preference of the trans forms over the cis isomers by at least 1 kcal/mol. The electrostatic interactions stabilize the cis conformations in all species except the cationic forms of the D,L-peptides, and it could further be shown that only the carbonyl group of X and the two end groups contribute significantly to the total electrostatic energy. One of the principal results of the experimental studies, i.e., the occurrence of 5-15% cis-proline in all the peptides with an uncharged C-terminus, was corroborated by our investigation of the cationic species. A detailed assessment of the electrostatic contribution to the total energy of the different conformations of H-Gly-L-Pro-OH indicates that the standard ECEPP parameters tend to overestimate the electrostatic interactions in aqueous solutions of the X-Pro dipeptides.

Notes: 1H-Nmr was used to measure the rate of cis-trans interconversion of X-Pro bonds in linear and cyclic oligopeptides. k(cis → trans) = 2.5 × 10-3 s-1 at 25°C was found for the zwitterionic form of H-Ala-Pro-OH, in good agreement with earlier measurements. Replacement of Ala by Phe, Tyr, or Trp resulted in a 10-fold slower interconversion rate, whereas after substitution of Ala by His or Glu, the rate decreased only slightly. Independent of the residues X, the interconversion rate was increased by a factor of ca. 20 when the peptide chain was elongated by addition of Ala to the C-terminal Pro. An additional increase by a factor of 6 was observed when going from the protected linear peptide CF3CO-Gly-Gly-Pro-Ala-OCH3 to the closely related cyclic compound c[-Gly-Gly-Pro-Gly-Ala-]. These data are evaluated with regard to their possible use in future studies on the role of X-Pro cis-trans isomerization in the kinetics of protein folding.

Notes: The program FANTOM (fast Newton-Raphson torsion angle energy minimizer) performs minimizations of the ECEPP/2 energy function for proteins with the Newton-Raphson method. It is implemented for use with conventional computer hardware. The torsion angles are chosen as independent variables. The first and second derivatives are calculated with a previously described rapid algorithm. For the matrix inversion a modified Cholesky factorization is used. A line search adjusts the step length and nonbonded interactions can be calculated with a cutoff. The following tests of the program are described: All local minima of the ECEPP/2 energy function for the amino acids glycine and alanine were determined. An exhaustive search by more than 16,000 independent energy minimizations was used to identify low-energy structures of Met-enkephalin, which were then compared with previously published structures of this pentapeptide. To investigate the use of FANTOM with disulfide bonds, it was applied with conotoxin. As an illustration of the intended primary use of the program, an energy refinement of the structure of the basic pancreatic trypsin inhibitor determined by nmr spectroscopy in solution is described.

Notes: The high-resolution 1H-nmr study of the ferrichrome cyclohexapeptides, in d6-DMSO solutions, has been extended to the amide NH spectral region. A total of ten diamagnetic analogues of ferrichrome that differ in the coordinated metal ion (Al3+, Ga3+ or Co3+), the primary structure, the nature of the bidentate hydroxamate moiety, or the isotope compositions (14N, 15N) have been investigated. The 3JαNH values reflect regiorous conformational isomorphism throughout the complete suite of analogues, quite independent of the residue occupancy of each site. Totally resolved amide multiplets have been obtained in most cases and the four-line (doublet of doublets) appearances of glycyl NH resonances has been observed for the first times; these data enabled stereospecific assignment and accurate analysis of the NH-CαH proton spin systems. The high resolution was made possible by the use of a suitable spectral deconvolution shceme at 360 MHz. The fine structure, extraordinarily well displayed in the 15N-peptide spectra, provides a series of parameter values whose consistency has been checked by computer simulation. Since the crystallographic structure for two of the ferric peptides is known to 0.002-Å resolution, a 3J vs θ correspondence could be confidently established. A Karplus curve was derived from the combined x-ray and nmr data: \documentclass{article}\pagestyle{empty}\begin{document}$$ ^3 J_{\alpha {\rm NH}} = 5.4\cos ^2 \theta - 1.3\cos \theta + 2.2{\rm Hz} $$\end{document} It is suggested that seriously nonplanar amides can exhibit 3JαNH values higher than predicted by the ferrichrome curve.

Notes: Why Pentose- and Not Hexose-Nucleid Acids? Part IV . ‘Homo-DNA’: 1H-, 13C-, 31P-, and 15N-NMR-Spectroscopic Investigation of ddGlc(A-A-A-A-A-T-T-T-T-T) in Aqueous SolutionFrom a comprehensive NMR structure analysis, it is concluded that the ‘homo-DNA’ oligonucleotide ddGlc(A-A-A-A-A-T-T-T-T-T) in 3 mM D2O solution (100 mM NaCl, 50 mM phosphate buffer, pH 7.0, T = 50°) forms a duplex of C2-symmetry, with its self-complementary oligonucleotide strands in antiparallel orientation. The 2′,3′-dideoxy-β-D-glucopyranosyl rings are in their most stable chair conformation, with all three substituents equatorial and with the adenine as well as the thymine bases in the anti-conformation. The base pairing is of the Watson-Crick type; this pairing mode (as opposed to the reverse-Hoogsteen mode) was deduced from the observation of inter strand NOEs between the adenine protons H—C(2) and the pyranose protons Hα-C(2′) of the sequentially succeeding thymidine nucleotides of the opposite strand, a correlation which discriminates between the Watson-Crick and the reverse-Hoogsteen pairing mode. The NOEs of the NH protons with either the adenine protons H—C(2) or H—C(8), that are normally used to identify the pairing mode in DNA duplexes, cannot be observed here, because the NH signals are very broad. This line broadening is primarily due to the fact that the exchange of the imino protons with the solvent is faster than for corresponding DNA duplexes.Computer-assisted modeling of the [ddGlc(A5-T5)]2 duplex with the program CONFOR [23], using the linear (idealized) homo-DNA single-strand conformation (α = -60°, β = 180°, γ = 60°, δ = 60°, ∊ = 180°, ζ = -60°, see [1] [3]) as the starting structure, resulted in two duplex models A and B (see Figs. 27-32, Scheme 9, and Table 4) which both contain quasi-linear double strands with the base-pairing axis inclined relative to the strand axes by ca. 60° and 45°, respectively, and with base-pair stacking distances of ca. 4.5 Å. While neither of the two models, taken separately, can satisfy all of the NMR constraints, the NMR data can be rationalized by the assumption that the observed duplex structure represents a dynamic equilibrium among conformers which relate to models A and B as their limiting structure. The required rapid equilibrium appears feasible, since the models A and B are interconvertible by two complementary 120° counter rotations around the α-axis and the γ-axis, respectively, of the phosphodiester backbone. The models A and B correspond to the two types of linear (idealized) single-strand backbone conformation derived previously by qualitative conformational analysis without and with allowance for gauche-trans-phosphodiester conformations, respectively [1] [3]. Refinement of the models A and B with the use of the program AMBER [27] by energy minimization in a water bath and molecular-dynamics simulations (2 ps, 300° K) resulted in two dynamic structures (Figs. 33 and 34, Table 4). These have roughly the same energy, closely resemble the starting structures A and B, and satisfy - as an ensemble - all of the NMR constraints without violating any van der Waals distances by more than 0.2 Å. Extensive fluctuations in base-pair distance and deviations from base-pair coplanarity, as well as the presence of water molecules in the cavities between some of the base pairs, were observed in both dynamic structures A and B, which, on the other hand, did not mutually interconvert within the short simulation time period used. These model properties, together with the conjectured equilibrium between the two structure types A and B, lead to the hypothesis of a homo-DNA duplex containing a ‘partially molten’ pairing core. This proposal could qualitatively account for a high rate of the NH exchange, as well as for part of the previously established [3] deficits in both enthalpic stabilization and entropic destabilization of homo-DNA duplexes relative to corresponding DNA duplexes. The phenomenon of the higher overall stability of homo-DNA duplexes vs. DNA duplexes (e.g, [ddGlc(A5-T5)]2, Tm = 59° vs. [d(A5-T5)]2, Tm = 33°, both at c ≈ 50 μM [3]) can then be seen as the result not only of a higher degree of conformational preorganization of the homo-DNA single strand toward the conformation of the duplex backbone [1] [3], but also of the entropic benefit of greater disorder in the central pairing zone of the homo-DNA duplex. This view of the structure of a homo-DNA duplex relates its characteristic properties to a central structural feature: the average base-pair distance in the models of homo-DNA is too large for regular base stacking (ca. 4.5 Å vs. ca. 3.5 Å in DNA). This difference in the distances between adjacent base pairs is a direct consequence of the quasi-linearity of the homo-DNA double strand as opposed to the right-handed twist of the helical DNA duplexes [1] [3], which is directly related to the specific conformational properties of pyranose rings as opposed to furanose rings [1]. Thus, the structural hypothesis derived from the NMR analysis of [ddGlc(A5-T5)]2 relates the conformational differences between homo-DNA and DNA directly to the sugar ring size, which is the essential constitutional difference between the two types of structure.The English footnotes to Figs. 1-34, Schemes 1-9, and Tables 1-4 provide an extension of this summary.