PowerNovo: de novo peptide sequencing via tandem mass spectrometry using an ensemble of transformer and BERT models.

The primary objective of analyzing the data obtained in a mass spectrometry-based proteomic experiment is peptide and protein identification, or correct assignment of the tandem mass spectrum to one amino acid sequence. Comparison of empirical fragment spectra with the theoretical predicted one or matching with the collected spectra library are commonly accepted strategies of proteins identification and defining of their amino acid sequences. Although these approaches are widely used and are appreciably efficient for the well-characterized model organisms or measured proteins, they cannot detect novel peptide sequences that have not been previously annotated or are rare. This study presents PowerNovo tool for de novo sequencing of proteins using tandem mass spectra acquired in a variety of types of mass analyzers and different fragmentation techniques. PowerNovo involves an ensemble of models for peptide sequencing: model for detecting regularities in tandem mass spectra, precursors, and fragment ions and a natural language processing model, which has a function of peptide sequence quality assessment and helps with reconstruction of noisy sequences. The results of testing showed that the performance of PowerNovo is comparable and even better than widely utilized PointNovo, DeepNovo, Casanovo, and Novor packages. Also, PowerNovo provides complete cycle of processing (pipeline) of mass spectrometry data and, along with predicting the peptide sequence, involves the peptide assembly and protein inference blocks.

Structural and functional characterization of cyclic pyrimidine-regulated anti-phage system.

3',5'-cyclic uridine monophosphate (cUMP) and 3',5'-cyclic cytidine monophosphate (cCMP) have been established as bacterial second messengers in the phage defense system, named pyrimidine cyclase system for anti-phage resistance (Pycsar). This system consists of a pyrimidine cyclase and a cyclic pyrimidine receptor protein. However, the molecular mechanism underlying cyclic pyrimidine synthesis and recognition remains unclear. Herein, we determine the crystal structures of a uridylate cyclase and a cytidylate cyclase, revealing the conserved residues for cUMP and cCMP production, respectively. In addition, a distinct zinc-finger motif of the uridylate cyclase is identified to confer substantial resistance against phage infections. Furthermore, structural characterization of cUMP receptor protein PycTIR provides clear picture of specific cUMP recognition and identifies a conserved N-terminal extension that mediates PycTIR oligomerization and activation. Overall, our results contribute to the understanding of cyclic pyrimidine-mediated bacterial defense.

Cloning and functional analysis of the juvenile hormone receptor gene CsMet in Coccinella septempunctata.

The potential role of the juvenile hormone receptor gene (methoprene-tolerant, Met) in reproduction of Coccinella septempunctata L. (Coleoptera: Coccinellidae)(Coleoptera: Coccinellidae), was investigated by cloning, analyzing expression profiles by quantitative real-time PCR, and via RNA interference (RNAi). CsMet encoded a 1518-bp open reading frames with a predicted protein product of 505 amino acids; the latter contained 2 Per-Arnt-Sim repeat profile at amino acid residues 30-83 and 102-175. CsMet was expressed in different C. septempunctata larvae developmental stages and was most highly expressed in third instar. CsMet expression in female adults gradually increased from 20 to 30 d, and expression levels at 25 and 30 d were significantly higher than levels at 1-15 d. CsMet expression in 20-d-old male adults was significantly higher than in males aged 1-15 d. CsMet expression levels in fat body tissues of male and female adults were significantly higher than expression in the head, thorax, and reproductive system. At 5 and 10 d after CsMet-dsRNA injection, CsMet expression was significantly lower than the controls by 75.05% and 58.38%, respectively. Ovary development and vitellogenesis in C. septempunctata injected with CsMet-dsRNA were significantly delayed and fewer mature eggs were produced. This study provides valuable information for the large-scale rearing of C. septempunctata.

The impact of exchanging the light and heavy chains on the structures of bovine ultralong antibodies.

The third complementary-determining regions of the heavy-chain (CDR3H) variable regions (VH) of some cattle antibodies are highly extended, consisting of 48 or more residues. These `ultralong' CDR3Hs form ß-ribbon stalks that protrude from the surface of the antibody with a disulfide cross-linked knob region at their apex that dominates antigen interactions over the other CDR loops. The structure of the Fab fragment of a naturally paired bovine ultralong antibody (D08), identified by single B-cell sequencing, has been determined to 1.6 Šresolution. By swapping the D08 native light chain with that of an unrelated antigen-unknown ultralong antibody, it is shown that interactions between the CDR3s of the variable domains potentially affect the fine positioning of the ultralong CDR3H; however, comparison with other crystallographic structures shows that crystalline packing is also a major contributor. It is concluded that, on balance, the exact positioning of ultralong CDR3H loops is most likely to be due to the constraints of crystal packing.

EFG-CS: Predicting chemical shifts from amino acid sequences with protein structure prediction using machine learning and deep learning models.

Nuclear magnetic resonance (NMR) crystallography is one of the main methods in structural biology for analyzing protein stereochemistry and structure. The chemical shift of the resonance frequency reflects the effect of the protons in a molecule producing distinct NMR signals in different chemical environments. Apprehending chemical shifts from NMR signals can be challenging since having an NMR structure does not necessarily provide all the required chemical shift information, making predictive models essential for accurately deducing chemical shifts, either from protein structures or, more ideally, directly from amino acid sequences. Here, we present EFG-CS, a web server that specializes in chemical shift prediction. EFG-CS employs a machine learning-based transfer prediction model for backbone atom chemical shift prediction, using ESMFold-predicted protein structures. Additionally, ESG-CS incorporates a graph neural network-based model to provide comprehensive side-chain atom chemical shift predictions. Our method demonstrated reliable performance in backbone atom prediction, achieving comparable accuracy levels with root mean square errors (RMSE) of 0.30 ppm for H, 0.22 ppm for Hα, 0.89 ppm for C, 0.89 ppm for Cα, 0.84 ppm for Cß, and 1.69 ppm for N. Moreover, our approach also showed predictive capabilities in side-chain atom chemical shift prediction achieving RMSE values of 0.71 ppm for Hß, 0.74-1.15 ppm for Hδ, and 0.58-0.94 ppm for Hγ, solely utilizing amino acid sequences without homology or feature curation. This work shows for the first time that generative AI protein models can predict NMR shifts nearly comparable to experimental models. This web server is freely available at https://biosig.lab.uq.edu.au/efg_cs, and the chemical shift prediction results can be downloaded in tabular format and visualized in 3D format.

Characterization of BrGH3A, a bovine rumen-derived glycoside hydrolase family 3 ß-glucosidase with a permuted domain arrangement.

The bovine rumen contains a large consortium of residential microbes that release a variety of digestive enzymes for feed degradation. However, the utilization of these microbial enzymes is still limited because these rumen microorganisms are mostly anaerobes and are thus unculturable. Therefore, we applied a sequence-based metagenomic approach to identify a novel 2,445-bp glycoside hydrolase family 3 ß-glucosidase gene known as BrGH3A from the metagenome of bovine ruminal fluid. BrGH3A ß-glucosidase is a 92-kDa polypeptide composed of 814 amino acid residues. Unlike most glycoside hydrolases in the same family, BrGH3A exhibited a permuted domain arrangement consisting of an (α/ß)6 sandwich domain, a fibronectin type III domain and a (ß/α)8 barrel domain. BrGH3A exhibited greater catalytic efficiency toward laminaribiose than cellobiose. We proposed that BrGH3A is an exo-acting ß-glucosidase from Spirochaetales bacteria that is possibly involved in the intracellular degradation of ß-1,3-/1,4-mixed linkage glucans that are present in grass cell walls. BrGH3A exhibits rich diversity in rumen hydrolytic enzymes and may represent a member of a new clan with a permuted domain topology within the large family.

Characterization of Synthetic Peptides by Mass Spectrometry.

In the quality control of synthetic peptides, mass spectroscopy (MS) serves as an optimal method for evaluating authenticity and integrity. Typically, the sequence of a synthetic peptide is already established, thereby directing the focus of analysis towards validating its identity and purity. This chapter outlines straightforward methodologies for conducting MS analyses specifically tailored for synthetic peptides.

Structural Characterization of Peptide Antibodies.

The role of proteins as very effective immunogens for the generation of antibodies is indisputable. Nevertheless, cases in which protein usage for antibody production is not feasible or convenient compelled the creation of a powerful alternative consisting of synthetic peptides. Synthetic peptides can be modified to obtain desired properties or conformation, tagged for purification, isotopically labeled for protein quantitation or conjugated to immunogens for antibody production. The antibodies that bind to these peptides represent an invaluable tool for biological research and discovery. To better understand the underlying mechanisms of antibody-antigen interaction, here, we present a pipeline developed by us to structurally classify immunoglobulin antigen binding sites and to infer key sequence residues and other variables that have a prominent role in each structural class.

Identification of a novel member of the genus Laulavirus (family Phenuiviridae) from the entomopathogenic ascomycete fungus Cordyceps javanica.

The virus family Phenuiviridae (order Hareavirales, comprising segmented negative-sense single stranded RNA viruses) has highly diverse members that are known to infect animals, plants, protozoans, and fungi. In this study, we identified a novel phenuivirus infecting a strain of the entomopathogenic fungus Cordyceps javanica isolated from a small brown plant hopper (Laodelphax striatellus), and this virus was tentatively named "Cordyceps javanica negative-strand RNA virus 1" (CjNRSV1). The CjNRSV1 genome consists of three negative-sense single stranded RNA segments (RNA1-3) with lengths of 7252, 2401, and 1117 nt, respectively. The 3'- and 5'-terminal regions of the RNA1, 2, and 3 segments have identical sequences, and the termini of the RNA segments are complementary to each other, reflecting a common characteristic of viruses in the order Hareavirales. RNA1 encodes a large protein (∼274 kDa) containing a conserved domain for the bunyavirus RNA-dependent RNA polymerase (RdRP) superfamily, with 57-80% identity to the RdRP encoded by phenuiviruses in the genus Laulavirus. RNA2 encodes a protein (∼79 kDa) showing sequence similarity (47-63% identity) to the movement protein (MP, a plant viral cell-to-cell movement protein)-like protein (MP-L) encoded by RNA2 of laulaviruses. RNA3 encodes a protein (∼28 kDa) with a conserved domain of the phenuivirid nucleocapsid protein superfamily. Phylogenetic analysis using the RdRPs of various phenuiviruses and other unclassified phenuiviruses showed CjNRSV1 to be grouped with established members of the genus Laulavirus. Our results suggest that CjNRSV1 is a novel fungus-infecting member of the genus Laulavirus in the family Phenuiviridae.

Complete genome analysis of a novel hypovirus in the phytopathogenic fungus Monilinia fructicola.

Monilinia fructicola is one of the most devastating fungal diseases of rosaceous fruit crops, both in the field and postharvest, causing significant yield losses. Here, we report the discovery of a novel positive single-stranded RNA virus, Monilinia fructicola hypovirus 3 (MfHV3), in a strain (hf-1) of the phytopathogenic fungus Monilinia fructicola. The complete genome of MfHV3 is 9259 nucleotides (nt) in length and contains a single large open reading frame (ORF) from nt position 462 to 8411. This ORF encodes a polyprotein with three conserved domains, namely UDP-glycosyltransferase, RNA-dependent RNA polymerase (RdRp), and DEAD-like helicase. The MfHV3 polyprotein shares the highest similarity with Colletotrichum camelliae hypovirus 1. Phylogenetic analysis indicated that MfHV3 clustered with members of the genus Betahypovirus within the family Hypoviridae. Taken together, the results of genomic organization comparisons, amino acid sequence alignments, and phylogenetic analysis convincingly show that MfHV3 is a new member of the genus Betahypovirus, family Hypoviridae.

Identification of heparin-binding amino acid residues in antibody HS4C3 with the potential to design antibodies against heparan sulfate domains.

Heparan sulfate (HS) is a linear polysaccharide with high structural and functional diversity. Detection and localization of HS in tissues can be performed using single chain variable fragment (scFv) antibodies. Although several anti-HS antibodies recognizing different sulfation motifs have been identified, little is known about their interaction with HS. In this study the interaction between the scFv antibody HS4C3 and heparin was investigated. Heparin-binding lysine and arginine residues were identified using a protect and label methodology. Site-directed mutagenesis was applied to further identify critical heparin-binding lysine/arginine residues using immunohistochemical and biochemical assays. In addition, computational docking of a heparin tetrasaccharide towards a 3-D homology model of HS4C3 was applied to identify potential heparin-binding sites. Of the 12 lysine and 15 arginine residues within the HS4C3 antibody, 6 and 9, respectively, were identified as heparin-binding. Most of these residues are located within one of the complementarity determining regions (CDR) or in their proximity. All basic amino acid residues in the CDR3 region of the heavy chain were involved in binding. Computational docking showed a heparin tetrasaccharide close to these regions. Mutagenesis of heparin-binding residues reduced or altered reactivity towards HS and heparin. Identification of heparin-binding arginine and lysine residues in HS4C3 allows for better understanding of the interaction with HS and creates a framework to rationally design antibodies targeting specific HS motifs.

Correlation between the Colony Phenotype and Amino Acid Sequence of the Variable Vaa Antigen in Clinical Isolates of Mycoplasma hominis.

A new Mycoplasma hominis phenotype forming mini-colonies (MC) on agar and distinct from the phenotype forming typical colonies (TC) not only in size, but also in morphology, growth rate, and resistance to adverse factors, has been previously identified. In this study, the phenotype of colonies was determined and a comparative analysis of the amino acid sequence of the main variable antigen Vaa of the laboratory strain N-34 and seven clinical isolates of M. hominis was performed. It is demonstrated that the amino acid sequence of Vaa in clinical isolates forming TC (similar to the laboratory strain N-34) is entirely analogous to that of laboratory strain. Clinical isolates forming MC carry amino acid substitutions in the variable C-terminal region of Vaa, which can contribute to adhesion to eukaryotic cells and immune evasion. The connection between colony phenotype and amino acid sequence of Vaa is established.

Protein SUMOylation promotes cAMP-independent EPAC1 activation.

Protein SUMOylation is a prevalent stress-response posttranslational modification crucial for maintaining cellular homeostasis. Herein, we report that protein SUMOylation modulates cellular signaling mediated by cAMP, an ancient and universal stress-response second messenger. We identify K561 as a primary SUMOylation site in exchange protein directly activated by cAMP (EPAC1) via site-specific mapping of SUMOylation using mass spectrometry. Sequence and site-directed mutagenesis analyses reveal that a functional SUMO-interacting motif in EPAC1 is required for the binding of SUMO-conjugating enzyme UBC9, formation of EPAC1 nuclear condensate, and EPAC1 cellular SUMOylation. Heat shock-induced SUMO modification of EPAC1 promotes Rap1/2 activation in a cAMP-independent manner. Structural modeling and molecular dynamics simulation studies demonstrate that SUMO substituent on K561 of EPAC1 promotes Rap1 interaction by increasing the buried surface area between the SUMOylated receptor and its effector. Our studies identify a functional SUMOylation site in EPAC1 and unveil a novel mechanism in which SUMOylation of EPAC1 leads to its autonomous activation. The findings of SUMOylation-mediated activation of EPAC1 not only provide new insights into our understanding of cellular regulation of EPAC1 but also will open up a new field of experimentation concerning the cross-talk between cAMP/EPAC1 signaling and protein SUMOylation, two major cellular stress response pathways, during cellular homeostasis.

Differential potassium channel inhibitory activities of a novel thermostable degradation peptide BmKcug1a-P1 from scorpion medicinal material and its N-terminal truncated/restored peptides.

Thermally stable full-length scorpion toxin peptides and partially degraded peptides with complete disulfide bond pairing are valuable natural peptide resources in traditional Chinese scorpion medicinal material. However, their pharmacological activities are largely unknown. This study discovered BmKcug1a-P1, a novel N-terminal degraded peptide, in this medicinal material. BmKcug1a-P1 inhibited hKv1.2 and hKv1.3 potassium channels with IC50 values of 2.12 ± 0.27 µM and 1.54 ± 0.28 µM, respectively. To investigate the influence of N-terminal amino acid loss on the potassium channel inhibiting activities, three analogs (i.e., full-length BmKcug1a, BmKcug1a-P1-D2 and BmKcug1a-P1-D4) of BmKcug1a-P1 were prepared, and their potassium channel inhibiting activities on hKv1.3 channel were verified by whole-cell patch clamp technique. Interestingly, the potassium channel inhibiting activity of full-length BmKcug1a on the hKv1.3 channel was significantly improved compared to its N-terminal degraded form (BmKcug1a-P1), while the activities of two truncated analogs (i.e., BmKcug1a-P1-D2 and BmKcug1a-P1-D4) were similar to that of BmKcug1a-P1. Extensive alanine-scanning experiments identified the bonding interface (including two key functional residues, Asn30 and Arg34) of BmKcug1a-P1. Structural and functional dissection further elucidated whether N-terminal residues of the peptide are located at the bonding interface is important in determining whether the N-terminus significantly influences the potassium channel inhibiting activity of the peptide. Altogether, this research identified a novel N-terminal degraded active peptide, BmKcug1a-P1, from traditional Chinese scorpion medicinal material and elucidated how the N-terminus of peptides influences their potassium channel inhibiting activity, contributing to the functional identification and molecular truncation optimization of full-length and degraded peptides from traditional Chinese scorpion medicinal material Buthus martensii Karsch.

Insights into Actin Isoform-Specific Interactions with Myosin via Computational Analysis.

Actin, which plays a crucial role in cellular structure and function, interacts with various binding proteins, notably myosin. In mammals, actin is composed of six isoforms that exhibit high levels of sequence conservation and structural similarity overall. As a result, the selection of actin isoforms was considered unimportant in structural studies of their binding with myosin. However, recent high-resolution structural research discovered subtle structural differences in the N-terminus of actin isoforms, suggesting the possibility that each actin isoform may engage in specific interactions with myosin isoforms. In this study, we aimed to explore this possibility, particularly by understanding the influence of different actin isoforms on the interaction with myosin 7A. First, we compared the reported actomyosin structures utilizing the same type of actin isoforms as the high-resolution filamentous skeletal α-actin (3.5 Å) structure elucidated using cryo-EM. Through this comparison, we confirmed that the diversity of myosin isoforms leads to differences in interaction with the actin N-terminus, and that loop 2 of the myosin actin-binding sites directly interacts with the actin N-terminus. Subsequently, with the aid of multiple sequence alignment, we observed significant variations in the length of loop 2 across different myosin isoforms. We predicted that these length differences in loop 2 would likely result in structural variations that would affect the interaction with the actin N-terminus. For myosin 7A, loop 2 was found to be very short, and protein complex predictions using skeletal α-actin confirmed an interaction between loop 2 and the actin N-terminus. The prediction indicated that the positively charged residues present in loop 2 electrostatically interact with the acidic patch residues D24 and D25 of actin subdomain 1, whereas interaction with the actin N-terminus beyond this was not observed. Additionally, analyses of the actomyosin-7A prediction models generated using various actin isoforms consistently yielded the same results regardless of the type of actin isoform employed. The results of this study suggest that the subtle structural differences in the N-terminus of actin isoforms are unlikely to influence the binding structure with short loop 2 myosin 7A. Our findings are expected to provide a deeper understanding for future high-resolution structural binding studies of actin and myosin.

A Bicyclic Analog of the Linear Peptide Arodyn Is a Potent and Selective Kappa Opioid Receptor Antagonist.

Kappa opioid receptor (KOR) antagonists have potential therapeutic applications in the treatment of stress-induced relapse to substance abuse and mood disorders. The dynorphin A analog arodyn (Ac[Phe1,2,3,Arg4,D-Ala8]dynorphin A-(1-11)-NH2) exhibits potent and selective kappa opioid receptor antagonism. Multiple cyclizations in longer peptides, such as dynorphin and its analogs, can extend the conformational constraint to additional regions of the peptide beyond what is typically constrained by a single cyclization. Here, we report the design, synthesis, and pharmacological evaluation of a bicyclic arodyn analog with two constraints in the opioid peptide sequence. The peptide, designed based on structure-activity relationships of monocyclic arodyn analogs, was synthesized by solid-phase peptide synthesis and cyclized by sequential ring-closing metathesis (RCM) in the C- and N-terminal sequences. Molecular modeling studies suggest similar interactions of key aromatic and basic residues in the bicyclic peptide with KOR as found in the cryoEM structure of KOR-bound dynorphin, despite substantial differences in the backbone conformations of the two peptides. The bicyclic peptide's affinities at KOR and mu opioid receptors (MOR) were determined in radioligand binding assays, and its KOR antagonism was determined in the [35S]GTPγS assay in KOR-expressing cells. The bicyclic analog retains KOR affinity and selectivity (Ki = 26 nM, 97-fold selectivity over MOR) similar to arodyn and exhibits potent KOR antagonism in the dynorphin-stimulated [35S]GTPγS assay. This bicyclic peptide represents a promising advance in preparing cyclic opioid peptide ligands and opens avenues for the rational design of additional bicyclic opioid peptide analogs.

Exploring the Structurally Conserved Regions and Functional Significance in Bacterial N-Terminal Nucleophile (Ntn) Amide-Hydrolases.

The initial adoption of penicillin as an antibiotic marked the start of exploring other compounds essential for pharmaceuticals, yet resistance to penicillins and their side effects has compromised their efficacy. The N-terminal nucleophile (Ntn) amide-hydrolases S45 family plays a key role in catalyzing amide bond hydrolysis in various compounds, including antibiotics like penicillin and cephalosporin. This study comprehensively analyzes the structural and functional traits of the bacterial N-terminal nucleophile (Ntn) amide-hydrolases S45 family, covering penicillin G acylases, cephalosporin acylases, and D-succinylase. Utilizing structural bioinformatics tools and sequence analysis, the investigation delineates structurally conserved regions (SCRs) and substrate binding site variations among these enzymes. Notably, sixteen SCRs crucial for substrate interaction are identified solely through sequence analysis, emphasizing the significance of sequence data in characterizing functionally relevant regions. These findings introduce a novel approach for identifying targets to enhance the biocatalytic properties of N-terminal nucleophile (Ntn) amide-hydrolases, while facilitating the development of more accurate three-dimensional models, particularly for enzymes lacking structural data. Overall, this research advances our understanding of structure-function relationships in bacterial N-terminal nucleophile (Ntn) amide-hydrolases, providing insights into strategies for optimizing their enzymatic capabilities.

TriplEP-CPP: Algorithm for Predicting the Properties of Peptide Sequences.

Advancements in medicine and pharmacology have led to the development of systems that deliver biologically active molecules inside cells, increasing drug concentrations at target sites. This improves effectiveness and duration of action and reduces side effects on healthy tissues. Cell-penetrating peptides (CPPs) show promise in this area. While traditional medicinal chemistry methods have been used to develop CPPs, machine learning techniques can speed up and reduce costs in the search for new peptides. A predictive algorithm based on machine learning models was created to identify novel CPP sequences using molecular descriptors using a combination of algorithms like k-nearest neighbors, gradient boosting, and random forest. Some potential CPPs were found and tested for cytotoxicity and penetrating ability. A new low-toxicity CPP was discovered from the Rhopilema esculentum venom proteome through this study.

C(P)XCG Proteins of Haloferax volcanii with Predicted Zinc Finger Domains: The Majority Bind Zinc, but Several Do Not.

In recent years, interest in very small proteins (µ-proteins) has increased significantly, and they were found to fulfill important functions in all prokaryotic and eukaryotic species. The halophilic archaeon Haloferax volcanii encodes about 400 µ-proteins of less than 70 amino acids, 49 of which contain at least two C(P)XCG motifs and are, thus, predicted zinc finger proteins. The determination of the NMR solution structure of HVO_2753 revealed that only one of two predicted zinc fingers actually bound zinc, while a second one was metal-free. Therefore, the aim of the current study was the homologous production of additional C(P)XCG proteins and the quantification of their zinc content. Attempts to produce 31 proteins failed, underscoring the particular difficulties of working with µ-proteins. In total, 14 proteins could be produced and purified, and the zinc content was determined. Only nine proteins complexed zinc, while five proteins were zinc-free. Three of the latter could be analyzed using ESI-MS and were found to contain another metal, most likely cobalt or nickel. Therefore, at least in haloarchaea, the variability of predicted C(P)XCG zinc finger motifs is higher than anticipated, and they can be metal-free, bind zinc, or bind another metal. Notably, AlphaFold2 cannot correctly predict whether or not the four cysteines have the tetrahedral configuration that is a prerequisite for metal binding.

Overview and Evolution of Insect Fibroin Heavy Chain (FibH).

The FibH gene, crucial for silk spinning in insects, encodes a protein that significantly influences silk fiber mechanics. Due to its large size and repetitive sequences, limited known sequences of insect FibH impede comprehensive understanding. Here, we analyzed 114 complete FibH gene sequences from Lepidoptera (71 moths, 24 butterflies) and 13 Trichoptera, revealing single-copy FibH in most species, with 2-3 copies in Hesperinae and Heteropterinae (subfamily of skippers). All FibH genes are structured with two exons and one intron (39-45 bp), with the second exon being notably longer. Moths exhibit higher GC content in FibH compared to butterflies and Trichoptera. The FibH composition varies among species, with moths and butterflies favoring Ala, Gly, Ser, Pro, Gln, and Asn, while Trichoptera FibH is enriched in Gly, Ser, and Arg, and has less Ala. Unique to Trichoptera FibH are Tyr, Val, Arg, and Trp, whereas Lepidoptera FibH is marked by polyAla (polyalanine), polySer (polyserine), and the hexapeptide GAGSGA. A phylogenetic analysis suggests that Lepidoptera FibH evolved from Trichoptera, with skipper FibH evolving from Papilionoidea. This study substantially expands the FibH repertoire, providing a foundation for the development of artificial silk.