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

Snake Venom: A Promising Source of Neurotoxins Targeting Voltage-Gated Potassium Channels

The venom derived from various sources of snakes represents a vast collection of predominantly protein-based toxins that exhibit a wide range of biological actions, including but not limited to inflammation, pain, cytotoxicity, cardiotoxicity, and neurotoxicity. The venom of a particular snake species is composed of several toxins, while the venoms of around 600 venomous snake species collectively encompass a substantial reservoir of pharmacologically intriguing compounds. Despite extensive research efforts, a significant portion of snake venoms remains uncharacterized. Recent findings have demonstrated the potential application of neurotoxins derived from snake venom in selectively targeting voltage-gated potassium channels (Kv). These neurotoxins include BPTI-Kunitz polypeptides, PLA2 neurotoxins, CRISPs, SVSPs, and various others. This study provides a comprehensive analysis of the existing literature on the significance of Kv channels in various tissues, highlighting their crucial role as proteins susceptible to modulation by diverse snake venoms. These toxins have demonstrated potential as valuable pharmacological resources and research tools for investigating the structural and functional characteristics of Kv channels.

Structure and selectivity engineering of the M1 muscarinic receptor toxin complex

Muscarinic toxins (MTs) are natural toxins produced by mamba snakes that primarily bind to muscarinic acetylcholine receptors (MAChRs) and modulate their function. Despite their similar primary and tertiary structures, MTs show distinct binding selectivity toward different MAChRs. The molecular details of how MTs distinguish MAChRs are not well understood. Here, we present the crystal structure of M₁AChR in complex with MT7, a subtype-selective anti-M₁AChR snake venom toxin. The structure reveals the molecular basis of the extreme subtype specificity of MT7 for M₁AChR and the mechanism by which it regulates receptor function. Through in vitro engineering of MT7 finger regions that was guided by the structure, we have converted the selectivity from M₁AChR toward M₂AChR, suggesting that the three-finger fold is a promising scaffold for developing G protein-coupled receptor modulators.

External release of entropy by synchronized movements of local secondary structures drives folding of a small, disulfide-bonded protein

A crucial mechanism to the formation of native, fully functional, 3D structures from local secondary structures is unraveled in this study. Through the introduction of various amino acid substitutions at four canonical β-turns in a three-fingered protein, Toxin α from Naja nigricollis, we found that the release of internal entropy to the external environment through the globally synchronized movements of local substructures plays a crucial role. Throughout the folding process, the folding species were saturated with internal entropy so that intermediates accumulated at the equilibrium state. Their relief from the equilibrium state was accomplished by the formation of a critical disulfide bridge, which could guide the synchronized movement of one of the peripheral secondary structure. This secondary structure collided with a core central structure, which flanked another peripheral secondary structure. This collision displaced the internal thermal fluctuations from the first peripheral structure to the second peripheral structure, where the displaced thermal fluctuations were ultimately released as entropy. Two protein folding processes that acted in succession were identified as the means to establish the flow of thermal fluctuations. The first process was the time-consuming assembly process, where stochastic combinations of colliding, native-like, secondary structures provided candidate structures for the folded protein. The second process was the activation process to establish the global mutual relationships of the native protein in the selected candidate. This activation process was initiated and propagated by a positive feedback process between efficient entropy release and well-packed local structures, which moved in synchronization. The molecular mechanism suggested by this experiment was assessed with a well-defined 3D structure of erabutoxin b because one of the turns that played a critical role in folding was shared with erabutoxin b.

A workflow for in silico design of hIL-10 and ebvIL-10 inhibitors using well-known miniprotein scaffolds

The over-expression of immune-suppressors such as IL-10 is a crucial landmark in both tumor progression, and latent viral and parasite infection. IL-10 is a multifunctional protein. Besides its immune-cell suppressive function, it also promotes B-cell tumorigenesis of lymphomas and melanoma. Human pathogens like unicellular parasites and viruses that remain latent inside B cells promote the over-expression of hIL-10 upon infection, which inhibits cell-mediated immune surveillance, and at the same time mediates B cell proliferation. The B-cell specific oncogenic latent virus Epstein-Barr virus (EBV) encodes a viral homologue of hIL-10 (ebvIL-10), expressed during lytic viral proliferation. Once expressed, ebvIL-10 inhibits cell-mediated immune surveillance, assuring EBV re-infection. During long-term latency, EBV-infected B cells over-express hIL-10 to assure B-cell proliferation, occasionally inducing EBV-mediated lymphomas. The amino acid sequences of hIL-10 and ebvIL-10 are more than 80% identical and thus have a very similar tridimensional structure. Based on their published crystallographic structures bound to their human receptor IL10R1, we report a structure-based design of hIL-10 and ebvIL-10 inhibitors based on 3 loops from IL10R1 that establish specific hydrogen bonds with the two IL10s. We have grafted these loops onto a permissible loop in three well-known miniprotein scaffolds-the Conus snail toxin MVIIA, the plant-derived trypsin inhibitor EETI, and the human appetite modulator AgRP. Our computational workflow described in detail below was invigorated by the negative and positive controls implemented, and therefore paves the way for future in vitro and in vivo validation assays of the IL-10 inhibitors engineered.

Soluble Immune Response Suppressor (SIRS): Reassessing the immunosuppressant potential of an elusive peptide

A previously studied immunosuppressive cytokine, Soluble Immune Response Suppressor (SIRS), may have relevance to current studies of immune suppression in a variety of human disease states. Despite extensive efforts using experimental models, mainly in mice, much remains to be discovered as to how autoimmune cells in mice and humans escape normal regulation and, conversely, how tumor cells evade evoking an immune response. It is the contention of this commentary that the literature pre-2000 contain results that might inform current studies. The broadly immunosuppressive protein, SIRS, was studied extensively from the 1970s to 1990s and culminated in the determination of the n-terminal 21mer sequence of this 15kDa protein which had high homology to the short neurotoxins from sea snakes, that are canonical members of the three finger neurotoxin superfamily (3FTx). It was not until 2007 that the prophylactic administration of the synthetic N-terminal peptide of the SIRS 21mer, identical to the published sequence, was reported to inhibit or delay the development of two autoimmune diseases in mice: experimental allergic encephalomyelitis (EAE) and type I diabetes (T1D). These findings were consistent with other studies of the 3FTx superfamily as important probes in the study of mammalian pharmacology. It is the perspective of this commentary that SIRS, SIRS peptide and the anti-peptide mAb, represent useful, pharmacologically-active probes for the study of the immune response as well as in the potential treatment of autoimmune, inflammatory diseases and cancer.

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.

Structural model of ligand-G protein-coupled receptor (GPCR) complex based on experimental double mutant cycle data: MT7 snake toxin bound to dimeric hM1 muscarinic receptor

The snake toxin MT7 is a potent and specific allosteric modulator of the human M1 muscarinic receptor (hM1). We previously characterized by mutagenesis experiments the functional determinants of the MT7-hM1 receptor interaction (Fruchart-Gaillard, C., Mourier, G., Marquer, C., Stura, E., Birdsall, N. J., and Servent, D. (2008) Mol. Pharmacol. 74, 1554-1563) and more recently collected evidence indicating that MT7 may bind to a dimeric form of hM1 (Marquer, C., Fruchart-Gaillard, C., Mourier, G., Grandjean, O., Girard, E., le Maire, M., Brown, S., and Servent, D. (2010) Biol. Cell 102, 409-420). To structurally characterize the MT7-hM1 complex, we adopted a strategy combining double mutant cycle experiments and molecular modeling calculations. First, thirty-three ligand-receptor proximities were identified from the analysis of sixty-one double mutant binding affinities. Several toxin residues that are more than 25 Å apart still contact the same residues on the receptor. As a consequence, attempts to satisfy all the restraints by docking the toxin onto a single receptor failed. The toxin was then positioned onto two receptors during five independent flexible docking simulations. The different possible ligand and receptor extracellular loop conformations were described by performing simulations in explicit solvent. All the docking calculations converged to the same conformation of the MT7-hM1 dimer complex, satisfying the experimental restraints and in which (i) the toxin interacts with the extracellular side of the receptor, (ii) the tips of MT7 loops II and III contact one hM1 protomer, whereas the tip of loop I binds to the other protomer, and (iii) the hM1 dimeric interface involves the transmembrane helices TM6 and TM7. These results structurally support the high affinity and selectivity of the MT7-hM1 interaction and highlight the atypical mode of interaction of this allosteric ligand on its G protein-coupled receptor target.

Conformational exchange is critical for the productivity of an oxidative folding intermediate with buried free cysteines

Much has been learned about the folding of proteins from comparative studies of the folding of proteins that are related in sequence and structure. Observation of the effects of mutations helps account for sequence-specific properties and large variations in folding rates observed in homologous proteins, which are not explained by structure-derived descriptions. The folding kinetics of variants of a β-stranded protein, toxin α from Naja nigricollis, depends on the length of their loop lk1. These proteins, named Tox60, Tox61, and Tox62, contain four disulfide bonds. We show that their oxidative refolding pathways are similar. Differences in these pathways are restricted to the last step of the reaction, that is, the closure of the last disulfide. At this step, two species of three-disulfide intermediates are observed: intermediate C lacking the B3 disulfide and intermediate D lacking the B2 disulfide. Surprisingly, D is the most productive intermediate for Tox61 despite the low accessibility of its free cysteines. However, in the case of Tox62, its conversion efficiency drops by 2 orders of magnitude and C becomes the most productive intermediate. NMR was used in order to study the structural dynamics of each of these intermediates. Both three-disulfide intermediates of Tox61 exist in two forms, exchanging on the 1- to 100-ms scale. One of these forms is structurally very close to the native Tox61, whereas the other is always significantly more flexible on a picosecond-to-nanosecond timescale. On the other hand, in the case of Tox62, the three-disulfide intermediates only show a native-like structure. The higher conformational heterogeneity of Tox61 intermediate D allows an increased accessibility of its free cysteines to oxidative agents, which explains its faster native disulfide formation. Thus, residue deletion in loop lk1 probably abrogates stabilizing intramolecular interactions, creates conformational heterogeneity, and increases the folding rate of Tox60 and Tox61 compared to Tox62.

Crystal structure of Natratoxin, a novel snake secreted phospholipaseA2 neurotoxin from Naja atra venom inhibiting A-type K+ currents

Snake secreted phospholipasesA2 (sPLA2s) are widely used as pharmacological tools to investigate their role in diverse pathophysiological processes. Some members of snake venom sPLA2s have been found to block voltage-activated K(+) channels (K(v) channels). However, most studies involved in their effects on ion channels were indirectly performed on motor nerve terminals while few studies were directly done on native neurons. Here, a novel snake sPLA2 peptide neurotoxin, Natratoxin, composed of 119 amino acid residues and purified from Naja atra venom was reported. It was characterized using whole-cell patch-clamp in acutely dissociated rat dorsal root ganglion (DRG) neurons. It was found to effectively inhibit A-type K(+) currents and cause alterations of channel gating characters, such as the shifts of steady-state activation and inactivation curves to hyperpolarization direction and changes of V(1/2) and slope factor. Therefore, Natratoxin was suggested to be a gating modifier of K(v) channel. In addition, this inhibitory effect was found to be independent of its enzymatic activity. These results suggested that the toxin enacted its inhibitory effect by binding to K(v) channel. To further elucidate the structural basis for this electrophysiological phenomenon, we determined the crystal structure of Natratoxin at 2.2 A resolution by molecular replacement method and refined to an R-factor of 0.190. The observed overall fold has a different structural organization from other K(+) channel inhibitors in animal toxins. Compared with other K(v) channel inhibitors, a similar putative functional surface in its C-terminal was revealed to contribute to protein-protein interaction in such a blocking effect. Our results demonstrated that the spatial distribution of key amino acid residues matters most in the recognition of this toxin towards its channel target rather than its type of fold.

Deduction of functional peptide motifs in scorpion toxins

Scorpion toxins are important physiological probes for characterizing ion channels. Molecular databases have limited functional annotation of scorpion toxins. Their function can be inferred by searching for conserved motifs in sequence signature databases that are derived statistically but are not necessarily biologically relevant. Mutation studies provide biological information on residues and positions important for structure-function relationship but are not normally used for extraction of binding motifs. 3D structure analyses also aid in the extraction of peptide motifs in which non-contiguous residues are clustered spatially. Here we present new, functionally relevant peptide motifs for ion channels, derived from the analyses of scorpion toxin native and mutant peptides.

Two-dimensional measurement of proton T1rho relaxation in unlabeled proteins: mobility changes in alpha-bungarotoxin upon binding of an acetylcholine receptor peptide

A method for the measurement of proton T(1)(rho) relaxation times in unlabeled proteins is described using a variable spin-lock pulse after the initial nonselective 90 degrees excitation in a HOHAHA pulse sequence. The experiment is applied to alpha-bungarotoxin (alpha-BTX) and its complex with a 25-residue peptide derived from the acetylcholine receptor (AChR) alpha-subunit. A good correlation between high T(1)(rho) values and increased local motion is revealed. In the free form, toxin residues associated with receptor binding according to the NMR structure of the alpha-BTX complex with an AChR peptide and the model for alpha-BTX with the AChR [Samson, A. O., et al. (2002) Neuron 35, 319-332] display high mobility. When the AChR peptide binds, a decrease in the relaxation times and the level of motion of residues involved in binding of the receptor alpha-subunit is exhibited, while residues implicated in binding gamma- and delta-subunits retain their mobility. In addition, the quantitative T(1)(rho) measurements enable us to corroborate the mapping of boundaries of the AChR determinant strongly interacting with the toxin [Samson, A. O., et al. (2001) Biochemistry 40, 5464-5473] and can similarly be applied to other protein complexes in which peptides represent one of the two interacting proteins. The presented method is advantageous because of its simplicity, generality, and time efficiency and paves the way for future investigation of proton relaxation rates in small unlabeled proteins.

Diversity of folds in animal toxins acting on ion channels

Animal toxins acting on ion channels of excitable cells are principally highly potent short peptides that are present in limited amounts in the venoms of various unrelated species, such as scorpions, snakes, sea anemones, spiders, insects, marine cone snails and worms. These toxins have been used extensively as invaluable biochemical and pharmacological tools to characterize and discriminate between the various ion channel types that differ in ionic selectivity, structure and/or cell function. Alongside the huge molecular and functional diversity of ion channels, a no less impressive structural diversity of animal toxins has been indicated by the discovery of an increasing number of polypeptide folds that are able to target these ion channels. Indeed, it appears that these peptide toxins have evolved over time on the basis of clearly distinct architectural motifs, in order to adapt to different ion channel modulating strategies (pore blockers compared with gating modifiers). Herein, we provide an up-to-date overview of the various types of fold from animal toxins that act on ion channels selective for K+, Na+, Ca2+ or Cl- ions, with special emphasis on disulphide bridge frameworks and structural motifs associated with these peptide folds.

Testing the nematophagous biological control strainPaecilomyces lilacinus251 for paecilotoxin production

Paecilomyces lilacinus is a nematophagous fungus currently developed as a biological control agent. In order to evaluate potential toxin production, culture extract and concentrated culture supernatant of P. lilacinus strain 251 were tested against Gram-negative and Gram-positive bacteria. High-performance liquid chromatography analysis was carried out to compare the chromatograms of P. lilacinus strain 251 with the chromatogram of known paecilotoxin. It was found that the 251 strain of P. lilacinus did not produce detectable levels of paecilotoxin or other toxins with antimicrobial activity.