Open Access Research article

Structure of α-conotoxin BuIA: influences of disulfide connectivity on structural dynamics

Ai-Hua Jin1, Hemma Brandstaetter1, Simon T Nevin2, Chia Chia Tan1, Richard J Clark1, David J Adams2, Paul F Alewood1, David J Craik1 and Norelle L Daly1*

Author Affiliations

1 Institute for Molecular Bioscience, Australian Research Council Special Research Centre for Functional and Applied Genomics, The University of Queensland, Brisbane QLD 4072, Australia

2 School of Biomedical Sciences, The University of Queensland, Brisbane QLD 4072, Australia

For all author emails, please log on.

BMC Structural Biology 2007, 7:28  doi:10.1186/1472-6807-7-28

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6807/7/28


Received:3 November 2006
Accepted:20 April 2007
Published:20 April 2007

© 2007 Jin et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

α-Conotoxins have exciting therapeutic potential based on their high selectivity and affinity for nicotinic acetylcholine receptors. The spacing between the cysteine residues in α-conotoxins is variable, leading to the classification of sub-families. BuIA is the only α-conotoxin containing a 4/4 cysteine spacing and thus it is of significant interest to examine the structure of this conotoxin.

Results

In the current study we show the native globular disulfide connectivity of BuIA displays multiple conformations in solution whereas the non-native ribbon isomer has a single well-defined conformation. Despite having multiple conformations in solution the globular form of BuIA displays activity at the nicotinic acetylcholine receptor, contrasting with the lack of activity of the structurally well-defined ribbon isomer.

Conclusion

These findings are opposite to the general trends observed for α-conotoxins where the native isomers have well-defined structures and the ribbon isomers are generally disordered. This study thus highlights the influence of the disulfide connectivity of BuIA on the dynamics of the three-dimensional structure.

Background

Nicotinic acetylcholine receptors (nAChRs) play a key role in the central and peripheral nervous system and are involved in neuronal growth and plasticity, development, learning, memory and pain sensation [1-5]. α-Conotoxins, found in the venoms of marine molluscs belonging to the genus Conus, act as antagonists at nAChRs and consequently have a range of potential therapeutic applications [6-8].

Members of the α-conotoxin family inhibit either muscle, neuronal or both types of nAChRs and prevent channel opening [6,9,10]. They comprise 12 to 19 amino acids, including four highly conserved cysteine residues that form two disulfide bonds. Native α-conotoxins have a disulfide connectivity linking the first and the third, and the second and the fourth cysteine residues (CysI–CysIII and CysII–CysIV), which is referred to as the "globular" isomer. The alternative connectivities (CysI–CysIV and CysII–CysIII) and (CysI–CysII and CysIII–CysIV) are referred to as the ribbon and beads forms, respectively and are shown in Fig. 1. The beads form has so far not been reported for naturally occurring conotoxins, but the ribbon connectivity is characteristic of the χ-conotoxins [6,11,12]. The segments between cysteine residues define two loops in the peptide backbone, also shown in Fig. 1.

thumbnailFigure 1. Disulfide isomers of BuIA. (A) A schematic representation of the ribbon, globular and beads disulfide connectivities of BuIA, a conotoxin having the unique 4/4 loop spacing. (B) HPLC traces of ribbon and globular BuIA. Analytical HPLC was carried out using a C18 Phenomenex analytical column (250 × 4.6 mm) with 90% acetonitrile and 0.043% trifluroacetic acid as the eluting solvent B and 1% B per minute as the gradient. Ribbon BuIA elutes earlier than the globular form but has a slightly broadened peak under the conditions utilized.

The spacing between the cysteine residues in α-conotoxins is variable, leading to their classification into sub-families [6]. Table 1 shows the sequences of selected α-conotoxins with different cysteine spacings. The original α-conotoxins discovered in the 1980's, such as the paralytic toxin GI, have a 3/5 spacing, with three amino acids in the first and five in the second loop [13]. They compete mainly for nAChRs of the muscle endplate. Members of the "non-classical" α-conotoxins display 4/X spacings (where X is 3, 4, 6, or 7) and target mainly neuronal nAChRs, with the exception of EI, which is muscle specific [14,15]. In addition to the variations in cysteine spacing there is a significant degree of sequence variation and it is this diversity that results in the exquisite selectivity of the α-conotoxins. The subtype selectivity of selected α-conotoxins is illustrated in Table 1.

Table 1. Sequences and activities of selected α-conotoxins

BuIA is the only α-conotoxin with the unique 4/4 spacing (see Table 1 and Fig. 1) and is the first conotoxin identified in the piscivorous Conus bullatus species. Azam et al. [16] used conserved features of the α-conotoxin signal sequence and 3'-untranslated sequence region of α-conotoxin prepropeptide genes to clone the gene coding for BuIA. A synthetic version of the corresponding 13-amino acid mature peptide with the globular disulfide connectivity (Cys2–Cys8, Cys3–Cys13) is active against α3- and α6-containing nAChR subtypes, whereas α2- and α4-containing nAChRs are blocked less effectively. Strikingly, BuIA kinetically discriminates between β2- and β4-containing receptors, as the off-rates are rapid for β2-subunit, but very slow for β4-containing nAChRs. This feature appears to extend across different mammalian species, as similar results for mouse, rat and human nAChRs were obtained [16]. Thus the novel α-conotoxin BuIA appears to be a valuable probe to distinguish among nAChRs containing different α and β subunits.

Understanding the structure-activity relationships of α-conotoxins has been a major focus of conotoxin research in recent years and, as a result the structures of many of the known α-conotoxins have been determined. Conotoxins are particularly suitable for structural analysis with NMR spectroscopy because of their small size, solubility and generally well-defined conformations, and the majority of α-conotoxin structures have been determined using this technique [17,18]. Despite the range of sequences and spacings between the cysteine residues a well-defined consensus structural motif has been found for the α-conotoxins, with the major features being the restraints imposed by the conserved globular disulfide connectivity and a helical region centred around CysII [18]. In the current study we have examined the structures and activities of the globular isomer and the non-native ribbon isomer of BuIA with a Cys2–Cys13, Cys3–Cys8 connectivity and show that the globular isomer is structurally heterogeneous while the non-native ribbon isomer has a well-defined conformation.

Results

Orthogonal syntheses of globular and ribbon BuIA were achieved using solid phase peptide synthesis via Boc chemistry. Purification via HPLC yielded a single product for each isomer whose mass was determined using electrospray ionization mass spectrometry. Both isomers have an observed mass of 1311.6 Da (calculated mass 1311.5 Da). HPLC traces for the purified products are shown in Fig. 1B, from which it is clear that the ribbon form elutes earlier than the globular form, indicating that its hydrophobic residues are buried to a greater extent. Interestingly, under the HPLC conditions used the ribbon form has a slightly broader peak than the globular form despite having a more defined structure as described below. When prepared without selective protection of the cysteine residues BuIA folds predominately into the globular isomer in a 4:2:1 ratio (globular:ribbon:beads), consistent with this being the native form in venom, as is the case with all other α-conotoxins. For the oxidation profile of BuIA without selective protection see Additional file 1.

Additional File 1. Oxidative folding of BuIA. The data provided represent the oxidation of conotoxin BuIA without selective protection of the cysteine residues.

Format: DOC Size: 43KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

The globular and ribbon isomers purified from oxidation without selective protection were analyzed by NMR spectroscopy and found to be identical to their respective isomers isolated following selective formation of the disulfide bonds. However, the TOCSY and NOESY spectra revealed considerable differences between the two isomers, as illustrated in the upper panels of Fig. 2. The non-native ribbon form with the Cys2–Cys13, Cys3–Cys8 disulfide connectivity displayed spectra consistent with a well-defined structure with a single conformation. Specifically, it yielded spectra with good dispersion of chemical shifts in the amide region ranging from 7.7 to 9.1 ppm and exactly the expected number of TOCSY spin systems. In contrast, the spectra of globular BuIA displayed more than the expected number of TOCSY cross-peaks, indicating the likely presence of conformational heterogeneity. The possibility that the extra cross-peaks were due to impurities was eliminated because only a single sharp HPLC peak was observed and only a single mass was detected by mass spectrometry.

thumbnailFigure 2. NMR analysis of ribbon and globular BuIA. The upper panels show the amide regions of the TOCSY spectra at 600 MHz of globular (A) and ribbon BuIA (B) at 282 K in 95% H2O/5% 2H2O at pH 3. The one letter code and the residue numbers are used for labeling. For globular BuIA (A) stretches of assigned residues are labeled, as a or b to indicate conformers A and B. The lower panels (C & D) show a secondary shift analysis of conformers A and B from globular BuIA and ribbon BuIA compared to AuIB. For each residue, represented by the one letter code, the differences between the actual chemical shifts for the α-protons and their respective random coil shifts are graphed. (C) Native AuIB (black bars), conformer A (white bars), conformer B (grey bars). (D) Ribbon AuIB (black bars), ribbon BuIA (white bars).

The sequential assignment of the spectra of the ribbon and globular forms was undertaken using a combination of TOCSY and NOESY spectra [19]. For the ribbon isomer the assignment was straightforward. By contrast, the presence of multiple sets of peaks due to three conformers (referred to here as conformers A, B and C), a lack of detection of sufficient sequential NOESY cross-peaks and peak overlap prevented complete sequential assignment of all of the globular BuIA conformers. Nonetheless, with the aid of the 900 MHz spectra conformers A and B were assigned, as indicated in Fig 2A. The first six residues of conformer C were assigned but the proline residue at position 6 has weak spectral peaks and has an αH shift overlapped with several other residues and thus assignment of Pro7 and sequentially following residues was not possible. A range of pH values (3, 4 and 5.5), temperatures 280–310 K and acetonitrile co-solvent concentrations (20% or 40% v/v) were used to determine if solution conditions influenced the presence of multiple conformations but similar ratios of isomers were present in all the conditions used.

Conformers A and C both display sharp peaks, in contrast to conformer B, which has significant broadening. Based on peak integrals the ratio of the isomers is approximately 3:2:1 (A:B:C). Because peaks from conformers A and C in the same solution are sharp, the broadened peaks of conformer B presumably do not reflect aggregation, but are likely to be the result of exchange broadening. Thus, not only does the globular isomer lead to at least three conformers in slow exchange on the NMR time scale, but one set of peaks probably reflects the presence of two or more conformers in intermediate to fast exchange. Overall this isomer of BuIA appears to be much more conformationally heterogeneous than other conotoxins studied to date.

Inspection of the sequential cross-peaks associated with the proline residues indicated that the conformational heterogeneity observed for globular BuIA originates from cis-trans isomerisation of the two proline residues. Conformer A has both proline residues in a trans conformation, consistent with the report of Chi et al. [20] whereas conformer B has a trans Pro6 and cis Pro7 based on strong NOEs observed between the α- and δ-protons of Thr5 and Pro6 and a strong αH-αH NOE between Pro6 and Pro7 respectively. The lack of a complete assignment for conformer C and peak overlap prevented determination of the proline configurations for this conformation. In the ribbon form of BuIA, which adopts an ordered structure, the same proline conformation pattern as conformer B was observed, namely trans Pro6 and cis Pro7.

Analysis of the secondary αH chemical shifts of ribbon BuIA indicates a lack of regular secondary structure. There is no continuous stretch of residues exhibiting negative deviations that would be indicative of helices or positive deviations indicative of extended structures, such as β-sheets. Panel D of Fig. 2 shows a comparison of the αH secondary shifts of ribbon BuIA with the ribbon isomer of AuIB [21]. The trends are similar, but in general the secondary shifts of ribbon BuIA deviate from random coil to a greater extent than that observed for ribbon AuIB consistent with a well-defined structure for ribbon BuIA.

Despite the lack of complete assignment of all the conformations of globular BuIA an analysis of the chemical shifts was very informative. The secondary shifts of conformer A are very similar to the related α-conotoxin AuIB [21] as shown in Fig. 2C, indicating that the native conotoxin fold is present [18]. In contrast, the secondary shifts of conformer B differ significantly from those of AuIB (Fig. 2D). While our work was in progress the structure of the native conformation (i.e. conformer A) of BuIA was reported by Chi et al and confirms that the consensus structure of the α-conotoxins is maintained [20]. It was not possible to directly compare our findings with those of Chi et al because chemical shifts or spectra were not provided in the study but consistent with our study, Chi et al reported the presence of additional peaks that they presumed arose from a cis conformation in the globular isomer. In the present study, we have shown there is more than one additional conformation present. The ribbon isomer, which is the focus of our study, was not examined in the Chi et al study [20].

Given the single conformation observed for ribbon BuIA it was of interest to determine the three dimensional structure of this isomer. Structure calculations were carried out with 96 NOE distance restraints (comprising 12 long-range, 18 medium-range, 50 sequential and 16 intraresidual distance restraints) and six dihedral angle restraints, using a simulated annealing protocol in CNS. The ϕ-angle restraints, derived from αN coupling constants, were -65 ± 15° for Ala9 and Leu11, -120 ± 15° for Thr5, and -120° ± 30° for Cys8. The ϕ-angle for Cys13 was restrained to 100 ± 80° based on the intraresidual Hα-HN NOE being clearly weaker than the NOE between HN and the Hα of the preceding residue. The χ1-angle of Cys8 was restrained to 180 ± 30° based on the αβ coupling constants and NOE intensities. Hydrogen bonds were not included in the calculations but measurement of the exchange rates for the amide protons of ribbon BuIA revealed that 3 of the 11 amides, namely Val10, Leu11 and Tyr12, are slowly exchanging. Their amide proton signals were still detected in NMR spectra recorded 30 minutes after dissolution of the sample in 2H2O, but all amide protons had exchanged after 45 minutes in 2H2O. This result suggests that Val10, Leu11 and Tyr12 are more protected from the solvent than are the other residues.

Fig. 3 shows a stereoview of the superposition of the 20 lowest energy structures of ribbon BuIA. Despite its small size it adopts a well-defined three-dimensional structure. The ensemble of structures aligns well over the whole molecule with a pairwise backbone RMSD of 0.36 ± 0.12 Å and a RMSD of 1.1 ± 0.49 Å for all heavy atoms, and had good structural and energetic statistics, as shown in Table 2. No NOE violations higher than 0.3 Å and no dihedral angle restraint violations greater than 3° were observed within these 20 structures. In the Ramachandran plot, showing the favored ϕ- and ψ-angle combinations for residues in proteins, all residues of ribbon BuIA (except Pro and Gly) lie in the most favored and additionally allowed regions.

Table 2. NMR and refinement statistics for ribbon BuIA

thumbnailFigure 3. Three-dimensional structure of the ribbon isomer of BuIA. (top) Stereoview of the 20 NMR-derived lowest energy structures of ribbon BuIA. For clarity only the backbone (N, Cα, C') is shown and the two disulfide bonds are omitted. The structures are superimposed over the whole molecule. N refers to the N-terminus; C refers to the C-terminus. Residues were labeled with the one letter code. (bottom) Solvent accessible surface of ribbon BuIA. The two views are rotated by 180° about the vertical axis. Hydrophobic residues (Pro, Ala, Val, Leu, Tyr) are green, hydrophilic residues (Ser, Thr) are blue, Gly is grey and the Cys residues are yellow. Selected residues are labeled with the one letter code.

Using the program PROMOTIF [22] the 20 lowest energy structures of ribbon BuIA were examined to identify consensus structural elements. In agreement with the prediction from the secondary shifts, ribbon BuIA shows no α-helices or β-sheets. The structure comprises an inverse γ-turn for residues 2–4 but no regular secondary structure is present in the remainder of the molecule. One of the two disulfide bonds stabilizing the three dimensional structure, namely that formed by Cys3 and Cys8, overlays very well in the 20 lowest energy structures. In 17 of 20 structures the torsion angles of this disulfide bond are representative of a right-handed hook [23]. By contrast, the disulfide bond joining Cys2 to Cys13 is classified as right-handed hook in 9 out of 20 structures but not recognized as a standard disulfide type in the other structures. Calculation of the hydrogen bonds using MolMol [24] revealed that in 13 of 20 structures the amide proton of Leu11 is bonded to the carbonyl of Ala9. This hydrogen bond is consistent with the relatively slow exchange of the amide proton of Leu11, but Val10 and Tyr12 were also slowly exchanging and no hydrogen bond acceptors were apparent from the calculated structures.

Using the program InsightII the solvent accessible surface of ribbon BuIA was calculated and is shown in Fig. 3, from which the large hydrophobic surface formed by Ala, Val, Leu, Tyr and Pro residues is clear. The hydrophobic patch extends across a wide stretch of the whole molecule, as it comprises only a few hydrophilic residues, namely Ser and Thr. Strikingly, there are no charged residues in the molecule. Tyr12 protrudes conspicuously, making the overall shape of ribbon BuIA non-spherical. In contrast to many other α-conotoxins the cysteine residues are rather exposed to the solvent, presumably as a result of the ribbon disulfide connectivity and the fact that BuIA is only 13 residues in length.

The biological activities of the BuIA globular and ribbon analogues, as well as a form where the cysteine residues were alkylated, were examined in an electrophysiological assay by measuring the % inhibition of ACh-evoked current amplitude at the α3β2 and α3β4 nAChR subtypes expressed in Xenopus oocytes (Fig. 4A). The alkylated form of BuIA was analyzed by mass spectrometry and had an observed molecular mass (1571.9 Da) consistent with the calculated mass (1571.64 Da). Membrane currents were evoked with 100 μM ACh and concentration-response curves were fitted to determine IC50 values and Hill coefficients. The concentration-response curves for inhibition of α3β2 and α3β4 nAChR subtypes by globular BuIA gave IC50 values of 4.8 ± 0.4 nM (nH = 1.3) and 59.1 ± 2.3 nM (nH = 1.2), respectively, whereas almost no inhibition was observed with 1 μM ribbon BuIA at either nAChR subtype (n = 4–7; Fig. 4B). Furthermore, the globular BuIA isomer is approximately ten-fold more potent at β2- compared to β4-subunit containing nAChR subtypes. These results for synthetic globular BuIA are consistent with that reported previously [16]. Alkylated BuIA was inactive at α3β2 and α3β4 subtypes when tested at 10 nM and 100 nM respectively. Increasing the concentration of alkylated BuIA to 10 μM also failed to inhibit nAChRs.

thumbnailFigure 4. Effect of globular and ribbon isomers of α-conotoxin BuIA at nAChR subtypes expressed in Xenopus oocytes. (A)Representative ACh-evoked currents mediated by α3β2 and α3β4 nAChR subtypes obtained in the absence (control) and presence of 10 nM and 100 nM BuIA-globular, respectively. Complete recovery was observed after 10 min washout (broken line). (B)Concentration-response curves for the inhibition of α3β2 and α3β4 nAChRs by globular (filled symbols) and ribbon (open symbols) isomers of BuIA. Best fit of the data gave IC50 values of 4.8 ± 0.4 nM (nH = 1.3) and 59.1 ± 2.3 nM (nH = 1.2) for α3β2 and α3β4 nAChRs, respectively, whereas the ribbon isomer of BuIA exhibited no inhibition at 1 μM. Responses are shown as a percentage of ACh (100 μM)-induced peak current amplitude after a 5 min incubation of the BuIA isomers with respect to control (ACh alone). Error bars are SEM with n = 4–7 for each data point.

Discussion

As BuIA is the only α-conotoxin reported so far with a 4/4 spacing its three dimensional structure is of significant interest. In the current study we examined the structural features and determined the biological activity of the globular (CysI–CysIII, CysII–CysIV) and ribbon (CysI–CysIV, CysII–CysIII) disulfide isomers of BuIA. These studies show that the general trend observed for α-conotoxins where the native globular form has a well-defined conformation and the ribbon isomer is disordered is reversed for BuIA.

NMR analysis of the globular isomer of BuIA revealed more cross-peaks than expected for a 13 residue peptide. At least three conformations are present in the two-dimensional spectra of globular BuIA and as native α-conotoxin structures determined to date generally display well-defined structures, the significant heterogeneity differs from usual findings for α-conotoxins. Surprisingly, this conformational heterogeneity is not seen in the ribbon isomer, which forms a single well-defined conformation. The observation that the HPLC trace of the ribbon reveals a broad peak, but the structure is well-defined is interesting and unusual. One possibly is that the different conditions under which the HPLC and NMR are carried out are responsible for these differences. The more hydrophobic environment used for the HPLC analysis coupled with the presence of TFA could potentially affect how the ribbon isomer behaves on RP-HPLC. However, recording NMR spectra of globular BuIA in the presence of varying amounts of acetonitrile co-solvent did not influence the conformational heterogeneity, making the sharper peak observed by HPLC analysis for globular BuIA compared to ribbon BuIA a puzzling result.

The conformational heterogeneity of globular BuIA originates form cis/trans isomerisation of the pair of adjacent proline residues in the middle of the sequence. The proline configurations for the globular isomers and ribbon BuIA are highlighted in Fig. 5. Conformer A has both proline residues in a trans conformation, consistent with the trend observed for other native α-conotoxins, whereas conformer B is broadened and displays a trans/cis configuration for Pro6 and Pro7 respectively. The conformations of the proline residues in conformer C were not determined because of the lack of a complete assignment. Ribbon BuIA displays the same proline geometry as conformer B despite the difference in disulfide connectivity. The BuIA sequence can clearly accommodate a range of geometries of the proline residues in contrast to what has generally been observed in previous structural analyses of α-conotoxins [17,18]. One exception is α-conotoxin GI that has previously been shown to have two interconverting conformations [25] but this differs from BuIA where there are clearly multiple conformations. Other peptides are known to contain similar proline motifs to BuIA, including SFTI-1[26], and in that case the Pro-Pro motif has a cis/trans configuration as shown in Fig. 5. This highlights the complexity that can be associated with the geometry of proline residues in small cysteine-rich peptides and how changes in structure and sequence can result in different geometries being favored.

thumbnailFigure 5. Proline geometries of the BuIA isomers and the sunflower trypsin inhibitor SFTI-1 [26]. The proline residues with trans geometry are marked as t and those with cis as c. The three-dimensional structures of a classic α-conotoxin, ribbon BuIA and SFTI-1 are shown to highlight the structural differences between the various peptides.

The structurally heterogeneous globular isomer of BuIA is active on α3-containing neuronal nAChRs, with a preference for receptors comprising β2- rather than β4-subunits. This specificity is quite similar to some of the 4/7 α-conotoxins such as MII [20]. The structure of conformer A of globular BuIA, determined by Chi et al [20], has been compared to that of MII and it appears that the two-turn helical motif present in both MII and BuIA is important for binding to the α3-subunit of neuronal nAChRs [20]. Recent studies on the binding of α-conotoxins to the acetylcholine binding protein, a homolog of the nAChR extracellular domain, have shown that the rigid consensus structure found in solution does not change significantly upon binding [27-29]. Thus, as ribbon BuIA does not contain the consensus structural motif it consequently exhibits no inhibitory activity at nAChRs up to 1 μM. It appears that the well-defined conformation of ribbon BuIA prevents it from adopting an appropriate conformation for binding to the nAChR.

We have also shown here that an analogue of BuIA with the cysteine residues reduced and alkylated is inactive at the nAChR. This is consistent with studies on ImI where the reduced and alkylated form had a significant decrease in activity at the α7 receptor [30]. Given that the reduced and alkylated form of BuIA is likely to be devoid of structure based on the strong influence of disulfide bonds on the structures of small cysteine-rich peptides, the lack of activity of this peptide indicates that it cannot adopt an appropriate conformation for binding to the nAChR. Thus, the linear sequence alone is not sufficient for activity and a propensity to form a helical structure in solution appears to be correlated with activity at the nAChR.

The ribbon isomer of BuIA, despite having a well-defined structure, is not the preferred isomer formed during air oxidation using standard conditions applicable to other conotoxins. The fact that the folding pathway favors the globular disulfide connectivity indicates that it is energetically more favorable, even though multiple conformations are present in solution. Although the ribbon isomer is less favoured during oxidation, as is often the case for non-native isomers, it displays significantly different properties to both non-native and native conotoxin ribbon isomers. The ribbon isomers of both α-AuIB and α-GI exhibit greater conformational flexibility than their globular counterparts [21,31], in contrast to ribbon BuIA, indicating that the unique sequence of BuIA stabilizes the ribbon isomer. Structures of globular and ribbon AuIB are shown in Fig. 6. Although a single structure is shown for the ribbon isomer of AuIB it displays considerable conformational heterogeneity relative to the globular isomer. Furthermore, the ribbon isomer of BuIA lacks any regular secondary structure, in contrast to the χ-conotoxins that naturally contain a ribbon connectivity. The χ-conotoxin MrIA consists of a β-hairpin connected by an inverse γ-turn [32], as shown in Fig. 6.

thumbnailFigure 6. Backbone trace of the structures of ribbon BuIA (A), globular AuIB (B), ribbon AuIB (C), and ribbon MrIA (D). Disulfide bonds are shown as balls and sticks. The helices are shown with thickened ribbons and β-strands are indicated with arrows.

Another similarity between the ribbon isomers of α-conotoxins is related to the geometry of the proline residues [6,18]. The inactive ribbon isomers of BuIA and GI [31] both contain a cis proline (Pro7 and Pro5 respectively) in contrast to the trans Pro seen in the native disulfide connectivity for all α-conotoxin structures determined to date [18]. In one conformation of globular BuIA (conformer B, Fig. 2) Pro7 is also in a cis geometry, but the alternative trans conformation may be stabilized upon binding to the receptor. The implications for the presence of the cis proline in activity have not yet been elucidated.

Conclusion

The native globular isomer of BuIA, although being the favored isomer during oxidative refolding, has multiple conformations in solution, unlike the majority of native α-conotoxins. In addition, the structure of the ribbon isomer is stabilized whereas all other non-native isomers display disordered structures. Overall our results highlight the influence of the disulfide connectivity of BuIA on the dynamics of the three-dimensional structure. This information is potentially important for on-going efforts to understand the structure-activity relationships of this valuable class of peptides.

Methods

Materials

Protected N-α-Boc-L-amino acids and reagents used for peptide chain assembly were purchased from NovaBiochem. Methylbenzhydrylamine resin was from Peptide Institute Inc. (Osaka, Japan). p-Cresol and p-thiocresol were purchased from Aldrich (Sydney, Australia). HBTU was purchased form Richelieu Biotechnologies (Quebec, Canada). Anhydrous HF was from Matheson Gas (BOC Gases, Melbourne, Australia). All of the solvents used for HPLC and peptide synthesis were purchased from Labscan, (Bangkok, Thailand).

Peptide synthesis

BuIA analogues were manually assembled using a stepwise in situ neutralization protocol for Boc chemistry [33]. The four Cys residues were differentially protected in pairs corresponding to the globular and ribbon isomers with HF-labile methylbenzyl groups and HF-resistant acetamidomethyl groups. An additional synthesis was also done without selective protection of the cysteine residues. Other amino acid side-chain protection was as follows: Tyr (BrZ), Thr (Bzl) and Ser (Bzl). All the syntheses were carried out on a 0.25 mmol scale using HBTU (0.5 M/DMF) as the activation reagent, and DIEA as neutralizing reagent on MBHA resin. Neat trifluoroacetic acid was used as the amine-deprotecting reagent. Deprotection and cleavage were carried out in HF:p-cresol:p-thiocresol (9:0.5:0.5) for 2 h. Methylbenzyl groups attached to cysteines were removed during the HF cleavage and the oxidation of the first disulfide bond was performed in 0.1 M NH4HCO3 oxidation buffer at pH 8 at a concentration of 0.1 mg/mL. The one pair disulfide bonded peptides were purified by HPLC and lyophilized prior to the formation of the second disulfide bond. Simultaneous deprotection and oxidation of acetamidomethyl protected cysteine residues was performed in 80% methanol at low pH. Iodine was added and the reaction mixture was stirred for 5 min under argon. Methanol and iodine were removed by rotary evaporation. The resultant mixture was diluted 10-fold and purified by RP-HPLC. The peptide synthesized without selective protection of the cysteine residues was oxidized in 0.1 M NH4HCO3 pH 8, and 1 mM reduced glutathione. Linear reduced BuIA was alkylated with iodoacetamide in 0.1 M ammonium acetate buffer for testing in biological assays.

RP-HPLC and mass spectrometry

Preparative HPLC were performed on Waters 600E solvent delivery system. Data were collected via a 484-absorbance detector from Waters at 230 nm. Analytical HPLC was performed on Shimadzu LC-2010 instruments. Preparative HPLC was carried out using a Phenomenex C18 column (250 × 21.20 mm, 10 μm). Analytical HPLC was performed on a Phenomenex C18 (250 × 4.6 mm, 10 μm) column. 90% acetonitrile and 0.043% trifluoroacetic acid was used as the eluting solvent B. 0.05% trifluoroacetic acid was used as the eluting solvent A. 1%B per minute was used as the linear gradient at standard column flow rates (8 mL/min for preparative HPLC and 1 mL/min for analytical HPLC). Electrospray mass spectra were acquired on a Micromass LCT (Manchester, UK), a liquid chromatography-orthogonal acceleration reflecting TOF mass spectrometer, coupled to an Agilent HPLC system. Samples (~5 μL) were injected into a moving solvent (70 μL/min, 70% acetonitrile/formic acid 0.05% in water) coupled directly to the electrospray ionisation source. Full scan mass spectra were acquired over the mass range of 400–1600 Da with a scan step size of 0.2 Da. Molecular masses were derived from the observed m/z values.

NMR spectroscopy

Although X-ray crystallography has been used in a few cases to solve conotoxin structures [34,35] NMR has been the predominant technique for structurally characterizing this class of peptides[17]. Samples for 1H NMR measurements contained ~1 mM peptide in 95% H2O/5% D2O (v/v) at pH 3, 4 and 5.5. Spectra on ribbon BuIA were recorded at 282 or 290 K on Bruker AVANCE-500, 600 and 900 spectrometers whereas a range of temperatures between 280 and 310 K were used for globular BuIA to determine if temperature influenced the multiple conformations present. 2D NMR spectra were recorded in phase-sensitive mode using time-proportional phase incrementation for quadrature detection in the t1 dimension. The 2D experiments consisted of a TOCSY using a MLEV-17 spin lock sequence with a mixing time of 80 ms, DQF-COSY, ECOSY and NOESY with mixing times of 100–250 ms. Solvent suppression was achieved using a modified WATERGATE sequence. Spectra were acquired over 6024 Hz with 4096 complex data points in F2 and 512 increments in the F1 dimension. Identification of slowly exchanging amide protons was achieved by acquisition of one-dimensional and TOCSY spectra immediately following dissolution of a fully protonated peptide in 2H23JHN-Hα coupling constants were measured from a one-dimensional spectrum or from the DQF-COSY spectrum.

Spectra were processed on a Silicon Graphics Indigo workstation using XWINNMR (Bruker) software. The t1 dimension was zero-filled to 1024 real data points, and 90° phase-shifted sine bell window functions were applied prior to Fourier transformation. Chemical shifts were referenced to internal 2, 2-dimethyl-2-silapentane-5-sulfonate.

Structure calculations

Preliminary structures of ribbon BuIA were calculated using a torsion angle simulated annealing protocol within the program DYANA [36]. Final structures were calculated using simulated annealing and energy minimization protocols within CNS version 1.1 [37]. The starting structures were generated using random (ϕ,ψ) dihedral angles and energy-minimized to produce structures with the correct local geometry. A set of 50 structures was generated by a torsion angle simulated annealing protocol [38]. This protocol involved a high-temperature phase comprising 4000 steps of 0.015 ps of torsion angle dynamics, a cooling phase with 4000 steps of 0.015 ps of torsion angle dynamics during which the temperature was lowered to 0 K, and finally an energy minimization phase comprising 500 steps of Powell minimization. Structures consistent with restraints were subjected to further molecular dynamics and energy minimization in a water shell, as described by Linge and Nilges [39]. The refinement in explicit water involved the following steps. First, heating to 500 K via steps of 100 K, each comprising 50 steps of 0.005 ps of Cartesian dynamics. Second, 2500 steps of 0.005 ps of Cartesian dynamics at 500 K before a cooling phase where the temperature was lowered in steps of 100 K, each comprising 2500 steps of 0.005 ps of Cartesian dynamics. Finally, the structures were minimized with 2000 steps of Powell minimization. Structures were analyzed using PROMOTIF [22] and PROCHECK-NMR [40].

Electrophysiology

RNA was prepared from linearized cDNA encoding the rat α3,β2 and β4 nAChR subunits (provided by Dr. J. Patrick, Baylor College of Medicine, Houston, TX) and transcribed using a mMESSAGE mMACHINE™ high yield capped RNA transcription kit. Xenopus laevis were anesthetized with oocytes surgically removed and placed in OR2 buffer (82 mM NaCl, 2 mM KCl, 1 mM MgCl2 and 5 mM HEPES at pH 7.4) with 3 mg/ml collagenase (Sigma) for 1–2 hrs at room temperature then stored in ND96 buffer (96 mM NaCl, 2 mM KCl, 1.8 mM CaCl2, 1 mM MgCl2 and 5 mM HEPES at pH 7.4) supplemented with 50 mg/l gentamycin and 5 mM pyruvic acid. Oocytes were injected with 5 ng of cRNA and kept at 18°C in ND96 buffer for 2–5 days before recording.

Membrane currents were recorded from Xenopus oocytes using the two electrode (virtual ground circuit) voltage clamp technique with an automated workstation with eight channels in parallel, and an on-line analysis (OpusXpress™ 6000A workstation, Molecular Devices, Union City, CA). Both the voltage recording and current injecting electrodes were pulled from borosilicate glass (GC150T-15; Harvard Apparatus, Edenbridge, UK) and filled with 3 M KCl with resistances between 0.2 – 1.5 MΩ. All recordings were conducted at room temperature (20 – 23°C) using a bath solution of ND96 as described above. During recordings, the oocytes were perfused continuously at a rate of 1.5 ml/min. Acetylcholine (100 μM) was applied for 2 s at 5 ml/min with 10 min washout periods between successive applications with 5 min incubations of the BuIA analogues. Cells were voltage clamped at a holding potential of -80 mV. Data were sampled at 500 Hz and low pass filtered at 200 Hz. Peak current amplitude was measured before and after the incubation of BuIA analogues, using the empirical Hill equation to calculate the IC50 and Hill coefficient (nH) of the concentration-response curves. The data were fitted by a non-linear least squares algorithm (SPSS inc. Chicago, IL) with pooled data represented as mean ± SEM.

Abbreviations

HPLC, high performance liquid chromatography; NMR, nuclear magnetic resonance; ACh, acetylcholine; nAChRs, nicotinic acetylcholine receptors.

Authors' contributions

A-HJ designed experiments and carried out synthetic studies. HB carried out structural studies and helped draft the manuscript. STN designed and carried out biological assays. CCT carried out synthetic and folding studies. RJC designed and carried out synthetic and folding studies. DJA designed biological experiments and revised the manuscript. PFA conceived study and designed synthetic experiments. DJC designed structural studies and revised the manuscript. NLD designed and carried out structural studies and drafted and revised the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

DJC and PFA are Senior Research Affiliates of the Australian Research Council (ARC) Special Research Centre for Functional and Applied Genomics, which provided infrastructure support for the work in this project. DJC is an ARC Professorial Fellow. NLD is an NHMRC Industry Fellow. This work was supported by a grant from the ARC. This work was presented in part at the IST World Congress in Glasgow July 2006.

References

  1. Role LW, Berg DK: Nicotinic receptors in the development and modulation of CNS synapses.

    Neuron 1996, 16(6):1077-1085. PubMed Abstract | Publisher Full Text OpenURL

  2. Jones S, Sudweeks S, Yakel JL: Nicotinic receptors in the brain: correlating physiology with function.

    Trends Neurosci 1999, 22(12):555-561. PubMed Abstract | Publisher Full Text OpenURL

  3. Mazurov A, Hauser T, Miller CH: Selective alpha7 nicotinic acetylcholine receptor ligands.

    Curr Med Chem 2006, 13(13):1567-1584. PubMed Abstract | Publisher Full Text OpenURL

  4. Dutertre S, Lewis RJ: Toxin insights into nicotinic acetylcholine receptors.

    Biochem Pharmacol 2006, 72(6):661-670. PubMed Abstract | Publisher Full Text OpenURL

  5. Lustig LR: Nicotinic acetylcholine receptor structure and function in the efferent auditory system.

    Anat Rec A Discov Mol Cell Evol Biol 2006, 288(4):424-434. PubMed Abstract | Publisher Full Text OpenURL

  6. McIntosh JM, Gardner S, Luo S, Garrett JE, Yoshikami D: Conus peptides: novel probes for nicotinic acetylcholine receptor structure and function.

    Eur J Pharmacol 2000, 393(1-3):205-208. PubMed Abstract | Publisher Full Text OpenURL

  7. Adams DJ, Alewood PF, Craik DJ, Drinkwater RD, Lewis RJ: Conotoxins and their potential pharmaceutical applications.

    Drug Dev Res 1999, 46:219-234. Publisher Full Text OpenURL

  8. Livett BG, Sandall DW, Keays D, Down J, Gayler KR, Satkunanathan N, Khalil Z: Therapeutic applications of conotoxins that target the neuronal nicotinic acetylcholine receptor.

    Toxicon 2006, 48(7):810-829. PubMed Abstract | Publisher Full Text OpenURL

  9. McIntosh JM, Santos AD, Olivera BM: Conus peptides targeted to specific nicotinic acetylcholine receptor subtypes.

    Ann Rev Biochem 1999, 68:59-88. PubMed Abstract | Publisher Full Text OpenURL

  10. Dutton JL, Craik DJ: alpha-Conotoxins: nicotinic acetylcholine receptor antagonists as pharmacological tools and potential drug leads.

    Curr Med Chem 2001, 8(4):327-344. PubMed Abstract | Publisher Full Text OpenURL

  11. Balaji RA, Ohtake A, Sato K, Gopalakrishnakone P, Kini RM, Seow KT, Bay BH: lambda-conotoxins, a new family of conotoxins with unique disulfide pattern and protein folding. Isolation and characterization from the venom of Conus marmoreus.

    J Biol Chem 2000, 275(50):39516-39522. PubMed Abstract | Publisher Full Text OpenURL

  12. Sharpe IA, Gehrmann J, Loughnan ML, Thomas L, Adams DA, Atkins A, Palant E, Craik DJ, Adams DJ, Alewood PF, Lewis RJ: Two new classes of conopeptides inhibit the alpha1-adrenoceptor and noradrenaline transporter.

    Nat Neurosci 2001, 4(9):902-907. PubMed Abstract | Publisher Full Text OpenURL

  13. Gray WR, Luque A, Olivera BM, Barrett J, Cruz LJ: Peptide toxins from Conus geographus venom.

    J Biol Chem 1981, 256(10):4734-4740. PubMed Abstract | Publisher Full Text OpenURL

  14. Park KH, Suk JE, Jacobsen R, Gray WR, McIntosh JM, Han KH: Solution conformation of alpha-conotoxin EI, a neuromuscular toxin specific for the alpha 1/delta subunit interface of torpedo nicotinic acetylcholine receptor.

    J Biol Chem 2001, 276(52):49028-49033. PubMed Abstract | Publisher Full Text OpenURL

  15. Martinez JS, Olivera BM, Gray WR, Craig AG, Groebe DR, Abramson SN, McIntosh JM: alpha-Conotoxin EI, a new nicotinic acetylcholine receptor antagonist with novel selectivity.

    Biochemistry 1995, 34(44):14519-14526. PubMed Abstract | Publisher Full Text OpenURL

  16. Azam L, Dowell C, Watkins M, Stitzel JA, Olivera BM, McIntosh JM: alpha-conotoxin BuIA, a novel peptide from Conus bullatus, distinguishes among neuronal nicotinic acetylcholine receptors.

    J Biol Chem 2005, 280(1):80-87. PubMed Abstract | Publisher Full Text OpenURL

  17. Marx UC, Daly NL, Craik DJ: NMR of conotoxins: structural features and an analysis of chemical shifts of post-translationally modified amino acids.

    Magn Reson Chem 2006, 44 Spec No:S41-50. PubMed Abstract | Publisher Full Text OpenURL

  18. Millard EL, Daly NL, Craik DJ: Structure-activity relationships of alpha-conotoxins targeting neuronal nicotinic acetylcholine receptors.

    Eur J Biochem 2004, 271(12):2320-2326. PubMed Abstract | Publisher Full Text OpenURL

  19. Wüthrich K: NMR of Proteins and Nucleic Acids. New York , Wiley-Interscience; 1986.

  20. Chi SW, Kim DH, Olivera BM, McIntosh JM, Han KH: NMR structure determination of alpha-conotoxin BuIA, a novel neuronal nicotinic acetylcholine receptor antagonist with an unusual 4/4 disulfide scaffold.

    Biochem Biophys Res Commun 2006, 349(4):1228-1234. PubMed Abstract | Publisher Full Text OpenURL

  21. Dutton JL, Bansal PS, Hogg RC, Adams DJ, Alewood PF, Craik DJ: A new level of conotoxin diversity, a non-native disulfide bond connectivity in alpha-conotoxin AuIB reduces structural definition but increases biological activity.

    J Biol Chem 2002, 277(50):48849-48857. PubMed Abstract | Publisher Full Text OpenURL

  22. Hutchinson EG, Thornton JM: PROMOTIF - A program to identify and analyze structural motifs in proteins.

    Protein Sci 1996, 5:212-220. PubMed Abstract | Publisher Full Text OpenURL

  23. Richardson JS: The anatomy and taxonomy of protein structure.

    Adv Protein Chem 1981, 34:167-339. PubMed Abstract OpenURL

  24. Koradi R, Billeter M, Wüthrich K: MOLMOL: a program for display and analysis of macromolecular structures.

    J Mol Graph 1996, 14(1):51-55. PubMed Abstract | Publisher Full Text OpenURL

  25. Maslennikov IV, Sobol AG, Gladky KV, Lugovskoy AA, Ostrovsky AG, Tsetlin VI, Ivanov VT, Arseniev AS: Two distinct structures of alpha-conotoxin GI in aqueous solution.

    Eur J Biochem 1998, 254(2):238-247. PubMed Abstract | Publisher Full Text OpenURL

  26. Luckett S, Garcia RS, Barker JJ, Konarev AV, Shewry PR, Clarke AR, Brady RL: High-resolution structure of a potent, cyclic proteinase inhibitor from sunflower seeds.

    J Mol Biol 1999, 290(2):525-533. PubMed Abstract | Publisher Full Text OpenURL

  27. Celie PHN, Kasheverov IE, Mordvintsev DY, Hogg RC, van Nierop P, van Elk R, van Rossum-Fikkert SE, Zhmak MN, Bertrand D, Tsetlin V, Sixma TK, Smit AB: Crystal structure of nicotinic acetylcholine receptor homolog AChBP in complex with an alpha-conotoxin PnIA variant.

    Nat Struct Mol Biol 2005, 12(7):582-588. PubMed Abstract | Publisher Full Text OpenURL

  28. Hansen SB, Sulzenbacher G, Huxford T, Marchot P, Taylor P, Bourne Y: Structures of Aplysia AChBP complexes with nicotinic agonists and antagonists reveal distinctive binding interfaces and conformations.

    Embo J 2005, 24(20):3635-3646. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Ulens C, Hogg RC, Celie PH, Bertrand D, Tsetlin V, Smit AB, Sixma TK: Structural determinants of selective alpha-conotoxin binding to a nicotinic acetylcholine receptor homolog AChBP.

    Proc Natl Acad Sci U S A 2006, 103(10):3615-3620. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Servent D, Thanh HL, Antil S, Bertrand D, Corringer PJ, Changeux JP, Menez A: Functional determinants by which snake and cone snail toxins block the alpha 7 neuronal nicotinic acetylcholine receptors.

    J Physiol Paris 1998, 92(2):107-111. PubMed Abstract | Publisher Full Text OpenURL

  31. Gehrmann J, Alewood PF, Craik DJ: Structure determination of the three disulfide bond isomers for a-conotoxin GI: A model for the role of disulfide bonds in structural stability.

    J Mol Biol 1998, 278:401-415. PubMed Abstract | Publisher Full Text OpenURL

  32. Nilsson KP, Lovelace ES, Caesar CE, Tynngard N, Alewood PF, Johansson HM, Sharpe IA, Lewis RJ, Daly NL, Craik DJ: Solution structure of chi-conopeptide MrIA, a modulator of the human norepinephrine transporter.

    Biopolymers 2005, 80(6):815-823. PubMed Abstract | Publisher Full Text OpenURL

  33. Schnölzer M, Alewood P, Jones A, Alewood D, Kent SBH: In situ neutralization in Boc-chemistry solid phase peptide synthesis.

    Int J Pept Protein Res 1992, 40:180-193. PubMed Abstract OpenURL

  34. Hu SH, Gehrmann J, Guddat LW, Alewood PF, Craik DJ, Martin JL: The 1.1 A crystal structure of the neuronal acetylcholine receptor antagonist, a-conotoxin PnIA from Conus pennaceus.

    Structure 1996, 4:417-423. PubMed Abstract | Publisher Full Text OpenURL

  35. Hu SH, Gehrmann J, Alewood PF, Craik DJ, Martin JL: Crystal structure at 1.1 A resolution of alpha-conotoxin PnIB: comparison with alpha-conotoxins PnIA and GI.

    Biochemistry 1997, 36(38):11323-11330. PubMed Abstract | Publisher Full Text OpenURL

  36. Guntert P, Mumenthaler C, Wüthrich K: Torsion angle dynamics for NMR structure calculation with the new program DYANA.

    J Mol Biol 1997, 273(1):283-298. PubMed Abstract | Publisher Full Text OpenURL

  37. Brünger AT, Adams PD, Rice LM: New applications of simulated annealing in X-ray crystallography and solution NMR.

    Structure 1997, 5(3):325-336. PubMed Abstract | Publisher Full Text OpenURL

  38. Rice LM, Brünger AT: Torsion angle dynamics: reduced variable conformational sampling enhances crystallographic structure refinement.

    Proteins 1994, 19(4):277-290. PubMed Abstract | Publisher Full Text OpenURL

  39. Linge JP, Nilges M: Influence of non-bonded parameters on the quality of NMR structures: a new force field for NMR structure calculation.

    J Biomol NMR 1999, 13(1):51-59. PubMed Abstract | Publisher Full Text OpenURL

  40. Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM: AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR.

    J Biomol NMR 1996, 8(4):477-486. PubMed Abstract | Publisher Full Text OpenURL

  41. McIntosh JM, Yoshikami D, Mahe E, Nielsen DB, Rivier JE, Gray WR, Olivera BM: A nicotinic acetylcholine receptor ligand of unique specificity, alpha-conotoxin ImI.

    J Biol Chem 1994, 269(24):16733-16739. PubMed Abstract | Publisher Full Text OpenURL

  42. Fainzilber M, Hasson A, Oren R, Burlingame AL, Gordon D, Spira ME, Zlotkin E: New mollusc-specific alpha-conotoxins block Aplysia neuronal acetylcholine receptors.

    Biochemistry 1994, 33(32):9523-9529. PubMed Abstract | Publisher Full Text OpenURL

  43. Luo S, Kulak JM, Cartier GE, Jacobsen RB, Yoshikami D, Olivera BM, McIntosh JM: alpha-conotoxin AuIB selectively blocks alpha3 beta4 nicotinic acetylcholine receptors and nicotine-evoked norepinephrine release.

    J Neurosci 1998, 18(21):8571-8579. PubMed Abstract | Publisher Full Text OpenURL