ExPASy: The proteomics server for in-depth protein knowledge and analysis

Insights into the Unusual Activity of a Novel Homospermidine Synthase with a Promising Application to Produce Spermidine

Spermidine is a naturally occurring polyamine with multiple biological activities and potential food and agricultural applications. However, sustainable and scalable spermidine production has not yet been attained. In this study, a homospermidine synthase (HSS) from Pseudomonas frederiksbergensis (PfHSS) capable of catalyzing the synthesis of spermidine from 1,3-diaminopropane and putrescine was identified based on multiple sequence alignment using Blastochloris viridis HSS (BvHSS) as a template. The optimal reaction pH and temperature for purified PfHSS were determined to be 8.5 and 45 °C, respectively, and K+ was able to promote the enzyme activity. Further analysis of the structural and functional relationships through molecular docking and molecular dynamics simulation indicates that glutamic acid at position 359 is the essential residue for the enzyme-catalyzed synthesis of spermidine. The whole-cell catalytic reaction yielded 1321.4 mg/L spermidine and 678.2 mg/L of homospermidine. This study presents a novel, promising, and sustainable biological method for producing spermidine.

Genome-Wide Identification of Bilberry WRKY Transcription Factors: Go Wild and Duplicate

WRKY transcription factor genes compose an important family of transcriptional regulators that are present in several plant species. According to previous studies, these genes can also perform important roles in bilberry (Vaccinium myrtillus L.) metabolism, making it essential to deepen our understanding of fruit ripening regulation and anthocyanin biosynthesis. In this context, the detailed characterization of these proteins will provide a comprehensive view of the functional features of VmWRKY genes in different plant organs and in response to different intensities of light. In this study, the investigation of the complete genome of the bilberry identified 76 VmWRKY genes that were evaluated and distributed in all twelve chromosomes. The proteins encoded by these genes were classified into four groups (I, II, III, and IV) based on their conserved domains and zinc finger domain types. Fifteen pairs of VmWRKY genes in segmental duplication and four pairs in tandem duplication were detected. A cis element analysis showed that all promoters of the VmWRKY genes contain at least one potential cis stress-response element. Differential expression analysis of RNA-seq data revealed that VmWRKY genes from bilberry show preferential or specific expression in samples. These findings provide an overview of the functional characterization of these proteins in bilberry.

Identification, characterization, and agglutinating activity of a novel C-type lectin domain family 3 member B (CLEC3B) discovered in golden pompano, Trachinotus ovatus

The lectins are a large family of carbohydrate-binding proteins that play important roles in the innate immune response of various organisms. Although C-type lectin domain family 3 member B (CLEC3B), an important member of C-type lectin, has been well documented in humans and several other higher vertebrates, little is currently known about this molecule in economically important marine fish species. In this study, through transcriptomic and BLAST screening, a novel CLEC3B gene was identified in the golden pompano (Trachinotus ovatus). The T. ovatus CLEC3B (ToCLEC3B) was subsequently characterized by bioinformatic analysis and compared with those reported in other species. In addition, the expression patterns of ToCLEC3B in different tissues under normal condition and at different times post pathogen challenge were assessed. Furthermore, the agglutinating activity of ToCLEC3B with and without Ca2+ against different bacteria and blood cells of donor species were verified using the recombinant T. ovatus CLEC3B (rToCLEC3B). Our results demonstrated that ToCLEC3B is a Ca2+-dependent galactose-binding lectin with a single copy of carbohydrate recognition domain (CRD). Similar to CLEC3B reported in other species, the CRD domain of ToCLEC3B consists of two α-helices, six β-sheets, and four loops, forming two Ca2+- and a galactose-binding sites. According to the phylogenetic analysis, the ToCLEC3B was highly similar (similarity at 95.00%) to that of its relative, the greater amberjack (Seriola dumerili). The expression of ToCLEC3B was detected in all tissues examined under normal condition and was significantly up-regulated by injection of pathogenic microbes. In addition, the rToCLEC3B exhibited strong agglutinating activity against different bacteria and blood cells of donor species in the presence of Ca2+. Our results indicate that ToCLEC3B is a constitutive and inducible acute-phase immune factor in the host's innate immune response of T. ovatus.

Identification and expression analysis of the SWEET genes in radish reveal their potential functions in reproductive organ development

BACKGROUND

Sugars produced by photosynthesis provide energy for biological activities and the skeletons for macromolecules; they also perform multiple physiological functions in plants. Sugar transport across plasma membranes mediated by the Sugar Will Eventually be Exported Transporter (SWEET) genes substantially affects these processes. However, the evolutionary dynamics and function of the SWEET genes are largely unknown in radish, an important Brassicaceae species.

METHODS AND RESULTS

Genome-wide identification and analysis of the RsSWEET genes from the recently updated radish reference genome was conducted using bioinformatics methods. The tissue-specific expression was analyzed using public RNA-seq data, and the expression levels in the bud, stamens, pistils, pericarps and seeds at 15 and 30 days after flowering (DAF) were determined by RT‒qPCR. Thirty-seven RsSWEET genes were identified and named according to their Arabidopsis homologous. They are unevenly distributed across the nine radish chromosomes and were further divided into four clades by phylogenetic analysis. There are 5-7 transmembrane domains and at least one MtN3_slv domain in the RsSWEETs. RNA-seq and RT‒qPCR revealed that the RsSWEETs exhibit higher expression levels in the reproductive organs, indicating that these genes might play vital roles in reproductive organ development. RsSWEET15.1 was found to be especially expressed in siliques according to the RNA-seq data, and the RT‒qPCR results further confirmed that it was most highly expressed levels in the seeds at 30 DAF, followed by the pericarp at 15 DAF, indicating that it is involved in seed growth and development.

CONCLUSIONS

This study suggests that the RsSWEET genes play vital roles in reproductive organ development and provides a theoretical basis for the future functional analysis of RsSWEETs in radish.

Comparison of IGLV2-14 light chain sequences of patients with AL amyloidosis or multiple myeloma

Light chain amyloidosis (AL) is one of the most common forms of systemic amyloidosis and is caused by the deposition of insoluble fibrils derived from misfolded and aggregated immunoglobulin light chains (LC). To uncover the causes leading to this aggregation, we compared AL LC sequences with those of patients with the related disease multiple myeloma (MM), which do not aggregate in insoluble fibrils in vivo. IGLV2-14 is one of the most common AL-associated IGLV subfamilies. Here, we analysed IGLV2-14 LC sequences of 13 AL and eight MM patients in detail. We found that AL-associated LCs presented a lower median mutation count (7.0 vs. 11.5 in MM; P = 0.045), as well as an overall composition of less charged amino acids than MM LCs. However, we did not find a mutation that was present in ≥ 50% of the AL and not in the MM sequences. Furthermore, we did not find a significant difference in the isoelectric point (pI) in general, suggesting similar stability of the LCs in AL and MM. However, the subgroup of patients without a detectable heavy chain stood out. Surprisingly, they are characterized by an increase in mutation count (median 7.0 vs. 5.5) and pI (median 7.82 vs. 6.44, P = 0.043). In conclusion, our data suggest that the amount of mutations and the introduction of charges play a crucial role in AL fibril formation, as well as the absence or presence of a potential heavy chain binding partner.

Genome-wide characterization of integrin (ITG) gene family and their expression profiling in half-smooth tongue sole (Cynoglossus semilaevis) upon Vibrio anguillarum infection

Integrins (ITGs) are transmembrane heterodimer receptors with ITGα subunit and ITGβ subunit, participating in various physiological processes, including immunity. At present, systematic research on ITGs in teleost is scarce, especially in half-smooth tongue sole (Cynoglossus semilaevis). In this study, a set of 28 ITG genes in half-smooth tongue sole have been identified and characterized. The phylogenetic analysis showed that ITGα and ITGβ subunits were respectively classified into five and two clusters, consistent with previous studies. The selection pressure analysis indicated that most of ITG genes were under purifying selection, except for ITGα11b and ITGαL with positive selection. The expression profiles of eight selected ITG genes, including ITGα1, ITGα5, ITGα8, ITGα11, ITGβ1, ITGβ2, ITGβ3, and ITGβ8, were analyzed in healthy tissues and after infection with Vibrio anguillarum, revealed their implications in immune response. The study provided a comprehensive characterization and expression analysis of ITG genes in half-smooth tongue sole, setting a solid foundation for further functional studies and promising potential in disease control.

Genome-wide identification of trihelix transcription factors in the apple genome in silico

Trihelix transcription factors are involved in the growth and development of plants, as well as various stress responses. In the study presented, we identified 37 trihelix family genes in the apple genome (MdTH). The trihelix genes were located on 13 chromosomes. Phylogenetic analysis of these MdTH and the trihelix genes of other species divided them into six subfamilies: GT-1, GT-2, SH4, SIP1, GTγ, and GTδ. The genes of different groups significantly diverged in their gene structure and conserved functional domains. Cis-element analysis showed that promoter sequences of MdTH genes contained light response elements, phytohormone response elements, and stress-related cis-elements. The expression pattern analysis results demonstrated that MdTH were regulated by drought, salinity, as well as high and low temperatures. MdTH4 and MdTH24 were highly regulated by soil salinity, MdTH4-by drought. MdTH30 showed high expression under low temperature; MdTH8, MdTH20, and MdTH36-under high temperature.

Systematic evolution of bZIP transcription factors in Malvales and functional exploration of AsbZIP14 and AsbZIP41 in Aquilaria sinensis

Introduction

Agarwood, the dark-brown resin produced by Aquilaria trees, has been widely used as incense, spice, perfume or traditional medicine and 2-(2-phenethyl) chromones (PECs) are the key markers responsible for agarwood formation. But the biosynthesis and regulatory mechanism of PECs were still not illuminated. The transcription factor of basic leucine zipper (bZIP) presented the pivotal regulatory roles in various secondary metabolites biosynthesis in plants, which might also contribute to regulate PECs biosynthesis. However, molecular evolution and function of bZIP are rarely reported in Malvales plants, especially in Aquilaria trees.

Methods and results

Here, 1,150 bZIPs were comprehensively identified from twelve Malvales and model species genomes and the evolutionary process were subsequently analyzed. Duplication types and collinearity indicated that bZIP is an ancient or conserved TF family and recent whole genome duplication drove its evolution. Interesting is that fewer bZIPs in A. sinensis than that species also experienced two genome duplication events in Malvales. 62 AsbZIPs were divided into 13 subfamilies and gene structures, conservative domains, motifs, cis-elements, and nearby genes of AsbZIPs were further characterized. Seven AsbZIPs in subfamily D were significantly regulated by ethylene and agarwood inducer. As the typical representation of subfamily D, AsbZIP14 and AsbZIP41 were localized in nuclear and potentially regulated PECs biosynthesis by activating or suppressing type III polyketide synthases (PKSs) genes expression via interaction with the AsPKS promoters.

Discussion

Our results provide a basis for molecular evolution of bZIP gene family in Malvales and facilitate the understanding the potential functions of AsbZIP in regulating 2-(2-phenethyl) chromone biosynthesis and agarwood formation.

Design of an alternate antibody fragment format that can be produced in the cytoplasm of Escherichia coli

With increased accessibility and tissue penetration, smaller antibody formats such as antibody fragments (Fab) and single chain variable fragments (scFv) show potential as effective and low-cost choices to full-length antibodies. These formats derived from the modular architecture of antibodies could prove to be game changers for certain therapeutic and diagnostic applications. Microbial hosts have shown tremendous promise as production hosts for antibody fragment formats. However, low target protein yields coupled with the complexity of protein folding result in production limitations. Here, we report an alternative antibody fragment format 'FabH3' designed to overcome some key bottlenecks associated with the folding and production of Fabs. The FabH3 molecule is based on the Fab format with the constant domains replaced by engineered immunoglobulin G1 (IgG1) CH3 domains capable of heterodimerization based on the electrostatic steering approach. We show that this alternative antibody fragment format can be efficiently produced in the cytoplasm of E. coli using the catalyzed disulfide-bond formation system (CyDisCo) in a natively folded state with higher soluble yields than its Fab counterpart and a comparable binding affinity against the target antigen.

Predicting the changes in neutralizing antibody interaction with G protein derived from Bangladesh isolates of Nipah virus: molecular dynamics based approach

The infectious Nipah virus (NiV) is categorized into NiV-M (Malaysia) and NiV-B (Bangladesh) groups based on its genome comparison, pathogenicity, and mortality rate. The development of therapeutic molecules has used NiV-M-derived data in multiple studies than NiV-B. In continuation with this, the protein level investigation is also less explored to understand the interaction with therapeutic neutralizing antibodies for NiV-B. So, this study focuses on understanding the impact of NiV-B-specific mutations on the interaction of therapeutic neutralizing antibodies with the G protein. The population-based comparative analysis of NiV-B G protein sequences with NiV-M sequence identified twenty-six mutations. These predominantly polar mutations were then used to model the mutant protein (G_MT). In a comparative study, the G protein G_MT and reference protein G_WT (Malaysian origin) were subjected to a protein docking with neutralizing human monoclonal antibody HENV26. The binding affinity and the free binding energy of the glycoprotein in complex with G-WT and G_MT were calculated using PRODIGY and MM/PBSA tools respectively. Based on the PRODIGY report, G-WT showed stronger binding (-13.8 kcal/mol) compared to that of the G_MT (-9.0 kcal/mol) with the HENV26 antibody. The stability of the complexes was evaluated using MM/PBSA which showed higher binding energy with HENV26 for G_WT (-75.11 kcal/mol) in contrast to G_MT (-41.66 kcal/mol). The results indicate that the mutant G protein has a reduced ability to bind to neutralizing antibodies, resulting in a decreased effectiveness against strains carrying these mutations.Communicated by Ramaswamy H. Sarma.

Identification and characterization of a novel chromosomal aminoglycoside 3'-O-phosphotransferase, APH(3')-Id, from Kluyvera intermedia DW18 isolated from the sewage of an animal farm

Background

Aminoglycosides, as important clinical antimicrobials, are used as second-line drugs for treating multidrug-resistant tuberculosis or combined with β-lactam drugs for treating severe infections such as sepsis. Aminoglycoside-modifying enzyme (AME) is the most important mechanism of aminoglycoside resistance and deserves more attention.

Methods

The bacterium Kluyvera intermedia DW18 was isolated from the sewage of an animal farm using the conventional method. The agar dilution method was used to determine the minimum inhibitory concentrations (MICs) of antimicrobials. A novel resistance gene was cloned, and the enzyme was expressed. The kinetic parameters were measured by a SpectraMax M5 multifunctional microplate reader. Bioinformatic analysis was performed to reveal the genetic context of the aph(3')-Id gene and its phylogenetic relationship with other AMEs.

Results

A novel aminoglycoside 3'-O-phosphotransferase gene designated aph(3')-Id was identified in K. intermedia DW18 and shared the highest amino acid identity of 77.49% with the functionally characterized aminoglycoside 3'-O-phosphotransferase APH(3')-Ia. The recombinant plasmid carrying the novel resistance gene (pMD19-aph(3')-Id/E. coli DH5α) showed 1,024-, 512-, 128- and 16-fold increased MIC levels for kanamycin, ribostamycin, paromomycin and neomycin, respectively, compared with the reference strain DH5α. APH(3')-Id showed the highest catalytic efficiency for ribostamycin [kcat/Km of (4.96 ± 1.63) × 105 M-1/s-1], followed by paromomycin [kcat/Km of (2.18 ± 0.21) × 105 M-1/s-1], neomycin [kcat/Km of (1.73 ± 0.20) × 105 M-1/s-1], and kanamycin [kcat/Km of (1.10 ± 0.18) × 105 M-1/s-1]. Three conserved functional domains of the aminoglycoside phosphotransferase family and ten amino acid residues responsible for the phosphorylation of kanamycin were found in the amino acid sequence of APH(3')-Id. No mobile genetic element (MGE) was discovered surrounding the aph(3')-Id gene.

Conclusion

In this work, a novel aminoglycoside 3'-O-phosphotransferase gene designated aph(3')-Id encoded in the chromosome of the environmental isolate Kluyvera intermedia DW18 was identified and characterized. These findings will help clinicians select effective antimicrobials to treat infections caused by pathogens with this kind of resistance gene.

CasDinG is a 5'-3' dsDNA and RNA/DNA helicase with three accessory domains essential for type IV CRISPR immunity

CRISPR-associated DinG protein (CasDinG) is essential to type IV-A CRISPR function. Here, we demonstrate that CasDinG from Pseudomonas aeruginosa strain 83 is an ATP-dependent 5'-3' DNA translocase that unwinds double-stranded (ds)DNA and RNA/DNA hybrids. The crystal structure of CasDinG reveals a superfamily 2 helicase core of two RecA-like domains with three accessory domains (N-terminal, arch, and vestigial FeS). To examine the in vivo function of these domains, we identified the preferred PAM sequence for the type IV-A system (5'-GNAWN-3' on the 5'-side of the target) with a plasmid library and performed plasmid clearance assays with domain deletion mutants. Plasmid clearance assays demonstrated that all three domains are essential for type IV-A immunity. Protein expression and biochemical assays suggested the vFeS domain is needed for protein stability and the arch for helicase activity. However, deletion of the N-terminal domain did not impair ATPase, ssDNA binding, or helicase activities, indicating a role distinct from canonical helicase activities that structure prediction tools suggest involves interaction with dsDNA. This work demonstrates CasDinG helicase activity is essential for type IV-A CRISPR immunity as well as the yet undetermined activity of the CasDinG N-terminal domain.

Genome-Wide Identification and Expression Profiling Analysis of SWEET Family Genes Involved in Fruit Development in Plum (Prunus salicina Lindl)

SWEETs (sugars will eventually be exported transporters) play a vital role in longer-distance sugar transportation, and thus control carbon flow and energy metabolism in plants. SWEET genes have been identified in various plant species, but their functions in fruit development remain uncharacterized. Here, we isolated 15 putative PsSWEETs from the Prunus salicina genome. For further analysis, comprehensive bioinformatics methods were applied to determine the gene structure, chromosome distribution, phylogeny, cis-acting regulatory elements, and expression profiles of PsSWEETs. qRT-PCR analysis suggested that these SWEETs might have diverse functions in the development of plum fruit. The relative expression levels of PsSWEET1 and PsSWEET9 were obviously higher in ripened fruit than the ones in other developmental stages, suggesting their possible roles in the transport and accumulation of sugars in plum fruit. Positive correlations were found between the expression level of PsSWEET3/10/13 and the content of sucrose, and the expression level of PsSWEET2 and the content of fructose, respectively, during the development of 'Furongli' fruit, suggesting their possible roles in the accumulation of sucrose and fructose. The current study investigated the initial genomic characterization and expression patterns of the SWEET gene family in plum, which could provide a foundation for the further understanding of the functional analysis of the SWEET gene family.

Two conformations of the Tom20 preprotein receptor in the TOM holo complex

The TOM complex is the main entry point for precursor proteins (preproteins) into mitochondria. Preproteins containing targeting sequences are recognized by the TOM complex and imported into mitochondria. We have determined the structure of the TOM core complex from Neurospora crassa by single-particle electron cryomicroscopy at 3.3 Å resolution, showing its interaction with a bound preprotein at 4 Å resolution, and of the TOM holo complex including the Tom20 receptor at 6 to 7 Å resolution. TOM is a transmembrane complex consisting of two β-barrels, three receptor subunits, and three short transmembrane subunits. Tom20 has a transmembrane helix and a receptor domain on the cytoplasmic side. We propose that Tom20 acts as a dynamic gatekeeper, guiding preproteins into the pores of the TOM complex. We analyze the interactions of Tom20 with other TOM subunits, present insights into the structure of the TOM holo complex, and suggest a translocation mechanism.

DichroIDP: a method for analyses of intrinsically disordered proteins using circular dichroism spectroscopy

Intrinsically disordered proteins (IDPs) are comprised of significant numbers of residues that form neither helix, sheet, nor any other canonical type of secondary structure. They play important roles in a broad range of biological processes, such as molecular recognition and signalling, largely due to their chameleon-like ability to change structure from unordered when free in solution to ordered when bound to partner molecules. Circular dichroism (CD) spectroscopy is a widely-used method for characterising protein secondary structures, but analyses of IDPs using CD spectroscopy have suffered because the methods and reference datasets used for the empirical determination of secondary structures do not contain adequate representations of unordered structures. This work describes the creation, validation and testing of a standalone Windows-based application, DichroIDP, and a new reference dataset, IDP175, which is suitable for analyses of proteins containing significant amounts of disordered structure. DichroIDP enables secondary structure determinations of IDPs and proteins containing intrinsically disordered regions.

Analysis of the Tomato mTERF Gene Family and Study of the Stress Resistance Function of SLmTERF-13

Mitochondrial transcription termination factor (mTERF) is a DNA-binding protein that is encoded by nuclear genes, ultimately functions in mitochondria and can affect gene expression. By combining with mitochondrial nucleic acids, mTERF regulates the replication, transcription and translation of mitochondrial genes and plays an important role in the response of plants to abiotic stress. However, there are few studies on mTERF genes in tomato, which limits the in-depth study and utilization of mTERF family genes in tomato stress resistance regulation. In this study, a total of 28 mTERF gene family members were obtained through genome-wide mining and identification of the tomato mTERF gene family. Bioinformatics analysis showed that all members of the family contained environmental stress or hormone response elements. Gene expression pattern analysis showed that the selected genes had different responses to drought, high salt and low temperature stress. Most of the genes played key roles under drought and salt stress, and the response patterns were more similar. The VIGS method was used to silence the SLmTERF13 gene, which was significantly upregulated under drought and salt stress, and it was found that the resistance ability of silenced plants was decreased under both kinds of stress, indicating that the SLmTERF13 gene was involved in the regulation of the tomato abiotic stress response. These results provide important insights for further evolutionary studies and contribute to a better understanding of the role of the mTERF genes in tomato growth and development and abiotic stress response, which will ultimately play a role in future studies of tomato gene function.

Genome-Wide Analysis of Haemonchus contortus Proteases and Protease Inhibitors Using Advanced Informatics Provides Insights into Parasite Biology and Host-Parasite Interactions

Biodiversity within the animal kingdom is associated with extensive molecular diversity. The expansion of genomic, transcriptomic and proteomic data sets for invertebrate groups and species with unique biological traits necessitates reliable in silico tools for the accurate identification and annotation of molecules and molecular groups. However, conventional tools are inadequate for lesser-known organismal groups, such as eukaryotic pathogens (parasites), so that improved approaches are urgently needed. Here, we established a combined sequence- and structure-based workflow system to harness well-curated publicly available data sets and resources to identify, classify and annotate proteases and protease inhibitors of a highly pathogenic parasitic roundworm (nematode) of global relevance, called Haemonchus contortus (barber's pole worm). This workflow performed markedly better than conventional, sequence-based classification and annotation alone and allowed the first genome-wide characterisation of protease and protease inhibitor genes and gene products in this worm. In total, we identified 790 genes encoding 860 proteases and protease inhibitors representing 83 gene families. The proteins inferred included 280 metallo-, 145 cysteine, 142 serine, 121 aspartic and 81 "mixed" proteases as well as 91 protease inhibitors, all of which had marked physicochemical diversity and inferred involvements in >400 biological processes or pathways. A detailed investigation revealed a remarkable expansion of some protease or inhibitor gene families, which are likely linked to parasitism (e.g., host-parasite interactions, immunomodulation and blood-feeding) and exhibit stage- or sex-specific transcription profiles. This investigation provides a solid foundation for detailed explorations of the structures and functions of proteases and protease inhibitors of H. contortus and related nematodes, and it could assist in the discovery of new drug or vaccine targets against infections or diseases.

Spinocerebellar ataxia 38: structure-function analysis shows ELOVL5 G230V is proteotoxic, conformationally altered and a mutational hotspot

Fatty acid elongase ELOVL5 is part of a protein family of multipass transmembrane proteins that reside in the endoplasmic reticulum where they regulate long-chain fatty acid elongation. A missense variant (c.689G>T p.Gly230Val) in ELOVL5 causes Spinocerebellar Ataxia subtype 38 (SCA38), a neurodegenerative disorder characterized by autosomal dominant inheritance, cerebellar Purkinje cell demise and adult-onset ataxia. Having previously showed aberrant accumulation of p.G230V in the Golgi complex, here we further investigated the pathogenic mechanisms triggered by p.G230V, integrating functional studies with bioinformatic analyses of protein sequence and structure. Biochemical analysis showed that p.G230V enzymatic activity was normal. In contrast, SCA38-derived fibroblasts showed reduced expression of ELOVL5, Golgi complex enlargement and increased proteasomal degradation with respect to controls. By heterologous overexpression, p.G230V was significantly more active than wild-type ELOVL5 in triggering the unfolded protein response and in decreasing viability in mouse cortical neurons. By homology modelling, we generated native and p.G230V protein structures whose superposition revealed a shift in Loop 6 in p.G230V that altered a highly conserved intramolecular disulphide bond. The conformation of this bond, connecting Loop 2 and Loop 6, appears to be elongase-specific. Alteration of this intramolecular interaction was also observed when comparing wild-type ELOVL4 and the p.W246G variant which causes SCA34. We demonstrate by sequence and structure analyses that ELOVL5 p.G230V and ELOVL4 p.W246G are position-equivalent missense variants. We conclude that SCA38 is a conformational disease and propose combined loss of function by mislocalization and gain of toxic function by ER/Golgi stress as early events in SCA38 pathogenesis.

Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces

The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Exploring the role of detoxification genes in the resistance of Bursaphelenchus xylophilus to different exogenous nematicidal substances using transcriptomic analyses

Bursaphelenchus xylophilus (Pine wood nematode, PWN) has become a worldwide forest disease due to its rapid infection ability, high lethality and difficulty in control. The main means of countering B. xylophilus is currently chemical control, but nematicides can present problems such as environmental pollution and drug resistance. The development of novel environmentally-friendly nematicides has thus become a focus of recent research. In this study, BxUGT3 and BxUGT34, which might be related to detoxification, were investigated by comparing transcriptomic and WGCNA approaches. Three other genes with a similar expression pattern, BxUGT13, BxUGT14, and BxUGT16, were found by gene family analysis. Further bioassays and qPCR assays confirmed that these five genes showed significant changes in transcript levels upon exposure to α-pinene and carvone, demonstrating that they respond to exogenous nematicidal substances. Finally, RNAi and bioassays showed that B. xylophilus with silenced BxUGT16 had increased mortality in the face of α-pinene and carvone stress, suggesting that BxUGT16 plays an important role in detoxification. Taken together, this study used novel molecular research methods, explored the detoxification mechanism of B. xylophilus at a transcriptomic level, and revealed a molecular target for the development of novel biopesticides.

Mapping the Human Chondroitin Sulfate Glycoproteome Reveals an Unexpected Correlation Between Glycan Sulfation and Attachment Site Characteristics

Chondroitin sulfate proteoglycans (CSPGs) control key events in human health and disease and are composed of chondroitin sulfate (CS) polysaccharide(s) attached to different core proteins. Detailed information on the biological effects of site-specific CS structures is scarce as the polysaccharides are typically released from their core proteins prior to analysis. Here we present a novel glycoproteomic approach for site-specific sequencing of CS modifications from human urine. Software-assisted and manual analysis revealed that certain core proteins carried CS with abundant sulfate modifications, while others carried CS with lower levels of sulfation. Inspection of the amino acid sequences surrounding the attachment sites indicated that the acidity of the attachment site motifs increased the levels of CS sulfation, and statistical analysis confirmed this relationship. However, not only the acidity but also the sequence and characteristics of specific amino acids in the proximity of the serine glycosylation site correlated with the degree of sulfation. These results demonstrate attachment site-specific characteristics of CS polysaccharides of CSPGs in human urine and indicate that this novel method may assist in elucidating the biosynthesis and functional roles of CSPGs in cellular physiology.

Accelerating the Screening of Small Peptide Ligands by Combining Peptide-Protein Docking and Machine Learning

This research introduces a novel pipeline that couples machine learning (ML), and molecular docking for accelerating the process of small peptide ligand screening through the prediction of peptide-protein docking. Eight ML algorithms were analyzed for their potential. Notably, Light Gradient Boosting Machine (LightGBM), despite having comparable F1-score and accuracy to its counterparts, showcased superior computational efficiency. LightGBM was used to classify peptide-protein docking performance of the entire tetrapeptide library of 160,000 peptide ligands against four viral envelope proteins. The library was classified into two groups, 'better performers' and 'worse performers'. By training the LightGBM algorithm on just 1% of the tetrapeptide library, we successfully classified the remaining 99%with an accuracy range of 0.81-0.85 and an F1-score between 0.58-0.67. Three different molecular docking software were used to prove that the process is not software dependent. With an adjustable probability threshold (from 0.5 to 0.95), the process could be accelerated by a factor of at least 10-fold and still get 90-95% concurrence with the method without ML. This study validates the efficiency of machine learning coupled to molecular docking in rapidly identifying top peptides without relying on high-performance computing power, making it an effective tool for screening potential bioactive compounds.

Integrative determination of the atomic structure of mutant huntingtin exon 1 fibrils from Huntington's disease

Neurodegeneration in Huntington's disease (HD) is accompanied by the aggregation of fragments of the mutant huntingtin protein, a biomarker of disease progression. A particular pathogenic role has been attributed to the aggregation-prone huntingtin exon 1 (HttEx1) fragment, whose polyglutamine (polyQ) segment is expanded. Unlike amyloid fibrils from Parkinson's and Alzheimer's diseases, the atomic-level structure of HttEx1 fibrils has remained unknown, limiting diagnostic and treatment efforts. We present and analyze the structure of fibrils formed by polyQ peptides and polyQ-expanded HttEx1. Atomic-resolution perspectives are enabled by an integrative analysis and unrestrained all-atom molecular dynamics (MD) simulations incorporating experimental data from electron microscopy (EM), solid-state NMR, and other techniques. Visualizing the HttEx1 subdomains in atomic detail helps explaining the biological properties of these protein aggregates, as well as paves the way for targeting them for detection and degradation.

The glycoside hydrolase gene family profile and microbial function of Debaryomyces hansenii Y4 during South-road dark tea fermentation

Microbes are crucial to the quality formation of Sichuan South-road Dark Tea (SSDT) during pile-fermentation, but their mechanism of action has not yet been elucidated. Here, the glycoside hydrolase (GH) gene family and microbial function of Debaryomyces hansenii Y4 during solid-state fermentation were analyzed, and the results showed that many GH genes being distributed in comparatively abundant GH17, GH18, GH76, GH31, GH47, and GH2 were discovered in D. hansenii. They encoded beta-galactosidase, alpha-D-galactoside galactohydrolase, alpha-xylosidase, mannosidase, etc., and most of the GHs were located in the exocellular space and participated in the degradation of polysaccharides and oligosaccharides. D. hansenii Y4 could develop the mellow mouthfeel and "reddish brown" factors of SSDT via increasing the levels of water extracts, soluble sugars and amino acids but decreasing the tea polyphenols and caffeine levels, combined with altering the levels of thearubiins and brown index. It may facilitate the isomerization between epicatechin gallate and catechin gallate. Moreover, the expression levels of DEHA2G24860g (Beta-galactosidase gene) and DEHA2G08602g (Mannan endo-1,6-alpha-mannosidase DFG5 gene) were sharply up-regulated in fermentative anaphase, and they were significantly and negatively correlated with epicatechin content, especially, the expression of DEHA2G08602g was significantly and negatively correlated with catechin gallate level. It was hypothesized that D. hansenii Y4 is likely to be an important functional microbe targeting carbohydrate destruction and catechin transformation during SSDT pile-fermentation, with DEHA2G08602g as a key thermotolerant functional gene.

Modulation of the conformation, fibrillation, and fibril morphologies of human brain α-, β-, and γ-syn proteins by the disaccharide chemical chaperone trehalose

Human α-, β-, and γ-synuclein (syn) are natively unfolded proteins present in the brain. Deposition of aggregated α-syn in Lewy bodies is associated with Parkinson's disease (PD) and γ-syn is known to be involved in both neurodegeneration and breast cancer. At physiological pH, while α-syn has the highest propensity for fibrillation followed by γ-syn, β-syn does not form any fibrils. Fibril formation in these proteins could be modulated by protein structure stabilizing osmolytes such as trehalose which has an exceptional stabilizing effect for globular proteins. We present a comprehensive study of the effect of trehalose on the conformation, aggregation, and fibril morphology of α-, β-, and γ-syn proteins. Rather than stabilizing the intrinsically disordered state of the synucleins, trehalose accelerates the rate of fibril formation by forming aggregation-competent partially folded intermediate structures. Fibril morphologies are also strongly dependent on the concentration of trehalose with ≤ 0.4M favoring the formation of mature fibrils in α-, and γ-syn with no effect on the fibrillation of β-syn. At ≥ 0.8M, trehalose promotes the formation of smaller aggregates that are more cytotoxic. Live cell imaging of preformed aggregates of a labeled A90C α-syn shows their rapid internalization into neural cells which could be useful in reducing the load of aggregated species of α-syn. The findings throw light on the differential effect of trehalose on the conformation and aggregation of disordered synuclein proteins with respect to globular proteins and could help in understanding the effect of osmolytes on intrinsically disordered proteins under cellular stress conditions.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Gene Family in Litchi chinensis

NAC proteins play an essential role in the growth and development of litchi, especially during reproductive development. However, a comprehensive analysis of the litchi NAC gene family is currently absent. Based on information from the litchi genome, we found that the 112 NAC genes of litchi show an uneven distribution on the chromosomes. Phylogenetic and conserved structural domain analyses indicated that different types of variability were exhibited in the family of litchi NACs (LcNACs). Gene covariance analysis showed that the LcNACs showed better similarity in the same genus than with Arabidopsis. We further investigated the differential expression patterns of LcNACs in buds and rudimentary leaves of litchi. qRT-PCR results implied that they were involved in the process. Profiling of LcNAC promoter elements in litchi showed that they were extensively involved in light response, phytohormone regulation, abiotic stress response, and plant growth and development processes. This study provides new insights into the identification, structural characterization, tissue-specific expression analysis, and promoter response elements of LcNACs. It reveals the characteristics of the LcNACs and lays the foundation for the subsequent understanding of its biological functions and molecular regulatory mechanisms.

Genome-Wide Identification and Analysis of the Plant Cysteine Oxidase (PCO) Gene Family in Brassica napus and Its Role in Abiotic Stress Response

Plant Cysteine Oxidase (PCO) is a plant O₂-sensing enzyme catalyzing the oxidation of cysteine to Cys-sulfinic acid at the N-termini of target proteins. To better understand the Brassica napus PCO gene family, PCO genes in B. napus and related species were analyzed. In this study, 20, 7 and 8 PCO genes were identified in Brassica napus, Brassica rapa and Brassica oleracea, respectively. According to phylogenetic analysis, the PCOs were divided into five groups: PCO1, PCO2, PCO3, PCO4 and PCO5. Gene organization and motif distribution analysis suggested that the PCO gene family was relatively conserved during evolution. According to the public expression data, PCO genes were expressed in different tissues at different developmental stages. Moreover, qRT-PCR data showed that most of the Bna/Bra/BoPCO5 members were expressed in leaves, roots, flowers and siliques, suggesting an important role in both vegetative and reproductive development. Expression of BnaPCO was induced by various abiotic stress, especially waterlogging stress, which was consistent with the result of cis-element analysis. In this study, the PCO gene family of Brassicaceae was analyzed for the first time, which contributes to a comprehensive understanding of the origin and evolution of PCO genes in Brassicaceae and the function of BnaPCO in abiotic stress responses.

Genome-wide characterization of L-aspartate oxidase genes in wheat and their potential roles in the responses to wheat disease and abiotic stresses

L-aspartate oxidase (AO) is the first enzyme in NAD⁺ biosynthesis and is widely distributed in plants, animals, and microorganisms. Recently, AO family members have been reported in several plants, including Arabidopsis thaliana and Zea mays. Research on AO in these plants has revealed that AO plays important roles in plant growth, development, and biotic stresses; however, the nature and functions of AO proteins in wheat are still unclear. In this study, nine AO genes were identified in the wheat genome via sequence alignment and conserved protein domain analysis. These nine wheat AO genes (TaAOs) were distributed on chromosomes 2, 5, and 6 of sub-genomes A, B, and D. Analysis of the phylogenetic relationships, conserved motifs, and gene structure showed that the nine TaAOs were clustered into three groups, and the TaAOs in each group had similar conserved motifs and gene structure. Meanwhile, the subcellular localization analysis of transient expression mediated by Agrobacterium tumetioniens indicated that TaAO3-6D was localized to chloroplasts. Prediction of cis-elements indicated that a large number of cis-elements involved in responses to ABA, SA, and antioxidants/electrophiles, as well as photoregulatory responses, were found in TaAO promoters, which suggests that the expression of TaAOs may be regulated by these factors. Finally, transcriptome and real-time PCR analysis showed that the expression of TaAOs belonging to Group III was strongly induced in wheat infected by F. graminearum during anthesis, while the expression of TaAOs belonging to Group I was heavily suppressed. Additionally, the inducible expression of TaAOs belonging to Group III during anthesis in wheat spikelets infected by F. graminearum was repressed by ABA. Finally, expression of almost all TaAOs was induced by exposure to cold treatment. These results indicate that TaAOs may participate in the response of wheat to F. graminearum infection and cold stress, and ABA may play a negative role in this process. This study lays a foundation for further investigation of TaAO genes and provides novel insights into their biological functions.

Molecular and Functional Characterisation of a Novel Intragenic 12q24.21 Deletion Resulting in MED13L Haploinsufficiency Syndrome

Background and Objectives: Heterozygous pathogenic variants in the MED13L gene cause impaired intellectual development and distinctive facial features with or without cardiac defects (MIM #616789). This complex neurodevelopmental disorder is characterised by various phenotypic features, including plagiocephaly, strabismus, clubfoot, poor speech, and developmental delay. The aim of this study was to evaluate the clinical significance and consequences of a novel heterozygous intragenic MED13L deletion in a proband with clinical features of a MED13L-related disorder through extensive clinical, molecular, and functional characterisation. Materials and Methods: Combined comparative genomic hybridisation and single-nucleotide polymorphism array (SNP-CGH) was used to identify the changes in the proband's gDNA sequence (DECIPHER #430183). Intragenic MED13L deletion was specified via quantitative polymerase chain reaction (qPCR) and Sanger sequencing of the proband's cDNA sample. Western blot and bioinformatics analyses were used to investigate the consequences of this copy number variant (CNV) at the protein level. CRISPR-Cas9 technology was used for a MED13L-gene-silencing experiment in a culture of the control individual's skin fibroblasts. After the MED13L-gene-editing experiment, subsequent functional fibroblast culture analyses were performed. Results: The analysis of the proband's cDNA sample allowed for specifying the regions of the breakpoints and identifying the heterozygous deletion that spanned exons 3 to 10 of MED13L, which has not been reported previously. In silico, the deletion was predicted to result in a truncated protein NP_056150.1:p.(Val104Glyfs*5), partly altering the Med13_N domain and losing the MedPIWI and Med13_C domains. After MED13L gene editing was performed, reduced cell viability; an accelerated aging process; and inhibition of the RB1, E2F1, and CCNC gene expression were found to exist. Conclusions: Based on these findings, heterozygous intragenic 12q24.21 deletion in the affected individual resulted in MED13L haploinsufficiency due to the premature termination of protein translation, therefore leading to MED13L haploinsufficiency syndrome.

Biomolecular condensate drives polymerization and bundling of the bacterial tubulin FtsZ to regulate cell division

Cell division is spatiotemporally precisely regulated, but the underlying mechanisms are incompletely understood. In the social bacterium Myxococcus xanthus, the PomX/PomY/PomZ proteins form a single megadalton-sized complex that directly positions and stimulates cytokinetic ring formation by the tubulin homolog FtsZ. Here, we study the structure and mechanism of this complex in vitro and in vivo. We demonstrate that PomY forms liquid-like biomolecular condensates by phase separation, while PomX self-assembles into filaments generating a single large cellular structure. The PomX structure enriches PomY, thereby guaranteeing the formation of precisely one PomY condensate per cell through surface-assisted condensation. In vitro, PomY condensates selectively enrich FtsZ and nucleate GTP-dependent FtsZ polymerization and bundle FtsZ filaments, suggesting a cell division site positioning mechanism in which the single PomY condensate enriches FtsZ to guide FtsZ-ring formation and division. This mechanism shares features with microtubule nucleation by biomolecular condensates in eukaryotes, supporting this mechanism's ancient origin.

Cold Cas: reevaluating the occurrence of CRISPR/Cas systems in Mycobacteriaceae

Bacterial CRISPR/Cas systems target foreign genetic elements such as phages and regulate gene expression by some pathogens, even in the host. The system is a marker for evolutionary history and has been used for inferences in Mycobacterium tuberculosis for 30 years. However, knowledge about mycobacterial CRISPR/Cas systems remains limited. It is believed that Type III-A Cas systems are exclusive to Mycobacterium canettii and the M. tuberculosis complex (MTBC) of organisms and that very few of the >200 diverse species of non-tuberculous mycobacteria (NTM) possess any CRISPR/Cas system. This study sought unreported CRISPR/Cas loci across NTM to better understand mycobacterial evolution, particularly in species phylogenetically near the MTBC. An analysis of available mycobacterial genomes revealed that Cas systems are widespread across Mycobacteriaceae and that some species contain multiple types. The phylogeny of Cas loci shows scattered presence in many NTM, with variation even within species, suggesting gains/losses of these loci occur frequently. Cas Type III-A systems were identified in pathogenic Mycobacterium heckeshornense and the geological environmental isolate Mycobacterium SM1. In summary, mycobacterial CRISPR/Cas systems are numerous, Type III-A systems are unreliable as markers for MTBC evolution, and mycobacterial horizontal gene transfer appears to be a frequent source of genetic variation.

Emergence of ADC-5 Cephalosporinase in environmental Acinetobacter baumannii from a German tank milk with a novel Sequence Type

Bacteria resistant to antibiotics arguably pose the greatest threat to human health in the twenty-first century. One such bacterium that typifies antibiotic resistance is Acinetobacter baumannii . Frequently, hospital strains of A. baumannii display multidrug resistant (MDR) or extensively drug resistant (XDR) phenotypes, often requiring the use of last resort antibiotics for treatment. In addition to hospital settings, A. baumannii has been isolated from many highly divergent sources including wastewater treatment plant effluent, soil, and agricultural run-off with global distribution. However, such isolates remain poorly characterized. In this study, we characterized a strain of A. baumannii, AB341-IK15, isolated from bulk tank milk in Germany that demonstrated resistance to ceftazidime and intermediate resistance to ceftriaxone and piperacillin/tazobactam. Further genetic characterization identified an ADC-5 cephalosporinase, first incidence in an environmental isolate; and an OXA-408 oxacillinase that may contribute to this phenotype. Interestingly, AB341-IK15 is of a novel sequence type. This research underscores the importance of studying isolates of A. baumannii of non-clinical origin to understand the antibiotic resistance and virulence potential of environmental isolates of A. baumannii as well to understand the diversity of this species.

Design of a Multi-Epitope Vaccine against Tuberculosis from Mycobacterium tuberculosis PE_PGRS49 and PE_PGRS56 Proteins by Reverse Vaccinology

Tuberculosis is a disease caused by Mycobacterium tuberculosis, representing the second leading cause of death by an infectious agent worldwide. The available vaccine against this disease has insufficient coverage and variable efficacy, accounting for a high number of cases worldwide. In fact, an estimated third of the world's population has a latent infection. Therefore, developing new vaccines is crucial to preventing it. In this study, the highly antigenic PE_PGRS49 and PE_PGRS56 proteins were analyzed. These proteins were used for predicting T- and B-cell epitopes and for human leukocyte antigen (HLA) protein binding efficiency. Epitopes GGAGGNGSLSS, FAGAGGQGGLGG, GIGGGTQSATGLG (PE_PGRS49), and GTGWNGGKGDTG (PE_PGRS56) were selected based on their best physicochemical, antigenic, non-allergenic, and non-toxic properties and coupled to HLA I and HLA II structures for in silico assays. A construct with an adjuvant (RS09) plus each epitope joined by GPGPG linkers was designed, and the stability of the HLA-coupled construct was further evaluated by molecular dynamics simulations. Although experimental and in vivo studies are still necessary to ensure its protective effect against the disease, this study shows that the vaccine construct is dynamically stable and potentially effective against tuberculosis.

Comparative mitochondrial genomics of American nectar-feeding long-nosed bats Leptonycteris spp. with insights into the phylogeny of the family Phyllostomidae

Among leaf-nosed bats (family Phyllostomidae), the genus Leptonycteris (subfamily Glossophaginae), contains three migratory and obligate nectar-dwelling species of great ecological and economic importance; the greater long-nosed bat L. nivalis, the lesser long-nosed bat L. yerbabuenae, and the southern long-nosed bat L. curasoae. According to the IUCN, the three species are categorized as 'vulnerable', 'endangered', and 'near threatened', respectively. In this study, we assembled and characterized in detail the mitochondrial genome of Leptonycteris spp. and examined the phylogenetic position of this genus in the family Phyllostomidae based on protein coding genes (PCGs). The mitogenomes of L. nivalis, L. curasoae, and L. yerbabuenae are 16,708, 16,758, and 16,729 bp in length and each encode 13 PCGs, 22 transfer RNA genes, 2 rRNA genes, and a putative control region (CR). Mitochondrial gene order is identical to that reported before for the family Phyllostomidae. All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that is missing the DHU arm in the three species. All PCGs are exposed to purifying selection with atp8 experiencing the most relaxed purifying selection as the ω ratio was higher for this gene compared to the other PCGs in each species. The CR of each species contains three functional domains: extended termination associated sequence (ETAS), Central, and conserved sequence block domain (CSB). A phylomitogenomic analysis revealed that Leptonycteris is monophyletic and most closely related to the genus Glossophaga. The analysis also supported the monophyly of the family Glossophaginae in the speciose family Phyllostomidae. The mitochondria characterization of these species provides relevant information to develop molecular markers for conservation purposes.

Laterosporulin25: A probiotically produced, novel defensin-like bacteriocin and its immunogenic properties

Although multiple vaccines have been developed against infectious diseases, the rapid emergence of new pathogens develops an urgent need for novel strategies to combat infectious diseases. Antimicrobial peptides (AMPs) are excellent agents to fight against infectious diseases having unique multiple mechanisms of action against various pathogens. Apart from the direct applications, AMPs can also be developed as subunit vaccines or could be used as a highly immunogenic carrier protein with highly antigenic but non-immunogenic antigens. Here in the present study, we have identified a novel defensin-like bacteriocin, laterosporulin25 (LS25) upon genome mining of Brevibacillus laterosporus DSM25, a probiotic bacterial strain. By using immunoinformatic tools, we have studied the immunogenic and physiochemical properties of LS25. LS25 is characterized as defensin-like bacteriocin, having 51 amino acids and a molecular weight of 5862.7 Da. The modeled tertiary structure of LS25 is docked with TLR3 and TLR4-MD2 complex to confirm the facilitation of induced immune response that is further validated using molecular dynamics simulations and In-silico immune stimulations. Overall, detailed immunoinformatics analysis suggested LS25 as a potential candidate to be used as an adjuvant or carrier protein for subunit vaccine development, however, further in-vitro and in-vivo experiments are essential to validate its potential.

Genome-Wide Identification and Characterization of the TLP Gene Family in Phyllostachys edulis and Association with Witches' Broom Disease Resistance in Bamboo

Thaumatin-like proteins (TLPs) are pathogenesis-related proteins with pivotal roles in plant defense mechanisms. In this study, various bioinformatics and RNA-seq methods were used to analyze the biotic and abiotic stress responses of the TLP family in Phyllostachys edulis. Overall, 81 TLP genes were identified in P. edulis; 166 TLPs from four plant species were divided into three groups and ten subclasses, with genetic covariance observed between these species. Subcellular localization in silico studies indicated that TLPs were primarily distributed in the extracellular. Analysis of the upstream sequences of TLPs demonstrated the presence of cis-acting elements related to disease defense, environmental stress, and hormonal responses. Multiple sequence alignment demonstrated that most TLPs possessed five conserved REDDD amino acid sequences with only a few amino acid residue differences. RNA-seq analysis of P. edulis responses to Aciculosporium take, the pathogenic fungus that causes witches' broom disease, showed that P. edulis TLPs (PeTLPs) were expressed in different organs, with the highest expression in buds. PeTLPs responded to both abscisic acid and salicylic acid stress. These PeTLP expression patterns were consistent with their gene and protein structures. Collectively, our findings provide a basis for further comprehensive analyses of the genes related to witches' broom in P. edulis.

The LiaFSR-LiaX System Mediates Resistance of Enterococcus faecium to Peptide Antibiotics and to Aureocin A53- and Enterocin L50-Like Bacteriocins

Multidrug-resistant Enterococcus faecium strains are currently a leading cause of difficult-to-treat nosocomial infections. The emerging resistance of enterococci to last-resort antibiotics, such as daptomycin, prompts a search for alternative antimicrobials. Aureocin A53- and enterocin L50-like bacteriocins are potent antimicrobial agents that form daptomycin-like cationic complexes and have a similar cell envelope-targeting mechanism of action, suggesting their potential as next-generation antibiotics. However, to ensure their safe use, the mechanisms of resistance to these bacteriocins and cross-resistance to antibiotics need to be well understood. Here, we investigated the genetic basis of E. faecium's resistance to aureocin A53- and enterocin L50-like bacteriocins and compared it with that to antibiotics. First, we selected spontaneous mutants resistant to the bacteriocin BHT-B and identified adaptive mutations in the liaFSR-liaX genes encoding the LiaFSR stress response regulatory system and the daptomycin-sensing protein LiaX, respectively. We then demonstrated that a gain-of-function mutation in liaR increases the expression of liaFSR, liaXYZ, cell wall remodeling-associated genes, and hypothetical genes involved in protection against various antimicrobials. Finally, we showed that adaptive mutations or overexpression of liaSR or liaR alone results in cross-resistance to other aureocin A53- and enterocin L50-like bacteriocins, as well as antibiotics targeting specific components of the cell envelope (daptomycin, ramoplanin, gramicidin) or ribosomes (kanamycin and gentamicin). Based on the obtained results, we concluded that activation of the LiaFSR-mediated stress response confers resistance to peptide antibiotics and bacteriocins via a cascade of reactions, eventually leading to cell envelope remodeling. IMPORTANCE Pathogenic enterococci carry virulence factors and a considerable resistome, which makes them one of the most serious and steadily increasing causes of hospital epidemiological risks. Accordingly, Enterococcus faecium is classified into a top-priority ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) group of six highly virulent and multidrug-resistant (MDR) bacterial pathogens for which novel antimicrobial agents need to be developed urgently. Alternative measures, such as the use of bacteriocins, separately or in combination with other antimicrobial agents (e.g., antibiotics), could be a potential solution, especially since several international health agencies recommend and support the development of such interventions. Nevertheless, in order to exploit their efficacy, more basic research on the mechanisms of cell killing and the development of resistance to bacteriocins is needed. The present study fills some of the knowledge gaps regarding the genetic basis of the development of resistance to potent antienterococcal bacteriocins, pointing out the common and divergent features regarding the cross-resistance to antibiotics.

Genome-wide analysis revealed the stepwise origin and functional diversification of HSDs from lower to higher plant species

Hydroxysteroid dehydrogenase (HSDs) is an oil-body sterol protein (steroleosin) with an NADP(H) binding domain that belongs to the short-chain dehydrogenase/reductase (SDR) superfamily. There are numerous studies on the characterization of HSDs in plants. However, thus far, the evolutionary differentiation and divergence analysis of these genes remain to be explored. The current study used an integrated method to elucidate the sequential evolution of HSDs in 64 sequenced plant genomes. Analyses were conducted on their origins, distribution, duplication, evolutionary paths, domain functions, motif composition, properties, and cis-elements. Results indicate that except for algae, HSD1 was widely distributed in plant species ranging from lower to higher plants, while HSD5 was restricted to terrestrial plants, and HSD2 was identified in fewer monocots and several dicot plants. Phylogenetic analysis of HSD proteins revealed that monocotyledonous HSD1 in moss and ferns appeared closest to the outgroup, V. carteri HSD-like, M. musculus HSD1, and H. sapiens HSD1. These data support the hypothesis that HSD1 originated in bryophytes and then in non-vascular and vascular plants, followed by HSD5 only in land plants. Gene structure analysis suggests that HSDs in plant species came up with a fixed number of six exons, and the intron phase was primarily 0, 1, 0, 0, and 0. Similarly, duplication analysis revealed that segmental duplications were the main reason for HSDs in plant species. Physicochemical properties suggest that dicotyledonous HSD1s and HSD5s were mainly acidic. The monocotyledonous HSD1s and HSD2s and the dicotyledonous HSD2s, HSD3s, HSD4s, and HSD6s were mainly basic, implying that HSDs in plants may have a variety of functions. Cis-regulatory elements and expression analysis revealed that HSDs in plants might have roles in several abiotic stresses. Due to the high expression of HSD1s and HSD5s in seeds, these HSDs in plants may have roles in fatty acid accumulation and degradation.

Spike-Independent Infection of Human Coronavirus 229E in Bat Cells

Bats are the reservoir for numerous human pathogens, including coronaviruses. Despite many coronaviruses having descended from bat ancestors, little is known about virus-host interactions and broader evolutionary history involving bats. Studies have largely focused on the zoonotic potential of coronaviruses with few infection experiments conducted in bat cells. To determine genetic changes derived from replication in bat cells and possibly identify potential novel evolutionary pathways for zoonotic virus emergence, we serially passaged six human 229E isolates in a newly established Rhinolophus lepidus (horseshoe bat) kidney cell line. Here, we observed extensive deletions within the spike and open reading frame 4 (ORF4) genes of five 229E viruses after passaging in bat cells. As a result, spike protein expression and infectivity of human cells was lost in 5 of 6 viruses, but the capability to infect bat cells was maintained. Only viruses that expressed the spike protein could be neutralized by 229E spike-specific antibodies in human cells, whereas there was no neutralizing effect on viruses that did not express the spike protein inoculated on bat cells. However, one isolate acquired an early stop codon, abrogating spike expression but maintaining infection in bat cells. After passaging this isolate in human cells, spike expression was restored due to acquisition of nucleotide insertions among virus subpopulations. Spike-independent infection of human coronavirus 229E may provide an alternative mechanism for viral maintenance in bats that does not rely on the compatibility of viral surface proteins and known cellular entry receptors. IMPORTANCE Many viruses, including coronaviruses, originated from bats. Yet, we know little about how these viruses switch between hosts and enter human populations. Coronaviruses have succeeded in establishing in humans at least five times, including endemic coronaviruses and the recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In an approach to identify requirements for host switches, we established a bat cell line and adapted human coronavirus 229E viruses by serial passage. The resulting viruses lost their spike protein but maintained the ability to infect bat cells, but not human cells. Maintenance of 229E viruses in bat cells appears to be independent of a canonical spike receptor match, which in turn might facilitate cross-species transmission in bats.

Genome-Wide Identification of B-Box Gene Family and Expression Analysis Suggest Its Roles in Responses to Cercospora Leaf Spot in Sugar Beet (Beta Vulgaris L.)

The B-box (BBX) protein, which is a zinc-finger protein containing one or two B-box domains, plays a crucial role in the growth and development of plants. Plant B-box genes are generally involved in morphogenesis, the growth of floral organs, and various life activities in response to stress. In this study, the sugar beet B-box genes (hereafter referred to as BvBBXs) were identified by searching the homologous sequences of the Arabidopsis thaliana B-box gene family. The gene structure, protein physicochemical properties, and phylogenetic analysis of these genes were systematically analyzed. In this study, 17 B-box gene family members were identified from the sugar beet genome. A B-box domain can be found in all sugar beet BBX proteins. BvBBXs encode 135 to 517 amino acids with a theoretical isoelectric point of 4.12 to 6.70. Chromosome localization studies revealed that BvBBXs were dispersed across nine sugar beet chromosomes except chromosomes 5 and 7. The sugar beet BBX gene family was divided into five subfamilies using phylogenetic analysis. The gene architectures of subfamily members on the same evolutionary tree branch are quite similar. Light, hormonal, and stress-related cis-acting elements can be found in the promoter region of BvBBXs. The BvBBX gene family was differently expressed in sugar beet following Cercospora leaf spot infection, according to RT-qPCR data. It is shown that the BvBBX gene family may influence how the plant reacts to a pathogen infection.

Design and Evolution of Enhanced Peptide-Peptide Ligation for Modular Transglutaminase Assembly

Robust and precise tools are needed to enhance the functionality and resilience of synthetic nanoarchitectures. Here, we have employed directed evolution and rational design to build a fast-acting molecular superglue from a bacterial adhesion protein. We have generated the SnoopLigase2 coupling system, a genetically encoded route for efficient transamidation between SnoopTag2 and DogTag2 peptides. Each peptide was selected for rapid reaction by phage display screening. The optimized set allows more than 99% completion and is compatible with diverse buffers, pH values, and temperatures, accelerating the reaction over 1000-fold. SnoopLigase2 directs a specific reaction in the mammalian secretory pathway, allowing covalent display on the plasma membrane. Transglutaminase 2 (TG2) has a network of interactions and substrates amidst the mammalian cell surface and extracellular matrix. We expressed a modified TG2 with resistance to oxidative inactivation and minimal self-reactivity. SnoopLigase2 enables TG2 functionalization with transforming growth factor alpha (TGFα) in routes that would be impossible through genetic fusion. The TG2:TGFα conjugate retained transamidase activity, stably anchored TGFα for signal activation in the extracellular environment, and reprogrammed cell behavior. This modular toolbox should create new opportunities for molecular assembly, both for novel biomaterials and complex cellular environments.

Classification of porcine reproductive and respiratory syndrome clinical impact in Ontario sow herds using machine learning approaches

Since the early 1990s, porcine reproductive and respiratory syndrome (PRRS) virus outbreaks have been reported across various parts of North America, Europe, and Asia. The incursion of PRRS virus (PRRSV) in swine herds could result in various clinical manifestations, resulting in a substantial impact on the incidence of respiratory morbidity, reproductive loss, and mortality. Veterinary experts, among others, regularly analyze the PRRSV open reading frame-5 (ORF-5) for prognostic purposes to assess the risk of severe clinical outcomes. In this study, we explored if predictive modeling techniques could be used to identify the severity of typical clinical signs observed during PRRS outbreaks in sow herds. Our study aimed to evaluate four baseline machine learning (ML) algorithms: logistic regression (LR) with ridge and lasso regularization techniques, random forest (RF), k-nearest neighbor (KNN), and support vector machine (SVM), for the clinical impact classification of ORF-5 sequences and demographic data into high impact and low impact categories. First, baseline classifiers were evaluated using different input representations of ORF-5 nucleotides, amino acid sequences, and demographic data using a 10-fold cross-validation technique. Then, we designed a consensus voting ensemble approach to aggregate the different types of input representations for genetic and demographic data for classifying clinical impact. In this study, we observed that: (a) for abortion and pre-weaning mortality (PWM), different classifiers gained improvement over baseline accuracy, which showed the plausible presence of both genotypic-phenotypic and demographic-phenotypic relationships, (b) for sow mortality (SM), no baseline classifier successfully established such linkages using either genetic or demographic input data, (c) baseline classifiers showed good performance with a moderate variance of the performance metrics, due to high-class overlap and the small dataset size used for training, and (d) the use of consensus voting ensemble techniques helped to make the predictions more robust and stabilized the performance evaluation metrics, but overall accuracy did not substantially improve the diagnostic metrics over baseline classifiers.

A Phenylboronic Acid-Based Transition State Analogue Yields Nanomolar Inhibition of Mandelate Racemase

Mandelate racemase (MR) catalyzes the Mg²⁺-dependent interconversion of (R)- and (S)-mandelate by stabilizing the altered substrate in the transition state (TS) by ∼26 kcal/mol. The enzyme has been employed as a model to explore the limits to which the free energy of TS stabilization may be captured by TS analogues to effect strong binding. Herein, we determined the thermodynamic parameters accompanying binding of a series of bromo-, chloro-, and fluoro-substituted phenylboronic acids (PBAs) by MR and found that binding was predominately driven by favorable entropy changes. 3,4-Dichloro-PBA was discovered to be the most potent inhibitor yet identified for MR, binding with a K_d^app value of 11 ± 2 nM and exceeding the binding of the substrate by ∼72,000-fold. The ΔC_p value accompanying binding (-488 ± 18 cal·mol^-1 K^-1) suggested that dispersion forces contribute significantly to the binding. The pH-dependence of the inhibition revealed that MR preferentially binds the anionic, tetrahedral form of 3,4-dichloro-PBA with a pH-independent K_i value of 5.7 ± 0.5 nM, which was consistent with the observed upfield shift of the ¹¹B NMR signal. The linear free energy relationship between log(k_cat/K_m) and log(1/K_i) for wild-type and 11 MR variants binding 3,4-dichloro-PBA had a slope of 0.8 ± 0.2, indicating that MR recognizes the inhibitor as an analogue of the TS. Hence, halogen substitution may be utilized to capture additional free energy of TS stabilization arising from dispersion forces to enhance the binding of boronic acid inhibitors by MR.

Functional Diversity of Homologous Oxidoreductases-Tuning of Substrate Specificity by a FAD-Stacking Residue for Iron Acquisition and Flavodoxin Reduction

Although bacterial thioredoxin reductase-like ferredoxin/flavodoxin NAD(P)⁺ oxidoreductases (FNRs) are similar in terms of primary sequences and structures, they participate in diverse biological processes by catalyzing a range of different redox reactions. Many of the reactions are critical for the growth, survival of, and infection by pathogens, and insight into the structural basis for substrate preference, specificity, and reaction kinetics is crucial for the detailed understanding of these redox pathways. Bacillus cereus (Bc) encodes three FNR paralogs, two of which have assigned distinct biological functions in bacillithiol disulfide reduction and flavodoxin (Fld) reduction. Bc FNR2, the endogenous reductase of the Fld-like protein NrdI, belongs to a distinct phylogenetic cluster of homologous oxidoreductases containing a conserved His residue stacking the FAD cofactor. In this study, we have assigned a function to FNR1, in which the His residue is replaced by a conserved Val, in the reduction of the heme-degrading monooxygenase IsdG, ultimately facilitating the release of iron in an important iron acquisition pathway. The Bc IsdG structure was solved, and IsdG-FNR1 interactions were proposed through protein-protein docking. Mutational studies and bioinformatics analyses confirmed the importance of the conserved FAD-stacking residues on the respective reaction rates, proposing a division of FNRs into four functionally unique sequence similarity clusters likely related to the nature of this residue.

The Role of Asn11 in Catalysis by Triosephosphate Isomerase

Four catalytic amino acids at triosephosphate isomerase (TIM) are highly conserved: N11, K13, H95, and E167. Asparagine 11 is the last of these to be characterized in mutagenesis studies. The ND2 side chain atom of N11 is hydrogen bonded to the O-1 hydroxyl of enzyme-bound dihydroxyacetone phosphate (DHAP), and it sits in an extended chain of hydrogen-bonded side chains that includes T75' from the second subunit. The N11A variants of wild-type TIM from Trypanosoma brucei brucei (TbbTIM) and Leishmania mexicana (LmTIM) undergo dissociation from the dimer to monomer under our assay conditions. Values of Kas = 8 × 103 and 1 × 106 M-1, respectively, were determined for the conversion of monomeric N11A TbbTIM and LmTIM into their homodimers. The N11A substitution at the variant of LmTIM previously stabilized by the E65Q substitution gives the N11A/E65Q variant that is stable to dissociation under our assay conditions. The X-ray crystal structure of N11A/E65Q LmTIM shows an active site that is essentially superimposable on that for wild-type TbbTIM, which also has a glutamine at position 65. A comparison of the kinetic parameters for E65Q LmTIM and N11A/E65Q LmTIM-catalyzed reactions of (R)-glyceraldehyde 3-phosphate (GAP) and (DHAP) shows that the N11A substitution results in a (13-14)-fold decrease in kcat/Km for substrate isomerization and a similar decrease in kcat for DHAP but only a 2-fold decrease in kcat for GAP.

Cathepsins in oral diseases: mechanisms and therapeutic implications

Cathepsins are a type of lysosomal globulin hydrolase and are crucial for many physiological processes, including the resorption of bone matrix, innate immunity, apoptosis, proliferation, metastasis, autophagy, and angiogenesis. Findings regarding their functions in human physiological processes and disorders have drawn extensive attention. In this review, we will focus on the relationship between cathepsins and oral diseases. We highlight the structural and functional properties of cathepsins related to oral diseases, as well as the regulatory mechanisms in tissue and cells and their therapeutic uses. Elucidating the associated mechanism between cathepsins and oral diseases is thought to be a promising strategy for the treatment of oral diseases and may be a starting point for further studies at the molecular level.

Genome-wide identification and analysis of the invertase gene family in tobacco (Nicotiana tabacum) reveals NtNINV10 participating the sugar metabolism

Sucrose (Suc) is directly associated with plant growth and development as well as tolerance to various stresses. Invertase (INV) enzymes played important role in sucrose metabolism by irreversibly catalyzing Suc degradation. However, genome-wide identification and function of individual members of the INV gene family in Nicotiana tabacum have not been conducted. In this report, 36 non-redundant NtINV family members were identified in Nicotiana tabacum including 20 alkaline/neutral INV genes (NtNINV1-20), 4 vacuolar INV genes (NtVINV1-4), and 12 cell wall INV isoforms (NtCWINV1-12). A comprehensive analysis based on the biochemical characteristics, the exon-intron structures, the chromosomal location and the evolutionary analysis revealed the conservation and the divergence of NtINVs. For the evolution of the NtINV gene, fragment duplication and purification selection were major factors. Besides, our analysis revealed that NtINV could be regulated by miRNAs and cis-regulatory elements of transcription factors associated with multiple stress responses. In addition, 3D structure analysis has provided evidence for the differentiation between the NINV and VINV. The expression patterns in diverse tissues and under various stresses were investigated, and qRT-PCR experiments were conducted to confirm the expression patterns. Results revealed that changes in NtNINV10 expression level were induced by leaf development, drought and salinity stresses. Further examination revealed that the NtNINV10-GFP fusion protein was located in the cell membrane. Furthermore, inhibition of the expression of NtNINV10 gene decreased the glucose and fructose in tobacco leaves. Overall, we have identified possible NtINV genes functioned in leaf development and tolerance to environmental stresses in tobacco. These findings provide a better understanding of the NtINV gene family and establish the basis for future research.

Identification of C2H2 zinc finger genes through genome-wide association study and functional analyses of LkZFPs in response to stresses in Larix kaempferi

BACKGROUND

C2H2 zinc finger proteins (C2H2-ZFPs), one of the largest transcription factors, play a variety of roles in plant development and growth as well as stress response. While, the evolutionary history and expression profile of the C2H2-ZFP genes in Larix kaempferi (LkZFPs) have not been reported so far.

RESULTS

In this study, the whole genome of the LkZFPs was identified and characterized, including physicochemical properties, phylogenetic relationships, conservative motifs, the promoter cis-elements and Gene Ontology (GO) annotation. We identified 47 LkZFPs and divided them into four subfamilies based on phylogenetic analysis and conserved motifs. Subcellular localization prediction showed that most of the LkZFPs were located in the nucleus. Promoter cis-element analysis suggested that the LkZFPs may be involved in the regulation of stress responses. Moreover, Real-time quantitative PCR (RT-qPCR) results showed that Q-type LkZFP genes were involved in the response to abiotic stress, such as salt, drought and hormone stresses. Subcellular localization results showed that LkZFP7 and LkZFP37 were located in the nucleus, LkZFP32 was located in both cytoplasm and nucleus.

CONCLUSION

The identification and functional analysis of LkZFPs suggested that some LkZFP genes might play important roles in coping with both biological and abiotic stresses. These results could further increase understanding of the function of the LkZFPs, and provide some research direction and theoretical support.

Over-expression of CcMYB24, encoding a R2R3-MYB transcription factor from a high-leaf-number mutant of Cymbidium, increases the number of leaves in Arabidopsis

Ornamental foliage plants have long been cultivated for their attractive leaves. Variation in leaf traits of ornamental foliage plants is one of the goals in breeding. MYB transcription factors regulate many aspects of leaf development, and thus influence morphological traits of leaves. However, little is known about the function of MYB transcription factors in leaf development of Cymbidium, one of the most economically important ornamental plants in the world. In the present study, a MYB transcription factor, CcMYB24, was identified and the corresponding gene cloned from a new orchid mutant, TRIR-2, which produces more leaves than control plants. The CcMYB24 showed a higher expression level in 'TRIR-2' than in control plants, and the protein was located in the nucleus. The sequence of CcMYB24 showed a high similarity with RAX2-like genes which belong to the R2R3-MYB gene family in other Cymbidium plants. Overexpression of CcMYB24 resulted in a phenotype with an increased number of leaves, elevated chlorophyll content, and decreased contents of carotenoids and flavonoids in Arabidopsis. These results provide functional evidence for the role of CcMYB24 in promoting the production of leaves in 'TRIR-2'. Understanding the role of CcMYB24 in Cymbidium will be beneficial for the molecular breeding of ornamental foliage plants.