Analysis of proton chemical shifts in regular secondary structure of proteins

Using NMR Chemical Shifts to Determine Residue-Specific Secondary Structure Populations for Intrinsically Disordered Proteins

Protein disorder is a pervasive phenomenon in biology and a natural consequence of polymer evolution that facilitates cell signaling by organizing sites for posttranslational modifications and protein-protein interactions into arrays of short linear motifs that can be rearranged by RNA splicing. Disordered proteins are missing the long-range nonpolar interactions that form tertiary structures, but they often contain regions with residual secondary structure that are stabilized by protein binding. NMR spectroscopy is uniquely suited to detect residual secondary structure in a disordered protein and it can provide atomic resolution data on the structure and dynamics of disordered protein interaction sites. Here we describe how backbone chemical shifts are used for assigning residual secondary structure in disordered proteins and discuss some of the tools available for estimating secondary structure populations with a focus on disordered proteins containing different levels of alpha helical secondary structure which are stabilized by protein binding.

Preparation of fibril nuclei of beta-amyloid peptides in reverse micelles

Protofibrils of beta-amyloid peptides formed by fibril nuclei incubated in reverse micelles.

Predicting the redox state and secondary structure of cysteine residues using multi-dimensional classification analysis of NMR chemical shifts

A tool for predicting the redox state and secondary structure of cysteine residues using multi-dimensional analyses of different combinations of nuclear magnetic resonance (NMR) chemical shifts has been developed. A data set of cysteine [Formula: see text], (13)C(α), (13)C(β), (1)H(α), (1)H(N), and (15)N(H) chemical shifts was created, classified according to redox state and secondary structure, using a library of 540 re-referenced BioMagResBank (BMRB) entries. Multi-dimensional analyses of three, four, five, and six chemical shifts were used to derive rules for predicting the structural states of cysteine residues. The results from 60 BMRB entries containing 122 cysteines showed that four-dimensional analysis of the C(α), C(β), H(α), and N(H) chemical shifts had the highest prediction accuracy of 100 and 95.9 % for the redox state and secondary structure, respectively. The prediction of secondary structure using 3D, 5D, and 6D analyses had the accuracy of ~90 %, suggesting that H(N) and [Formula: see text] chemical shifts may be noisy and made the discrimination worse. A web server (6DCSi) was established to enable users to submit NMR chemical shifts, either in BMRB or key-in formats, for prediction. 6DCSi displays predictions using sets of 3, 4, 5, and 6 chemical shifts, which shows their consistency and allows users to draw their own conclusions. This web-based tool can be used to rapidly obtain structural information regarding cysteine residues directly from experimental NMR data.

Structure of ring-shaped Aβ₄₂ oligomers determined by conformational selection

The oligomerization of amyloid beta (Aβ) peptides into soluble non-fibrillar species plays a critical role in the pathogenesis of Alzheimer’s disease. However, it has been challenging to characterize the tertiary and quaternary structures of Aβ peptides due to their disordered nature and high aggregation propensity. In this work, replica exchange molecular dynamics simulations were used to explore the conformational space of Aβ₄₂ monomer. Among the most populated transient states, we identified a particular conformation which was able to generate ring-shaped pentamers and hexamers, when docked onto itself. The structures of these aggregates were stable during microsecond all-atom MD simulations in explicit solvent. In addition to high resolution models of these oligomers, this study provides support for the conformational selection mechanism of Aβ peptide self-assembly.

CSI 3.0: a web server for identifying secondary and super-secondary structure in proteins using NMR chemical shifts

The Chemical Shift Index or CSI 3.0 (http://csi3.wishartlab.com) is a web server designed to accurately identify the location of secondary and super-secondary structures in protein chains using only nuclear magnetic resonance (NMR) backbone chemical shifts and their corresponding protein sequence data. Unlike earlier versions of CSI, which only identified three types of secondary structure (helix, β-strand and coil), CSI 3.0 now identifies total of 11 types of secondary and super-secondary structures, including helices, β-strands, coil regions, five common β-turns (type I, II, I′, II′ and VIII), β hairpins as well as interior and edge β-strands. CSI 3.0 accepts experimental NMR chemical shift data in multiple formats (NMR Star 2.1, NMR Star 3.1 and SHIFTY) and generates colorful CSI plots (bar graphs) and secondary/super-secondary structure assignments. The output can be readily used as constraints for structure determination and refinement or the images may be used for presentations and publications. CSI 3.0 uses a pipeline of several well-tested, previously published programs to identify the secondary and super-secondary structures in protein chains. Comparisons with secondary and super-secondary structure assignments made via standard coordinate analysis programs such as DSSP, STRIDE and VADAR on high-resolution protein structures solved by X-ray and NMR show >90% agreement between those made with CSI 3.0.

Insights into the molecular flexibility of θ-defensins by NMR relaxation analysis

θ-Defensins are mammalian cyclic peptides that have antimicrobial activity and show potential as stable scaffolds for peptide-based drug design. The cyclic cystine ladder structural motif of θ-defensins has been characterized using NMR spectroscopy and is important for their structure and stability. However, the effect of the pronounced elongated topology of θ-defensins on their molecular motion is not yet understood. Studies of molecular motion by NMR relaxation measurements have been facilitated by the recent development of a semirecombinant method for producing cyclic peptides that allows for isotopic labeling. Here we have undertaken a multifield (15)N NMR relaxation analysis of the anti-HIV θ-defensin, HTD-2, and interpreted the experimental data using various models of overall and internal molecular motion. We found that it was necessary to apply a model that includes internal motion to account for the variations in the experimental T1 and NOE data at different backbone amide sites in the peptide. Although an isotropic model with internal motion was the simplest model that provided a satisfactory fit with the experimental data, we cannot exclude the possibility that overall motion is anisotropic, especially considering the strikingly elongated topology of θ-defensins. The presence of flexible side chains, self-association, interactions with solvent, and internal motions are all potential contributors to the observed relaxation data. Internal motion consistent with the constraints imposed by the cyclic cystine ladder was observed in that the order parameters, S(2), show that residues in the turns are more flexible than those in the β-sheet. This study provides insights into the dynamics of θ-defensins and information that might be useful in their application as scaffolds in drug design.

Type I and II β-turns prediction using NMR chemical shifts

A method for predicting type I and II β-turns using nuclear magnetic resonance (NMR) chemical shifts is proposed. Isolated β-turn chemical-shift data were collected from 1,798 protein chains. One-dimensional statistical analyses on chemical-shift data of three classes β-turn (type I, II, and VIII) showed different distributions at four positions, (i) to (i + 3). Considering the central two residues of type I β-turns, the mean values of Cο, Cα, H(N), and N(H) chemical shifts were generally (i + 1) > (i + 2). The mean values of Cβ and Hα chemical shifts were (i + 1) < (i + 2). The distributions of the central two residues in type II and VIII β-turns were also distinguishable by trends of chemical shift values. Two-dimensional cluster analyses on chemical-shift data show positional distributions more clearly. Based on these propensities of chemical shift classified as a function of position, rules were derived using scoring matrices for four consecutive residues to predict type I and II β-turns. The proposed method achieves an overall prediction accuracy of 83.2 and 84.2% with the Matthews correlation coefficient values of 0.317 and 0.632 for type I and II β-turns, indicating that its higher accuracy for type II turn prediction. The results show that it is feasible to use NMR chemical shifts to predict the β-turn types in proteins. The proposed method can be incorporated into other chemical-shift based protein secondary structure prediction methods.

Design of supramolecular amino acids to template peptide folding

We describe in this manuscript a new design of supramolecular amino acids (SAAs) for directing peptide folding in organic environments. The incorporated supramolecular motif has a strong driving force to dimerize in a sequence- and orientation-specific manner. By introducing such SAAs into the primary sequence of peptides, the specific and directional dimerization of the supramolecular units should facilitate the folding of the peptides. Our approach may provide a general strategy to program secondary structures in organic media.

NightShift: NMR shift inference by general hybrid model training--a framework for NMR chemical shift prediction

Background

NMR chemical shift prediction plays an important role in various applications in computational biology. Among others, structure determination, structure optimization, and the scoring of docking results can profit from efficient and accurate chemical shift estimation from a three-dimensional model.

A variety of NMR chemical shift prediction approaches have been presented in the past, but nearly all of these rely on laborious manual data set preparation and the training itself is not automatized, making retraining the model, e.g., if new data is made available, or testing new models a time-consuming manual chore.

Results

In this work, we present the framework NightShift (NMR Shift Inference by General Hybrid Model Training), which enables automated data set generation as well as model training and evaluation of protein NMR chemical shift prediction.

In addition to this main result – the NightShift framework itself – we describe the resulting, automatically generated, data set and, as a proof-of-concept, a random forest model called Spinster that was built using the pipeline.

Conclusion

By demonstrating that the performance of the automatically generated predictors is at least en par with the state of the art, we conclude that automated data set and predictor generation is well-suited for the design of NMR chemical shift estimators.

The framework can be downloaded from https://bitbucket.org/akdehof/nightshift. It requires the open source Biochemical Algorithms Library (BALL), and is available under the conditions of the GNU Lesser General Public License (LGPL). We additionally offer a browser-based user interface to our NightShift instance employing the Galaxy framework via https://ballaxy.bioinf.uni-sb.de/.

Combined nuclear magnetic resonance spectroscopy and molecular dynamics study of growth hormone releasing hexapeptide GHRP-6 and a cyclic analogue

The Growth Hormone Releasing Hexapeptide, GHRP-6 was the first of a family of synthetic peptides that enhance the release of the Growth Hormone by the pituitary gland in a dose-dependent manner. Since its discovery, it has been used as a benchmark and starting point in numerous researches aiming to obtain new drugs. Complete resonance assignment of GHRP-6 NMR spectra in both open and cyclic forms are reported, showing some differences to random coil chemical shifts. Connectivities observed in the ROESY spectra indicate spatial proximity between the aromatic residues side-chains in both molecules, as well as between residues DPhe5 and Lys6 sidechains. An ensemble of 10 structures was generated for each one of the molecules, showing RMSD values indicative of nonrandom structures. Molecular Dynamics simulations, both with and without explicit solvent, were carried out for GHRP-6 and its cyclic analogue. Conformational analysis performed on the trajectories showed a nonrandom structure with a well preserved backbone. The presence of geometrical patterns resembling those typical of π-π interactions in both peptides, suggest that this kind of interactions may be relevant for the biological activity of GHRP-6. Same conclusion can be drawn from the spatial proximity of residues DPhe5 and Lys6 sidechains.

Conformational distributions of denatured and unstructured proteins are similar to those of 20 × 20 blocked dipeptides

Understanding intrinsic conformational preferences of amino-acids in unfolded proteins is important for elucidating the underlying principles of their stability and re-folding on biological timescales. Here, to investigate the neighbor interaction effects on the conformational propensities of amino-acids, we carried out (1)H NMR experiments for a comprehensive set of blocked dipeptides and measured the scalar coupling constants between alpha protons and amide protons as well as their chemical shifts. Detailed inspection of these NMR properties shows that, irrespective of amino-acid side-chain properties, the distributions of the measured coupling constants and chemical shifts of the dipeptides are comparatively narrow, indicating small variances of their conformation distributions. They are further compared with those of blocked amino-acids (Ac-X-NHMe), oligopeptides (Ac-GGXGG-NH(2)), and native (lysozyme), denatured (lysozyme and outer membrane protein X from Escherichia coli), unstructured (Domain 2 of the protein 5A of Hepatitis C virus), and intrinsically disordered (hNlg3cyt: intracellular domain of human NL3) proteins. These comparative investigations suggest that the conformational preferences and local solvation environments of the blocked dipeptides are quite similar to not only those of other short oligopeptides but also those of denatured and natively unfolded proteins.

Accurate measurement of alpha proton chemical shifts of excited protein states by relaxation dispersion NMR spectroscopy

Carr-Purcell-Meiboom-Gill relaxation dispersion NMR spectroscopy can provide detailed information about low populated, invisible states of protein molecules, including backbone chemical shifts of the invisible conformer and bond vector orientations that can be used as structural constraints. Notably, the measurement of 1Halpha chemical shifts in excited protein states has not been possible to date because, in the absence of suitable labeling, the homonuclear proton scalar coupling network in side chains of proteins leads to a significant degradation in the performance of proton-based relaxation dispersion experiments. Here we have overcome this problem through a labeling scheme in which proteins are prepared with U-2H glucose and 50% D2O/50% H2O that results in deuteration levels of between 50-88% at the Cbeta carbon. Effects from residual 1Halpha-1Hbeta scalar couplings can be suppressed through a new NMR experiment that is presented here. The utility of the methodology is demonstrated on a ligand binding exchanging system and it is shown that 1Halpha chemical shifts extracted from dispersion profiles are, on average, accurate to 0.03 ppm, an order of magnitude better than they can be predicted from structure using a database approach. The ability to measure 1Halpha chemical shifts of invisible conformers is particularly important because such shifts are sensitive to both secondary and tertiary structure. Thus, the methodology presented is a valuable addition to a growing list of experiments for characterizing excited protein states that are difficult to study using the traditional techniques of structural biology.

Relationship between chemical shift value and accessible surface area for all amino acid atoms

BACKGROUND

Chemical shifts obtained from NMR experiments are an important tool in determining secondary, even tertiary, protein structure. The main repository for chemical shift data is the BioMagResBank, which provides NMR-STAR files with this type of information. However, it is not trivial to link this information to available coordinate data from the PDB for non-backbone atoms due to atom and chain naming differences, as well as sequence numbering changes.

RESULTS

We here describe the analysis of a consistent set of chemical shift and coordinate data, in which we focus on the relationship between the per-atom solvent accessible surface area (ASA) in the reported coordinates and their reported chemical shift value. The data is available online on http://www.ebi.ac.uk/pdbe/docs/NMR/shiftAnalysis/index.html.

CONCLUSION

Atoms with zero per-atom ASA have a significantly larger chemical shift dispersion and often have a different chemical shift distribution compared to those that are solvent accessible. With higher per-atom ASA, the chemical shift values also tend towards random coil values. The per-atom ASA, although not the determinant of the chemical shift, thus provides a way to directly correlate chemical shift information to the atomic coordinates.

Structure and dynamics of the Abeta(21-30) peptide from the interplay of NMR experiments and molecular simulations

We combine molecular dynamics simulations and new high-field NMR experiments to describe the solution structure of the Abeta(21-30) peptide fragment that may be relevant for understanding structural mechanisms related to Alzheimer's disease. By using two different empirical force-field combinations, we provide predictions of the three-bond scalar coupling constants ((3)J(H(N)H(alpha))), chemical-shift values, (13)C relaxation parameters, and rotating-frame nuclear Overhauser effect spectroscopy (ROESY) crosspeaks that can then be compared directly to the same observables measured in the corresponding NMR experiment of Abeta(21-30). We find robust prediction of the (13)C relaxation parameters and medium-range ROESY crosspeaks by using new generation TIP4P-Ew water and Amber ff99SB protein force fields, in which the NMR validates that the simulation yields both a structurally and dynamically correct ensemble over the entire Abeta(21-30) peptide. Analysis of the simulated ensemble shows that all medium-range ROE restraints are not satisfied simultaneously and demonstrates the structural diversity of the Abeta(21-30) conformations more completely than when determined from the experimental medium-range ROE restraints alone. We find that the structural ensemble of the Abeta(21-30) peptide involves a majority population (approximately 60%) of unstructured conformers, lacking any secondary structure or persistent hydrogen-bonding networks. However, the remaining minority population contains a substantial percentage of conformers with a beta-turn centered at Val24 and Gly25, as well as evidence of the Asp23 to Lys28 salt bridge important to the fibril structure. This study sets the stage for robust theoretical work on Abeta(1-40) and Abeta(1-42), for which collection of detailed NMR data on the monomer will be more challenging because of aggregation and fibril formation on experimental timescales at physiological conditions. In addition, we believe that the interplay of modern molecular simulation and high-quality NMR experiments has reached a fruitful stage for characterizing structural ensembles of disordered peptides and proteins in general.

Application of the random coil index to studying protein flexibility

Protein flexibility lies at the heart of many protein-ligand binding events and enzymatic activities. However, the experimental measurement of protein motions is often difficult, tedious and error-prone. As a result, there is a considerable interest in developing simpler and faster ways of quantifying protein flexibility. Recently, we described a method, called Random Coil Index (RCI), which appears to be able to quantitatively estimate model-free order parameters and flexibility in protein structural ensembles using only backbone chemical shifts. Because of its potential utility, we have undertaken a more detailed investigation of the RCI method in an attempt to ascertain its underlying principles, its general utility, its sensitivity to chemical shift errors, its sensitivity to data completeness, its applicability to other proteins, and its general strengths and weaknesses. Overall, we find that the RCI method is very robust and that it represents a useful addition to traditional methods of studying protein flexibility. We have implemented many of the findings and refinements reported here into a web server that allows facile, automated predictions of model-free order parameters, MD RMSF and NMR RMSD values directly from backbone 1H, 13C and 15N chemical shift assignments. The server is available at http://wishart.biology.ualberta.ca/rci.

Spectral fitting for signal assignment and structural analysis of uniformly 13C-labeled solid proteins by simulated annealing based on chemical shifts and spin dynamics

We describe an approach for the signal assignment and structural analysis with a suite of two-dimensional (13)C-(13)C magic-angle-spinning solid-state NMR spectra of uniformly (13)C-labeled peptides and proteins. We directly fit the calculated spectra to experimental ones by simulated annealing in restrained molecular dynamics program CNS as a function of atomic coordinates. The spectra are calculated from the conformation dependent chemical shift obtained with SHIFTX and the cross-peak intensities computed for recoupled dipolar interactions. This method was applied to a membrane-bound 14-residue peptide, mastoparan-X. The obtained C', C(alpha) and C(beta) chemical shifts agreed with those reported previously at the precisions of 0.2, 0.7 and 0.4 ppm, respectively. This spectral fitting program also provides backbone dihedral angles with a precision of about 50 degrees from the spectra even with resonance overlaps. The restraints on the angles were improved by applying protein database program TALOS to the obtained chemical shifts. The peptide structure provided by these restraints was consistent with the reported structure at the backbone RMSD of about 1 A.

Empirical isotropic chemical shift surfaces

A list of proteins is given for which spatial structures, with a resolution better than 2.5 A, are known from entries in the Protein Data Bank (PDB) and isotropic chemical shift (ICS) values are known from the RefDB database related to the Biological Magnetic Resonance Bank (BMRB) database. The structures chosen provide, with unknown uncertainties, dihedral angles phi and psi characterizing the backbone structure of the residues. The joint use of experimental ICSs of the same residues within the proteins, again with mostly unknown uncertainties, and ab initio ICS(phi,psi) surfaces obtained for the model peptides For-(L-Ala)(n)-NH(2), with n = 1, 3, and 5, resulted in so-called empirical ICS(phi,psi) surfaces for all major nuclei of the 20 naturally occurring alpha-amino acids. Out of the many empirical surfaces determined, it is the 13C(alpha)-1H(alpha) ICS(phi,psi) surface which seems to be most promising for identifying major secondary structure types, alpha-helix, beta-strand, left-handed helix (alpha(D)), and polyproline-II. Detailed tests suggest that Ala is a good model for many naturally occurring alpha-amino acids. Two-dimensional empirical 13C(alpha)-1H(alpha) ICS(phi,psi) correlation plots, obtained so far only from computations on small peptide models, suggest the utility of the experimental information contained therein and thus they should provide useful constraints for structure determinations of proteins.

Protein backbone chemical shifts predicted from searching a database for torsion angle and sequence homology

Chemical shifts of nuclei in or attached to a protein backbone are exquisitely sensitive to their local environment. A computer program, SPARTA, is described that uses this correlation with local structure to predict protein backbone chemical shifts, given an input three-dimensional structure, by searching a newly generated database for triplets of adjacent residues that provide the best match in phi/psi/chi(1 )torsion angles and sequence similarity to the query triplet of interest. The database contains (15)N, (1)H(N), (1)H(alpha), (13)C(alpha), (13)C(beta) and (13)C' chemical shifts for 200 proteins for which a high resolution X-ray (< or =2.4 A) structure is available. The relative importance of the weighting factors for the phi/psi/chi(1) angles and sequence similarity was optimized empirically. The weighted, average secondary shifts of the central residues in the 20 best-matching triplets, after inclusion of nearest neighbor, ring current, and hydrogen bonding effects, are used to predict chemical shifts for the protein of known structure. Validation shows good agreement between the SPARTA-predicted and experimental shifts, with standard deviations of 2.52, 0.51, 0.27, 0.98, 1.07 and 1.08 ppm for (15)N, (1)H(N), (1)H(alpha), (13)C(alpha), (13)C(beta) and (13)C', respectively, including outliers.

A Bayesian-probability-based method for assigning protein backbone dihedral angles based on chemical shifts and local sequences

Chemical shifts contain substantial information about protein local conformations. We present a method to assign individual protein backbone dihedral angles into specific regions on the Ramachandran map based on the amino acid sequences and the chemical shifts of backbone atoms of tripeptide segments. The method uses a scoring function derived from the Bayesian probability for the central residue of a query tripeptide segment to have a particular conformation. The Ramachandran map is partitioned into representative regions at two levels of resolution. The lower resolution partitioning is equivalent to the conventional definitions of different secondary structure regions on the map. At the higher resolution level, the alpha and beta regions are further divided into subregions. Predictions are attempted at both levels of resolution. We compared our method with TALOS using the original TALOS database, and obtained comparable results. Although TALOS may produce the best results with currently available databases which are much enlarged, the Bayesian-probability-based approach can provide a quantitative measure for the reliability of predictions.

PREDITOR: a web server for predicting protein torsion angle restraints

Every year between 500 and 1000 peptide and protein structures are determined by NMR and deposited into the Protein Data Bank. However, the process of NMR structure determination continues to be a manually intensive and time-consuming task. One of the most tedious and error-prone aspects of this process involves the determination of torsion angle restraints including phi, psi, omega and chi angles. Most methods require many days of additional experiments, painstaking measurements or complex calculations. Here we wish to describe a web server, called PREDITOR, which greatly accelerates and simplifies this task. PREDITOR accepts sequence and/or chemical shift data as input and generates torsion angle predictions (with predicted errors) for phi, psi, omega and chi-1 angles. PREDITOR combines sequence alignment methods with advanced chemical shift analysis techniques to generate its torsion angle predictions. The method is fast (<40 s per protein) and accurate, with 88% of phi/psi predictions being within 30° of the correct values, 84% of chi-1 predictions being correct and 99.97% of omega angles being correct. PREDITOR is 35 times faster and up to 20% more accurate than any existing method. PREDITOR also provides accurate assessments of the torsion angle errors so that the torsion angle constraints can be readily fed into standard structure refinement programs, such as CNS, XPLOR, AMBER and CYANA. Other unique features to PREDITOR include dihedral angle prediction via PDB structure mapping, automated chemical shift re-referencing (to improve accuracy), prediction of proline cis/trans states and a simple user interface. The PREDITOR website is located at: .

Accurate prediction of protein torsion angles using chemical shifts and sequence homology

Torsion angle restraints are frequently used in the determination and refinement of protein structures by NMR. These restraints may be obtained by J coupling, cross-correlation measurements, nuclear Overhauser effects (NOEs) or secondary chemical shifts. Currently most backbone (phi/psi) torsion angles are determined using a combination of J(HNHalpha) couplings and chemical shift measurements while most side-chain (chi1) angles and cis/trans peptide bond angles (omega) are determined via NOEs. The dependency on multiple experimental (and computational) methods to obtain different torsion angle restraints is both time-consuming and error prone. The situation could be greatly improved if the determination of all torsion angles (phi, psi, chi and omega) could be made via a single type of measurement (i.e. chemical shifts). Here we describe a program, called SHIFTOR, that is able to accurately predict a large number of protein torsion angles (phi, psi, omega, chi1) using only 1H, 13C and 15N chemical shift assignments as input. Overall, the program is 100x faster and its predictions are approximately 20% better than existing methods. The program is also capable of predicting chi1 angles with 81% accuracy and omega angles with 100% accuracy. SHIFTOR exploits many of the recent developments and observations regarding chemical shift dependencies as well as using information in the Protein Databank to improve the quality of its shift-derived torsion angle predictions. SHIFTOR is available as a freely accessible web server at http://wishart.biology.ualberta.ca/shiftor.

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.

Chemical shifts provide fold populations and register of beta hairpins and beta sheets

A detailed analysis of peptide backbone amide (H(N)) and H alpha chemical shifts reveals a consistent pattern for beta hairpins and three-stranded beta sheets. The H alpha's at non-hydrogen-bonded strand positions are inwardly directed and shifted downfield approximately 1 ppm due largely to an anisotropy contribution from the cross-strand amide function. The secondary structure associated H alpha shift deviations for the H-bonded strand positions are also positive but much smaller (0.1-0.3 ppm) and the turn residues display negative H alpha chemical shift deviations (CSDs). The pattern of (H(N)) shift deviations is an even better indicator of both hairpin formation and register, with the cross-strand H-bonded sites shifted downfield (also by approximately 1 ppm) and with diagnostic values for the first turn residue and the first strand position following the turn. These empirical observations, initially made for [2:2]/[2:4]-type-I' and -II' hairpins, are rationalized and can be extended to the analysis of other turns, hairpin classes ([3:5], [4:4]/[4:6]), and three-stranded peptide beta-sheet models. The H alpha's at non-hydrogen-bonded sites and (H(N))'s in the intervening H-bonded sites provide the largest and most dependable measures of hairpin structuring and can be used for melting studies; however the intrinsic temperature dependence of (H(N)) shifts deviations needs to reflect the extent of solvent sequestration in the folded state. Several observations made in the course of this study provide insights into beta-sheet folding mechanisms: (1) The magnitude of the (H(N)) shifts suggests that the cross-strand H-bonds in peptide hairpins are as short as those in protein beta sheets. (2) Even L-Pro-Gly turns, which are frequently used in unfolded controls for hairpin peptides, can support hairpin populations in aqueous fluoroalcohol media. (3) The good correlation between hairpin population estimates from cross-strand H-bonded (H(N)) shift deviations, H alpha shift deviations, and structuring shifts at the turn locus implies that hairpin folding transitions approximate two-state behavior.

Toward direct determination of conformations of protein building units from multidimensional NMR experiments VI. Chemical shift analysis of his to gain 3D structure and protonation state information

NMR--chemical shift structure correlations were investigated by using GIAO-RB3LYP/6-311++G(2d,2p) formalism. Geometries and chemical shifts (CSI values) of 103 different conformers of N'-formyl-L-histidinamide were determined including both neutral and charged protonation forms. Correlations between amino acid torsional angle values and chemical shifts were investigated for the first time for an aromatic and polar amino acid residue whose side chain may carry different charges. Linear correlation coefficients of a significant level were determined between chemical shifts and dihedral angles for CSI[(1)H(alpha)]/phi, CSI[(13)C(alpha)]/phi, and CSI[(13)C(alpha)]/psi. Protonation of the imidazole ring induces the upfield shift of CSI[(13)C(alpha)] for positively charged histidines and an opposite effect for the negative residue. We investigated the correspondence of theoretical and experimental (13)C(alpha), (13)C(beta), and (1)H(alpha) chemical shifts and the nine basic conformational building units characteristic for proteins. These three chemical shift values allow the identification of conformational building units at 80% accuracy. These results enable the prediction of additional regular secondary structural elements (e.g., polyProlineII, inverse gamma-turns) and loops beyond the assignment of chemical shifts to alpha-helices and beta-pleated sheets. Moreover, the location of the His residue can be further specified in a beta-sheet. It is possible to determine whether the appropriate residue is located at the middle or in a first/last beta-strand within a beta-sheet based on calculated CSI values. Thus, the attractive idea of establishing local residue specific backbone folding parameters in peptides and proteins by employing chemical shift information (e.g., (1)H(alpha) and (13)C(alpha)) obtained from selected heteronuclear correlation NMR experiments (e.g., 2D-HSQC) is reinforced.

Electronic properties of amino acid side chains: quantum mechanics calculation of substituent effects

BACKGROUND

Electronic properties of amino acid side chains such as inductive and field effects have not been characterized in any detail. Quantum mechanics (QM) calculations and fundamental equations that account for substituent effects may provide insight into these important properties. PM3 analysis of electron distribution and polarizability was used to derive quantitative scales that describe steric factors, inductive effects, resonance effects, and field effects of amino acid side chains.

RESULTS

These studies revealed that: (1) different semiempirical QM methods yield similar results for the electronic effects of side chain groups, (2) polarizability, which reflects molecular deformability, represents steric factors in electronic terms, and (3) inductive effects contribute to the propensity of an amino acid for alpha-helices.

CONCLUSION

The data provide initial characterization of the substituent effects of amino acid side chains and suggest that these properties affect electron density along the peptide backbone.

Linear analysis of carbon-13 chemical shift differences and its application to the detection and correction of errors in referencing and spin system identifications

Statistical analysis reveals that the set of differences between the secondary shifts of the alpha- and beta-carbons for residues i of a protein (Deltadelta13C(alpha)i - Deltadelta13C(beta)i) provides the means to detect and correct referencing errors for 1H and 13C nuclei within a given dataset. In a correctly referenced protein dataset, linear regression plots of Deltadelta13C(alpha)i, Deltadelta13C(beta)i, or Deltadelta1H(alpha)i vs. (Deltadelta13C(alpha)i - Deltadelta13C(beta)i) pass through the origin from two directions, the helix-to-coil and strand-to-coil directions. Thus, linear analysis of chemical shifts (LACS) can be used to detect referencing errors and to recalibrate the 1H and 13C chemical shift scales if needed. The analysis requires only that the signals be identified with distinct residue types (intra-residue spin systems). LACS allows errors in calibration to be detected and corrected in advance of sequence-specific assignments and secondary structure determinations. Signals that do not fit the linear model (outliers) deserve scrutiny since they could represent errors in identifying signals with a particular residue, or interesting features such as a cis-peptide bond. LACS provides the basis for the automated detection of such features and for testing reassignment hypotheses. Early detection and correction of errors in referencing and spin system identifications can improve the speed and accuracy of chemical shift assignments and secondary structure determinations. We have used LACS to create a database of offset-corrected chemical shifts corresponding to nearly 1800 BMRB entries: 300 with and 1500 without corresponding three-dimensional (3D) structures. This database can serve as a resource for future analysis of the effects of amino acid sequence and protein secondary and tertiary structure on NMR chemical shifts.

Protein chemical shifts arising from alpha-helices and beta-sheets depend on solvent exposure

The NMR chemical shifts of certain atomic nuclei in proteins ((1)H(alpha),(13)C(alpha), and (13)C(beta)) depend sensitively on whether or not the amino acid residue is part of a secondary structure (alpha-helix, beta-sheet), and if so, whether it is helix or sheet. The physical origin of the different chemical shifts of atomic nuclei in alpha-helices versus beta-sheets is a problem of long standing. We report that the chemical shift contributions arising from secondary structure (secondary structure shifts) depend strongly on the extent of exposure to solvent. This behavior is observed for (1)H(alpha), (13)C(alpha), and (13)C(beta) (sheet), but not for(13)C(beta) (helix), whose secondary structure shifts are small. When random coil values are subtracted from the chemical shifts of all(1)H(alpha) nuclei (Pro residues excluded) and the residual chemical shifts are summed to plot the mean values against solvent exposure, the results give a funnel-shaped curve that approaches a small value at full-solvent exposure. When chemical shifts are plotted instead against E(local), the electrostatic contribution to conformational energy produced by local dipole-dipole interactions, a well characterized dependence of (1)H(alpha) chemical shifts on E(local) is found. The slope of this plot varies with both the type of amino acid and the extent of solvent exposure. These results indicate that secondary structure shifts are produced chiefly by the electric field of the protein, which is screened by water dipoles at residues in contact with solvent.

Secondary structural effects on protein NMR chemical shifts

For an amino acid in protein, its chemical shift, delta(phi, psi)(s), is expressed as a function of its backbone torsion angles (phi and psi) and secondary state (s): delta(phi, psi)(s=deltaphi, psi)_coil+Deltadelta(phi, psi)_s), where delta(phi, psi)(coil) represents its chemical shift at coil state (s=coil); Delta delta(phi, psi)(s) (s=sheet or helix) is herein defined as secondary structural effect correction factor, which are quantitatively determined from Residue-specific Secondary Structure Shielding Surface (RSS) for (13)CO, (13)Calpha, (13)Cbeta,(1)Halpha, (15)N, and (1)HN nuclei. The secondary structural effect correction factors defined in this study differ from those in earlier investigations by separating out the backbone conformational effects. As a consequence, their magnitudes are significantly smaller than those earlier reported. The present Delta delta(phi, psi)(sheet) and Delta delta(phi, psi)(helix) were found varying little with backbone conformation and the 20 amino acids, specifically for (13)CO, (13)Calpha, and (1)Halpha nuclei. This study also carries out some useful investigations on other chemical shift prediction approaches - the traditional shielding surfaces, SHIFTS, SHIFTX, PROSHIFT, and identifies some unexpected shortcomings with these methods. It provides some useful insights into understanding protein chemical shifts and suggests a new route to improving chemical shifts prediction. The RSS surfaces were incorporated into the program PRSI.

Structure of the His44 → Ala Single Point Mutant of the Distal Finger Motif of HIV-1 Nucleocapsid Protein: A Combined NMR, Molecular Dynamics Simulation, and Fluorescence Study

The nucleocapsid protein (NCp7) of human immunodeficiency virus type 1 (HIV-1) contains two highly conserved CCHC zinc fingers that strongly bind Zn(2+) through coordination of one His and three Cys residues. It has been suggested that NCp7 function is conformation specific since substitution of any of the zinc coordinating residues in the zinc finger motifs leads to subsequent loss of viral infectivity. To further determine the structural requirements necessary for this specific conformation, we investigated by (1)H 2D NMR and molecular dynamics simulations the structure of the distal finger motif of NCp7 in which the zinc coordinating amino acid, His 44, was substituted by a noncoordinating Ala residue. While the fold of the N-terminal part of this mutated peptide was similar to that of the native peptide, an increased lability and significant conformational changes were observed in the vicinity of the His-to-Ala mutation. Moreover, molecular dynamics simulations suggested a mechanism by which the variant peptide can bind zinc ion even though one zinc-coordinating amino acid was lacking. Using the fluorescence of the naturally occurring Trp37 residue, the binding affinity of the variant peptide to the (TG)(3) model oligonucleotide was found to be decreased by about 2 orders of magnitude with respect with the native peptide. Modeling of the DNA:NCp7 complex using structures of the variant peptide suggests that the residues forming a hydrophobic cleft in the native protein are improperly oriented for efficient DNA binding by the variant peptide.

Quantitative analysis of conformational exchange contributions to 1H-15N multiple-quantum relaxation using field-dependent measurements. Time scale and structural characterization of exchange in a calmodulin C-terminal domain mutant

Multiple-quantum spin relaxation is a sensitive probe for correlated conformational exchange dynamics on microsecond to millisecond time scales in biomolecules. We measured differential 1H-15N multiple-quantum relaxation rates for the backbone amide groups of the E140Q mutant of the C-terminal domain of calmodulin at three static magnetic field strengths. The differential multiple-quantum relaxation rates range between -88.7 and 92.7 s(-1), and the mean and standard deviation are 7.0 +/- 24 s(-1), at a static magnetic field strength of 14.1 T. Together with values of the 1H and 15N chemical shift anisotropies (CSA) determined separately, the field-dependent data enable separation of the different contributions from dipolar-dipolar, CSA-CSA, and conformational exchange cross-correlated relaxation mechanisms to the differential multiple-quantum relaxation rates. The procedure yields precise quantitative information on the dominant conformational exchange contributions observed in this protein. The field-dependent differences between double- and zero-quantum relaxation rates directly benchmark the rates of conformational exchange, showing that these are fast on the chemical shift time scale for the large majority of residues in the protein. Further analysis of the differential 1H-15N multiple-quantum relaxation rates using previously determined exchange rate constants and populations, obtained from 15N off-resonance rotating-frame relaxation data, enables extraction of the product of the chemical shift differences between the resonance frequencies of the 1H and 15N spins in the exchanging conformations, deltasigma(H)deltasigma(N). Thus, information on the 1H chemical shift differences is obtained, while circumventing complications associated with direct measurements of conformational exchange effects on 1H single-quantum coherences in nondeuterated proteins. The method significantly increases the information content available for structural interpretation of the conformational exchange process, partly because deltasigma(H)deltasigma(N) is a signed quantity, and partly because two chemical shifts are probed simultaneously. The present results support the hypothesis that the exchange in the calcium-loaded state of the E140Q mutant involves conformations similar to those of the wild-type apo (closed) and calcium-loaded (open) states.

Asymmetric -helicity loss within a peptide adsorbed onto charged colloidal substrates

A combination of circular dichroism and solution 1H NMR spectroscopy provides a localized description of the distribution of alpha -helical structure within the capped peptide DDDDAAAAARRRR (4DAR5) in aqueous solution and adsorbed onto anionic and cationic colloidal substrates. The adsorption-induced conformational changes are different from those observed upon heating 4DAR5 in solution, in which case the alanine segment remains largely alpha -helical and the transition to a coil structure propagates from the termini. Adsorption is driven by electrostatic complementarity, which places the charged peptide segment adjacent to the substrate of opposite charge. A similar pattern of alpha -helicity loss is observed whether the peptide is adsorbed onto anionic or cationic colloidal silica, despite inverse orientations; significant alpha -helicity loss occurs within the central alanine segment and the terminal arginine segment, whereas alpha -helicity is retained in the aspartate segment. This pattern of adsorption-induced conformational change illustrates the complex and subtle balance among the intramolecular and intermolecular factors that influence the conformations of adsorbed peptides and proteins.

Conformation of secondary amides. A predictive algorithm that correlates DFT-calculated structures and experimental proton chemical shifts

The magnetic deshielding caused by the amido group on CON-CHalpha protons of secondary amides can easily be correlated with DFT-based structures at the B3LYP/6-31G level of theory via a novel algorithm that refines previous models, such as the classical McConnell equation. The shift is given by delta = a + 2.16 cos2(alpha - 35)/d, where alpha denotes the virtual dihedral angle resulting from linking the carbonyl and the alpha-carbons and d is the distance (A) between the shifted proton and the carbonyl oxygen. Notably, in this equation a is a parameter that can be optimized for different solvents, namely, CDCl3, DMSO-d6, and D2O. For the development of these correlations, the preferential conformation of amides is taken from the optimized structures in the gas phase obtained at the DFT level. The deshielding on anti and gauche protons in both rotamers of (Z)-acetamides and E/Z isomers of formamides has been evaluated. This methodology has proved to be highly reliable, allowing us to discard ab initio or DFT conformational arrangements when shifts calculated by the above-mentioned equation differ from the experimental values. Thus, the anti disposition between the CHalpha proton and the N-H bond appears to be the more stable conformation of simple amides. For amides bearing only one proton at Calpha, a local syn minimum can equally be characterized. The rotational barriers around the CON-alkyl bond along with the pyramidalization of the amido group have also been reassessed. As the conformation is taken away from anti or local syn minima, the nonplanarity of the amido group appears to increase.

Probing the conformational and dynamical effects of O-glycosylation within the immunodominant region of a MUC1 peptide tumor antigen

MUC1 mucin is a large transmembrane glycoprotein, the extracellular domain of which is formed by a repeating 20 amino acid sequence, GVTSAPDTRPAPGSTAPPAH. In normal breast epithelial cells, the extracellular domain is densely covered with highly branched complex carbohydrate structures. However, in neoplastic breast tissue, the extracellular domain is under-glycosylated, resulting in the exposure of a highly immunogenic core peptide epitope (PDTRP in bold above), as well as in the exposure of normally cryptic core Tn (GalNAc), STn (sialyl alpha2-6 GalNAc) and TF (Gal beta1-3 GalNAc) carbohydrates. Here, we report the results of 1H NMR structural studies, natural abundance 13C NMR relaxation measurements and distance-restrained MD simulations designed to probe the structural and dynamical effects of Tn-glycosylation within the PDTRP core peptide epitope. Two synthetic peptides were studied: a nine-residue MUC1 peptide of the sequence, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6-Arg7-Pro8-Ala9, and a Tn-glycosylated version of this peptide, Thr1-Ser2-Ala3-Pro4-Asp5-Thr6(alphaGalNAc)-Arg7-Pro8-Ala9. The results of these studies show that a type I beta-turn conformation is adopted by residues PDTR within the PDTRP region of the unglycosylated MUC1 sequence. The existence of a similar beta-turn within the PDTRP core peptide epitope of the under-glycosylated cancer-associated MUC1 mucin protein might explain the immunodominance of this region in vivo, as the presence of defined secondary structure within peptide epitope regions has been correlated with increased immunogenicity in other systems. Our results have also shown that Tn glycosylation at the central threonine within the PDTRP core epitope region shifts the conformational equilibrium away from the type I beta-turn conformation and toward a more rigid and extended state. The significance of these results are discussed in relation to the possible roles that peptide epitope secondary structure and glycosylation state may play in MUC1 tumor immunogenicity.

Hydrogen exchange, core modules, and new designed proteins

A strategy for design of new proteins that mimic folding properties of native proteins is based on peptides modeled on the slow exchange cores of natural proteins. We have synthesized peptides, called core modules, that correspond to the elements of secondary structure that carry the very slowest exchanging amides in a protein. The expectation is that, if soluble in water, core modules will form conformational ensembles that favor native-like structure. Core modules modeled on natural bovine pancreatic trypsin inhibitor have been shown by NMR studies to meet this expectation. The next step toward production of a native state mimic is to further shift the conformational bias of a core module toward more ordered structure by promoting module-module interactions that are mutually stabilizing. For this, two core modules were incorporated into a single molecule by means of a long cross-link. From a panel of several two-module peptides, one very promising lead emerged; it is called BetaCore. BetaCore is monomeric in water and forms a new fold composed of a four-stranded, antiparallel beta-sheet. The single, dominant conformation of BetaCore is characterized by various NMR experiments. Here we compare the individual core module to the two-module BetaCore and discuss the progressive stabilization of intramodule structure and the formation of new intermodule interactions.

Backbone structure confirmation and side chain conformation refinement of a bradykinin mimic BKM-824 by comparing calculated (1)H, (13)C and (19)F chemical shifts with experiment

Calculated and experimental (1)H, (13)C and (19)F chemical shifts were compared in BKM-824, a cyclic bradykinin antagonist mimic, c[Ava(1)-Igl(2)-Ser(3)-DF5F(4)-Oic(5)-Arg(6)] (Ava=5-aminovaleric acid, Igl=alpha-(2-indanyl)glycine, DF5F=pentafluorophenylalanine, Oic=(2S,3aS,7aS)-octahydroindole-2-carboxylic acid). The conformation of BKM-824 has been studied earlier by NMR spectroscopy (M. Miskolzie et al., J. Biomolec. Struct. Dyn. 17, 947-955 (2000)). All NMR structures have qualitatively the same backbone structure but there is considerable variation in the side chain conformations. We have carried out quantum mechanical optimization for three representative NMR structures at the B3LYP/6-31G* level, constraining the backbone dihedral angles at their NMR structure values, followed by NMR chemical shift calculations at the optimized structures with the 6-311G** basis set. There is an intramolecular hydrogen bond at Ser(3) in the optimized structures. The experimental (13)C chemical shifts at five C(alpha) positions as well as at the Cbeta, Cgamma and Cdelta position of Ava(1), which forms part of the backbone, are well reproduced by the calculations, confirming the NMR backbone structure. A comparison between the calculated and experimental H(beta) chemical shifts in Igl(2) shows that the dominant conformation at this residue is gauche. Changes of proton chemical shifts with the scan of the chi(1) angle in DF5F(4) suggest that chi(1)180 degrees. The calculated (1)H and (13)C chemical shifts are in good agreement with experiment at the rigid residue Oic(5). None of the models gives accurate results for Arg(6), presumably because of its positive charge. Our study indicates that calculated NMR shifts can be used as additional constraints in conjunction with NMR data to determine protein conformations. However, to be computationally effective, a database of chemical shifts in small peptide fragments should be precalculated.

BetaCore, a designed water soluble four-stranded antiparallel beta-sheet protein

BetaCore is a designed approximately 50-residue protein in which two BPTI-derived core modules, CM I and CM II, are connected by a 22-atom cross-link. At low temperature and pH 3, homo- and heteronuclear NMR data report a dominant folded ('f') conformation with well-dispersed chemical shifts, i, i+1 periodicity, numerous long-range NOEs, and slowed amide hydrogen isotope exchange patterns that is a four-stranded antiparallel beta-sheet with nonsymmetrical and specific association of CM I and CM II. BetaCore 'f' conformations undergo reversible, global, moderately cooperative, non-two-state thermal transitions to an equilibrium ensemble of unfolded 'u' conformations. There is a significant energy barrier between 'f' and 'u' conformations. This is the first designed four-stranded antiparallel beta-sheet that folds in water.

Cation-pi interaction in a folded polypeptide

Cationic and aromatic side chains from protein residues interact to stabilize tertiary structure. The stabilization energy originates in part from electrostatic attraction between the cation, and regions of high electron density in pi-orbitals of the aromatic group, leading to the name cation-pi interaction. The lysine and tyrosine containing peptide, N-acetyl-Pro-Pro-Lys-Tyr-Asp-Lys-NH(2), has near uv CD characteristic of tyrosine in a structured environment. Nuclear Overhauser effect (NOE), coupling constant, and ring current chemical shift constraints obtained with (1)H NMR confirm that the peptide (t6p) folds. Simulated annealing consistent with all NMR constraints produces a 40-structure ensemble for t6p with potential energies within one standard deviation of the lowest value observed. Calculated binding energies indicate that cation-pi and cation-phenolic OH interactions exists between the Lys3 and Tyr4 side chains in most of the structures. The t6p peptide in solution is a model for these interactions in a protein. A perturbing electric field from the cationic ground state charge intermingles the excited states of the aromatic group. This intermingling effect may provide a cation-pi signature effect in the tyrosine spectroscopy. The absorption and CD for the lowest energy electronic transitions of the tyrosine phenol were computed for the ensemble. Red-shifted peak energy and hypochromicity in the absorbance band, and decreasing rotational strength, correlates with increasing binding energy of the complex indicating the cation-pi spectroscopic signature. The ensemble average spectroscopic signature effects in t6p are small and in agreement with observation.