Use of chemical shifts and coupling constants in nuclear magnetic resonance structural studies on peptides and proteins

Formerly degenerate seventh zinc finger domain from transcription factor ZNF711 rehabilitated by experimental NMR structure

Domain Z7 of nuclear transcription factor ZNF711 has the consensus last metal-ligand H23 found in odd-numbered zinc-fingers of this protein replaced by a phenylalanine. Ever since the discovery of ZNF711 it has been thought that Z7 is probably non-functional because of the H23F substitution. The presence of H26 three positions downstream prompted us to examine if this histidine could substitute as the last metal ligand. The Z7 domain adopts a stable tertiary structure upon metal binding. The NMR structure of Zn2+-bound Z7 shows the classical ββα-fold of CCHH zinc fingers. Mutagenesis and pH titration experiments indicate that H26 is not involved in metal binding and that Z7 has a tridentate metal-binding site comprised of only residues C3, C6, and H19. By contrast, an F23H mutation that introduces a histidine in the consensus position forms a tetradentate ligand. The structure of the WT Z7 is stable causing restricted ring-flipping of phenyalanines 10 and 23. Dynamics are increased with either the H26A or F23H substitutions and aromatic ring rotation is no longer hindered in the two mutants. The mutations have only small effects on the Kd values for Zn2+ and Co2+ and retain the high thermal stability of the WT domain above 80 °C. Like two previously reported designed zinc fingers with the last ligand replaced by water, the WT Z7 domain is catalytically active, hydrolyzing 4-nitophenyl acetate. We discuss the implications of naturally occurring tridentate zinc fingers for cancer mutations and drug targeting of notoriously undruggable transcription factors. Our findings that Z7 can fold with only a subset of three metal ligands suggests the recent view that most everything about protein structure can be predicted through homology modeling might be premature for at least the resilient and versatile zinc-finger motif.

Stretching Peptides to Generate Small Molecule β-Strand Mimics

Advances in the modulation of protein-protein interactions (PPIs) enable both characterization of PPI networks that govern diseases and design of therapeutics and probes. The shallow protein surfaces that dominate PPIs are challenging to target using standard methods, and approaches for accessing extended backbone structures are limited. Here, we incorporate a rigid, linear, diyne brace between side chains at the i to i+2 positions to generate a family of low-molecular-weight, extended-backbone peptide macrocycles. NMR and density functional theory studies show that these stretched peptides adopt stable, rigid conformations in solution and can be tuned to explore extended peptide conformational space. The diyne brace is formed in excellent conversions (>95%) and amenable to high-throughput synthesis. The minimalist structure-inducing tripeptide core (<300 Da) is amenable to further synthetic elaboration. Diyne-braced inhibitors of bacterial type 1 signal peptidase demonstrate the utility of the technique.

NMR structure of the Bacillus cereus hemolysin II C-terminal domain reveals a novel fold

In addition to multiple virulence factors, Bacillus cereus a pathogen that causes food poisoning and life-threatening wound infections, secretes the pore-forming toxin hemolysin II (HlyII). The HlyII toxin has a unique 94 amino acid C-terminal domain (HlyIIC). HlyIIC exhibits splitting of NMR resonances due to cis/trans isomerization of a single proline near the C-terminus. To overcome heterogeneity, we solved the structure of P405M-HlyIIC, a mutant that exclusively stabilizes the trans state. The NMR structure of HlyIIC reveals a novel fold, consisting of two subdomains αA-β1-β2 and β3-β4-αB-β5, that come together in a barrel-like structure. The barrel core is fastened by three layers of hydrophobic residues. The barrel end opposite the HlyIIC-core has a positively charged surface, that by binding negatively charged moieties on cellular membranes, may play a role in target-cell surface recognition or stabilization of the heptameric pore complex. In the WT domain, dynamic flexibility occurs at the N-terminus and the first α-helix that connects the HlyIIC domain to the HlyII-core structure. In the destabilizing P405M mutant, increased flexibility is evident throughout the first subdomain, suggesting that the HlyIIC structure may have arisen through gene fusion.

NMR proton chemical shift prediction of C·C mismatches in B-DNA

Accurate prediction of DNA chemical shifts facilitates resonance assignment and allows recognition of different conformational features. Based on the nearest neighbor model and base pair replacement approach, we have determined a set of triplet chemical shift values and correction factors for predicting the proton chemical shifts of B-DNA containing an internal C·C mismatch. Our results provide a reliable chemical shift prediction with an accuracy of 0.07 ppm for non-labile protons and 0.09 ppm for labile protons. In addition, we have also shown that the correction factors for C·C mismatches can be used interchangeably with those for T·T mismatches. As a result, we have generalized a set of correction factors for predicting the flanking residue chemical shifts of pyrimidine·pyrimidine mismatches.

NMR assignments for the cis and trans forms of the hemolysin II C-terminal domain

Pathogenic bacteria secrete pore-forming toxins (PFTs) to selectively defend against immune cells and to break through cellular barriers in the host. Understanding how PFTs attack cell membranes is not only essential for therapeutic intervention but for designing agents to deliver drugs to specific human cell subtypes, for example in anti-cancer or anti-viral therapies. Many toxins contain accessory domains that help recognize specific molecular epitopes on the membranes of target cells, including proteins, carbohydrates, and lipids. Here we report NMR assignments for the 94-residue 10 kDa C-terminal accessory domain of Bacillus cereus hemolysin II, HlyIIC, that has no known structural or functional homologues. The HlyIIC domain exists in a dynamic equilibrium due to cis/trans isomerization of its Gly86-Pro87 peptide bond. The cis and trans forms are about equally populated and are in slow exchange on the NMR timescale, giving rise to separate signals for approximately half of the residues in the domain. Assignments for the cis and trans forms were achieved with the aid of a P87M mutant that stabilizes the trans form, and separate sequential walks for the two forms in 3D NMR spectra of the wild-type HlyIIC. Based on backbone chemical shifts, the domain has a α1-α2-β1-β2-β3-β4-α3-β5 order of secondary structure elements. The last strand in the trans form and in the P87M mutant is shortened near Pro87 compared to the cis form. Both cis/trans isomerization and the P87M mutation cause large chemical shift changes throughout HlyIIC, suggesting that the proline is important in stabilizing the structure of the domain. The NMR assignments pave the way for solving the structures of the multiple conformational forms of HlyIIC and establishing their mechanism of interconversion.

NMR assignments for the telokin-like domain of bacteriophage P22 coat protein

The bacteriophage P22 virion is assembled from identical coat protein monomers in a complex reaction that is generally conserved among tailed, double-stranded DNA bacteriophages and viruses. Many coat proteins of dsDNA viruses have structures based on the HK97 fold, but in some viruses and phages there are additional domains. In the P22 coat protein, a "telokin-like" domain was recently identified, whose structure has not yet been characterized at high-resolution. Two recently published low-resolution cryo-EM reconstructions suggest markedly different folds for the telokin-like domain that lead to alternative conclusions about its function in capsid assembly and stability. Here we report (1)H, (15)N, and (13)C NMR resonance assignments for the telokin-like domain. The secondary structure predicted from the chemical shift values obtained in this work shows significant discrepancies from both cryo-EM models but agrees better with one of the models. In particular, the functionally important "D-loop" in one model shows chemical shifts and solvent exchange protection more consistent with β-sheet structure. Our work will set the basis for a high-resolution NMR structure determination of the telokin-like domain that will help improve the cryo-EM models, and in turn lead to a better understanding of how coat protein monomers assemble into the icosahedral capsids required for virulence.

Quantum Mechanical Study of Vicinal J Spin-Spin Coupling Constants for the Protein Backbone

We have performed densisty functional theory (DFT) calculations of vicinal J coupling constants involving the backbone torsional angle for the protein GB3 using our recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach (Xiao He et al. J. Phys. Chem. B 2009, 113, 10380-10388). Interestingly, the calculated values based on an NMR structure are more accurate than those based on a high-resolution X-ray strucure because the NMR structure was refined using a large number of residual dipolar couplings (RDCs) whereas the hydrogen atoms were added into the X-ray structure in idealized positions, confirming that the postioning of the hydrogen atoms relative to the backbone atoms is important to the accuracy of J coupling constant prediction. By comparing three Karplus equations, our results have demonstrated that hydrogen bonding, substituent and electrostatic effects could have significant impacts on vicinal J couplings even though they depend mostly on the intervening dihedral angles. The root-mean-square deviations (RMSDs) of the calculated (3)J(H(N),H(α)), (3)J(H(N),C(β)), (3)J(H(N),C') values based on the NMR structure are 0.52, 0.25, and 0.35 Hz, respectively, after taking the dynamic effect into consideration. The excellent accuracy demonstrates that our AF-QM/MM approach is a useful tool to study the relationship between J coupling constants and the structure and dynamics of proteins.

NMR proton chemical shift prediction of T·T mismatches in B-DNA duplexes

A proton chemical shift prediction scheme for B-DNA duplexes containing a T·T mismatch has been established. The scheme employs a set of T·T mismatch triplet chemical shift values, 5'- and 3'-correction factors extracted from reference sequences, and also the B-DNA chemical shift values predicted by Altona et al. The prediction scheme was tested by eight B-DNA duplexes containing T·T mismatches. Based on 560 sets of predicted and experimental proton chemical shift values, the overall prediction accuracy for non-labile protons was determined to be 0.07 ppm with an excellent correlation coefficient of 0.9996. In addition, the prediction accuracy for 96 sets of labile protons was found to be 0.22 ppm with a correlation coefficient of 0.9961. The prediction scheme developed herein can facilitate resonance assignments of B-DNA duplexes containing T·T mismatches and be generalized for the chemical shift prediction of other DNA mismatches. Our chemical shift data will also be useful for establishing structure-chemical shift information in B-DNA containing mismatches.

13C structuring shifts for the analysis of model β-hairpins and β-sheets in proteins: diagnostic shifts appear only at the cross-strand H-bonded residues

The present studies have shown that (13)C=O, (13)C(α) and (13)C(β) of H-bonded strand residues in β-hairpins provide additional probes for quantitating the extent of folding in β-hairpins and other β-sheet models. Large differences in the structuring shifts (CSDs) of these (13)C sites in H-bonded versus non-H-bonded sites are observed: the differences between H-bonded and non-H-bonded sites are greater than 1.2 ppm for all three (13)C probes. This prompts us to suggest that efforts to determine the extent of hairpin folding from (13)C shifts should be based exclusively on the observation at the cross-strand H-bonded sites. Furthermore, the statistics suggest the (13)C' and (13)C(β) CSDs will provide the best differentiation with 100%-folded CSD values approaching -2.6 and +3 ppm, respectively, for the H-bonded sites. These conclusions can be extended to edge-strands of protein β-sheets. Our survey of reported (13)C shifts in β-proteins indicates that some of the currently employed random coil values need to be adjusted, particularly for ionization-induced effects.

Structural characterization by NMR of a double phosphorylated chimeric peptide vaccine for treatment of Alzheimer's disease

Rational design of peptide vaccines becomes important for the treatment of some diseases such as Alzheimer's disease (AD) and related disorders. In this study, as part of a larger effort to explore correlations of structure and activity, we attempt to characterize the doubly phosphorylated chimeric peptide vaccine targeting a hyperphosphorylated epitope of the Tau protein. The 28-mer linear chimeric peptide consists of the double phosphorylated B cell epitope Tau₂₂₉₋₂₃₇[pThr231/pSer235] and the immunomodulatory T cell epitope Ag85B₂₄₁₋₂₅₅ originating from the well-known antigen Ag85B of the Mycobacterium tuberculosis, linked by a four amino acid sequence -GPSL-. NMR chemical shift analysis of our construct demonstrated that the synthesized peptide is essentially unfolded with a tendency to form a β-turn due to the linker. In conclusion, the -GPSL- unit presumably connects the two parts of the vaccine without transferring any structural information from one part to the other. Therefore, the double phosphorylated epitope of the Tau peptide is flexible and accessible.

High-resolution characterization of intrinsic disorder in proteins: expanding the suite of (13)C-detected NMR spectroscopy experiments to determine key observables

Order in disorder: The characterization of intrinsically disordered proteins by NMR spectroscopy is a necessity on the one hand and a continuous challenge on the other. We propose two experiments that provide diagnostic parameters to monitor the degree of unfolding of a polypeptide. The test was performed on the yeast Cox17 protein, known to gain its function through maturation from an intrinsically disordered state (see figure).

Nuclear magnetic resonance structure and dynamics of the response regulator Sma0114 from Sinorhizobium meliloti

Receiver domains control intracellular responses triggered by signal transduction in bacterial two-component systems. Here, we report the solution nuclear magnetic resonance structure and dynamics of Sma0114 from the bacterium Sinorhizobium meliloti, the first such characterization of a receiver domain from the HWE-kinase family of two-component systems. The structure of Sma0114 adopts a prototypical α(5)/β(5) Rossman fold but has features that set it apart from other receiver domains. The fourth β-strand of Sma0114 houses a PFxFATGY sequence motif, common to many HWE-kinase-associated receiver domains. This sequence motif in Sma0114 may substitute for the conserved Y-T coupling mechanism, which propagates conformational transitions in the 455 (α4-β5-α5) faces of receiver domains, to prime them for binding downstream effectors once they become activated by phosphorylation. In addition, the fourth α-helix of the consensus 455 face in Sma0114 is replaced with a segment that shows high flexibility on the pico- to nanosecond time scale by (15)N relaxation data. Secondary structure prediction analysis suggests that the absence of helix α4 may be a conserved property of the HWE-kinase-associated family of receiver domains to which Sma0114 belongs. In spite of these differences, Sma0114 has a conserved active site, binds divalent metal ions such as Mg(2+) and Ca(2+) that are required for phosphorylation, and exhibits micro- to millisecond active-site dynamics similar to those of other receiver domains. Taken together, our results suggest that Sma0114 has a conserved active site but differs from typical receiver domains in the structure of the 455 face that is used to effect signal transduction following activation.

NMR assignments for the Sinorhizobium meliloti response regulator Sma0114

Response regulators are terminal ends of bacterial two-component systems that undergo extensive structural reorganization in response to phosphoryl transfer from their cognate histidine kinases. The response regulator encoded by the gene sma0114 of Sinorhizobium meliloti is a part of a unique class of two-component systems that employ HWE histidine kinases. The distinct features of Sma0114 include a PFxFATGY motif that houses the conserved threonine in the "Y-T coupling" conformational switch which mediates output response through downstream protein-protein interactions, and the replacement of the conserved phenylalanine/tyrosine in Y-T coupling by a leucine. Here we present (1)H, (15)N, and (13)C NMR assignments for Sma0114. We identify the secondary structure of the protein based on TALOS chemical shift analysis, (3)J(HNHα) coupling constants and hydrogen-deuterium exchange. The secondary structure determined by NMR is in good agreement with that predicted from the sequence. Both methods suggest that Sma0114 differs from standard CheY-like folds by missing the fourth α-helix. Our initial NMR characterization of Sma0114 paves the way to a full investigation of the structure and dynamics of this response regulator.

β-Sheet 13C structuring shifts appear only at the H-bonded sites of hairpins

The (13)C chemical shifts measured for designed β-hairpins indicate that the structuring shifts (upfield for Cα and C', downfield for Cβ) previously reported as diagnostic for β-structuring in proteins appear only at the H-bonded strand residues. The resulting periodicity of structuring shift magnitudes is not, however, a consequence of H-bonding status; rather, it reflects a previously unrecognized alternation in the backbone torsion angles of β-strands. This feature of hairpins is also likely to be present in proteins. The study provides reference values for the expectation shifts for (13)C sites in β-structures that should prove useful in the characterization of the folding equilibria of β-sheet models.

Dynamic alpha-helix structure of micelle-bound human amylin

Amylin is an endocrine hormone that regulates metabolism. In patients afflicted with type 2 diabetes, amylin is found in fibrillar deposits in the pancreas. Membranes are thought to facilitate the aggregation of amylin, and membrane-bound oligomers may be responsible for the islet beta-cell toxicity that develops during type 2 diabetes. To better understand the structural basis for the interactions between amylin and membranes, we determined the NMR structure of human amylin bound to SDS micelles. The first four residues in the structure are constrained to form a hairpin loop by the single disulfide bond in amylin. The last nine residues near the C terminus are unfolded. The core of the structure is an alpha-helix that runs from about residues 5-28. A distortion or kink near residues 18-22 introduces pliancy in the angle between the N- and C-terminal segments of the alpha-helix. Mobility, as determined by (15)N relaxation experiments, increases from the N to the C terminus and is strongly correlated with the accessibility of the polypeptide to spin probes in the solution phase. The spin probe data suggest that the segment between residues 5 and 17 is positioned within the hydrophobic lipid environment, whereas the amyloidogenic segment between residues 20 and 29 is at the interface between the lipid and solvent. This orientation may direct the aggregation of amylin on membranes, whereas coupling between the two segments may mediate the transition to a toxic structure.

Structural and dynamic characterization of intrinsically disordered human securin by NMR spectroscopy

Understanding the molecular action of securin, the inhibitor of separase in mitosis, is of immense theoretical and biomedical importance. The residue-level structural description of an intrinsically disordered protein of this length (202 amino acids, containing 24 prolines), however, represents a particular challenge. Here we combined (1)H-detected and (13)C-detected protonless NMR experiments to achieve full assignment of securin's backbone amide resonances. Chemical shifts, (15)N relaxation rates (R(1), R(2), (1)H-(15)N NOEs), (1)H exchange rates with the solvent (CLEANEX-PM), and (1)H-(15)N residual dipolar couplings were determined along the entire length of the protein. This analysis showed that securin is not entirely disordered, but segregates into a largely disordered N-terminal half and a C-terminal half with transient segmental order, within which the segment D(150)-F(159) has a significant helical tendency and segments E(113)-S(127) and W(174)-L(178) also show a significant deviation from random-coil behavior. These results, in combination with bioinformatic and biochemical data on the securin/separase interaction, shed light on the inhibitory action of securin on separase.

Crystal polymorphism of protein GB1 examined by solid-state NMR spectroscopy and X-ray diffraction

The study of micro- or nanocrystalline proteins by magic-angle spinning (MAS) solid-state NMR (SSNMR) gives atomic-resolution insight into structure in cases when single crystals cannot be obtained for diffraction studies. Subtle differences in the local chemical environment around the protein, including the characteristics of the cosolvent and the buffer, determine whether a protein will form single crystals. The impact of these small changes in formulation is also evident in the SSNMR spectra; however, the changes lead only to correspondingly subtle changes in the spectra. Here, we demonstrate that several formulations of GB1 microcrystals yield very high quality SSNMR spectra, although only a subset of conditions enable growth of single crystals. We have characterized these polymorphs by X-ray powder diffraction and assigned the SSNMR spectra. Assignments of the 13C and 15N SSNMR chemical shifts confirm that the backbone structure is conserved, indicative of a common protein fold, but side chain chemical shifts are changed on the surface of the protein, in a manner dependent upon crystal packing and electrostatic interactions with salt in the mother liquor. Our results demonstrate the ability of SSNMR to reveal minor structural differences among crystal polymorphs. This ability has potential practical utility for studying the formulation chemistry of industrial and therapeutic proteins, as well as for deriving fundamental insights into the phenomenon of single-crystal growth.

13C-detection in RNA bases: revealing structure-chemical shift relationships

The chemical shifts of the unprotonated carbons in the proton-deficient nucleobases of RNA are rarely reported, despite the valuable information that they contain about base-pairing and base-stacking. We have developed 13C-detected 2D-experiments to identify the unprotonated 13C in the RNA bases and have assigned all the base nuclei of uniformly 13C,15N-labeled HIV-2 TAR-RNA. The 13C chemical shift distributions revealed perturbations correlated with the base-pairing and base-stacking properties of all four base-types. From this work, we conclude that the information contained in the chemical shift perturbations within the base rings can provide valuable restraint information for solving RNA structures, especially in conformational averaged regions, where NOE-based information is not available.

NMR study and molecular dynamics simulations of optimized β-hairpin fragments of protein G

The stability and structure of several beta-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the beta-turn conformation. Finally, we show that the effects of combining these single mutations are additive, suggesting that independent stabilizing interactions can be isolated and evaluated in a simple model system. Our results indicate that the structure and stability of this beta-hairpin peptide can be modulated in numerous ways and thus contributes toward a more complete understanding of this important model beta-hairpin as well as to the folding and stability of larger peptides and proteins.

Proton chemical shift prediction of A.A mismatches in B-DNA duplexes

A proton chemical shift prediction method has been developed for double helical DNAs containing A.A mismatches. This method makes use of the chemical shift prediction scheme for normal B-DNA duplexes developed by Altona and co-workers and a set of A.A mismatch triplet chemical shift values and corrections factors extracted from reference sequences. The triplet values are used for predicting chemical shifts of A.A mismatches whereas the normal B-DNA chemical shifts and correction factors are used for the flanking residues of A.A mismatches. Both 5'- and 3'-correction factors have been determined from the chemical shift differences upon replacing the A.A mismatch in a duplex with an A.T base pair. Based on 560 sets of predicted and experimental chemical shifts, the overall prediction accuracy for various types of protons has been determined to be 0.07 ppm with an excellent correlation coefficient of 0.9996.

Predicting 13Calpha chemical shifts for validation of protein structures

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

Macrocyclic beta-sheet peptides that mimic protein quaternary structure through intermolecular beta-sheet interactions

This paper reports the design, synthesis, and characterization of a family of cyclic peptides that mimic protein quaternary structure through beta-sheet interactions. These peptides are 54-membered-ring macrocycles comprising an extended heptapeptide beta-strand, two Hao beta-strand mimics [JACS 2000, 122, 7654] joined by one additional alpha-amino acid, and two delta-linked ornithine beta-turn mimics [JACS 2003, 125, 876]. Peptide 3a, as the representative of these cyclic peptides, contains a heptapeptide sequence (TSFTYTS) adapted from the dimerization interface of protein NuG2 [PDB ID: 1mio]. 1H NMR studies of aqueous solutions of peptide 3a show a partially folded monomer in slow exchange with a strongly folded oligomer. NOE studies clearly show that the peptide self-associates through edge-to-edge beta-sheet dimerization. Pulsed-field gradient (PFG) NMR diffusion coefficient measurements and analytical ultracentrifugation (AUC) studies establish that the oligomer is a tetramer. Collectively, these experiments suggest a model in which cyclic peptide 3a oligomerizes to form a dimer of beta-sheet dimers. In this tetrameric beta-sheet sandwich, the macrocyclic peptide 3a is folded to form a beta-sheet, the beta-sheet is dimerized through edge-to-edge interactions, and this dimer is further dimerized through hydrophobic face-to-face interactions involving the Phe and Tyr groups. Further studies of peptides 3b-3n, which are homologues of peptide 3a with 1-6 variations in the heptapeptide sequence, elucidate the importance of the heptapeptide sequence in the folding and oligomerization of this family of cyclic peptides. Studies of peptides 3b-3g show that aromatic residues across from Hao improve folding of the peptide, while studies of peptides 3h-3n indicate that hydrophobic residues at positions R3 and R5 of the heptapeptide sequence are important in oligomerization.

Oxygen as a Paramagnetic Probe of Clustering and Solvent Exposure in Folded and Unfolded States of an SH3 Domain

The N-terminal SH3 domain of the Drosophila modular protein Drk undergoes slow exchange between a folded (Fexch) and highly populated unfolded (Uexch) state under nondenaturing buffer conditions, enabling both Fexch and Uexch states to be simultaneously monitored. The addition of dissolved oxygen, equilibrated to a partial pressure of either 30 atm or 60 atm, provides the means to study solvent exposure with atomic resolution via 13C NMR paramagnetic shifts in 1H,13C HSQC (heteronuclear single quantum coherence) spectra. Absolute differences in these paramagnetic shifts between the Fexch and Uexch states allow the discrimination of regions of the protein which undergo change in solvent exposure upon unfolding. Contact with dissolved oxygen for both the Fexch and Uexch states could also be assessed through 13C paramagnetic shifts which were normalized based on the corresponding paramagnetic shifts seen in the free amino acids. In the Fexch state, the 13C nuclei belonging to the hydrophobic core of the protein exhibited very weak normalized paramagnetic shifts while those with greater solvent accessible surface area exhibited significantly larger normalized shifts. The Uexch state displayed less varied 13C paramagnetic shifts although distinct regions of protection from solvent exposure could be identified by a lack of such shifts. These regions, which included Phe9, Thr12, Ala13, Lys21, Thr22, Ile24, Ile27, and Arg38, overlapped with those found to have residual nativelike and non-native structures in previous studies and in some cases provided novel information. Thus, the paramagnetic shifts from dissolved oxygen are highly useful in the study of a transient structure or clustering in disordered systems, where conventional NMR measurements (couplings, chemical shift deviations from random coil values, and NOEs) may give little information.

Explicit factorization of external coordinates in constrained statistical mechanics models

If a macromolecule is described by curvilinear coordinates or rigid constraints are imposed, the equilibrium probability density that must be sampled in Monte Carlo simulations includes the determinants of different mass-metric tensors. In this work, the authors explicitly write the determinant of the mass-metric tensor G and of the reduced mass-metric tensor g, for any molecule, general internal coordinates and arbitrary constraints, as a product of two functions; one depending only on the external coordinates that describe the overall translation and rotation of the system, and the other only on the internal coordinates. This work extends previous results in the literature, proving with full generality that one may integrate out the external coordinates and perform Monte Carlo simulations in the internal conformational space of macromolecules.

Adsorption Mechanism at the Molecular Level between Polymers and Uremic Octapeptide by the 2D1H NMR Technique

To remove uremic octapeptide from the blood stream of uremic patients, various modified polyacylamide cross-linked absorbents were prepared. Adsorption experiments showed these absorbents have significant differences in adsorption capacity to the target peptide. In this paper, two-dimension proton nuclear magnetic resonance (2D 1H NMR) spectroscopy was used to investigate the interaction mechanism between the peptide and the adsorbents. Because of the insolubility of the absorbent, some soluble linear polymers with the same functional groups as the absorbents were employed as the model adsorbents in 2D 1H NMR. The preferred binding site for the peptide and polymers was identified to be at the C-terminal carboxyl group of the octapeptide via chemical shift perturbation effects. In this study, we found that hydrogen bonding, electrostatic, and hydrophobic interactions all play a role in the interaction force but had different contributions. Especially, the great chemical shift changes of the aromatic amino acid residues (Trp) during the interaction between butyl-modified polyacrylamide and octapeptide suggested the hydrophobic interaction, incorporated with the electrostatic force, played an important role in the binding reaction in aqueous solutions. This information not only rationally explained the results of the adsorption experiments, but also identified the effective binding site and mechanism, and shall provide a structural basis for designing better affinity-type adsorbents for the target peptide.

A simple method to predict protein flexibility using secondary chemical shifts

Protein motions play a critical role in many biological processes, such as enzyme catalysis, allosteric regulation, antigen-antibody interactions, and protein-DNA binding. NMR spectroscopy occupies a unique place among methods for investigating protein dynamics due to its ability to provide site-specific information about protein motions over a large range of time scales. However, most NMR methods require a detailed knowledge of the 3D structure and/or the collection of additional experimental data (NOEs, T1, T2, etc.) to accurately measure protein dynamics. Here we present a simple method based on chemical shift data that allows accurate, quantitative, site-specific mapping of protein backbone mobility without the need of a three-dimensional structure or the collection and analysis of NMR relaxation data. Further, we show that this chemical shift method is able to quantitatively predict per-residue RMSD values (from both MD simulations and NMR structural ensembles) as well as model-free backbone order parameters.

Energetic characterization of short helical polyalanine peptides in water: analysis of 13C=O chemical shift data

Measured at 2 degrees C in water, NMR chemical shifts of (13)C=O labeled central alanine residues of peptides W-Lys(5)-(t)L(3)-Ala(n)-(t)L(3)-Lys(5)NH(2), n = 9, 11, 13, 15, 19 and W-Lys(5)-(t)L(3)-a-Ala(n)-A-Inp-(t)L(2)-Lys(5)NH(2) (a = D-Ala; (t)L = tert-leucine; Inp = 4-carboxypiperidine) are used to assign jt(L) and ct(L), the N- and C-terminal (t)L capping parameters and length-dependent values for w(Ala)(n), the alanine helical propensity for Ala(n) peptides. These parameters allow Lifson-Roig characterization of the stabilities of Ala(n)() helices in water. To facilitate chemical shift characterization, different (13)C/(12)C ratios are incorporated into specific Ala sites to code up to six residue sites per peptide. Large left/right chemical shift anisotropies are intrinsic to helical polyalanines, and a correcting L-R-based model is introduced. Capping parameters jt(L) = ct(L) lie in the range of 0.3 to 0.5; the (t)L residues are thus moderately helix-destabilizing. For helical conformations of lengths shorter than eight residues, assigned values for w(Ala) approach 1.0 but increase monotonically with length to a value of 1.59 for w(Ala)(19).

Bayesian estimation of Karplus parameters and torsion angles from three-bond scalar couplings constants

We apply Bayesian inference to analyze three-bond scalar coupling constants in an objective and consistent way. The Karplus curve and a Gaussian error law are used to model scalar coupling measurements. By applying Bayes' theorem, we obtain a probability distribution for all unknowns, i.e., the torsion angles, the Karplus parameters, and the standard deviation of the Gaussian. We infer all these unknowns from scalar coupling data using Markov chain Monte Carlo sampling and analytically derive a probability distribution that only involves the torsion angles.

Determination of Membrane Protein Structure and Dynamics by Magic-Angle-Spinning Solid-State NMR Spectroscopy†

It is shown that molecular structure and dynamics of a uniformly labeled membrane protein can be studied under magic-angle-spinning conditions. For this purpose, dipolar recoupling experiments are combined with novel through-bond correlation schemes that probe mobile protein segments. These NMR schemes are demonstrated on a uniformly [13C,15N] variant of the 52-residue polypeptide phospholamban. When reconstituted in lipid bilayers, the NMR data are consistent with an alpha-helical trans-membrane segment and a cytoplasmic domain that exhibits a high degree of structural disorder.

Water-solubilized, cap-stabilized, helical polyalanines: calibration standards for NMR and CD analyses

NMR and CD studies are reported for two length series of solubilized, spaced, highly helical polyalanines that are N-capped by the optimal helix stabilizer (beta)Asp-Hel and C-capped by beta-aminoalanine beta and that are studied in water at 2 degrees C, pH 1-8. NMR analysis yields a structural characterization of the peptide Ac(beta)AspHelAla(8)betaNH(2) and selected members of one (beta)AspHelAla(n)beta series. At pH > 4.5 the (beta)AspHel cap provides a preorganized triad of carboxylate anion and two amide residues that is complementary to the helical polyalanine N-terminus. The C-terminal beta-aminoalanine assumes a helix-stabilizing conformation consistent with literature precedents. H(N)CO NMR experiments applied to capped, uniformly (13)C- and (15)N-labeled Ala(8) and Ala(12) peptides define Ala(n) hydrogen bonding signatures as alpha-helical without detectable 3(10) character. Relative NH-->ND exchange rates yield site protection factors PF(i) that define uniquely high fractional helicities FH for the peptide Ala(n) regions. These Ala(n) calibration series, studied in water and lacking helix-stabilizing tertiary structure, yield the first (13)C NMR chemical shifts, (3)J(HNH)(alpha) coupling constants, and CD ellipticities [theta(Molar)](lambda,n) characteristic of a fully helical alanine within an Ala(n) context. CD data are used to assign parameters X and [theta](lambda,infinity), required for rigorous calculation of FH values from CD ellipticities.

GROMACS: Fast, flexible, and free

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

NMR structure of the C-terminal domain of SecA in the free state

SecA is an integral component of the prokaryotic Sec preprotein secretory translocase system. We report here the solution NMR structure of a fragment corresponding to the C-terminal domain of Escherichia coli SecA. In the presence of Zn2+, the fragment adopts a shortened version of the classic betabetaalpha zinc finger fold. The isolated C-terminal domain shows substantial differences from the X-ray structure of a homologous SecA domain bound to the chaperone-like cofactor SecB. The differences between the structures of the free and bound forms suggest that binding to SecB causes a perturbation of the C-terminal domain's intrinsically favored betabetaalpha fold.

Insights into Partially Folded or Unfolded States of Metalloproteins from Nuclear Magnetic Resonance

Nuclear magnetic resonance (NMR) provides detailed insights into the conformational features of unfolded and partially folded proteins. In the case of metalloproteins, special attention should be devoted to the characterization of the properties of the metal binding sites, and specific approaches need to be developed depending on the nature of the metal ion and its coordination environment. At the same time, metal-based NMR parameters may help in getting a better picture of the average structural properties of the metalloprotein. A critical evaluation of the limits of applicability of paramagnetic effects for solution structure determination in partially folded or unfolded proteins is presented. The coupling between NMR characterization of structure and dynamic of the polypeptide chain and of the metal environment provides insights into the stabilizing role of metal ions in metalloproteins. The overall approach is illustrated for some case examples of increasing flexibility obtained far from native conditions for cytochrome c and superoxide dismutase, two metalloproteins that have been extensively studied in our lab and whose misfolded forms may be relevant for important biological processes.

NMR assignments for a helical 40 kDa membrane protein

Backbone nuclear magnetic resonance (NMR) assignments were achieved for diacylglycerol kinase (DAGK) in detergent micelles. DAGK is a homotrimeric integral membrane protein comprised of 121 residue subunits, each having three transmembrane segments. Assignments were made using TROSY-based pulse sequences. DAGK was found to be an almost exclusively helical protein. This work points to the feasibility of both solving the structure of DAGK using solution NMR methods and using NMR as a primary tool in structural studies of other helical integral membrane proteins of similar size and complexity.

NMR Characterization of Native Liquid Spider Dragline Silk from Nephila edulis

Solid spider dragline silk is well-known for its mechanical properties. Nonetheless a detailed picture of the spinning process is lacking. Here we report NMR studies on the liquid silk within the wide sac of the major ampullate (m.a.) gland from the spider Nephila edulis. The resolution in the NMR spectra is shown to be significantly improved by the application of magic-angle spinning (MAS). From the narrow width of the resonance lines and the chemical shifts observed, it is concluded that the silk protein within the wide sac of the m.a. gland is dynamically disordered throughout the molecule in the sense that each amino acid of a given type senses an identical environment, on average. The NMR data obtained are consistent with an isotropic liquid phase.

Mapping the Electronic Structure of the Blue Copper Site in Plastocyanin by NMR Relaxation

The biological function of metalloproteins stems from the electronic and geometric structures of their active sites. Thus, in blue copper proteins such as plastocyanins, an unusual electronic structure of the metal site is believed to contribute to the rapid, long-range electron-transfer reactivity that characterizes these proteins. To clarify this structure-function relationship, numerous quantum chemical calculations of the electronic structure of the blue copper proteins have been made. However, the obtained structures depend strongly on the applied model. Experimental approaches based on ENDOR spectroscopy and X-ray absorption have also been used to elucidate the electronic structure of the blue copper site. Still, the determination of the electronic structure relies on a calibration with quantum chemical calculations, performed on small model complexes. Here we present an approach that allows a direct experimental mapping of the electron spin delocalization in paramagnetic metalloproteins using oxidized plastocyanin from Anabaena variabilis as an example. The approach utilizes the longitudinal paramagnetic relaxation of protons close to the metal site and relies on the dependence of these relaxations on the spatial distribution of the unpaired electron of the metal ion. Surprisingly it is found that the unpaired electron of the copper ion in plastocyanin is less delocalized than predicted by most of the quantum chemical calculations.

DSA: a new internal standard for NMR studies in aqueous solution

[structure: see text] The widely used internal standard for NMR studies in aqueous solution DSS (sodium 4,4-dimethyl-4-silapentane-1-sulfonate) can interact with cationic peptides, diminishing its value for such studies. This paper introduces DSA (4,4-dimethyl-4-silapentane-1-ammonium trifluoroacetate) as a new internal standard that does not suffer from this problem.

NMR identification of left-handed polyproline type II helices

NMR characteristics of a model left-handed 3(1)-helical peptide are reported in this study. With temperature and sequence corrections on the predicted random coil (15)N chemical shifts, a significant (15)N chemical shift deviation is observed for the model 3(1) peptide. The (15)N chemical shift differences also correlate well with the molar ellipticities (at 220 nm) of the CD spectra at different temperatures, indicating that the (15)N chemical shift is a sensitive probe for 3(1)-helices. The average (3)J(HNalpha) and (1)J(CalphaHalpha) values of the model peptide are determined to be 6.5 and 142.6 Hz, respectively, which are consistent with the values calculated from the geometry of 3(1)-helices. With careful measurements of amide (15)N chemical shifts and incorporating temperature and sequence effect corrections, the (15)N chemical shifts can be used together with (3)J(HNalpha) and (1)J(CalphaHalpha) to differentiate 3(1)-helices from random coils with high confidence. Based on the observed NMR characteristics, a strategy is developed for probing left-handed 3(1)-helical structures from other secondary structures.

NMR structure of two novel polyethylene glycol conjugates of the human growth hormone-releasing factor, hGRF(1-29)-NH2

Two novel mono-PEGylated derivatives of hGRF(1-29)-NH(2) [human growth hormone-releasing factor, fragment 1-29] have been synthesized by regio-specific conjugation of Lys(12) or Lys(21) to a monomethoxy-PEG(5000) chain (compounds Lys(12)PEG-GRF and Lys(21)PEG-GRF). The PEG moiety has been covalently linked at the amino group of a norleucine residue via a carbamate bond. The Lys(12)PEG-GRF regioisomer was found to be slightly less active in vitro than both the unmodified peptide and Lys(21)PEG-GRF. To assess whether the differences in the biological activity of the PEGylated analogues could be related to conformational rearrangements induced by the PEG moiety, the structure of these PEGylated derivatives has been worked out (TFE solution) by means of NMR spectroscopy and molecular dynamics. Secondary structure shifts, hydrogen/deuterium exchange kinetics, temperature coefficients of amide protons, and NOE-based molecular models point out that hGRF(1-29)-NH(2), Lys(21)PEG-GRF and Lys(12)PEG-GRF share a remarkably similar pattern of secondary structure. All three compounds adopt an alpha-helix conformation which spans the whole length of the molecule, and which becomes increasingly rigid on going from the N-terminus to the C-terminus. Residues Lys(12) and Lys(21) are enclosed in all the compounds considered into well-defined alpha-helical domains, indicating that PEGylation either at Lys(12) or Lys(21) does not alter the tendency of the peptide to adopt a stable alpha-helix conformation, nor does it induce appreciable conformational mobility in the proximity of the PEGylation sites. No significant variation of the amphiphilic organization of the alpha-helix is observed among the three peptides. Therefore, the different biological activities observed for the PEGylated analogues are not due to conformational effects, but are rather due to sterical hindrance effects. The relationship between the biological activitiy of the mono-PEGylated derivatives and sterical hindrance is discussed in terms of the topology of interaction between hGRF(1-29)-NH(2) and its receptor.

An empirical correlation between secondary structure content and averaged chemical shifts in proteins

It is shown that the averaged chemical shift (ACS) of a particular nucleus in the protein backbone empirically correlates well to its secondary structure content (SSC). Chemical shift values of more than 200 proteins obtained from the Biological Magnetic Resonance Bank are used to calculate ACS values, and the SSC is estimated from the corresponding three-dimensional coordinates obtained from the Protein Data Bank. ACS values of (1)H(alpha) show the highest correlation to helical and sheet structure content (correlation coefficient of 0.80 and 0.75, respectively); (1)H(N) exhibits less reliability (0.65 for both sheet and helix), whereas such correlations are poor for the heteronuclei. SSC estimated using this correlation shows a good agreement with the conventional chemical shift index-based approach for a set of proteins that only have chemical shift information but no NMR or x-ray determined three-dimensional structure. These results suggest that even chemical shifts averaged over the entire protein retain significant information about the secondary structure. Thus, the correlation between ACS and SSC can be used to estimate secondary structure content and to monitor large-scale secondary structural changes in protein, as in folding studies.

Conformational study of linear and cyclic peptides corresponding to the 276-284 epitope region of HSV gD-1

The results of conformational analysis of linear and cyclic peptides from the 276SALLEDPVG(284) sequence of glycoprotein D of Herpes simplex virus are presented. The epitope peptides were synthesized by SPPS and on resin cyclization was applied for preparation of cyclic compounds. Circular dichroism spectroscopy, Fourier-transform infrared spectroscopy and nuclear magnetic resonance (NMR) were used to determine of the solution structure of both linear and cyclic peptides. The results indicated that the cyclopeptides containing the core of the epitope (DPVG) as a part of the cycle have more stable beta-turn structure than the linear peptides or the cyclic analogues, where the core motif is not a part of the cycle. NMR study of H-SALLc(EDPVGK)-NH(2) confirm presence of a type I beta-turn structure which includes the DPVG epitope core.

3D TEDOR NMR experiments for the simultaneous measurement of multiple carbon-nitrogen distances in uniformly (13)C,(15)N-labeled solids

We describe three-dimensional magic-angle-spinning NMR experiments for the simultaneous measurement of multiple carbon-nitrogen distances in uniformly (13)C,(15)N-labeled solids. The approaches employ transferred echo double resonance (TEDOR) for (13)C-(15)N coherence transfer and (15)N and (13)C frequency labeling for site-specific resolution, and build on several previous 3D TEDOR techniques. The novel feature of the 3D TEDOR pulse sequences presented here is that they are specifically designed to circumvent the detrimental effects of homonuclear (13)C-(13)C J-couplings on the measurement of weak (13)C-(15)N dipolar couplings. In particular, homonuclear J-couplings lead to two undesirable effects: (i) they generate anti-phase and multiple-quantum (MQ) spin coherences, which lead to spurious cross-peaks and phase-twisted lines in the 2D (15)N-(13)C correlation spectra, and thus degrade the spectral resolution and prohibit the extraction of reliable cross-peak intensities, and (ii) they significantly reduce cross-peak intensities for strongly J-coupled (13)C sites (e.g., CO and C(alpha)). The first experiment employs z-filter periods to suppress the anti-phase and MQ coherences and generates 2D spectra with purely absorptive peaks for all TEDOR mixing times. The second approach uses band-selective (13)C pulses to refocus J-couplings between (13)C spins within the selective pulse bandwidth and (13)C spins outside the bandwidth. The internuclear distances are extracted by using a simple analytical model, which accounts explicitly for multiple spin-spin couplings contributing to cross-peak buildup. The experiments are demonstrated in two U-(13)C,(15)N-labeled peptides, N-acetyl-L-Val-L-Leu (N-ac-VL) and N-formyl-L-Met-L-Leu-L-Phe (N-f-MLF), where 20 and 26 (13)C-(15)N distances up to approximately 5-6 A were measured, respectively. Of the measured distances, 10 in N-ac-VL and 13 in N-f-MLF are greater than 3 A and provide valuable structural constraints.

Comparative structural connectivity spectra analysis (CoSCoSA) models of steroid binding to the corticosteroid binding globulin

A three-dimensional quantitative spectrometric data-activity relationship (3D-QSDAR) model was developed that is built by combining NMR spectral information with structural information in a 3D-connectivity matrix. The 3D-connectivity matrix is built by displaying all possible carbon-to-carbon connections with their assigned carbon NMR chemical shifts and distances between the carbons. Selected 2D (13)C-(13)C COrrelation SpectroscopY (COSY) (through-bond nearest neighbors) and selected theoretical 2D (13)C-(13)C distance connectivity spectral slices from the 3D-connectivity matrix to produce a relationship among the spectral patterns for 30 steroids binding to corticosteroid binding globulin. We call this technique a comparative structural connectivity spectra analysis (CoSCoSA) modeling. A CoSCoSA principal component linear regression model based on the combination of (13)C-(13)C COSY and (13)C-(13)C distance spectra principal components (PCs) had an r(2) of 0.96 and a leave-one-out (LOO) cross-validation q(2) of 0.92. A CoSCoSA parallel distributed artificial neural network (PD-ANN) model based on the combination of (13)C-(13)C COSY and (13)C-(13)C distance spectra had an r(2) of 0.96, a leave-three-out q(3)(2) of 0.78, and a leave-ten-out q(10)(2) of 0.73. CoSCoSA modeling attempts to uniquely combine the quantum mechanics information from the NMR chemical shifts with internal molecular atom-to-atom distances into an accurate modeling technique. The CoSCoSA modeling technique has the flexibility and accuracy to outperform the cross-validated variance q(2) of previously published quantitative structure-activity relationship (QSAR), quantitative spectral data-activity relationship (QSDAR), self-organizing map (SOM), and electrotopological state (E-state) models.

Principles of mucin architecture: structural studies on synthetic glycopeptides bearing clustered mono-, di-, tri-, and hexasaccharide glycodomains

The structural characteristics of a mucin glycopeptide motif derived from the N-terminal fragment STTAV of the cell surface glycoprotein CD43 have been investigated by NMR. In this study, a series of molecules prepared by total synthesis were examined, consisting of the peptide itself, three glycopeptides having clustered sites of alpha-O-glycosylation on the serine and threonine side chains with the Tn, TF, and STF carbohydrate antigens, respectively, and one with the beta-O-linked TF antigen. Additionally, a glycopeptide having the sequence SSSAVAV, triglycosylated with the Le(y) epitope, was investigated. NMR data for the tri-STF-STTAV glycopeptide were used to solve the structure of this construct through restrained molecular dynamics calculations. The calculations revealed a defined conformation for the glycopeptide core rooted in the interaction of the peptide and the first N-acetylgalactosamine residue. The similarity of the NMR data for each of the alpha-O-linked glycopeptides demonstrates that this structure persists for each construct and that the mode of attachment of the first sugar and the peptide is paramount in establishing the organization of the core. The core provides a common framework on which a variety of glycans may be displayed. Remarkably, while there is a profound organizational effect on the peptide backbone with the alpha-linked glycans, attachment via a beta-linkage has little apparent consequence.

Stabilization of the helical structure of Y2-selective analogues of neuropeptide Y by lactam bridges

The importance of helical structure in an analogue of NPY selective for the Y2 receptor, Ac[Leu28,31]NPY24-36, has been investigated by introducing a lactam bridge between positions 28 and 32. The resulting analogue, Ac-cyclo28/32[Ala24,Lys28,Leu31,Glu32]NPY24-36, is a potent Y2-selective agonist. Structural analysis by NMR shows that this analogue forms a helical structure in a 40% trifluoroethanol/water mixture, whereas in water only the region around the lactam bridge (Lys28-Glu32) adopts helical-like structure, with both N- and C-termini being poorly defined. The observation of well-defined helical structure in aqueous TFE contrasts with that reported for a similar analogue, Ac-cyclo28/32[Lys28,Glu32]NPY25-36 (Rist et al. FEBS Lett. 1996, 394, 169-173), which consisted of a hairpin-like structure that brought the N- and C-termini into proximity. We have therefore determined the structures of this analogue, as well as those of Ac-cyclo28/32[Ala24,Lys28,Leu31,Glu32]NPY24-36 and Ac-cyclo28/32[Ala24,Lys28,Glu32]NPY24-36, under identical solution conditions (30% TFE/H2O mixture at 308 K) and find essentially the same helical structure in all three peptides. These findings support the proposal that these Y2-selective analogues adopt a helical structure when bound to the Y2 receptor.

Molecular modelling in structural biology

Molecular modelling is a powerful methodology for analysing the three dimensional structure of biological macromolecules. There are many ways in which molecular modelling methods have been used to address problems in structural biology. It is not widely appreciated that modelling methods are often an integral component of structure determination by NMR spectroscopy and X-ray crystallography. In this review we consider some of the numerous ways in which modelling can be used to interpret and rationalise experimental data and in constructing hypotheses that can be tested by experiment. Genome sequencing projects are producing a vast wealth of data describing the protein coding regions of the genome under study. However, only a minority of the protein sequences thus identified will have a clear sequence homology to a known protein. In such cases valuable three-dimensional models of the protein coding sequence can be constructed by homology modelling methods. Threading methods, which used specialised schemes to relate protein sequences to a library of known structures, have been shown to be able to identify the likely protein fold even in cases where there is no clear sequence homology. The number of protein sequences that cannot be assigned to a structural class by homology or threading methods, simply because they belong to a previously unidentified protein folding class, will decrease in the future as collaborative efforts in systematic structure determination begin to develop. For this reason, modelling methods are likely to become increasingly useful in the near future. The role of the blind prediction contests, such as the Critical Assessment of techniques for protein Structure Prediction (CASP), will be briefly discussed. Methods for modelling protein-ligand and protein-protein complexes are also described and examples of their applications given.

Solution structure of Pisum sativum defensin 1 by high resolution NMR: plant defensins, identical backbone with different mechanisms of action

Pisum sativum defensin 1 (Psd1) is a 46 amino acid residue plant defensin isolated from seeds of pea. The three-dimensional structure in solution of Psd1 was determined by two-dimensional NMR data recorded at 600 MHz. Experimental restraints were used for structure calculation using CNS and torsion-angle molecular dynamics. The 20 lowest energy structures were selected and further subjected to minimization, giving a root-mean-square deviation of 0.78(+/- 0.22) A in the backbone and 1.91(+/-0.60) A for over all atoms of the molecule. The protein has a globular fold with a triple-stranded antiparalell beta-sheet and an alpha-helix (from residue Asn17 to Leu27). Psd1 presents the so called "cysteine stabilized alpha/beta motif" and presents identical three-dimensional topology in the backbone with other defensins and neurotoxins. Comparison of the electrostatic surface potential among proteins with high three-dimensional (selected using the softwares TOP and DALI) topology gave insights into the mode of action of Psd1. The surface topologies between proteins that present antifungal activity or sodium channel inhibiting activity are different. On the other hand the surface topology presents several common features with potassium channel inhibitors, suggesting that Psd1 presents this activity. Other common features with potassium channel inhibitors were found including the presence of a lysine residue essential for inhibitory activity. The identity of Psd1 in primary sequence is not enough to infer a mechanism of action, in contrast with the strategy proposed here.