Search

The calculation of DFT (density functional theory) chemical shifts have become an important technique for the verification of a proposed structure. An easily calculated metric developed for proton and carbon chemical shifts of natural products and organic compounds, the calculated chemical shift index (CCSI), has been developed, which uses the deviation of each pair of calculated and experimental chemical shifts. The mean absolute deviation (MAD), which is commonly used as the goodness of fit metric for DFT calculated chemical shifts, can conceal large deviations in the calculated data. A classification strategy is also proposed for the CCSI to highlight when further assessment of the NMR data is required.

In the present study, for the first time 1-amide 4-sulfonamide-1,2,3-triazole scaffolds were synthesized by using an azide-alkyne Huisgen cycloaddition reaction. The target products were obtained in moderate to good yields (45–75%) by using catalytic CuI and green system H2O/EtOH. The easy availability of the inexpensive starting materials, avoiding isolation and handling of hazardous organic azides and mild reaction conditions make this method a valuable tool for generating functionalized 1,2,3-triazole derivatives. The unambiguous characterization of synthesized compounds was accomplished by using various spectroscopic techniques such as 1H NMR, 13C NMR, and FT-IR. The information regarding optimized geometry, were obtained by applying DFT/B3LYP-6-31G(d) method. The electrophilicity index, 1H and 13C chemical shift values, lithium and sodium ion affinities of the desired product 3b have been also calculated by the mentioned method. As a whole, the calculated results were found in close agreement to that of experimental data. The studies revealed that the compound 3b possesses good Li+ and Na+ affinity and cation π interaction plays a vital role in the complexation of 3b. For the first time, nucleus–independent chemical shift index was used to confirm the cation π interaction of 3b.

SARS-CoV-2 is the novel coronavirus that is the causative agent of COVID-19, a sometimes-lethal respiratory infection responsible for a world-wide pandemic. The envelope (E) protein, one of four structural proteins encoded in the viral genome, is a 75-residue integral membrane protein whose transmembrane domain exhibits ion channel activity and whose cytoplasmic domain participates in protein-protein interactions. These activities contribute to several aspects of the viral replication-cycle, including virion assembly, budding, release, and pathogenesis. Here, we describe the structure and dynamics of full-length SARS-CoV-2 E protein in hexadecylphosphocholine micelles by NMR spectroscopy. We also characterized its interactions with four putative ion channel inhibitors. The chemical shift index and dipolar wave plots establish that E protein consists of a long transmembrane helix (residues 8–43) and a short cytoplasmic helix (residues 53–60) connected by a complex linker that exhibits some internal mobility. The conformations of the N-terminal transmembrane domain and the C-terminal cytoplasmic domain are unaffected by truncation from the intact protein. The chemical shift perturbations of E protein spectra induced by the addition of the inhibitors demonstrate that the N-terminal region (residues 6–18) is the principal binding site. The binding affinity of the inhibitors to E protein in micelles correlates with their antiviral potency in Vero E6 cells: HMA ≈ EIPA > DMA >> Amiloride, suggesting that bulky hydrophobic groups in the 5’ position of the amiloride pyrazine ring play essential roles in binding to E protein and in antiviral activity. An N15A mutation increased the production of virus-like particles, induced significant chemical shift changes from residues in the inhibitor binding site, and abolished HMA binding, suggesting that Asn15 plays a key role in maintaining the protein conformation near the binding site. These studies provide the foundation for complete structure determination of E protein and for structure-based drug discovery targeting this protein., Author summary The novel coronavirus SARS-CoV-2, the causative agent of the world-wide pandemic of COVID-19, has become one of the greatest threats to human health. While rapid progress has been made in the development of vaccines, drug discovery has lagged, partly due to the lack of atomic-resolution structures of the free and drug-bound forms of the viral proteins. The SARS-CoV-2 envelope (E) protein, with its multiple activities that contribute to viral replication, is widely regarded as a potential target for COVID-19 treatment. As structural information is essential for drug discovery, we established an efficient sample preparation system for biochemical and structural studies of intact full-length SARS-CoV-2 E protein and characterized its structure and dynamics. We also characterized the interactions of amilorides with specific E protein residues and correlated this with their antiviral activity during viral replication. The binding affinity of the amilorides to E protein correlated with their antiviral potency, suggesting that E protein is indeed the likely target of their antiviral activity. We found that residue asparagine15 plays an important role in maintaining the conformation of the amiloride binding site, providing molecular guidance for the design of inhibitors targeting E protein.

Human islet amyloid polypeptide (hIAPP) aggregates into fibrils through oligomers that have been postulated to contain α-helices as well as β-sheets., Human islet amyloid polypeptide (hIAPP) aggregates into fibrils through oligomers that have been postulated to contain α-helices as well as β-sheets. We employ a site-specific isotope labeling strategy that is capable of detecting changes in dihedral angles when used in conjunction with 2D IR spectroscopy. The method is analogous to the chemical shift index used in NMR spectroscopy for assigning protein secondary structure. We introduce isotope labels at two neighbouring residues, which results in an increased intensity and positive frequency shift if those residues are α-helical versus a negative frequency shift in β-sheets and turns. The 2D IR dihedral index approach is demonstrated for hIAPP in micelles for which the polypeptide structure is known, using pairs of 13C18O isotope labels L12A13 and L16V17, along with single labeled control experiments. Applying the approach to aggregation experiments performed in buffer, we show that about 27–38% of hIAPP peptides adopt an α-helix secondary structure in the monomeric state at L12A13, prior to aggregation, but not at L16V17 residues. At L16V17, the kinetics are described solely by the monomer and fiber conformations, but at L12A13 the kinetics exhibit a third state that is created by an oligomeric intermediate. Control experiments performed with a single isotope label at A13 exhibit two-state kinetics, indicating that a previously unknown change in dihedral angle occurs at L12A13 as hIAPP transitions from the intermediate to fiber structures. We propose a mechanism for aggregation, in which helices seed oligomer formation via structures analogous to leucine rich repeat proteins.

α-Conotoxins preferentially antagonize muscle and neuronal nicotinic acetylcholine receptors (nAChRs). Native α-conotoxins have two disulfide links, CI–CIII and CII–CIV, and owing to the inherent properties of disulfide bonds, α-conotoxins have been systematically engineered to improve their chemical and biological properties. In this study, we explored the possibility of simplifying the disulfide framework of α-conotoxins Vc1.1, BuIA, ImI, and AuIB, by introducing [C2H,C8F] modification to the CI–CIII bond. We therefore explored the possibility of using hydrophobic packing of standard amino acid side chains to replace disulfide bonds as an alternative strategy to nonnatural amino acid cross-links. The impact of CI–CIII disulfide bond replacement on the conformation of the α-conotoxins was investigated using molecular dynamics (MD) simulations and nuclear magnetic resonance chemical shift index study. Two-electrode voltage clamp techniques and MD simulations were used to study the impact of disulfide bond...

Squamous cell carcinoma antigen recognized by T-cells 3 (SART3) is an essential recycling factor in pre-mRNA splicing, which is required for association of U4/U6 small nuclear ribonucleoprotein (snRNP). SART3 contains two RNA recognition motifs (RRMs), and they are responsible for the tertiary interaction with U6 small nuclear RNA. Despite the importance of structural studies for understanding complicate U4/U6 snRNP recycling mechanism, only a few of them have been performed for SART3. Here, the structure of SART3 RRM2 was characterized by heteronuclear multi-dimensional nuclear magnetic resonance experiments. Nearly complete 1H, 15N, and 13C chemical shifts of the backbone residues of RRM2 were assigned. In addition, the secondary structure of RRM2 were predicted by the chemical shift index and TALOS+ analyses, and the results showed that RRM2 forms a “β1-α1-β2-β3-α2-β4-β5” structure, where β4 is not common in the canonical RRM domain structures. Our results will provide structural basis for investigation of SART3-mediated U4/U6 snRNP complex formation.

The first report on crystal and molecular structure of 3,6-diiodo-9-ethyl-9H-carbazole is presented. Experimental room-temperature X-ray and 13C chemical shift studies were supported by advanced theoretical calculations using density functional theory (DFT). The 13C nuclear magnetic shieldings were predicted at the non-relativistic and relativistic level of theory using the zeroth-order regular approximation (ZORA). Theoretical relativistic calculations of chemical shifts of carbons C3 and C6, directly bonded to iodine atoms, produced a reasonable agreement with experiment (initial deviation from experiment of 44.3 dropped to 4.25 ppm). The changes in ring aromatic character via simple harmonic oscillator model of aromaticity (HOMA) and the nucleus independent chemical shift indexes (NICS) calculations were estimated. A good linear correlation between experimental and theoretically predicted structural and NMR parameters was observed.

Predicting secondary structure of protein through assigned backbone chemical shifts has been used widely because of its convenience and flexibility. In spite of its usefulness, chemical shift based analysis has some defects including isotopic shifts and solvent interaction. Here, it is shown that corrected chemical shift analysis for secondary structure of protein. It is included chemical shift correction through consideration of deuterium isotopic effect and calculate chemical shift index using probability-based methods. Enhanced method was applied successfully to one of the proteins from Mycobacterium tuberculosis. It is suggested that correction of chemical shift analysis could increase accuracy of secondary structure prediction of protein and small molecule in solution.

Conversion of Hydrogenobacter thermophilus cytochrome c(552) into a b-type cytochrome by mutagenesis of both heme-binding cysteines to alanines significantly reduces the stability of the protein (Tomlinson, E. J., and Ferguson, S. J. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 5156-5160). To understand the effects of this change on the structure and dynamics of the protein, hetero-nuclear (15)N-edited NMR techniques have been used to characterize this b-type variant. The backbone (15)N, (1)H(N), and (1)H(alpha), and (1)H(beta) resonances of the protein have been assigned. Analysis of (3)J(HN)alpha coupling constants, nuclear Overhauser enhancement intensities, and chemical shift index data demonstrates that the four alpha-helices present in the wild-type protein are retained in the b-type variant. Comparison of the chemical shifts for the b-type and wild-type proteins indicates that the tertiary structures of the two proteins are closely similar. Some subtle differences are, however, observed for residues in the N-terminal region and in the vicinity of the heme-binding pocket. Hydrogen exchange studies show that there are 25 backbone amide protons that exchange very slowly in the b-type variant and confirm that the fluctuations within the b-type protein are of a similar extent to those in the wild-type protein. These data demonstrate the notable retention of the native secondary structure and tertiary fold despite the absence of covalent linkages between the heme group and the protein.

The Bas-Congo virus (BASV) is the first rhabdovirus associated with a human outbreak of acute hemorrhagic fever. The single-stranded, negative-sense RNA genome of BASV contains the five core genes present in all rhabdoviral genomes plus an additional three genes, annotated U1, U2, and U3, with weak (

A common obstacle to NMR studies of proteins is sample preparation. In many cases, proteins targeted for NMR studies are poorly expressed and/or expressed in insoluble forms. Here, we describe a novel approach to overcome these problems. In the protein S tag-intein (PSTI) technology, two tandem 92-residue N-terminal domains of protein S (PrS(2)) from Myxococcus xanthus is fused at the N-terminal end of a protein to enhance its expression and solubility. Using intein technology, the isotope-labeled PrS(2)-tag is replaced with non-isotope labeled PrS(2)-tag, silencing the NMR signals from PrS(2)-tag in isotope-filtered (1)H-detected NMR experiments. This method was applied to the E. coli ribosome binding factor A (RbfA), which aggregates and precipitates in the absence of a solubilization tag unless the C-terminal 25-residue segment is deleted (RbfAΔ25). Using the PrS(2)-tag, full-length well-behaved RbfA samples could be successfully prepared for NMR studies. PrS(2) (non-labeled)-tagged RbfA (isotope-labeled) was produced with the use of the intein approach. The well-resolved TROSY-HSQC spectrum of full-length PrS(2)-tagged RbfA superimposes with the TROSY-HSQC spectrum of RbfAΔ25, indicating that PrS(2)-tag does not affect the structure of the protein to which it is fused. Using a smaller PrS-tag, consisting of a single N-terminal domain of protein S, triple resonance experiments were performed, and most of the backbone (1)H, (15)N and (13)C resonance assignments for full-length E. coli RbfA were determined. Analysis of these chemical shift data with the Chemical Shift Index and heteronuclear (1)H-(15)N NOE measurements reveal the dynamic nature of the C-terminal segment of the full-length RbfA protein, which could not be inferred using the truncated RbfAΔ25 construct. CS-Rosetta calculations also demonstrate that the core structure of full-length RbfA is similar to that of the RbfAΔ25 construct.

A cDNA encoding a cysteine protease inhibitor, cystatin was cloned from sesame (Sesamum indicum L.) seed. This clone was constructed into an expression vector and expressed in E. coli and purified to homogeneous. The recombinant sesame cystatin (SiCYS) showed effectively inhibitory activity toward C1 cysteine proteases. In order to unravel its inhibitory action from structural point of view, multidimensional heteronuclear NMR techniques were used to characterize the structure of SiCYS. The full (1)H, (15)N, and (13)C resonances of SiCYS were assigned. The secondary structure of SiCYS was identified by using the assigned chemical shifts of (1)H(α), (13)C(α), (13)C(β), and (13)CO through the consensus chemical shift index (CSI). The results of CSI analysis of SiCYS suggest eight β-strands (residues 33-46, 51-61, 63-75, 80-87, 150-155, 157-169, 172-183, and 192-195) and two α-helices (residues 16-30, and 120-135).

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., (13)C-(13)C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (-0.75) commensurate to the control (-0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.

The 18.5-kDa myelin basic protein (MBP), the most abundant isoform in human adult myelin, is a multifunctional, intrinsically disordered protein that maintains compact assembly of the sheath. Solution NMR spectroscopy and a hydrophobic moment analysis of MBP's amino-acid sequence have previously revealed three regions with high propensity to form strongly amphipathic α-helices. These regions, located in the central, N- and C-terminal parts of the protein, have been shown to play a role in the interactions of MBP with cytoskeletal proteins, Src homology 3-domain-containing proteins, Ca(2+)-activated calmodulin (Ca(2+)-CaM), and myelin-mimetic membrane bilayers. Here, we have further characterized the structure-function relationship of these three domains. We constructed three recombinant peptides derived from the 18.5-kDa murine MBP: (A22-K56), (S72-S107), and (S133-S159) (which are denoted α1, α2, and α3, respectively). We used a variety of biophysical methods (circular dichroism spectroscopy, isothermal titration calorimetry, transmission electron microscopy, fluorimetry, and solution NMR spectroscopy and chemical shift index analysis) to characterize the interactions of these peptides with actin and Ca(2+)-CaM. Our results show that all three peptides can adopt α-helical structure inherently even in aqueous solution. Both α1- and α3-peptides showed strong binding with Ca(2+)-CaM, and both adopted an α-helical conformation upon interaction, but the binding of the α3-peptide appeared to be more dynamic. Only the α1-peptide exhibited actin polymerization and bundling activity, and the addition of Ca(2+)-CaM resulted in depolymerization of actin that had been polymerized by α1. The results of this study proved that there is an N-terminal binding domain in MBP for Ca(2+)-CaM (in addition to the primary site located in the C-terminus), and that it is sufficient for CaM-induced actin depolymerization. These three domains of MBP represent molecular recognition fragments with multiple roles in both membrane- and protein-association.

The most refined nucleus-independent chemical shift index (NICS(0)(πzz)) and the extra cyclic resonance energies (ECREs), based on the block localized wave function (BLW) method, show that the aromaticity of all azines is like that of benzene. The same is true for aza-naphthalenes relative to naphthalene. The lower relative energies of isomers with vicinal N's are due to the weakness of NN bonds rather than to reduced aromaticity.

The LAH4 family of histidine-rich peptides exhibits potent antimicrobial and DNA transfection activities, both of which require interactions with cellular membranes. The bilayer association of the peptides has been shown to be strongly pH-dependent, with in-planar alignments under acidic conditions and transmembrane orientations when the histidines are discharged. Therefore, we investigated the pH- and temperature-dependent conformations of LAH4 in DPC micellar solutions and in a TFE/PBS solvent mixture. In the presence of detergent and at pH 4.1, LAH4 adopts helical conformations between residues 9 and 24 concomitantly with a high hydrophobic moment. At pH 6.1, a helix-loop-helix structure forms with a hinge encompassing residues His10–Ala13. The data suggest that the high density of histidine residues and the resulting electrostatic repulsion lead to both a decrease in the pK values of the histidines and a less stable α-helical conformation of this region. The hinged structure at pH 6.1 facilitates membrane anchoring and insertion. At pH 7.8, the histidines are uncharged and an extended helical conformation including residues 4–21 is again obtained. LAH4 thus exhibits a high degree of conformational plasticity. The structures provide a stroboscopic view of the conformational changes that occur during membrane insertion, and are discussed in the context of antimicrobial activity and DNA transfection.

The study of the interaction of ghrelin (1), the endogenous ligand for the GH secretagogues receptor (GHS-R1a), and des-acyl ghrelin (2) with the GHS-R1a by NMR using living cells is presented, using GHS-R1a stably transfected cell lines (CHO and HEK 293) and wild type cells. Therefore, the interaction of 1 and 2 with the GHS-R1a receptor has been performed using quasi-physiological conditions. Ghrelin (1), showed a higher number of residues affected by chemical shift perturbation (CSP) or chemical shift exchange (CSE) effects: Ser3, Phe4, Leu5, Val12, Gln13/Gln14, Lys16/Lys19, Glu17 and Lys24 were much more affected in 1 than in des-acyl ghrelin (2). The chemical shift index CSI values indicated the presence of a possible α-helical region between Glu8 and Lys20 for ghrelin (1). After analysing the NMR data, two possible structures have arisen, which present different proline rotamers: the EEZE and the EZEZ conformers, at positions Pro7, Pro21, Pro22 and Pro27, respectively, keeping a left-handed α-helix from Glu8 to Lys20. These experimental evidences might imply that the GHS-R1a receptor is acting as a prolyl-cis/trans isomerase. © 2010 Elsevier Ltd. All rights reserved.

Estimation of secondary structure in polypeptides is important for studying their structure, folding and dynamics. In NMR spectroscopy, such information is generally obtained after sequence specific resonance assignments are completed. We present here a new methodology for assignment of secondary structure type to spin systems in proteins directly from NMR spectra, without prior knowledge of resonance assignments. The methodology, named Combination of Shifts for Secondary Structure Identification in Proteins (CSSI-PRO), involves detection of specific linear combination of backbone (1)H(alpha) and (13)C' chemical shifts in a two-dimensional (2D) NMR experiment based on G-matrix Fourier transform (GFT) NMR spectroscopy. Such linear combinations of shifts facilitate editing of residues belonging to alpha-helical/beta-strand regions into distinct spectral regions nearly independent of the amino acid type, thereby allowing the estimation of overall secondary structure content of the protein. Comparison of the predicted secondary structure content with those estimated based on their respective 3D structures and/or the method of Chemical Shift Index for 237 proteins gives a correlation of more than 90% and an overall rmsd of 7.0%, which is comparable to other biophysical techniques used for structural characterization of proteins. Taken together, this methodology has a wide range of applications in NMR spectroscopy such as rapid protein structure determination, monitoring conformational changes in protein-folding/ligand-binding studies and automated resonance assignment.

Angiogenin is a potent inducer of blood vessel formation and is overexpressed in many cancers. It is a member of the pancreatic RNase superfamily. Angiogenin has 33% sequence identity with RNase A. All RNases family proteins except angiogenin and turtle RNase have four disulfide bonds. Angiogenin has three disulfide bonds, [C27-C82], [C40-C93] and [C58-C108], respectively. To examine the role of disulfide bond [C40-C93] near the nuclear localization site, recombinant des [40-93] mutant by replacement of cysteines at positions 40 and 93 with serines was cloned, expressed, and purified. CD spectrum of des [40-93] was similar to that of wild type and implies that the overall folding of two forms is similar. 2D and 3D NMR experiments for 15 N and/or 13 C-isotope labeled des [40-93] were performed and the backbone resonance assignment was done. Based on the result of backbone assignment and chemical shift index (CSI), we confirmed that deletion of [C40-C93] disulfide bond affected the local structural stability of bovine angiogenin near the nuclear localization site and ribonucleolytic active site. Based on this assignment, studies on structure and dynamics of des [40-93] mutant will be performed.

F 1 F o -ATP synthase is a large multiprotein complex, including at least 10 subunits in the membrane-bound F o -sector. One of these F o proteins is subunit e (Su e), involved in the stable dimerization of F 1 F o -ATP synthase, and required for the establishment of normal cristae membrane architecture. As a step toward enabling structure-function studies of the F o -sector, the Su e transmembrane region was structurally characterized in micelles. Based on a series of NMR and CD (circular dichroism) studies, a structural model of the Su e/micelle complex was constructed, indicating Su e is largely helical, and emerges from the micelle with Arg20 near the phosphate head groups. Su e only adopts this folded conformation in the context of the micelle, and is essentially disordered in DMSO, water or trifluoroethanol/ water. Within the micelle the C-terminal Ala10-Arg20 stretch is helical, while the region N-terminal may be transiently helical, based on negative CSI (chemical shift index) values. The Ala10-Arg20 helix contains the G 14 XXXG 18 motif, which has been proposed to play an important role in dimer formation with another protein from the F o -sector. The Gly on the C-terminal end of this motif (Glyl8) is slightly more mobile than the more buried Gly 14, based on NMR order parameter measurements (Gly 14 S 2 = 0.950; Glyl8 S 2 = 0.895). Only one Su e transmembrane peptide is bound per micelle, and micelles are 22-23 A in diameter, composed of 51 ±4 dodecylphosphocholine detergent molecules. Although there is no evidence for Su e homodimerization via the transmembrane domain, potentially synergistic roles for N-terminal (membrane) and C-terminal (soluble) domain interactions may still occur. Furthermore, the presence of a buried charged residue (Arg7) suggests there may be interactions with other F o -sector protein(s) that stabilize this charge, and possibly drive the folding of the N-terminal 9 residues of the transmembrane domain.

The determination of crystal structures of the troponin complex (Takeda et al. 2003. Nature. 424:35–41; Vinogradova et al. 2005. Proc. Natl. Acad. Sci. USA. 102:5038–5043) has advanced knowledge of the regulation of muscle contraction at the molecular level. However, there are domains important for actin binding that are not visualized. We present evidence that the C-terminal region of troponin I (TnI residues 135–182) is flexible in solution and has no stable secondary structure. We use NMR spectroscopy to observe the backbone dynamics of skeletal [2H, 13C, 15N]-TnI in the troponin complex in the presence of Ca2+ or EGTA/Mg2+. Residues in this region give stronger signals than the remainder of TnI, and chemical shift index values indicate little secondary structure, suggesting a very flexible region. This is confirmed by NMR relaxation measurements. Unlike TnC and other regions of TnI in the complex, the C-terminal region of TnI is not affected by Ca2+ binding. Relaxation measurements and reduced spectral density analysis are consistent with the C-terminal region of TnI being a tethered domain connected to the rest of the troponin complex by a flexible linker, residues 137–146, followed by a collapsed region with at most nascent secondary structure.

CH(alpha) residual dipolar couplings (Deltardc's) were measured for the oxidized cytochrome b562 from Escherichia coli as a result of its partial self-orientation in high magnetic fields due to the anisotropy of the overall magnetic susceptibility tensor. Both the low spin iron (III) heme and the four-helix bundle fold contribute to the magnetic anisotropy tensor. CH(alpha) Deltardc's, which span a larger range than the analogous NH values (already available in the literature) sample large space variations at variance with NH Deltardc's, which are largely isooriented within alpha helices. The whole structure is now significantly refined with the chemical shift index and CH(alpha) Deltardc's. The latter are particularly useful also in defining the molecular magnetic anisotropy parameters. It is shown here that the backbone folding can be conveniently and accurately determined using backbone restraints only, which include NOEs, hydrogen bonds, residual dipolar couplings, pseudocontact shifts, and chemical shift index. All these restraints are easily and quickly determined from the backbone assignment. The calculated backbone structure is comparable to that obtained by using also side chain restraint. Furthermore, the structure obtained with backbone only restraints is, in its whole, very similar to that obtained with the complete set of restraints. The paramagnetism based restraints are shown to be absolutely relevant, especially for Deltardc's.

The wild-type (wt) N-terminal 23-residue fusion peptide (FP) of the human immunodeficiency virus (HIV) fusion protein gp41 and its V2E mutant have been studied by nuclear magnetic resonance (NMR) spectroscopy in dodecylphosphocholine (DPC) micelles as membrane mimics. A number of NMR techniques have been used. Pulsed field-gradient diffusion measurements in DPC and in 4:1 DPC/sodium dodecylsulfate mixed micelles showed that there is no major difference between the partition coefficients of the fusogenic wt peptide and the V2E mutant in these micelles, indicating that there is no correlation between the activity of the fusion peptides and their membrane affinities. The nuclear Overhauser enhancement (NOE) patterns and the chemical shift index for these two peptides indicated that both FP are in an α helical conformation between the Ile4 to Leu12 or to Ala15 region. Simulated annealing showed that the helical region extends from Ile4 to Met19. The two FPs share similar conformational characteristics, indicating that the conformation of the FP is not an important factor determining its activity. The spin-label studies, utilizing spin labels 5- and 16-doxystearic acids in the DPC micelles, provided clear indication that the wt FP inserts its N-terminus into the micelles while the V2E mutant does not insert into the micelles. The conclusion from the spin-label results is corroborated by deuterium amide proton exchange experiments. The correlation between the oblique insertion of the FP and its fusogenic activity is in excellent agreement with results from our molecular dynamics simulation and from other previous studies.

Random coil chemical shifts are commonly used to detect protein secondary structural elements in chemical shift index (CSI) calculations. Though this technique is widely used and seems reliable for folded proteins, the choice of reference random coil chemical shift values can significantly alter the outcome of secondary structure estimation. In order to evaluate these effects, we present a comparison of secondary structure content calculated using CSI, based on five different reference random coil chemical shift value sets, to that derived from three-dimensional structures. Our results show that none of the reference random coil data sets chosen for evaluation fully reproduces the actual secondary structures. Among the reference values generally available to date, most tend to be good estimators only of helices. Based on our evaluation, we recommend the experimental values measured by Schwarzinger et al.(2000), and statistical values obtained by Lukin et al. (1997), as good estimators of both helical and sheet content.

Fibroblast growth factors (FGFs) regulate a wide range of important cellular processes. The biological activities of FGFs are mediated by cell surface receptors (FGFRs). In the present study for the first time we report the cloning, expression, and characterization of the ligand (FGF)-binding D2 domain of human FGFR2. D2 domain is expressed in Escherichia coli in high yields (10 mg/L) as inclusion bodies. The protein is recovered by dissolving the inclusion bodies in 8 M urea and subsequently refolding on nickel affinity column. The protein is purified (to approximately 97% purity) to homogeneity using heparin-Sepharose affinity column. Far-UV circular dichroism data and chemical shift index plot based on 1H-alpha, 13C-alpha, 13C-beta, and 13carbonyl group chemical shifts suggest that D2 domain is an all beta-sheet protein consisting of 9 beta-strands. Isothermal titration calorimetry and equilibrium urea unfolding experiments show that recombinant D2 domain is in a biologically active conformation and binds strongly to its ligand (FGF) and to the heparin analog, sucrose octasulfate (SOS). Using a variety of triple resonance NMR experiments, complete assignment of 1H, 15N, and 13C resonances in D2 domain has been accomplished. The findings of the present study not only pave way for an in-depth investigation of the molecular mechanism(s) underlying the activation of FGF signaling but also provide avenues for the rational design of potent inhibitors against FGF-mediated pathogenesis.

All cytokines belonging to the interleukin-6 (IL-6)-type family of cytokines utilize receptors that have a modular build of several immunoglobulin-like and fibronectin type III-like domains. Characteristic of these receptors is a cytokine receptor homology region consisting of two such fibronectin domains defined by a set of four conserved cysteines and a tryptophan-serine-X-tryptophan-serine sequence motif. On target cells, interleukin-6 first binds to its specific receptor and subsequently to a homodimer of the signal transducer protein gp130. The interleukin-6 receptor consists of three extracellular domains. The N-terminal immunoglobulin-like domain is not involved in ligand binding, whereas the third membrane proximal fibronectin-like domain accounts for more than 90% of the binding energy to IL-6. Here, the key residues of this fibronectin-like domain involved in the interaction with IL-6 are described. Chemical shift mapping data with 15N-labeled IL-6R-D3 and unlabeled IL-6 coupled with recent structural data clearly reveal the epitope within the IL-6R-D3 responsible for mediating the high affinity interaction with its cognate cytokine.

A large number of bioactive peptides have been isolated from amphibian skin secretions. These peptides have a variety of actions including antibiotic and anticancer activities and the inhibition of neuronal nitric oxide synthase. We have investigated the structure-activity relationship of citropin 1.1, a broad-spectrum antibiotic and anticancer agent that also causes inhibition of neuronal nitric oxide synthase, by making a number of synthetically modified analogues. Citropin 1.1 has been shown previously to form an amphipathic alpha-helix in aqueous trifluoroethanol. The results of the structure-activity studies indicate the terminal residues are important for bacterial activity and increasing the overall positive charge, while maintaining an amphipathic distribution of residues, increases activity against Gram-negative organisms. Anticancer activity generally mirrors antibiotic activity suggesting a common mechanism of action. The N-terminal residues are important for inhibition of neuronal nitric oxide synthase, as is an overall positive charge greater than three. The structure of one of the more active synthetic modifications (A4K14-citropin 1.1) was determined in aqueous trifluoroethanol, showing that this peptide also forms an amphipathic alpha-helix.

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.

Bleomycin (Bm), a 1.4 kDa glycopeptide excreted by Streptomyces verticillus, is a natural antibacterial compound used in therapy as antineoplastic drug. To counteract its biological activity, cells have developed several resistance mechanisms, one of these based on proteins able to tightly bind Bm. In this paper, the interaction of Zn(2+)-Bm with the Streptoalloteichus hindustanus Bm resistance protein (ShBle) has been investigated by solution state NMR. Sequential nOe and chemical shift index have shown that the fold of the protein (in absence or presence of Bm) is identical to the previously published X-ray structure. The dimeric nature of ShBle is confirmed by the diffusion tensor as determined by NMR relaxation data. Using isotope filtered nOe experiment, intermolecular nOes between Bm and ShBle have been observed as used for modeling. While the interaction of the Bm metal binding site with ShBle appears to be uniquely defined, several conformations of the bithiazole moieties are compatible with the NMR data. Binding of Bm also induces changes of the local dynamics (stretch N85-G91), as shown by (15)N relaxation data. These results are discussed in the context of several Bm analogues able to interact with ShBle and of the recently published X-rays structures.

The importance of protein chemical shift values for the determination of three-dimensional protein structure has increased in recent years because of the large databases of protein structures with assigned chemical shift data. These databases have allowed the investigation of the quantitative relationship between chemical shift values obtained by liquid state NMR spectroscopy and the three-dimensional structure of proteins. A neural network was trained to predict the 1H, 13C, and 15N of proteins using their three-dimensional structure as well as experimental conditions as input parameters. It achieves root mean square deviations of 0.3 ppm for hydrogen, 1.3 ppm for carbon, and 2.6 ppm for nitrogen chemical shifts. The model reflects important influences of the covalent structure as well as of the conformation not only for backbone atoms (as, e.g., the chemical shift index) but also for side-chain nuclei. For protein models with a RMSD smaller than 5 A a correlation of the RMSD and the r.m.s. deviation between the predicted and the experimental chemical shift is obtained. Thus the method has the potential to not only support the assignment process of proteins but also help with the validation and the refinement of three-dimensional structural proposals. It is freely available for academic users at the PROSHIFT server: www.jens-meiler.de/proshift.html

The alpha subunit of tryptophan synthase (alphaTS) from S. typhimurium belongs to the triosephosphate isomerase (TIM) or the (beta/alpha)(8) barrel fold, one of the most common structures in biology. To test the conservation of the global fold in the isolated Escherichia coli homolog, we have obtained a majority of the backbone assignments for the 29-kD alphaTS by using standard heteronuclear multidimensional NMR methods on uniformly (15)N- and (15)N/(13)C-labeled protein and on protein selectively (15)N-labeled at key hydrophobic residues. The secondary structure mapped by chemical shift index, nuclear Overhauser enhancements (NOEs), and hydrogen-deuterium (H-D) exchange, and several abnormal chemical shifts are consistent with the conservation of the global TIM barrel fold of the isolated E. coli alphaTS. Because most of the amide protons that are slow to exchange with solvent correspond to the beta-sheet residues, the beta-barrel is likely to play an important role in stabilizing the previously detected folding intermediates for E. coli alphaTS. A similar combination of uniform and selective labeling can be extended to other TIM barrel proteins to obtain insight into the role of the motif in stabilizing what appear to be common partially folded forms.

The present work describes the dynamics of the apo form of cytochrome b(562), a small soluble protein consisting of 106 amino acid residues [Itagaki, E., and Hager, L. P. (1966) J. Biol. Chem. 241, 3687-3695]. The presence of exchange in the millisecond time scale is demonstrated for the last part of helix IV (residues 95-105 in the holo form). The chemical shift index analysis [Wishart, D. S., and Sykes, B. D. (1994) J. Biomol. NMR 4, 171-180] based on H(alpha), C(alpha), C(beta), and C' chemical shifts suggests a larger helical content than shown in the NMR structure based on NOEs. These results indicate the presence of helical-like conformations participating in the exchange process. This hypothesis is consistent with amide deuterium exchange rates and the presence of some hydrogen bonds identified from amide chemical shift temperature coefficients [Baxter, N. J., and Williamson, M. P. (1997) J. Biomol. NMR 9, 359-369]. (15)N relaxation indicates limited mobility for the amide protons of this part of the helix in the picosecond time scale. A 30 ns stochastic dynamics simulation shows small fluctuations around the helical conformation on this time scale. These fluctuations, however, do not result in a significant decrease of the calculated order parameters which are consistent with the experimental (15)N relaxation data. These results resolve an apparent discrepancy in the NMR structures between the disorder observed in helix IV due to a lack of NOEs and the secondary structure predictions based on H(alpha) chemical shifts [Feng, Y., Wand, A. J., and Sligar, S. G. (1994) Struct. Biol. 1, 30-35].

A computational approach has been developed to assess the power of paramagnetism-based backbone constraints with respect to the determination of the tertiary structure, once the secondary structure elements are known. This is part of the general assessment of paramagnetism-based constraints which are known to be relevant when used in conjunction with all classical constraints. The paramagnetism-based constraints here investigated are the pseudocontact shifts, the residual dipolar couplings due to self-orientation of the metalloprotein in high magnetic fields, and the cross correlation between dipolar relaxation and Curie relaxation. The relative constraints are generated by back-calculation from a known structure. The elements of secondary structure are supposed to be obtained from chemical shift index. The problem of the reciprocal orientation of the helices is addressed. It is shown that the correct fold can be obtained depending on the length of the α-helical stretches with respect to the length of the non helical segments connecting the α-helices. For example, the correct fold is straightforwardly obtained for the four-helix bundle protein cytochrome b 562, while the double EF-hand motif of calbindin D9k is hardly obtained without ambiguity. In cases like calbindin D9k, the availability of datasets from different metal ions is helpful, whereas less important is the location of the metal ion with respect to the secondary structure elements.

The intrinsically disordered 18.5 kDa classic isoform of MBP (myelin basic protein) interacts with Fyn kinase during oligodendrocyte development and myelination. It does so primarily via a central proline-rich SH3 (Src homology 3) ligand (T92–R104, murine 18.5 kDa MBP sequence numbering) that is part of a molecular switch due to its high degree of conservation and modification by MAP (mitogen-activated protein) and other kinases, especially at residues T92 and T95. Here, we show using co-transfection experiments of an early developmental oligodendroglial cell line (N19) that an MBP segment upstream of the primary ligand is involved in MBP–Fyn–SH3 association in cellula. Using solution NMR spectroscopy in vitro, we define this segment to comprise MBP residues (T62–L68), and demonstrate further that residues (V83–P93) are the predominant SH3-target, assessed by the degree of chemical shift change upon titration. We show by chemical shift index analysis that there is no formation of local poly-proline type II structure in the proline-rich segment upon binding, and by NOE (nuclear Overhauser effect) and relaxation measurements that MBP remains dynamic even while complexed with Fyn–SH3. The association is a new example first of a non-canonical SH3-domain interaction and second of a fuzzy MBP complex., MBP interacts with Fyn kinase during oligodendrocyte development and myelination. We show that there is no binding-induced PPII formation in the primary ligand segment, and that a region upstream is required for in vitro interaction.

Protein chemical shifts have long been used by NMR spectroscopists to assist with secondary structure assignment and to provide useful distance and torsion angle constraint data for structure determination. One of the most widely used methods for secondary structure identification is called the Chemical Shift Index (CSI). The CSI method uses a simple digital chemical shift filter to locate secondary structures along the protein chain using backbone 13C and 1H chemical shifts. While the CSI method is simple to use and easy to implement, it is only about 75–80 % accurate. Here we describe a significantly improved version of the CSI (2.0) that uses machine-learning techniques to combine all six backbone chemical shifts (13Cα, 13Cβ, 13C, 15N, 1HN, 1Hα) with sequence-derived features to perform far more accurate secondary structure identification. Our tests indicate that CSI 2.0 achieved an average identification accuracy (Q3) of 90.56 % for a training set of 181 proteins in a repeated tenfold cross-validation and 89.35 % for a test set of 59 proteins. This represents a significant improvement over other state-of-the-art chemical shift-based methods. In particular, the level of performance of CSI 2.0 is equal to that of standard methods, such as DSSP and STRIDE, used to identify secondary structures via 3D coordinate data. This suggests that CSI 2.0 could be used both in providing accurate NMR constraint data in the early stages of protein structure determination as well as in defining secondary structure locations in the final protein model(s). A CSI 2.0 web server ( http://csi.wishartlab.com ) is available for submitting the input queries for secondary structure identification.