|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 283, Issue 15, 9805-9813, April 11, 2008
Alanine Scanning Mutagenesis of the Prototypic Cyclotide Reveals a Cluster of Residues Essential for Bioactivity*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Cyclotides display a wide range of biological activities. In addition to their uterotonic activity (12), hemolytic (6), anti-HIV4 (13), anti-bacterial (14), anti-tumor (15), anti-fouling (16), and insecticidal (17) activities have been reported. The mechanism(s) of action of these activities are unknown but may involve interactions with membranes (18, 19), possibly via self-association of cyclotides to form multimeric species (20). Reported structures of cyclotides show the presence of a surface-exposed hydrophobic patch (2, 21) that may be involved in such interactions. Their natural function in plants appears to be as host defense agents against insect pests based on their potent inhibition of growth and development of caterpillar larvae of the cotton budworm and cotton bollworm (17, 22).
Fig. 1 shows the three-dimensional structure and amino acid sequence of kalata B1 and highlights the knotted arrangement of the three conserved disulfide bonds that are characteristic of all cyclotides. The cystine knot and cyclic backbone account for the exceptional resistance of the cyclotides to thermal, chemical, and enzymatic degradation compared with other peptides of similar size (23). This chemical and biological resilience has led to suggestions of the use of cyclotides as protein-engineering scaffolds for stabilizing bioactive peptide epitopes in the optimization of peptide drugs (24-27).
The sequences of currently known cyclotides reveal subsets of conserved and variable residues. In describing cyclotides it is convenient to discuss their structures in terms of the six backbone loops between successive cysteine residues. The cyclotides may be regarded as a natural combinatorial template in which substitutions occur in these loops, which protrude from the cystine knot core (Fig. 1). Loop 1 is the most highly conserved and contains a structurally important glutamate residue that forms a network of hydrogen bonds to the backbone of the first two residues of loop 3 (28). Cyclotide sequences have been divided into the Möbius and bracelet subclasses (3), which differ from one another mostly in loops 2, 3, and 5, but the characteristic discriminating feature is a cis-Pro peptide bond in loop 5 of the Möbius subfamily, so named for its conceptual 180° twist of the circular backbone.
|
A better appreciation of the structure-activity relationships of the cyclotides also has the potential to help uncover the mechanism(s) mediating their various bioactivities. Hence, in the current study we chemically synthesized a complete suite of kalata B1 mutants, in which all residues, except for the six conserved cysteines, were individually replaced with an alanine. We analyzed these alanine mutants for structural and functional perturbations relative to native kalata B1 and determined their activities against insects and human cells. This approach has provided significant insights into the structural constraints of the cyclotide framework and the functional elements for bioactivity and has provided a more detailed understanding of the possible self-association process of these remarkable peptides.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
NMR Spectroscopy—NMR solutions for assignment and structure determination were prepared in H2O:D2O (95:5) and adjusted to pH 5.0-5.5 (or 3.0-3.5 for mutant E7A) with 0.1 M NaOH and 0.1 M HCl. Samples for amide temperature coefficients were prepared in 100 mM sodium phosphate buffered at pH 5.5 with 5% D2O adjusted to pH 5.0-5.5 (or 3.0-3.5 for mutant E7A). Chemical shifts and amide temperature coefficients for mutant E7A were compared with the wild-type peptide under similar conditions. One- and two-dimensional spectra were acquired on ARX 500-MHz or 600-MHz spectrometers at 25 °C.
Structure Calculations—Preliminary structures were calculated using a torsion angle-simulated annealing protocol within the program DYANA (31). Final structures were calculated using simulated annealing and energy minimization protocols within CNS version 1.1 (32). 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 (33, 34). Structures consistent with restraints were subjected to further molecular dynamics and energy minimization in a water shell, as described by Linge and Nilges (35). Structures were analyzed using PROMOTIF (36) and PROCHECK-NMR (37). The Ramachandran statistics for E7A are 75.8% of residues in the most favored and 24.2% in the additionally allowed. For T27A the Ramachandran statistics are 87.4% in the most favored and 12.6% in the additionally allowed. Hemolytic Assays—Stock solutions of each peptide were prepared and standardized using UV spectroscopy (
= 6410 cm-1-1 at 215 nm) at 350 µM ± 3% then serially diluted in triplicate. Human erythrocytes were washed with repeated centrifugation in phosphate-buffered saline (pH 7.4) to remove serum contaminants. A 0.25% v/v stock solution was prepared from the cell pellet, and 100 µl of this stock was added to each diluted peptide solution aliquot and then incubated at 37 °C for 1 h. After centrifugation of intact cells the percentage of hemolysis of the supernatant was measured by visual absorption spectroscopy (
= 405 nm).
Insecticidal Assays—Adult Drosophila melanogaster were allowed to oviposit overnight on plates (7 g of agar/7 g of sucrose/260 ml of water/70 ml of apple juice), and the resulting eggs were collected and washed with repeated centrifugation. The volume of eggs for each trial could only be approximated to within ± 50 µl, so the wild-type peptide was included in each batch of assays on alanine mutants as a reference positive control. The egg pellet was diluted with 50% glycerol and aliquoted into vials containing 3 ml of basic larval media (10 g of agar/80 g of cornmeal/200 g of yeast/850 ml of water/100 ml of treacle/11 ml of 10% w/v nipogen/5 ml of propionic acid). Glycerol was added to increase the viscosity to assist in homogeneous dispersion of eggs throughout the suspension. After 24-h incubation test compounds were added as UV absorption standardized solutions to give final peptide concentrations of 230, 100, or 0 ppm (or 80, 35, and 0 µM, respectively). Larvae were then allowed to mature and pupate, and the mass of total flies emerging was measured over 2-3 days.
Stability Assays—Enzymatic digests and acid hydrolysis conditions were identical to those described previously (23). Briefly, 1 mg/ml solutions of WT-kB1, G6A, or E7A were prepared in 100 mM ammonium bicarbonate buffer (pH 8.1). The enzymes used in these studies were trypsin, chymotrypsin, and proteinase K. Enzymatic digestions (50:1 substrate-to-enzyme ratio) were carried out at 37 °C with aliquots taken at set times up to 24 h. DnaK fragment 2 was used as a control to ensure that the enzymes were active. With this control, and because we were primarily interested in quantifying short term degradation, if present, no corrections for autodigestion of the enzymes were required. Acid hydrolysis reactions (1 mg/ml peptide solution in 1 M HCl) were performed at two temperatures, 37 °C and 80 °C. Samples were analyzed by liquid chromatography/mass spectrometry on an Agilent Series 1100 HPLC coupled to a QSTAR mass spectrometer (Applied Biosystems, Foster City, CA) equipped with an electrospray ionization source. Data were collected, processed, and analyzed using the Analyst QS1.1 software package. Thermal stability was assessed by monitoring the quality of one-dimensional NMR spectra of compounds in water at pH 5.0-5.5 for G6A and 3.0-3.5 for E7A on cycling from 25 °C to 80 °C and back to 25 °C.
Metal Binding—Titration of 0.4 mM kalata B1, T27A, T16A, and E7A with MnCl2 was performed at 40 °C, pH 5.0. One-dimensional and total correlation spectroscopy spectra (80-ms mixing time) were acquired at concentrations of 0 and 2.0 mM MnCl2 on an ARX 500-MHz spectrometer. The paramagnetic relaxation enhancement induced by Mn2+ was qualitatively characterized by relative cross-peak intensity in the total correlation spectroscopy spectra,
![]() |
cross-peak intensities in the spectra of the sample with and without Mn2+. | RESULTS |
|---|
|
|
|---|
|
-proton chemical shifts of peptides are sensitive indicators of local structure (40) and were used to highlight regions of perturbed structure. The majority of alanine mutants had no significant
-H chemical shift perturbations compared with the wild-type peptide, as shown in Fig. 2A. A small subset of mutants showed some minor localized shift perturbations as shown in Fig. 2B. The alanine mutants in loop 2 account for the majority of these slightly perturbed structures. For mutant E7A several expected peaks were missing from the spectra at pH 5.0-5.5 (presumably due to broadening) but were present at pH 3.0-3.5, indicating that this mutant is less able to form a rigid secondary structure at the higher pH. Overall though, the chemical shift analysis showed that generally no major structural changes are induced by Ala substitution at any of the individual non-Cys residues and that the overall cyclotide fold is maintained in all mutants. Two alanine mutants were selected for full structure determination. Mutant E7A was selected to investigate the effect of removal of the internal hydrogen-bonding network from the glutamate side chain to the backbone amides of loop 3. Mutant T27A was selected as a relatively unperturbed example based on the chemical shift data. The three-dimensional structures of E7A and T27A were solved by simulated annealing using distance and angle restraints derived from two-dimensional NMR spectra. Simulated annealing calculations gave ensembles of 20 low energy structures with backbone root mean square deviations of 0.5 Å for E7A and 0.5 Å for T27A as summarized in Table 2. The structures of E7A and T27A are shown in Fig. 3, and they confirm that the overall fold of the mutants does not substantially differ from the wild-type peptide.
|
|
3 mg per test compound by using triplicates at two peptide concentrations. Triplicates were found to vary by <10%. Wild-type kalata B1 and alanine mutants gave reproducible dose response curves with an approximate LD50 of 150 ppm (52 µM) for the wild-type peptide. The mass of emerging adult flies varied inversely with the mass of eggs aliquoted. The activity of the alanine mutants relative to the wild-type peptide is reported in Table 1, where it is clear that substitution of a specific set of individual residues leads to a loss of activity. Thus, despite the fact that no Ala substitution causes a major structural change, there are functional effects of individual residues. The relative hemolysis of human erythrocytes caused by each mutant is also shown in Table 1, which indicates that residues in loops 1 and 3 along with a few other isolated residues strongly affect activity. Interestingly, a comparison of the results from both bioassays indicates that hemolytic activity and insecticidal activity depend on a similar subset of residues.
Stability Analysis—The stability of the G6A and E7A mutants was assessed by measuring their resistance to proteolysis and to acid hydrolysis. Mass spectra were examined for the presence of intact peptides and/or degradation products, and the percentage of intact peptide remaining was plotted against time. The two alanine mutants studied were completely impervious to enzymatic digestion by the three endoproteinases trypsin, chymotrypsin, and proteinase K, as was the prototypic cyclotide kalata B1. A linear control peptide degraded rapidly, with <5% remaining after 5 min under the same conditions. Acid hydrolysis experiments were conducted at 37 °C and 80 °C. The former represents likely conditions in the human stomach (pH
1, 37 °C), whereas the higher temperature was used to force the hydrolysis to occur and allow comparison of the stability of the mutants to the wild-type cyclotide. At 37 °C, kalata B1, and the two mutants were extremely stable, with no degradation products observed for up to 24 h. At 80 °C, the peptides degraded within 90 min and showed similar degradation curves from 0-90 min (data not shown).
Thermal stability was assessed by monitoring the one-dimensional NMR spectra of peptides on cycling from 25 °C to 80 °C and back to 25 °C. The spectra of both mutants were unchanged on return to the original temperature, but E7A lost more amide peaks at higher temperature than the wild-type peptide and G6A (data not shown). Overall, the data suggest that the structures of the mutants are highly resistant to extremes of temperature, pH, and enzyme conditions, but the E7A mutant is slightly more susceptible to thermal unfolding and is more flexible.
|
70% attenuation) of loop 1 residues. By contrast, E7A showed significantly less broadening (supplemental Fig. S1). It is clear that removal of the E7 carboxyl group reduces the effects of Mn2+, implicating it as the primary binding site. | DISCUSSION |
|---|
|
|
|---|
Chemical solid phase peptide synthesis was used to generate the 23 non-cysteine alanine mutants. Backbone cyclization was achieved using a modified form of native chemical ligation in which the linear peptide is synthesized with a C-terminal thioester that is repeatedly exchanged with cysteine thiol side chains, culminating in the N-terminal cysteine thiol that undergoes an irreversible exchange with the N-terminal amine to form the linking backbone amide (43). Exploitation of this mechanism permitted designing the synthesis with any linear precursor that terminated in an N-terminal cysteine. This allowed the synthesis of the alanine mutants to start at multiple different positions in the sequence, allowing the resin from a bulk synthesis of unmodified cyclotide sequence to be split at the appropriate point in the chain assembly process to produce a set of alanine mutants from the final loop(s).
The mutants were cyclized and subjected to oxidative folding in ammonium bicarbonate/2-propanol buffer, with the yields summarized in Table 1. All of the mutants formed a native fold in this buffer, with the exception of W23A and P24A, which involve substitutions in loop 5. There was a variation in yield of the correctly folded product relative to the wild-type peptide; peptides with lower yields had mutations primarily concentrated in loops 2 and 5, along with an isolated example of P3A in loop 6 having a low yield. The tight turns of loops 2 and 5 and their close packing in the structure of the native peptide explain the sensitivity of these regions to modification. Mutant G12A had a particularly low yield that is likely associated with its broadened chromatographic peaks during purification, suggesting that this mutant may be conformationally flexible. Introduction of a side chain at this position appears to disfavor the positive
-angle present for this highly conserved glycine that is present in almost half of known cyclotide sequences (44).
The chemical shifts of the
-protons of the alanine mutants relative to the wild-type peptide gave useful insights into the structural resilience of the cyclotide framework. The majority of mutants had no significant perturbations in shift, showing that the framework is highly tolerant to point substitutions. A small subset of mutants experienced minor local shift perturbations in the modified region (Fig. 2B) generally in mutants that also displayed lowered oxidative folding yields. The most significant changes in chemical shifts occurred in loop 2, with loops 5 and 6 displaying the lowest variation in shifts. Combined with the folding yield observations this finding indicates that future sequence modification attempts may be more successful if they concentrate on loops 3 or 6. Loop 3 represents the most diverse loop in terms of length in native cyclotide sequences and is extended into a short 310 helix in the bracelet cyclotides (28).
The glutamic acid in loop 1 is almost completely conserved in the cyclotides, only conservatively substituted with aspartic acid in one known sequence (45), and is involved in a network of hydrogen bonds that stabilize the structure (28). Analysis of the NMR structure of mutant E7A was initially hindered by the poor spectra and an atypical pH sensitivity of this compound. Numerous amide resonances were not detected at pH
5.5, indicating that this mutant undergoes a conformational exchange process not seen in the wild-type peptide. These resonances were more readily detected at lower pH, and the structure determined at pH 3.5 is shown in Fig. 3 (left structures). Despite the loss of the network of hydrogen bonds involving E7 the overall fold of E7A is maintained. Thus, the most significant structural effects of mutation of this conserved Glu are related to an increased flexibility of the structure at higher pH. This effect is similar to the loss of stability previously seen in acyclic permutants of kalata B1 (46, 47). The importance of the glutamic acid residue has also been recently highlighted by the observation that methylation of this side chain resulted in a 48-fold decrease in cytotoxicity activity for the bracelet cyclotide cycloviolacin O2 (48). This finding confirms the critical importance of this residue in both subfamilies of cyclotides. The structure of mutant T27A, determined as an example of a mutant that retained its bioactivity, showed no significant difference to the wild-type peptide structure, as shown in Fig. 3 (right structures). This finding supports the use of chemical shifts as a marker for a lack of structural change for the other mutants also.
|
We were at first very surprised to see that the surface-exposed hydrophobic patch of kalata B1 shown in Fig. 4c is clearly distinct from the residues that are crucial for bioactivity. Spin label studies recently demonstrated that this surface-exposed hydrophobic patch essentially forms the membrane binding face of the cyclotide kalata B1, as shown in Fig. 4d (42). This finding, together with the mutagenesis results from this work clearly suggests that more than simple membrane interactions must mediate bioactivity. The high level of sensitivity of bioactivity to changes in hydrophilic rather than hydrophobic residues indicates that the activity of kalata B1 is not dependent on gross physical characteristics and that the peptide does not act as a simple detergent. This extreme sensitivity of activity to loss of hydrophilic residues (e.g. T8A) or to small steric changes (e.g. G6A) and the clustering of critical residues instead indicate that a specific recognition process takes place that does not directly involve the membrane binding face.
|
The hydrophobic and putative self-recognition faces are approximately orthogonal to each other in Fig. 4. This is consistent with a self-association model in which four cyclotide molecules come together, with the activity-influencing hydrophilic surfaces facing inward and the hydrophobic surfaces oriented together on one face of the tetramer (Fig. 5). Such a reversible oligomerization could improve the absorption and dispersion of cyclotides in biological systems by allowing both the solubility of monomers in aqueous environments and of oligomers in the membrane, switching the effective hydrophobicity of the peptides.
The dramatic effect of the loss of structural rigidity on bioactivity seen in mutants E7A and G12A is also consistent with a self-association mechanism, where the conformation of the peptide that is capable of self-recognition is thrown into dynamic equilibrium with other conformations. In a tetramer dependent on multiple self-recognition interfaces any reduction in the average lifetime of each self-association face would have a multiplicative effect on the lifetime of the entire oligomer. Mutants that displayed a loss of structural rigidity as judged by NMR (E7A and G12A) displayed a dramatic loss of bioactivity. This is consistent with the reduced rigidity and total loss of activity in a series of acyclic permutants (46) and capped acyclic permutants (47).
Analysis of the stability of two inactive mutants was undertaken to determine if the resilience of the scaffold could be retained while eliminating activities that would be undesirable in pharmaceutical applications. Two mutants on loop 1 were selected, G6A because of the presumably steric nature of its interference with activity and E7A because of the possible structural contribution of this residue in the wild-type peptide. Mutant G6A was found to display similar resistance to proteolytic and thermal degradation and to extremes of pH as the wild-type peptide, consistent with a precise steric clash removing hemolytic activity in this mutant. Mutant E7A was also found to have comparable resistance to the wild-type peptide, indicating that the cyclic backbone and cystine knot alone are sufficient structural features to confer extreme stability. The resilience of both inactive mutants strengthens the potential of the cyclotide framework for stabilization of bioactive peptide epitopes for therapeutic applications.
Previously, it has been reported that the divalent cation Mn2+ binds to kalata B1 (41). The presence of a metal binding site may be functionally significant, because it may facilitate the interaction of kalata B1 (and other cyclotides) to membrane surfaces, particularly anionic membrane surfaces (42). A representative group of mutants, which were selected based on their activity and the presence/absence of a carboxyl group, were tested for Mn2+ binding and compared with kalata B1. Kalata B1 and the two carboxyl-containing mutants, T16A (an inactive mutant) and T27A (an active mutant), showed very similar responses to Mn2+. Significant broadening of the amide signals from loop 1 and complete broadening of the amide signals from loop 3 was observed, which is consistent with previous reports and identifies the metal binding face of kalata B1 (42). At similar concentrations of peptide and Mn2+, the amide signals around the corresponding region for E7A showed significantly less broadening, confirming that the carboxyl of Glu-7 is responsible for metal binding. This finding is consistent with the fact that carboxyl groups have been reported to coordinate Mn2+ binding in proteins (51); removal of the carboxyl in the E7A mutant quenches the binding. Because the inactive mutant T16A binds Mn2+, it appears that metal binding may not be an important factor for the activity of the cyclotides.
The current study has clearly shown that the cyclotide framework can accommodate sequence variations and maintain a well defined structure with the native fold. Given this conservation of the overall fold in the kalata B1 mutants, it appears likely that it is the side chains that have a greater role in maintaining the activity of the cyclotides. Such a scenario has also been observed for other disulfide-rich peptides such as the
-conotoxins where the native
-helical fold is maintained despite the sequence diversity throughout the family, but the selectivity and specificity for the nicotinic acetylcholine receptor varies enormously based on the side chains that are available to interact with the receptor subunits (52, 53).
In conclusion, significant advances have been made in the understanding of the adaptability of the cyclotide framework to modifications. Selected mutants lacking undesirable hemolytic activity were found to retain the characteristic stability of the framework, verifying its suitability for pharmacological applications. The conserved hydrogen-bonding network of the loop 1 glutamate to loop 3 was found to play a critical functional role, and the activity of kalata B1 was found to depend on a surface patch of predominantly hydrophilic residues. This information will be valuable in future experiments aimed at exploiting the cyclic cystine knot scaffold in drug design applications.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Fig. S1. ![]()
1 Supported by an Australian Postgraduate award. ![]()
2 An National Health and Medical Research Council Industry Fellow. ![]()
3 To whom correspondence should be addressed: Tel.: 61-7-3346-2019; Fax: 61-7-3346-2029; E-mail: d.craik{at}imb.uq.edu.au.
4 The abbreviations used are: HIV, human immunodeficiency virus; HPLC, high-performance liquid chromatography; kB1, kalata B1; E7A, [E7A]-kalata B1; T27A, [T27A]-kalata B1. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. K. Wang, S.-H. Hu, J. L. Martin, T. Sjogren, J. Hajdu, L. Bohlin, P. Claeson, U. Goransson, K. J. Rosengren, J. Tang, et al. Combined X-ray and NMR Analysis of the Stability of the Cyclotide Cystine Knot Fold That Underpins Its Insecticidal Activity and Potential Use as a Drug Scaffold J. Biol. Chem., April 17, 2009; 284(16): 10672 - 10683. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| All ASBMB Journals | Molecular and Cellular Proteomics |
| Journal of Lipid Research | ASBMB Today |