Use of chemical shifts in macromolecular structure determination

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.

Protein structure determination from NMR chemical shifts

NMR spectroscopy plays a major role in the determination of the structures and dynamics of proteins and other biological macromolecules. Chemical shifts are the most readily and accurately measurable NMR parameters, and they reflect with great specificity the conformations of native and nonnative states of proteins. We show, using 11 examples of proteins representative of the major structural classes and containing up to 123 residues, that it is possible to use chemical shifts as structural restraints in combination with a conventional molecular mechanics force field to determine the conformations of proteins at a resolution of 2 angstroms or better. This strategy should be widely applicable and, subject to further development, will enable quantitative structural analysis to be carried out to address a range of complex biological problems not accessible to current structural techniques.

A new model for chemical shifts of amide hydrogens in proteins

We propose a new computational model to predict amide proton chemical shifts in proteins. In addition to the ring-current, susceptibility and electrostatic effects of earlier models, we add a hydrogen-bonding term based on density functional calculations of model peptide-peptide and peptide-water systems. Both distance and angular terms are included, and the results are rationalized in terms of natural bond orbital analysis of the interactions. Comparison to observed shifts for 15 proteins shows a significant improvement over existing structure-shift correlations. These additions are incorporated in a new version of the SHIFTS program.

A simple genetic algorithm for the optimization of multidomain protein homology models driven by NMR residual dipolar coupling and small angle X-ray scattering data

Most proteins comprise several domains and/or participate in functional complexes. Owing to ongoing structural genomic projects, it is likely that it will soon be possible to predict, with reasonable accuracy, the conserved regions of most structural domains. Under these circumstances, it will be important to have methods, based on simple-to-acquire experimental data, that allow to build and refine structures of multi-domain proteins or of protein complexes from homology models of the individual domains/proteins. It has been recently shown that small angle X-ray scattering (SAXS) and NMR residual dipolar coupling (RDC) data can be combined to determine the architecture of such objects when the X-ray structures of the domains are known and can be considered as rigid objects. We developed a simple genetic algorithm to achieve the same goal, but by using homology models of the domains considered as deformable objects. We applied it to two model systems, an S1KH bi-domain of the NusA protein and the gammaS-crystallin protein. Despite its simplicity our algorithm is able to generate good solutions when driven by SAXS and RDC data.

Use of 13Calpha chemical shifts in protein structure determination

A physics-based method aimed at determining protein structures by using NOE-derived distances together with observed and computed 13C chemical shifts is proposed. The approach makes use of 13Calpha chemical shifts, computed at the density functional level of theory, to obtain torsional constraints for all backbone and side-chain torsional angles without making a priori use of the occupancy of any region of the Ramachandran map by the amino acid residues. The torsional constraints are not fixed but are changed dynamically in each step of the procedure, following an iterative self-consistent approach intended to identify a set of conformations for which the computed 13Calpha chemical shifts match the experimental ones. A test is carried out on a 76-amino acid, all-alpha-helical protein; namely, the Bacillus subtilis acyl carrier protein. It is shown that, starting from randomly generated conformations, the final protein models are more accurate than an existing NMR-derived structure model of this protein, in terms of both the agreement between predicted and observed 13Calpha chemical shifts and some stereochemical quality indicators, and of similar accuracy as one of the protein models solved at a high level of resolution. The results provide evidence that this methodology can be used not only for structure determination but also for additional protein structure refinement of NMR-derived models deposited in the Protein Data Bank.

PWWP module of human hepatoma-derived growth factor forms a domain-swapped dimer with much higher affinity for heparin

Hepatoma-derived growth factor (hHDGF)-related proteins (HRPs) comprise a new growth factor family sharing a highly conserved and ordered N-terminal PWWP module (residues 1-100, previously referred to as a HATH domain) and a variable disordered C-terminal domain. We have shown that the PWWP module is responsible for heparin binding and have solved its structure in solution. Here, we show that under physiological conditions, both the PWWP module and hHDGF can form dimers. Surface plasmon resonance (SPR) studies revealed that the PWWP dimer binds to heparin with affinity that is two orders of magnitude higher (K(d)=13 nM) than that of the monomeric PWWP module (K(d)=1.2 microM). The monomer-dimer equilibrium properties and NMR structural data together suggest that the PWWP dimer is formed through a domain-swapping mechanism. The domain-swapped PWWP dimer structures were calculated on the basis of the NMR data. The results show that the two PWWP protomers exchange their N-terminal hairpin to form a domain-swapped dimer. The two monomers in a dimer are linked by the long flexible L2 loops, a feature supported by NMR relaxation data for the monomer and dimer. The enhanced heparin-binding affinity of the dimer can be rationalized in the framework of the dimer structure.

Identification of a Collapsed Intermediate with Non-native Long-range Interactions on the Folding Pathway of a Pair of Fyn SH3 Domain Mutants by NMR Relaxation Dispersion Spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central beta-sheet is at least partially formed, but without interactions between the NH2- and COOH-terminal beta-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F<-->I<-->U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (< 25 degrees C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F<-->I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.

The mechanism of discrimination between cognate and non-specific DNA by dimeric b/HLH/LZ transcription factors

The Myc/Max/Mad proteins are basic region-helix-loop-helix-leucine zipper (b/HLH/LZ) transcription factors that regulate the transcription of numerous genes involved in cell growth and proliferation. The Max protein is the obligate heterodimeric partner of the Myc and Mad proteins. Heterodimerization and DNA binding to target gene promoters are mediated by the b/HLH/LZ domains. Max can also form a homodimeric b/HLH/LZ. The enhanced expression of Myc and binding to promoters of target genes contribute to almost every aspect of tumor biology. However, the detailed mechanism by which dimeric and heterodimeric b/HLH/LZs discriminate cognate DNA (E-Box: CACGTG) from non-specific sequences in the target gene promoters is still unknown. Here, we use the Max b/HLH/LZ homodimer as a model for this class of transcription factors in the characterization and understanding of the mechanism of discrimination between the E-Box and non-specific DNA sequences. We report the characterization of a cognate and a non-specific Max b/HLH/LZ/DNA complex by EMSA, CD and NMR. Our results support a detailed mechanism by which dimeric b/HLH/LZ transcription factors can discriminate E-Box sequences from non-specific DNA. The mechanism proceeds via the conformational selection of fitting b/HLH/LZ homodimers with the basic region only partially helical. Next, the basic region undergoes a DNA-assisted folding or induced-fit. It is this step that provides the discrimination by stabilizing and destabilizing the alpha-helical conformation of the basic region in the cognate and non-specific complexes, respectively. This leads to a low affinity complex with a higher probability of being dissociated and hence to discrimination. A description of the side-chains and nucleotides proposed to be involved in the discrimination process is provided.

Mapping protein-protein interaction by 13C'-detected heteronuclear NMR spectroscopy

The copper-mediated protein-protein interaction between yeast Atx1 and Ccc2 has been examined by protonless heteronuclear NMR and compared with the already available (1)H-(15)N HSQC information. The observed chemical shift variations are analyzed with respect to the actual solution structure, available through intermolecular NOEs. The advantage of using the CON-IPAP spectrum with respect to the (1)H-(15)N HSQC resides in the increased number of signals observed, including those of prolines. CBCACO-IPAP experiments allow us to focus on the interaction region and on side-chain carbonyls, while a newly designed CEN-IPAP experiment on side-chains of lysines. An attempt is made to rationalize the chemical shift variations on the basis of the structural data involving the interface between the proteins and the nearby regions. It is here proposed that protonless (13)C direct-detection NMR is a useful complement to (1)H based NMR spectroscopy for monitoring protein-protein and protein-ligand interactions.

Determination of the human type I interferon receptor binding site on human interferon-alpha2 by cross saturation and an NMR-based model of the complex

Type I interferons (IFNs) are a family of homologous helical cytokines that exhibit pleiotropic effects on a wide variety of cell types, including antiviral activity and antibacterial, antiprozoal, immunomodulatory, and cell growth regulatory functions. Consequently, IFNs are the human proteins most widely used in the treatment of several kinds of cancer, hepatitis C, and multiple sclerosis. All type I IFNs bind to a cell surface receptor consisting of two subunits, IFNAR1 and IFNAR2, associating upon binding of interferon. The structure of the extracellular domain of IFNAR2 (R2-EC) was solved recently. Here we study the complex and the binding interface of IFNalpha2 with R2-EC using multidimensional NMR techniques. NMR shows that IFNalpha2 does not undergo significant structural changes upon binding to its receptor, suggesting a lock-and-key mechanism for binding. Cross saturation experiments were used to determine the receptor binding site upon IFNalpha2. The NMR data and previously published mutagenesis data were used to derive a docking model of the complex with an RMSD of 1 Angstrom, and its well-defined orientation between IFNalpha2 and R2-EC and the structural quality greatly improve upon previously suggested models. The relative ligand-receptor orientation is believed to be important for interferon signaling and possibly one of the parameters that distinguish the different IFN I subtypes. This structural information provides important insight into interferon signaling processes and may allow improvement in the development of therapeutically used IFNs and IFN-like molecules.

NMR analysis of Caenorhabditis elegans FLP-18 neuropeptides: implications for NPR-1 activation

Phe-Met-Arg-Phe-NH2 (FMRFamide)-like peptides (FLPs) are the largest neuropeptide family in animals, particularly invertebrates. FLPs are characterized by a C-N-terminal gradient of decreasing amino acid conservation. Neuropeptide receptor 1 (NPR-1) is a G-protein coupled receptor (GPCR), which has been shown to be a strong regulator of foraging behavior and aggregation responses in Caenorhabditis elegans. Recently, ligands for NPR-1 were identified as neuropeptides coded by the precursor genes flp-18 and flp-21 in C. elegans. The flp-18 gene encodes eight FLPs including DFDGAMPGVLRF-NH2 and EMPGVLRF-NH2. These peptides exhibit considerably different activities on NPR-1, with the longer one showing a lower potency. We have used nuclear magnetic resonance and biological activity to investigate structural features that may explain these activity differences. Our data demonstrate that long-range electrostatic interactions exist between N-terminal aspartates and the C-terminal penultimate arginine as well as N-terminal hydrogen-bonding interactions that form transient loops within DFDGAMPGVLRF-NH2. We hypothesize that these loops, along with peptide charge, diminish the activity of this peptide on NPR-1 relative to that of EMPGVLRF-NH2. These results provide some insight into the large amino acid diversity in FLPs.

Probabilistic approach to determining unbiased random-coil carbon-13 chemical shift values from the protein chemical shift database

We describe a probabilistic model for deriving, from the database of assigned chemical shifts, a set of random coil chemical shift values that are "unbiased" insofar as contributions from detectable secondary structure have been minimized (RCCSu). We have used this approach to derive a set of RCCSu values for 13Calpha and 13Cbeta for 17 of the 20 standard amino acid residue types by taking advantage of the known opposite conformational dependence of these parameters. We present a second probabilistic approach that utilizes the maximum entropy principle to analyze the database of 13Calpha and 13Cbeta chemical shifts considered separately; this approach yielded a second set of random coil chemical shifts (RCCSmax-ent). Both new approaches analyze the chemical shift database without reference to known structure. Prior approaches have used either the chemical shifts of small peptides assumed to model the random coil state (RCCSpeptide) or statistical analysis of chemical shifts associated with structure not in helical or strand conformation (RCCSstruct-stat). We show that the RCCSmax-ent values are strikingly similar to published RCCSpeptide and RCCSstruct-stat values. By contrast, the RCCSu values differ significantly from both published types of random coil chemical shift values. The differences (RCCSpeptide - RCCSu) for individual residue types show a correlation with known intrinsic conformational propensities. These results suggest that random coil chemical shift values from both prior approaches are biased by conformational preferences. RCCSu values appear to be consistent with the current concept of the "random coil" as the state in which the geometry of the polypeptide ensemble samples the allowed region of (phi, psi)-space in the absence of any dominant stabilizing interactions and thus represent an improved basis for the detection of secondary structure. Coupled with the growing database of chemical shifts, this probabilistic approach makes it possible to refine relationships among chemical shifts, their conformational propensities, and their dependence on pH, temperature, or neighboring residue type.

Improving the accuracy of protein secondary structure prediction using structural alignment

Background

The accuracy of protein secondary structure prediction has steadily improved over the past 30 years. Now many secondary structure prediction methods routinely achieve an accuracy (Q3) of about 75%. We believe this accuracy could be further improved by including structure (as opposed to sequence) database comparisons as part of the prediction process. Indeed, given the large size of the Protein Data Bank (>35,000 sequences), the probability of a newly identified sequence having a structural homologue is actually quite high.

Results

We have developed a method that performs structure-based sequence alignments as part of the secondary structure prediction process. By mapping the structure of a known homologue (sequence ID >25%) onto the query protein's sequence, it is possible to predict at least a portion of that query protein's secondary structure. By integrating this structural alignment approach with conventional (sequence-based) secondary structure methods and then combining it with a "jury-of-experts" system to generate a consensus result, it is possible to attain very high prediction accuracy. Using a sequence-unique test set of 1644 proteins from EVA, this new method achieves an average Q3 score of 81.3%. Extensive testing indicates this is approximately 4–5% better than any other method currently available. Assessments using non sequence-unique test sets (typical of those used in proteome annotation or structural genomics) indicate that this new method can achieve a Q3 score approaching 88%.

Conclusion

By using both sequence and structure databases and by exploiting the latest techniques in machine learning it is possible to routinely predict protein secondary structure with an accuracy well above 80%. A program and web server, called PROTEUS, that performs these secondary structure predictions is accessible at . For high throughput or batch sequence analyses, the PROTEUS programs, databases (and server) can be downloaded and run locally.

Proton character of the peptide unit in the Ca2+-binding sites of calcium pump

The effect of forming calcium pump structures in biological systems on the proton character of the peptide unit has been studied theoretically using the density-functional theory calculations with a large basis set. One acetic acid, one acetate, and three acetamide molecules as well as the modeling peptide unit (MPU) have been employed to mimic the amino acid residues forming the Ca2+-binding sites. To highlight the limiting case of the Ca2+-binding effect on the proton property and the proton countertransport possibility in the direction opposite to the ion, the MPU bounded by the bare or the hydrated Ca2+ has also been investigated. The natural bond orbital (NBO) analysis indicates that the increase of the p-character of the (N-H) sigma orbital results in weakening of the N-H bond which is lengthened when a Ca2+ ion is introduced to the MPU. Calculated NMR shielding sigma(H1) of the MPU shifts upfield upon the Ca2+ ion combination, which reveals the donating of the electron from the amide H as represented by the increase of the calculated positive natural charge for amide H of the MPU. Moreover, the proton affinities (PA) and gas-phase basicities (GB) for the amide nitrogen active site of the MPU are reduced; that is, the acidity of the amide hydrogen gets stronger because of the influence of the Ca2+ ion. To prove the transport possibility of the N-H proton in the direction opposite to the Ca2+ ion along the N-H...O=C hydrogen bond in the helical peptide linkage, NH3 and H2O are used here to assist the dissociation of the amide H of the MPU, and the calculated results show the notable decrease of the deprotonation energies compared to that of the case without this assistance. Moreover, calculated results also reveal that the variation of the quantities discussed here for amide H of the MPU gets smaller when the acidity of Ca2+ ion decreases. Ionization states of the acidic residues forming the Ca2+-binding sites may influence the activity of the amide H of the MPU and further affect the transport tendency of the peptide unit proton in the direction opposite to Ca2+.

Structure of tightly membrane-bound mastoparan-X, a G-protein-activating peptide, determined by solid-state NMR

The structure of mastoparan-X (MP-X), a G-protein activating peptide from wasp venom, in the state tightly bound to anionic phospholipid bilayers was determined by solid-state NMR spectroscopy. Carbon-13 and nitrogen-15 NMR signals of uniformly labeled MP-X were completely assigned by multidimensional intraresidue C-C, N-CalphaCbeta, and N-Calpha-C', and interresidue Calpha-CalphaCbeta, N-CalphaCbeta, and N-C'-Calpha correlation experiments. The backbone torsion angles were predicted from the chemical shifts of 13C', 13Calpha, 13Cbeta, and 15N signals with the aid of protein NMR database programs. In addition, two 13C-13C and three 13C-15N distances between backbone nuclei were precisely measured by rotational resonance and REDOR experiments, respectively. The backbone structure of MP-X was determined from the 26 dihedral angle restraints and five distances with an average root-mean-square deviation of 0.6 A. Peptide MP-X in the bilayer-bound state formed an amphiphilic alpha-helix for residues Trp3-Leu14 and adopted an extended conformation for Asn2. This membrane-bound conformation is discussed in relation to the peptide's activities to form pores in membranes and to activate G-proteins. This study demonstrates the power of multidimensional solid-state NMR of uniformly isotope-labeled molecules and distance measurements for determining the structures of peptides bound to lipid membranes.

The inverted chevron plot measured by NMR relaxation reveals a native-like unfolding intermediate in acyl-CoA binding protein

The folding kinetics of bovine acyl-CoA binding protein was studied by 15N relaxation dispersion measurements under equilibrium conditions. Relaxation dispersion profiles were measured at several concentrations of guanidine hydrochloride (GuHCl). The unfolding rate constant (k(u)) was determined under conditions favoring folding, for which the folding rate constant (k(f)) dominates the relaxation in stopped-flow kinetic measurements. Conversely, k(f) was determined under conditions favoring unfolding, for which k(u) dominates stopped-flow data. The rates determined by NMR therefore complement those from stopped-flow kinetics and define an "inverted chevron" plot. The combination of NMR relaxation and stopped-flow kinetic measurements allowed determination of k(f) and k(u) in the range from 0.48 M GuHCl to 1.28 M GuHCl. Individually, the stopped-flow and NMR data fit two-state models for folding. However, although the values of k(f) determined by the two methods agree, the values of k(u) do not. As a result, a combined analysis of all data does not comply with a two-state model but indicates that an unfolding intermediate exists on the native side of the dominant energy barrier. The denaturant and temperature dependencies of the chemical shifts and k(u) indicate that the intermediate state is structurally similar to the native state. Equilibrium unfolding monitored by optical spectroscopy corroborate these conclusions. The temperature dependence of the chemical shifts identifies regions of the protein that are selectively destabilized in the intermediate. These results illustrate the power of combining stopped-flow kinetics and NMR spectroscopy to analyze protein folding.

Secondary structure and dynamics of micelle bound beta- and gamma-synuclein

We have used solution state NMR spectroscopy to characterize the secondary structure and backbone dynamics of the proteins beta- and gamma-synuclein in their detergent micelle-bound conformations. Comparison of the results with those previously obtained for the Parkinson's disease-linked protein alpha-synuclein shows that structural differences between the three homologous synuclein family members are directly related to variations in their primary amino acid sequences. An 11-residue deletion in the lipid-binding domain of beta-synuclein leads to the destabilization of an entire segment of the micelle-bound helical structure containing the deletion site. The acidic C-terminal tail region of gamma-synuclein, which displays extensive sequence divergence, is more highly disordered than the corresponding regions in the other two family members. The observed structural differences are likely to mediate functional variations between the three proteins, with differences between alpha- and beta-synuclein expected to revolve around their lipid interactions, while differences in gamma-synuclein function are expected to result from different protein-protein interactions mediated by its unique C-terminal tail.

Genetic algorithms as a tool for helix design – computational and experimental studies on prion protein helix 1

Evolutionary computing is a general optimization mechanism successfully implemented for a variety of numeric problems in a variety of fields, including structural biology. We here present an evolutionary approach to optimize helix stability in peptides and proteins employing the AGADIR energy function for helix stability as scoring function. With the ability to apply masks determining positions, which are to remain constant or fixed to a certain class of amino acids, our algorithm is capable of developing stable helical scaffolds containing a wide variety of structural and functional amino acid patterns. The algorithm showed good convergence behaviour in all tested cases and can be parameterized in a wide variety of ways. We have applied our algorithm for the optimization of the stability of prion protein helix 1, a structural element of the prion protein which is thought to play a crucial role in the conformational transition from the cellular to the pathogenic form of the prion protein, and which therefore poses an interesting target for pharmacological as well as genetic engineering approaches to counter the as of yet uncurable prion diseases. NMR spectroscopic investigations of selected stabilizing and destabilizing mutations found by our algorithm could demonstrate its ability to create stabilized variants of secondary structure elements.

NMR study of the tetrameric KcsA potassium channel in detergent micelles

Nuclear magnetic resonance (NMR) studies of large membrane-associated proteins are limited by the difficulties in preparation of stable protein-detergent mixed micelles and by line broadening, which is typical of these macroassemblies. We have used the 68-kDa homotetrameric KcsA, a thermostable N-terminal deletion mutant of a bacterial potassium channel from Streptomyces lividans, as a model system for applying NMR methods to membrane proteins. Optimization of measurement conditions enabled us to perform the backbone assignment of KcsA in SDS micelles and establish its secondary structure, which was found to closely agree with the KcsA crystal structure. The C-terminal cytoplasmic domain, absent in the original structure, contains a 14-residue helix that could participate in tetramerization by forming an intersubunit four-helix bundle. A quantitative estimate of cross- relaxation between detergent and KcsA backbone amide protons, together with relaxation and light scattering data, suggests SDS-KcsA mixed micelles form an oblate spheroid with approximately 180 SDS molecules per channel. K(+) ions bind to the micelle-solubilized channel with a K(D) of 3 +/- 0.5 mM, resulting in chemical shift changes in the selectivity filter. Related pH-induced changes in chemical shift along the "outer" transmembrane helix and the cytoplasmic membrane interface hint at a possible structural explanation for the observed pH-gating of the potassium channel.

Solution structure of the severe acute respiratory syndrome-coronavirus heptad repeat 2 domain in the prefusion state

The envelope glycoprotein, termed the spike protein, of severe acute respiratory syndrome coronavirus (SARS-CoV) is known to mediate viral entry. Similar to other class 1 viral fusion proteins, the heptad repeat regions of SARS-CoV spike are thought to undergo conformational changes from a prefusion form to a subsequent post-fusion form that enables fusion of the viral and host membranes. Recently, the structure of a post-fusion form of SARS-CoV spike, which consists of isolated domains of heptad repeats 1 and 2 (HR1 and HR2), has been determined by x-ray crystallography. To date there is no structural information for the prefusion conformations of SARS-CoV HR1 and HR2. In this work we present the NMR structure of the HR2 domain (residues 1141-1193) from SARS-CoV (termed S2-HR2) in the presence of the co-solvent trifluoroethanol. We find that in the absence of HR1, S2-HR2 forms a coiled coil symmetric trimer with a complex molecular mass of 18 kDa. The S2-HR2 structure, which is the first example of the prefusion form of coronavirus envelope, supports the current model of viral membrane fusion and gives insight into the design of structure-based antagonists of SARS.

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.

Characterizing a Partially Folded Intermediate of the Villin Headpiece Domain Under Non-denaturing Conditions: Contribution of His41 to the pH-dependent Stability of the N-terminal Subdomain

The contribution of interactions involving the imidazole ring of His41 to the pH-dependent stability of the villin headpiece (HP67) N-terminal subdomain has been investigated by nuclear magnetic resonance (NMR) spin relaxation. NMR-derived backbone N-H order parameters (S2) for wild-type (WT) HP67 and H41Y HP67 indicate that reduced conformational flexibility of the N-terminal subdomain in WT HP67 is due to intramolecular interactions with the His41 imidazole ring. These interactions, together with desolvation effects, contribute to significantly depress the pKa of the buried imidazole ring in the native state. 15N R1rho relaxation dispersion data indicate that WT HP67 populates a partially folded intermediate state that is 10.9 kJ mol(-1) higher in free energy than the native state under non-denaturing conditions at neutral pH. The partially folded intermediate is characterized as having an unfolded N-terminal subdomain while the C-terminal subdomain retains a native-like fold. Although the majority of the residues in the N-terminal subdomain sample a random-coil distribution of conformations, deviations of backbone amide 1H and 15N chemical shifts from canonical random-coil values for residues within 5A of the His41 imidazole ring indicate that a significant degree of residual structure is maintained in the partially folded ensemble. The pH-dependence of exchange broadening is consistent with a linear three-state exchange model whereby unfolding of the N-terminal subdomain is coupled to titration of His41 in the partially folded intermediate with a pKa,I=5.69+/-0.07. Although maintenance of residual interactions with the imidazole ring in the unfolded N-terminal subdomain appears to reduce pKa,I compared to model histidine compounds, protonation of His41 disrupts these interactions and reduces the difference in free energy between the native state and partially folded intermediate under acidic conditions. In addition, chemical shift changes for residues Lys70-Phe76 in the C-terminal subdomain suggest that the HP67 actin binding site is disrupted upon unfolding of the N-terminal subdomain, providing a potential mechanism for regulating the villin-dependent bundling of actin filaments.

Multistate Folding of the Villin Headpiece Domain

The villin headpiece (HP67) is a 67 residue, monomeric protein derived from the C-terminal domain of villin. Wild-type HP67 (WT HP67) is the smallest fragment of villin that retains strong in vitro actin-binding activity. WT HP67 is made up of two subdomains, which form a tightly packed interface. The C-terminal subdomain of WT HP67, denoted HP35, is rich in helical structure, folds in isolation, and has been widely used as a model system for folding studies. In contrast, very little is known about the folding of the intact villin headpiece domain. Here, NMR, CD and H/2H amide exchange measurements are used to follow the pH, thermal and urea-induced unfolding of WT HP67 and a mutant (HP67 H41Y) in which a buried conserved histidine in the N-terminal subdomain, His41, has been mutated to Tyr. Although most small proteins display two-state equilibrium unfolding, the results presented here demonstrate that unfolding of the villin headpiece is a multistate process. The presence of a folded N-terminal subdomain is shown to stabilize the C-terminal subdomain, increasing the midpoints of the thermal and urea-induced unfolding transitions and increasing protection factors for H/2H exchange. Histidine 41 has been shown to act as a pH-dependent switch in wild-type HP67: the N-terminal subdomain is unfolded when His41 is protonated, while the C-terminal subdomain remains folded irrespective of the protonation state of His41. Mutation of His41 to Tyr eliminates the segmental pH-dependent unfolding of the headpiece. The mutation stabilizes both domains, but folding is still multistate, indicating that His41 is not solely responsible for the unusual equilibrium unfolding behavior of villin headpiece domain.

Proposal for sets of se NMR chemical shifts in planar and perpendicular orientations of aryl group and the applications

The orientational effect of p-YC(6)H(4) (Ar) on delta(Se) is elucidated for ArSeR, based on experimental and theoretical investigations. Sets of delta(Se) are proposed for pl and pd employing 9-(arylselanyl)anthracenes (1) and 1-(arylselanyl)anthraquinones (2), respectively, where Se-C(R) in ArSeR is on the Ar plane in pl and perpendicular to the plane in pd. Absolute magnetic shielding tensors of Se (sigma(Se)) are calculated for ArSeR (R = H, Me, and Ph), assuming pl and pd, with the DFT-GIAO method. Observed characters are well reproduced by the total shielding tensors (sigma(t)(Se)). The paramagnetic terms (sigma(P)(Se)) are governed by sigma(P)(Se)(xx) + sigma(P)(Se)(yy), where the direction of n(P)(Se) is set to the z-axis. The mechanisms of the orientational effect are established both for pl and pd. Sets of delta(Se: 1) and delta(Se: 2) act as the standards for pl and pd, respectively, when delta(Se) of ArSeR are analyzed based on the orientational effect.

Solution structure of the N-terminal A domain of the human voltage-gated Ca2+channel beta4a subunit

Ca2+ channel beta subunits regulate trafficking and gating (opening and closing) of voltage-dependent Ca2+ channel alpha1 subunits. Based on primary sequence comparisons, they are thought to be modular structures composed of five domains (A-E) that are related to the large family of membrane associated guanylate-kinase (MAGUK) proteins. The crystal structures of the beta subunit core, B-D, domains have recently been reported; however, very little is known about the structures of the A and E domains. The N-terminal A domain is a hypervariable region that differs among the four subtypes of Ca2+ channel beta subunits (beta1-beta4). Furthermore, this domain undergoes alternative splicing to create multiple N-terminal structures within a given gene class that have distinct effects on gating. We have solved the solution structure of the A domain of the human beta4a subunit, a splice variant that we have shown previously to have alpha1 subunit subtype-specific effects on Ca2+ channel trafficking and gating.

Solution structure of (gamma)S-crystallin by molecular fragment replacement NMR

The solution structure of murine gammaS-crystallin (gammaS) has been determined by multidimensional triple resonance NMR spectroscopy, using restraints derived from two sets of dipolar couplings, recorded in different alignment media, and supplemented by a small number of NOE distance restraints. gammaS consists of two topologically similar domains, arranged with an approximate twofold symmetry, and each domain shows close structural homology to closely related (approximately 50% sequence identity) domains found in other members of the gamma-crystallin family. Each domain consists of two four-strand "Greek key" beta-sheets. Although the domains are tightly anchored to one another by the hydrophobic surfaces of the two inner Greek key motifs, the N-arm, the interdomain linker and several turn regions show unexpected flexibility and disorder in solution. This may contribute entropic stabilization to the protein in solution, but may also indicate nucleation sites for unfolding or other structural transitions. The method used for solving the gammaS structure relies on the recently introduced molecular fragment replacement method, which capitalizes on the large database of protein structures previously solved by X-ray crystallography and NMR.

Estimation of protein secondary structure content directly from NMR spectra using an improved empirical correlation with averaged chemical shift

We have recently shown that the averaged chemical shift (ACS) of a nucleus in the protein backbone correlates well empirically to its secondary structure content (SSC). This allows the estimation of SSC directly from the NMR spectrum without the time intensive process of chemical shift assignment. Here, we present an empirical correlation that accounts both for contributions to the relevant protein and chemical shift databases made subsequent to the original analysis, and for missing or inconsistently referenced resonances. Our results affirm that this method provides a significant tool for initial structural prediction from NMR data prior to complete chemical shift assignment.

Correlated dynamics of consecutive residues reveal transient and cooperative unfolding of secondary structure in proteins

Nuclear spin relaxation is a powerful method for studying molecular dynamics at atomic resolution. Recent methods development in biomolecular NMR spectroscopy has enabled detailed investigations of molecular dynamics that are critical for biological function, with prominent examples addressing allostery, enzyme catalysis, and protein folding. Dynamic processes with similar correlation times are often detected in multiple locations of the molecule, raising the question of whether the underlying motions are correlated (corresponding to concerted fluctuations involving many atoms distributed across extended regions of the molecule) or uncorrelated (corresponding to independent fluctuations involving few atoms in localized regions). Here, we have used (13)C(alpha)(i - 1)/(13)C(alpha)(i) differential multiple-quantum spin relaxation to provide direct evidence for correlated dynamics of consecutive amino acid residues in the protein sequence. By monitoring overlapping pairs of residues (i - 1 and i, i and i + 1, etc.), we identified correlated motions that extend through continuous segments of the sequence. We detected significant correlated conformational transitions in the native state of the E140Q mutant of the calmodulin C-terminal domain. Previous work has shown that this domain exchanges between two major conformational states that resemble the functionally relevant open and closed states of the WT protein, with a mean correlation time of approximately 20 micros. The present results reveal that an entire alpha-helix undergoes partial unraveling in a transient and cooperative manner.

A NMR Study of the Interaction of a Three-domain Construct of ATP7A with Copper(I) and Copper(I)-HAH1

ATP7A is a P-type ATPase involved in copper(I) homeostasis in humans. It possesses a long N-terminal tail protruding into the cytosol and containing six copper(I)-binding domains, which are individually folded and capable of binding one copper(I) ion. ATP7A receives copper from a soluble protein, the metallochaperone HAH1. The exact role and interplay of the six soluble domains is still quite unclear, as it has been extensively demonstrated that they are strongly redundant with respect to copper(I) transport in vivo. In the present work, a three-domain (fourth to sixth, MNK456) construct has been investigated in solution by NMR, in the absence and presence of copper(I). In addition, the interaction of MNK456 with copper(I)-HAH1 has been studied. It is proposed that the fourth domain is the preferential site for the initial interaction with the partner. A significant dependence of the overall domain dynamics on the metallation state and on the presence of HAH1 is observed. This dependence could constitute the molecular mechanism to trigger copper(I) translocation and/or ATP7A relocalization from the trans-Golgi network to the plasmatic membrane.

Secondary structure and dynamics of an intrinsically unstructured linker domain

The transient secondary structure and dynamics of an intrinsically unstructured linker domain from the 70 kDa subunit of human replication protein A was investigated using solution state NMR. Stable secondary structure, inferred from large secondary chemical shifts, was observed for a segment of the intrinsically unstructured linker domain when it is attached to an N-terminal protein interaction domain. Results from NMR relaxation experiments showed the rotational diffusion for this segment of the intrinsically unstructured linker domain to be correlated with the N-terminal protein interaction domain. When the N-terminal domain is removed, the stable secondary structure is lost and faster rotational diffusion is observed. The large secondary chemical shifts were used to calculate phi and psi dihedral angles and these dihedral angles were used to build a backbone structural model. Restrained molecular dynamics were performed on this new structure using the chemical shift based dihedral angles and a single NOE distance as restraints. In the resulting family of structures a large, solvent exposed loop was observed for the segment of the intrinsically unstructured linker domain that had large secondary chemical shifts.

Structure-function studies of the functional and binding epitope of the human 37 kDa laminin receptor precursor protein

Expression of the 37 kDa laminin receptor precursor protein (37LRP) correlates directly with increased invasiveness and the metastatic potential of tumors. The 37LRP matures to a 67 kDa protein which facilitates the binding of cancer cells to basement membranes. The palindrome peptide sequence LMWWML, corresponding to the 173-178-residue stretch of the human 37LRP sequence, has been identified as the laminin-1-binding site. Peptides from 37LRP of species that contain this palindrome-bind laminin-1 with high affinity. Nuclear magnetic resonance (NMR) conformational studies have been undertaken on a synthetic 15-residue peptide (KGAHSVGLMWWMLAR) containing the palindrome to establish the structural basis of this activity. To further correlate the structural data with laminin-1-binding function, analogous structural studies were conducted for a similar peptide (RGKHSIGLIWYLLAR) lacking the palindrome, originating from 37LRP sequence of Saccharomyces cerevisiae and exhibiting low laminin-1-binding affinity. Finally, in vitro cell invasion assays were performed to investigate the possibility that the laminin-1-binding affinity of the peptides influences their inhibitory activity.

Comparison of structure and dynamics of micelle-bound human alpha-synuclein and Parkinson disease variants

Three point mutations (A30P, E46K, and A53T) as well as gene triplication genetically link the 140-residue protein alpha-synuclein (aS) to the development of Parkinson disease. Here, the structure and dynamics of micelle-bound aS(A30P) and aS(A53T) are described and compared with wild-type aS, in addition to describing the aS-micelle interaction. A53T is sensed only by directly adjacent residues and leaves the backbone structure and dynamics indistinguishable from the wild type. A30P interrupts one helix turn (Val26-Ala29) and destabilizes the preceding one. A shift in helix register following A30P disturbs the canonical succession of polar and hydrophobic residues for at least two turns. The shortened helix-N adopts a slightly higher helical content and is less bent, indicating that strain was present in the micelle-bound helix. In the vicinity of the A30P-induced perturbations, the underlying micelle environment has rearranged, but nevertheless all aS variants maintain similar interrelationships with the micelle. Moreover, aS-micelle immersion correlates well with fast and slow aS backbone dynamics, allowing a rare insight into protein-micelle interplay.

Microsecond Protein Dynamics Measured by 13Cα Rotating-Frame Spin Relaxation

NMR spin relaxation in the rotating frame (R1rho) is a unique method for atomic-resolution characterization of conformational (chemical) exchange processes occurring on the microsecond timescale. We present a rotating-frame 13C(alpha) relaxation dispersion experiment for measuring conformational dynamics in uniformly 13C-labeled proteins. The experiment was validated by using the E140Q mutant of the C-terminal fragment of calmodulin, which exhibits significant conformational exchange between two major conformations, as gauged from previous 15N and 1H relaxation studies. Consistent with previous work, the present 13C(alpha) R1rho experiment detects conformational-exchange dynamics throughout the protein. The average correlation time of <tau(ex)>=25+/-8 micros is in excellent agreement with those determined previously from 1H and 15N R1rho relaxation data: <tau(ex)>=19+/-7 and 21+/-3 micros, respectively. The extracted chemical-shift differences between the exchanging states reveal significant fluctuations in dihedral angles within single regions of Ramachandran phi-psi space, that were not identified from the 1H and 15N relaxation data. The present results underscore the advantage of using several types of nuclei to probe exchange dynamics in biomolecules.

Microsecond timescale backbone conformational dynamics in ubiquitin studied with NMR R1rho relaxation experiments

NMR spin relaxation experiments are used to characterize the dynamics of the backbone of ubiquitin. Chemical exchange processes affecting residues Ile 23, Asn 25, Thr 55, and Val 70 are characterized using on- and off-resonance rotating-frame 15N R1rho relaxation experiments to have a kinetic exchange rate constant of 25,000 sec(-1) at 280 K. The exchange process affecting residues 23, 25, and 55 appears to result from disruption of N-cap hydrogen bonds of the alpha-helix and possibly from repacking of the side chain of Ile 23. Chemical exchange processes affecting other residues on the surface of ubiquitin are identified using 1H-15N multiple quantum relaxation experiments. These residues are located near or at the regions known to interact with various enzymes of the ubiquitin-dependent protein degradation pathway.

Solution structure and dynamics of integral membrane proteins by NMR: a case study involving the enzyme PagP

Solution NMR spectroscopy is rapidly becoming an important technique for the study of membrane protein structure and dynamics. NMR experiments on large perdeuterated proteins typically exploit the favorable relaxation properties of backbone amide (15)N-(1)H groups to obtain sequence-specific chemical shift assignments, structural restraints, and a wide range of dynamics information. These methods have proven successful in the study of the outer membrane enzyme, PagP, not only for obtaining the global fold of the protein but also for characterizing in detail the conformational fluctuations that are critical to its activity. NMR methods can also be extended to take advantage of slowly relaxing methyl groups, providing additional probes of structure and dynamics at side chain positions. The current work on PagP demonstrates how solution NMR can provide a unique atomic resolution description of the dynamic processes that are key to the function of many membrane protein systems.

Off-resonance rotating-frame amide proton spin relaxation experiments measuring microsecond chemical exchange in proteins

NMR spin relaxation in the rotating frame (R(1rho)) is a unique method for atomic-resolution characterization of conformational (chemical) exchange processes occurring on the microsecond time scale. Here, we use amide 1H off-resonance R(1rho) relaxation experiments to determine exchange parameters for processes that are significantly faster than those that can be probed using 15N or 13C relaxation. The new pulse sequence is validated using the E140Q mutant of the C-terminal domain of calmodulin, which exhibits significant conformational exchange contributions to the transverse relaxation rates. The 1H off-resonance R(1rho) data sample the entire relaxation dispersion profiles for the large majority of residues in this protein, which exchanges between conformations with a time constant of approximately 20 micros. This is in contrast to the case for 15N, where additional laboratory-frame relaxation data are required to determine the exchange parameters reliably. Experiments were performed on uniformly 15N-enriched samples that were either highly enriched in 2H or fully protonated. In the latter case, dipolar cross-relaxation with aliphatic protons were effectively decoupled to first order using a selective inversion pulse. Deuterated and protonated samples gave the same results, within experimental errors. The use of deuterated samples increases the sensitivity towards exchange contributions to the 1H transverse relaxation rates, since dipolar relaxation is greatly reduced. The exchange correlation times determined from the present 1H off-resonance R(1rho) experiments are in excellent agreement with those determined previously using a combination of 15N laboratory-frame and off-resonance R(1rho) relaxation data, with average values of [see text] and 21 +/- 3 micros , respectively.

Effect of urea on peptide conformation in water: molecular dynamics and experimental characterization

Molecular dynamics simulations of a ribonuclease A C-peptide analog and a sequence variant were performed in water at 277 and 300 K and in 8 M urea to clarify the molecular denaturation mechanism induced by urea and the early events in protein unfolding. Spectroscopic characterization of the peptides showed that the C-peptide analog had a high alpha-helical content, which was not the case for the variant. In the simulations, interdependent side-chain interactions were responsible for the high stability of the alpha-helical C-peptide analog in the different solvents. The other peptide displayed alpha-helical unwinding that propagated cooperatively toward the N-terminal. The conformations sampled by the peptides depended on their sequence and on the solvent. The ability of water molecules to form hydrogen bonds to the peptide as well as the hydrogen bond lifetimes increased in the presence of urea, whereas water mobility was reduced near the peptide. Urea accumulated in excess around the peptide, to which it formed long-lived hydrogen bonds. The unfolding mechanisms induced by thermal denaturation and by urea are of a different nature, with urea-aqueous solutions providing a better peptide solvation than pure water. Our results suggest that the effect of urea on the chemical denaturation process involves both the direct and indirect mechanisms.

Protein energetic conformational analysis from NMR chemical shifts (PECAN) and its use in determining secondary structural elements

We present an energy model that combines information from the amino acid sequence of a protein and available NMR chemical shifts for the purposes of identifying low energy conformations and determining elements of secondary structure. The model ("PECAN", Protein Energetic Conformational Analysis from NMR chemical shifts) optimizes a combination of sequence information and residue-specific statistical energy function to yield energetic descriptions most favorable to predicting secondary structure. Compared to prior methods for secondary structure determination, PECAN provides increased accuracy and range, particularly in regions of extended structure. Moreover, PECAN uses the energetics to identify residues located at the boundaries between regions of predicted secondary structure that may not fit the stringent secondary structure class definitions. The energy model offers insights into the local energetic patterns that underlie conformational preferences. For example, it shows that the information content for defining secondary structure is localized about a residue and reaches a maximum when two residues on either side are considered. The current release of the PECAN software determines the well-defined regions of secondary structure in novel proteins with assigned chemical shifts with an overall accuracy of 90%, which is close to the practical limit of achievable accuracy in classifying the states.

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.

Position dependence of the 13C chemical shifts of alpha-helical model peptides. Fingerprint of the 20 naturally occurring amino acids

The position dependence of the (13)C chemical shifts was investigated at the density functional level for alpha-helical model peptides represented by the sequence Ac-(Ala)(i)-X-(Ala)(j)-NH(2), where X represents any of the 20 naturally occurring amino acids, with 0 < or = i < or = 8 and i + j = 8. Adoption of the locally dense basis 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 20 naturally occurring amino acids in alpha-helices, there is (1) significant variability of the computed (13)C shielding as a function of both the guest residue (X) and the position along the sequence; for example, at the N terminus, the (13)C(alpha) and (13)C(beta) shieldings exhibit a uniform pattern of variation with respect to both the central or the C-terminal positions; (2) good agreement between computed and observed (13)C(alpha) and (13)C(beta) chemical shifts in the interior of the helix, with correlation coefficients of 0.98 and 0.99, respectively; for (13)C(alpha) chemical shifts, computed in the middle of the helix, only five residues, namely Asn, Asp, Ser, Thr, and Leu, exhibit chemical shifts beyond the observed standard deviation; and (3) better agreement for four of these residues (Asn, Asp, Ser, and Thr) only for the computed values of the (13)C(alpha) chemical shifts at the N terminus. The results indicate that (13)C(beta), but not (13)C(beta), chemical shifts are sensitive enough to reflect the propensities of some amino acids for specific positions within an alpha-helix, relative to the N and C termini of peptides and proteins.

NMR chemical shift and relaxation measurements provide evidence for the coupled folding and binding of the p53 transactivation domain

The interaction between the acidic transactivation domain of the human tumor suppressor protein p53 (p53TAD) and the 70 kDa subunit of human replication protein A (hRPA70) was investigated using heteronuclear magnetic resonance spectroscopy. A 1H-15N heteronuclear single quantum coherence (HSQC) titration experiment was performed on a 15N-labeled fragment of hRPA70, containing the N-terminal 168 residues (hRPA701-168) and p53TAD. HRPA701-168 residues important for binding were identified and found to be localized to a prominent basic cleft. This binding site overlapped with a previously identified single-stranded DNA-binding site, suggesting that a competitive binding mechanism may regulate the formation of p53TAD-hRPA70 complex. The amide 1H and 15N chemical shifts of an uniformly 15N-labeled sample of p53TAD were also monitored before and after the addition of unlabeled hRPA701-168. In the presence of unlabeled hRPA701-168, resonance lineshapes increased and corresponding intensity reductions were observed for specific p53TAD residues. The largest intensity reductions were observed for p53TAD residues 42-56. Minimal binding was observed between p53TAD and a mutant form of hRPA701-168, where the basic cleft residue R41 was changed to a glutamic acid (R41E), demonstrating that ionic interactions play an important role in specifying the binding interface. The region of p53TAD most affected by binding hRPA701-168 was found to have some residual alpha helical and beta strand structure; however, this structure was not stabilized by binding hRPA701-168. 15N relaxation experiments were performed to monitor changes in backbone dynamics of p53TAD when bound to hRPA701-168. Large changes in both the transverse (R2) and rotating frame (R1) relaxation rates were observed for a subset of the p53TAD residues that had 1H-15N HSQC resonance intensity reductions during the complex formation. The folding of p53TAD upon complex formation is suggested by the pattern of changes observed for both R2 and R1. A model that couples the formation of a weak encounter complex between p53TAD and hRPA701-168 to the folding of p53TAD is discussed in the context of a functional role for the p53-hRPA70 complex in DNA repair.

Solution structure of synthetic penaeidin-4 with structural and functional comparisons with penaeidin-3

Antimicrobial peptide structure has direct implications for the complexity of functions and mechanisms of action. The penaeidin antimicrobial peptide family from shrimp is divided into multiple class designations based on primary structure. The penaeidin classes are not only characterized by variability in primary sequence but also by variation in target specificity and effectiveness. Whereas class 4 exhibits low isoform diversity within species and is highly conserved between species, the primary sequence of penaeidin class 3 is less conserved between species and exhibits considerable isoform diversity within species. All penaeidins, regardless of class or species, are composed of two dramatically different domains: an unconstrained proline-rich domain and a disulfide bond-stabilized cysteine-rich domain. The proline-rich domain varies in length and is generally less conserved, whereas the spacing and specific residue content of the cysteine-rich domain is more conserved. The structure of the synthetic penaeidin class 4 (PEN4-1) from Litopenaeus setiferus was analyzed using several approaches, including chemical mapping of disulfide bonds, circular dichroism analysis of secondary structural characteristics, and complete characterization of the solution structure of the peptide by proton NMR. L. setiferus PEN4-1 was then compared with the previously characterized structure of penaeidin class 3 from Litopenaeus vannamei. Moreover, the specificity of these antimicrobial peptides was examined through direct comparison of activity against a panel of microbes. The penaeidin classes differ in microbial target specificity, which correlates to variability in specific domain sequence. However, the tertiary structure of the cysteine-rich domain and indeed the overall structure of penaeidins are conserved across classes.

A simple method to adjust inconsistently referenced 13C and 15N chemical shift assignments of proteins

Inconsistent 13C and 15N chemical shift referencing is a continuing problem associated with protein chemical shift assignments deposited in BioMagResBank (BMRB). Here we describe a simple and robust approach that can quantitatively determine the 13C and 15N referencing offsets solely from chemical shift assignment data and independently of 3D coordinate data. This novel structure-independent approach permitted the assessment and determination of 13C and 15N reference offsets for all protein entries deposited in the BMRB. Tests on 452 proteins with known 3D structures show that this structure-independent approach yields 13C and 15N referencing offsets that exhibit excellent agreement with those calculated on the basis of 3D structures. Furthermore, this protocol appears to improve the accuracy of chemical shift-derived secondary structural identification, and has been formally incorporated into a computer program called PSSI (http//www.pronmr.com).

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.

How Does an Amide-N Chemical Shift Tensor Vary in Peptides?

This study addresses a void in the existing literature on the amide-(15)N chemical shift anisotropy (CSA) tensor of peptides: a systematic investigation of how the tensor varies in different peptides. Amide-(15)N CSA tensors for several dipeptides are obtained using quantum chemical calculations, as well as for a series of model Ala-X and X-Ala sequences in both α-helical and β-sheet conformations (where X is one of the naturally occurring amino acids). The calculated values show a significant variation in both isolated and extended peptide structures. Hydrogen bonding at both the carbonyl group and the N-H bond of the peptide plane is shown to affect the principal values of the tensor. Calculations on model peptides indicate that the amide-(15)N CSA tensor is dependent on atoms located within a distance of five bonds. Consequently, the tensor of a given peptide residue is unaffected by residues other than those adjacent to it, which implies that the amide-(15)N CSA tensor should be considered in the context of tripeptide sequences. This further suggests that the amide-(15)N CSA tensor of the second residue of a given tripeptide sequence may be extrapolated to the same sequence in any other polypeptide or protein, given the same backbone conformation and intermolecular environment. These conclusions will facilitate future NMR structural studies of proteins.

The NMR solution structure of a mutant of the Max b/HLH/LZ free of DNA: insights into the specific and reversible DNA binding mechanism of dimeric transcription factors

Basic region-helix1-loop-helix2-leucine zipper (b/H(1)LH(2)/LZ) transcription factors bind specific DNA sequence in their target gene promoters as dimers. Max, a b/H(1)LH(2)/LZ transcription factor, is the obligate heterodimeric partner of the related b/H(1)LH(2)/LZ proteins of the Myc and Mad families. These heterodimers specifically bind E-box DNA sequence (CACGTG) to activate (e.g. c-Myc/Max) and repress (e.g. Mad1/Max) transcription. Max can also homodimerize and bind E-box sequences in c-Myc target gene promoters. While the X-ray structure of the Max b/H(1)LH(2)/LZ/DNA complex and that of others have been reported, the precise sequence of events leading to the reversible and specific binding of these important transcription factors is still largely unknown. In order to provide insights into the DNA binding mechanism, we have solved the NMR solution structure of a covalently homodimerized version of a Max b/H(1)LH(2)/LZ protein with two stabilizing mutations in the LZ, and characterized its backbone dynamics from (15)N spin-relaxation measurements in the absence of DNA. Apart from minor differences in the pitch of the LZ, possibly resulting from the mutations in the construct, we observe that the packing of the helices in the H(1)LH(2) domain is almost identical to that of the two crystal structures, indicating that no important conformational change in these helices occurs upon DNA binding. Conversely to the crystal structures of the DNA complexes, the first 14 residues of the basic region are found to be mostly unfolded while the loop is observed to be flexible. This indicates that these domains undergo conformational changes upon DNA binding. On the other hand, we find the last four residues of the basic region form a persistent helical turn contiguous to H(1). In addition, we provide evidence of the existence of internal motions in the backbone of H(1) that are of larger amplitude and longer time-scale (nanoseconds) than the ones in the H(2) and LZ domain. Most interestingly, we note that conformers in the ensemble of calculated structures have highly conserved basic residues (located in the persistent helical turn of the basic region and in the loop) known to be important for specific binding in a conformation that matches that of the DNA-bound state. These partially prefolded conformers can directly fit into the major groove of DNA and as such are proposed to lie on the pathway leading to the reversible and specific DNA binding. In these conformers, the conserved basic side-chains form a cluster that elevates the local electrostatic potential and could provide the necessary driving force for the generation of the internal motions localized in the H(1) and therefore link structural determinants with the DNA binding function. Overall, our results suggests that the Max homodimeric b/H(1)LH(2)/LZ can rapidly and preferentially bind DNA sequence through transient and partially prefolded states and subsequently, adopt the fully helical bound state in a DNA-assisted mechanism or induced-fit.

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.

Peptaibol zervamicin IIb structure and dynamics refinement from transhydrogen bond J couplings

Zervamicin IIB (Zrv-IIB) is a channel-forming peptaibol antibiotic of fungal origin. The measured transhydrogen bond (3h)J(NC') couplings in methanol solution heaving average value of -0.41 Hz indicate that the stability of the Zrv-IIB helix in this milieu is comparable to the stability of helices in globular proteins. The N-terminus of the peptide forms an alpha-helix, whereas 3(10)-helical hydrogen bonds stabilize the C-terminus. However, two weak transhydrogen bond peaks are observed in a long-range HNCO spectrum for HN Aib(12). Energy calculations using the Empirical Conformation Energy Program for Peptides (ECEPP)/2 force field and the implicit solvent model show that the middle of the peptide helix accommodates a bifurcated hydrogen bond that is simultaneously formed between HN Aib(12) and CO Leu(8) and CO Aib(9). Several lowered (3h)J(NC') on a polar face of the helix correlate with the conformational exchange process observed earlier and imply dynamic distortions of a hydrogen bond pattern with the predominant population of a properly folded helical structure. The refined structure of Zrv-IIB on the basis of the observed hydrogen bond pattern has a small ( approximately 20 degrees ) angle of helix bending that is virtually identical to the angle of bending in dodecylphosphocholine (DPC) micelles, indicating the stability of a hinge region in different environments. NMR parameters ((1)HN chemical shifts and transpeptide bond (1)J(NC') couplings) sensitive to hydrogen bonding along with the solvent accessible surface area of carbonyl oxygens indicate a large polar patch on the convex side of the helix formed by three exposed backbone carbonyls of Aib(7), Aib(9), and Hyp(10) and polar side chains of Hyp(10), Gln(11), and Hyp(13). The unique structural features, high helix stability and the enhanced polar patch, set apart Zrv-IIB from other peptaibols (for example, alamethicin) and possibly underlie its biological and physiological properties.

Polyproline II helix conformation in a proline-rich environment: a theoretical study

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP(II) structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro)(3)-X-(Pro)(3)-Gly-Tyr-NH(2), with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro)(3)-Ala-Ala-(Pro)(3)-Gly-Tyr-NH(2) (PAAP), by using a theoretical approach that includes a solvent effect as well as cis <--> trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since (13)C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles (phi,psi) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged (13)C chemical shifts for the guest residues in the PXP peptide (X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.