PredictProtein - Predicting Protein Structure and Function for 29 Years

MicroProteinDB: A database to provide knowledge on sequences, structures and function of ncRNA-derived microproteins

Omics-based technologies have revolutionized our comprehension of microproteins encoded by ncRNAs, revealing their abundant presence and pivotal roles within complex functional landscapes. Here, we developed MicroProteinDB (http://bio-bigdata.hrbmu.edu.cn/MicroProteinDB), which offers and visualizes the extensive knowledge to aid retrieval and analysis of computationally predicted and experimentally validated microproteins originating from various ncRNA types. Employing prediction algorithms grounded in diverse deep learning approaches, MicroProteinDB comprehensively documents the fundamental physicochemical properties, secondary and tertiary structures, interactions with functional proteins, family domains, and inter-species conservation of microproteins. With five major analytical modules, it will serve as a valuable knowledge for investigating ncRNA-derived microproteins.

The RHCE gene encodes the chicken blood system I

BACKGROUND

There are 13 known chicken blood systems, which were originally detected by agglutination of red blood cells by specific alloantisera. The genomic region or specific gene responsible has been identified for four of these systems (A, B, D and E). We determined the identity of the gene responsible for the chicken blood system I, using DNA from multiple birds with known chicken I blood system serology, 600K and 54K single nucleotide polymorphism (SNP) data, and lowpass sequence information.

RESULTS

The gene responsible for the chicken I blood system was identified as RHCE, which is also one of the genes responsible for the highly polymorphic human Rh blood group locus, for which maternal/fetal antigenic differences can result in fetal hemolytic anemia with fetal mortality. We identified 17 unique RHCE haplotypes in the chicken, with six haplotypes corresponding to known I system serological alleles. We also detected deletions in the RHCE gene that encompass more than 6000 bp and that are predicted to remove its last seven exons.

CONCLUSIONS

RHCE is the gene responsible for the chicken I blood system. This is the fifth chicken blood system for which the responsible gene and gene variants are known. With rapid DNA-based testing now available, the impact of I blood system variation on response against disease, general immune function, and animal production can be investigated in greater detail.

ASH1L guards cis-regulatory elements against cyclobutane pyrimidine dimer induction

The histone methyltransferase ASH1L, first discovered for its role in transcription, has been shown to accelerate the removal of ultraviolet (UV) light-induced cyclobutane pyrimidine dimers (CPDs) by nucleotide excision repair. Previous reports demonstrated that CPD excision is most efficient at transcriptional regulatory elements, including enhancers, relative to other genomic sites. Therefore, we analyzed DNA damage maps in ASH1L-proficient and ASH1L-deficient cells to understand how ASH1L controls enhancer stability. This comparison showed that ASH1L protects enhancer sequences against the induction of CPDs besides stimulating repair activity. ASH1L reduces CPD formation at C-containing but not at TT dinucleotides, and no protection occurs against pyrimidine-(6,4)-pyrimidone photoproducts or cisplatin crosslinks. The diminished CPD induction extends to gene promoters but excludes retrotransposons. This guardian role against CPDs in regulatory elements is associated with the presence of H3K4me3 and H3K27ac histone marks, which are known to interact with the PHD and BRD motifs of ASH1L, respectively. Molecular dynamics simulations identified a DNA-binding AT hook of ASH1L that alters the distance and dihedral angle between neighboring C nucleotides to disfavor dimerization. The loss of this protection results in a higher frequency of C->T transitions at enhancers of skin cancers carrying ASH1L mutations compared to ASH1L-intact counterparts.

Genomic exploration of the fermented meat isolate Staphylococcus shinii IMDO-S216 with a focus on competitiveness-enhancing secondary metabolites

BACKGROUND

Staphylococcus shinii appears as an umbrella species encompassing several strains of Staphylococcus pseudoxylosus and Staphylococcus xylosus. Given its phylogenetic closeness to S. xylosus, S. shinii can be found in similar ecological niches, including the microbiota of fermented meats where the species may contribute to colour and flavour development. In addition to these conventional functionalities, a biopreservation potential based on the production of antagonistic compounds may be available. Such potential, however, remains largely unexplored in contrast to the large body of research that is available on the biopreservative properties of lactic acid bacteria. The present study outlines the exploration of the genetic basis of competitiveness and antimicrobial activity of a fermented meat isolate, S. shinii IMDO-S216. To this end, its genome was sequenced, de novo assembled, and annotated.

RESULTS

The genome contained a single circular chromosome and eight plasmid replicons. Focus of the genomic exploration was on secondary metabolite biosynthetic gene clusters coding for ribosomally synthesized and posttranslationally modified peptides. One complete cluster was coding for a bacteriocin, namely lactococcin 972; the genes coding for the pre-bacteriocin, the ATP-binding cassette transporter, and the immunity protein were also identified. Five other complete clusters were identified, possibly functioning as competitiveness factors. These clusters were found to be involved in various responses such as membrane fluidity, iron intake from the medium, a quorum sensing system, and decreased sensitivity to antimicrobial peptides and competing microorganisms. The presence of these clusters was equally studied among a selection of multiple Staphylococcus species to assess their prevalence in closely-related organisms.

CONCLUSIONS

Such factors possibly translate in an improved adaptation and competitiveness of S. shinii IMDO-S216 which are, in turn, likely to improve its fitness in a fermented meat matrix.

Cereulide production capacities and genetic properties of 31 emetic Bacillus cereus group strains

The highly potent toxin cereulide is a frequent cause of foodborne intoxications. This extremely resistant toxin is produced by Bacillus cereus group strains carrying the plasmid encoded cesHPTABCD gene cluster. It is known that the capacities to produce cereulide vary greatly between different strains but the genetic background of these variations is not clear. In this study, cereulide production capacities were associated with genetic characteristics. For this, cereulide levels in cultures of 31 strains were determined after incubation in tryptic soy broth for 24 h at 24 °C, 30 °C and 37 °C. Whole genome sequencing based data were used for an in-depth characterization of gene sequences related to cereulide production. The taxonomy, population structure and phylogenetic relationships of the strains were evaluated based on average nucleotide identity, multi-locus sequence typing (MLST), core genome MLST and single nucleotide polymorphism analyses. Despite a limited strain number, the approach of a genome wide association study (GWAS) was tested to link genetic variation with cereulide quantities. Our study confirms strain-dependent differences in cereulide production. For most strains, these differences were not explainable by sequence variations in the cesHPTABCD gene cluster or the regulatory genes abrB, spo0A, codY and pagRBc. Likewise, the population structure and phylogeny of the tested strains did not comprehensively reflect the cereulide production capacities. GWAS yielded first hints for associated proteins, while their possible effect on cereulide synthesis remains to be further investigated.

Experimental and computational approaches for membrane protein insertion and topology determination

Membrane proteins play pivotal roles in a wide array of cellular processes and constitute approximately a quarter of the protein-coding genes across all organisms. Despite their ubiquity and biological significance, our understanding of these proteins remains notably less comprehensive compared to their soluble counterparts. This disparity in knowledge can be attributed, in part, to the inherent challenges associated with employing specialized techniques for the investigation of membrane protein insertion and topology. This review will center on a discussion of molecular biology methodologies and computational prediction tools designed to elucidate the insertion and topology of helical membrane proteins.

Altered socio-affective communication and amygdala development in mice with protocadherin10-deficient interneurons

Autism spectrum disorder (ASD) is a group of neurodevelopmental conditions associated with deficits in social interaction and communication, together with repetitive behaviours. The cell adhesion molecule protocadherin10 (PCDH10) is linked to ASD in humans. Pcdh10 is expressed in the nervous system during embryonic and early postnatal development and is important for neural circuit formation. In mice, strong expression of Pcdh10 in the ganglionic eminences and in the basolateral complex (BLC) of the amygdala was observed at mid and late embryonic stages, respectively. Both inhibitory and excitatory neurons expressed Pcdh10 in the BLC at perinatal stages and vocalization-related genes were enriched in Pcdh10-expressing neurons in adult mice. An epitope-tagged Pcdh10-HAV5 mouse line revealed endogenous interactions of PCDH10 with synaptic proteins in the young postnatal telencephalon. Nuanced socio-affective communication changes in call emission rates, acoustic features and call subtype clustering were primarily observed in heterozygous pups of a conditional knockout (cKO) with selective deletion of Pcdh10 in Gsh2-lineage interneurons. These changes were less prominent in heterozygous ubiquitous Pcdh10 KO pups, suggesting that altered anxiety levels associated with Gsh2-lineage interneuron functioning might drive the behavioural effects. Together, loss of Pcdh10 specifically in interneurons contributes to behavioural alterations in socio-affective communication with relevance to ASD.

Signaling Transduction Pathways and G-Protein-Coupled Receptors in Different Stages of the Embryonic Diapause Termination Process in Artemia

Artemia is a widely distributed small aquatic crustacean, renowned for its ability to enter a state of embryonic diapause. The embryonic diapause termination (EDT) is closely linked to environmental cues, but the precise underlying mechanisms remain elusive. In this study, ATAC-seq and RNA-seq sequencing techniques were employed to explore the gene expression profiles in Artemia cysts 30 min after EDT. These profiles were compared with those during diapause and 5 h after EDT. The regulatory mechanisms governing the EDT process were analyzed through Gene Ontology (GO) enrichment analysis of differentially expressed genes. Furthermore, the active G-protein-coupled receptors (GPCRs) were identified through structural analysis. The results unveiled that the signaling transduction during EDT primarily hinges on GPCRs and the cell surface receptor signaling pathway, but distinct genes are involved across different stages. Hormone-mediated signaling pathways and the tachykinin receptor signaling pathway exhibited heightened activity in the '0-30 min' group, whereas the Wnt signaling pathway manifested its function solely in the '30 min-5 h' group. These results imply a complete divergence in the mechanisms of signal regulation during these two stages. Moreover, through structural analysis, five GPCRs operating at different stages of EDT were identified. These findings provide valuable insights into the signal regulation mechanisms governing Artemia diapause.

Metabolic Investigation and Auxiliary Enzyme Modelization of the Pyrrocidine Pathway Allow Rationalization of Paracyclophane-Decahydrofluorene Formation

Fungal paracyclophane-decahydrofluorene-containing natural products are complex polycyclic metabolites derived from similar hybrid PKS-NRPS pathways. Herein we studied the biosynthesis of pyrrocidines, one representative of this family, by gene inactivation in the producer Sarocladium zeae coupled to thorough metabolic analysis and molecular modeling of key enzymes. We characterized nine pyrrocidines and analogues as well as in mutants a variety of accumulating metabolites with new structures including rare cis-decalin, cytochalasan, and fused 6/15/5 macrocycles. This diversity highlights the extraordinary plasticity of the pyrrocidine biosynthetic gene cluster. From accumulating metabolites, we delineated the scenario of pyrrocidine biosynthesis. The ring A of the decahydrofluorene is installed by PrcB, a membrane-bound cyclizing isomerase, on a PKS-NRPS-derived pyrrolidone precursor. Docking experiments in PrcB allowed us to characterize the active site suggesting a mechanism triggered by arginine-mediated deprotonation at the terminal methyl of the substrate. Next, two integral membrane proteins, PrcD and PrcE, each predicted as a four-helix bundle, perform hydroxylation of the pyrrolidone ring and paracyclophane formation, respectively. Modelization of PrcE highlights a topological homology with vitamin K oxido-reductase and the presence of a disulfide bond. Our results suggest a previously unsuspected coupling mechanism via a transient loss of aromaticity of tyrosine residue to form the strained paracyclophane motif. Finally, the lipocalin-like protein PrcX drives the exo-cycloaddition yielding ring B and C of the decahydrofluorene to afford pyrrocidine A, which is transformed by a reductase PrcI to form pyrrocidine B. These insights will greatly facilitate the microbial production of pyrrocidine analogues by synthetic biology.

Substrate Selection Criteria in Regulated Intramembrane Proteolysis

Alzheimer's disease is the most common form of dementia encountered in an aging population. Characteristic amyloid deposits of Aβ peptides in the brain are generated through cleavage of amyloid precursor protein (APP) by γ-secretase, an intramembrane protease. Cryo-EM structures of substrate γ-secretase complexes revealed details of the process, but how substrates are recognized and enter the catalytic site is still largely ignored. γ-Secretase cleaves a diverse range of substrate sequences without a common consensus sequence, but strikingly, single point mutations within the transmembrane domain (TMD) of specific substrates may greatly affect cleavage efficiencies. Previously, conformational flexibility was hypothesized to be the main criterion for substrate selection. Here we review the 3D structure and dynamics of several γ-secretase substrate TMDs and compare them with mutants shown to affect the cleavage efficiency. In addition, we present structural and dynamic data on ITGB1, a known nonsubstrate of γ-secretase. A comparison of biophysical details between these TMDs and changes generated by introducing crucial mutations allowed us to unravel common principles that differ between substrates and nonsubstrates. We identified three motifs in the investigated substrates: a highly flexible transmembrane domain, a destabilization of the cleavage region, and a basic signature at the end of the transmembrane helix. None of these appears to be exclusive. While conformational flexibility on its own may increase cleavage efficiency in well-known substrates like APP or Notch1, our data suggest that the three motifs seem to be rather variably combined to determine whether a transmembrane helix is efficiently recognized as a γ-secretase substrate.

Nucleolar protein TAAP1/C22orf46 confers pro-survival signaling in non-small cell lung cancer

Tumor cells subvert immune surveillance or lytic stress by harnessing inhibitory signals. Hence, bispecific antibodies have been developed to direct CTLs to the tumor site and foster immune-dependent cytotoxicity. Although applied with success, T cell-based immunotherapies are not universally effective partially because of the expression of pro-survival factors by tumor cells protecting them from apoptosis. Here, we report a CRISPR/Cas9 screen in human non-small cell lung cancer cells designed to identify genes that confer tumors with the ability to evade the cytotoxic effects of CD8+ T lymphocytes engaged by bispecific antibodies. We show that the gene C22orf46 facilitates pro-survival signals and that tumor cells devoid of C22orf46 expression exhibit increased susceptibility to T cell-induced apoptosis and stress by genotoxic agents. Although annotated as a non-coding gene, we demonstrate that C22orf46 encodes a nucleolar protein, hereafter referred to as "Tumor Apoptosis Associated Protein 1," up-regulated in lung cancer, which displays remote homologies to the BH domain containing Bcl-2 family of apoptosis regulators. Collectively, the findings establish TAAP1/C22orf46 as a pro-survival oncogene with implications to therapy.

DeepSS2GO: protein function prediction from secondary structure

Predicting protein function is crucial for understanding biological life processes, preventing diseases and developing new drug targets. In recent years, methods based on sequence, structure and biological networks for protein function annotation have been extensively researched. Although obtaining a protein in three-dimensional structure through experimental or computational methods enhances the accuracy of function prediction, the sheer volume of proteins sequenced by high-throughput technologies presents a significant challenge. To address this issue, we introduce a deep neural network model DeepSS2GO (Secondary Structure to Gene Ontology). It is a predictor incorporating secondary structure features along with primary sequence and homology information. The algorithm expertly combines the speed of sequence-based information with the accuracy of structure-based features while streamlining the redundant data in primary sequences and bypassing the time-consuming challenges of tertiary structure analysis. The results show that the prediction performance surpasses state-of-the-art algorithms. It has the ability to predict key functions by effectively utilizing secondary structure information, rather than broadly predicting general Gene Ontology terms. Additionally, DeepSS2GO predicts five times faster than advanced algorithms, making it highly applicable to massive sequencing data. The source code and trained models are available at https://github.com/orca233/DeepSS2GO.

A phosphatase gene is linked to nectar dihydroxyacetone accumulation in mānuka (Leptospermum scoparium)

Floral nectar composition beyond common sugars shows great diversity but contributing genetic factors are generally unknown. Mānuka (Leptospermum scoparium) is renowned for the antimicrobial compound methylglyoxal in its derived honey, which originates from the precursor, dihydroxyacetone (DHA), accumulating in the nectar. Although this nectar trait is highly variable, genetic contribution to the trait is unclear. Therefore, we investigated key gene(s) and genomic regions underpinning this trait. We used RNAseq analysis to identify nectary-associated genes differentially expressed between high and low nectar DHA genotypes. We also used a mānuka high-density linkage map and quantitative trait loci (QTL) mapping population, supported by an improved genome assembly, to reveal genetic regions associated with nectar DHA content. Expression and QTL analyses both pointed to the involvement of a phosphatase gene, LsSgpp2. The expression pattern of LsSgpp2 correlated with nectar DHA accumulation, and it co-located with a QTL on chromosome 4. The identification of three QTLs, some of the first reported for a plant nectar trait, indicates polygenic control of DHA content. We have established plant genetics as a key influence on DHA accumulation. The data suggest the hypothesis of LsSGPP2 releasing DHA from DHA-phosphate and variability in LsSgpp2 gene expression contributing to the trait variability.

Insights into the inner workings of transformer models for protein function prediction

MOTIVATION

We explored how explainable artificial intelligence (XAI) can help to shed light into the inner workings of neural networks for protein function prediction, by extending the widely used XAI method of integrated gradients such that latent representations inside of transformer models, which were finetuned to Gene Ontology term and Enzyme Commission number prediction, can be inspected too.

RESULTS

The approach enabled us to identify amino acids in the sequences that the transformers pay particular attention to, and to show that these relevant sequence parts reflect expectations from biology and chemistry, both in the embedding layer and inside of the model, where we identified transformer heads with a statistically significant correspondence of attribution maps with ground truth sequence annotations (e.g. transmembrane regions, active sites) across many proteins.

AVAILABILITY AND IMPLEMENTATION

Source code can be accessed at https://github.com/markuswenzel/xai-proteins.

Bioinformatics approach for prediction and analysis of the Non-Structural Protein 4B (NSP4B) of the Zika virus

BACKGROUND

The Nonstructural Protein (NSP) 4B of Zika virus of 251 amino acids from (ZIKV/Human/POLG_ZIKVF) with accession number (A0A024B7W1), Induces the production of Endoplasmic Reticulum ER-derived membrane vesicles, which are the sites of viral replication. To understand the physical basis of how proteins fold in nature and to solve the challenge of protein structure prediction, Ab-initio and comparative modeling are crucial tools.

RESULTS

The systematic in silico technique, ThreaDom, had only predicted one domain (4 - 190) of NSP4B. I-TASSER, and Alphafold were ranked as the best servers for full-length 3-D protein structure predictions of NSP4B, where the predicted models were evaluated quantitatively using benchmarked metrics including C-score (-3.43), TM-score (0.77949), RMSD (2.73), and Z-score (1.561). The functional and protein binding motifs were realized using motif databases, secondary and surface accessibility predictions combined with Post-Translational Modification Sites (PTMs) prediction. Two highly conserved protein-binding motifs (Flavi NS4B and Bacillus papRprotein), together with three (PTMs) (Casein Kinase II, Myristyl site, and ASN-Glycosylation site) were predicted utilizing the Motif scan and Scanprosite servers. These patterns and PTMs were associated with NSP4B's role in triggering the development of the viral replication complex and its participation in the localization of NS3 and NS5 on the membrane. Only one hit from Structural Classification of Protein (SCOP) matched the protein sequence at positions 10 to 397 and was categorized six-hairpin glycosidases superfamily according to CATH (Class, Architecture, Topology, and Homology). Integrating this NSP4B information with the templates' SCOP and CATH annotations achieves it easier to attribute structure-function/evolution links to both previously known and recently discovered protein structures.

Euglena's atypical respiratory chain adapts to the discoidal cristae and flexible metabolism

Euglena gracilis, a model organism of the eukaryotic supergroup Discoba harbouring also clinically important parasitic species, possesses diverse metabolic strategies and an atypical electron transport chain. While structures of the electron transport chain complexes and supercomplexes of most other eukaryotic clades have been reported, no similar structure is currently available for Discoba, limiting the understandings of its core metabolism and leaving a gap in the evolutionary tree of eukaryotic bioenergetics. Here, we report high-resolution cryo-EM structures of Euglena's respirasome I + III2 + IV and supercomplex III2 + IV2. A previously unreported fatty acid synthesis domain locates on the tip of complex I's peripheral arm, providing a clear picture of its atypical subunit composition identified previously. Individual complexes are re-arranged in the respirasome to adapt to the non-uniform membrane curvature of the discoidal cristae. Furthermore, Euglena's conformationally rigid complex I is deactivated by restricting ubiquinone's access to its substrate tunnel. Our findings provide structural insights for therapeutic developments against euglenozoan parasite infections.

Efficient encoding of large antigenic spaces by epitope prioritization with Dolphyn

We investigate a relatively underexplored component of the gut-immune axis by profiling the antibody response to gut phages using Phage Immunoprecipitation Sequencing (PhIP-Seq). To cover large antigenic spaces, we develop Dolphyn, a method that uses machine learning to select peptides from protein sets and compresses the proteome through epitope-stitching. Dolphyn compresses the size of a peptide library by 78% compared to traditional tiling, increasing the antibody-reactive peptides from 10% to 31%. We find that the immune system develops antibodies to human gut bacteria-infecting viruses, particularly E.coli-infecting Myoviridae. Cost-effective PhIP-Seq libraries designed with Dolphyn enable the assessment of a wider range of proteins in a single experiment, thus facilitating the study of the gut-immune axis.

Construction and characterization of a temperature-sensitive pRC4 replicon for Rhodococcus and Gordonia

Temperature-sensitive plasmids are useful for genome engineering and several synthetic biology applications. There are only limited reports on temperature-sensitive plasmids for Rhodococcus and none for Gordonia. Here, we report the construction of a temperature-sensitive pRC4 replicon that is functional in Rhodococcus and Gordonia. The amino acid residues were predicted for the temperature-sensitive phenotype in the pRC4 replicon using in silico methods and molecular simulation of the DNA-binding replication protein with the origin of replication. The amino acid residues were mutated, and the temperature-sensitive phenotype was validated in Gordonia sp. IITR100. Similar results were also observed in Rhodococcus erythropolis, suggesting that the temperature-sensitive phenotype was exhibited across genera.

Serpentoviruses Exhibit Diverse Organization and ORF Composition with Evidence of Recombination

Serpentoviruses are a subfamily of positive sense RNA viruses in the order Nidovirales, family Tobaniviridae, associated with respiratory disease in multiple clades of reptiles. While the broadest viral diversity is reported from captive pythons, other reptiles, including colubrid snakes, turtles, and lizards of captive and free-ranging origin are also known hosts. To better define serpentoviral diversity, eleven novel serpentovirus genomes were sequenced with an Illumina MiSeq and, when necessary, completed with other Sanger sequencing methods. The novel serpentoviral genomes, along with 57 other previously published serpentovirus genomes, were analyzed alongside four outgroup genomes. Genomic analyses included identifying unique genome templates for each serpentovirus clade, as well as analysis of coded protein composition, potential protein function, protein glycosylation sites, differences in phylogenetic history between open-reading frames, and recombination. Serpentoviral genomes contained diverse protein compositions. In addition to the fundamental structural spike, matrix, and nucleoprotein proteins required for virion formation, serpentovirus genomes also included 20 previously uncharacterized proteins. The uncharacterized proteins were homologous to a number of previously characterized proteins, including enzymes, transcription factors, scaffolding, viral resistance, and apoptosis-related proteins. Evidence for recombination was detected in multiple instances in genomes from both captive and free-ranging snakes. These results show serpentovirus as a diverse clade of viruses with genomes that code for a wide diversity of proteins potentially enhanced by recombination events.

PanEffect: a pan-genome visualization tool for variant effects in maize

SUMMARY

Understanding the effects of genetic variants is crucial for accurately predicting traits and functional outcomes. Recent approaches have utilized artificial intelligence and protein language models to score all possible missense variant effects at the proteome level for a single genome, but a reliable tool is needed to explore these effects at the pan-genome level. To address this gap, we introduce a new tool called PanEffect. We implemented PanEffect at MaizeGDB to enable a comprehensive examination of the potential effects of coding variants across 50 maize genomes. The tool allows users to visualize over 550 million possible amino acid substitutions in the B73 maize reference genome and to observe the effects of the 2.3 million natural variations in the maize pan-genome. Each variant effect score, calculated from the Evolutionary Scale Modeling (ESM) protein language model, shows the log-likelihood ratio difference between B73 and all variants in the pan-genome. These scores are shown using heatmaps spanning benign outcomes to potential functional consequences. In addition, PanEffect displays secondary structures and functional domains along with the variant effects, offering additional functional and structural context. Using PanEffect, researchers now have a platform to explore protein variants and identify genetic targets for crop enhancement.

AVAILABILITY AND IMPLEMENTATION

The PanEffect code is freely available on GitHub (https://github.com/Maize-Genetics-and-Genomics-Database/PanEffect). A maize implementation of PanEffect and underlying datasets are available at MaizeGDB (https://www.maizegdb.org/effect/maize/).

Babesia bovis RON2 binds to bovine erythrocytes through a highly conserved epitope

B. bovis invasion of bovine erythrocytes requires tight junction formation involving AMA-1/RON2 complex interaction. RON2 has been considered a vaccine candidate since antibodies targeting the protein can inhibit parasite invasion of target cells; however, the mechanism controlling B. bovis RON2 interaction with red blood cells is not yet fully understood. This study was thus aimed at identifying B. bovis RON2 protein regions associated with interaction with bovine erythrocytes. Natural selection analysis of the ron2 gene identified predominantly negative selection signals in the C-terminal region. Interestingly, protein-cell and competition assays highlighted the RON2-C region's role in peptide 42918-mediated erythrocyte binding, probably to a sialoglycoprotein receptor. This peptide (1218SFIMVKPPALHCVLKPVETL1237) lies within an intrinsically disordered region of the RON2 secondary structure flanked by two helical residues. The study provides, for the first time, valuable insights into RON2's role in interaction with its target cells. Future studies are required for studying the peptide's potential as an anti-B. bovis vaccine component.

Advancements and hurdles in the development of a vaccine for triple-negative breast cancer: A comprehensive review of multi-omics and immunomics strategies

Triple-Negative Breast Cancer (TNBC) presents a significant challenge in oncology due to its aggressive behavior and limited therapeutic options. This review explores the potential of immunotherapy, particularly vaccine-based approaches, in addressing TNBC. It delves into the role of immunoinformatics in creating effective vaccines against TNBC. The review first underscores the distinct attributes of TNBC and the importance of tumor antigens in vaccine development. It then elaborates on antigen detection techniques such as exome sequencing, HLA typing, and RNA sequencing, which are instrumental in identifying TNBC-specific antigens and selecting vaccine candidates. The discussion then shifts to the in-silico vaccine development process, encompassing antigen selection, epitope prediction, and rational vaccine design. This process merges computational simulations with immunological insights. The role of Artificial Intelligence (AI) in expediting the prediction of antigens and epitopes is also emphasized. The review concludes by encapsulating how Immunoinformatics can augment the design of TNBC vaccines, integrating tumor antigens, advanced detection methods, in-silico strategies, and AI-driven insights to advance TNBC immunotherapy. This could potentially pave the way for more targeted and efficacious treatments.

Multiple circulating forms of neprilysin detected with novel epitope-directed monoclonal antibodies

Neprilysin (NEP) is an emerging biomarker for various diseases including heart failure (HF). However, major inter-assay inconsistency in the reported concentrations of circulating NEP and uncertainty with respect to its correlations with type and severity of disease are in part attributed to poorly characterized antibodies supplied in commercial ELISA kits. Validated antibodies with well-defined binding footprints are critical for understanding the biological and clinical context of NEP immunoassay data. To achieve this, we applied in silico epitope prediction and rational peptide selection to generate monoclonal antibodies (mAbs) against spatially distant sites on NEP. One of the selected epitopes contained published N-linked glycosylation sites at N285 and N294. The best antibody pair, mAb 17E11 and 31E1 (glycosylation-sensitive), were characterized by surface plasmon resonance, isotyping, epitope mapping, and western blotting. A validated two-site sandwich NEP ELISA with a limit of detection of 2.15 pg/ml and working range of 13.1-8000 pg/ml was developed with these mAbs. Western analysis using a validated commercial polyclonal antibody (PE pAb) and our mAbs revealed that non-HF and HF plasma NEP circulates as a heterogenous mix of moieties that possibly reflect proteolytic processing, post-translational modifications and homo-dimerization. Both our mAbs detected a ~ 33 kDa NEP fragment which was not apparent with PE pAb, as well as a common ~ 57-60 kDa moiety. These antibodies exhibit different affinities for the various NEP targets. Immunoassay results are dependent on NEP epitopes variably detected by the antibody pairs used, explaining the current discordant NEP measurements derived from different ELISA kits.

DescribePROT in 2023: more, higher-quality and experimental annotations and improved data download options

The DescribePROT database of amino acid-level descriptors of protein structures and functions was substantially expanded since its release in 2020. This expansion includes substantial increase in the size, scope, and quality of the underlying data, the addition of experimental structural information, the inclusion of new data download options, and an upgraded graphical interface. DescribePROT currently covers 19 structural and functional descriptors for proteins in 273 reference proteomes generated by 11 accurate and complementary predictive tools. Users can search our resource in multiple ways, interact with the data using the graphical interface, and download data at various scales including individual proteins, entire proteomes, and whole database. The annotations in DescribePROT are useful for a broad spectrum of studies that include investigations of protein structure and function, development and validation of predictive tools, and to support efforts in understanding molecular underpinnings of diseases and development of therapeutics. DescribePROT can be freely accessed at http://biomine.cs.vcu.edu/servers/DESCRIBEPROT/.

Prop3D: A flexible, Python-based platform for machine learning with protein structural properties and biophysical data

BACKGROUND

Machine learning (ML) has a rich history in structural bioinformatics, and modern approaches, such as deep learning, are revolutionizing our knowledge of the subtle relationships between biomolecular sequence, structure, function, dynamics and evolution. As with any advance that rests upon statistical learning approaches, the recent progress in biomolecular sciences is enabled by the availability of vast volumes of sufficiently-variable data. To be useful, such data must be well-structured, machine-readable, intelligible and manipulable. These and related requirements pose challenges that become especially acute at the computational scales typical in ML. Furthermore, in structural bioinformatics such data generally relate to protein three-dimensional (3D) structures, which are inherently more complex than sequence-based data. A significant and recurring challenge concerns the creation of large, high-quality, openly-accessible datasets that can be used for specific training and benchmarking tasks in ML pipelines for predictive modeling projects, along with reproducible splits for training and testing.

RESULTS

Here, we report 'Prop3D', a platform that allows for the creation, sharing and extensible reuse of libraries of protein domains, featurized with biophysical and evolutionary properties that can range from detailed, atomically-resolved physicochemical quantities (e.g., electrostatics) to coarser, residue-level features (e.g., phylogenetic conservation). As a community resource, we also supply a 'Prop3D-20sf' protein dataset, obtained by applying our approach to CATH . We have developed and deployed the Prop3D framework, both in the cloud and on local HPC resources, to systematically and reproducibly create comprehensive datasets via the Highly Scalable Data Service ( HSDS ). Our datasets are freely accessible via a public HSDS instance, or they can be used with accompanying Python wrappers for popular ML frameworks.

CONCLUSION

Prop3D and its associated Prop3D-20sf dataset can be of broad utility in at least three ways. Firstly, the Prop3D workflow code can be customized and deployed on various cloud-based compute platforms, with scalability achieved largely by saving the results to distributed HDF5 files via HSDS . Secondly, the linked Prop3D-20sf dataset provides a hand-crafted, already-featurized dataset of protein domains for 20 highly-populated CATH families; importantly, provision of this pre-computed resource can aid the more efficient development (and reproducible deployment) of ML pipelines. Thirdly, Prop3D-20sf's construction explicitly takes into account (in creating datasets and data-splits) the enigma of 'data leakage', stemming from the evolutionary relationships between proteins.

Beyond Blast: Enabling Microbiologists to Better Extract Literature, Taxonomic Distributions and Gene Neighborhood Information for Protein Families

Capturing the published corpus of information on all members of a given protein family should be an essential step in any study focusing on specific members of that said family. Using a previously gathered dataset of more than 280 references mentioning a member of the DUF34 (NIF3/Ngg1-interacting Factor 3), we evaluated the efficiency of different databases and search tools, and devised a workflow that experimentalists can use to capture the most published information on members of a protein family in the least amount of time. To complement this workflow, web-based platforms allowing for the exploration of protein family members across sequenced genomes or for the analysis of gene neighborhood information were reviewed for their versatility and ease of use. Recommendations that can be used for experimentalist users, as well as educators, are provided and integrated within a customized, publicly accessible Wiki.

Insights on the nuclear shuttling of H2A-H2B histone chaperones

All cellular processes that involve the unwinding of DNA also lead to the systematic shuttling of histones. Histone shuttling across the nuclear membrane is facilitated by a class of proteins known as - histone chaperones. Histone chaperones are classified based on their binding to H3/H4 histones or H2A/H2B histones. During the shuttling process, two types of signals - NLS and NES are recognized by the nuclear transport proteins. However, this is the nuclear transport protein and the mechanism of signal recognition by the protein is still unknown. Thus, in this piece of work, the NLS and NES signals are predicted on important H2A/H2B binding histone chaperones. In addition, cellular localization and potential DNA binding regions of histone chaperones are predicted. Mapping of predicted regions on the histone chaperone's structure suggested that the critical binding regions mainly lie on the disordered region of the histone chaperones. NLS and NES are present in the N- and C-terminal of the histone chaperones. Most histone chaperones contain bipartiate NLS signals. This article sheds light on the crucial aspect that in addition of being directly engaged in nucleosome synthesis and disassembly in vivo, histone chaperone also performs various specific roles via histone binding activity.

Recombinant p40 Protein Promotes Expression of Occludin in HaCaT Keratinocytes: A Brief Communication

The ability of epithelial barriers to perform as the first defense line against external damage derives from tight junctions, protein complexes that block microorganisms through the paracellular space. Indeed, disturbances of barrier permeability caused by bacterial metabolites and other inflammatory stimuli are the consequence of changes in protein expression in these complexes. Postbiotics, molecules derived from bacteria with beneficial effects on the host, improve barrier function through the activation of survival pathways in epithelial cells. Lacticaseibacillus rhamnosus GG secretes the muramidase p40, which protects intestinal barriers through an EGFR-dependent pathway. In this work, we cloned, expressed, and purified the recombinant p40 protein from L. rhamnosus GR-1 to evaluate its effect on cell viability, cell cytotoxicity, TEER, and protein levels of tight junctions, as well as EGFR activation via Western blot on HaCaT keratinocytes subjected to LPS. We found a novel mutation at residue 368 that does not change the structure of p40. Our protein also reduces the LPS-induced increase in cell cytotoxicity when it is added prior to this stimulus. Furthermore, although LPS did not cause changes in barrier function, p40 increased TEER and occludin expression in HaCaT, but unlike previous work with p40 from LGG, we found that recombinant p40 did not activate EGFR. This suggests that recombinant p40 enhances epithelial barrier function through distinct signaling pathways.

Computational redesign of taxane-10β-hydroxylase for de novo biosynthesis of a key paclitaxel intermediate

Paclitaxel (Taxol®) is the most popular anticancer diterpenoid predominantly present in Taxus. The core skeleton of paclitaxel is highly modified, but researches on the cytochrome P450s involved in post-modification process remain exceedingly limited. Herein, the taxane-10β-hydroxylase (T10βH) from Taxus cuspidata, which is the third post-modification enzyme that catalyzes the conversion of taxadiene-5α-yl-acetate (T5OAc) to taxadiene-5α-yl-acetoxy-10β-ol (T10OH), was investigated in Escherichia coli by combining computation-assisted protein engineering and metabolic engineering. The variant of T10βH, M3 (I75F/L226K/S345V), exhibited a remarkable 9.5-fold increase in protein expression, accompanied by respective 1.3-fold and 2.1-fold improvements in turnover frequency (TOF) and total turnover number (TTN). Upon integration into the engineered strain, the variant M3 resulted in a substantial enhancement in T10OH production from 0.97 to 2.23 mg/L. Ultimately, the titer of T10OH reached 3.89 mg/L by fed-batch culture in a 5-L bioreactor, representing the highest level reported so far for the microbial de novo synthesis of this key paclitaxel intermediate. This study can serve as a valuable reference for further investigation of other P450s associated with the artificial biosynthesis of paclitaxel and other terpenoids. KEY POINTS: • The T10βH from T. cuspidata was expressed and engineered in E. coli unprecedentedly. • The expression and activity of T10βH were improved through protein engineering. • De novo biosynthesis of T10OH was achieved in E. coli with a titer of 3.89 mg/L.

Discovering the deep evolutionary roots of serum amyloid A protein family

Deep evolutionary origin of the conserved animal serum amyloid A (SAA) apolipoprotein family leading to yet unknown highly similar SAA-like sequences occurring in certain bacterial genomes is demonstrated in this contribution. Horizontal gene transfer event of corresponding genes between gut bacteria and non-vertebrate animals was discovered in the reconstructed phylogenetic tree obtained with maximum likelihood and neighbor-joining methods, respectively. This detailed phylogeny based on totally 128 complete sequences comprised diverse serum amyloid A isoforms from various animal vertebrate and non-vertebrate phyla and also corresponding genes coding for highly similar proteins from animal gut bacteria. Typical largely conserved sequence motifs and a peculiar structural fold consisting mainly of four α-helices in a bundle within all reconstructed clades of the SAA protein family are discussed with respect to their supposed biological functions in various organisms that contain corresponding genes.

Microproteins-Discovery, structure, and function

Advances in proteogenomic technologies have revealed hundreds to thousands of translated small open reading frames (sORFs) that encode microproteins in genomes across evolutionary space. While many microproteins have now been shown to play critical roles in biology and human disease, a majority of recently identified microproteins have little or no experimental evidence regarding their functionality. Computational tools have some limitations for analysis of short, poorly conserved microprotein sequences, so additional approaches are needed to determine the role of each member of this recently discovered polypeptide class. A currently underexplored avenue in the study of microproteins is structure prediction and determination, which delivers a depth of functional information. In this review, we provide a brief overview of microprotein discovery methods, then examine examples of microprotein structures (and, conversely, intrinsic disorder) that have been experimentally determined using crystallography, cryo-electron microscopy, and NMR, which provide insight into their molecular functions and mechanisms. Additionally, we discuss examples of predicted microprotein structures that have provided insight or context regarding their function. Analysis of microprotein structure at the angstrom level, and confirmation of predicted structures, therefore, has potential to identify translated microproteins that are of biological importance and to provide molecular mechanism for their in vivo roles.

Harnessing Attenuation-Related Mutations of Viral Genomes: Development of a Serological Assay to Differentiate between Capripoxvirus-Infected and -Vaccinated Animals

Sheeppox, goatpox, and lumpy skin disease caused by the sheeppox virus (SPPV), goatpox virus (GTPV), and lumpy skin disease virus (LSDV), respectively, are diseases that affect millions of ruminants and many low-income households in endemic countries, leading to great economic losses for the ruminant industry. The three viruses are members of the Capripoxvirus genus of the Poxviridae family. Live attenuated vaccines remain the only efficient means for controlling capripox diseases. However, serological tools have not been available to differentiate infected from vaccinated animals (DIVA), though crucial for proper disease surveillance, control, and eradication efforts. We analysed the sequences of variola virus B22R homologue gene for SPPV, GTPV, and LSDV and observed significant differences between field and vaccine strains in all three capripoxvirus species, resulting in the truncation and absence of the B22R protein in major vaccines within each of the viral species. We selected and expressed a protein fragment present in wildtype viruses but absent in selected vaccine strains of all three species, taking advantage of these alterations in the B22R gene. An indirect ELISA (iELISA) developed using this protein fragment was evaluated on well-characterized sera from vaccinated, naturally and experimentally infected, and negative cattle and sheep. The developed wildtype-specific capripox DIVA iELISA showed >99% sensitivity and specificity for serum collected from animals infected with the wildtype virus. To the best of our knowledge, this is the first wildtype-specific, DIVA-capable iELISA for poxvirus diseases exploiting changes in nucleotide sequence alterations in vaccine strains.

The DNA-binding induced (de)AMPylation activity of a Coxiella burnetii Fic enzyme targets Histone H3

The intracellular bacterial pathogen Coxiella burnetii evades the host response by secreting effector proteins that aid in establishing a replication-friendly niche. Bacterial filamentation induced by cyclic AMP (Fic) enzymes can act as effectors by covalently modifying target proteins with the posttranslational AMPylation by transferring adenosine monophosphate (AMP) from adenosine triphosphate (ATP) to a hydroxyl-containing side chain. Here we identify the gene product of C. burnetii CBU_0822, termed C. burnetii Fic 2 (CbFic2), to AMPylate host cell histone H3 at serine 10 and serine 28. We show that CbFic2 acts as a bifunctional enzyme, both capable of AMPylation as well as deAMPylation, and is regulated by the binding of DNA via a C-terminal helix-turn-helix domain. We propose that CbFic2 performs AMPylation in its monomeric state, switching to a deAMPylating dimer upon DNA binding. This study unveils reversible histone modification by a specific enzyme of a pathogenic bacterium.

A putative lipase affects Pseudomonas aeruginosa biofilm matrix production

Pseudomonas aeruginosa is an opportunistic pathogen that is widely known for infecting patients with underlying conditions. This species often survives antibiotic therapy by forming biofilms, in which the cells produce a protective extracellular matrix. P. aeruginosa also produces virulence factors that enhance its ability to cause disease. One signaling pathway that influences virulence is the nitrogen-related phosphotransferase system (Nitro-PTS), which consists of an initial phosphotransferase, PtsP, a phosphocarrier, PtsO, and a terminal phosphate receptor, PtsN. The physiological role of the Nitro-PTS in P. aeruginosa is poorly understood. However, PtsN, when deprived of its upstream phosphotransfer proteins, has an antagonistic effect on biofilm formation. We thus conducted a transposon mutagenesis screen in an unphosphorylated-PtsN (i.e., ∆ptsP) background to identify downstream proteins with unacknowledged roles in PtsN-mediated biofilm suppression. We found an unstudied gene, PA14_04030, whose disruption restored biofilm production. This gene encodes a predicted phospholipase with signature alpha/beta hydrolase folds and a lipase signature motif with an active-site Ser residue. Hence, we renamed the gene bipL, for biofilm-impacting phospholipase. Deletion of bipL in a ∆ptsP background increased biofilm formation, supporting the idea that BipL is responsible for reducing biofilm formation in strains with unphosphorylated PtsN. Moreover, substituting the putative catalytic Ser for Ala phenocopied bipL deletion, indicating that this residue is important for the biofilm-suppressive activity of BipL in vivo. As our preliminary data suggest that BipL is a lipase, we performed lipidomics to detect changes in the lipid profile due to bipL deletion and found changes in some lipid species. IMPORTANCE Biofilm formation by bacteria occurs when cells secrete an extracellular matrix that holds them together and shields them from environmental insults. Biofilms of bacterial opportunistic human pathogens such as Pseudomonas aeruginosa pose a substantial challenge to clinical antimicrobial therapy. Hence, a more complete knowledge about the bacterial factors that influence and regulate production of the biofilm matrix is one key to formulate more effective therapeutic strategies. In this study, we screen for factors that are important for reducing biofilm matrix production in certain genetic backgrounds. We unexpectedly found a gene encoding a putative lipase enzyme and showed that its predicted catalytic site is important for its ability to reduce biofilm formation. Our findings suggest that lipase enzymes have previously uncharacterized functions in biofilm matrix regulation.

An Orphan VrgG Auxiliary Module Related to the Type VI Secretion Systems from Pseudomonas ogarae F113 Mediates Bacterial Killing

The model rhizobacterium Pseudomonas ogarae F113, a relevant plant growth-promoting bacterium, encodes three different Type VI secretion systems (T6SS) in its genome. In silico analysis of its genome revealed the presence of a genetic auxiliary module containing a gene encoding an orphan VgrG protein (VgrG5a) that is not genetically linked to any T6SS structural cluster, but is associated with genes encoding putative T6SS-related proteins: a possible adaptor Tap protein, followed by a putative effector, Tfe8, and its putative cognate immunity protein, Tfi8. The bioinformatic analysis of the VgrG5a auxiliary module has revealed that this cluster is only present in several subgroups of the P. fluorescens complex of species. An analysis of the mutants affecting the vgrG5a and tfe8 genes has shown that the module is involved in bacterial killing. To test whether Tfe8/Tfi8 constitute an effector-immunity pair, the genes encoding Tfe8 and Tfi8 were cloned and expressed in E. coli, showing that the ectopic expression of tfe8 affected growth. The growth defect was suppressed by tfi8 ectopic expression. These results indicate that Tfe8 is a bacterial killing effector, while Tfi8 is its cognate immunity protein. The Tfe8 protein sequence presents homology to the proteins of the MATE family involved in drug extrusion. The Tfe8 effector is a membrane protein with 10 to 12 transmembrane domains that could destabilize the membranes of target cells by the formation of pores, revealing the importance of these effectors for bacterial interaction. Tfe8 represents a novel type of a T6SS effector present in pseudomonads.

Potential for host-symbiont communication via neurotransmitters and neuromodulators in an aneural animal, the marine sponge Amphimedon queenslandica

Interkingdom signalling within a holobiont allows host and symbionts to communicate and to regulate each other's physiological and developmental states. Here we show that a suite of signalling molecules that function as neurotransmitters and neuromodulators in most animals with nervous systems, specifically dopamine and trace amines, are produced exclusively by the bacterial symbionts of the demosponge Amphimedon queenslandica. Although sponges do not possess a nervous system, A. queenslandica expresses rhodopsin class G-protein-coupled receptors that are structurally similar to dopamine and trace amine receptors. When sponge larvae, which express these receptors, are exposed to agonists and antagonists of bilaterian dopamine and trace amine receptors, we observe marked changes in larval phototactic swimming behaviour, consistent with the sponge being competent to recognise and respond to symbiont-derived trace amine signals. These results indicate that monoamines synthesised by bacterial symbionts may be able to influence the physiology of the host sponge.

Combined prediction and design reveals the target recognition mechanism of an intrinsically disordered protein interaction domain

An increasing number of protein interaction domains and their targets are being found to be intrinsically disordered proteins (IDPs). The corresponding target recognition mechanisms are mostly elusive because of challenges in performing detailed structural analysis of highly dynamic IDP-IDP complexes. Here, we show that by combining recently developed computational approaches with experiments, the structure of the complex between the intrinsically disordered C-terminal domain (CTD) of protein 4.1G and its target IDP region in NuMA can be dissected at high resolution. First, we carry out systematic mutational scanning using dihydrofolate reductase-based protein complementarity analysis to identify essential interaction regions and key residues. The results are found to be highly consistent with an α/β-type complex structure predicted by AlphaFold2 (AF2). We then design mutants based on the predicted structure using a deep learning protein sequence design method. The solved crystal structure of one mutant presents the same core structure as predicted by AF2. Further computational prediction and experimental assessment indicate that the well-defined core structure is conserved across complexes of 4.1G CTD with other potential targets. Thus, we reveal that an intrinsically disordered protein interaction domain uses an α/β-type structure module formed through synergistic folding to recognize broad IDP targets. Moreover, we show that computational prediction and experiment can be jointly applied to segregate true IDP regions from the core structural domains of IDP-IDP complexes and to uncover the structure-dependent mechanisms of some otherwise elusive IDP-IDP interactions.

Identification of an Antimicrobial Peptide from the Venom of the Trinidad Thick-Tailed Scorpion Tityus trinitatis with Potent Activity against ESKAPE Pathogens and Clostridioides difficile

Envenomation by the Trinidad thick-tailed scorpion Tityus trinitatis may result in fatal myocarditis and there is a high incidence of acute pancreatitis among survivors. Peptidomic analysis (reversed-phase HPLC followed by MALDI-TOF mass spectrometry and automated Edman degradation) of T. trinitatis venom led to the isolation and characterization of three peptides with antimicrobial activity. Their primary structures were established asTtAP-1 (FLGSLFSIGSKLLPGVFKLFSRKKQ.NH2), TtAP-2 (IFGMIPGLIGGLISAFK.NH2) and TtAP-3 (FFSLIPSLIGGLVSAIK.NH2). In addition, potassium channel and sodium channel toxins, present in the venom in high abundance, were identified by CID-MS/MS sequence analysis. TtAP-1 was the most potent against a range of clinically relevant Gram-positive and Gram-negative aerobes and against the anaerobe Clostridioides difficile (MIC = 3.1-12.5 µg/mL). At a concentration of 1× MIC, TtAP-1 produced rapid cell death (<15 min against Acinetobacter baumannii and Staphylococcus aureus). The therapeutic potential of TtAP-1 as an anti-infective agent is limited by its high hemolytic activity (LC50 = 18 µg/mL against mouse erythrocytes) but the peptide constitutes a template for the design of analogs that maintain the high bactericidal activity against ESKAPE pathogens but are less toxic to human cells. It is suggested that the antimicrobial peptides in the scorpion venom facilitate the action of the neurotoxins by increasing the membrane permeability of cells from either prey or predator.

Serine protease 35 regulates the fibroblast matrisome in response to hyperosmotic stress

Hyperosmotic stress occurs in several diseases, but its long-term effects are largely unknown. We used sorbitol-treated human fibroblasts in 3D culture to study the consequences of hyperosmotic stress in the skin. Sorbitol regulated many genes, which help cells cope with the stress condition. The most robustly regulated gene encodes serine protease 35 (PRSS35). Its regulation by hyperosmotic stress was dependent on the kinases p38 and JNK and the transcription factors NFAT5 and ATF2. We identified different collagens and collagen-associated proteins as putative PRSS35 binding partners. This is functionally important because PRSS35 affected the extracellular matrix proteome, which limited cell proliferation. The in vivo relevance of these findings is reflected by the coexpression of PRSS35 and its binding partners in human skin wounds, where hyperosmotic stress occurs as a consequence of excessive water loss. These results identify PRSS35 as a key regulator of the matrisome under hyperosmotic stress conditions.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance

PURPOSE

Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis.

METHODS

Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins.

RESULTS

Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis.

CONCLUSION

This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

Bioinformatic analysis of SIRT7 sequence and structure

Sirtuins are highly conserved proteins that perform very important functions in different cellular processes. Notably, SIRT7 is the least studied human sirtuin, but it is known to be involved in a wide variety of processes in both health and disease. In this way, SIRT7 activity-regulating molecules could be beneficial for the treatment of relevant diseases such as cardiovascular and bone diseases, where SIRT7 levels are reduced, or obesity and cancer, where they are increased. In this work, using bioinformatic methods, the sequence and structure of SIRT7 orthologs in a wide variety of organisms were analyzed. Thus, the catalytic domain was found to be quite conserved (83.23% identity) and key residues such as D118, Y119, R120, D170, H187, N189, C198, C225, C228, V273, G298, F239 and V237 were identified. Furthermore, a phylogenetic tree was constructed where SIRT7 orthologs from mammals, birds, reptiles, amphibians, fish, insects, and arachnids were found to cluster in different groups. Finally, predicted three-dimensional structures showed a classic structure of the central catalytic region of most sirtuins, while the flanking N- and C-terminal regions were unique to each phylogenetic group. All this helps to understand a little more how SIRT7 works and gives clues for the future design and development of small molecules that benefit human and animal health.Communicated by Ramaswamy H. Sarma.

A unique mRNA decapping complex in trypanosomes

Removal of the mRNA 5' cap primes transcripts for degradation and is central for regulating gene expression in eukaryotes. The canonical decapping enzyme Dcp2 is stringently controlled by assembly into a dynamic multi-protein complex together with the 5'-3'exoribonuclease Xrn1. Kinetoplastida lack Dcp2 orthologues but instead rely on the ApaH-like phosphatase ALPH1 for decapping. ALPH1 is composed of a catalytic domain flanked by C- and N-terminal extensions. We show that T. brucei ALPH1 is dimeric in vitro and functions within a complex composed of the trypanosome Xrn1 ortholog XRNA and four proteins unique to Kinetoplastida, including two RNA-binding proteins and a CMGC-family protein kinase. All ALPH1-associated proteins share a unique and dynamic localization to a structure at the posterior pole of the cell, anterior to the microtubule plus ends. XRNA affinity capture in T. cruzi recapitulates this interaction network. The ALPH1 N-terminus is not required for viability in culture, but essential for posterior pole localization. The C-terminus, in contrast, is required for localization to all RNA granule types, as well as for dimerization and interactions with XRNA and the CMGC kinase, suggesting possible regulatory mechanisms. Most significantly, the trypanosome decapping complex has a unique composition, differentiating the process from opisthokonts.

The intrinsically disordered cytoplasmic tail of a dendrite branching receptor uses two distinct mechanisms to regulate the actin cytoskeleton

Dendrite morphogenesis is essential for neural circuit formation, yet the molecular mechanisms underlying complex dendrite branching remain elusive. Previous studies on the highly branched Caenorhabditis elegans PVD sensory neuron identified a membrane co-receptor complex that links extracellular signals to intracellular actin remodeling machinery, promoting high-order dendrite branching. In this complex, the claudin-like transmembrane protein HPO-30 recruits the WAVE regulatory complex (WRC) to dendrite branching sites, stimulating the Arp2/3 complex to polymerize actin. We report here our biochemical and structural analysis of this interaction, revealing that the intracellular domain (ICD) of HPO-30 is intrinsically disordered and employs two distinct mechanisms to regulate the actin cytoskeleton. First, HPO-30 ICD binding to the WRC requires dimerization and involves the entire ICD sequence, rather than a short linear peptide motif. This interaction enhances WRC activation by the GTPase Rac1. Second, HPO-30 ICD directly binds to the sides and barbed end of actin filaments. Binding to the barbed end requires ICD dimerization and inhibits both actin polymerization and depolymerization, resembling the actin capping protein CapZ. These dual functions provide an intriguing model of how membrane proteins can integrate distinct mechanisms to fine-tune local actin dynamics.

Crystal structure of Prp16 in complex with ADP

DEAH-box helicases play a crucial role in pre-mRNA splicing as they are responsible for major rearrangements of the spliceosome and are involved in various quality-ensuring steps. Prp16 is the driving force during spliceosomal catalysis, remodeling the C state into the C* state. Here, the first crystal structure of Prp16 from Chaetomium thermophilum in complex with ADP is reported at 1.9 Å resolution. Comparison with the other spliceosomal DEAH-box helicases Prp2, Prp22 and Prp43 reveals an overall identical domain architecture. The β-hairpin, which is a structural element of the RecA2 domain, exhibits a unique position, punctuating its flexibility. Analysis of cryo-EM models of spliceosomal complexes containing Prp16 reveals that these models show Prp16 in its nucleotide-free state, rendering the model presented here the first structure of Prp16 in complex with a nucleotide.

Coordination of the N-Terminal Heme in the Non-Classical Peroxidase from Escherichia coli

The non-classical bacterial peroxidase from Escherichia coli, YhjA, is proposed to deal with peroxidative stress in the periplasm when the bacterium is exposed to anoxic environments, defending it from hydrogen peroxide and allowing it to thrive under those conditions. This enzyme has a predicted transmembrane helix and is proposed to receive electrons from the quinol pool in an electron transfer pathway involving two hemes (NT and E) to accomplish the reduction of hydrogen peroxide in the periplasm at the third heme (P). Compared with classical bacterial peroxidases, these enzymes have an additional N-terminal domain binding the NT heme. In the absence of a structure of this protein, several residues (M82, M125 and H134) were mutated to identify the axial ligand of the NT heme. Spectroscopic data demonstrate differences only between the YhjA and YhjA M125A variant. In the YhjA M125A variant, the NT heme is high-spin with a lower reduction potential than in the wild-type. Thermostability was studied by circular dichroism, demonstrating that YhjA M125A is thermodynamically more unstable than YhjA, with a lower T_M (43 °C vs. 50 °C). These data also corroborate the structural model of this enzyme. The axial ligand of the NT heme was validated to be M125, and mutation of this residue was proven to affect the spectroscopic, kinetic, and thermodynamic properties of YhjA.

In silico characterization of the psilocybin biosynthesis pathway

Nearly all mushrooms of the Psilocybe genus contain the natural product psilocybin, which is a psychoactive alkaloid derived from l-tryptophan. Considering their use in ancient times, as well as their psychedelic properties, these mushrooms have re-emerged with psychotherapeutic potential for treating depression, which has triggered increased pharmaceutical interest. However, the psilocybin biosynthesis pathway was only recently defined and, as such, little exists in the way of structural data. Accordingly, the aim of this study was to structurally characterize this pathway by generating homology models for the four Psilocybe cubensis enzymes involved in psilocybin biosynthesis (PsiD, a decarboxylase; PsiH, a monooxygenase; PsiK, a phosphotransferase; PsiM, a methyltransferase). Following initial model generation and alignment with the identified structural templates, repeated refinement of the models was carried out using secondary structure prediction, geometry evaluation, energy minimization, and molecular dynamics simulations in water. The final models were then evaluated using molecular docking interactions with their substrates, i.e., psilocybin precursors (l-tryptophan, tryptamine, 4-hydroxytryptamine, and norbaeocystin/baeocystin), all of which generated feasible binding modes for the expected biotransformation. Further plausibility of the psilocybin → aeruginascin, 4-hydroxytryptamine → norpsilocin, and tryptamine → N,N-dimethyltryptamine conversions, all mediated by the generated model for PsiM, suggests valid routes of formation for these key secondary metabolites. The structural characterization of these enzymes and their binding modes which emerged from this study can lead to a better understanding of psilocybin synthesis, thereby paving the way for the development of novel substrates and selective inhibitors, as well as improved biotechnological manipulation and production of psilocybin in vitro.

Structures of Tetrahymena thermophila respiratory megacomplexes on the tubular mitochondrial cristae

Tetrahymena thermophila, a classic ciliate model organism, has been shown to possess tubular mitochondrial cristae and highly divergent electron transport chain involving four transmembrane protein complexes (I-IV). Here we report cryo-EM structures of its ~8 MDa megacomplex IV₂+ (I + III₂+ II)₂, as well as a ~ 10.6 MDa megacomplex (IV₂ + I + III₂+ II)₂ at lower resolution. In megacomplex IV₂+ (I + III₂+ II)₂, each CIV₂ protomer associates one copy of supercomplex I + III₂ and one copy of CII, forming a half ring-shaped architecture that adapts to the membrane curvature of mitochondrial cristae. Megacomplex (IV₂+ I + III₂+ II)₂ defines the relative position between neighbouring half rings and maintains the proximity between CIV₂ and CIII₂ cytochrome c binding sites. Our findings expand the current understanding of divergence in eukaryotic electron transport chain organization and how it is related to mitochondrial morphology.

Structure and function characterization of the α-L-arabinofuranosidase from the white-rot fungus Trametes hirsuta

α-L-Arabinofuranosidases (Abfs) play a crucial role in the degradation of hemicelluloses, especially arabinoxylans (AX). Most of the available characterized Abfs are from bacteria, while fungi, as natural decomposers, contain Abfs with little attention given. An arabinofuranosidase (ThAbf1), belonging to the glycoside hydrolase 51 (GH51) family, from the genome of the white-rot fungus Trametes hirsuta, was recombinantly expressed, characterized, and functionally determined. The general biochemical properties showed that the optimal conditions for ThAbf1 were pH 6.0 and 50°C. In substrate kinetics assays, ThAbf1 preferred small fragment arabinoxylo-oligosaccharides (AXOS) and could surprisingly hydrolyze di-substituted 2³,3³-di-L-arabinofuranosyl-xylotriose (A^2,3XX). It also synergized with commercial xylanase (XYL) and increased the saccharification efficiency of arabinoxylan. The crystal structure of ThAbf1 indicated the presence of an adjacent cavity next to the catalytic pocket which led to the ability of ThAbf1 to degrade di-substituted AXOS. The narrow binding pocket prevents ThAbf1 from binding larger substrates. These findings have strengthened our understanding of the catalytic mechanism of GH51 family Abfs and provided a theoretical foundation for the development of more efficient and versatile Abfs to accelerate the degradation and biotransformation of hemicellulose in biomass. KEY POINTS: • ThAbf1 from Trametes hirsuta degraded di-substituted arabinoxylo-oligosaccharide. • ThAbf1 performed detailed biochemical characterization and kinetics. • ThAbf1 structure has been obtained to illustrate the substrate specificity.

Protein-Gene Orthology in Baculoviridae: An Exhaustive Analysis to Redefine the Ancestrally Common Coding Sequences

Baculoviruses are entomopathogens that carry large, double-stranded circular DNA genomes and infect insect larvae of Lepidoptera, Hymenoptera and Diptera, with applications in the biological control of agricultural pests, in the production of recombinant proteins and as viral vectors for various purposes in mammals. These viruses have a variable genetic composition that differs between species, with some sequences shared by all known members, and others that are lineage-specific or unique to isolates. Based on the analysis of nearly 300 sequenced genomes, a thorough bioinformatic investigation was conducted on all the baculoviral protein coding sequences, characterizing their orthology and phylogeny. This analysis confirmed the 38 protein coding sequences currently considered as core genes, while also identifying novel coding sequences as candidates to join this set. Accordingly, homology was found among all the major occlusion body proteins, thus proposing that the polyhedrin, granulin and CUN085 genes be considered as the 39th core gene of Baculoviridae.

Bioinformatics characteristics and expression analysis of TLR3 and its adaptor protein TRIF in Largemouth bass (Micropterus salmoides) upon Flavobacterium columnare infection

TLR3 and TRIF (adaptor protein for TLR3) are vital to the MyD88-independent pathway mediated by Toll-like receptors (TLRs). In order to identify the role of TLR3 and TRIF in Micropterus salmoides, the Ms_TLR3 and Ms_TRIF (Ms: abbreviation for M. salmoides) were cloned and characterized in this study. The open reading frames (ORFs) of Ms_TLR3 and Ms_TRIF genes were 2736 bp and 1791 bp in length, encoding 911 and 596 amino acids, respectively. The protein structure of Ms_TLR3 includes a signal peptide, 18 LRR-related domains, a low complexity region, a transmembrane region, and a TIR domain. However, only a TIR domain and a coiled coil domain were found in Ms_TRIF. Both Ms_TLR3 and Ms_TRIF showed the highest homology to that of M. dolomieu. Ms_TLR3 and Ms_TRIF showed similar expression patterns in various tissues, with the highest expression level in the head kidney. After stimulation of Flavobacterium columnare, the mRNA expressions of Ms_TLR3 and Ms_TRIF were significantly up-regulated at 1 dpi in the gill, spleen and head kidney, and at 6 hpi in the trunk kidney. Furthermore, morphological changes in the gills of largemouth bass challenged with F. columnare suggested that F. columnare infection can destroy the gill filament. Taken together, Ms_TLR3 and Ms_TRIF are indeed involved in F. columnare infection and the subsequent immune response in largemouth bass. Moreover, Ms_TLR3 and Ms_TRIF might respectively play their potential roles in mucosal (mainly in the gill) and systemic (mainly in the head kidney) immune response to bacterial infection.