Snails In Silico: A Review of Computational Studies on the Conopeptides

Force Field Analysis Software and Tools (FFAST): Assessing Machine Learning Force Fields under the Microscope

As the sophistication of machine learning force fields (MLFF) increases to match the complexity of extended molecules and materials, so does the need for tools to properly analyze and assess the practical performance of MLFFs. To go beyond average error metrics and into a complete picture of a model's applicability and limitations, we developed FFAST (force field analysis software and tools): a cross-platform software package designed to gain detailed insights into a model's performance and limitations, complete with an easy-to-use graphical user interface. The software allows the user to gauge the performance of any molecular force field,─such as popular state-of-the-art MLFF models, ─ on various popular data set types, providing general prediction error overviews, outlier detection mechanisms, atom-projected errors, and more. It has a 3D visualizer to find and picture problematic configurations, atoms, or clusters in a large data set. In this paper, the example of the MACE and NequIP models is used on two data sets of interest [stachyose and docosahexaenoic acid (DHA)]─to illustrate the use cases of the software. With this, it was found that carbons and oxygens involved in or near glycosidic bonds inside the stachyose molecule present increased prediction errors. In addition, prediction errors on DHA rise as the molecule folds, especially for the carboxylic group at the edge of the molecule. We emphasize the need for a systematic assessment of MLFF models for ensuring their successful application to the study of dynamics of molecules and materials.

Inferring Pathways of Oxidative Folding from Prefolding Free Energy Landscapes of Disulfide-Rich Toxins

Short, cysteine-rich peptides can exist in stable or metastable structural ensembles due to the number of possible patterns of formation of their disulfide bonds. One interesting subset of this peptide group is the conotoxins, which are produced by aquatic snails in the family Conidae. The μ conotoxins, which are antagonists and blockers of the voltage-gated sodium channel, exist in a folding spectrum: on one end of the spectrum are more hirudin-like folders, which form disulfide bonds and then reshuffle them, leading to an ensemble of kinetically trapped isomers, and on the other end are more BPTI-like folders, which form the native disulfide bonds one by one in a particular order, leading to a preponderance of conformations existing in a single stable state. In this Article, we employ the composite diffusion map approach to study the unified free energy surface of prefolding μ-conotoxin equilibrium. We identify the two most important nonlinear collective modes of the unified folding landscape and demonstrate that in the absence of their disulfides, the conotoxins can be thought of as largely disordered polymers. A small increase in the number of hydrophobic residues in the protein shifts the free energy landscape toward hydrophobically collapsed coil conformations responsible for cysteine proximity in hirudin-like folders, compared to semiextended coil conformations with more distal cysteines in BPTI-like folders. Overall, this work sheds important light on the folding processes and free energy landscapes of cysteine-rich peptides and demonstrates the extent to which sequence and length contribute to these landscapes.

Bibliometric Review of the Literature on Cone Snail Peptide Toxins from 2000 to 2022

The venom of marine cone snails is mainly composed of peptide toxins called conopeptides, among which conotoxins represent those that are disulfide-rich. Publications on conopeptides frequently state that conopeptides attract considerable interest for their potent and selective activity, but there has been no analysis yet that formally quantifies the popularity of the field. We fill this gap here by providing a bibliometric analysis of the literature on cone snail toxins from 2000 to 2022. Our analysis of 3028 research articles and 393 reviews revealed that research in the conopeptide field is indeed prolific, with an average of 130 research articles per year. The data show that the research is typically carried out collaboratively and worldwide, and that discoveries are truly a community-based effort. An analysis of the keywords provided with each article revealed research trends, their evolution over the studied period, and important milestones. The most employed keywords are related to pharmacology and medicinal chemistry. In 2004, the trend in keywords changed, with the pivotal event of that year being the approval by the FDA of the first peptide toxin drug, ziconotide, a conopeptide, for the treatment of intractable pain. The corresponding research article is among the top ten most cited articles in the conopeptide literature. From the time of that article, medicinal chemistry aiming at engineering conopeptides to treat neuropathic pain ramped up, as seen by an increased focus on topological modifications (e.g., cyclization), electrophysiology, and structural biology.

Cone snail species off the Brazilian coast and their venoms: a review and update

The genus Conus includes over 900 species of marine invertebrates known as cone snails, whose venoms are among the most powerful described so far. This potency is mainly due to the concerted action of hundreds of small bioactive peptides named conopeptides, which target different ion channels and membrane receptors and thus interfere with crucial physiological processes. By swiftly harpooning and injecting their prey and predators with such deadly cocktails, the slow-moving cone snails guarantee their survival in the harsh, competitive marine environment. Each cone snail species produces a unique venom, as the mature sequences of conopeptides from the venoms of different species share very little identity. This biochemical diversity, added to the numerous species and conopeptides contained in their venoms, results in an immense biotechnological and therapeutic potential, still largely unexplored. That is especially true regarding the bioprospection of the venoms of cone snail species found off the Brazilian coast - a region widely known for its biodiversity. Of the 31 species described in this region so far, only four - Conus cancellatus, Conus regius, Conus villepinii, and Conus ermineus - have had their venoms partially characterized, and, although many bioactive molecules have been identified, only a few have been actually isolated and studied. In addition to providing an overview on all the cone snail species found off the Brazilian coast to date, this review compiles the information on the structural and pharmacological features of conopeptides and other molecules identified in the venoms of the four aforementioned species, paving the way for future studies.

High-Throughput Prediction and Design of Novel Conopeptides for Biomedical Research and Development

Cone snail venoms have been considered a valuable treasure for international scientists and businessmen, mainly due to their pharmacological applications in development of marine drugs for treatment of various human diseases. To date, around 800 Conus species are recorded, and each of them produces over 1,000 venom peptides (termed as conopeptides or conotoxins). This reflects the high diversity and complexity of cone snails, although most of their venoms are still uncharacterized. Advanced multiomics (such as genomics, transcriptomics, and proteomics) approaches have been recently developed to mine diverse Conus venom samples, with the main aim to predict and identify potentially interesting conopeptides in an efficient way. Some bioinformatics techniques have been applied to predict and design novel conopeptide sequences, related targets, and their binding modes. This review provides an overview of current knowledge on the high diversity of conopeptides and multiomics advances in high-throughput prediction of novel conopeptide sequences, as well as molecular modeling and design of potential drugs based on the predicted or validated interactions between these toxins and their molecular targets.

Absolute Structure Determination and Kv1.5 Ion Channel Inhibition Activities of New Debromoaplysiatoxin Analogues

Potassium channel Kv1.5 has been considered a key target for new treatments of atrial tachyarrhythmias, with few side effects. Four new debromoaplysiatoxin analogues with a 6/6/12 fused ring system were isolated from marine cyanobacterium Lyngbya sp. Their planar structures were elucidated by HRESIMS, 1D and 2D NMR. The absolute configuration of oscillatoxin J (1) was determined by single-crystal X-ray diffraction, and the absolute configurations of oscillatoxin K (2), oscillatoxin L (3) and oscillatoxin M (4) were confirmed on the basis of GIAO NMR shift calculation followed by DP4 analysis. The current study confirmed the absolute configuration of the pivotal chiral positions (7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R) at traditional ATXs with 6/12/6 tricyclic ring system. Compound 1, 2 and 4 exhibited blocking activities against Kv1.5 with IC50 values of 2.61 ± 0.91 µM, 3.86 ± 1.03 µM and 3.79 ± 1.01 µM, respectively. However, compound 3 exhibited a minimum effect on Kv1.5 at 10 µM. Furthermore, all of these new debromoaplysiatoxin analogs displayed no apparent activity in a brine shrimp toxicity assay.

Structural and dynamical effects of targeted mutations on μO-Conotoxin MfVIA: Molecular simulation studies

Conotoxins are a group of cysteine-rich, neurotoxic peptides isolated from the venom of marine cone snails. MfVIA is a member of the μO-conotoxin family, and acts as an inhibitor of subtype 1.8 voltage-gated sodium ion channels (Na_V1.8). The unique selectivity of MfVIA as an inhibitor of Na_V1.8 makes it an ideal peptide for elucidation of the physiological functions of this voltage-gated ion channel. Previous experimental studies of point mutations of MfVIA showed that the double mutant [E5K,E8K] exhibited greater activity at Na_V1.8 relative to the wild-type toxin. The present study employs molecular dynamics (MD) simulations to examine the effects of various mutations at these key residues (E5 and E8) on the structure and dynamics of MfVIA. Five double mutants were studied, in which the positions 5 and 8 residues were mutated to amino acids with a range of different physicochemical properties, namely [E5A,E8A], [E5D,E8D], [E5F,E8F], [E5K,E8K], and [E5R,E8R]. Except for [E5D,E8D], all of the mutants tend to show decreased contacts at the N-terminus owing to the loss of the R1-E5 salt bridge relative to that of the wild-type, which subsequently lead to greater exposure and flexibility of the N-terminus for most of the mutant peptides studied, potentially rendering them more able to interact with other species, including Na_V1.8. Molecular docking studies of the peptides to Na_V1.8 via different binding mechanisms suggest that the [E5R, E8R] mutant may be especially worthy of further investigation owing to its predicted binding mode, which differs markedly from those of the other peptides in this study.

Neurobiological activity of conotoxins via sodium channel modulation

Conotoxins (CnTX) are bioactive peptides produced by marine molluscs belonging to Conus genus. The biochemical structure of these venomous peptides is characterized by a low number of amino acids linked with disulfide bonds formed by a high degree of post-translational modifications and glycosylation steps which increase the diversity and rate of evolution of these molecules. CnTX different isoforms are known to target ion channels and, in particular, voltage-gated sodium (Na⁺) channels (Na_v channels). These are transmembrane proteins fundamental in excitable cells for generating the depolarization of plasma membrane potential known as action potential which propagates electrical signals in muscles and nerves for physiological functions. Disorders in Na_v channel activity have been shown to induce neurological pathologies and pain states. Here, we describe the current knowledge of CnTX isoform modulation of the Na_v channel activity, the mechanism of action and the potential therapeutic use of these toxins in counteracting neurological dysfunctions.

Graph-Directed Approach for Downselecting Toxins for Experimental Structure Determination

Conotoxins are short, cysteine-rich peptides of great interest as novel therapeutic leads and of great concern as lethal biological agents due to their high affinity and specificity for various receptors involved in neuromuscular transmission. Currently, of the approximately 6000 known conotoxin sequences, only about 3% have associated structural characterization, which leads to a bottleneck in rapid high-throughput screening (HTS) for identification of potential leads or threats. In this work, we combine a graph-based approach with homology modeling to expand the library of conotoxin structures and to identify those conotoxin sequences that are of the greatest value for experimental structural characterization. The latter would allow for the rapid expansion of the known structural space for generating high quality template-based models. Our approach generalizes to other evolutionarily-related, short, cysteine-rich venoms of interest. Overall, we present and validate an approach for venom structure modeling and experimental guidance and employ it to produce a 290%-larger library of approximate conotoxin structures for HTS. We also provide a set of ranked conotoxin sequences for experimental structure determination to further expand this library.

ConoMode, a database for conopeptide binding modes

ConoMode is a database for complex three-dimensional (3D) structures of conopeptides binding with their target proteins. Conopeptides, a large family of peptides from the venom of marine snails of the Conus genus, have exceptionally diverse sequences, and their high specificity to block ion channels makes them crucial as drug leads and tools for physiological studies. ConoMode is a specialized archive for the collection of 3D coordinate data for the conopeptides and their binding target proteins from published literature and the Protein Data Bank. These 3D structures can be determined using experimental methods such as X-ray crystallography and electron microscopy and computational methods including docking, homology modeling and molecular dynamics simulations. The binding modes for the conopeptides determined using computational modeling must be validated based on experimental data. The 3D coordinate data from ConoMode can be searched, visualized, downloaded and uploaded. Currently, ConoMode manages 19 conopeptide sequences (from 10 Conus species), 15 protein sequences and 37 3D structures. ConoMode utilizes a modern technical framework to provide a good user experience on mobile devices with touch interaction features. Furthermore, the database is fully optimized for unstructured data and flexible data models.

Database URL: http://conomode.qnlm.ac/conomode/conomode/index

Conotoxins: Chemistry and Biology

The venom of the marine predatory cone snails (genus Conus) has evolved for prey capture and defense, providing the basis for survival and rapid diversification of the now estimated 750+ species. A typical Conus venom contains hundreds to thousands of bioactive peptides known as conotoxins. These mostly disulfide-rich and well-structured peptides act on a wide range of targets such as ion channels, G protein-coupled receptors, transporters, and enzymes. Conotoxins are of interest to neuroscientists as well as drug developers due to their exquisite potency and selectivity, not just against prey but also mammalian targets, thereby providing a rich source of molecular probes and therapeutic leads. The rise of integrated venomics has accelerated conotoxin discovery with now well over 10,000 conotoxin sequences published. However, their structural and pharmacological characterization lags considerably behind. In this review, we highlight the diversity of new conotoxins uncovered since 2014, their three-dimensional structures and folds, novel chemical approaches to their syntheses, and their value as pharmacological tools to unravel complex biology. Additionally, we discuss challenges and future directions for the field.

Conotoxin Diversity in the Venom Gland Transcriptome of the Magician's Cone, Pionoconus magus

The transcriptomes of the venom glands of two individuals of the magician’s cone, Pionoconus magus, from Okinawa (Japan) were sequenced, assembled, and annotated. In addition, RNA-seq raw reads available at the SRA database from one additional specimen of P. magus from the Philippines were also assembled and annotated. The total numbers of identified conotoxin precursors and hormones per specimen were 118, 112, and 93. The three individuals shared only five identical sequences whereas the two specimens from Okinawa had 30 sequences in common. The total number of distinct conotoxin precursors and hormones for P. magus was 275, and were assigned to 53 conotoxin precursor and hormone superfamilies, two of which were new based on their divergent signal region. The superfamilies that had the highest number of precursors were M (42), O1 (34), T (27), A (18), O2 (17), and F (13), accounting for 55% of the total diversity. The D superfamily, previously thought to be exclusive of vermivorous cones was found in P. magus and contained a highly divergent mature region. Similarly, the A superfamily alpha 4/3 was found in P. magus despite the fact that it was previously postulated to be almost exclusive of the genus Rhombiconus. Differential expression analyses of P. magus compared to Chelyconus ermineus, the only fish-hunting cone from the Atlantic Ocean revealed that M and A2 superfamilies appeared to be more expressed in the former whereas the O2 superfamily was more expressed in the latter.

Transcriptomic Analysis of Marine Gastropod Hemifusus tuba Provides Novel Insights into Conotoxin Genes

The marine gastropod Hemifusus tuba is served as a luxury food in Asian countries and used in traditional Chinese medicine to treat lumbago and deafness. The lack of genomic data on H. tuba is a barrier to aquaculture development and functional characteristics of potential bioactive molecules are poorly understood. In the present study, we used high-throughput sequencing technologies to generate the first transcriptomic database of H. tuba. A total of 41 unique conopeptides were retrieved from 44 unigenes, containing 6-cysteine frameworks belonging to four superfamilies. Duplication of mature regions and alternative splicing were also found in some of the conopeptides, and the de novo assembly identified a total of 76,306 transcripts with an average length of 824.6 nt, of which including 75,620 (99.1%) were annotated. In addition, simple sequence repeats (SSRs) detection identified 14,000 unigenes containing 20,735 SSRs, among which, 23 polymorphic SSRs were screened. Thirteen of these markers could be amplified in Hemifusus ternatanus and seven in Rapana venosa. This study provides reports of conopeptide genes in Buccinidae for the first time as well as genomic resources for further drug development, gene discovery and population resource studies of this species.

Structural and Functional Analyses of Cone Snail Toxins

Cone snails are marine gastropod mollusks with one of the most powerful venoms in nature. The toxins, named conotoxins, must act quickly on the cone snails´ prey due to the fact that snails are extremely slow, reducing their hunting capability. Therefore, the characteristics of conotoxins have become the object of investigation, and as a result medicines have been developed or are in the trialing process. Conotoxins interact with transmembrane proteins, showing specificity and potency. They target ion channels and ionotropic receptors with greater regularity, and when interaction occurs, there is immediate physiological decompensation. In this review we aimed to evaluate the structural features of conotoxins and the relationship with their target types.