The 13 C a chemical shifts for 16,299 residues from 213 conformations of four proteins (experimentally determined by X-ray crystallography and Nuclear Magnetic Resonance methods) were computed by using a combination of approaches that includes, but is not limited to, the use of density functional theory. Initially, a validation test of this methodology was carried out by a detailed examination of the correlation between computed and observed 13 C a chemical shifts of 10,564 (of the 16,299) residues from 139 conformations of the human protein ubiquitin. The results of this validation test on ubiquitin show agreement with conclusions derived from computation of the chemical shifts at the ab initio Hartree-Fock level. Further, application of this methodology to 5,735 residues from 74 conformations of the three remaining proteins that differ in their number of amino acid residues, sequence and three-dimensional structure, together with a new scoring function, namely the conformationally averaged root-mean-square-deviation, enables us to: (a) offer a criterion for an accurate assessment of the quality of NMRderived protein conformations; (b) examine whether X-ray or NMR-solved structures are better representations of the observed 13 C a chemical shifts in solution; (c) provide evidence indicating that the proposed methodology is more accurate than automated predictors for validation of protein structures; (d) shed light as to whether the agreement between computed and observed 13 C a chemical shifts is influenced by the identity of an amino acid residue or its location in the sequence; and (e) provide evidence confirming the presence of dynamics for proteins in solution, and hence showing that an ensemble of conformations is a better representation of the structure in solution than any single conformation.

The dependence of the (13)C chemical shift on side-chain orientation was investigated at the density functional level for a two-strand antiparallel beta-sheet model peptide represented by the amino acid sequence Ac-(Ala)(3)-X-(Ala)(12)-NH(2) where X represents any of the 17 naturally occurring amino acids, i.e., not including alanine, glycine and proline. The dihedral angles adopted for the backbone were taken from, and fixed at, observed experimental values of an antiparallel beta-sheet. We carried out a cluster analysis of the ensembles of conformations generated by considering the side-chain dihedral angles for each residue X as variables, and use them to compute the (13)C chemical shifts at the density functional theory level. It is shown that the adoption of the locally-dense basis set approach for the quantum chemical calculations enabled us to reduce the length of the chemical-shift calculations while maintaining good accuracy of the results. For the 17 naturally occurring amino acids in an antiparallel beta-sheet, there is (i) good agreement between computed and observed (13)C(alpha) and (13)C(beta) chemical shifts, with correlation coefficients of 0.95 and 0.99, respectively; (ii) significant variability of the computed (13)C(alpha) and (13)C(beta) chemical shifts as a function of chi(1) for all amino acid residues except Ser; and (iii) a smaller, although significant, dependence of the computed (13)C(alpha) chemical shifts on chi(xi) (with xi > or = 2) compared to chi(1) for eleven out of seventeen residues. Our results suggest that predicted (13)C(alpha) and (13)C(beta) chemical shifts, based only on backbone (phi,psi) dihedral angles from high-resolution X-ray structure data or from NMR-derived models, may differ significantly from those observed in solution if the dihedral-angle preferences for the side chains are not taken into account.

A theoretical study to identify the conformational preferences of lysine-based oligopeptides has been carried out. The solvation free energy and free energy of ionization of the oligopeptides have been calculated by using a fast multigrid boundary element method that considers the coupling between the conformation of the molecule and the ionization equilibria explicitly, at a given pH value. It has been found experimentally that isolated alanine and lysine residues have somewhat small intrinsic helix-forming tendencies; however, results from these simulations indicate that conformations containing right-handed alpha-helical turns are energetically favorable at low values of pH for lysine-based oligopeptides. Also, unusual patterns of interactions among lysine side chains with large hydrophobic contacts and close proximity (5-6 A) between charged NH3+ groups are observed. Similar arrangements of charged groups have been seen for lysine and arginine residues in experimentally determined structures of proteins available from the Protein Data Bank. The lowest-free-energy conformation of the sequence Ac-(LYS)6-NMe from these simulations showed large pKalpha shifts for some of the NH3+ groups of the lysine residues. Such large effects are not observed in the lowest-energy conformations of oligopeptide sequences with two, three, or four lysine residues. Calculations on the sequence Ac-LYS-(ALA)4-LYS-NMe also reveal low-energy alpha-helical conformations with interactions of one of the LYS side chains with the helix backbone in an arrangement quite similar to the one described recently by (Proc. Natl. Acad. Sci. U.S.A. 93:4025-4029). The results of this study provide a sound basis with which to discuss the nature of the interactions, such as hydrophobicity, charge-charge interaction, and solvent polarization effects, that stabilize right-handed alpha-helical conformations.