Protein identification and analysis tools in the ExPASy server

Genome-wide identification and analysis of the NF-Y transcription factor family reveal its potential roles in tobacco (Nicotiana tabacum L.)

Nuclear Factor Y (NF-Y) represents a group of transcription factors commonly present in higher eukaryotes, typically consisting of three subunits: NF-YA, NF-YB, and NF-YC. They play crucial roles in the embryonic development, photosynthesis, flowering, abiotic stress responses, and other essential processes in plants. To better understand the genome-wide NF-Y domain-containing proteins, the protein physicochemical properties, chromosomal localization, synteny, phylogenetic relationships, genomic structure, promoter cis-elements, and protein interaction network of NtNF-Ys in tobacco (Nicotiana tabacum L.) were systematically analyzed. In this study, we identified 58 NtNF-Ys in tobacco, respectively, and divided into three subfamilies corresponding to their phylogenetic relationships. Their tissue specificity and expression pattern analyses for leaf development, drought and saline-alkali stress, and ABA response were carried out using RNA-seq or qRT-PCR. These findings illuminate the role of NtNF-Ys in regulating plant leaf development, drought and saline-alkali stress tolerance, and ABA response. This study offers new insights to enhance our understanding of the roles of NtNF-Ys and identify potential genes involved in leaf development, as well as drought and saline-alkali stress tolerance of plants.

KNUTS-DB - a data-driven knowledge database for NUTS

OBJECTIVES

Databases specifying the nutrients and allergens present in multi-ingredient foods are required to explore the effect of food consumption on health outcomes accurately. The phytochemicals found in tree nuts have been associated with antioxidant, anti-inflammatory, antiproliferative, antiviral, chemo preventive and hypercholesterolaemic actions, all of which are known to affect the initiation and progression of several pathogenic processes. We have developed KNUTS-DB - a data-driven knowledge database for nuts - containing information on the chemical nutrients, molecular targets, pathways, disease associations, and clinically defined food allergen properties of peanuts and tree nuts.

METHODS

The database includes data sets associated with extremely rich and diverse metadata on almonds, cashew, pecan, walnut, pistachio, peanut and walnut. Additionally, the database allows users to perform pathway and GO-Term enrichment analysis for user-defined chemical nutrients. The results of the analysis are presented using heatmaps and network plots.

RESULTS

The database can be searched using options such as similarity search, interaction values, type of allergens, and specific nut types. KNUTS-DB offers researchers a unique platform to explore nuts' chemical composition and to investigate associations between nut composition and health outcomes - providing deeper insights into the molecular mechanisms.

CONCLUSIONS

KNUTS-DB, the first of its kind, is freely available to all users via https://allergypred.charite.de/KNutsDB/ without any login or registration.

Isoelectric point, net charge and amino acid analysis of experimentally validated therapeutic antibodies

The isoelectric point (pI) of an antibody is known to affect its non-specific interactions and repulsive self-interactions. However, analytical outcomes for the pI of a large number of therapeutic antibodies remain unexplored. In this study, we explored the pI and net charge of variable heavy (VH), variable light (VL), CDR (complementarity determining regions) and whole IgG on a large number of therapeutic antibodies, additionally amino acids distribution in the CDR regions were also analyzed. A total of 708 experimentally validated antibodies from the Thera-SAbDab database were analyzed in this study. Analysis of the antibody dataset showed that the pI of the whole IgG sequence is between 5 and 9, while the majority was in the intermediate range between 7 and 9 (86.7%). The charge had a wide range from - 10 to 12, with the majority falling between the charges 2-6 (53.4%). However, the combined pI score of the CDRs of light chains (60%) as well as for the heavy chains (67%) was observed in the range of 4-6. The amino acid composition analysis of CDR regions revealed that most of the amino acids in the light chain are uncharged-polar (46.3%) followed by hydrophobic-aliphatic (28.4%), while in the heavy chain; it is hydrophobic-aliphatic (35.2%) followed by uncharged-polar (24.6%). In conclusion, the pI and net charge analysis of therapeutic antibodies are crucial for understanding pharmacokinetic properties. Moreover, amino acid composition of the light and heavy chain CDR regions has a significant impact on the pI and charge of the entire IgG antibody.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-025-00356-y.

Aptamer- vs Fab-Conjugated Liposomes: A Comparative Study in Targeting Acute Myeloid Leukemia Cells

Acute myeloid leukemia (AML) is a hematologic malignancy characterized by uncontrolled proliferation of abnormal myeloid cells with a generally poor prognosis despite advancements in chemotherapy and stem cell transplantation. To enhance therapeutic efficacy and minimize systemic toxicity, we designed liposomal nanoparticles functionalized with two distinct targeting ligands, a DNA aptamer or fragment-antigen-binding (Fab) antibody, targeting the surface marker transmembrane glycoprotein CD33 antigen (CD33) on AML cells. Aptamer- and Fab-conjugated liposomes (Apt-Lipm and Fab-Lipm, respectively) were prepared and tested for cellular uptake by CD33-positive AML cell lines. Comparative studies revealed that Fab-Lipm exhibited significantly superior binding affinity, targeting efficiency, and cellular uptake compared with Apt-Lipm. Furthermore, we demonstrated the intracellular distribution and endocytic pathways of Fab-Lipm during the cellular uptake. This comparative study of aptamer- and Fab-conjugated liposomes suggests that the Fab-conjugated liposomal system offers enhanced precision in targeting AML cells for the development of effective therapeutic strategies against hematologic malignancies.

De novo Biosynthesis of Caffeic Acid and Chlorogenic Acid in Escherichia coli via Enzyme Engineering and Pathway Engineering

Caffeic acid (CA) and chlorogenic acid (CGA) have diverse health benefits, including hemostatic, antioxidant, and antiinflammatory, highlighting their potential for medical applications. However, the absence of high-performance production strains increases production costs, limiting their wider application. In this study, we engineered Escherichia coli for the de novo production of CA and CGA. To improve production, a highly efficient mutant tyrosine ammonia-lyase from Rhodotorula taiwanensis (RtTALT415M/Y458F) was identified using genome mining and protein engineering. By engineering the tyrosine biosynthetic pathway through the deletion of pheA and tyrR, along with the overexpression of aroGfbr and tyrAfbr, we developed an engineered E. coli strain, CA11, which produced 6.36 g/L of CA with a yield of 0.06 g/g glucose and a productivity of 0.18 g/L/h. This represents the highest titer reported for microbial synthesis of CA using glucose as the sole carbon source in E. coli. Based on strain CA11, we further developed strain CGA13, with optimized replicons, promoters, and ribosome-binding sites, which produced 1.53 g/L of CGA in fed-batch fermentation, highlighting its potential for industrial-scale production.

Development of a multi-epitope vaccine from outer membrane proteins and identification of novel drug targets against Francisella tularensis: an In Silico approach

Background

Francisella tularensis is a category A potential thread agent, making the development of vaccines and countermeasures a high priority. Therefore, identifying new vaccine candidates and novel drug targets is essential for addressing this significant public health concern.

Methods

This study presents an in silico analysis of two strategies against F. tularensis infection: the development of a multi-epitope vaccine (MEV) and the identification of novel drug targets. Outer membrane proteins (OMPs) were predicted using subcellular localization tools and immunogenicity was evaluated using a reverse vaccinology pipeline. Epitopes from these OMPs were combined to create candidate MEV for prophylactic protection. Concurrently, cytoplasmic proteins were subjected to rigorous analysis to identify potential novel drug targets.

Results

Of 1,921 proteins, we identified 12 promising protein vaccine candidates from F. tularensis OMPs and proposed a multi-epitope vaccine (MEV) designed using seven immunodominant epitopes derived from four of these OMPs, including two hypothetical proteins (WP_003026145.1 and WP_003029346.1), an OmpA family protein (WP_003020808.1), and PD40 (WP_003021546.1). In addition, we proposed 10 novel drug targets for F. tularensis: Asp-tRNA (Asn)/Glu-tRNA (Gln) amidotransferase subunit GatC (WP_003017413.1), NAD(P)-binding protein (WP_042522581.1), 30S ribosomal protein S16 (WP_003023081.1), Class I SAM-dependent methyltransferase (WP_003022345.1), haloacid dehalogenase (WP_003014157.1), uroporphyrinogen-III synthase (WP_003022220.1), and four hypothetical proteins (WP_003017784.1, WP_003020080.1, WP_003020066.1, and WP_003022350.1).

Conclusion

This study designed an MEV and proposed novel drug targets to address tularemia, offering broad protection against various F. tularensis strains. MEV, with favorable physicochemical properties, showed strong potential through molecular docking and dynamic simulations. Immune simulations suggest that it may elicit robust responses against pathogens. The identification of novel drug targets can lead to the discovery of new antimicrobial agents. However, further in vitro and in vivo studies are required to validate their efficacy and capability.

Immunoinformatics method to design universal multi-epitope nanoparticle vaccine for TGEV S protein

Porcine transmissible enteritis virus (TGEV) is a fatal pathogen affecting newborn piglets, presenting a significant challenge to global intensive pig farming biosecurity due to its ongoing mutation. There is still a lack of effective vaccines to combat this virus, Vaccination has long been considered the most effective way to overcome infectious diseases, however, traditional vaccines cannot be brought to market quickly enough to deal with rapid mutations and emerging viruses. Therefore, this study addresses this gap by using immunoinformatics methods and ferritin nanoparticle delivery system to build a platform for rapid research and development of porcine coronavirus vaccine, designing a candidate nanoparticle vaccine that targets the TGEV S protein. To this end, multiple servers and strict screening criteria were used to analyze the S protein, and 3 CTL dominant epitopes, 3 Th dominant epitopes, and 6 B cell dominant epitopes were obtained. The candidate nanoparticle vaccine was constructed by incorporating ferritin sequences through the C-terminus after they were tandemly linked in a certain order using a flexible linker. Further experimental analyses showed that the designed candidate nanoparticle vaccine possessed relatively high antigenicity, immunogenicity, non-allergenicity, non-transmembrane proteins, suitable physicochemical properties, and high solubility upon overexpression. Tertiary structure modeling and disulfide engineering ensured conformational similarity to natural proteins and high stability. Additionally, the model predicted 6 Linear Epitopes and 6 Discontinuous Epitopes for B-cell conformational epitopes. Docking with TLR-3 and TLR-4 molecules shows a large number of interacting hydrogen-bonded amino acid residues and hydrophobically interacting amino acid residues. Immunomimetic assays show high levels of immunoglobulin, T-lymphocyte and IFN-γ secretion and may elicit specific immune responses. Through computerized cloning, the candidate nanoparticle vaccine can be efficiently expressed in the E. coli K12 expression system, aligning with future large-scale industrial production strategies. Overall, the results indicate that the constructed candidate nanoparticle vaccine can be effectively expressed and may be able to induce a strong immune response, which is expected to be an ideal candidate vaccine against TGEV.

Design of a peptide-based vaccine against human respiratory syncytial virus using a reverse vaccinology approach: evaluation of immunogenicity, antigenicity, allergenicity, and toxicity

Background

Attempts to develop an hRSV vaccine have faced safety and efficacy challenges, with only three FDA-approved vaccines (Moderna's Mresvia, Pfizer's Abrysvo, and GSK's Arexvy) available. These vaccines are limited to individuals over 60 years, require boosters, and only reduce disease severity without clearing the infection. Therefore, we employed a reverse vaccinology approach in this study to identify the most promising antigenic epitopes capable of eliciting a robust and protective immune response.

Methodology

This study employed computational techniques to design a novel multi-epitope vaccine targeting hRSV. Using bioinformatics tools, candidate epitopes were identified from conserved viral proteins (F and G glycoproteins), assessing their immunogenicity, antigenicity, and allergenicity. Key tools included ExPASy, ProtParam, VaxiJen v2.0, AllergenFP v1.0, AllerTOP v2.0, NetCTL v1.2, IEDB, and Toxin-Pred. The vaccine construct was assessed for stability and toxicity through in silico analyses. We then characterized its kinetic properties, evaluated its structural integrity, and analyzed its interactions with Toll-like receptors (TLRs) using molecular docking, modeling, and refinement with AlphaFold3 and ClusPro.

Results

The designed constructs showed strong antigenicity (0.5996 for F-based and 0.6048 for G-based vaccine), non-allergenicity, and stability (instability index <40). Among these, most amino acids were in the extracellular domain of the construct. Molecular docking and dynamics simulations indicated strong binding interactions with TLR1 and TLR4 and minimal RMSF fluctuations, which ensured structural stability. Strong humoral and cellular responses were suggested by in silico immune simulation demonstrating robust immune activation, with high levels of IgG, IgM, IL-2, and IFN-γ. The physical and chemical analyses revealed that the majority of amino acids from the F and G proteins were located in the extracellular domain of the construct. The presence of signal peptide cleavage sites in both glycoprotein components further facilitates antigen presentation to the immune system.

Conclusions

This study presents a promising peptide-based vaccine candidate against hRSV that can effectively engage the immune system, showing strong immunogenicity and antigenicity. Future in vitro and in vivo studies are essential to evaluate the ability of the multi-epitope vaccine candidate to stimulate both humoral and cell-mediated immune responses and to assess its efficacy and safety profile.

Strictly G-Specific Alginate Lyase Aly7Sa for Efficient Preparation of Unsaturated Guluronate Oligosaccharides

Alginate is a commercially valuable polysaccharide consisting of β-d-mannuronate (M) and its C5 epimer, α-l-guluronate (G). Alginate lyases are efficient tools for the degradation of alginate and the preparation of oligosaccharides. In this research, an endolytic alginate lyase Aly7Sa with strict G specificity was expressed and characterized with the optimum reaction conditions at 30 °C and pH 6.5. The main degradation products of Aly7Sa for alginate were trisaccharide to octasaccharide, and those of PolyG were disaccharide to heptasaccharide. By utilizing HPAEC-PAD/MS and NMR methods, we identified the structure of products obtained from alginate. Interestingly, the trisaccharide to hexasaccharide products of Aly7Sa contained only unsaturated guluronate oligosaccharides, which were different from all of the characterized G-specific alginate lyases. The absence of oligosaccharide products with M residues demonstrated the strict G specificity of Aly7Sa. The targeted preparation was carried out based on the regular oligosaccharide pattern of Aly7Sa. By single-step purification employing gel-permeation chromatography, 4.8 mg of ΔGG, 6.8 mg of ΔGGG, and 3.7 mg of ΔGGGG were obtained with 100 mg of alginate as substrate. The strictly G-specific alginate lyase Aly7Sa provided an efficient tool for the preparation of unsaturated guluronate oligosaccharides, and the unique substrate specificity of the enzyme could also serve the research and development of alginate.

Genome-wide Identification, Characterization, and Expression Analysis of NHX Genes in Phaseolus vulgaris L. under Salt Stress: An In Silico Approach

BACKGROUND

Climate change is among the major triggering agents of abiotic stresses (e.g., saline stress), culminating in a vulnerability of common bean production systems. In recent decades, important research has identified and characterized genes that can mitigate the adverse effects caused by salt stress; among them, the Na+/H+ antiporters (NHXs) gene stands out. The NHX genes are widely distributed in all organisms and play significant roles in osmotic regulation in plants under salt stress conditions. Genome-wide identification of NHX genes has been performed in several plant species but not in Phaseolus vulgaris L.

METHODS

This study aimed to identify and characterize NHX genes in P. vulgaris L. using a genome-wide analysis approach conducted in silico. The common bean genome revealed nine putative PvNHX genes, and their subcellular localization, phylogenetic relationship, cis-regulatory elements, conserved motifs identification, chromosomal location, expression patterns, and interaction networks were analyzed.

RESULTS

Promoter analysis suggested that PvNHX genes shared hormone-related elements and were light-responsive and stress-responsive. Seven PvNHX genes were under the regulation of five microRNA (miRNA) families. RNA-seq analysis revealed that most PvNHX genes were expressed in response to salt stress. Currently, the most assertive strategy to confront these adversities is to use the information generated by sequencing plants to identify candidate genes that can be introgressed to improve programs in producing resilient cultures.

CONCLUSION

These results can provide valuable information for future studies on the functional mechanism of PvNHX genes in common beans in response to salt stress.

The complete mitochondrial genome analysis of Haemaphysalis hystricis Supino, 1897 (Ixodida: Ixodidae) and its phylogenetic implications

In order to study the sequence characteristics, gene order, and codon usage of the mitochondrial genome of Haemaphysalis hystricis, and to explore its phylogenetic relationship, a total of 36 H. hystricis isolated from dogs were used as sample in this study. The mitochondrial genome of a H. hystricis was amplified with several pairs of specific primers by PCR, and was sequenced by first generation sequencing. The mitochondrial genome of H. hystricis was 14,719 bp in size, and it contained 37 genes including 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and AT-rich region. Each PCG sequence had different lengths, the sequence longest and shortest gene were nad5 (1,652 bp) and atp8 (155 bp), respectively, among the 13 PCGs. All PCGs used ATN as their initiation codon, 10 of 13 PCGs used TAN as their termination codon, and 3 of which had incomplete termination codon (TA/T). Most of the 22 tRNAs with different sizes could form the classical cloverleaf structures expect for tRNA-Ala, tRNA-Ser1, tRNA-Ser2, and tRNA-Glu, and there were base mismatch (U-U and U-G) in all the 22 tRNAs sequences. Two rRNAs, namely rrnL and rrnS, had different lengths, rrnL located between tRNA-Leu1 and tRNA-Val, and rrnS located between tRNA-Val and tRNA-Ile, respectively. Two AT (D-loop) control areas with different lengths were in the mitochondrial genome, the NCRL was located between tRNA-Leu2 and tRNA-Cys, and the NCRS was located between rrnS and tRNA-Ile. The complete mitochondrial genome sequence of H. hystricis was AT preferences, and the gene order is the same as that of other Haemaphysalis family ticks. However, phylogenetic analysis showed that H. hystricis was most closely related to Haemaphysalis longicornis among the selected ticks. The mitochondrial genome not only enriches the genome database, provides more novel genetic markers for identifying tick species, and studying its molecular epidemiology, population genetics, systematics, but also have implications for the diagnosis, prevention, and control of ticks and tick-borne diseases in animals and humans.

Structural-functional analyses of the huntingtin/HAP40 complex in Drosophila and humans

Huntington's disease (HD) is a neurodegenerative disorder caused by an abnormal CAG expansion in the Huntingtin (HTT) gene. Given its simple genetic cause but complex pathogenic mechanisms, interest in targeting HTT for HD treatment is growing, necessitating a clear understanding of HTT regulation. HTT protein primarily exists in a core complex with HAP40, forming a highly ordered structure with two large globular domains connected by a bridge. We previously demonstrated that HAP40 is conserved in Drosophila, controls HTT's function, protein stability, and levels, and is a potential modifier of HD pathogenesis, supporting its central role in HTT regulation. Here, we showed that HTT synergizes with HAP40 to induce novel gain-of-function effects in Drosophila when overexpressed. Protein modeling revealed that despite their prominent evolutionary and sequence divergence, the fly and human HTT-HAP40 complexes share a high degree of structural similarity. Protein-contact maps and molecular simulations showed that HAP40 preferentially binds to HTT's C-terminal domain in both complexes. By examining the interfacial contacts between HTT and HAP40 in fly and human complexes, we identified ten conserved bonds that are important for HAP40's affinity for HTT. Finally, we showed that the conserved N-terminal BΦ motif in HAP40 is not essential for HTT binding but important for HAP40's functions. Through the structural-functional analyses of the fly and human HTT-HAP40 complexes, our results support that the structural similarity underlies the functional conservation of the two complexes from these evolutionarily distant species and further uncover novel insight into HAP40 regulation and its interaction with HTT.

Genome-wide identification and expression analysis of C3H gene family in melon

Zinc finger protein (ZFP) represent a significant class of transcription factors in plants, involved in various functions, including tissue development, signal transduction, and responses to both biotic and abiotic stresses. ZFPs are categorized into 10 distinct subfamilies, among which the C3H gene family is recognized as a functionally significant group of transcription factors.To date, no studies have been reported regarding the C3H gene family in melon (Cucumis melo). In this study, 38 CmC3H genes were identified in the melon genome, and these genes are unevenly distributed across the 12 chromosomes. Phylogenetic analysis classified the C3H family members into four groups, with significant differences observed in sequence, protein motifs, and gene structure among CmC3H genes within the same group. The CmC3H family contains one pair of segmentally duplicated genes and shares 20, 7, 39, and 38 pairs of homologous C3H genes with Arabidopsis thaliana, rice (Oryza sativa), cucumber (Cucumis sativus), and watermelon (Citrullus lanatus), respectively. Promoter region analysis revealed a high abundance of cis-elements associated with growth and development, hormone regulation, and stress responses. Expression profiling revealed that CmC3H family members exhibit significant tissue-specific expression patterns. Quantitative PCR analysis indicated that six genes (CmC3H4, CmC3H7, CmC3H13, CmC3H24, CmC3H33, and CmC3H38) may play roles in melon's drought stress resistance. Heavy metal lead stress appears to suppress the expression of CmC3H genes. The genes CmC3H24 and CmC3H33 may be involved in regulating melon's resistance to Fusarium wilt infection. CmC3H11 and CmC3H21 can be considered as the key candidate genes for improving the melon's ability to resist both biotic and abiotic stresses.This study provides preliminary insights into the expression profiles of CmC3H genes under drought stress, heavy metal lead stress, and Fusarium wilt infection, offering a theoretical foundation for the molecular mechanisms underlying melon improvement and stress resistance.

Identification and characterization of the elusive protein backbone of the immuno-dominant and species-specific Em2(G11) metacestode antigen of Echinococcus multilocularis

Alveolar echinococcosis (AE) caused by Echinococcus multilocularis, is a severe zoonotic disease in humans. One of the major metacestode antigens of E. multilocularis is the Em2 or Em2(G11) native purified antigen. The Em2 antigen is used for the serological and histopathological diagnosis of AE in humans and plays an important role in parasite-host interactions. As the Em2(G11) antigen is a mucin-type and glycosylated protein, the protein backbone has not been identified yet. We have targeted the protein backbone identification through mass spectrometry (LC-MS/MS) analysis of the Em2(G11) antigen. As a result, we evidenced that the Em2(G11) antigen consists of 33 unique protein candidates of which the most abundant was ''EmuJ_001105600.1''. This protein (889 amino acids) had 427 predicted glycosylation sites. Amino acid composition comparison was in agreement with earlier studies and further confirmed the candidate of interest as the most likely Em2(G11) protein backbone. NCBI BLAST revealed no other known protein homologues in related Echinococcus species nor helminths. After successfully producing this protein recombinantly (Em2rec), a monoclonal antibody (mAbEm2rec) was raised against it. Immunohistochemical stainings of liver tissue sections of AE patients showed that the mAbEm2rec reacts specifically with E. multilocularis antigens solely after deglycosylation with an O-glycosidase cocktail. Similarly, in ELISA, the mAbEm2rec recognized the recombinant and native antigens of E. multilocularis after deglycosylation. These results reveal the nature of this highly glycosylated and specific protein, where mucins are covering the proteomic backbone. For antibody detection in human patients, the native Em2(G11) antigen was superior compared to the Em2rec antigen, indicating the importance of glycosylated epitopes in this immuno-dominant antigen. Of note is the second most abundant protein in the Em2(G11) antigen, namely phosphoenolpyruvate carboxykinase (PEPCK; EmuJ_000292700.1). PEPCK is known to play an important part in the metabolic pathway of gluconeogenesis in E. multilocularis. However, whether this co-eluted protein has any functional importance in the parasite-host interplay of nutrients, growth, and diagnostic significance, is not explored. By combining various approaches, we were able to uncover and confirm the protein backbone of the diagnostic Em2(G11) antigen of E. multilocularis.

In Silico Functional Annotation and Structural Characterization of Hypothetical Proteins in Bacillus paralicheniformis and Bacillus subtilis Isolated from Honey

Bacillus species are ubiquitous and survive in competitive microbial communities under adverse environmental conditions. Bacillus paralicheniformis and Bacillus subtilis obtained from honey revealed a significant proportion of proteins within their genomes as uncharacterized hypothetical proteins (HPs). A total of 1007 HP sequences were evaluated, resulting in the successful annotation of 56 HPs by assigning specific functions to them. A systematic in silico approach, integrating a range of bioinformatics tools and databases to annotate functions, characterize physicochemical properties, determine subcellular localization, and study protein-protein interactions, was used. Homology and de novo models were generated for the HPs, coupled with iterative remodeling and molecular dynamics (MD) simulations. HPs having significant roles in sporulation, biofilm formation, motility, ion transportation, regulation of metabolic processes, DNA repair, replication, and transcription were identified. Classical MD simulations of globular and transducer membrane proteins, along with postprocessing analyses, refined our structural predictions and provided deeper insights into the stability and functional dynamics of the protein structures under physiological conditions. Moreover, we observed a correlation between the percentage of α helix, β sheet, and coil structures in globular proteins and transducer membrane proteins. The integration of iterative loop modeling, MD simulations, and Dictionary of Secondary Structure in Proteins analysis further validated our predicted models and facilitated the identification of regions critical for protein function, thereby enhancing the overall reliability and robustness of our functional annotations. Furthermore, annotation of these hypothetical proteins aids in identifying novel proteins within bacterial cells, ultimately contributing to a deeper understanding of bacterial cell biology and their use for biotechnological purposes.

Identification and characterisation of a novel EhOrc1/Cdc6 from the human pathogen Entamoeba histolytica: an in silico approach

The onset of a pre-replication complex on origin commences DNA replication. The Origin recognition complex (Orc), Cell division cycle protein 6 (Cdc6), and the minichromosome maintenance (Mcm) replicative helicase, along with Chromatin licensing and DNA replication factor 1 (Cdt1), make up the pre-replication complex in eukaryotes. Eukaryotic Orc is made up of six subunits, designated Orc1-6 while monomeric Cdc6 has sequence similarity with Orc1. However, Orc has remained unexplored in the protozoan parasite Entamoeba histolytica. Here we report a single functional Orc1/Cdc6 protein in E. histolytica. Its structural and functional aspects have been highlighted by a detailed in silico analysis that reflects physicochemical characteristics, predictive 3D structure modelling, protein-protein interaction studies, molecular docking and simulation. This in silico study provides insight into EhOrc1/Cdc6 and points out that E. histolytica carries pre-replication machinery that is less complex than higher eukaryotes and closer to archaea. Additionally, it lays the groundwork for future investigations into the methods of origin recognition, and anomalies of cell cycle observed in this enigmatic parasite.Communicated by Ramaswamy H. Sarma.

Whole-Genome Identification of the Flax Fatty Acid Desaturase Gene Family and Functional Analysis of the LuFAD2.1 Gene Under Cold Stress Conditions

Fatty acid desaturase (FAD) is essential for plant growth and development and plant defence response. Although flax (Linum usitatissimum L.) is an important oil and fibre crop, but its FAD gene remains understudied. This study identified 43 LuFAD genes in the flax genome. The phylogenetic analysis divided the FAD genes into seven subfamilies. LuFAD is unevenly distributed on 15 chromosomes, and fragment duplication is the only driving force for the amplification of the LuFAD gene family. In the LuFAD gene promoter region, most elements respond to plant hormones (MeJA, ABA) and abiotic stresses (anaerobic and low temperature). The expression pattern analysis showed that the temporal and spatial expression patterns of all LuFAD genes in different tissues and the response patterns to abiotic stresses (heat and salt) were identified. Subcellular localisation showed that all LuFAD2-GFP were expressed in the endoplasmic reticulum membrane. RT-qPCR analysis revealed that LuFAD2 was significantly upregulated under cold, salt and drought stress, and its overexpression in Arabidopsis thaliana enhanced cold tolerance genes and reduced ROS accumulation. This study offers key insights into the FAD gene family's role in flax development and stress adaptation.

Uncovering the structural stability of Magnaporthe oryzae effectors: a secretome-wide in silico analysis

Rice blast, caused by the ascomycete fungus Magnaporthe oryzae, is a deadly disease and a major threat to global food security. The pathogen secretes small proteinaceous effectors, virulence factors, inside the host to manipulate and perturb the host immune system, allowing the pathogen to colonize and establish a successful infection. While the molecular functions of several effectors are characterized, very little is known about the structural stability of these effectors. We analyzed a total of 554 small secretory proteins (SSPs) from the M. oryzae secretome to decipher key features of intrinsic disorder (ID) and the structural dynamics of the selected putative effectors through thorough and systematic in silico studies. Our results suggest that out of the total SSPs, 66% were predicted as effector proteins, released either into the apoplast or cytoplasm of the host cell. Of these, 68% were found to be intrinsically disordered effector proteins (IDEPs). Among the six distinct classes of disordered effectors, we observed peculiar relationships between the localization of several effectors in the apoplast or cytoplasm and the degree of disorder. We determined the degree of structural disorder and its impact on protein foldability across all the putative small secretory effector proteins from the blast pathogen, further validated by molecular dynamics simulation studies. This study provides definite clues toward unraveling the mystery behind the importance of structural distortions in effectors and their impact on plant-pathogen interactions. The study of these dynamical segments may help identify new effectors as well.Communicated by Ramaswamy H. Sarma.

Advancements in the conservation of the conformational epitope of membrane protein immunogens

Generating antibodies targeting native membrane proteins presents various challenges because these proteins are often embedded in the lipid bilayer, possess various extracellular and intracellular domains, and undergo post-translational modifications. These properties of MPs make it challenging to preserve their stable native conformations for immunization or antibody generation outside of the membranes. In addition, MPs are often hydrophobic due to their membrane-spanning regions, making them difficult to solubilize and purify in their native form. Therefore, employing purified MPs for immunogen preparation may result in denaturation or the loss of native structure, rendering them inadequate for producing antibodies recognizing native conformations. Despite these obstacles, various new approaches have emerged to address these problems. We outline recent advancements in designing and preparing immunogens to produce antibodies targeting MPs. Strategies outlined here are relevant for producing antibodies for research, diagnostics, and therapies and designing immunogens for vaccination purposes.

Development of a multi-epitope vaccine candidate to combat SARS-CoV-2 and dengue virus co-infection through an immunoinformatic approach

Background

Although the SARS-CoV-2 and dengue viruses seriously endanger human health, there is presently no vaccine that can stop a person from contracting both viruses at the same time. In this study, four antigens from SARS-CoV-2 and dengue virus were tested for immunogenicity, antigenicity, allergenicity, and toxicity and chosen to predict dominant T- and B-cell epitopes.

Methods

For designing a multi-epitope vaccine, the sequences were retrieved, and using bioinformatics and immunoinformatics, the physicochemical and immunological properties, as well as secondary structures, of the vaccine were predicted and studied. Additionally, the three-dimensional structure was estimated, improved upon, and confirmed using bioinformatics methods before being docked with TLR-2 and TLR-4. Eight helper T-cell lymphocyte (HTL) epitopes, ten cytotoxic T-cell lymphocyte (CTL) epitopes, nine B-cell epitopes, and TLR agonists were used to create a new multi-epitope vaccine. Furthermore, according to the immunological stimulation hypothesis, the vaccine could stimulate T and B cells to create large quantities of Th1 cytokines and antibodies.

Results

The study indicates that the developed vaccine is a favorable vaccine candidate with antigenicity, immunogenicity, non-toxicity, and non-allergenicity properties. The vaccine construct was made up of 460 amino acids, had an MW of 49391.51 Da, a theoretical pI of 9.86, and the formula C2203H3433N643O618S18, a lipid index of 39.84, a GRAVY of -0.473, an aliphatic index of 63.80, and an instability index of 39.84, which classifies the protein to be stable.

Conclusion

The acquired data showed that both vaccine designs had a considerable chance of preventing the co-infection of SARS-CoV-2 and dengue virus and that they demonstrate good results following in-silico testing. Furthermore, the vaccine may be an effective strategy in preventing SARS-CoV-2 and dengue since it can cause noticeably high levels of Th1 cytokines and antibodies.

Evolutionary and functional analysis of the DIR gene family in Moso bamboo: insights into rapid shoot growth and stress responses

Dirigent (DIR) proteins are key regulators of lignin and lignan biosynthesis and play critical roles in plant hormone responses, abiotic stress tolerance, and growth and development. This study identified and characterized 47 PeDIR genes in Moso bamboo, classifying them into three groups. Phylogenetic and comparative analyses revealed strong evolutionary conservation, with the Moso bamboo PeDIR genes being most closely related to those in rice and maize. DIR proteins within each subfamily exhibited high conservation in motif composition, domain structure, and 3D configuration. Subcellular localization and protein interaction studies further elucidated PeDIR gene functions. Specifically, PeDIR02 primarily localized to the cell membrane and was shown to be unable to form homodimers in yeast two-hybrid (Y2H) assays. Transcriptome and expression analyses revealed the involvement of PeDIR genes in rapid shoot growth, indicating roles in lignin biosynthesis and cell wall modification. Transcriptome and qRT-PCR data also demonstrated the responsiveness of these genes to hormones and abiotic stresses, such as drought and salinity. This study constructed the first comprehensive regulatory network between transcription factors (TFs) and PeDIR genes, identifying ERF, DOF, and MYB TFs as key synergistic regulators of PeDIR gene expression.

Global Gac/Rsm regulatory system activates the biosynthesis of mupirocin by controlling the MupR/I quorum sensing system in Pseudomonas sp. NCIMB 10586

The biosynthesis of mupirocin, a clinically significant antibiotic produced by Pseudomonas sp. NCIMB 10586, is activated by the N-acyl homoserine lactone (AHL) MupR/I quorum sensing (QS) system. However, to date, limited research has focused on the influence of global regulators such as the GacS/A two-component system (TCS) on the MupR/I QS system or mupirocin biosynthesis. In this study, we characterized the regulatory components of the Gac/Rsm transduction system in the mupirocin-producing model strain NCIMB 10586 and investigated their interconnection with the MupR/I QS circuit and subsequent mupirocin biosynthesis. The production of mupirocin was hampered by either gacS inactivation, gacA inactivation, or the double-mutant of the sRNAs ( RsmY and RsmZ). Similarly, the expressions of mupR and mupI, and AHL synthesis significantly decreased in gacS, gacA, or rsmY/Z mutants, indicating that the GacS/A system stimulates mupirocin biosynthesis via the MupR/I QS system. Five CsrA family proteins, RsmA/E/I/F/N, were found in strain NCIMB 10586, and the single and multiple mutants of rsmA/E/I/F/N showed different phenotypes with respect to mupirocin production. Our results revealed that mupirocin biosynthesis was likely to be negatively regulated by RsmA/E/I, but positively regulated by RsmF. Additionally, the RsmF protein was shown to interact with the 5' leader of mupR mRNA. In summary, the Gac/Rsm system positively regulates the biosynthesis of mupirocin mainly through the MupR/I QS system, and the model of the regulatory mechanism is proposed. The elucidation of the Gac/Rsm-MupR/I regulatory pathway could help devise ways for improving mupirocin production through genetic engineering.IMPORTANCEThe Gac/Rsm regulatory system plays a global regulatory role in bacterial physiology and metabolism, including secondary metabolism. Mupirocin is a clinically important antibiotic, produced by Pseudomonas sp. NCIMB 10586, whose biosynthesis is activated by the MupR/I quorum sensing system. Global regulators have important impacts on the gene expression of secondary metabolic gene clusters and QS genes, and the GacS/A two-component system is one of the main regulators across Pseudomonas species, which significantly influences antibiotic production. Our study presented that the expressions of QS genes and mup gene cluster were downregulated in gacS, gacA, or rsmY/Z mutants compared to the wild-type. The inactivation of rsmA/E/I/F/N in NCIMB 10586, encoding CsrA family proteins, showed different regulatory traits of mupirocin production, in which the RsmF protein could interact with the 5' UTR region of mupR mRNA. These findings provide the understanding of the regulatory role of Gac/Rsm on mupirocin biosynthesis and mupR/I QS system and lay foundations for further improving mupirocin production.

Frataxin Traps Low Abundance Quaternary Structure to Stimulate Human Fe-S Cluster Biosynthesis

Iron-sulfur clusters are essential protein cofactors synthesized in human mitochondria by an NFS1-ISD11-ACP-ISCU2-FXN assembly complex. Surprisingly, researchers have discovered three distinct quaternary structures for cysteine desulfurase subcomplexes, which display similar interactions between NFS1-ISD11-ACP protomeric units but dramatically different dimeric interfaces between the protomers. Although the role of these different architectures is unclear, possible functions include regulating activity and promoting the biosynthesis of distinct sulfur-containing biomolecules. Here, crystallography, native ion-mobility mass spectrometry, and chromatography methods reveal the Fe-S assembly subcomplex exists as an equilibrium mixture of these different quaternary structures. Isotope labeling and native mass spectrometry experiments show that the NFS1-ISD11-ACP complexes disassemble into protomers, which can then undergo exchange reactions and dimerize to reform native complexes. Single crystals isolated in distinct architectures have the same activity profile and activation by the Friedreich's ataxia (FRDA) protein frataxin (FXN) when rinsed and dissolved in assay buffer. These results suggest FXN functions as a "molecular lock" and shifts the equilibrium toward one of the architectures to stimulate the cysteine desulfurase activity and promote iron-sulfur cluster biosynthesis. An NFS1-designed variant similarly shifts the equilibrium and partially replaces FXN in activating the complex. We propose that eukaryotic cysteine desulfurases are unusual members of the morpheein class of enzymes that control their activity through their oligomeric state. Overall, the findings support architectural switching as a regulatory mechanism linked to FXN activation of the human Fe-S cluster biosynthetic complex and provide new opportunities for therapeutic interventions of the fatal neurodegenerative disease FRDA.

Evaluation of Synthetic Peptides from Schistosoma mansoni ATP Diphosphohydrolase 1: In Silico Approaches for Characterization and Prospective Application in Diagnosis of Schistosomiasis

Schistosomiasis is the infection caused by Schistosoma mansoni and constitutes a worldwide public health problem. The parasitological recommended method and serological methods can be used for the detection of eggs and antibodies, respectively. However, both have limitations, especially in low endemicity areas. Thus, new approaches for the diagnosis of schistosomiasis are essential. In this study, a six-amino acid peptide and derived sequences from SmATPDase1 were synthesized for the evaluation of immunogenicity. SmATPDase1 is included in a protein group in S. mansoni tegument; therefore, its peptides could be potential candidates for diagnostic antigens. In the hypothetical SmATPDase1 three-dimensional structure, peptides are located in a region exposed and accessible to antibody binding. In addition, peptide amino acid sequences are conserved in the most relevant Schistosoma species and have low identity with human NTPDases isoforms. Swiss mice immunization resulted in significant anti-peptide polyclonal antibodies production, which recognized a 63 kDa protein in tegument and adult worm preparations. By immunofluorescence microscopy, polyclonal antibodies also identified this enzyme in cercariae. Sera of infected animals presented high seropositivity in ELISA-peptides, with an area under curve (AUC) greater than 0.96 for all peptides. In mice with low parasite burden, we observed a seropositivity AUC > 0.9. Reactivity in the prepatent period exhibited AUC values greater than 0.94 for all peptides. Anti-P1425 monoclonal antibodies were successfully produced, and mAbs recognized the integral protein in ELISA and Western blots. The data indicate that peptides from SmATPDase1 are potential biomarkers for schistosomiasis, and anti-peptide antibodies are interesting tools for the detection of the infection.

Identification and expression analysis of the RBOH gene family of Isatis indigotica Fort. and the potential regulation mechanism of RBOH gene on H2O2 under salt stress

KEY MESSAGE

RBOH gene may regulate the resistance of Isatis indigotica Fort. to salt stress by mediating the production of H2O2. RBOH gene plays an important role in plant growth and development, abiotic and biotic stress response, and hormone signalling. However, studies on RBOH gene expression and molecular mechanism of Isatis indigotica Fort. under salt stress have not been reported. This study identified 10 genes of the I. indigotica RBOH gene family (IiRBOH) and divided them into five subfamilies (I-V). Genes within the same class show conserved structural features and similar amino acid sequences. Analysis of CRE suggested that IiRBOH genes might play roles in growth and development, metabolism, hormone regulation, and stress response. Two physiological indicators of I. indigotica treated with salt for different days were detected. It was found that the content of H2O2 in the I. indigotica tissue first increased, then decreased and increased again. The catalase activity also showed a trend of first increasing and then decreasing. The qRT-PCR results showed that these IiRBOH genes showed different expression patterns in response to salt stress, and some of these genes may be involved in the resistance of I. indigotica to salt stress. Through RT-PCR analysis and screening on the PlantCARE website, it was found that IiRBOHA and IiRBOHC not only possess W-box CRE but also exhibit high expression under salt stress. Y1H experiments were conducted with the WRKY genes predicted by phylogenetic analysis to regulate salt stress potentially, and it was discovered that IiWRKY6 and IiWRKY54 can directly activate the transcription of the IiRBOHA gene promoter. This study preliminarily explored the mechanism by which the RBOH gene in I. indigotica mediates H2O2 to resist salt stress, thus laying a foundation for further research on the biological functions of the RBOH gene in I. indigotica.

Sam-Sam Association Between EphA2 and SASH1: In Silico Studies of Cancer-Linked Mutations

Recently, SASH1 has emerged as a novel protein interactor of a few Eph tyrosine kinase receptors like EphA2. These interactions involve the first N-terminal Sam (sterile alpha motif) domain of SASH1 (SASH1-Sam1) and the Sam domain of Eph receptors. Currently, the functional meaning of the SASH1-Sam1/EphA2-Sam complex is unknown, but EphA2 is a well-established and crucial player in cancer onset and progression. Thus, herein, to investigate a possible correlation between the formation of the SASH1-Sam1/EphA2-Sam complex and EphA2 activity in cancer, cancer-linked mutations in SASH1-Sam1 were deeply analyzed. Our research plan relied first on searching the COSMIC database for cancer-related SASH1 variants carrying missense mutations in the Sam1 domain and then, through a variety of bioinformatic tools and molecular dynamic simulations, studying how these mutations could affect the stability of SASH1-Sam1 alone, leading eventually to a defective fold. Next, through docking studies, with the support of AlphaFold2 structure predictions, we investigated if/how mutations in SASH1-Sam1 could affect binding to EphA2-Sam. Our study, apart from presenting a solid multistep research protocol to analyze structural consequences related to cancer-associated protein variants with the support of cutting-edge artificial intelligence tools, suggests a few mutations that could more likely modulate the interaction between SASH1-Sam1 and EphA2-Sam.

Design of a Multi-Epitope Vaccine against the Glycoproteins of Newcastle Disease Virus by Using an Immunoinformatics Approach

Newcastle disease (ND) causes major economic losses in poultry farming in Madagascar and many other countries. Previously, vaccines based on attenuated or inactivated Newcastle disease viruses (NDV) have been effective against this disease. However, their efficacy has declined due to viral mutations over time. To address this, two new multi-epitope vaccines (MEV) have been designed using immunoinformatics methods. First, 26 conserved epitopes from the fusion protein and 22 from the hemagglutinin-neuraminidase protein of 12 NDV strains isolated in Madagascar were selected to design the MEV. These epitopes were fused with specific linkers. Additionally, the adjuvant Avian Beta-Defensins-1 and the 6xHis tag were added to the N- and C-terminal ends of the vaccine formulations, respectively. The antigenicity, allergenicity, solubility, and physicochemical properties of the designed MEV were evaluated. Their three-dimensional structures were also modeled. Molecular docking studies and dynamic simulations were then conducted with the chicken Toll-like receptor 7 (TLR7) to assess the binding affinity of the MEV with this receptor. Finally, an immunological simulation was carried out to assess the ability of the candidate vaccine to induce an effective immune response. Through immunoinformatics analysis, both MEVs developed in this study were found to be highly antigenic, nonallergenic, and physicochemically stable. In addition, they showed significant interaction with the TLR7 receptor. They also have the capacity to trigger immune responses and promote the formation of memory cells following immunization. Therefore, these vaccines represent promising candidates for the control of ND. As this is an immunoinformatics study based on in silico methods, both in vitro and in vivo experiments are required to confirm and extend these results.

Crystallographic mining driven computer-guided approach to identify the ASK1 inhibitor likely to perturb the catalytic region

The pathological levels of reactive oxygen species (ROS) and oxidative stress has been recognized as a critical driver for inflammatory disorders. Apoptosis signal-regulating kinase 1 (ASK1) has been reported to be activated by intracellular ROS and its inhibition leads to a down regulation of p38-and JNK-dependent signaling. ASK1 inhibitors are reported to have the potential to treat clinically important inflammatory pathologies including liver, pulmonary and renal disorders. In view of its biological and pathological significance, inhibition of ASK1 with small molecules has been pursued as an attractive strategy to combat human diseases such as non-alcoholic steatohepatitis (NASH). Despite several ASK1 inhibitors being developed, the failure in Phase 3 clinical trials of most advanced candidate selonsertib's, underscores to discover therapeutic agents with diverse chemical moiety. Here, by using structural pharmacophore and enumeration strategy on mining co-crystals of ASK1, different scaffolds were generated to enhance the chemical diversity keeping the critical molecular interaction in the catalytic site intact. A total of 15,772 compounds were generated from diverse chemical scaffolds and were evaluated using a virtual screening pipeline. Based on docking and MM-GBSA scores, a lead candidate, S3C-1-D424 was identified from top hits. A comparative molecular dynamics simulations (MD) of APO, Selonsertib and shortlisted potential candidates combined with pharmacokinetics profiling and thermodynamic analysis, demonstrating their suitability as potential ASK1 inhibitors to explore further for establishment towards hit-to-lead campaign.Communicated by Ramaswamy H. Sarma.

A microaerobically induced small heat shock protein contributes to Rhizobium leguminosarum/Pisum sativum symbiosis and interacts with a wide range of bacteroid proteins

During the establishment of the symbiosis with legume plants, rhizobia are exposed to hostile physical and chemical microenvironments to which adaptations are required. Stress response proteins including small heat shock proteins (sHSPs) were previously shown to be differentially regulated in bacteroids induced by Rhizobium leguminosarum bv. viciae UPM791 in different hosts. In this work, we undertook a functional analysis of the host-dependent sHSP RLV_1399. A rlv_1399-deleted mutant strain was impaired in the symbiotic performance with peas but not with lentil plants. Expression of rlv_1399 gene was induced under microaerobic conditions in a FnrN-dependent manner consistent with the presence of an anaerobox in its regulatory region. Overexpression of this sHSP improves the viability of bacterial cultures following exposure to hydrogen peroxide and to cationic nodule-specific cysteine-rich (NCR) antimicrobial peptides. Co-purification experiments have identified proteins related to nitrogenase synthesis, stress response, carbon and nitrogen metabolism, and to other relevant cellular functions as potential substrates for RLV_1399 in pea bacteroids. These results, along with the presence of unusually high number of copies of shsp genes in rhizobial genomes, indicate that sHSPs might play a relevant role in the adaptation of the bacteria against stress conditions inside their host.IMPORTANCEThe identification and analysis of the mechanisms involved in host-dependent bacterial stress response is important to develop optimal Rhizobium/legume combinations to maximize nitrogen fixation for inoculant development and might have also applications to extend nitrogen fixation to other crops. The data presented in this work indicate that sHSP RLV_1399 is part of the bacterial stress response to face specific stress conditions offered by each legume host. The identification of a wide diversity of sHSP potential targets reveals the potential of this protein to protect essential bacteroid functions. The finding that nitrogenase is the most abundant RLV_1399 substrate suggests that this protein is required to obtain an optimal nitrogen-fixing symbiosis.

Genome-wide identification and comparative analysis of the AP2/ERF gene family in Prunus dulcis and Prunus tenella: expression of PdAP2/ERF genes under freezing stress during dormancy

The AP2/ERF (APETALA2/ethylene responsive factor) transcription factor family, one of the largest in plants, plays a crucial role in regulating various biological processes, including plant growth and development, hormone signaling, and stress response. This study identified 114 and 116 AP2/ERF genes in the genomes of 'Wanfeng' almond (Prunus dulcis) and 'Yumin' wild dwarf almond (Prunus tenella), respectively. These genes were categorized into five subfamilies: AP2, DREB, ERF, RAV, and Soloist. The PdAP2/ERF and PtAP2/ERF members both demonstrated high conservation in protein motifs and gene structures. Members of both families were unevenly distributed across eight chromosomes, with 30 and 27 pairs of segmental duplications and 15 and 18 pairs of tandem repeated genes, respectively. The promoter regions of PdAP2/ERF and PtAP2/ERF family members contained numerous important cis-elements related to growth and development, hormone regulation, and stress response. Expression pattern analysis revealed that PdAP2/ERF family members exhibited responsive characteristics under freezing stress at different temperatures in perennial dormant branches. Quantitative fluorescence analysis indicated that PdAP2/ERF genes might be more intensely expressed in the phloem of perennial dormant branches of almond, with the opposite trend observed in the xylem. This study compared the characteristics of PdAP2/ERF and PtAP2/ERF gene family members and initially explored the expression patterns of PdAP2/ERF genes in the phloem and xylem of perennial dormant branches. The findings provide a theoretical foundation for future research on almond improvement and breeding, as well as the molecular mechanisms underlying resistance to freezing stress.

In silico design of a multi-epitope vaccine against Mycobacterium avium subspecies paratuberculosis

The widespread chronic enteritis known as Paratuberculosis (PTB) or Johne's disease (JD) is caused by Mycobacterium avium subspecies paratuberculosis (MAP), posing a significant threat to global public health. Given the challenges associated with PTB or JD, the development and application of vaccines are potentially important for disease control. The aim of this study was to design a multi-epitope vaccine against MAP. A total of 198 MAP genomes were analyzed using pan-genome and reverse vaccinology approaches. B-cell and T-cell epitope analysis was performed on the selected promising cross-protective antigens followed by selection of epitopes with high antigenicity, no allergenicity, and no toxicity for the design of the vaccine. The designed vaccine was evaluated through molecular dynamics simulations, molecular docking, and immunological simulations. The results revealed the identification of five promising cross-protective antigens. In total, 10 B-cell epitopes, 10 HTL epitopes, and 9 CTL epitopes were selected for the design of the vaccine. Both the vaccine candidate and the vaccine-TLR4 complex demonstrated considerable stability in molecular dynamics simulations. Molecular docking studies confirmed that the vaccine candidate successfully interacted with TLR4. Immunological simulations showed an increase in both B-cell and T-cell populations after vaccination. Additionally, the vaccine candidate exhibited a codon adaptability index of 1.0 and a GC content of 53.64%, indicating strong potential for successful expression in Escherichia coli. This research developed a multi-epitope vaccine targeting MAP through pan-genomes and reverse vaccinology methods, offering innovative strategies for creating effective vaccines against MAP.

In silico design of an epitope-based vaccine ensemble for fasliolopsiasis

Introduction

Fasciolopsiasis, a food-borne intestinal disease is most common in Asia and the Indian subcontinent. Pigs are the reservoir host, and fasciolopsiasis is most widespread in locations where pigs are reared and aquatic plants are widely consumed. Human infection has been most commonly documented in China, Bangladesh, Southeast Asia, and parts of India. It predominates in school-age children, and significant worm burdens are not uncommon. The causal organism is Fasciolopsis buski, a giant intestinal fluke that infects humans and causes diarrhoea, fever, ascites, and intestinal blockage. The increasing prevalence of medication resistance and the necessity for an effective vaccination make controlling these diseases challenging.

Methods

Over the last decade, we have achieved major advances in our understanding of intestinal fluke biology by in-depth interrogation and analysis of evolving F. buski omics datasets. The creation of large omics datasets for F. buski by our group has accelerated the discovery of key molecules involved in intestinal fluke biology, toxicity, and virulence that can be targeted for vaccine development. Finding successful vaccination antigen combinations from these huge number of genes/proteins in the available omics datasets is the key in combating these neglected tropical diseases. In the present study, we developed an in silico workflow to select antigens for composing a chimeric vaccine, which could be a significant technique for developing a fasciolopsiasis vaccine that prevents the parasite from causing serious harm.

Results and discussion

This chimeric vaccine can now be tested experimentally and compared to other vaccine candidates to determine its potential influence on human health. Although the results are encouraging, additional validation is needed both in vivo and in vitro. Considering the extensive genetic data available for intestinal flukes that has expanded with technological advancements, we may need to reassess our methods and suggest a more sophisticated technique in the future for identifying vaccine molecules.

Comparative analysis of the WRKY gene family between Chimonanthus praecox and C. salicifolius

Gene duplications provide evolutionary potentials for generating novel functions. Chimonanthus praecox and C. salicifolius are closely related species from Calycantaceae, Magnoliids. In this study, we compared the WRKY gene family from C. praecox and C. salicifolius, and predicted the potential gene function through gene expression patterns to explore the evolution of orthologous and paralogous gene pairs. A total of 73 and 85 WRKY genes were identified and analyzed from the whole genome sequencing of C. praecox and C. salicifolius. Based on the phylogenetic analysis, CpWRKY and CsWRKY genes were clustered into three groups (Group I、II、III) and 5 subgroups (Group IIa、IIb、IIc、IId、IIe). In C. praecox and C. salicifolius, we identified thirty-six and fifty-four pairs of WRKY segmental duplicated genes, respectively, along with two and three pairs of tandem duplicates, indicating that segmental duplication plays a crucial role in the evolution of Chimonanthus WRKY gene family. Most WRKY duplication gene pairs originated from segmental duplications before the first whole genome duplication (WGD), highlighting this period as a significant source of genetic diversity and functionality for the WRKY family. The analysis of WRKY gene expression levels suggests that CsWRKY18 and CsWRKY68 may promote the growth of the roots in C. salicifolius. Comparisons of expression profiles between species revealed that five orthologous gene pairs presented identical expression trends, indicating functional conservation and absence of neo-functionalization or sub-functionalization. However, most orthologous gene pairs exhibit differences in expression patterns, suggesting that they have undergone functional divergence. This functional differentiation may be due to the different selective pressures faced by C. praecox and C. salicifolius during their speciation processes. This study provided detailed information on the WRKY gene family from C. praecox and C. salicifolius, and a new insight for studying gene duplication and function evolution.

Innovative epitopes in Staphylococcal Protein-A an immuno-informatics approach to combat MDR-MRSA infections

Background

Methicillin-resistant Staphylococcus aureus (MRSA) poses a significant challenge in clinical environments due to its resistance to standard antibiotics. Staphylococcal Protein A (SpA), a crucial virulence factor of MRSA, undermines host immune responses, making it an attractive target for vaccine development. This study aimed to identify potential epitopes within SpA that could elicit robust immune responses, ultimately contributing to the combat against multidrug-resistant (MDR) MRSA.

Methods

The SpA protein sequence was retrieved from the UniProt database, with various bioinformatics tools employed for epitope prediction. T-cell epitopes were identified using the Tepitool server, focusing on high-affinity interactions with prevalent human leukocyte antigens (HLAs). B-cell epitopes were predicted using the BepiPred tool. Predicted epitopes underwent docking studies with HLA molecules to evaluate binding properties. In-silico analyses confirmed the antigenicity, promiscuity, and non-glycosylated nature of the selected epitopes.

Results

Several T and B cell epitopes within SpA were identified, showcasing high binding affinities and extensive population coverage. A multi-epitope vaccine construct, linked by synthetic linkers and an adjuvant, was modelled, refined, and validated through various bioinformatics assessments. The vaccine candidate was subsequently docked with Toll-like receptor 4 (TLR-4) to evaluate its potential for immunogenicity.

Conclusion

This study lays the groundwork for developing epitope-based vaccines targeting SpA in MRSA, identifying promising candidates for experimental validation and contributing to innovative immunotherapeutic strategies against MRSA infections.

Genome-wide identification of CAMTA gene family in teak (Tectona grandis) and functional characterization of TgCAMTA1 and TgCAMTA3 in cold tolerance

BACKGROUND

Calmodulin-binding transcription activator (CAMTA) proteins play significant roles in signal transduction, growth and development, as well as abiotic stress responses, in plants. Understanding their involvement in the low-temperature stress response of teak is vital for revealing cold resistance mechanisms.

RESULTS

Through bioinformatics analysis, the CAMTA gene family in teak was examined, and six CAMTA genes were identified in teak. The encoded proteins were predicted to be located in the nucleus and exhibited hydrophilic properties, with molecular weights ranging from 103.4 to 123.3 kDa and isoelectric points ranging from 5.49 to 7.55. On the basis of protein sequence homology, the CAMTA family could be divided into three subgroups. Domain and 3D structure analyses demonstrated that all the TgCAMTA proteins contained the typical CAMTA domain with the CaMBD binding domain, which was exposed on the surface. Expression analysis of different tissues revealed the expression of TgCAMTA genes in teak roots, stems, leaves, flowers, fruits, and branches. Furthermore, the promoter region contained various cis-acting elements related to light, hormone, and abiotic stress responses. After low-temperature stress treatment, different expression patterns of TgCAMTAs were observed in teak roots, stems, and leaves, with TgCAMTA1 showing the highest expression level in leaves compared with stems. Transgenic lines carrying the TgCAMTA1/3 promoter::GUS construct cold stress induction of TgCAMTA1/3 genes revealed the presence of multiple low-temperature responsive cis-acting elements in the TgCAMTA1/3 promoter region. Subcellular localization analysis indicated that these genes were functional predominantly in the nucleus. Compared with wild-type Arabidopsis, TgCAMTA1/3-overexpressing Arabidopsis presented increased tolerance to freezing stress, with increased expression of AtCOR genes. Moreover, under low-temperature conditions, TgCAMTA3-overexpressing Arabidopsis presented significantly elevated expression levels of genes related to the CBF signaling pathway, including AtCBF1/2/3.

CONCLUSIONS

Our findings add significantly to the existing knowledge regarding cold stress tolerance and help elucidate cold response mechanisms in teak.

Molecular Evolution of Paralogous Cold Shock Proteins in E. coli: A Study of Asymmetric Divergence and Protein Functional Networks

Nine homologous Cold Shock Proteins (Csps) have been recognized in the E.coli Cold Shock Domain gene family. These Csps function as RNA chaperones. This study aims to establish the evolutionary relationships among these genes by identifying and classifying their paralogous counterparts. It focuses on the physicochemical, structural, and functional analysis of the genes to explore the phylogeny of the Csp gene family. Computational tools were employed for protein molecular modeling, conformational analysis, functional studies, and duplication-divergence assessments. The research also examined amino acid conservation, protein mutations, domain-motif patterns, and evolutionary residue communities to better understand residual interactions, evolutionary coupling, and co-evolution. H33, M5, W11 and F53 residues were highly conserved within the protein family. It was further seen that residues M5, G17, G58, G61, P62, A64, V67 were intolerant to any kind of mutation whereas G3, D40, G41, Y42, S44, T54, T68, S69 were most tolerable towards substitutions. The study of residue communities displayed that the strongest residue coupling was observed in N13, F18, S27, F31, and W11. It was observed that all the gene pairs except CspF/CspH had new motifs generated over time. It was ascertained that all the gene pairs underwent asymmetric expression divergence after duplication. The Ka/ Ks ratio also revealed that all residues undertook neutral and purifying selection pressure. New functions were seen to develop in gene pairs evident from generation of new motifs. The discovery of new motifs and functions in Csps highlights their adaptive versatility, crucial for E. coli's resilience to environmental stressors and valuable for understanding bacterial stress response mechanisms. These findings will pave the way for future investigations into Csp evolution, with potential applications in microbial ecology and antimicrobial strategy development.

Investigating novel Streptomyces bacteriophage endolysins as potential antimicrobial agents

As antibiotic resistance has become a major global threat, the World Health Organization (WHO) has urgently called for alternative strategies for control of bacterial infections. Endolysin, a phage-encoded protein, can degrade bacterial peptidoglycan (PG) and disrupt bacterial growth. According to the WHO, there are only three endolysin products currently in clinical phase development. In this study, we explore novel endolysins from Streptomyces phages as only a few of them have been experimentally characterized. Using several bioinformatics tools, we identified nine different functional domain combinations from 250 Streptomyces phages putative endolysins. LazerLemon gp35 (CHAP; LL35lys), Nabi gp26 (amidase; Nb26lys), and Tribute gp42 (PGRP/amidase; Tb42lys) were selected for experimental studies. We hypothesized that (i) the proteins of interest will have the ability to degrade purified PG, and (ii) the proteins will have potential antimicrobial activity against bacteria from families of importance in antibiotic resistance, such as ESKAPE safe relatives (Enterococcus raffinosus, Staphylococcus epidermidis, Klebsiella aerogenes, Acinetobacter baylyi, Pseudomonas putida, and Escherichia coli). LL35lys, Nb26lys, and Tb42lys exhibit PG-degrading activity on zymography and hydrolysis assay. The enzymes (100 µg/mL) can reduce PG turbidity to 32%-40%. The killing assay suggests that Tb42lys has a broader range (E. coli, P. putida, A. baylyi and K. aerogenes). While Nb26lys better attacks Gram-negative than -positive bacteria, LL35lys can only reduce the growth of the Gram-positive ESKAPE strains but does so effectively with a low MIC90 of 2 µg/mL. A higher concentration (≥300 µg/mL) of Nb26lys is needed to inhibit P. putida and K. aerogenes. From 250 putative endolysins, bioinformatic methods were used to select three putative endolysins for cloning and study: LL35lys, Nb26lys, and Tb42lys. All have shown PG-degrading activity, a critical function of endolysin. With a low MIC, LL35lys shows activity for the Gram-positive ESKAPE strains, while Nb26lys and Tb42lys are active against the Gram negatives. Therefore, endolysins from Streptomyces phages have potential as possible antimicrobial agents against ESKAPE bacteria.

IMPORTANCE

As antibiotic resistance has become a major global threat, the World Health Organization (WHO) has urgently called for alternative strategies for control of bacterial infections. Endolysin, a phage-encoded protein, can degrade bacterial peptidoglycan in the bacterial cell wall and disrupt bacterial growth. According to the WHO, there are only three endolysin products currently in clinical phase development. In this study we explored novel endolysins from Streptomyces phages as only a few of them have been experimentally characterized. Using several bioinformatics tools, we identified nine different combinations of functional enzymatic domain types from 250 Streptomyces bacteriophages possible endolysins. From these, three potential endolysins were selected for experimental characterization. All three showed positive results in degrading cell wall material and disrupting bacterial growth, indicating their potential as possible antimicrobial agents.

Identification and Characterization of a Protease Producing Bacillus cereus Strain From Tannery Waste for Efficient Dehairing of Goat Skin

Environmental pollution has been a significant concern for the last few years. The leather industry significantly contributes to the economy but is one of Bangladesh's most prominent polluting industries. It is also responsible for several severe diseases such as cancer, lung diseases, and heart diseases of leather workers because they use bleaching agents and chemicals, and these have numerous adverse effects on human health. The study was aimed at isolating, identifying, and molecularly characterizing bacteria that produce protease enzymes that are highly capable of dehairing goat hide. Several attempts were made to isolate and identify protease-producing bacterial strains from different soil samples of tannery wastes. Initially, a total of four isolates were selected from tannery soil. After the different phenotypic and morphological characterization, one isolate (BS2) showed Gram-positive, rod-shaped, and spore-forming characteristics and could produce novel hair-degrading protease enzymes. The growth profile and protease activity of the bacteria at 37°C were observed; a basal level of extracellular protease increased with time. The enzyme's proteolytic activity was measured, and the unit of enzyme activity of the enzyme sample was 18.1. The ExPASy server (ProtParam) tool was used to estimate the physicochemical characteristics of the proteins and found molecular weight (MW) (7375.94 Da), aliphatic index (71.56), instability index (II, 80.63), Grand Average of Hydropathy (GRAVY) (-0.231), and isoelectric point (11.41). The protein-protein interactions (PPI) networks were generated using the Search Tool for the Retrieval of Interacting Genes (STRING) database and Cytoscape software. The PSIPRED v.4.0 and SAVES v.6.0 programs were used to determine the secondary and three-dimensional assembly of the selected protein. They found alpha helix (16, 25.00%), extended strand (6, 9.38%), beta-turn (5, 7.81%), and random coil (37, 57.81%). DNA isolation and purification were carried out, and the purity ratio was ~2.17 at 260 and 280 nm. Polymerase chain reaction (PCR) for amplifying the 16S rRNA gene was conducted, and the isolate was authentically recognized as Bacillus cereus (BS2) based on morphological, biochemical, and molecular analyses. The quantitative assessment has shown that 40 mL of culture centrifugation could dehair 2 × 1 cm of goat leather sample in 9 h. This potential bacterial strain can be used in the leather industry as an ecofriendly alternative to chemical dehairing, which can reduce environmental pollution and the risk of severe diseases among leather industry workers.

Investigating Protein-Protein Interactions of Autophagy-Involved TNIP1

Myriad proteins are involved in the process of autophagy, which they participate in via their protein-protein interactions (PPI). Herein we outline a methodology for examining such interactions utilizing the case of intrinsically disordered protein (IDP) TNIP1 and its interaction with linear M1-linked polyubiquitin. This includes methods for recombinant production, purification, immuno-identification, and analysis of an IDP associated with autophagy, its ordered binding partner, and means of quantitatively analyzing their interaction.

Design and evaluation of a multi-epitope DNA vaccine against HPV16

Cervical cancer, among the deadliest cancers affecting women globally, primarily arises from persistent infection with high-risk human papillomavirus (HPV). To effectively combat persistent infection and prevent the progression of precancerous lesions into malignancy, a therapeutic HPV vaccine is under development. This study utilized an immunoinformatics approach to predict epitopes of cytotoxic T lymphocytes (CTLs) and helper T lymphocytes (HTLs) using the E6 and E7 oncoproteins of the HPV16 strain as target antigens. Subsequently, through meticulous selection of T-cell epitopes and other necessary elements, a multi-epitope vaccine was constructed, exhibiting good immunogenic, physicochemical, and structural characteristics. Furthermore, in silico simulations showed that the vaccine not only interacted well with toll-like receptors (TLR2/TLR3/TLR4), but also induced a strong innate and adaptive immune response characterized by elevated Th1-type cytokines, such as interferon-gamma (IFN-γ) and interleukin-2 (IL2). Additionally, our study investigated the effects of different immunization intervals on immune responses, aiming to optimize a time-efficient immunization program. In animal model experiments, the vaccine exhibited robust immunogenic, therapeutic, and prophylactic effects. Administered thrice, it consistently induced the expansion of specific CD4 and CD8 T cells, resulting in substantial cytokines release and increased proliferation of memory T cell subsets in splenic cells. Overall, our findings support the potential of this multi-epitope vaccine in combating HPV16 infection and signify its candidacy for future HPV vaccine development.

Design and in silico analysis of a novel peptide-based multiepitope vaccine against glioblastoma multiforme by targeting tumor-associated macrophage

CD204 is a distinct indicator for tumor-associated macrophages (TAMs) in glioma. Evidence indicates that CD204-positive TAMs are involved in the aggressive behavior of various types of cancers. This study was conducted to develop a new and effective peptide-based vaccine for GBM, specifically targeting CD204. Epitopes of the target protein were identified using NetMHCpan 4.1a, NetMHCIIpan-4.0, and ABCpred tools. Subsequently, the predicted epitopes were evaluated using bioinformatics tools to assess their antigenicity, non-allergenicity, immunogenicity, non-toxicity, and potential to stimulate the production of IL-4 and IFN-γ in HTL epitopes. Selected T-cell epitopes demonstrated a robust binding affinity with the particular HLA alleles. Finally, four HTL epitopes, three CTL epitopes, and two B-cell epitopes, jointed via linkers and adjuvant, were used for the final vaccine construct design. Analysis disclosed that the developed vaccine demonstrated robust antigenic properties while proving soluble, stable, non-toxic, and non-allergenic. Additionally, molecular docking studies and molecular dynamics simulations confirmed a robust correlation between the designed vaccine and TLR-2 and TLR-4 immune receptors. The molecular docking results demonstrated a strong interaction between the newly developed vaccine and TLR2 (-895.1 kcal/mol) and TLR4 (-881.0 kcal/mol) receptors. During the simulation, the vaccine-TLR2 and vaccine-TLR4 complexes exhibited binding energies of -113.41 and -106.61 kcal/mol, respectively. Analysis by different bioinformatic tools indicated the potential of the designed vaccine in immune stimulation and a significant elevation in IgG and IgM antibodies, T-helper cells, T-cytotoxic cells, INF-γ, IL-2, and IL-4. Research findings show that the newly designed multi-epitope vaccine is promising in providing long-term immunity against GBM and offers a promising therapeutic alternative.

Bioinformatics analysis and alternative polyadenylation in Heat Shock Proteins 70 (HSP70) family members

OBJECTIVE

The Heat Shock Protein 70 (HSP70) family is a highly conserved group of molecular chaperones essential for maintaining cellular homeostasis. These proteins are necessary for protein folding, assembly, and degradation and involve cell recovery from stress conditions. HSP70 proteins are upregulated in response to heat shock, oxidative stress, and pathogenic infections. Their primary role is preventing protein aggregation, refolding misfolded proteins, and targeted degradation of irreparably damaged proteins. Given their involvement in fundamental cellular processes and stress responses, HSP70 proteins are critical for cell survival and modulating disease outcomes in cancer, neurodegeneration, and other pathologies. The present study aims to understand domain architecture, physicochemical properties, phosphorylation, ubiquitination, and alternative polyadenylation site prediction in various HSP70 members.

METHOD

SMART and InterProScan software were used for domain analysis. EXPASY Protparam, NetPhos 3.1 server DTU, and MUbisiDa were used for physicochemical analysis, phosphorylation, and ubiquitination site analysis, respectively. Alternative polyadenylation was studied using the EST database.

RESULT

Domain analysis shows that coiled-coil and nucleotide-binding domains are present in some of the HSP70 members. Five HSP70 family members have alternate polyadenylation sites in their 3'UTR.

CONCLUSION

The present work has provided valuable insights into their structure, functions, interactome, and polyadenylation patterns. Studying their therapeutic potential in diseases like cancer can be helpful.

StaPep: An Open-Source Toolkit for Structure Prediction, Feature Extraction, and Rational Design of Hydrocarbon-Stapled Peptides

All-hydrocarbon stapled peptides, with their covalent side-chain constraints, provide enhanced proteolytic stability and membrane permeability, making them superior to linear peptides. However, tools for extracting structural and physicochemical descriptors to predict the properties of hydrocarbon-stapled peptides are lacking. To address this, we present StaPep, a Python-based toolkit for generating 3D structures and calculating 21 features for hydrocarbon-stapled peptides. StaPep supports peptides containing two non-standard amino acids (norleucine and 2-aminoisobutyric acid) and six non-natural anchoring residues (S3, S5, S8, R3, R5, and R8), with customization options for other non-standard amino acids. We showcase StaPep's utility through three case studies. The first generates 3D structures of these peptides with a mean RMSD of 1.62 ± 0.86, offering essential structural insights for drug design and biological activity prediction. The second develops machine learning models based on calculated molecular features to differentiate between membrane-permeable and non-permeable stapled peptides, achieving an AUC of 0.93. The third constructs regression models to predict the antimicrobial activity of stapled peptides against Escherichia coli, with a Pearson correlation of 0.84. StaPep's pipeline spans data retrieval, structure generation, feature calculation, and machine learning modeling for hydrocarbon-stapled peptides. The source codes and data set are freely available on Github: https://github.com/dahuilangda/stapep_package.

Genome-Wide Identification and Analysis of BrTCP Transcription Factor Family Genes Involved in Cold Stress Tolerance in Winter Rapeseed (Brassica rapa L.)

TCP transcription factors are important during plant growth and stress responses. However, their role in the cold stress response of Brassica rapa L. remains poorly understood. In this research, we identified the TCPs gene family in B. rapa to learn the features of the BrTCP gene family, functionally annotating the interacting proteins of TCP4 and analyzing their expression levels. Our results illustrated the presence of 19 members of the BrTCPs family in B. rapa, exhibiting molecular weights ranging from 27,367.45 to 59,433.64 Da. All identified proteins were classified as unstable, with isoelectric points ranging from 5.5 to 9.48. Subcellular localization forecasted that TCP proteins were all positioned in the nucleus. The BrTCP gene structure is relatively simple, with only seven members possessing introns, and none of the members contain UTR regions. BrTCPs comprise hormone-, light-, and stress-responsive elements. We found that the frequency of photoresponsive elements was greatest in the promoter region, suggesting that BrTCP genes are regulated by light signals and function synergistically with plant growth and development. In addition, five candidate interaction proteins of BrTCP4 were identified using yeast two-hybrid screening. RNA-Seq and q-PCR analyses of the interacting genes revealed differential expression of BrTCP family genes across various tissues following cold stress. Significant responses were observed under low-temperature stress, drought stress, and rehydration treatment, suggesting that these genes play crucial roles as regulators of the molecular network mechanisms responding to stress. This study enhances our understanding of the BrTCP family and provides significant insights into the stress tolerance mechanisms of B. rapa.

The Identification of AMT Family Genes and Their Expression, Function, and Regulation in Chenopodium quinoa

Quinoa (Chenopodium quinoa) is an Andean allotetraploid pseudocereal crop with higher protein content and balanced amino acid composition in the seeds. Ammonium (NH4+), a direct source of organic nitrogen assimilation, mainly transported by specific transmembrane ammonium transporters (AMTs), plays important roles in the development, yield, and quality of crops. Many AMTs and their functions have been identified in major crops; however, no systematic analyses of AMTs and their regulatory networks, which is important to increase the yield and protein accumulation in the seeds of quinoa, have been performed to date. In this study, the CqAMTs were identified, followed by the quantification of the gene expression, while the regulatory networks were predicted based on weighted gene co-expression network analysis (WGCNA), with the putative transcriptional factors (TFs) having binding sites on the promoters of CqAMTs, nitrate transporters (CqNRTs), and glutamine-synthases (CqGSs), as well as the putative TF expression being correlated with the phenotypes and activities of GSs, glutamate synthase (GOGAT), nitrite reductase (NiR), and nitrate reductase (NR) of quinoa roots. The results showed a total of 12 members of the CqAMT family with varying expressions in different organs and in the same organs at different developmental stages. Complementation expression analyses in the triple mep1/2/3 mutant of yeast showed that except for CqAMT2.2b, 11/12 CqAMTs restored the uptake of NH4+ in the host yeast. CqAMT1.2a was found to mainly locate on the cell membrane, while TFs (e.g., CqNLPs, CqG2Ls, B3 TFs, CqbHLHs, CqZFs, CqMYBs, CqNF-YA/YB/YC, CqNACs, and CqWRKY) were predicted to be predominantly involved in the regulation, transportation, and assimilation of nitrogen. These results provide the functions of CqAMTs and their possible regulatory networks, which will lead to improved nitrogen use efficiency (NUE) in quinoa as well as other major crops.

Dissecting Cytophagalysin: Structural and Biochemical Studies of a Bacterial Pappalysin-Family Metallopeptidase

Cytophaga is a genus of Gram-negative bacteria occurring in soil and the gut microbiome. It is closely related to pathogenic Flavobacterium spp. that cause severe diseases in fish. Cytophaga strain L43-1 secretes cytophagalysin (CPL1), a 137 kDa peptidase with reported collagenolytic and gelatinolytic activity. We performed highly-confident structure prediction calculations for CPL1, which identified 11 segments and domains, including a signal peptide for secretion, a prosegment (PS) for latency, a metallopeptidase (MP)-like catalytic domain (CD), and eight immunoglobulin (Ig)-like domains (D3-D10). In addition, two short linkers were found at the D8-D9 and D9-D10 junctions, and the structure would be crosslinked by four disulfide bonds. The CPL1 CD was found closest to ulilysin from Methanosarcina acetivorans, which assigns CPL1 to the lower-pappalysin family within the metzincin clan of MPs. Based on the structure predictions, we aimed to produce constructs spanning the full-length enzyme, as well as PS+CD, PS+CD+D3, and PS+CD+D3+D4. However, we were successful only with the latter three constructs. We could activate recombinant CPL1 by PS removal employing trypsin, and found that both zymogen and mature CPL1 were active in gelatin zymography and against a fluorogenic gelatin variant. This activity was ablated in a mutant, in which the catalytic glutamate described for lower pappalyins and other metzincins was replaced by alanine, and by a broad-spectrum metal chelator. Overall, these results proved that our recombinant CPL1 is a functional active MP, thus supporting the conclusions derived from the structure predictions.

Developing a universal multi-epitope protein vaccine candidate for enhanced borna virus pandemic preparedness

Introduction

Borna disease virus 1 (BoDV-1) is an emerging zoonotic RNA virus that can cause severe acute encephalitis with high mortality. Currently, there are no effective countermeasures, and the potential risk of a future outbreak requires urgent attention. To address this challenge, the complete genome sequence of BoDV-1 was utilized, and immunoinformatics was applied to identify antigenic peptides suitable for vaccine development.

Methods

Immunoinformatics and antigenicity-focused protein screening were employed to predict B-cell linear epitopes, B-cell conformational epitopes, and cytotoxic T lymphocyte (CTL) epitopes. Only overlapping epitopes with antigenicity greater than 1 and non-toxic, non-allergenic properties were selected for subsequent vaccine construction. The epitopes were linked using GPGPG linkers, incorporating β-defensins at the N-terminus to enhance immune response, and incorporating Hit-6 at the C-terminus to improve protein solubility and aid in protein purification. Computational tools were used to predict the immunogenicity, physicochemical properties, and structural stability of the vaccine. Molecular docking was performed to predict the stability and dynamics of the vaccine in complex with Toll-like receptor 4 (TLR-4) and major histocompatibility complex I (MHC I) receptors. The vaccine construct was cloned through in silico restriction to create a plasmid for expression in a suitable host.

Results

Among the six BoDV-1 proteins analyzed, five exhibited high antigenicity scores. From these, eight non-toxic, non-allergenic overlapping epitopes with antigenicity scores greater than 1 were selected for vaccine development. Computational predictions indicated favorable immunogenicity, physicochemical properties, and structural stability. Molecular docking analysis showed that the vaccine remained stable in complex with TLR-4 and MHC I receptors, suggesting strong potential for immune recognition. A plasmid construct was successfully generated, providing a foundation for the experimental validation of vaccines in future pandemic scenarios.

Discussion

These findings demonstrate the potential of the immunoinformatics-designed multi-epitope vaccines for the prevention and treatment of BoDV-1. Relevant preparations were made in advance for possible future outbreaks and could be quickly utilized for experimental verification.

Molecular dynamics simulations revealed structural differences among TTHA1873-DNA interaction

This study aims to investigate the comparative binding pattern of TTHA1873 and its mutants (R55A and R138A) with DNA through molecular docking and molecular dynamics (MD) simulations. The docking results suggests that the Wild type (WT-TTHA1873), R55A, R138A and double mutant R55A/R138A having docking scores of -225.80 kcal/mol, -209.81 kcal/mol, -197.53 kcal/mol, -195.55 kcal/mol respectively and WT-TTHA1873 has more significant binding capability with DNA in comparison to mutants. The MD analysis revealed that the WT-TTHA1873 demonstrated stable interactions with DNA and exhibited a reduced conformational space compared to the mutants. By examining the atomic interactions, it was observed that significant variations in the hydrogen bonding pattern between WT-TTHA1873 and its mutants while interacting with DNA resulted in structural anomalies in the mutants and differences in DNA-binding specificity. The calculated binding free energies imply more stability of the WT-TTHA1873-DNA complex, while the mutants showed lesser binding affinity toward its interacting partner, double-stranded DNA. It is apparent that substituting single mutation R55A and R138A on TTHA1873 abolishes their DNA-binding ability. The present study portrays the critical role of R55 and R138 from TTHA1873 as likely involved in DNA binding.

The conformation of the nSrc specificity-determining loop in the Src SH3 domain is modulated by a WX conserved sequence motif found in SH3 domains

Cellular signaling networks are modulated by multiple protein-protein interaction domains that coordinate extracellular inputs and processes to regulate cellular processes. Several of these domains recognize short linear motifs, or SLiMs, which are often highly conserved and are closely regulated. One such domain, the Src homology 3 (SH3) domain, typically recognizes proline-rich SLiMs and is one of the most abundant SLiM-binding domains in the human proteome. These domains are often described as quite versatile, and indeed, SH3 domains can bind ligands in opposite orientations dependent on target sequence. Furthermore, recent work has identified diverse modes of binding for SH3 domains and a wide variety of sequence motifs that are recognized by various domains. Specificity is often attributed to the RT and nSrc loops near the peptide-binding cleft in this domain family, particularly for Class I binding, which is defined as RT and nSrc loop interactions with the N-terminus of the ligand. Here, we used the Src and Abl SH3 domains as a model to further investigate the role of the RT and nSrc loops in SH3 specificity. We created chimeric domains with both the RT and nSrc loop sequences swapped between these SH3 domains, and used fluorescence anisotropy assays to test how relative binding affinities were affected for Src SH3- and Abl SH3-specific ligands. We also used Alphafold-Multimer to model our SH3:peptide complexes in combination with molecular dynamics simulations. We identified a position that contributes to the nSrc loop conformation in Src SH3, the amino acid immediately following a highly conserved Trp that creates a hydrophobic pocket critical for SH3 ligand recognition. We defined this as the WX motif, where X = Trp for Src and Cys for Abl. A broad importance of this position for modulating nSrc loop conformation in SH3 domains is suggested by analyses of previously deposited SH3 structures, multiple sequence alignment of SH3 domains in the human proteome, and our biochemical and computational data of mutant Src and Abl SH3 domains. Overall, our work uses experimental approaches and structural modeling to better understand specificity determinants in SH3 domains.

Unraveling the interaction between a glycolytic regulator protein EhPpdk and an anaphase promoting complex protein EhApc10: yeast two hybrid screening, in vitro binding assays and molecular simulation study

The anaphase promoting complex (APC or cyclosome) is a major ubiquitin ligase that coordinates mitotic and G1 progression, acting as a major regulator of chromosome segregation. While the human APC contains fourteen subunits, it is yet to be explored in the pathogen Entamoeba histolytica. Our study reveals the existence of a single functional Apc10 homolog in E. histolytica, which acts as a processivity factor of ubiquitin ligase activity in human. A cDNA library generated from HM1:IMSS strain of E. histolytica was screened for interaction partners of EhApc10 in yeast two hybrid study. The novel interactor, a glycolytic enzyme, pyruvate phosphate dikinase (Ppdk) was found to interact with EhApc10 and further validated by in vitro assay. A comprehensive in silico study has emphasized the structural and functional aspects, encompassing physicochemical traits, predictive 3D structure modelling, validation of EhApc10-EhPpdk interaction through molecular docking and simulation. The interplay between a cell cycle protein and a glycolytic enzyme highlights the connection between cellular metabolism and the cell cycle regulatory mechanism. The study serves as the groundwork for future research on the non-mitotic role of APC beyond cell cycle.