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.

Chitlac is a biocompatible modified polysaccharide composed of a chitosan backbone to which lactitol moieties have been chemically inserted via a reductive N-alkylation reaction with lactose. The physical-chemical and biological properties of Chitlac that have been already reported in the literature suggest a high accessibility of terminal galactose in the lactitol side chain. This finding may account for its biocompatibility which makes it extremely interesting for the production of biomaterials. The average structure and the dynamics of the side chains of Chitlac have been studied by means of NMR (nuclear Overhauser effect and nuclear relaxation) and molecular dynamics to ascertain this hypothesis. A complete assignment of the (1)H and (13)C NMR signals of the modified polysaccharide has been accomplished together with the determination of the apparent pKa values of the primary and secondary amines (6.69 and 5.87, respectively). NMR and MD indicated a high mobility of Chitlac side chains with comparable average internuclear distances between the two techniques. It was found that the highly flexible lactitol side chain in Chitlac can adopt two distinct conformations differing in the orientation with respect to the polysaccharide chain: a folded conformation, with the galactose ring parallel to the main chain, and an extended conformation, where the lactitol points away from the chitosan backbone. In both cases, the side chain resulted to be highly hydrated and fully immersed in the solvent.

No abstract

No abstract

Cysteines possess a unique property among the 20 naturally occurring amino acids: it can be present in proteins in either the reduced or oxidized form, and can regulate the activity of some proteins. Consequently, to augment our previous treatment of the other types of residues, the 13Cα and 13Cβ chemical shifts of 837 cysteines in disulfide-bonded cystine from a set of seven non-redundant proteins, determined by X-ray crystallography and NMR spectroscopy, were computed at the DFT level of theory. Our results indicate that the errors between observed and computed 13Cα chemical shifts of such oxidized cysteines can be attributed to several effects such as: (a) the quality of the NMR-determined models, as evaluated by the conformational-average (ca) rmsd value; (b) the existence of high B-factor or crystal-packing effects for the X-ray-determined structures; (c) the dynamics of the disulfide bonds in solution; and (d) the differences in the experimental conditions under which the observed 13Cα chemical shifts and the protein models were determined by either X-ray crystallography or NMR-spectroscopy. These quantum-chemical-based calculations indicate the existence of two, almost non-overlapped, basins for the oxidized and reduced –SH 13Cβ, but not for the 13Cα, chemical shifts, in good agreement with the observation of 375 13Cα and 337 13Cβ resonances from 132 proteins by Sharma and Rajarathnam (2000). Overall, our results indicate that explicit consideration of the disulfide bonds is a necessary condition for an accurate prediction of 13Cα and 13Cβ chemical shifts of cysteines in cystines.

No abstract