Motions and structural variability within toxins: implication for their use as scaffolds for protein engineering

Animal toxins are small proteins built on the basis of a few disulfide bonded frameworks. Because of their high variability in sequence and biologic function, these proteins are now used as templates for protein engineering. Here we report the extensive characterization of the structure and dynamics of two toxin folds, the "three-finger" fold and the short alpha/beta scorpion fold found in snake and scorpion venoms, respectively. These two folds have a very different architecture; the short alpha/beta scorpion fold is highly compact, whereas the "three-finger" fold is a beta structure presenting large flexible loops. First, the crystal structure of the snake toxin alpha was solved at 1.8-A resolution. Then, long molecular dynamics simulations (10 ns) in water boxes of the snake toxin alpha and the scorpion charybdotoxin were performed, starting either from the crystal or the solution structure. For both proteins, the crystal structure is stabilized by more hydrogen bonds than the solution structure, and the trajectory starting from the X-ray structure is more stable than the trajectory started from the NMR structure. The trajectories started from the X-ray structure are in agreement with the experimental NMR and X-ray data about the protein dynamics. Both proteins exhibit fast motions with an amplitude correlated to their secondary structure. In contrast, slower motions are essentially only observed in toxin alpha. The regions submitted to rare motions during the simulations are those that exhibit millisecond time-scale motions. Lastly, the structural variations within each fold family are described. The localization and the amplitude of these variations suggest that the regions presenting large-scale motions should be those tolerant to large insertions or deletions.     

Identification and characterization of a taxon-specific three-finger toxin from the venom of the Green Vinesnake (Oxybelis fulgidus; family Colubridae)

Snake venoms contain a variety of protein and peptide toxins, and the three-finger toxins (3FTxs) are among the best characterized family of venom proteins. The compact nature and highly conserved molecular fold of 3FTxs, together with their abundance in many venoms, has contributed to their utility in structure-function studies. Although many target the nicotinic acetylcholine receptor of vertebrate skeletal muscle, often binding with nanomolar K_ds, several non-conventional 3FTxs show pronounced taxon-specific neurotoxic effects. Here we describe the purification and characterization of fulgimotoxin, a monomeric 3FTx from the venom of Oxybelis fulgidus, a neotropical rear-fanged snake. Fulgimotoxin retains the canonical 5 disulfides of the non-conventional 3FTxs and is highly neurotoxic to lizards; however, mice are unaffected, demonstrating that this toxin is taxon-specific in its effects. Analysis of structural features of fulgimotoxin and other colubrid venom 3FTxs indicate the presence of a “colubrid toxin motif” (CYTLY) and a second conserved segment (WAVK) found in Boiga and Oxybelis taxon-specific 3FTxs, both in loop II. Because specific residues in loop II conventional α-neurotoxic 3FTxs are intimately associated with receptor binding, we hypothesize that this loop, with its highly conserved substitutions, confers taxon-specific neurotoxicity. These findings underscore the importance of rear-fanged snake venoms for understanding the evolution of toxin molecules and demonstrate that even among well-characterized toxin families, novel structural and functional motifs may be found.     

Conserved structural determinants in three-fingered protein domains

The three-dimensional structures of some components of snake venoms forming so-called 'three-fingered protein' domains (TFPDs) are similar to those of the ectodomains of activin, bone morphogenetic protein and transforming growth factor-beta receptors, and to a variety of proteins encoded by the Ly6 and Plaur genes. The analysis of sequences of diverse snake toxins, various ectodomains of the receptors that bind activin and other cytokines, and numerous gene products encoded by the Ly6 and Plaur families of genes has revealed that they differ considerably from each other. The sequences of TFPDs may consist of up to six disulfide bonds, three of which have the same highly conserved topology. These three disulfide bridges and an asparagine residue in the C-terminal part of TFPDs are essential for the TFPD-like fold. Analyses of the three-dimensional structures of diverse TFPDs have revealed that the three highly conserved disulfides impose a major stabilizing contribution to the TFPD-like fold, in both TFPDs contained in some snake venoms and ectodomains of several cellular receptors, whereas the three remaining disulfide bonds impose specific geometrical constraints in the three fingers of some TFPDs.     

Solution structure of long neurotoxin NTX-1 from the venom of Naja naja oxiana by 2D-NMR spectroscopy.

The NMR solution structures of NTX-1 (PDB code 1W6B and BMRB 6288), a long neurotoxin isolated from the venom of Naja naja oxiana, and the molecular dynamics simulation of these structures are reported. Calculations are based on 1114 NOEs, 19 hydrogen bonds, 19 dihedral angle restraints and secondary chemical shifts derived from 1H to 13C HSQC spectrum. Similar to other long neurotoxins, the three-finger like structure shows a double and a triple stranded beta-sheet as well as some flexible regions, particularly at the tip of loop II and the C-terminal tail. The solution NMR and molecular dynamics simulated structures are in good agreement with root mean square deviation values of 0.23 and 1 A for residues involved in beta-sheet regions, respectively. The overall fold in the NMR structure is similar to that of the X-ray crystallography, although some differences exist in loop I and the tip of loop II. The most functionally important residues are located at the tip of loop II and it appears that the mobility and the local structure in this region modulate the binding of NTX-1 and other long neurotoxins to the nicotinic acetylcholine receptor.     

Chemical engineering of a three-fingered toxin with anti-alpha7 neuronal acetylcholine receptor activity

Though it possesses four disulfide bonds the three-fingered fold is amenable to chemical synthesis, using a Fmoc-based method. Thus, we synthesized a three-fingered curaremimetic toxin from snake with high yield and showed that the synthetic and native toxins have the same structural and biological properties. Both were characterized by the same 2D NMR spectra, identical high binding affinity (K(d) = 22 +/- 5 pM) for the muscular acetylcholine receptor (AChR) and identical low affinity (K(d) = 2.0 +/- 0.4 microM) for alpha7 neuronal AchR. Then, we engineered an additional loop cyclized by a fifth disulfide bond at the tip of the central finger. This loop is normally present in longer snake toxins that bind with high affinity (K(d) = 1-5 nM) to alpha7 neuronal AchR. Not only did the chimera toxin still bind with the same high affinity to the muscular AchR but also it displayed a 20-fold higher affinity (K(d) = 100 nM) for the neuronal alpha7 AchR, as compared with the parental short-chain toxin. This result demonstrates that the engineered loop contributes, at least in part, to the high affinity of long-chain toxins for alpha7 neuronal receptors. That three-fingered proteins with four or five disulfide bonds are amenable to chemical synthesis opens new perspectives for engineering new activities on this fold.     

Solution NMR Structures of Pyrenophora tritici-repentis ToxB and Its Inactive Homolog Reveal Potential Determinants of Toxin Activity

Pyrenophora tritici-repentis Ptr ToxB (ToxB) is a proteinaceous host-selective toxin produced by Pyrenophora tritici-repentis (P. tritici-repentis), a plant pathogenic fungus that causes the disease tan spot of wheat. One feature that distinguishes ToxB from other host-selective toxins is that it has naturally occurring homologs in non-pathogenic P. tritici-repentis isolates that lack toxic activity. There are no high-resolution structures for any of the ToxB homologs, or for any protein with >30% sequence identity, and therefore what underlies activity remains an open question. Here, we present the NMR structures of ToxB and its inactive homolog Ptr toxb. Both proteins adopt a β-sandwich fold comprising three strands in each half that are bridged together by two disulfide bonds. The inactive toxb, however, shows higher flexibility localized to the sequence-divergent β-sandwich half. The absence of toxic activity is attributed to a more open structure in the vicinity of one disulfide bond, higher flexibility, and residue differences in an exposed loop that likely impacts interaction with putative targets. We propose that activity is regulated by perturbations in a putative active site loop and changes in dynamics distant from the site of activity. Interestingly, the new structures identify AvrPiz-t, a secreted avirulence protein produced by the rice blast fungus, as a structural homolog to ToxB. This homology suggests that fungal proteins involved in either disease susceptibility such as ToxB or resistance such as AvrPiz-t may have a common evolutionary origin.     

Selective Removal of Individual Disulfide Bonds within a Potato Type II Serine Proteinase Inhibitor from Nicotiana alata Reveals Differential Stabilization of the Reactive-Site Loop

The 53-amino-acid trypsin inhibitor 1 from Nicotiana alata (T1) belongs to the potato type II family also known as the PinII family of proteinase inhibitors, one of the major families of canonical proteinase inhibitors. T1 contains four disulfide bonds, two of which (C4-C41 and C8-C37) stabilize the reactive-site loop. To investigate the influence of these two disulfide bonds on the structure and function of potato II inhibitors, we constructed two variants of T1, C4A/C41A-T1 and C8A/C37A-T1, in which these two disulfide bonds were individually removed and replaced by alanine residues. Trypsin inhibition assays show that wild-type T1 has a K_i of < 5 nM, C4A/C41A-T1 has a weaker K_i of ∼ 350 nM, and the potency of the C8A/C37A variant is further decreased to a K_i of ∼ 1.8 μM. To assess the influence of the disulfide bonds on the structure of T1, we determined the structure and dynamics of both disulfide variants by NMR spectroscopy. The structure of C4A/C41A-T1 and the amplitude of intrinsic flexibility in the reactive-site loop resemble that of the wild-type protein closely, despite the lack of the C4-C41 disulfide bond, whereas the timescale of motions is markedly decreased. The rescue of the structure despite loss of a disulfide bond is due to a previously unrecognized network of interactions, which stabilizes the structure of the reactive-site loop in the region of the missing disulfide bond, while allowing intrinsic motions on a fast (picosecond-nanosecond) timescale. In contrast, no comparable interactions are present around the C8-C37 disulfide bond. Consequently, the reactive-site loop becomes disordered and highly flexible in the structure of C8A/C37A-T1, making it unable to bind to trypsin. Thus, the reactive-site loop of T1 is stabilized differently by the C8-C37 and C4-C41 disulfide bonds. The C8-C37 disulfide bond is essential for the inhibitory activity of T1, whereas the C4-C41 disulfide bond is not as critical for maintaining the three-dimensional structure and function of the molecule but is responsible for maintaining flexibility of the reactive-site loop on a microsecond-nanosecond timescale.     

The role of conserved Glu residue on cyclotide stability and activity: a structural and functional study of kalata B12, a naturally occurring Glu to Asp mutant

Cyclotides are a family of plant defense proteins with a unique cyclic backbone and cystine knot. Their remarkable stability under harsh thermal, enzymatic, and chemical conditions, combined with their range of bioactivities, including anti-HIV activity, underpins their potential as protein drug scaffolds. The vast majority of cyclotides possess a conserved glutamate residue in loop 1 of the sequence that is involved in a structurally important network of hydrogen bonds to an adjacent loop (loop 3). A single native cyclotide sequence, kalata B12, has been discovered that has an aspartic acid in this otherwise conserved position. Previous studies have determined that methylation of the glutamate or substitution with alanine abolishes the membrane disrupting activity that is characteristic of the family. To further understand the role of this conserved structural feature, we studied the folding, structure, stability, and activity of the natural aspartic acid variant kalata B12 and compared it to the prototypical cyclotide kalata B1, along with its glutamate to alanine or aspartate mutants. We show that the overall fold of kalata B12 is similar to the structure of other cyclotides, confirming that the cyclotide framework is robust and tolerant to substitution, although the structure appears to be more flexible than other cyclotides. Modification of the glutamate in kalata B1 or replacing the aspartate in kalata B12 with a glutamate reduces the efficiency of oxidative folding relative to the native peptides. The bioactivity of all modified glutamate cyclotides is abolished, suggesting an important functional role of this conserved residue. Overall, this study shows that the presence of a glutamic acid in loop 1 of the cyclotides improves stability and is essential for the membrane disrupting activity of cyclotides.     

Engineering of three-finger fold toxins creates ligands with original pharmacological profiles for muscarinic and adrenergic receptors

Protein engineering approaches are often a combination of rational design and directed evolution using display technologies. Here, we test "loop grafting," a rational design method, on three-finger fold proteins. These small reticulated proteins have exceptional affinity and specificity for their diverse molecular targets, display protease-resistance, and are highly stable and poorly immunogenic. The wealth of structural knowledge makes them good candidates for protein engineering of new functionality. Our goal is to enhance the efficacy of these mini-proteins by modifying their pharmacological properties in order to extend their use in imaging, diagnostics and therapeutic applications. Using the interaction of three-finger fold toxins with muscarinic and adrenergic receptors as a model, chimeric toxins have been engineered by substituting loops on toxin MT7 by those from toxin MT1. The pharmacological impact of these grafts was examined using binding experiments on muscarinic receptors M1 and M4 and on the α(1A)-adrenoceptor. Some of the designed chimeric proteins have impressive gain of function on certain receptor subtypes achieving an original selectivity profile with high affinity for muscarinic receptor M1 and α(1A)-adrenoceptor. Structure-function analysis supported by crystallographic data for MT1 and two chimeras permits a molecular based interpretation of these gains and details the merits of this protein engineering technique. The results obtained shed light on how loop permutation can be used to design new three-finger proteins with original pharmacological profiles.     

Evidence for the existence of a common ancestor of scorpion toxins affecting ion channels

All scorpion toxins from different 30 species are simply reviewed. A new classification system of scorpion toxins is first proposed: scorpion toxins are classified into three families (long-chain scorpion toxins with 4 disulfide bridges, short-chain scorpion toxins with 3 disulfide bridges, and intermediate-type scorpion toxins with 3 or 4 disulfide bridges). Intermediate-type scorpion toxins provide a strong proof for the conclusion that channel toxins from scorpion venoms evolve from a common ancestor. Common organization of precursor nucleotides and genomic sequence, similar 3-dimensional structure, and the existence of intermediate type scorpion toxins and functionally intercrossing scorpion toxins show that all scorpion toxins affecting ion channels evolve from the common ancestor, which produce millions of scorpion toxins with function-diversity.     

Disulfide Bonds Play a Critical Role in the Structure and Function of the Receptor-binding Domain of the SARS-CoV-2 Spike Antigen

The current coronavirus pandemic is exerting a tremendously detrimental impact on global health. The Spike proteins of coronaviruses, responsible for cell receptor binding and viral internalization, possess multiple and frequently conserved disulfide bonds raising the question about their role in these proteins. Here, we present a detailed structural and functional investigation of the disulfide bonds of the SARS-CoV-2 Spike receptor-binding domain (RBD). Molecular dynamics simulations of the RBD predict increased flexibility of the surface loops when the four disulfide bonds of the domain are reduced. This flexibility is particularly prominent for the disulfide bond-containing surface loop (residues 456–490) that participates in the formation of the interaction surface with the Spike cell receptor ACE2. In vitro, disulfide bond reducing agents affect the RBD secondary structure, lower its melting temperature from 52 °C to 36–39 °C and decrease its binding affinity to ACE2 by two orders of magnitude at 37 °C. Consistent with these in vitro findings, the reducing agents tris(2-carboxyethyl)phosphine (TCEP) and dithiothreitol (DTT) were able to inhibit viral replication at low millimolar levels in cell-based assays. Our research demonstrates the mechanism by which the disulfide bonds contribute to the molecular structure of the RBD of the Spike protein, allowing the RBD to execute its viral function.     

Location of the disulfide bonds in human plasma prekallikrein: the presence of four novel apple domains in the amino-terminal portion of the molecule

The location of 16 of the 18 disulfide bonds in human plasma prekallikrein was determined by amino acid sequence analysis of cystinyl peptides produced by chemical and enzymatic digestions. A unique structure, named the apple domain, was established for each of the four tandem repeats in the amino-terminal portion of the molecule. The apple domains (90 or 91 amino acids) contain 3 highly conserved disulfide bonds linking the first and sixth, second and fifth, and third and fourth half-cystine residues present in each repeat. The fourth tandem repeat contains an extra disulfide bond that forms a second small loop within the apple domain. The carboxyl-terminal portion of plasma prekallikrein containing the catalytic region of the molecule was found to have disulfide bonds located in positions similar to those of other serine proteases.     

Structure of variant 2 scorpion toxin from Centruroides sculpturatus Ewing

Centruroides sculpturatus Ewing variant 2 toxin (CsE-v2) is a neurotoxin isolated from the venom of a scorpion native to the Arizona desert. The structure of CsE-v2 was solved in two different crystal forms using a combination of molecular replacement and multiple isomorphous replacement techniques. Crystals of CsE-v2 display a temperature-dependent, reversible-phase transition near room temperature. At lower temperature the space group changes from P3(2)21 to P3(1)21 with an approximate doubling of the C-axis. The small-cell structure, which has one molecule per asymmetric unit, has an R factor of 0.229 at 2.8 A resolution. The large-cell structure has two molecules per asymmetric unit and was refined at 2.2 A resolution to an R factor of 0.255. CsE-v2 is a rigid, compact structure with four intrachain disulfide bonds. The structure is similar to other long-chain beta neurotoxins, and the largest differences occur in the last six residues. The high-resolution structure of CsE-v2 corrects an error in the reported C-terminal sequence; the terminal tripeptide sequence is Ser 64-Cys 65-Ser 66 rather than Ser 64-Ser 65-Cys 66. Comparison of CsE-v2 with long-chain alpha toxins reveals four insertions and one deletion, as well as additional residues at the N and C termini. Structural alignment of alpha and beta toxins suggests that the primary distinguishing feature between the two classes is the length of the loop between the second and third strands in a three-strand beta sheet. The shorter loop in alpha toxins exposes a critical lysine side chain, whereas the longer loop in beta toxins buries the corresponding basic residue (either arginine or lysine).     

Role of individual cysteine residues and disulfide bonds in the structure and function of Aspergillus ribonucleolytic toxin restrictocin.

Restrictocin, produced by the fungus Aspergillus restrictus, belongs to the group of ribonucleolytic toxins called ribotoxins. It specifically cleaves a single phosphodiester bond in a conserved stem and loop structure in the 28S rRNA of large ribosomal subunit and potently inhibits eukaryotic protein synthesis. Restrictocin contains 149 amino acid residues and includes four cysteines at positions 5, 75, 131, and 147. These cysteine residues are involved in the formation of two disulfide bonds, one between Cys 5 and Cys 147 and another between Cys 75 and Cys 131. In the current study, all four cysteine residues were changed to alanine individually and in different combinations by site-directed mutagenesis so as to remove one or both the disulfides. The mutants were expressed and purified from Escherichia coli. Removal of any cysteine or any one of the disulfide bonds individually did not affect the ability of the toxin to specifically cleave the 28S rRNA or to inhibit protein synthesis in vitro. However, the toxin without both disulfide bonds completely lost both ribonucleolytic and protein synthesis inhibition activities. The active mutants, containing only one disulfide bond, exhibited relatively high susceptibility to trypsin digestion. Thus, none of the four cysteine residues is directly involved in restrictocin catalysis; however, the presence of any one of the two disulfide bonds is absolutely essential and sufficient to maintain the enzymatically active conformation of restrictocin. For maintenance of the unique stability displayed by the native toxin, both disulfide bonds are required.     

Structure of the Sensor Domain of Mycobacterium tuberculosis PknH Receptor Kinase Reveals a Conserved Binding Cleft

Since their discovery over 20 years ago, eukaryotic-like transmembrane receptor Ser/Thr protein kinases (STPKs) have been shown to play critical roles in the virulence, growth, persistence, and reactivation of many bacteria. Information regarding the signals transmitted by these proteins, however, remains scarce. To enhance understanding of the basis for STPK receptor signaling, we determined the 1.7-Å-resolution crystal structure of the extracellular sensor domain of the Mycobacterium tuberculosis receptor STPK, PknH (Rv1266c). The PknH sensor domain adopts an unanticipated fold containing two intramolecular disulfide bonds and a large hydrophobic and polar cleft. The residues lining the cleft and those surrounding the disulfide bonds are conserved. These results suggest that PknH binds a small-molecule ligand that signals by changing the location or quaternary structure of the kinase domain.     

Synthetic peptide antigens derived from long-chain alpha-neurotoxins: Immunogenicity effect against elapid venoms

Three-finger toxins (3FTXs), especially α-neurotoxins, are the most poorly neutralized elapid snake toxins by current antivenoms. In this work, the conserved structural similarity and motif arrangements of long-chain α-neurotoxins led us to design peptides with consensus sequences. Eight long-chain α-neurotoxins (also known as Type II) were used to generate a consensus sequence from which two peptides were chemically synthesized, LCP1 and LCP2. Rabbit sera raised against them were able to generate partially-neutralizing antibodies, which delayed mice mortality in neutralization assays against Naja haje, Dendrospis polylepis and Ophiophagus hannah venoms.     

Folding studies of Cox17 reveal an important interplay of cysteine oxidation and copper binding

Cox17 is a key mitochondrial copper chaperone involved in the assembly of cytochrome c oxidase (COX). The NMR solution structure of the oxidized apoCox17 isoform consists of a coiled-coil conformation stabilized by two disulfide bonds involving Cys(26)/Cys(57) and Cys(36)/Cys(47). This appears to be a conserved tertiary fold of a class of proteins, localized within the mitochondrial intermembrane space, that contain a twin Cys-x(9)-Cys sequence motif. An isomerization of one disulfide bond from Cys(26)/Cys(57) to Cys(24)/Cys(57) is required prior to Cu(I) binding to form the Cu(1)Cox17 complex. Upon further oxidation of the apo-protein, a form with three disulfide bonds is obtained. The reduction of all disulfide bonds provides a molten globule form that can convert to an additional conformer capable of binding up to four Cu(I) ions in a polycopper cluster. This form of the protein is oligomeric. These properties are framed within a complete model of mitochondrial import and COX assembly.     

Plethodontid modulating factor, a hypervariable salamander courtship pheromone in the three-finger protein superfamily

The soluble members of the three-finger protein superfamily all share a relatively simple 'three-finger' structure, yet perform radically different functions. Plethodontid modulating factor (PMF), a pheromone protein produced by the lungless salamander, Plethodon shermani, is a new and unusual member of this group. It affects female receptivity when delivered to the female's nares during courtship. As with other plethodontid pheromone genes, PMF is hyperexpressed in a specialized male mental (chin) gland. Unlike other plethodontid pheromone genes, however, PMF is also expressed at low levels in the skin, liver, intestine and kidneys of both sexes. The PMF sequences obtained from all tissue types were highly variable, with 103 unique haplotypes identified which averaged 35% sequence dissimilarity (range 1-60%) at the protein level. Despite this variation, however, all PMF sequences contained a conserved approximately 20-amino-acid secretion signal sequence and a pattern of eight cysteines that is also found in cytotoxins and short neurotoxins from snake venoms, as well as xenoxins from Xenopus. Although they share a common cysteine pattern, PMF isoforms differ from other three-finger proteins in: (a) amino-acid composition outside of the conserved motif; (b) length of the three distinguishing 'fingers'; (c) net charge at neutral pH. Whereas most three-finger proteins have a net positive charge at pH 7.0, PMF has a high net negative charge at neutral pH (pI range of most PMFs 3.5-4.0). Sequence comparisons suggest that PMF belongs to a distinct multigene subfamily within the three-finger protein superfamily.