Sequence Conservation Research Articles

A simple BLASTn-based approach generates novel insights into the regulation and biological function of type I toxin-antitoxins.

Bacterial chromosomal type I toxin-antitoxin systems consist of a small protein, typically under 60 amino acids, and a small RNA (sRNA) that represses toxin translation. These gene pairs have gained attention over the last decade for their contribution to antibiotic persistence and phage tolerance in bacteria. However, biological functions for many remain elusive as gene deletions often fail to produce an observable phenotype. For many pairs, it is still unknown when the toxin and/or antitoxin gene are natively expressed within the bacterium. We examined sequence conservation of three type I toxin-antitoxin systems, tisB/istR-1, shoB/ohsC, and zor/orz, in over 2,000 Escherichia coli strains, including pathogenic and commensal isolates. Using our custom database, we found that these gene pairs are widespread across E. coli and have expression potential via BLASTn. We identified an alternative, dominant sequence variant of TisB and confirmed that it is toxic upon overproduction. Additionally, analyses revealed a highly conserved sequence in the zorO mRNA untranslated region that is required for full toxicity. We further noted that over 30% of E. coli genomes contain an orz antitoxin gene only and confirmed its expression in a representative strain: the first confirmed report of a type I antitoxin without its cognate toxin. Our results add to our understanding of these systems, and our methodology is applicable for other type I loci to identify critical regulatory and functional features.IMPORTANCEChromosomal type I toxin-antitoxins are a class of genes that have gained increasing attention over the last decade for their roles in antibiotic persistence which may contribute to therapeutic failures. However, the control of many of these genes and when they function have remained elusive. We demonstrate that a simple genetic conservation-based approach utilizing free, publicly available data yields known and novel insights into the regulation and function of three chromosomal type I toxin-antitoxins in Escherichia coli. This study also provides a framework for how this approach could be applied to other genes of interest.

Phylogenetic insight into ABCE gene subfamily in plants.

ATP-BINDING CASSETTE SUBFAMILY E MEMBER (ABCE) proteins are one of the most conserved proteins across eukaryotes and archaea. Yeast and most animals possess a single ABCE gene encoding the critical translational factor ABCE1. In several plant species, including Arabidopsis thaliana and Oryza sativa, two or more ABCE gene copies have been identified, however information related to plant ABCE gene family is still missing. In this study we retrieved ABCE gene sequences of 76 plant species from public genome databases and comprehensively analyzed them with the reference to A. thaliana ABCE2 gene (AtABCE2). Using bioinformatic approach we assessed the conservation and phylogeny of plant ABCEs. In addition, we performed haplotype analysis of AtABCE2 and its paralogue AtABCE1 using genomic sequences of 1,135 A. thaliana ecotypes. Plant ABCE proteins showed overall high sequence conservation, sharing at least 78% of amino acid sequence identity with AtABCE2. We found that over half of the selected species have two to eight ABCE genes, suggesting that in plants ABCE genes can be classified as a low-copy gene family, rather than a single-copy gene family. The phylogenetic trees of ABCE protein sequences and the corresponding coding sequences demonstrated that Brassicaceae and Poaceae families have independently undergone lineage-specific split of the ancestral ABCE gene. Other plant species have gained ABCE gene copies through more recent duplication events. We also noticed that ploidy level but not ancient whole genome duplications experienced by a species impacts ABCE gene family size. Deeper analysis of AtABCE2 and AtABCE1 from 1,135 A. thaliana ecotypes revealed four and 35 non-synonymous SNPs, respectively. The lower natural variation in AtABCE2 compared to AtABCE1 is in consistence with its crucial role for plant viability. Overall, while the sequence of the ABCE protein family is highly conserved in the plant kingdom, many plants have evolved to have more than one copy of this essential translational factor.

Molecular epidemiology of the Old World screwworm fly (OWSF) in Iraq, and the genetic structure of various OWSF populations worldwide

Despite being endemic in Iraq, no reports have been published in the past 10 years to update the molecular epidemiology of the Old World screwworm fly (OWSF), Chrysomya bezziana, in this country. In the present study, 130 sheep from 10 Iraqi governorates were found infected with C. bezziana larvae, whose identities were PCR-confirmed based on the cytochrome b (Cytb) gene, and 23 isolates from various tested governorates were successfully sequenced. Although most isolates (n = 20) belonged to the common haplotype circulating in Iraq, two new haplotypes were detected. Significant changes in OWSF epidemiology in Iraq were also suggested, since infestations were detected, for the first time, in Nineveh governorate. Isolates of the present study were combined to those previously published from Iraq and worldwide, collected after searching the GenBank, and various genetic and population structure analyses were conducted. These isolates displayed a great statistically significant value when tested for the purifying (negative) selection, suggesting the limited occurrence of genetic variations, which was evidenced by the high sequence conservation (C = 0.937) value detected. A few isolates from Africa were revealed during our search, and clustered in a separate lineage other than that of the Asian isolates. The latter displayed different genetic variation patterns when compared. For example, isolates from geographically separate regions, e.g., the Gulf Arab countries and South-Eastern Asia had marked genetic differences. On the other hand, isolates from regions with close geographic proximity (the Gulf Arab countries and Iran) had limited genetic subdivision. This is not the case when comparing isolates from 10 islands in the Indonesian Archipelago. Populations from Sumatra and Sumba were isolated and displayed high genetic variations toward the other populations. On the contrary, populations from Sulawesi, Lombok and Sumbawa displayed limited genetic variations. This is particularly important, since it can help detecting the dynamics of establishing the sterile insect technique over various regions as an effective control strategy against the OWSFs.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation.

Intrinsically disordered regions (IDRs) are critical for cellular function yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. Here, we focus on the first two domains of coronavirus nucleocapsid proteins: the disordered N-terminal domain (NTD) and the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system for exploring the interplay between an IDR adjacent to a folded domain and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved dynamic and disordered RNA:protein complex is found across nucleocapsid protein orthologs despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

The pan-plastome of Prunus mume: insights into Prunus diversity, phylogeny, and domestication history.

Prunus mume in the Rosaceae and commonly referred to as mei or Chinese plum is widely used as a traditional ornamental flowering plant and fruit tree in China. Although some population and genetic analyses have been conducted for this species, no extensive comparisons of genetic variation from plastomes have yet been investigated. We de novo assembled a total of 322 complete P. mume plastomes in this study and did a series of comparative analyses to better resolve pan-plastomic patterns of P. mume. To determine the phylogeny and domestication history of this species, we reconstructed the phylogenetic tree of Prunus genus, and resolved the population structure of P. mume. We also examined the nucleotide variation of P. mume to find potential DNA barcodes. The assembled plastomes exhibited a typical quadripartite structure and ranged from 157,871 bp to 158,213 bp in total size with a GC content ranging from 36.73 to 36.75%. A total of 112 unique genes were identified. Single nucleotide variants (SNVs) were the most common variants found among the plastomes, followed by nucleotide insertions/deletions (InDels), and block substitutions with the intergenic spacer (IGS) regions containing the greatest number of variants. From the pan-plastome data six well-supported genetic clusters were resolved using multiple different population structure analyses. The different cultivars were unevenly distributed among multiple clades. We also reconstructed a phylogeny for multiple species of Prunus to better understand genus level diversity and history from which a complex introgressive relationship between mei and other apricots/plums was resolved. This study constructed the pan-plastome of P. mume, which indicated the domestication of P. mume involved multiple genetic origins and possible matrilineal introgression from other species. The phylogenetic analysis in Prunus and the population structure of P. mume provide an important maternal history for Prunus and the groundwork for future studies on intergenomic sequence transfers, cytonuclear incompatibility, and conservation genetics.

KDM4B mutations in human cancers

Homologous recombination (HR) is essential for repair of DNA double-strand breaks (DSBs) and restart of stalled or collapsed replication forks. Most cancers are characterized by mutations in components of the DSB repair pathways. Redundant DSB repair pathways exist in eukaryotes from yeast to humans and recent evidence has shown that complete loss of HR function appears to be lethal. Recent evidence has also shown that cancer cells with mutations in one DSB repair pathway can be killed by inhibiting one or more parallel pathways, a strategy that is currently aggressively explored as a cancer therapy. KDM4B is a histone demethylase with pleiotropic functions, which participates in preparing DSBs for repair by contributing to chromatin remodeling. In this report we carried out a pan-cancer analysis of KDM4B mutations with the goal of understanding their distribution and interaction with other DSB genes. We find that although KDM4B mutations co-occur with DSB repair genes, most KDM4B mutations are not drivers or pathogenic. A sequence conservation analysis from yeast to humans shows that highly conserved residues are resistant to mutation. Finally, all mutations occur in a heterozygous state. A single mutation, R986L, was predicted to significantly affect protein structure using computational modeling. This analysis suggests that KDM4B makes contributions to DSB repair but is not a key player.

Role of Rabenosyn-5 and Rab5b in host cell cytosol uptake reveals conservation of endosomal transport in malaria parasites.

Vesicular trafficking, including secretion and endocytosis, plays fundamental roles in the unique biology of Plasmodium falciparum blood-stage parasites. Endocytosis of host cell cytosol (HCC) provides nutrients and room for parasite growth and is critical for the action of antimalarial drugs and parasite drug resistance. Previous work showed that PfVPS45 functions in endosomal transport of HCC to the parasite's food vacuole, raising the possibility that malaria parasites possess a canonical endolysosomal system. However, the seeming absence of VPS45-typical functional interactors such as rabenosyn 5 (Rbsn5) and the repurposing of Rab5 isoforms and other endolysosomal proteins for secretion in apicomplexans question this idea. Here, we identified a parasite Rbsn5-like protein and show that it functions with VPS45 in the endosomal transport of HCC. We also show that PfRab5b but not PfRab5a is involved in the same process. Inactivation of PfRbsn5L resulted in PI3P and PfRab5b decorated HCC-filled vesicles, typical for endosomal compartments. Overall, this indicates that despite the low sequence conservation of PfRbsn5L and the unusual N-terminal modification of PfRab5b, principles of endosomal transport in malaria parasite are similar to that of model organisms. Using a conditional double protein inactivation system, we further provide evidence that the PfKelch13 compartment, an unusual apicomplexa-specific endocytosis structure at the parasite plasma membrane, is connected upstream of the Rbsn5L/VPS45/Rab5b-dependent endosomal route. Altogether, this work indicates that HCC uptake consists of a highly parasite-specific part that feeds endocytosed material into an endosomal system containing more canonical elements, leading to the delivery of HCC to the food vacuole.

Genomic Analysis of Brassinosteroid Biosynthesis Gene Family Reveals Its Roles in Cotton Development across Gossypium Species.

Cotton is a globally significant economic crop. Brassinosteroids (BRs) are crucial to cotton development. This study systematically analyzed the BR synthase gene family in four cotton species and identified 60 BR genes: 20 in Gossypium hirsutum (GhBRs), 20 in G. barbadense (GbBRs), 10 in G. arboreum (GaBRs), and 10 in G. raimondii (GrBRs). The analysis was extended to chromosomal localization, evolutionary relationships, domain features, and cis-regulatory elements in the promoter regions of BR synthase genes. The results showed that the BR synthase genes were evenly distributed across different subgenomes and chromosomes. Bioinformatics analyses revealed high conservation of amino acid sequences, secondary structures, and conserved domains among the subfamily members, which is closely linked to their pivotal roles in the BR biosynthesis pathway. Cis-element distribution analysis of the BR synthase genes further underscored the complexity of BR gene expression regulation, which is influenced by multiple factors, including plant hormones, abiotic stress, and transcription factors. Expression profiling of GhBRs genes in various cotton tissues and developmental stages highlighted the key roles of GhROT3-1 and GhDET2-1 in fiber elongation and initiation, respectively. Protein-protein interactions and transcription factor analyses further elucidated the regulatory mechanisms of GhROT3-1 and GhDET2-1 in cotton growth and development. This study lays a theoretical foundation for understanding the role of the BR signaling pathway in cotton development, facilitating molecular breeding.

MEM-based pangenome indexing for k-mer queries.

Pangenomes are growing in number and size, thanks to the prevalence of high-quality long-read assemblies. However, current methods for studying sequence composition and conservation within pangenomes have limitations. Methods based on graph pangenomes require a computationally expensive multiple-alignment step, which can leave out some variation. Indexes based on k-mers and de Bruijn graphs are limited to answering questions at a specific substring length k. We present Maximal Exact Match Ordered (MEMO), a pangenome indexing method based on maximal exact matches (MEMs) between sequences. A single MEMO index can handle arbitrary-length queries over pangenomic windows. MEMO enables both queries that test k-mer presence/absence (membership queries) and that count the number of genomes containing k-mers in a window (conservation queries). MEMO's index for a pangenome of 89 human autosomal haplotypes fits in 2.04 GB, 8.8× smaller than a comparable KMC3 index and 11.4× smaller than a PanKmer index. MEMO indexes can be made smaller by sacrificing some counting resolution, with our decile-resolution HPRC index reaching 0.67 GB. MEMO can conduct a conservation query for 31-mers over the human leukocyte antigen locus in 13.89 seconds, 2.5x faster than other approaches. MEMO's small index size, lack of k-mer length dependence, and efficient queries make it a flexible tool for studying and visualizing substring conservation in pangenomes.

Bioinformatic prediction of proteins relevant to functions of the bacterial OLE ribonucleoprotein complex.

OLE (ornate, large, extremophilic) RNAs are members of a noncoding RNA class present in many Gram-positive, extremophilic bacteria. The large size, complex structure, and extensive sequence conservation of OLE RNAs are characteristics consistent with the hypothesis that they likely function as ribozymes. The OLE RNA representative from Halalkalibacterium halodurans is known to localize to the phospholipid membrane and requires at least three essential protein partners: OapA, OapB, and OapC. However, the precise biochemical functions of this unusual ribonucleoprotein (RNP) complex remain unknown. Genetic disruption of OLE RNA or its partners revealed that the complex is beneficial under diverse stress conditions. To search for additional links between OLE RNA and other cellular components, we used phylogenetic profiling to identify proteins that are either correlated or anticorrelated with the presence of OLE RNA in various bacterial species. This analysis revealed strong correlations between the essential protein-binding partners of OLE RNA and organisms that carry the ole gene. Similarly, proteins involved in sporulation are correlated, suggesting a potential role for the OLE RNP complex in spore formation. Intriguingly, the Mg2+ transporter MpfA is strongly anticorrelated with OLE RNA. Evidence indicates that MpfA is structurally related to OapA and therefore MpfA may serve as a functional replacement for some contributions otherwise performed by the OLE RNP complex in species that lack this device. Indeed, OLE RNAs might represent an ancient RNA class that enabled primitive organisms to sense and respond to major cellular stresses.IMPORTANCEOLE (ornate, large, extremophilic) RNAs were first reported nearly 20 years ago, and they represent one of the largest and most intricately folded noncoding RNA classes whose biochemical function remains to be established. Other RNAs with similar size, structural complexity, and extent of sequence conservation have proven to catalyze chemical transformations. Therefore, we speculate that OLE RNAs likewise operate as ribozymes and that they might catalyze a fundamental reaction that has persisted since the RNA World era-a time before the emergence of proteins in evolution. To seek additional clues regarding the function of OLE RNA, we undertook a computational effort to identify potential protein components of the OLE ribonucleoprotein (RNP) complex or other proteins that have functional links to this device. This analysis revealed known protein partners and several additional proteins that might be physically or functionally linked to the OLE RNP complex. Finally, we identified a Mg2+ transporter protein, MpfA, that strongly anticorrelates with the OLE RNP complex. This latter result suggests that MpfA might perform at least some functions that are like those carried out by the OLE RNP complex.

Genotyping Hepatitis B virus by Next-Generation Sequencing: Detection of Mixed Infections and Analysis of Sequence Conservation.

Our aim was to develop an accurate, highly sensitive method for HBV genotype determination and detection of genotype mixtures. We examined the preS and 5' end of the HBV X gene (5X) regions of the HBV genome using next-generation sequencing (NGS). The 1852 haplotypes obtained were subjected to genotyping via the Distance-Based discrimination method (DB Rule) using two sets of 95 reference sequences of genotypes A-H. In clinical samples from 125 patients, the main genotypes were A, D, F and H in Caucasian, B and C in Asian and A and E in Sub-Saharan patients. Genotype mixtures were identified in 28 (22.40%) cases, and potential intergenotypic recombination was observed in 29 (23.20%) cases. Furthermore, we evaluated sequence conservation among haplotypes classified into genotypes A, C, D, and E by computing the information content. The preS haplotypes exhibited limited shared conserved regions, whereas the 5X haplotypes revealed two groups of conserved regions across the genotypes assessed. In conclusion, we developed an NGS-based HBV genotyping method utilizing the DB Rule for genotype classification. We identified two regions conserved across different genotypes at 5X, offering promising targets for RNA interference-based antiviral therapies.

Transcriptomic Profiles of Long Noncoding RNAs and Their Target Protein-Coding Genes Reveals Speciation Adaptation on the Qinghai-Xizang (Tibet) Plateau in Orinus.

Long noncoding RNAs (lncRNAs) are RNA molecules longer than 200 nt, which lack the ability to encode proteins and are involved in multifarious growth, development, and regulatory processes in plants and mammals. However, the environmental-regulated expression profiles of lncRNAs in Orinus that may associated with their adaptation on the Qinghai-Xizang (Tibet) Plateau (QTP) have never been characterized. Here, we utilized transcriptomic sequencing data of two Orinus species (O. thoroldii and O. kokonoricus) to identify 1624 lncRNAs, including 1119 intergenic lncRNAs, 200 antisense lncRNAs, five intronic lncRNAs, and 300 sense lncRNAs. In addition, the evolutionary relationships of Orinus lncRNAs showed limited sequence conservation among 39 species, which implied that Orinus-specific lncRNAs contribute to speciation adaptation evolution. Furthermore, considering the cis-regulation mechanism, from 286 differentially expressed lncRNAs (DElncRNAs) and their nearby protein coding genes (PCGs) between O. thoroldii and O. kokonoricus, 128 lncRNA-PCG pairs were obtained in O. thoroldii, whereas 92 lncRNA-PCG pairs were obtained in O. kokonoricus. In addition, a total of 19 lncRNA-PCG pairs in O. thoroldii and 14 lncRNA-PCG pairs in O. kokonoricus were found to participate in different biological processes, indicating that the different expression profiles of DElncRNAs between O. thoroldii and O. kokonoricus were associated with their adaptation at different elevations on the QTP. We also found several pairs of DElncRNA nearby transcription factors (TFs), indicating that these DElncRNAs regulate the expression of TFs to aid O. thoroldii in adapting to the environment. Therefore, this work systematically identified a series of lncRNAs in Orinus, laying the groundwork for further exploration into the biological function of Orinus in environmental adaptation.

ZEPPI: Proteome-scale sequence-based evaluation of protein–protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence coevolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI can be implemented on a proteome-wide scale and is applied here to millions of structural models of dimeric complexes in the Escherichia coli and human interactomes found in the PrePPI database. PrePPI's scoring function is based primarily on the evaluation of protein-protein interfaces, and ZEPPI adds a new feature to this analysis through the incorporation of evolutionary information. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. As we discuss, the standard CAPRI scores used to evaluate docking models are based on model quality and not on the ability to give yes/no answers as to whether two proteins interact. ZEPPI is able to detect weak signals from PPI models that the CAPRI scores define as incorrect and, similarly, to identify potential PPIs defined as low confidence by the current PrePPI scoring function. A number of examples that illustrate how the combination of PrePPI and ZEPPI can yield functional hypotheses are provided.

A set of nematode rRNA cistron databases and a primer assessment tool to enable more flexible and comprehensive metabarcoding.

The ITS-2-rRNA has been particularly useful for nematode metabarcoding but does not resolve all phylogenetic relationships, and reference sequences are not available for many nematode species. This is a particular issue when metabarcoding complex communities such as wildlife parasites or terrestrial and aquatic free-living nematode communities. We have used markerDB to produce four databases of distinct regions of the rRNA cistron: the 18S rRNA gene, the 28S rRNA gene, the ITS-1 intergenic spacer and the region spanning ITS-1_5.8S_ITS-2. These databases comprise 2645, 254, 13,461 and 10,107 unique full-length sequences representing 1391, 204, 1837 and 1322 nematode species, respectively. The comparative analysis illustrates the complementary value but also reveals a better representation of Clade III, IV and V than Clade I and Clade II nematodes in each case. Although the ITS-1 database includes the largest number of unique full-length sequences, the 18S rRNA database provides the widest taxonomic coverage. We also developed PrimerTC, a tool to assess primer sequence conservation across any reference sequence database, and have applied it to evaluate a large number of previously published rRNA cistron primers. We identified sets of primers that currently provide the broadest taxonomic coverage for each rRNA marker across the nematode phylum. These new resources will facilitate more comprehensive metabarcoding of nematode communities using either short-read or long-read sequencing platforms. Further, PrimerTC is available as a simple WebApp to guide or assess PCR primer design for any genetic marker and/or taxonomic group beyond the nematode phylum.

Rev Protein Diversity in HIV-1 Group M Clades.

The HIV-1 Rev protein expressed in the early stage of virus replication is involved in the nuclear export of some forms of virus RNA. Naturally occurring polymorphisms in the Rev protein could influence its activity. The association between the genetic features of different virus variants and HIV infection pathogenesis has been discussed for many years. In this study, Rev diversity among HIV-1 group M clades was analyzed to note the signatures that could influence Rev activity and, subsequently, clinical characteristics. From the Los Alamos HIV Sequence Database, 4962 Rev sequences were downloaded and 26 clades in HIV-1 group M were analyzed for amino acid changes, conservation in consensus sequences, and the presence of clade-specific amino acid substitutions (CSSs) and the Wu-Kabat protein variability coefficient (WK). Subtypes G, CRF 02_AG, B, and A1 showed the largest amino acid changes and diversity. The mean conservation of the Rev protein was 80.8%. In consensus sequences, signatures that could influence Rev activity were detected. In 15 out of 26 consensus sequences, an insertion associated with the reduced export activity of the Rev protein, 95QSQGTET96, was identified. A total of 32 CSSs were found in 16 clades, wherein A6 had the 41Q substitution in the functionally significant region of Rev. The high values of WK coefficient in sites 51 and 82, located on the Rev interaction surface, indicate the susceptibility of these positions to evolutionary replacements. Thus, the noted signatures require further investigation.

Research progress of circRNA as a biomarker of osteoporosis.

Osteoporosis, as a chronic metabolic bone disease, has the characteristic of insidious disease progression, which often leads to relatively delayed disease diagnosis. Therefore, early screening for osteoporosis has become a major public health challenge. The latest research indicates that circRNA is widely involved in the regulation of bone metabolism and is closely related to the occurrence and development of osteoporosis. Based on its high degree of sequence conservation and stability, circRNA has the potential to become a new clinical biomarker. The study of biomarkers is generally based on body fluid samples or adjacent tissue samples, with blood being the most commonly used, which can be divided into sources such as serum, plasma, peripheral blood monocytes, and plasma exosomes. Therefore, this article aims to review the research status of circRNA as a biomarker of osteoporosis.

Int&in: A machine learning-based web server for active split site identification in inteins.

Inteins are proteins that excise themselves out of host proteins and ligate the flanking polypeptides in an auto-catalytic process called protein splicing. In nature, inteins are either contiguous or split. In the case of split inteins, the two fragments must first form a complex for the splicing to occur. Contiguous inteins have previously been artificially split in two fragments because split inteins allow for distinct applications than contiguous ones. Even naturally split inteins have been split at unnatural split sites to obtain fragments with reduced affinity for one another, which are useful to create conditional inteins or to study protein-protein interactions. So far, split sites in inteins have been heuristically identified. We developed Int&in, a web server freely available for academic research (https://intein.biologie.uni-freiburg.de) that runs a machine learning model using logistic regression to predict active and inactive split sites in inteins with high accuracy. The model was trained on a dataset of 126 split sites generated using the gp41-1, Npu DnaE and CL inteins and validated using 97 split sites extracted from the literature. Despite the limited data size, the model, which uses various protein structural features, as well as sequence conservation information, achieves an accuracy of 0.79 and 0.78 for the training and testing sets, respectively. We envision Int&in will facilitate the engineering of novel split inteins for applications in synthetic and cell biology.

Staphylococcal peroxidase inhibitor (SPIN): Residue-level investigation of the helical bundle domain

Myeloperoxidase is a critical component of the antibacterial arsenal of neutrophils, whereby it consumes H2O2 as an oxidant to convert halogen and pseudohalogen anions into cytotoxic hypohalous acids. Following phagocytosis by neutrophils, the human pathogen Staphylococcus aureus secretes a potent myeloperoxidase inhibitory protein, called SPIN, as part of its immune evasion repertoire. The matured S. aureus SPIN polypeptide consists of only 73 residues yet contains two functional domains: whereas the 60 residue C-terminal helical bundle domain is responsible for MPO binding, the 13 residue N-terminal domain is required to inhibit MPO. Previous studies have informed understanding of the SPIN N-terminal domain, but comparatively little is known about the helical domain insofar as the contribution of individual residues is concerned. To address this limitation, we carried out a residue-level structure/function investigation on the helical bundle domain of S. aureus SPIN. Using sequence conservation and existing structures of SPIN bound to human MPO as a guide, we selected residues L49, E50, H51, E52, Y55, and Y75 for interrogation by site-directed mutagenesis. We found that loss of L49 or E52 reduced SPIN activity by roughly an order of magnitude, but that loss of Y55 or H51 caused progressively greater loss of inhibitory potency. Direct binding studies by SPR showed that loss of inhibitory potency in these SPIN mutants resulted from a diminished initial interaction between the inhibitor and MPO. Together, our studies provide new insights into the structure/function relationships of SPIN and identify positions Y55 and H51 as critical determinants of SPIN function.

Gene Function is a Driver of Activin Signaling Pathway Evolution Following Whole-Genome Duplication in Rainbow Trout (Oncorhynchus mykiss).

The genomes of plant and animal species are influenced by ancestral whole-genome duplication (WGD) events, which have profound impacts on the regulation and function of gene networks. To gain insight into the consequences of WGD events, we characterized the sequence conservation and expression patterns of ohnologs in the highly duplicated activin receptor signaling pathway in rainbow trout (RBT). The RBT activin receptor signaling pathway is defined by tissue-specific expression of inhibitors and ligands and broad expression of receptors and Co-Smad signaling molecules. Signaling pathway ligands exhibited shared expression, while inhibitors and Smad signaling molecules primarily express a single dominant ohnolog. Our findings suggest that gene function influences ohnolog evolution following duplication of the activin signaling pathway in RBT.

A conserved interaction between the GRIP and finger domains in ENaC subunit is critical for proteolytic activation

The epithelial sodium channel (ENaC) mediates Na+ transport in several epithelia, including the lung and the aldosterone-sensitive distal nephron. ENaC is a heterotrimer comprising ⍺βγ or δβγ subunits. Though sharing similar structural folds as other ENaC/Degenerin family members, ENaC subunits have a unique extracellular GRIP (Gating Release of Inhibition by Proteolysis) domain, and the a and g subunits have polybasic sites that can be cleaved by proteases. Aldosterone promotes double cleavage of the α and γ subunits at these sites, leading to release of the first two b-strands of the GRIP domain (P1 and P2) and channel activation. Structural and sequence conservation suggest that the interaction between P1 and the channel are mediated by a Pro-Leu motif in P1 and a conserved Trp residue in helix α2 of the finger domain. We hypothesized that this interaction is critical for the inhibitory properties of the imbedded tract, defined by P1, and thus channel activation upon release. To test our hypothesis, we substituted the conserved Trp with Ala or Phe in each of the human ENaC subunits and performed two-electrode voltage clamp experiments in Xenopus oocytes to measure the effects on function. Compared to wild type channels, we found that ⍺W251F, βW218A, βW218F, γW229F, and δW396F had no effect on amiloride-sensitive currents. Interestingly, δW396A decreased amiloride-sensitive current by 20% (p=0.01). However, amiloride-sensitive currents were greatly enhanced for Ala mutations in subunits subject to cleavage: ⍺W251A (p=0.02) and γW229A (p=0.01). In line with channel activation, ⍺W251A weakened sensitivity to extracellular Na+ (Na+ self-inhibition), and γW229A completely lost its sensitivity to both extracellular Na+ and chymotrypsin. We also tested the effect of inhibitory peptides (α8 and γ11) corresponding to the imbedded inhibitory tracts (i.e., P1 of the GRIP domain) of the ⍺ and γ subunits. ⍺W251A was insensitive to 10 μM ⍺8 compared to wild type ENaC (85% inhibition), and γW229A was insensitive to 100 μM γ11 compared to wild type ENaC (40% inhibition). Collectively, our results provide evidence for the idea that the interaction between the Pro-Leu motif in GRIP domain’s P1 β-strand and the conserved Trp in the finger domain’s α2 helix is critical for the autoinhibitory properties of P1, and therefore activation upon P1 release after proteolysis of the ⍺ and γ subunits. R01 DK125439 Kashlan (PI) ENaC Regulation by Biliary Factors. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.