Protein Structure and Function Prediction Using I-TASSER

The role of N-terminal acetylation of COVID fusion peptides in the interactions with liquid-ordered lipid bilayers

The partitioning of viral fusion peptides in lipid membranes with varying order was investigated due to the fusion mechanism being a potential therapeutic approach. Using a planar bilayer model and advanced techniques such as neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), the structural aspects of peptide-lipid interactions were explored. The study focused on two target membranes: one forming a liquid-ordered domain and the other forming a liquid-disordered domain. Surprisingly, the COVID fusion peptide did not bind significantly to either membrane, as demonstrated by both QCM-D and NR data, suggesting negligible or no interaction with the bilayers. However, the acetylated COVID fusion peptide showed distinct behaviour, indicating a crucial role of N-terminal acetylation in binding to cholesterol-rich liquid-ordered domains. The acetylated peptide induced changes in the structure and thickness of the ordered bilayer with cholesterol whereas proteins and peptides commonly only bind to disordered phases. This study provides valuable insights into the mechanisms of viral membrane fusion and highlights the importance of acetylation in influencing peptide-lipid interactions, laying the groundwork for potential antiviral therapeutic strategies.

Improving Protein Secondary Structure Prediction by Deep Language Models and Transformer Networks

Protein secondary structure prediction is useful for many applications. It can be considered a language translation problem, that is, translating a sequence of 20 different amino acids into a sequence of secondary structure symbols (e.g., alpha helix, beta strand, and coil). Here, we develop a novel protein secondary structure predictor called TransPross based on the transformer network and attention mechanism widely used in natural language processing to directly extract the evolutionary information from the protein language (i.e., raw multiple sequence alignment [MSA] of a protein) to predict the secondary structure. The method is different from traditional methods that first generate a MSA and then calculate expert-curated statistical profiles from the MSA as input. The attention mechanism used by TransPross can effectively capture long-range residue-residue interactions in protein sequences to predict secondary structures. Benchmarked on several datasets, TransPross outperforms the state-of-art methods. Moreover, our experiment shows that the prediction accuracy of TransPross positively correlates with the depth of MSAs, and it is able to achieve the average prediction accuracy (i.e., Q3 score) above 80% for hard targets with few homologous sequences in their MSAs. TransPross is freely available at https://github.com/BioinfoMachineLearning/TransPro .

Epitopes mapping for identification of potential cross-reactive peptide against leptospirosis

Leptospira, the pathogenic helical spirochetes that cause leptospirosis, is an emerging zoonotic disease with effective dissemination tactics in the host and can infect humans and animals with moderate or severe illnesses. Thus, peptide-based vaccines may be the most effective strategy to manage the immune response against Leptospira to close these gaps. In the current investigation, highly immunogenic proteins from the proteome of Leptospira interorgan serogroup Icterohaemorrhagie serovar Lai strain 56601 were identified using immunoinformatic methods. It was discovered that the conserved and most immunogenic outer membrane Lepin protein was both antigenic and non-allergenic by testing 15 linear B-cells and the ten best T-cell (Helper-lymphocyte (HTL) with the most significant number of HLA-DR binding alleles and the eight cytotoxic T lymphocyte (CTL)) epitopes. Furthermore, a 3D structural model of CTL epitopes was created using the Pep-Fold3 platform. Using the Autodock 4.2 docking server, research was conducted to determine how well the top-ranked CTL peptide models attach to HLA-A*0201 (PDB ID: 4U6Y). With HLA-A*0201, the epitope SSGTGNLHV binds with a binding energy of -1.29 kcal/mol. Utilizing molecular dynamics modeling, the projected epitope-allele docked complex structure was optimized, and the stability of the complex system was assessed. Therefore, this epitope can trigger an immunological response and produce effective Leptospira vaccine candidates. Overall, this study offers a unique vaccination candidate and may encourage additional research into leptospirosis vaccines.Communicated by Ramaswamy H. Sarma.

The banana MaFLA27 confers cold tolerance partially through modulating cell wall remodeling

Fasciclin-like arabinogalactan proteins (FLAs) have been shown to improve plant tolerance to salt stress. However, their role in cold tolerance (CT) remains unclear. Here, we report that banana MaFLA27 positively regulates CT in Arabidopsis. MaFLA27-overexpression (OE) caused the upregulation of differentially expressed arabinogalactan proteins (AGPs) and genes involved in the biosynthesis of cellulose, lignin, and xylan, as well as the degradation of pectin and xyloglucan. Correspondingly, MaFLA27-OE plants exhibited increased cell wall thickness, enhanced cellulose lignin and starch granule content, elevated levels of partially homogalacturonans recognized by JIM5 and JIM7 antibodies, xyloglucan components recognized by CCRC-M39/104 and LM15 antibodies, LM14 antibody binding AGPs. In contrast, transgenic plants showed a decreased degree of pectin methyl-esterification and accumulated less reactive oxygen species after cold acclimation when compared to wild-type plants. A higher number of pectin methylesterases and cellulose and xylan biosynthesis genes were elevated after cold acclimation. Additionally, both Arabidopsis mutant cesa8 and cellulose inhibitor-treated plants displayed decreased freezing tolerance. Our data suggested that MaFLA27-OE in Arabidopsis may perceive and transmit low-temperature stress signals to the cellulose synthase complexes, activating cellulose synthesis and enhancing cold tolerance. These findings reveal a previously unreported cold-tolerance function of FLAs and highlight associated cell wall-mediated tolerance mechanisms.

Molecular and ultrastructural characteristics of virulent and attenuated vaccine strains of goose parvovirus LIV-22

The disease caused by goose parvovirus (GPV) affects young goslings and ducks and leads to substantial losses for farmers due to high mortality rates, reaching 70-100% in naive flocks. Here, we present the results of a study focusing on the historical virulent GPV LIV-22 strain, which was isolated in the USSR in 1972. An attenuated GPV LIV-22 vaccine strain that was generated by continuous passaging in goose embryonic fibroblasts was also studied. Phylogenetic analysis placed both GPV LIV-22 strains in the classical GPV group, close to the vaccine and low-pathogenic strains. However, several individual changes in the GPV LIV-22 VP1 gene highlight the uniqueness of the evolution and adaptation mechanism of GPV LIV-22 strains. Transmission electron microscopy (TEM) revealed severe ultrastructural changes in goose hepatocytes and enterocytes as early as 24-48 h postinfection, confirming abrupt GPV pathogenesis. This description of some of the essential characteristics of the GPV LIV-22 virulent and vaccine strain will be useful for studying GPV evolution and molecular pathogenesis.

GhLCYε-3 characterized as a lycopene cyclase gene responding to drought stress in cotton

Drought stress significantly affects crop productivity. Carotenoids are essential photosynthetic pigment for plants, bacteria, and algae, with signaling and antioxidant functions. Lutein is a crucial branch product in the carotenoid synthesis pathway, which effectively improves the stress tolerance of higher plants. lycopene cyclase, a central enzyme for lutein synthesis, holds great significance in regulating lutein production. This research establishes a correlation between lutein content and stress resistance by measuring the drought resistance and lutein content of various cotton materials. To identify which crucial genes are associated with lutein, the lycopene cyclase family (LCYs) was analyzed. The research found that LCYs form a highly conserved family divided into two subfamilies, LCY-ε (lycopene ε-cyclase) and LCY-β (lycopene β-cyclase). Most members of the LCY family contain photoresponsive elements and abscisic acid elements. qRT-PCR demonstrates showed that most genes responded positively to drought stress, and GhLCYε-3 was expressed significantly differently under drought stress. Virus-induced gene silencing (VIGS) assay showed that the content of GhLCYε-3 was significantly increased with MDA and PRO, and the contents of chlorophyll and lutein were significantly decreased in pYL156 plants. The decrease in GhLCYε-3 expression is speculated to lead to reduced lutein content in vivo, resulting in the accumulation of reactive oxygen species (ROS) and decreased drought tolerance. This research enriched the understanding of LCY gene family and lutein function, and provided a new reference for cotton planting in arid areas.

Synopsis

Lycopene cyclase plays an important role in enhancing the ability of scavenging ROS and drought resistance of plants.

Preparation of recombinant type I collagen (PF-I-80) and its functional characterization and biomedical applications in wound healing

This study evaluates the potential applications of recombinant PF-I-80 protein in regenerative medicine and the treatment of inflammatory diseases, focusing on its effects on cell migration, differentiation, and anti-inflammatory properties. Various in vitro assays were conducted, including scratch assays, Transwell experiments, RT-PCR and Western Blot to analyze gene and protein expression related to differentiation and inflammation, and immunofluorescence staining to observe cellular changes. The results indicated that PF-I-80 significantly promoted cell migration, highlighting its potential in tissue repair and regeneration. It also enhanced cell differentiation, demonstrating its applicability in tissue repair, and showed significant anti-inflammatory effects by reducing the expression of pro-inflammatory cytokines. In animal models, PF-I-80 notably reduced levels of inflammatory factors IL-1β and TNF-α, shortened the inflammatory phase, and accelerated wound healing. Additionally, PF-I-80 increased FGF-2 levels, which promoted the proliferation of endothelial and fibroblast cells and enhanced collagen synthesis. These in vitro and in vivo findings position PF-I-80 as a promising biomaterial for applications in regenerative medicine and inflammatory disease treatment.

The arms race of ray-finned fish against the derepression of LTR retroelements

Transposable elements (TEs) are dynamic components of eukaryotic genomes that are important in speciation. Evolutionary pressure has led to the evolution of silencing mechanisms, such as the involvement of Krüppel box-associated zinc finger proteins (KRAB-ZFPs) that regulate TE transcription via TRIM28-mediated NuRD complex recruitment in sarcopterygians. The transcriptional activity of genes that encode proteins of the NuRD complex was also detected in actinopterygians, even though genes encoding crucial proteins such as TRIM28 and KRAB-ZNF are missing in this evolutionary lineage. This study hypothesized that TRIM33 could serve as a replacement for TRIM28 in this lineage and suggested an evolutionary relationship between the sarcopterygian KRAB-ZFPs, the actinopterygian KRAB-like ZNFs and the fish N-terminal zinc-associated domain (FiNZ) ZNFs that are specific to cyprinids. Through coimmunoprecipitation analyses conducted in the basal teleost Anguilla marmorata and the cyprinid Danio rerio, we evaluated the interaction of TRIM33 with the KRAB-like and FiNZ ZFPs. Moreover, we explored the expression profiles of TEs and genes that encode proteins that are potentially involved in their silencing during zebrafish development. Our results may provide new evidence for more thoroughly understanding the evolutionary mechanisms underlying the adaptation of ray-finned fish.

Exploring the Potential of Biomimetic Peptides in Targeting Fibrillar and Filamentous Alpha-Synuclein-An In Silico and Experimental Approach to Parkinson's Disease

Alpha-synuclein (ASyn) is a protein that is known to play a critical role in Parkinson's disease (PD) due to its propensity for misfolding and aggregation. Furthermore, this process leads to oxidative stress and the formation of free radicals that cause neuronal damage. In this study, we have utilized a biomimetic approach to design new peptides derived from marine natural resources. The peptides were designed using a peptide scrambling approach where antioxidant moieties were combined with fibrillary inhibition motifs in order to design peptides that would have a dual targeting effect on ASyn misfolding. Of the 20 designed peptides, 12 were selected for examining binding interactions through molecular docking and molecular dynamics approaches, which revealed that the peptides were binding to the pre-NAC and NAC (non-amyloid component) domain residues such as Tyr39, Asn65, Gly86, and Ala85, among others. Because ASyn filaments derived from Lewy body dementia (LBD) have a different secondary structure compared to pathogenic ASyn fibrils, both forms were tested computationally. Five of those peptides were utilized for laboratory validation based on those results. The binding interactions with fibrils were confirmed using surface plasmon resonance studies, where EQALMPWIWYWKDPNGS, PYYYWKDPNGS, and PYYYWKELAQM showed higher binding. Secondary structural analyses revealed their ability to induce conformational changes in ASyn fibrils. Additionally, PYYYWKDPNGS and PYYYWKELAQM also demonstrated antioxidant properties. This study provides insight into the binding interactions of varying forms of ASyn implicated in PD. The peptides may be further investigated for mitigating fibrillation at the cellular level and may have the potential to target ASyn.

CRB2 Depletion Induces YAP Signaling and Disrupts Mechanosensing in Podocytes

Focal Segmental Glomerulosclerosis (FSGS) is a histologic lesion caused by a variety of injurious stimuli that lead to dysfunction/loss of glomerular visceral epithelial cells (i.e. podocytes). Pathogenic mutations in CRB2, encoding the type 1 transmembrane protein Crumb 2 Homolog Protein, have been shown to cause early-onset corticosteroid-resistant nephrotic syndrome (SRNS)/FSGS. Here, we identified a 2-generation East Asian kindred (DUK40595) with biopsy-proven SRNS/FSGS caused by a compound heterozygous mutation in CRB2 comprised of the previously described truncating mutation p.Gly1036_Alafs*43 and a rare 9-bp deletion mutation p.Leu1074_Asp1076del. Because compound heterozygous mutations involving the truncating p.Gly1036_Alafs*43 variant have been associated with reduced CRB2 expression in podocytes and autosomal recessive SRNS/FSGS, we sought to define the pathogenic effects of CRB2 deficiency in podocytes. We show that CRB2 knockdown induces YAP activity and target gene expression in podocytes. It upregulates YAP-mediated mechanosignaling and increases the density of focal adhesion and F-actin. Using Elastic Resonator Interference Stress Microscopy (ERISM), we demonstrate that CRB2 knockdown also enhances podocyte contractility in a substrate stiffness-dependent manner. The knockdown effect decreases with increasing substrate stiffness, indicating impaired mechanosensing in CRB2 knockdown cells at low substrate stiffness. While the mechanical activation of CRB2 knockdown cells is associated with increased YAP activity, the enhanced cell contractility is not significantly reduced by the selective YAP inhibitors K-975 and verteporfin, suggesting that multiple pathways may be involved in mechanosignaling downstream of CRB2. Taken together, these studies provide the first evidence that CRB2 deficiency may impair podocyte mechanotransduction via disruption of YAP signaling in podocytes.

Pangenome analysis of five representative Tropheryma whipplei strains following multiepitope-based vaccine design via immunoinformatic approaches

Whipple disease caused by Tropheryma whipplei a gram-positive bacterium is a systemic disorder that impacts not only the gastrointestinal tract but also the vascular system, joints, central nervous system, and cardiovascular system. Due to the lack of an approved vaccine, this study aimed to utilize immunoinformatic approaches to design multiepitope -based vaccine by utilizing the proteomes of five representative T. whipplei strains. The genomes initially comprised a total of 4,844 proteins ranging from 956 to 1012 proteins per strain. We collected 829 nonredundant lists of core proteins, that were shared among all the strains. Following subtractive proteomics, one extracellular protein, WP_033800108.1, a WhiB family transcriptional regulator, was selected for the chimeric-based multiepitope vaccine. Five immunodominant epitopes were retrieved from the WhiB family transcriptional regulator protein, indicating MHC-I and MHC-II with a global population coverage of 70.61%. The strong binding affinity, high solubility, nontoxicity, nonallergenic properties and high antigenicity scores make the selected epitopes more appropriate. Integration of the epitopes into a chimeric vaccine was carried out by applying appropriate adjuvant molecules and linkers, leading to the vaccine construct having enhanced immunogenicity and successfully eliciting both innate and adaptive immune responses. Moreover, the abilityof the vaccine to bind TLR4, a core innate immune receptor, was confirmed. Molecular dynamics simulations have also revealed the promising potential stability of the designed vaccine at 400 ns. In summary, we have designed a potential vaccine construct that has the ability not only to induce targeted immunogenicity for one strain but also for global T. whipplei strains. This study proposes a potential universal vaccine, reducing Whipple's disease risk and laying the groundwork for future research on multi-strain pathogens.

Characterization of a Xylosyltransferase from Panax notoginseng Catalyzing Ginsenoside 2'-O Glycosylation in the Biosynthesis of Notoginsenosides

Notoginsenosides are important bioactive compounds from Panax notoginseng (Burk.) F. H. Chen, most of which have xylose in their sugar chains. However, the xylosyltransferases involved in the generation of notoginsenosides remain poorly understood, posing a bottleneck for further study of the biosynthesis of notoginsenosides. In this work, a new xylosyltransferase gene, PnUGT57 (named UGT94BW1), was identified from P. notoginseng, which has a distinct sequence and could catalyze the 2'-O glycosylation of ginsenosides Rh1 and Rg1 to produce notoginsenosides R2 and R1, respectively. We first characterized the optimal conditions for the PnUGT57 activity and its enzymatic kinetic parameters, and then, molecular docking and site-directed mutagenesis were performed to elucidate the catalytic mechanism of PnUGT57. Combined with the results of site-directed mutagenesis, Glu26, Ser266, Glu267, Trp347, Ser348, and Glu352 in PnUGT57 were identified as the key residues involved in 2'-O glycosylation of C-6 O-Glc, and PnUGT57R175A and PnUGT57G237A could significantly improve the catalytic activity of PnUGT57. These findings not only provide a new xylosyltransferase gene for augmenting the plant xylosyltransferase database but also identify the pivotal sites and catalytic mechanism of the enzyme, which would provide reference for the modification and application of xylosyltransferases in the future.

Structural impact, ligand-protein interactions, and molecular phenotypic effects of TGF-β1 gene variants: In silico analysis with implications for idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic Pulmonary Fibrosis (IPF) is a chronic interstitial lung disease resulting in progressively deteriorating lung function. Transforming growth factor-β1 (TGF-β1) belongs to the TGF superfamily and exerts a profibrotic role in promoting lung fibrosis by facilitating fibroblast infiltration and activity, extracellular matrix deposition, and inhibition of collagen breakdown, thus promoting tissue remodelling and IPF.

MATERIALS AND METHODS

We evaluated the link between pathogenic TGF-β1 SNPs and IPF pathogenesis and the structure-activity functional consequences of those SNPs on the TGF-β1 protein. Several computational algorithms were merged to address the functional consequences of TGF-β1 gene mutations to protein stability, putative post-translational modification sites, ligand-protein interactions, and molecular phenotypic effects. These included FATHMM, POLYPHEN2, PROVEAN, and SIFT tools (identifying deleterious nsSNPs in the TGF-β1 gene), along with Pmut, PhD-SNP, SNAP, MutPred and the related TMHMM, MARCOIL, and DisProt algorithms (predicting structural disorders). INPS-MD was also used to evaluate the mutation-induced TGF-β1 protein's stability and MODPRED for recognition of post-translational TGF-β1 modification.

RESULTS

In total, 14 major pathogenic variants markedly impact the destabilization of the TGF-β1 protein, with most of these high-risk mutations associated with decreased stability of the TGF-β1 protein as per the I-Mutant, MUpro, and INPS-MD tools. R205W, R185W, R180Q, D86Y, and I300T variants were proposed to participate in the post-translational modifications, thus affecting affect protein-ligand interactions. Furthermore, at-risk genetic variants appear to target conserved regions in the alpha helices, random coils, and extracellular loops, resulting in a varied composition of amino acids, charge, hydrophobicity, and spatial architecture.

CONCLUSIONS

This study manuscript comprehensively analyzes gene variants within the TGF-β1 gene, offering novel insights into their structural and functional implications in interacting with target sites. This study is significant for the development of targeted therapeutic strategies and personalized treatment approaches for patients with inflammatory lung diseases such as IPF.

NEKL-4 regulates microtubule stability and mitochondrial health in ciliated neurons

Ciliopathies are often caused by defects in the ciliary microtubule core. Glutamylation is abundant in cilia, and its dysregulation may contribute to ciliopathies and neurodegeneration. Mutation of the deglutamylase CCP1 causes infantile-onset neurodegeneration. In C. elegans, ccpp-1 loss causes age-related ciliary degradation that is suppressed by a mutation in the conserved NEK10 homolog nekl-4. NEKL-4 is absent from cilia, yet it negatively regulates ciliary stability via an unknown, glutamylation-independent mechanism. We show that NEKL-4 was mitochondria-associated. Additionally, nekl-4 mutants had longer mitochondria, a higher baseline mitochondrial oxidation state, and suppressed ccpp-1∆ mutant lifespan extension in response to oxidative stress. A kinase-dead nekl-4(KD) mutant ectopically localized to ccpp-1∆ cilia and rescued degenerating microtubule doublet B-tubules. A nondegradable nekl-4(PEST∆) mutant resembled the ccpp-1∆ mutant with dye-filling defects and B-tubule breaks. The nekl-4(PEST∆) Dyf phenotype was suppressed by mutation in the depolymerizing kinesin-8 KLP-13/KIF19A. We conclude that NEKL-4 influences ciliary stability by activating ciliary kinesins and promoting mitochondrial homeostasis.

Binding structures of SERF1a with NT17-polyQ peptides of huntingtin exon 1 revealed by SEC-SWAXS, NMR and molecular simulation

The aberrant fibrillization of huntingtin exon 1 (Httex1) characterized by an expanded polyglutamine (polyQ) tract is a defining feature of Huntington's disease, a neurodegenerative disorder. Recent investigations underscore the involvement of a small EDRK-rich factor 1a (SERF1a) in promoting Httex1 fibrillization through interactions with its N terminus. By establishing an integrated approach with size-exclusion-column-based small- and wide-angle X-ray scattering (SEC-SWAXS), NMR, and molecular simulations using Rosetta, the analysis here reveals a tight binding of two NT17 fragments of Httex1 (comprising the initial 17 amino acids at the N terminus) to the N-terminal region of SERF1a. In contrast, examination of the complex structure of SERF1a with a coiled NT17-polyQ peptide (33 amino acids in total) indicates sparse contacts of the NT17 and polyQ segments with the N-terminal side of SERF1a. Furthermore, the integrated SEC-SWAXS and molecular-simulation analysis suggests that the coiled NT17 segment can transform into a helical conformation when associated with a polyQ segment exhibiting high helical content. Intriguingly, NT17-polyQ peptides with enhanced secondary structures display diminished interactions with SERF1a. This insight into the conformation-dependent binding of NT17 provides clues to a catalytic association mechanism underlying SERF1a's facilitation of Httext1 fibrillization.

A glimpse into the hidden world of the flexible C-terminal protein binding domains of human RAD52

Human RAD52 protein binds DNA and is involved in genomic stability maintenance and several forms of DNA repair, including homologous recombination and single-strand annealing. Despite its importance, there are very few structural details about the variability of the RAD52 ring size and the RAD52 C-terminal protein-protein interaction domains. Even recent attempts to employ cryogenic electron microscopy (cryoEM) methods on full-length yeast and human RAD52 do not reveal interpretable structures for the C-terminal half that contains the replication protein A (RPA) and RAD51 binding domains. In this study, we employed the monodisperse purification of two RAD52 deletion constructs and small angle X-ray scattering (SAXS) to construct a structural model that includes RAD52's RPA binding domain. This model is of interest to DNA repair specialists as well as for drug development against HR-deficient cancers.

Spectrum of Clinical Variability with SEPT9 Gene Mutation in Hereditary Neuralgic Amyotrophy: Understanding the Pathogenesis Using Molecular Dynamics Simulation Study

BACKGROUND

Hereditary Neuralgic Amyotrophy (HNA) is an autosomal dominant disorder characterized by episodes of severe pain and amyotrophy affecting the brachial plexus as well as other sites. Mutations in the SEPTIN9 gene have been identified as genetic abnormality for HNA. Although the genetic mutations are known, their pathogenesis for the causation of this disorder is not exactly elucidated.

OBJECTIVE

In this study, we have investigated the phenotypic and genetic features in a large pedigree with HNA.

METHODS

We report the clinical spectrum and genetic analysis of a family with 9 affected members. Clinical heterogeneity has been reported in the individuals having mutations in SEPTIN9 gene. After taking informed consent, we have done genetic analysis of 6 affected and 4 unaffected members of the family to identify the molecular abnormalities of SEPTIN9 gene.

RESULTS AND CONCLUSIONS

Genetic analysis has identified the presence of NM_001113491.2:p.Arg106Trp mutation in SEPTIN9 gene. The same mutation has been identified in 6 affected members of the family. Molecular simulation study has revealed that the mutation has significantly altered the conformation of septin-9 protein, thereby impairing the microtubule binding and bundling ability. Although the affected members shared a common recurrent mutation, they have a wide spectrum of clinical variability. This may be due to the variable penetrance of the mutation and different epigenetic influences in the family. This is the first genetically confirmed case series of HNA reported from India.

A single-cell atlas of the miracidium larva of Schistosoma mansoni reveals cell types, developmental pathways, and tissue architecture

Schistosoma mansoni is a parasitic flatworm that causes the major neglected tropical disease schistosomiasis. The miracidium is the first larval stage of the life cycle. It swims and infects a freshwater snail, transforms into a mother sporocyst, where its stem cells generate daughter sporocysts that give rise to human-infective cercariae larvae. To understand the miracidium at cellular and molecular levels, we created a whole-body atlas of its ~365 cells. Single-cell RNA sequencing identified 19 transcriptionally distinct cell clusters. In situ hybridisation of tissue-specific genes revealed that 93% of the cells in the larva are somatic (57% neural, 19% muscle, 13% epidermal or tegument, 2% parenchyma, and 2% protonephridia) and 7% are stem. Whereas neurons represent the most diverse somatic cell types, trajectory analysis of the two main stem cell populations indicates that one of them is the origin of the tegument lineage and the other likely contains pluripotent cells. Furthermore, unlike the somatic cells, each of these stem populations shows sex-biased transcriptional signatures suggesting a cell-type-specific gene dosage compensation for sex chromosome-linked loci. The miracidium represents a simple developmental stage with which to gain a fundamental understanding of the molecular biology and spatial architecture of schistosome cells.

Genome Mining of Pseudanabaena galeata CCNP1313 Indicates a New Scope in the Search for Antiproliferative and Antiviral Agents

Compounds derived from natural sources pave the way for novel drug development. Cyanobacteria is an ubiquitous phylum found in various habitats. The fitness of those microorganisms, within different biotopes, is partially dependent on secondary metabolite production. Their enhanced production under biotic/abiotic stress factors accounts for better survival rates of cells, and thereby cyanobacteria are as an enticing source of bioactive compounds. Previous studies have shown the potent activity of extracts and fractions from Pseudanabaena galeata (Böcher 1949) strain CCNP1313 against cancer cells and viruses. However, active agents remain unknown, as the selected peptides had no effect on the tested cell lines. Here, we present a bottom-up approach, pinpointing key structures involved in secondary metabolite production. Consisting of six replicons, a complete genome sequence of P. galeata strain CCNP1313 was found to carry genes for non-ribosomal peptide/polyketide synthetases embedded within chromosome spans (4.9 Mbp) and for a ribosomally synthesized peptide located on one of the plasmids (0.2 Mbp). Elucidation of metabolite synthesis pathways led to prediction of their structure. While none of the synthesis-predicted products were found in mass spectrometry analysis, unexplored synthetases are characterized by structural similarities to those producing potent bioactive compounds.

Frequency of Fusion Inhibitor Resistance Mutations Among Therapy-Naïve HIV Patients

Glycoprotein 41 (gp41) of the human immunodeficiency virus type 1 (HIV-1) protein plays a critical role in membrane fusion. Gp41 binds to proteins in the plasma membrane of CD4+ T cells, particularly the T-cell antigen receptor (TCR). These findings indicate that gp41 is involved in the assembly of HIV-1 at the plasma membrane of T cells and affects the stimulation of the TCR. To control HIV-1, new inhibitors were introduced to target the gp41 protein. However, mutations in this region might reduce their efficacy. The Gp41 region was amplified from the sera of 30 patients using nested polymerase chain reaction. The sequences were analyzed by bioinformatics tools to identify mutations and gp41 structural features. Subtyping and the interaction between fusion inhibitors and gp41 proteins were also examined. As the first report from Iran, docking analysis between fusion inhibitors and Iranian gp41 proteins showed that mutations in gp41 could not reduce the efficacy of the fusion inhibitors. Most of the patients were infected with CRF35-AD. Several post-modification positions, including glycosylation and phosphorylation sites, were identified in the gp41 protein. Our findings revealed no known multinational drug resistance to gp41 inhibitors; thus, fusion inhibitors can effectively inhibit HIV in Iranian patients. In addition, the present study introduced a new gp41 region (36-44 aa), which considerably influences the interactions between gp41 inhibitors and the gp41 protein. This region may play a pivotal role in suppressing gp41 inhibitors in CFR35-AD. Furthermore, gp41 can be considered a good target for subtyping analysis via the phylogenetic method.

Dominant B cell-T cell epitopes instigated robust immune response in-silico against Scrub Typhus

Scrub typhus, a potentially life-threatening infectious disease, is attributed to bacteria Orientia tsutsugamushi (O. tsutsugamushi). The transmission of this illness to humans occurs through the bite of infected chiggers, which are the larval forms of mites belonging to the genus Leptotrombidium. In this research, we developed a subunit vaccine specifically designed to target outer membrane proteins. Immunodominant cytotoxic T-lymphocytes (CTLs), B- lymphocytes (BCLs), and major histocompatibility complex (MHC)- II epitopes were identified using machine learning and bioinformatics approaches. These epitopes were arranged in different combinations with the help of suitable linkers like AAY, KK, GPGPG and adjuvant (cholera toxin B) that resulted in a vaccine construct. Physiochemical properties were assessed, where the predicted solubility (0.571) was higher than threshold value. Tertiary structure was predicted using I-TASSER web server and evaluated using Ramachandran plot (94 % residues in most favourable region) and z-score (-6.04), which had shown the structure to have good stability and residue arrangement. Molecular docking with immune receptors, Toll-like receptor (TLR)-2 and -4 showed good residue interaction with 13 and 5 hydrogen bonds respectively. Molecular dynamics simulations of receptor-ligand complex provided the idea about the strong interaction having 1.524751 × 10-5 eigenvalue. Amino acid sequence of vaccine was converted to nucleotide sequence and underwent codon optimization. The optimized codon sequence was used for in-silico cloning, which provided idea about the possibility of synthesis of vaccine using E. coli as host. Overall, this study provided a promising blueprint for a scrub typhus vaccine, although experimental validation is needed for confirmation. Furthermore, it is crucial to acknowledge that while bioinformatics provides valuable insights, in-vitro and in-vivo studies are imperative for a comprehensive evaluation of vaccine candidate. Thus, the integration of computational predictions with empirical research is essential to validate the efficacy, safety, and real-world applicability of the designed vaccine against Scrub Typhus. Nevertheless, the findings are good to carry forward for in-vitro and in-vivo investigations.

Structure elucidation of a multi-modular recombinant endoglucanase, AtGH9C-CBM3A-CBM3B from Acetivibrio thermocellus ATCC 27405 and its substrate binding analysis

Cellulases from GH9 family show endo-, exo- or processive endocellulase activity, but the reason behind the variation is unclear. A GH9 recombinant endoglucanase, AtGH9C-CBM3A-CBM3B from Acetivibrio thermocellus was structurally characterized for conformation, binding and dynamics assessment. Modeled AtGH9C-CBM3A-CBM3B depicted (α/α)6-barrel structure with Asp98, Asp101 and Glu489 acting as catalytic triad. CD results revealed 25.2 % α-helix, 18.4 % β-sheet and rest 56.4 % of random coils, corroborating with predictions from PSIPRED and SOPMA. MD simulation of AtGH9C-CBM3A-CBM3B bound cellotetraose showed structural stability and global compactness with lowered RMSD values (1.5 nm) as compared with only AtGH9C-CBM3A-CBM3B (1.8 nm) for 200 ns. Higher fluctuation in RMSF values in far-positioned CBM3B pointed to its redundancy in substrate binding. Docking studies showed maximum binding with cellotetraose (ΔG = -5.05 kcal/mol), with reduced affinity towards ligands with degree of polymerization (DP) lower (DP < 4) or higher than 4 (DP > 4). Processivity index displayed the enzyme to be processive with loop 3 (342-379 aa) possibly blocking the non-reducing end of cellulose chain, resulting in cellotetraose release. SAXS analysis of AtGH9C-CBM3A-CBM3B at 5 mg/mL displayed monodispersed state with fist-and-elbow shape in solution. Negative zeta potential of -24 mV at 5 mg/mL indicated stability and free from aggregation.

Study of key residues in MERS-CoV and SARS-CoV-2 main proteases for resistance against clinically applied inhibitors nirmatrelvir and ensitrelvir

The Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is an epidemic, zoonotically emerging pathogen initially reported in Saudi Arabia in 2012. MERS-CoV has the potential to mutate or recombine with other coronaviruses, thus acquiring the ability to efficiently spread among humans and become pandemic. Its high mortality rate of up to 35% and the absence of effective targeted therapies call for the development of antiviral drugs for this pathogen. Since the beginning of the SARS-CoV-2 pandemic, extensive research has focused on identifying protease inhibitors for the treatment of SARS-CoV-2. Our intention was therefore to assess whether these protease inhibitors are viable options for combating MERS-CoV. To that end, we used previously established protease assays to quantify inhibition of SARS-CoV-2, MERS-CoV and other main proteases. Nirmatrelvir inhibited several of these proteases, whereas ensitrelvir was less broadly active. To simulate nirmatrelvir's clinical use against MERS-CoV and subsequent resistance development, we applied a safe, surrogate virus-based system. Using the surrogate virus, we previously selected hallmark mutations of SARS-CoV-2-Mpro, such as T21I, M49L, S144A, E166A/K/V and L167F. In the current study, we selected a pool of MERS-CoV-Mpro mutants, characterized the resistance and modelled the steric effect of catalytic site mutants S142G, S142R, S147Y and A171S.

A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-β (Aβ42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the β-content of Aβ42; however, it increases the tendency of the monomeric Aβ42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aβ42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aβ42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

A common protein C inhibitor exosite partially controls the heparin induced activation and inhibition of serine proteases

Protein C inhibitor (PCI) maintains hemostasis by inhibiting both procoagulant and anticoagulant serine proteases, and plays important roles in coagulation, fibrinolysis, reproduction, and anti-angiogenesis. The reactive site loop of PCI traps and irreversibly inhibits the proteases like APC (activating protein C), thrombin (FIIa) and factor Xa (FXa). Previous studies on antithrombin (ATIII) had identified Tyr253 and Glu255 as functional exosites that interact and aid in the inhibition of factor IXa and FXa. Presence of exosite in PCI is not known, however a sequence comparison with the PCI from different vertebrate species and ATIII identified Glu239 to be absolutely conserved. PCI residues analogous to ATIII exosite residues were mutated to R238A and E239A. Purified variant PCI in the presence of heparin (10 μg/ml) showed a 2-4 fold decrease in the rate of inhibition of the proteases. However, the stoichiometry of inhibition of FIIa, APC, and FXa by native PCI, R238A and E239A variants were found to be close to 1.0, which also indicated the formation of stable complexes based on SDS-PAGE and western blot analysis with thrombin and APC. Our findings revealed the possible presence of an exosite in PCI that influences the protease inhibition rates.

Mutagenesis and functional analysis of the varicella-zoster virus portal protein

Herpesviruses replicate by cleaving concatemeric dsDNA into single genomic units that are packaged through an oligomeric portal present in preformed procapsids. In contrast to what is known about phage portal proteins, details concerning herpesvirus portal structure and function are not as well understood. A panel of 65 Varicella-Zoster virus (VZV) recombinant portal proteins with five amino acid in-frame insertions were generated by random transposon mutagenesis of the VZV portal gene, ORF54. Subsequently, 65 VZVLUC recombinant viruses (TNs) were generated via recombineering. Insertions were mapped to predicted portal domains (clip, wing, stem, wall, crown, and β-hairpin tunnel-loop) and recombinant viruses were characterized for plaque morphology, replication kinetics, pORF54 expression, and classified based on replication in non-complementing (ARPE19) or complementing (ARPE54C50) cell lines. The N- and C-termini were tolerant to insertion mutagenesis, as were certain clip sub-domains. The majority of mutants mapping to the wing, wall, β-hairpin tunnel loop, and stem domains were lethal. Elimination of the predicted ORF54 start codon revealed that the first 40 amino acids of the N-terminus were not required for viral replication. Stop codon insertions in the C-terminus showed that the last 100 amino acids were not required for viral replication. Lastly, a putative protease cleavage site was identified in the C-terminus of pORF54. Cleavage was likely orchestrated by a viral protease; however, processing was not required for DNA encapsidation and viral replication. The panel of recombinants should prove valuable in future studies to dissect mammalian portal structure and function.IMPORTANCEThough nucleoside analogs and a live-attenuated vaccine are currently available to treat some human herpesvirus family members, alternate methods of combating herpesvirus infection could include blocking viral replication at the DNA encapsidation stage. The approval of Letermovir provided proof of concept regarding the use of encapsidation inhibitors to treat herpesvirus infections in the clinic. We propose that small-molecule compounds could be employed to interrupt portal oligomerization, assembly into the capsid vertex, or affect portal function/dynamics. Targeting portal at any of these steps would result in disruption of viral DNA packaging and a decrease or absence of mature infectious herpesvirus particles. The oligomeric portals of herpesviruses are structurally conserved, and therefore, it may be possible to find a single compound capable of targeting portals from one or more of the herpesvirus subfamilies. Drug candidates from such a series would be effective against viruses resistant to the currently available antivirals.

KEGG orthology prediction of bacterial proteins using natural language processing

BACKGROUND

The advent of high-throughput technologies has led to an exponential increase in uncharacterized bacterial protein sequences, surpassing the capacity of manual curation. A large number of bacterial protein sequences remain unannotated by Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology, making it necessary to use auto annotation tools. These tools are now indispensable in the biological research landscape, bridging the gap between the vastness of unannotated sequences and meaningful biological insights.

RESULTS

In this work, we propose a novel pipeline for KEGG orthology annotation of bacterial protein sequences that uses natural language processing and deep learning. To assess the effectiveness of our pipeline, we conducted evaluations using the genomes of two randomly selected species from the KEGG database. In our evaluation, we obtain competitive results on precision, recall, and F1 score, with values of 0.948, 0.947, and 0.947, respectively.

CONCLUSIONS

Our experimental results suggest that our pipeline demonstrates performance comparable to traditional methods and excels in identifying distant relatives with low sequence identity. This demonstrates the potential of our pipeline to significantly improve the accuracy and comprehensiveness of KEGG orthology annotation, thereby advancing our understanding of functional relationships within biological systems.

KCD: A prediction web server of knowledge-based circular dichroism

We present a web server that predicts the far-UV circular dichroism (CD) spectra of proteins by utilizing their three-dimensional (3D) structures from the Protein Data Bank (PDB). The main algorithm is based on the classical theory of optical activity together with a set of atomic complex polarizabilities, which are obtained from the analysis of a series of synchrotron radiation CD spectra and their related 3D structures from the PDB. The results of our knowledge-based CD method (KCD) are in good agreement with measured spectra that could include the effect of D-amino acids. Our method also delivers some of the most accurate predictions, in comparison with the calculated spectra from well-established models. Specifically, using a metric of closeness based on normalized absolute deviations between experimental and calculated spectra, the mean values for a series of 57 test proteins give the following figures for such models: 0.26 KCD, 0.27 PDBMD2CD, 0.30 SESCA, and 0.47 DichroCalc. From another point of view, it is worth mentioning the remarkable capabilities of the recent approaches based on artificial intelligence, which can precisely predict the native structure of proteins. The structure of proteins, however, is flexible and can be modified by a diversity of environmental factors such as interactions with other molecules, mechanical stresses, variations of temperature, pH, or ionic strength. Experimental CD spectra together with reliable predictions can be utilized to assess eventual secondary structural changes. A similar kind of evaluation can be done for the case of an incomplete protein structure that has been reconstructed by using different approaches. The KCD method can be freely accessed from: https://kcd.cinvestav.mx/.

Template-Based Modelling of the Structure of Fungal Effector Proteins

The discovery of new fungal effector proteins is necessary to enable the screening of cultivars for disease resistance. Sequence-based bioinformatics methods have been used for this purpose, but only a limited number of functional effector proteins have been successfully predicted and subsequently validated experimentally. A significant obstacle is that many fungal effector proteins discovered so far lack sequence similarity or conserved sequence motifs. The availability of experimentally determined three-dimensional (3D) structures of a number of effector proteins has recently highlighted structural similarities amongst groups of sequence-dissimilar fungal effectors, enabling the search for similar structural folds amongst effector sequence candidates. We have applied template-based modelling to predict the 3D structures of candidate effector sequences obtained from bioinformatics predictions and the PHI-BASE database. Structural matches were found not only with ToxA- and MAX-like effector candidates but also with non-fungal effector-like proteins-including plant defensins and animal venoms-suggesting the broad conservation of ancestral structural folds amongst cytotoxic peptides from a diverse range of distant species. Accurate modelling of fungal effectors were achieved using RaptorX. The utility of predicted structures of effector proteins lies in the prediction of their interactions with plant receptors through molecular docking, which will improve the understanding of effector-plant interactions.

Effectiveness of Narciclasine in Suppressing the Inflammatory Response in Sepsis: Molecular Docking and In Silico Studies

Narciclasine is an alkaloid belonging to the Amaryllidaceae family which has been reported to have many beneficial properties. Especially its anticancer properties have been widely reported. Here, we have focused on its potential use in suppressing the inflammatory response in sepsis using in silico methods. Lipopolysaccharide (LPS) is an endotoxin which is present in the outer membrane of gram-negative bacteria and is a crucial player in the pathogenesis of gram-negative sepsis. Activation of toll-like receptor 4 (TLR4) signaling by LPS is an important event in the pathogenesis of gram-negative sepsis. This initiates a downstream signaling pathway comprising of several adaptor proteins such as toll/interleukin-1 receptor domain-containing adapter protein (TIRAP), myeloid differentiation primary response protein 88 (MyD88), interleukin-1 receptor-associated kinase (IRAK)-1, IRAK-4, interferon regulatory factor 3 (IRF-3), tumor necrosis factor receptor-associated factor 6 (TRAF-6) leading to nuclear factor kappa B (NF-κβ) activation resulting in elevated production of inflammatory cytokines such as tumor necrosis factor alpha (TNF-α) and interleukin (IL)-6. S100 calcium binding proteins A8/A9 (S100A8/A9) have been found to be an agonist of TLR4, and it amplifies the inflammatory response in sepsis. Molecular docking studies of narciclasine with target proteins associated with the LPS-TLR4 pathway showed that it has good binding affinity and stable interactions with the targets studied. Molecular dynamics (MD) simulation studies over 100 ns showed that most of the ligand-target complexes were stable. The structures of all the targets except TRAF-6 were retrieved from the Protein Data Bank (PDB) database. Homology modeling was done to predict the 3-dimensional structure of TRAF-6. MD simulation of narciclasine-TRAF-6 complex showed that the structure is stable. Metapocket was used for active site prediction in the target proteins. Toxicity analysis by admetSAR revealed that narciclasine was readily biodegradable and exhibited minimum toxicity. These results indicate that narciclasine has effective anti-inflammatory properties which could be useful in suppressing the inflammatory response in sepsis.

Advances in AI for Protein Structure Prediction: Implications for Cancer Drug Discovery and Development

Recent advancements in AI-driven technologies, particularly in protein structure prediction, are significantly reshaping the landscape of drug discovery and development. This review focuses on the question of how these technological breakthroughs, exemplified by AlphaFold2, are revolutionizing our understanding of protein structure and function changes underlying cancer and improve our approaches to counter them. By enhancing the precision and speed at which drug targets are identified and drug candidates can be designed and optimized, these technologies are streamlining the entire drug development process. We explore the use of AlphaFold2 in cancer drug development, scrutinizing its efficacy, limitations, and potential challenges. We also compare AlphaFold2 with other algorithms like ESMFold, explaining the diverse methodologies employed in this field and the practical effects of these differences for the application of specific algorithms. Additionally, we discuss the broader applications of these technologies, including the prediction of protein complex structures and the generative AI-driven design of novel proteins.

Construction of an aerolysin-based multi-epitope vaccine against Aeromonas hydrophila: an in silico machine learning and artificial intelligence-supported approach

Aeromonas hydrophila, a gram-negative coccobacillus bacterium, can cause various infections in humans, including septic arthritis, diarrhea (traveler's diarrhea), gastroenteritis, skin and wound infections, meningitis, fulminating septicemia, enterocolitis, peritonitis, and endocarditis. It frequently occurs in aquatic environments and readily contacts humans, leading to high infection rates. This bacterium has exhibited resistance to numerous commercial antibiotics, and no vaccine has yet been developed. Aiming to combat the alarmingly high infection rate, this study utilizes in silico techniques to design a multi-epitope vaccine (MEV) candidate against this bacterium based on its aerolysin toxin, which is the most toxic and highly conserved virulence factor among the Aeromonas species. After retrieval, aerolysin was processed for B-cell and T-cell epitope mapping. Once filtered for toxicity, antigenicity, allergenicity, and solubility, the chosen epitopes were combined with an adjuvant and specific linkers to create a vaccine construct. These linkers and the adjuvant enhance the MEV's ability to elicit robust immune responses. Analyses of the predicted and improved vaccine structure revealed that 75.5%, 19.8%, and 1.3% of its amino acids occupy the most favored, additional allowed, and generously allowed regions, respectively, while its ERRAT score reached nearly 70%. Docking simulations showed the MEV exhibiting the highest interaction and binding energies (-1,023.4 kcal/mol, -923.2 kcal/mol, and -988.3 kcal/mol) with TLR-4, MHC-I, and MHC-II receptors. Further molecular dynamics simulations demonstrated the docked complexes' remarkable stability and maximum interactions, i.e., uniform RMSD, fluctuated RMSF, and lowest binding net energy. In silico models also predict the vaccine will stimulate a variety of immunological pathways following administration. These analyses suggest the vaccine's efficacy in inducing robust immune responses against A. hydrophila. With high solubility and no predicted allergic responses or toxicity, it appears safe for administration in both healthy and A. hydrophila-infected individuals.

Bioinformatic analysis of human ZPR1 gene pathogenic exome mutations

Advanced sequencing technologies enable rapid detection of sequence variants, aiming to uncover the molecular foundations of human genetic disorders. The challenge lies in interpreting the influence of new exome variants that lead to diverse phenotypes. Our study introduces a detailed, multi-tiered method for assessing the impact of novel variants, particularly focusing on the zinc finger protein 1 (ZPR1) gene. Herein, we employed a combination of variant effect predictors, protein stability analyses, and the American College of Medical Genetics and Association of Molecular Pathology (ACMG/AMP) guidelines. Our structural analysis pinpoints specific amino acid residues in the ZPR1 zinc finger domains that are sensitive to changes, distinguishing between benign and disease-causing coding variants using rigorous in silico tools. We examined 223 germline ZPR1 exome variants, uncovering significant ethnic disparities in the frequency of heterozygous harmful ZPR1 variants, ranging from 0.04% in the Ashkenazi Jewish population to 0.34% in African/African Americans. Additionally, the discovery of three homozygous carriers in European and South Asian groups suggests a higher occurrence of ZPR1 variants in these demographics, meriting further exploration. This research provides insights into the prevalence and implications of amino acid substitutions in the ZPR1 protein.

Unveiling the structure and interactions of SOG1, a NAC domain transcription factor: An in-silico perspective

SOG1 is a crucial plant-specific NAC domain family transcription factor and functions as the central regulator of DNA damage response, acting downstream of ATM and ATR kinases. In this study, various in-silico approaches have been employed for the characterization of SOG1 transcription factor in a comparative manner with its orthologues from various plant species. Amino acid sequences of more than a hundred SOG1 or SOG1-like proteins were retrieved and their relationship was determined through phylogenetic and motif analyses. Various physiochemical properties and secondary structural components of SOG1 orthologues were determined in selective plant species including Arabidopsis thaliana, Oryza sativa, Amborella trichopoda, and Physcomitrella patens. Furthermore, fold recognition or threading and homology-based three-dimensional models of SOG1 were constructed followed by subsequent evaluation of quality and accuracy of the generated protein models. Finally, extensive DNA-Protein and Protein-Protein interaction studies were performed using the HADDOCK server to give an insight into the mechanism of how SOG1 binds with the promoter region of its target genes or interacts with other proteins to regulate the DNA damage responses in plants. Our docking analysis data have shown the molecular mechanism of SOG1's binding with 5'-CTT(N)7AAG-3' and 5'-(N)4GTCAA(N)4-3' consensus sequences present in the promoter region of its target genes. Moreover, SOG1 physically interacts and forms a thermodynamically stable complex with NAC103 and BRCA1 proteins, which possibly serve as coactivators or mediators in the transcription regulatory network of SOG1. Overall, our in-silico study will provide meaningful information regarding the structural and functional characterization of the SOG1 transcription factor.

Precision Nutrition Unveiled: Gene-Nutrient Interactions, Microbiota Dynamics, and Lifestyle Factors in Obesity Management

BACKGROUND

Obesity is a complex metabolic disorder that is associated with several diseases. Recently, precision nutrition (PN) has emerged as a tailored approach to provide individualised dietary recommendations.

AIM

This review discusses the major intrinsic and extrinsic components considered when applying PN during the management of obesity and common associated chronic conditions.

RESULTS

The review identified three main PN components: gene-nutrient interactions, intestinal microbiota, and lifestyle factors. Genetic makeup significantly contributes to inter-individual variations in dietary behaviours, with advanced genome sequencing and population genetics aiding in detecting gene variants associated with obesity. Additionally, PN-based host-microbiota evaluation emerges as an advanced therapeutic tool, impacting disease control and prevention. The gut microbiome's composition regulates diverse responses to nutritional recommendations. Several studies highlight PN's effectiveness in improving diet quality and enhancing adherence to physical activity among obese patients. PN is a key strategy for addressing obesity-related risk factors, encompassing dietary patterns, body weight, fat, blood lipids, glucose levels, and insulin resistance.

CONCLUSION

PN stands out as a feasible tool for effectively managing obesity, considering its ability to integrate genetic and lifestyle factors. The application of PN-based approaches not only improves current obesity conditions but also holds promise for preventing obesity and its associated complications in the long term.

NEKL-4 regulates microtubule stability and mitochondrial health in C. elegans ciliated neurons

Ciliopathies are often caused by defects in the ciliary microtubule core. Glutamylation is abundant in cilia, and its dysregulation may contribute to ciliopathies and neurodegeneration. Mutation of the deglutamylase CCP1 causes infantile-onset neurodegeneration. In C. elegans, ccpp-1 loss causes age-related ciliary degradation that is suppressed by mutation in the conserved NEK10 homolog nekl-4. NEKL-4 is absent from cilia, yet negatively regulates ciliary stability via an unknown, glutamylation-independent mechanism. We show that NEKL-4 was mitochondria-associated. nekl-4 mutants had longer mitochondria, a higher baseline mitochondrial oxidation state, and suppressed ccpp-1 mutant lifespan extension in response to oxidative stress. A kinase-dead nekl-4(KD) mutant ectopically localized to ccpp-1 cilia and rescued degenerating microtubule doublet B-tubules. A nondegradable nekl-4(PESTΔ) mutant resembled the ccpp-1 mutant with dye filling defects and B-tubule breaks. The nekl-4(PESTΔ) Dyf phenotype was suppressed by mutation in the depolymerizing kinesin-8 KLP-13/KIF19A. We conclude that NEKL-4 influences ciliary stability by activating ciliary kinesins and promoting mitochondrial homeostasis.

A viral insulin-like peptide inhibits IGF-1 receptor phosphorylation and regulates IGF1R gene expression

OBJECTIVE

The insulin/IGF superfamily is conserved across vertebrates and invertebrates. Our team has identified five viruses containing genes encoding viral insulin/IGF-1 like peptides (VILPs) closely resembling human insulin and IGF-1. This study aims to characterize the impact of Mandarin fish ranavirus (MFRV) and Lymphocystis disease virus-Sa (LCDV-Sa) VILPs on the insulin/IGF system for the first time.

METHODS

We chemically synthesized single chain (sc, IGF-1 like) and double chain (dc, insulin like) forms of MFRV and LCDV-Sa VILPs. Using cell lines overexpressing either human insulin receptor isoform A (IR-A), isoform B (IR-B) or IGF-1 receptor (IGF1R), and AML12 murine hepatocytes, we characterized receptor binding, insulin/IGF signaling. We further characterized the VILPs' effects of proliferation and IGF1R and IR gene expression, and compared them to native ligands. Additionally, we performed insulin tolerance test in CB57BL/6 J mice to examine in vivo effects of VILPs on blood glucose levels. Finally, we employed cryo-electron microscopy (cryoEM) to analyze the structure of scMFRV-VILP in complex with the IGF1R ectodomain.

RESULTS

VILPs can bind to human IR and IGF1R, stimulate receptor autophosphorylation and downstream signaling pathways. Notably, scMFRV-VILP exhibited a particularly strong affinity for IGF1R, with a mere 10-fold decrease compared to human IGF-1. At high concentrations, scMFRV-VILP selectively reduced IGF-1 stimulated IGF1R autophosphorylation and Erk phosphorylation (Ras/MAPK pathway), while leaving Akt phosphorylation (PI3K/Akt pathway) unaffected, indicating a potential biased inhibitory function. Prolonged exposure to MFRV-VILP led to a significant decrease in IGF1R gene expression in IGF1R overexpressing cells and AML12 hepatocytes. Furthermore, insulin tolerance test revealed scMFRV-VILP's sustained glucose-lowering effect compared to insulin and IGF-1. Finally, cryo-EM analysis revealed that scMFRV-VILP engages with IGF1R in a manner closely resembling IGF-1 binding, resulting in a highly analogous structure.

CONCLUSIONS

This study introduces MFRV and LCDV-Sa VILPs as novel members of the insulin/IGF superfamily. Particularly, scMFRV-VILP exhibits a biased inhibitory effect on IGF1R signaling at high concentrations, selectively inhibiting IGF-1 stimulated IGF1R autophosphorylation and Erk phosphorylation, without affecting Akt phosphorylation. In addition, MFRV-VILP specifically regulates IGF-1R gene expression and IGF1R protein levels without affecting IR. CryoEM analysis confirms that scMFRV-VILP' binding to IGF1R is mirroring the interaction pattern observed with IGF-1. These findings offer valuable insights into IGF1R action and inhibition, suggesting potential applications in development of IGF1R specific inhibitors and advancing long-lasting insulins.

Investigation and drug design for novel molecules from natural products as inhibitors for controlling multiple myeloma disease using in-silico tools

Multiple myeloma (MM) is a disease that causes plasma cell growth in the bone marrow and immune globulin buildup in blood and urine. Despite recent advances in MM therapy, many still die due to its high mortality rate. A study using computational simulations analyzed 100 natural ingredients from the SANC database to determine if they inhibited the IgH domain, a known cause of multiple myeloma. Natural component Diospyrin inhibited the IgH enzyme with the best binding energy of -10.3 kcal/mol and three carbon-hydrogen bonds, followed by Parviflorone F complex with a binding energy of -10.1 kcal/mol and two conventional-hydrogen bonds. As a result, the Molecular Dynamic simulation was used to test the stability of the two complexes. During the simulation, the Diospyrin molecule dissociated from the protein at roughly 67.5 ns, whereas the Parviflorone F molecule stayed attached to the protein throughout. The latter was the subject of the investigation. The analysis of the production run data revealed that the Parviflorone F molecule exhibits a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The analysis of the production run data revealed that the Parviflorone F molecule exhibited a variety of conformations within the binding pocket while keeping a relatively constant distance from the protein's center of mass. The root mean square deviation (RMSD) plots for both the protein and complex showed a stable and steady average value of 4.4 Å for the first 82 nanoseconds of manufacture. As a result, the average value increased to 8.3 Å. Furthermore, the components of the binding free energy, as computed by MM-GBSA, revealed that the mean binding energy of the Parviflorone F molecule was -23.88 kcal/mol. Finally, after analyzing all of the examination data, Parviflorone F was identified as a powerful inhibitor of the IgH domain and hence of the MM disease, which requires further in-vivo conformation.Communicated by Ramaswamy H. Sarma.

The Application of MD Simulation to Lead Identification, Vaccine Design, and Structural Studies in Combat against Leishmaniasis - A Review

Drug discovery, vaccine design, and protein interaction studies are rapidly moving toward the routine use of molecular dynamics simulations (MDS) and related methods. As a result of MDS, it is possible to gain insights into the dynamics and function of identified drug targets, antibody-antigen interactions, potential vaccine candidates, intrinsically disordered proteins, and essential proteins. The MDS appears to be used in all possible ways in combating diseases such as cancer, however, it has not been well documented as to how effectively it is applied to infectious diseases such as Leishmaniasis. As a result, this review aims to survey the application of MDS in combating leishmaniasis. We have systematically collected articles that illustrate the implementation of MDS in drug discovery, vaccine development, and structural studies related to Leishmaniasis. Of all the articles reviewed, we identified that only a limited number of studies focused on the development of vaccines against Leishmaniasis through MDS. Also, the PCA and FEL studies were not carried out in most of the studies. These two were globally accepted utilities to understand the conformational changes and hence it is recommended that this analysis should be taken up in similar approaches in the future.

Immunoinformatic Analysis of Leishmania Major gp46 Protein and Potential Targets for Vaccination against Leishmaniasis

BACKGROUND

Cutaneous leishmaniasis (CL) is a parasitic disease with a significant burden in the Old World countries.

OBJECTIVE

In the current study, some of the primary biochemical properties and IFN-γ inducing epitopes with specific binding capacity to human and mouse MHC alleles were predicted for Leishmania major gp46 antigenic protein.

METHODS

Several online servers were used to predict physico-chemical traits, allergenicity, antigenicity, transmembrane domain and signal peptide, subcellular localization, post-translational modifications (PTMs), secondary and tertiary structures, tertiary model refining with validations. Also, IEDB web server was used to predict mouse/human cytotoxic T-lymphocyte (CTL) and helper T-lymphocyte (HTL) epitopes.

RESULTS

The 33.25 kDa protein was stable, hydrophilic, antigenic, while non-allergenic, with enhanced thermotolerance and 45 PTM sites. The secondary structure encompassed a random coil, followed by extended strands and helices. Ramachandran-based analysis of the refined model showed 73.1%, 21.6%, 3.4% and 1.9% of residues in the most favored, additional allowed, generously-allowed and disallowed regions, respectively. Epitope screening demonstrated 4 HTL epitopes against seemingly protective HLA alleles, 5 HTL epitopes against the HLA reference set, 3 human CTL epitopes and a number of mouse MHC-restricted epitopes.

CONCLUSION

This paper provides insights into the bioinformatics characteristics of the L. major gp46 protein as a promising vaccine candidate.

Antimicrobial peptides grafted onto the surface of N-acetylcysteine-chitosan nanoparticles can revitalize drugs against clinical isolates of Mycobacterium tuberculosis

Tuberculosis is caused by Mycobacterium tuberculosis (MTB) and is the leading cause of death from infectious diseases in the World. The search for new antituberculosis drugs is a high priority, since several drug-resistant TB-strains have emerged. Many nanotechnology strategies are being explored to repurpose or revive drugs. An interesting approach is to graft antimicrobial peptides (AMPs) to antibiotic-loaded nanoparticles. The objective of the present work was to determine the anti-MTB activity of rifampicin-loaded N-acetylcysteine-chitosan-based nanoparticles (NPs), conjugated with the AMP Ctx(Ile21)-Ha; against clinical isolates (multi- and extensively-drug resistant) and the H37Rv strain. The modified chitosan and drug-loaded NPs were characterized with respect to their physicochemical stability and their antimycobacterial profile, which showed potent inhibition (MIC values <0.977 μg/mL) by the latter. Furthermore, their accumulation within macrophages and cytotoxicity were determined. To understand the possible mechanisms of action, an in silico study of the peptide against MTB membrane receptors was performed. The results presented herein demonstrate that antibiotic-loaded NPs grafted with an AMP can be a powerful tool for revitalizing drugs against multidrug-resistant M. tuberculosis strains, by launching multiple attacks against MTB. This approach could potentially serve as a novel treatment strategy for various long-term diseases requiring extended treatment periods.

Transglutaminase 2 has higher affinity for relaxed than for stretched fibronectin fibers

Transglutaminase 2 (TG2) plays a vital role in stabilizing extracellular matrix (ECM) proteins through enzymatic crosslinking during tissue growth, repair, and inflammation. TG2 also binds non-covalently to fibronectin (FN), an essential component of the ECM, facilitating cell adhesion, migration, proliferation, and survival. However, the interaction between TG2 and fibrillar FN remains poorly understood, as most studies have focused on soluble or surface-adsorbed FN or FN fragments, which differ in their conformations from insoluble FN fibers. Using a well-established in vitro FN fiber stretch assay, we discovered that the binding of a crosslinking enzyme to ECM fibers is mechano-regulated. TG2 binding to FN is tuned by the mechanical tension of FN fibers, whereby TG2 predominantly co-localizes to low-tension FN fibers, while fiber stretching reduces their affinity for TG2. This mechano-regulated binding relies on the proximity between the N-terminal β-sandwich and C-terminal β-barrels of TG2. Crosslinking mass spectrometry (XL-MS) revealed a novel TG2-FN synergy site within TG2's C-terminal β-barrels that interacts with FN regions located outside of the canonical gelatin binding domain, specifically FNI2 and FNIII14-15. Combining XL-MS distance restraints with molecular docking revealed the mechano-regulated binding mechanism between TG2 and modules FNI7-9 by which mechanical forces regulate TG2-FN interactions. This highlights a previously unrecognized role of TG2 as a tension sensor for FN fibers. This novel interaction mechanism has significant implications in physiology and mechanobiology, including how forces regulate cell adhesion, spreading, migration, phenotype modulation, depending on the tensional state of ECM fibers. Data are available via ProteomeXchange with identifier PXD043976.

Unmasking crucial residues in adipose triglyceride lipase for coactivation with comparative gene identification-58

Lipolysis is an essential metabolic process that releases unesterified fatty acids from neutral lipid stores to maintain energy homeostasis in living organisms. Adipose triglyceride lipase (ATGL) plays a key role in intracellular lipolysis and can be coactivated upon interaction with the protein comparative gene identification-58 (CGI-58). The underlying molecular mechanism of ATGL stimulation by CGI-58 is incompletely understood. Based on analysis of evolutionary conservation, we used site directed mutagenesis to study a C-terminally truncated variant and full-length mouse ATGL providing insights in the protein coactivation on a per-residue level. We identified the region from residues N209-N215 in ATGL as essential for coactivation by CGI-58. ATGL variants with amino acids exchanges in this region were still able to hydrolyze triacylglycerol at the basal level and to interact with CGI-58, yet could not be activated by CGI-58. Our studies also demonstrate that full-length mouse ATGL showed higher tolerance to specific single amino acid exchanges in the N209-N215 region upon CGI-58 coactivation compared to C-terminally truncated ATGL variants. The region is either directly involved in protein-protein interaction or essential for conformational changes required in the coactivation process. Three-dimensional models of the ATGL/CGI-58 complex with the artificial intelligence software AlphaFold demonstrated that a large surface area is involved in the protein-protein interaction. Mapping important amino acids for coactivation of both proteins, ATGL and CGI-58, onto the 3D model of the complex locates these essential amino acids at the predicted ATGL/CGI-58 interface thus strongly corroborating the significance of these residues in CGI-58-mediated coactivation of ATGL.

Characterization of terminal flowering cowpea (Vigna unguiculata (L.) Walp.) mutants obtained by induced mutagenesis digs out the loss-of-function of phosphatidylethanolamine-binding protein

Cowpea (Vigna unguiculata (L.) Walp) is one of the major food legume crops grown extensively in arid and semi-arid regions of the world. The determinate habit of cowpea has many advantages over the indeterminate and is well adapted to modern farming systems. Mutation breeding is an active research area to develop the determinate habit of cowpea. The present study aimed to develop new determinate habit mutants with terminal flowering (TFL) in locally well-adapted genetic backgrounds. Consequently, the seeds of popular cowpea cv P152 were irradiated with doses of gamma rays (200, 250, and, 300 Gy), and the M1 populations were grown. The M2 populations were produced from the M1 progenies and selected determinate mutants (TFLCM-1 and TFLCM-2) from the M2 generation (200 Gy) were forwarded up to the M5 generation to characterize the mutants and simultaneously they were crossed with P152 to develop a MutMap population. In the M5 generation, determinate mutants (80-81 days) were characterized by evaluating the TFL growth habit, longer peduncles (30.75-31.45 cm), erect pods (160°- 200°), number of pods per cluster (4-5 nos.), and early maturity. Further, sequencing analysis of the VuTFL1 gene in the determinate mutants and MutMap population revealed a single nucleotide transversion (A-T at 1196 bp) in the fourth exon and asparagine (N) to tyrosine (Y) amino acid change at the 143rd position of phosphatidylethanolamine-binding protein (PEBP). Notably, the loss of function PEPB with a higher confidence level modification of anti-parallel beta-sheets and destabilization of the protein secondary structure was observed in the mutant lines. Quantitative real-time PCR (qRT-PCR) analysis showed that the VuTFL1 gene was downregulated at the flowering stage in TFL mutants. Collectively, the insights garnered from this study affirm the effectiveness of induced mutation in modifying the plant's ideotype. The TFL mutants developed during this investigation have the potential to serve as a valuable resource for fostering determinate traits in future cowpea breeding programs and pave the way for mechanical harvesting.

Design of a novel multi-epitopes vaccine against Escherichia fergusonii: a pan-proteome based in- silico approach

Escherichia fergusonii a gram-negative rod-shaped bacterium in the Enterobacteriaceae family, infect humans, causing serious illnesses such as urinary tract infection, cystitis, biliary tract infection, pneumonia, meningitis, hemolytic uremic syndrome, and death. Initially treatable with penicillin, antibiotic misuse led to evolving resistance, including resistance to colistin, a last-resort drug. With no licensed vaccine, the study aimed to design a multi-epitope vaccine against E. fergusonii. The study started with the retrieval of the complete proteome of all known strains and proceeded to filter the surface exposed virulent proteins. Seventeen virulent proteins (4 extracellular, 4 outer membranes, 9 periplasmic) with desirable physicochemical properties were identified from the complete proteome of known strains. Further, these proteins were processed for B-cell and T-cell epitope mapping. Obtained epitopes were evaluated for antigenicity, allergenicity, solubility, MHC-binding, and toxicity and the filtered epitopes were fused by specific linkers and an adjuvant into a vaccine construct. Structure of the vaccine candidate was predicted and refined resulting in 78.1% amino acids in allowed regions and VERIFY3D score of 81%. Vaccine construct was docked with TLR-4, MHC-I, and MHC-II, showing binding energies of -1040.8 kcal/mol, -871.4 kcal/mol, and -1154.6 kcal/mol and maximum interactions. Further, molecular dynamic simulation of the docked complexes was carried out resulting in a significant stable nature of the docked complexes (high B-factor and deformability values, lower Eigen and high variance values) in terms of intermolecular binding conformation and interactions. The vaccine was also reported to stimulate a variety of immunological pathways after administration. In short, the designed vaccine revealed promising predictions about its immune protective potential against E. fergusonii infections however experimental validation is needed to validate the results.

Computer-Aided Multi-Epitope Based Vaccine Design Against Monkeypox Virus Surface Protein A30L: An Immunoinformatics Approach

Monkeypox, a viral zoonotic disease resembling smallpox, has emerged as a significant national epidemic primarily in Africa. Nevertheless, the recent global dissemination of this pathogen has engendered apprehension regarding its capacity to metamorphose into a sweeping pandemic. To effectively combat this menace, a multi-epitope vaccine has been meticulously engineered with the specific aim of targeting the cell envelope protein of Monkeypox virus (MPXV), thereby stimulating a potent immunological response while mitigating untoward effects. This new vaccine uses T-cell and B-cell epitopes from a highly antigenic, non-allergenic, non-toxic, conserved, and non-homologous A30L protein to provide protection against the virus. In order to ascertain the vaccine design with the utmost efficacy, protein-protein docking methodologies were employed to anticipate the intricate interactions with Toll-like receptors (TLR) 2, 3, 4, 6, and 8. This meticulous approach led the researchers to discern an optimal vaccine architecture, bolstered by affirmative prognostications derived from both molecular dynamics (MD) simulations and immune simulations. The current research findings indicate that the peptides ATHAAFEYSK, FFIVVATAAV, and MNSLSIFFV exhibited antigenic properties and were determined to be non-allergenic and non-toxic. Through the utilization of codon optimization and in-silico cloning techniques, our investigation revealed that the prospective vaccine exhibited a remarkable expression level within Escherichia coli. Moreover, upon conducting immune simulations, we observed the induction of a robust immune response characterized by elevated levels of both B-cell and T-cell mediated immunity. Moreover, as the initial prediction with in-silico techniques has yielded promising results these epitope-based vaccines can be recommended to in vitro and in silico studies to validate their immunogenic properties.

In silico identification and characterization of potential druggable targets among hypothetical proteins of Leptospira interrogans serovar Copenhageni: a comprehensive bioinformatics approach

Leptospirosis is one of the neglected zoonosis, affecting human and animal populations worldwide. Reliable effective therapeutics and concerns to look for more research into the molecular analysis of its genome is therefore needed. In the genomic pool of the Leptospira interrogans many hypothetical proteins are still uncharacterized. In the current research, we performed extensive in silico analysis to prioritize the potential hypothetical proteins of L. interrogans serovar Copenhageni via stepwise reducing the available hypothetical proteins (Total 3606) of the assembly to only 15, based on non-homologous to homosapien, essential, functional, virulent, cellular localization. Out of them, only two proteins WP_000898918.1 (Hypothetical Protein 1) & WP_001014594.1 (Hypothetical Protein 2) were found druggable and involved in protein-protein interaction network. The 3 D structures of these two target proteins were predicted via ab initio homology modeling followed by structures refinement and validation, as no structures were available till date. The analysis also revealed that the functional domains, families and protein-protein interacting partners identified in both proteins are crucial for the survival of the bacteria. The binding cavities were predicted for both the proteins through blind and specific protein-ligand docking with their respective ligands and inhibitors and were found to be in accordance with the druggable sites predicted by DoGSiteScorer. The docking interactions were found within the active functional domains for both the proteins while for Hypothetical Protein 2, the same residues were involved in interactions with Cytidine-5'-triphosphate in blind and specific docking. Furthermore, the simulations of molecular dynamics and free binding energy revealed the stable substrate binding and efficient binding energies, and were in accordance to our docking results. The work predicted two unique hypothetical proteins of L. interrogans as a potential druggable targets for designing of inhibitors for them.Communicated by Ramaswamy H. Sarma.

Impacts of Hydrophobic Mismatch on Antimicrobial Peptide Efficacy and Bilayer Permeabilization

Antimicrobial resistance continues to be a major threat to world health, with the continued emergence of resistant bacterial strains. Antimicrobial peptides have emerged as an attractive option for the development of novel antimicrobial compounds in part due to their ubiquity in nature and the general lack of resistance development to this class of molecules. In this work, we analyzed the antimicrobial peptide C18G and several truncated forms for efficacy and the underlying mechanistic effects of the sequence truncation. The peptides were screened for antimicrobial efficacy against several standard laboratory strains, and further analyzed using fluorescence spectroscopy to evaluate binding to model lipid membranes and bilayer disruption. The results show a clear correlation between the length of the peptide and the antimicrobial efficacy. Furthermore, there is a correlation between peptide length and the hydrophobic thickness of the bilayer, indicating that hydrophobic mismatch is likely a contributing factor to the loss of efficacy in shorter peptides.

Genome sequencing of cucumber mosaic virus (CMV) isolates infecting chilli and its interaction with host ferredoxin protein of different host for causing mosaic symptoms

Chilli (Capsicum annuum L.) is an important vegetable crop grown in the Indian sub-continent and is prone to viral infections under field conditions. During the field survey, leaf samples from chilli plants showing typical symptoms of disease caused by cucumber mosaic virus (CMV) such as mild mosaic, mottling and leaf distortion were collected. DAC-ELISA analysis confirmed the presence of CMV in 71 out of 100 samples, indicating its widespread prevalence in the region. Five CMV isolates, named Gu1, Gu2, BA, Ho, and Sal were mechanically inoculated onto cucumber and Nicotiana glutinosa plants to study their virulence. Inoculated plants expressed the characteristic symptoms of CMV such as chlorotic spots followed by mild mosaic and leaf distortion. Complete genomes of the five CMV isolates were amplified, cloned, and sequenced, revealing RNA1, RNA2, and RNA3 sequences with 3358, 3045, and 2220 nucleotides, respectively. Phylogenetic analysis classified the isolates as belonging to the CMV-IB subgroup, distinguishing them from subgroup IA and II CMV isolates. Recombination analysis showed intra and interspecific recombination in all the three RNA segments of these isolates. In silico protein-protein docking approach was used to decipher the mechanism behind the production of mosaic symptoms during the CMV-host interaction in 13 host plants. Analysis revealed that the production of mosaic symptoms could be due to the interaction between the coat protein (CP) of CMV and chloroplast ferredoxin proteins. Further, in silico prediction was validated in 13 host plants of CMV by mechanical sap inoculation. Twelve host plants produced systemic symptoms viz., chlorotic spot, chlorotic ringspot, chlorotic local lesion, mosaic and mild mosaic and one host plant, Solanum lycopersicum produced mosaic followed by shoestring symptoms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03777-8.

Genome-wide identification and expression analysis of Na+/H+antiporter (NHX) genes in tomato under salt stress

Plant Na +/H + antiporter (NHX) genes enhance salt tolerance by preventing excessive Na+ accumulation in the cytosol through partitioning of Na+ ions into vacuoles or extracellular transport across the plasma membrane. However, there is limited detailed information regarding the salt stress responsive SlNHXs in the most recent tomato genome. We investigated the role of this gene family's expression patterns in the open flower tissues under salt shock in Solanum lycopersicum using a genome-wide approach. A total of seven putative SlNHX genes located on chromosomes 1, 4, 6, and 10 were identified, but no ortholog of the NHX5 gene was identified in the tomato genome. Phylogenetic analysis revealed that these genes are divided into three different groups. SlNHX proteins with 10-12 transmembrane domains were hypothetically localized in vacuoles or cell membranes. Promoter analysis revealed that SlNHX6 and SlNHX8 are involved with the stress-related MeJA hormone in response to salt stress signaling. The structural motif analysis of SlNHX1, -2, -3, -4, and -6 proteins showed that they have highly conserved amiloride binding sites. The protein-protein network revealed that SlNHX7 and SlNHX8 interact physically with Salt Overly Sensitive (SOS) pathway proteins. Transcriptome analysis demonstrated that the SlNHX2 and SlNHX6 genes were substantially expressed in the open flower tissues. Moreover, quantitative PCR analysis indicated that all SlNHX genes, particularly SlNHX6 and SlNHX8, are significantly upregulated by salt shock in the open flower tissues. Our results provide an updated framework for future genetic research and development of breeding strategies against salt stress in the tomato.