Solution structure of the calcium channel antagonist omega-conotoxin GVIA.

The three-dimensional solution structure is reported for omega-conotoxin GVIA, which is a potent inhibitor of presynaptic calcium channels in vertebrate neuromuscular junctions. Structures were generated by a hybrid distance geometry and restrained molecular dynamics approach using interproton distance, torsion angle, and hydrogen-bonding constraints derived from 1H NMR data. Conformations of GVIA with low constraint violations converged to a common peptide fold. The secondary structure in the peptide is an antiparallel triple-stranded beta-sheet containing a beta-hairpin and three tight turns. The NMR data are consistent with the region of the peptide from residues S9 to C16 being more dynamic than the rest of the peptide. The peptide has an amphiphilic structure with a positively charged hydrophilic side and an opposite side that contains a small hydrophobic region. Residues that are thought to be important in binding and function are located on the hydrophilic face of the peptide.     

Synthesis and biological characterization of a series of analogues of ω-conotoxin gvia

The 27-residue polypeptide ω-conotoxin GVIA (ω-CgTx), from the venom of the cone shell Conus geographus, blocks N-type neuronal calcium channels. It contains three disulphide bridges. We reporte here the synthesis and biological characterization of a seires of analogues in which one disulphide has been replaced by substitution of appropriate Cys residues with Ser, viz. [Ser^1,16]-ω-CgTx, [Ser^8,19]-ω-CgTx, [Ser^15,26]-ω-CgTx, [Ser¹⁶]-ω-CgTx_8-27 and [Ser¹⁵]-ω-CgTx_1-19. All syntheses were conducted manually using either Boc or Fmoc methodology. Deprotected peptides were oxidized to their bridged forms using either aerial oxidation or aqueous dimethyl sulphoxide. Peptides were purified using RP-HPLC, and their purity and identity were checked by RP-HPLC, capillary electrophoresis and mass spectrometry. Inhibition of neuronal N-type calcium channels was assessed as the inhibition of the twitch responses of rat vas deferens stimualted with single electrical pulses at 20 second intervals. None of these analogues was biologically active, suggesting that the disulphides play an important role in maintaining biological activity.     

“Ensemble” iterative relaxation matrix approach: A new NMR refinement protocol applied to the solution structure of crambin

The structure in solution of crambin, a small protein of 46 residues, has been determined from 2D NMR data using an iterative relaxation matrix approach (IRMA) together with distance geometry, distance bound driven dynamics, molecular dynamics, and energy minimization. A new protocol based on an “ensemble” approach is proposed and compared to the more standard initial rate analysis approach and a “single structure” relaxation matrix approach. The effects of fast local motions are included and R-factor calculations are performed on NOE build-ups to describe the quality of agreement between theory and experiment. A new method for stereospecific assignment of prochiral groups, based on a comparison of theoretical and experimental NOE intensities, has been applied. The solution structure of crambin could be determined with a precision (rmsd from the average structure) of 0.7 Å on backbone atoms and 1.1 Å on all heavy atoms and is largely similar to the crystal structure with a small difference observed in the position of the side chain of Tyr-29 which is determined in solution by both J-coupling and NOE data. Regions of higher structural variability (suggesting higher mobility) are found hi the solution structure, in particular for the loop between the two helices (Gly-20 to Pro-22). © 1993 Wiley-Liss, Inc.     

Solid-phase Synthesis of ω-Agatoxin IVA,a P-type Calcium Channel Blocker

ω-Agatoxin IVA, isolated from the venom of funnel web spider Agelenopsis aperta, blocks potently and selectively P-type calcium channels. This toxin, composed of 48 amino acids and containing 8 cysteine residues, was synthesized by the solid-phase procedure. The Cys residues were protected by acetamidomethyl (Acm) groups which were removed by mercuric acetate. During treatment with mercuric acetate, a by-product was detected, involving modification of tryptophan residues by the Acm groups. This side reaction can be completely prevented by addition of an excess of tryptophan in the reaction medium during Acm deprotection. The resulting peptide was submitted to an oxidative refolding, in different conditions, in order to determine the most favourable protocol. After formation of the four disulphide bonds, the toxin was purified by successive preparative HPLC, on two different supports, and fully characterized by analytical HPLC, capillary electrophoresis, amino acid analysis, mass spectrometry and Edman degradation. It was found to block the P-type calcium channel with a similar biological potency as described for the natural product.     

Solution secondary structure of calcium-saturated troponin C monomer determined by multidimensional heteronuclear NMR spectroscopy.

The solution secondary structure of calcium-saturated skeletal troponin C (TnC) in the presence of 15% (v/v) trifluoroethanol (TFE), which has been shown to exist predominantly as a monomer (Slupsky CM, Kay CM, Reinach FC, Smillie LB, Sykes BD, 1995, Biochemistry 34, forthcoming), has been investigated using multidimensional heteronuclear nuclear magnetic resonance spectroscopy. The 1H, 15N, and 13C NMR chemical shift values for TnC in the presence of TFE are very similar to values obtained for calcium-saturated NTnC (residues 1-90 of skeletal TnC), calmodulin, and synthetic peptide homodimers. Moreover, the secondary structure elements of TnC are virtually identical to those obtained for calcium-saturated NTnC, calmodulin, and the synthetic peptide homodimers, suggesting that 15% (v/v) TFE minimally perturbs the secondary and tertiary structure of this stably folded protein. Comparison of the solution structure of calcium-saturated TnC with the X-ray crystal structure of half-saturated TnC reveals differences in the phi/psi angles of residue Glu 41 and in the linker between the two domains. Glu 41 has irregular phi/psi angles in the crystal structure, producing a kink in the B helix, whereas in calcium-saturated TnC, Glu 41 has helical phi/psi angles, resulting in a straight B helix. The linker between the N and C domains of calcium-saturated TnC is flexible in the solution structure.     

Solution structure of BmKK2, a new potassium channel blocker from the venom of chinese scorpion Buthus martensi Karsch

A natural K+ channel blocker, BmKK2 (a member of scorpion toxin subfamily alpha-KTx 14), which is composed of 31 amino acid residues and purified from the venom of the Chinese scorpion Buthus martensi Karsch, was characterized using whole-cell patch-clamp recording in rat hippocampal neurons. The three dimensional structure of BmKK2 was determined with two-dimensional NMR spectroscopy and molecular modelling techniques. In solution this toxin adopted a common alpha/beta-motif, but showed distinct local conformation in the loop between alpha-helix and beta-sheet in comparison with typical short-chain scorpion toxins (e.g., CTX and NTX). Also, the alpha helix is shorter and the beta-sheet element is smaller (each strand consisted only two residues). The unusual structural feature of BmKK2 was attributed to the shorter loop between the alpha-helix and beta-sheet and the presence of two consecutive Pro residues at position 21 and 22 in the loop. Moreover, two models of BmKK2/hKv1.3 channel and BmKK2/rSK2 channel complexes were simulated with docking calculations. The results demonstrated the existence of a alpha-mode binding between the toxin and the channels. The model of BmKK2/rSK2 channel complex exhibited favorable contacts both in electrostatic and hydrophobic, including a network of five hydrogen bonds and bigger interface containing seven pairs of inter-residue interactions. In contrast, the model of BmKK2/hKv1.3 channel complex, containing only three pairs of inter-residue interactions, exhibited poor contacts and smaller interface. The results well explained its lower activity towards Kv channel, and predicted that it may prefer a type of SK channel with a narrower entryway as its specific receptor.     

Oxidative folding of ω-conotoxin MVIIC: Effects of temperature and salt

Oxidative folding of o-conotoxin MVIIC, a highly basic 26-amino acid peptide with three disulfide bonds, predominantly gave two products with mismatched disulfide bonds in 0.1M NH₄OAc buffer (pH 7.7) at 21°C both in the presence and absence of redox reagents such as reduced and oxidized glutathione. A low reaction temperature (5°C) and a high salt concentration in buffer such as 2M (NH₄)₂SO₄ were necessary to obtain the correctly folded biologically active product. The folding reaction was found to proceed via a two-stage pathway of (I) the formation and (II) the rearrangement of the mismatched disulfide bonds. Both the reaction temperature and the salt strongly affected the equilibrium between mismatched and correctly formed disulfide bonds in the second stage. Such an effect of salts on the rearrangement reaction could be explained by anion binding at a low concentration and the salting out effect at a high concentration by analyzing the rank order of their effectiveness. The anion-binding effect was also confirmed by examining the folding of the tetra-acetylated peptide at the Lys side chains. CD study suggested that the yield of the biologically active product was correlated with its conformational change as functions of temperature and salt concentration. © 1996 John Wiley & Sons, Inc.     

Solution NMR resonance assignment strategies for β‐barrel membrane proteins

Membrane proteins in detergent micelles are large and dynamic complexes that present challenges for solution NMR investigations such as spectral overlap and line broadening. In this study, multiple methods are introduced to facilitate resonance assignment of β‐barrel membrane proteins using Opa₆₀ from Neisseria gonorrhoeae as a model system. Opa₆₀ is an eight‐stranded β‐barrel with long extracellular loops (～63% of the protein) that engage host receptors and induce engulfment of the bacterium. The NMR spectra of Opa₆₀ in detergent micelles exhibits significant spectral overlap and resonances corresponding to the loop regions had variable line widths, which interfered with a complete assignment of the protein. To assign the β‐barrel residues, trypsin cleavage was used to remove much of the extracellular loops while preserving the detergent solubilized β‐barrel. The removal of the loop resonances significantly improved the assignment of the Opa₆₀ β‐barrel region (97% of the resonances corresponding to the β‐barrel and periplasmic turns were assigned). For the loop resonance assignments, two strategies were implemented; modulating temperature and synthetic peptides. Lowering the temperature broadened many peaks beyond detection and simplified the spectra to only the most dynamic regions of the loops facilitating 27 loop resonances to be assigned. To further assign functionally important and unstructured regions of the extracellular loops, a synthetic 20 amino acid peptide was synthesized and had nearly complete spectral overlap with the full‐length protein allowing 17 loop resonances to be assigned. Collectively, these strategies are effective tools that may accelerate solution NMR structure determination of β‐barrel membrane proteins.     

Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

Solution structure of the phosphocarrier protein HPr from Bacillus subtilis by two-dimensional NMR spectroscopy.

The solution structure of the phosphocarrier protein, HPr, from Bacillus subtilis has been determined by analysis of two-dimensional (2D) NMR spectra acquired for the unphosphorylated form of the protein. Inverse-detected 2D (1H-15N) heteronuclear multiple quantum correlation nuclear Overhauser effect (HMQC NOESY) and homonuclear Hartmann-Hahn (HOHAHA) spectra utilizing 15N assignments (reported here) as well as previously published 1H assignments were used to identify cross-peaks that are not resolved in 2D homonuclear 1H spectra. Distance constraints derived from NOESY cross-peaks, hydrogen-bonding patterns derived from 1H-2H exchange experiments, and dihedral angle constraints derived from analysis of coupling constants were used for structure calculations using the variable target function algorithm, DIANA. The calculated models were refined by dynamical simulated annealing using the program X-PLOR. The resulting family of structures has a mean backbone rmsd of 0.63 A (N, C alpha, C', O atoms), excluding the segments containing residues 45-59 and 84-88. The structure is comprised of a four-stranded antiparallel beta-sheet with two antiparallel alpha-helices on one side of the sheet. The active-site His 15 residue serves as the N-cap of alpha-helix A, with its N delta 1 atom pointed toward the solvent to accept the phosphoryl group during the phosphotransfer reaction with enzyme I. The existence of a hydrogen bond between the side-chain oxygen atom of Tyr 37 and the amide proton of Ala 56 is suggested, which may account for the observed stabilization of the region that includes the beta-turn comprised of residues 37-40. If the beta alpha beta beta alpha beta (alpha) folding topology of HPr is considered with the peptide chain polarity reversed, the protein fold is identical to that described for another group of beta alpha beta beta alpha beta proteins that include acylphosphatase and the RNA-binding domains of the U1 snRNP A and hnRNP C proteins.     

Solution NMR structure determination of proteins revisited

This 'Perspective' bears on the present state of protein structure determination by NMR in solution. The focus is on a comparison of the infrastructure available for NMR structure determination when compared to protein crystal structure determination by X-ray diffraction. The main conclusion emerges that the unique potential of NMR to generate high resolution data also on dynamics, interactions and conformational equilibria has contributed to a lack of standard procedures for structure determination which would be readily amenable to improved efficiency by automation. To spark renewed discussion on the topic of NMR structure determination of proteins, procedural steps with high potential for improvement are identified.     

Automated screening for metabolites in complex mixtures using 2D COSY NMR spectroscopy

One of the greatest challenges in metabolomics is the rapid and unambiguous identification and quantification of metabolites in a biological sample. Although one-dimensional (1D) proton nuclear magnetic resonance (NMR) spectra can be acquired rapidly, they are complicated by severe peak overlap that can significantly hinder the automated identification and quantification of metabolites. Furthermore, it is currently not reasonable to assume that NMR spectra of pure metabolites are available a priori for every metabolite in a biological sample. In this paper we develop and report on tests of methods that assist in the automatic identification of metabolites using proton two-dimensional (2D) correlation spectroscopy (COSY) NMR. Given a database of 2D COSY spectra for the metabolites of interest, our methods provide a list sorted by a heuristic likelihood of the metabolites present in a sample that has been analyzed using 2D COSY NMR. Our models attempt to correct the displacement of the peaks that can occur from one sample to the next, due to pH, temperature and matrix effects, using a statistical and chemical model. The correction of one peak can result in an implied correction of others due to spin–spin coupling. Furthermore, these displacements are not independent: they depend on the relative position of functional groups in the molecule. We report experimental results using defined mixtures of amino acids as well as real complex biological samples that demonstrate that our methods can be very effective at automatically and rapidly identifying metabolites.     

NMR structure of duplex DNA containing the α‐OH‐PdG·dA base pair: A mutagenic intermediate of acrolein

Acrolein, a cell metabolic product and main component of cigarette smoke, reacts with DNA generating α‐OH‐PdG lesions, which have the ability to pair with dATP during replication thereby causing G to T transversions. We describe the solution structure of an 11‐mer DNA duplex containing the mutagenic α‐OH‐PdG·dA base pair intermediate, as determined by solution nuclear magnetic resonance (NMR) spectroscopy and retrained molecular dynamics (MD) simulations. The NMR data support a mostly regular right‐handed helix that is only perturbed at its center by the presence of the lesion. Undamaged residues of the duplex are in anti orientation, forming standard Watson‐Crick base pairs alignments. Duplication of proton signals at and near the damaged base pair reveals the presence of two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. The α‐OH‐PdG adduct assumes a syn conformation pairing to its partner dA base that is protonated at pH 6.6. The three‐dimensional structure obtained by restrained molecular dynamics simulations show hydrogen bond interactions that stabilize α‐OH‐PdG in a syn conformation and across the lesion containing base pair. We discuss the implications of the structures for the mutagenic bypass of acrolein lesions. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 391–401, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

1H NMR structure of an antifungal γ-thionin protein SIα1: Similarity to scorpion toxins

The three-dimensional structure of the Sorghum bicolor seed protein γ-thionin SIα1 has been determined by 2D ¹H nuclear magnetic resonance (NMR) spectroscopy. The secondary structure of this 47-residue antifungal protein with four disulphide bridges consists of a three-stranded antiparallel sheet and one helix. The helix is tethered to the sheet by two disulphide bridges which link two successive turns of the helix to alternate residues i, i + 2 in one strand. Possible binding sites for antifungal activity are discussed. The same fold has been observed previously in several scorpion toxins. Proteins 32:334–349, 1998. © 1998 Wiley-Liss, Inc.     

The solution structure of BmTx3B, a member of the scorpion toxin subfamily alpha-KTx 16

This article reports the solution structure of BmTx3B (alpha-KTx16.2), a potassium channel blocker belonging to the subfamily alpha-KTx16, purified from the venom of the Chinese scorpion Buthus martensi Karsch. In solution, BmTx3B assumes a typical CSalphabeta motif, with an alpha-helix connected to a triple-stranded beta-sheet by 3 disulfide bridges, which belongs to the first structural group of short-chain scorpion toxins. On the other hand, BmTx3B is quite different from other toxins (such as ChTx and AgTx2) of this group in terms of the electrostatic and hydrophobic surface distribution. The functional surface (beta-face) of the molecule is characterized by less basic residues (only 2: Lys28 and Arg35) and extra aromatic residues (Phe1, Phe9, Trp15, and Tyr37). The peptide shows a great preference for the Kca1.1 channel over the Kv channel (about a 10(3)-fold difference). The model of BmTx3B/Kca1.1 channel complex generated by docking and dynamic simulation reveals that the stable binding between the BmTx3B and Kca1.1 channel is favored by a number of aromatic pi-pi stacking interactions. The influences of these structural features on the kinetic behavior of the toxin binding to Kca1.1 channel are also discussed.     

Direct NOE refinement of biomolecular structures using 2D NMR data

Summary A set of computer programs called DINOSAUR has been developed, which allows the refinement of biomolecular structures directly from 2D NOE intensities. The NOE restraining potential implemented emphasises the weak intensities corresponding to larger distances which are more likely to determine the three-dimensional structure. An approximation based on a two-spin approximation is proposed for the gradient of the NOE intensities instead of the exact solution which is extremely time-consuming. The DINOSAUR routines have been implemented in various refinement programs (Distance bound Driven Dynamics, Molecular Dynamics and Energy Minimisation) and tested on an eight-residue model peptide.     

Assignments and structure determination of the catalytic domain of human fibroblast collagenase using 3D double and triple resonance NMR spectroscopy

We report here the backbone 1HN, 15N, 13C, 13CO, and 1H NMR assignmentsfor the catalytic domain of human fibroblast collagenase (HFC). Three independentassignment pathways (matching 1H, 13C, and 13CO resonances) were used to establishsequential connections. The connections using 13C resonances were obtained fromHNCOCA and HNCA experiments; 13CO connections were obtained from HNCO andHNCACO experiments. The sequential proton assignment pathway was established from a 3D(1H/15N) NOESY-HSQC experiment. Amino acid typing was accomplished using 13C and15N chemical shifts, specific labeling of 15N-Leu, and spin pattern recognition from DQF-COSY. The secondary structure was determined by analyzing the 3D (1H/15N) NOESY-HSQC. A preliminary NMR structure calculation of HFC was found to be in agreement withrecent X-ray structures of human fibroblast collagenase and human neutrophil collagenase aswell as similar to recent NMR structures of a highly homologous protein, stromelysin. Allthree helices were located; a five-stranded -sheet (four parallel strands, one antiparallelstrand) was also determined. -Sheet regions were identified by cross-stranddN and dNN connections and by strong intraresidue dN correlations, and were corroborated byobserving slow amide proton exchange. Chemical shift changes in a selectively 15N-labeledsample suggest that substantial structural changes occur in the active site cleft on the bindingof an inhibitor.     

Structural classification of αββ and ββα supersecondary structure units in proteins

We present a fully automatic structural classification of supersecondary structure units, consisting of two hydrogen-bonded β strands, preceded or followed by an α helix. The classification is performed on the spatial arrangement of the secondary structure elements, irrespective of the length and conformation of the intervening loops. The similarity of the arrangements is estimated by a structure alignment procedure that uses as similarity measure the root mean square deviation of superimposed backbone atoms. Applied to a set of 141 well-resolved nonhomologous protein structures, the classification yields 11 families of recurrent arrangements. In addition, fragments that are structurally intermediate between the families are found; they reveal the continuity of the classification. The analysis of the families shows that the α helix and β hairpin axes can adopt virtually all relative orientations, with, however, some preferable orientations; moreover, according to the orientation, preferences in the left/right handedness of the α–β connection are observed. These preferences can be explained by favorable side by side packing of the α helix and the β hairpin, local interactions in the region of the α–β connection or stabilizing environments in the parent protein. Furthermore, fold recognition procedures and structure prediction algorithms coupled to database-derived potentials suggest that the preferable nature of these arrangements does not imply their intrinsic stability. They usually accommodate a large number of sequences, of which only a subset is predicted to stabilize the motif. The motifs predicted as stable could correspond to nuclei formed at the very beginning of the folding process. Proteins 30:193–212, 1998. © 1998 Wiley-Liss, Inc.     

Determination of absolute configuration in 4‐aryl‐3,4‐dihydro‐2(1H)‐pyrimidones by high performance liquid chromatography and CD spectroscopy

The absolute configuration of three 4‐aryl‐3,4‐dihydro‐2(1H)‐pyrimidones (Biginelli compounds, DHPMs) was established by comparison of the typical circular dichroism (CD) spectra of individual enantiomers with reference samples of known absolute configuration. The enantiomers were obtained by semipreparative separation of racemic mixtures on a Chiralcel OD‐H chiral stationary phase. The method was used to establish the enantiopreference of various lipases in biocatalytic kinetic resolution experiments employing activated DHPM esters. Chirality 11:659–662, 1999. © 1999 Wiley‐Liss, Inc.