ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids

Deciphering the genetic impact of signal peptide missense CTLA-4 polymorphism with rheumatoid arthritis in the Indian population: A case-control and in silico studies

Cytoplasmic T Lymphocyte Antigen-4 (CTLA-4) gene encodes for a glycoprotein, expressed on activated T-cells to transfer an inhibitory signal to control T-cell activation and proliferation. Techniques coupled with Real-time Polymerase Chain Reaction (PCR) and High-Resolution Melting Analysis (HRMA) were used to screen a missense signal peptide polymorphism (CTLA-4 + 49 A/G rs231775) in the Indian population to detect its association with Rheumatoid Arthritis (RA). Further, the resulting outcome was confirmed by Sanger's sequencing technique, and genotype frequencies were calculated. In eukaryotic cells, the M domain of the Signal Recognition Particle (SRP-54) recognizes the N-terminal region of the Signal Peptide (SP) sequence. SP directs the polypeptide chain into the Sec-61 translocon of the Endoplasmic Reticulum (ER) for further protein modification. As the Single Nucleotide Polymorphism (SNP) rs231775 lies in the signal peptide region of CTLA-4, an in-silico study was also performed to predict the mRNA stability and SP-SRP protein interaction. From the study, it was observed that the genotype frequency of rs231775 SNP G/G homozygous dominant was significantly higher in RA patients than G/A heterozygous dominant and A/A homozygous recessive conditions (Odd Ratio (OR) = 2.0862; 95 % Confidence Interval (C.I) = 1.2584 to 3.4584; Relative Risk (RR) = 1.8507; p = 0.0044). Moreover, the rs231775 SNP G allele frequency was higher in RA than the control group G = 0.407 (40.7 %) vs 0.32 (32 %). In silico approaches of Protein-Protein docking and Molecular Dynamics (MD) simulation reveal CTLA-4 rs231775 SNP (G allele) has destabilized the SP-SRP protein complex, which may affect the translocation of CTLA-4 nascent polypeptide chains into the ER via activating Regulation of Aberrant Protein Production (RAPP) pathway.

Computer-aided engineering of stabilized fibroblast growth factor 21

FGF21 is an endocrine signaling protein belonging to the family of fibroblast growth factors (FGFs). It has emerged as a molecule of interest for treating various metabolic diseases due to its role in regulating glucogenesis and ketogenesis in the liver. However, FGF21 is prone to heat, proteolytic, and acid-mediated degradation, and its low molecular weight makes it susceptible to kidney clearance, significantly reducing its therapeutic potential. Protein engineering studies addressing these challenges have generally shown that increasing the thermostability of FGF21 led to improved pharmacokinetics. Here, we describe the computer-aided design and experimental characterization of FGF21 variants with enhanced melting temperature up to 15 °C, uncompromised efficacy at activation of MAPK/ERK signaling in Hep G2 cell culture, and ability to stimulate proliferation of Hep G2 and NIH 3T3 fibroblasts cells comparable with FGF21-WT. We propose that stabilizing the FGF21 molecule by rational design should be combined with other reported stabilization strategies to maximize the pharmaceutical potential of FGF21.

The biosynthesis of the odorant 2-methylisoborneol is compartmentalized inside a protein shell

Terpenoids are the largest class of natural products, found across all domains of life. One of the most abundant bacterial terpenoids is the volatile odorant 2-methylisoborneol (2-MIB), partially responsible for the earthy smell of soil and musty taste of contaminated water. Many bacterial 2-MIB biosynthetic gene clusters were thought to encode a conserved transcription factor, named EshA in the model soil bacterium Streptomyces griseus. Here, we revise the function of EshA, now referred to as Sg Enc, and show that it is a Family 2B encapsulin shell protein. Using cryo-electron microscopy, we find that Sg Enc forms an icosahedral protein shell and encapsulates 2-methylisoborneol synthase (2-MIBS) as a cargo protein. Sg Enc contains a cyclic adenosine monophosphate (cAMP) binding domain (CBD)-fold insertion and a unique metal-binding domain, both displayed on the shell exterior. We show that Sg Enc CBDs do not bind cAMP. We find that 2-MIBS cargo loading is mediated by an N-terminal disordered cargo-loading domain and that 2-MIBS activity and Sg Enc shell structure are not modulated by cAMP. Our work redefines the function of EshA and establishes Family 2B encapsulins as cargo-loaded protein nanocompartments involved in secondary metabolite biosynthesis.

SID-2 is a conserved extracellular vesicle protein that is not associated with environmental RNAi in parasitic nematodes

In the free-living nematode Caenorhabditis elegans, the transmembrane protein SID-2 imports double-stranded RNA into intestinal cells to trigger systemic RNA interference (RNAi), allowing organisms to sense and respond to environmental cues such as the presence of pathogens. This process, known as environmental RNAi, has not been observed in the most closely related parasites that are also within clade V. Previous sequence-based searches failed to identify sid-2 orthologues in available clade V parasite genomes. In this study, we identified sid-2 orthologues in these parasites using genome synteny and protein structure-based comparison, following identification of a SID-2 orthologue in extracellular vesicles from the murine intestinal parasitic nematode Heligmosomoides bakeri. Expression of GFP-tagged H. bakeri SID-2 in C. elegans showed similar localization to the intestinal apical membrane as seen for GFP-tagged C. elegans SID-2, and further showed mobility in intestinal cells in vesicle-like structures. We tested the capacity of H. bakeri SID-2 to functionally complement environmental RNAi in a C. elegans SID-2 null mutant and show that H. bakeri SID-2 does not rescue the phenotype in this context. Our work identifies SID-2 as a highly abundant EV protein whose ancestral function may be unrelated to environmental RNAi, and rather highlights an association with extracellular vesicles in nematodes.

Optogenetically engineered Septin-7 enhances immune cell infiltration of tumor spheroids

Chimeric antigen receptor T cell therapies have achieved great success in eradicating some liquid tumors, whereas the preclinical results in treating solid tumors have proven less decisive. One of the principal challenges in solid tumor treatment is the physical barrier composed of a dense extracellular matrix, which prevents immune cells from penetrating the tissue to attack intratumoral cancer cells. Here, we improve immune cell infiltration into solid tumors by manipulating septin-7 functions in cells. Using protein allosteric design, we reprogram the three-dimensional structure of septin-7 and insert a blue light-responsive light-oxygen-voltage-sensing domain 2 (LOV2), creating a light-controllable septin-7-LOV2 hybrid protein. Blue light inhibits septin-7 function in live cells, inducing extended cell protrusions and cell polarization, enhancing cell transmigration efficiency through confining spaces. We genetically edited human natural killer cell line (NK92) and mouse primary CD8+ T-cells expressing the engineered protein, and we demonstrated improved penetration and cytotoxicity against various tumor spheroid models. Our proposed strategy to enhance immune cell infiltration is compatible with other methodologies and therefore, could be used in combination to further improve cell-based immunotherapies against solid tumors.

First description of novel compound heterozygous mutations in HYCC1: clinical evaluations and molecular analysis in patient with hypomyelinating leukodystrophy-5 with retrospective view

Hypomyelinating leukodystrophy-5 (HLD5) is a rare autosomal recessive hypomyelination disorder characterized by congenital cataract, progressive neurologic impairment, and myelin deficiency in the central and peripheral nervous system, caused by mutations in the HYCC1 gene. Here we report a 23-year-old girl with HLD5 from unrelated families. Molecular analysis was performed using sequence screening of the HYCC1 gene. In addition, in silico prediction tools and molecular investigation were used to predict the structural effect of the mutations. Results showed a novel compound heterozygous mutation in the HYCC1 gene. Moreover, in silico tools and 3D structural modeling revealed that c.521C > A (p.Ala174Glu) and c.652C > G (p.Gln218Glu) mutations could affect the structure, stability, and conformational analyses in the N-ter domain of the Hyccin protein. We also, we compared the phenotype of our patient with those of previously reported cases with HLD5 syndrome and our findings indicate the absence of reliable genotype-phenotype correlations. To the best of our knowledge, this is the first report describing a Tunisian HLD5 patient with compound heterozygous mutations (c.521C > A (p.Ala174Glu) and c.652C > G (p.Gln218Glu)) in HYCC1 gene.

Cryo-EM Structure of raiA ncRNA From Clostridium Reveals a New RNA 3D Fold

Advancements in genome-wide sequence analysis have led to the discovery of numerous novel bacterial non-coding RNAs (ncRNAs). These ncRNAs have been categorized into various RNA families and classes based on their size, structure, function, and evolutionary relationships. One such ncRNA family, raiA, is notably abundant in the bacterial phyla Firmicutes and Actinobacteria and is remarkably well-conserved across many Gram-positive bacteria. In this study, we integrated cryo-electron microscopy single-particle analysis with computational modeling and biochemical techniques to elucidate the structural characteristics of raiA from Clostridium sp. CAG 138. Our findings reveal the globular 3D fold of raiA, providing valuable structural insights. This analysis paves the way for future investigations into the functional properties of raiA, potentially uncovering new regulatory mechanisms in bacterial ncRNAs.

Microbial Cell Factories in the Bioeconomy Era: From Discovery to Creation

Microbial cell factories (MCFs) are extensively used to produce a wide array of bioproducts, such as bioenergy, biochemical, food, nutrients, and pharmaceuticals, and have been regarded as the "chips" of biomanufacturing that will fuel the emerging bioeconomy era. Biotechnology advances have led to the screening, investigation, and engineering of an increasing number of microorganisms as diverse MCFs, which are the workhorses of biomanufacturing and help develop the bioeconomy. This review briefly summarizes the progress and strategies in the development of robust and efficient MCFs for sustainable and economic biomanufacturing. First, a comprehensive understanding of microbial chassis cells, including accurate genome sequences and corresponding annotations; metabolic and regulatory networks governing substances, energy, physiology, and information; and their similarity and uniqueness compared with those of other microorganisms, is needed. Moreover, the development and application of effective and efficient tools is crucial for engineering both model and nonmodel microbial chassis cells into efficient MCFs, including the identification and characterization of biological parts, as well as the design, synthesis, assembly, editing, and regulation of genes, circuits, and pathways. This review also highlights the necessity of integrating automation and artificial intelligence (AI) with biotechnology to facilitate the development of future customized artificial synthetic MCFs to expedite the industrialization process of biomanufacturing and the bioeconomy.

Interleukin-2 receptor α (IL-2Rα/CD25) shedding is differentially regulated by N- and O-glycosylation

The cytokine interleukin-2 (IL-2) is a critical regulator of immune responses, with an especially well-characterized role in regulating T-cell homeostasis. IL-2 signaling involves three distinct receptor subunits: the IL-2Rα (CD25), IL-2Rβ, and IL-2Rγ. The intracellular transduction of IL-2-induced signals is strictly dependent on IL-2Rβ and IL-2Rγ, while the IL-2Rα is not directly involved in signaling. Instead, it has the highest affinity towards IL-2 and is thus responsible for regulating the affinity of a cell for IL-2. In addition to the membrane-bound IL-2Rα, a soluble form of the receptor (sIL-2Rα) has been described, which is present in the blood of healthy individuals, increased under various pathological conditions, and able to bind IL-2 and thus modulate its function. The sIL-2Rα is generated by proteolytic cleavage of the membrane-bound receptor. Here, we analyze whether glycosylation of the IL-2Rα regulates its proteolysis. We find that constitutive IL-2Rα shedding is affected by glycosylation and discover distinct roles for N- and O-glycosylation. Furthermore, we show that induced shedding by the metalloproteases ADAM10 and ADAM17 is also differentially regulated by distinct types of glycans. Finally, we identify a specific role for an N-glycan at an exosite in ADAM17-mediated proteolysis that does not affect ADAM10, indicating distinct substrate recognition mechanisms. These results further the understanding of the mechanisms leading to sIL-2Rα generation, and thus offer the opportunity to specifically modulate the generation of the soluble receptor.

Structural basis for FN3K-mediated protein deglycation

Protein glycation is a universal, non-enzymatic modification that occurs when a sugar covalently attaches to a primary amine. These spontaneous modifications may have deleterious or regulatory effects on protein function, and their removal is mediated by the conserved metabolic kinase fructosamine-3-kinase (FN3K). Despite its crucial role in protein repair, we currently have a poor understanding of how FN3K engages or phosphorylates its substrates. By integrating structural biology and biochemistry, we elucidated the catalytic mechanism for FN3K-mediated protein deglycation. Our work identifies key amino acids required for binding and phosphorylating glycated substrates and reveals the molecular basis of an evolutionarily conserved protein repair pathway. Additional structural-functional studies revealed unique structural features of human FN3K as well as differences in the dimerization behavior and regulation of FN3K family members. Our findings improve our understanding of the structure of FN3K and its catalytic mechanism, which opens new avenues for therapeutically targeting FN3K.

Preclinical alternative drug discovery programs for monogenic rare diseases. Should small molecules or gene therapy be used? The case of hereditary spastic paraplegias

Patients diagnosed with rare diseases and their and families search desperately to organize drug discovery campaigns. Alternative models that differ from default paradigms offer real opportunities. There are, however, no clear guidelines for the development of such models, which reduces success rates and raises costs. We address the main challenges in making the discovery of new preclinical treatments more accessible, using rare hereditary paraplegia as a paradigmatic case. First, we discuss the necessary expertise, and the patients' clinical and genetic data. Then, we revisit gene therapy, de novo drug development, and drug repurposing, discussing their applicability. Moreover, we explore a pool of recommended in silico tools for pathogenic variant and protein structure prediction, virtual screening, and experimental validation methods, discussing their strengths and weaknesses. Finally, we focus on successful case applications.

Molecular Dynamic Simulations Reveal that Water-Soluble QTY-Variants of Glutamate Transporters EAA1, EAA2 and EAA3 Retain the Conformational Characteristics of Native Transporters

OBJECTIVE

Glutamate transporters play a crucial role in neurotransmitter homeostasis, but studying their structure and function is challenging due to their membrane-bound nature. This study aims to investigate whether water-soluble QTY-variants of glutamate transporters EAA1, EAA2 and EAA3 retain the conformational characteristics and dynamics of native membrane-bound transporters.

METHODS

Molecular dynamics simulations and comparative genomics were used to analyze the structural dynamics of both native transporters and their QTY-variants. Native transporters were simulated in lipid bilayers, while QTY-variants were simulated in aqueous solution. Lipid distortions, relative solvent accessibilities, and conformational changes were examined. Evolutionary conservation profiles were correlated with structural dynamics. Statistical analyses included multivariate analysis to account for confounding variables.

RESULTS

QTY-variants exhibited similar residue-wise conformational dynamics to their native counterparts, with correlation coefficients of 0.73 and 0.56 for EAA1 and EAA3, respectively (p < 0.001). Hydrophobic interactions of native helices correlated with water interactions of QTY- helices (rs = 0.4753, p < 0.001 for EAA1). QTY-variants underwent conformational changes resembling the outward-to-inward transition of native transporters.

CONCLUSIONS

Water-soluble QTY-variants retain key structural properties of native glutamate transporters and mimic aspects of native lipid interactions, including conformational flexibility. This research provides valuable insights into the conformational changes and molecular mechanisms of glutamate transport, potentially offering a new approach for studying membrane protein dynamics and drug interactions.

A side-by-side comparison of variant function measurements using deep mutational scanning and base editing

Variant annotation is a crucial objective in mammalian functional genomics. Deep Mutational Scanning (DMS) is a well-established method for annotating human gene variants, but CRISPR base editing (BE) is emerging as an alternative. However, questions remain about how well high-throughput base editing measurements can annotate variant function and the extent of downstream experimental validation required. This study presents the first direct comparison of DMS and BE in the same lab and cell line. Results indicate that focusing on the most likely edits and highest efficiency sgRNAs enhances the agreement between a "gold standard" DMS dataset and a BE screen. A simple filter for sgRNAs making single edits in their window could sufficiently annotate a large proportion of variants directly from sgRNA sequencing of large pools. When multi-edit guides are unavoidable, directly measuring the variants created in the pool, rather than sgRNA abundance, can recover high-quality variant annotation measurements in multiplexed pools. Taken together, our data show a surprising degree of correlation between base editor data and gold standard deep mutational scanning.

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

BACKGROUND

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

Molecular insights into the initiation step of the Rcs signaling pathway

The Rcs pathway is repressed by the inner membrane protein IgaA under non-stressed conditions. This repression is hypothesized to be relieved by the binding of the outer membrane-anchored RcsF to IgaA. However, the precise mechanism by which RcsF binding triggers the signaling remains unclear. Here, we present the 1.8 Å resolution crystal structure capturing the interaction between IgaA and RcsF. Our comparative structural analysis, examining both the bound and unbound states of the periplasmic domain of IgaA (IgaAp), highlights rotational flexibility within IgaAp. Conversely, the conformation of RcsF remains unchanged upon binding. Our in vivo and in vitro studies do not support the model of a stable complex involving RcsF, IgaAp, and RcsDp. Instead, we demonstrate that the elements beyond IgaAp play a role in the interaction between IgaA and RcsD. These findings collectively allow us to propose a potential mechanism for the signaling across the inner membrane through IgaA.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

Drug to genome to drug: a computational large-scale chemogenomics screening for novel drug candidates against sporotrichosis

Sporotrichosis is recognized as the predominant subcutaneous mycosis in South America, attributed to pathogenic species within the Sporothrix genus. Notably, in Brazil, Sporothrix brasiliensis emerges as the principal species, exhibiting significant sapronotic, zoonotic and enzootic epidemic potential. Consequently, the discovery of novel therapeutic agents for the treatment of sporotrichosis is imperative. The present study is dedicated to the repositioning of pharmaceuticals for sporotrichosis therapy. To achieve this goal, we designed a pipeline with the following steps: (a) compilation and preparation of Sporothrix genome data; (b) identification of orthologous proteins among the species; (c) identification of homologous proteins in publicly available drug-target databases; (d) selection of Sporothrix essential targets using validated genes from Saccharomyces cerevisiae; (e) molecular modeling studies; and (f) experimental validation of selected candidates. Based on this approach, we were able to prioritize eight drugs for in vitro experimental validation. Among the evaluated compounds, everolimus and bifonazole demonstrated minimum inhibitory concentration (MIC) values of 0.5 µg/mL and 4.0 µg/mL, respectively. Subsequently, molecular docking studies suggest that bifonazole and everolimus may target specific proteins within S. brasiliensis- namely, sterol 14-α-demethylase and serine/threonine-protein kinase TOR, respectively. These findings shed light on the potential binding affinities and binding modes of bifonazole and everolimus with their probable targets, providing a preliminary understanding of the antifungal mechanism of action of these compounds. In conclusion, our research advances the understanding of the therapeutic potential of bifonazole and everolimus, supporting their further investigation as antifungal agents for sporotrichosis in prospective hit-to-lead and preclinical investigations.

Bioinformatics Study on Site-Specific Variations of Eotaxin-3, a Key Chemokine in Eosinophilic Esophagitis (EoE)

Eotaxin-3 is a key chemokine with a relevant role in eosinophilic esophagitis, a rare chronic immune/antigen-mediated inflammatory disorder. Eotaxin-3 is a potent activator of eosinophil emergence and migration, which may lead to allergic airway inflammation. We investigated, using bioinformatics tools, the protein structure and the possible effects of the known variations reported in public databases. Following a procedure already established, we created a 3D model of the whole protein and modeled the structure of 105 protein variants due to known point mutations. The effects of the amino acid substitution at the level of impact on protein structure, stability, and possibly function were detected by the bioinformatics procedure and described in detail. A web application was implemented to browse the results of the analysis and visualize the 3D models, with the opportunity of downloading the models and analyzing them using their own software. Among 105 amino acid substitutions investigated, the study evidenced in 44 cases at least one change in any of the investigated structural parameters. Other six variations are also relevant, although a structural effect was not detected by our analysis, because they affected amino acids highly conserved, which suggests a possible function role. All these variations should be the object of particular attention, as they may induce a loss of functionality in the protein.

In Silico Functional and Structural Analysis of Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in Human Paired Box 4 Gene

In human genome, members of Paired box (PAX) transcription factor family are highly sequence-specific DNA-binding proteins. Among PAX gene family members, PAX4 gene has significant role in growth, proliferation, differentiation, and insulin secretion of pancreatic β-cells. Single nucleotide polymorphisms (SNPs) in PAX4 gene progress in the pathogenesis of various human diseases. Hence, the molecular mechanism of how these SNPs in PAX4 gene significantly progress diseases pathogenesis needs to be elucidated. For the reason, a series of bioinformatic analyzes were done to identify the SNPs of PAX4 gene that contribute in diseases pathogenesis. From the analyzes, 4145 SNPs (rsIDs) in PAX4 gene were obtained, where, 362 missense (8.73%), 169 synonymous (4.08%), and 2323 intron variants (56.04%). The rest SNPs were unspecified. Among the 362 missense variants, 118 nsSNPs were found as deleterious in SIFT analysis. Among those, 25 nsSNPs were most probably damaging and 23 were deleterious as observed in PolyPhen-2 and PROVEAN analyzes, respectively. Following all analyzes, 14 nsSNPs (rs149708455, rs115887120, rs147279315, rs35155575, rs370095957, rs373939873, rs145468905, rs121917718, rs2233580, rs3824004, rs372751660, rs369459316, rs375472849, rs372497946) were common and observed as deleterious, probably damaging, affective and diseases associated. Following structural analyzes, 11 nsSNPs guided proteins were found as most unstable and highly conserved. Among these, R20W, R39Q, R45Q, R60H, G65D, and A223D mutated proteins were highly harmful. Hence, the results from above-mentioned integrated comprehensive bioinformatic analyzes guide how different nsSNPs in PAX4 gene alter structural and functional characteristics of the protein that might progress diseases pathogenesis in human including type 2 diabetes.

Day-night gene expression reveals circadian gene disco as a candidate for diel-niche evolution in moths

Temporal ecological niche partitioning is an underappreciated driver of speciation. While insects have long been models for circadian biology, the genes and circuits that allow adaptive changes in diel-niches remain poorly understood. We compared gene expression in closely related day- and night-active non-model wild silk moths, with otherwise similar ecologies. Using an ortholog-based pipeline to compare RNA-Seq patterns across two moth species, we find over 25 pairs of gene orthologs showing differential expression. Notably, the gene disco, involved in circadian control, optic lobe and clock neuron development in Drosophila, shows robust adult circadian mRNA cycling in moth heads. Disco is highly conserved in moths and has additional zinc-finger domains with specific nocturnal and diurnal mutations. We propose disco as a candidate gene for the diversification of temporal diel-niche in moths.

A bacterial sense of touch: T4P retraction motor as a means of surface sensing by Pseudomonas aeruginosa PA14

Most microbial cells found in nature exist in matrix-covered, surface-attached communities known as biofilms. This mode of growth is initiated by the ability of the microbe to sense a surface on which to grow. The opportunistic pathogen Pseudomonas aeruginosa (Pa) PA14 utilizes a single polar flagellum and type 4 pili (T4P) to sense surfaces. For Pa, T4P-dependent "twitching" motility is characterized by effectively pulling the cell across a surface through a complex process of cooperative binding, pulling, and unbinding. T4P retraction is powered by hexameric ATPases. Pa cells that have engaged a surface increase production of the second messenger cyclic AMP (cAMP) over multiple generations via the Pil-Chp system. This rise in cAMP allows cells and their progeny to become better adapted for surface attachment and activates virulence pathways through the cAMP-binding transcription factor Vfr. While many studies have focused on mechanisms of T4P twitching and regulation of T4P production and function by the Pil-Chp system, the mechanism by which Pa senses and relays a surface-engagement signal to the cell is still an open question. Here we review the current state of the surface sensing literature for Pa, with a focus on T4P, and propose an integrated model of surface sensing whereby the retraction motor PilT senses and relays the signal to the Pil-Chp system via PilJ to drive cAMP production and adaptation to a surface lifestyle.

Metagenome-derived SusD-homologs affiliated with Bacteroidota bind to synthetic polymers

Starch utilization system (Sus)D-homologs are well known for their carbohydrate-binding capabilities and are part of the sus operon in microorganisms affiliated with the phylum Bacteroidota. Until now, SusD-like proteins have been characterized regarding their affinity toward natural polymers. In this study, three metagenomic SusD homologs (designated SusD1, SusD38489, and SusD70111) were identified and tested with respect to binding to natural and non-natural polymers. SusD1 and SusD38489 are cellulose-binding modules, while SusD70111 preferentially binds chitin. Employing translational fusion proteins with superfolder GFP (sfGFP), pull-down assays, and surface plasmon resonance (SPR) has provided evidence for binding to polyethylene terephthalate (PET) and other synthetic polymers. Structural analysis suggested that a Trp triad might be involved in protein adsorption. Mutation of these residues to Ala resulted in an impaired adsorption to microcrystalline cellulose (MC), but not so to PET and other synthetic polymers. We believe that the characterized SusDs, alongside the methods and considerations presented in this work, will aid further research regarding bioremediation of plastics.

IMPORTANCE

SusD1 and SusD38489 can be considered for further applications regarding their putative adsorption toward fossil-fuel based polymers. This is the first time that SusD homologs from the polysaccharide utilization loci (PUL), largely described for the phylum Bacteroidota, are characterized as synthetic polymer-binding proteins.

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

Computational exploration of SLC14A1 genetic variants through structure modeling, protein-ligand docking, and molecular dynamics simulation

The urea transporter UT-B1, encoded by the SLC14A1 gene, has been hypothesized to be a significant protein whose deficiency and dysfunction contribute to the pathogenesis of bladder cancer and many other diseases. Several studies reported the association of genetic alterations in the SLC14A1 (UT-B1) gene with bladder carcinogenesis, suggesting a need for thorough characterization of the UT-B1 protein's coding and non-coding variants. This study used various computational techniques to investigate the commonly occurring germ-line missense and non-coding SNPs (ncSNPs) of the SLC14A1 gene (UT-B1) for their structural, functional, and molecular implications for disease susceptibility and dysfunctionality. SLC14A1 missense variants, primarily identified from the ENSEMBL genome browser, were screened through twelve functionality prediction tools leading to two variants D280Y (predicted detrimental by maximum tools) and D280N (high global MAF) for rs1058396. Subsequently, the ConSurf and NetSurf tools revealed the D280 residue to be in a variable site and exposed on the protein surface. According to I-Mutant2.0 and MUpro, both variants are predicted to cause a significant effect on protein stability. Analysis of molecular docking anticipated these two variants to decrease the binding affinity of UT-B1 protein for the examined ligands to a significant extent. Molecular dynamics also disclosed the possible destabilization of the UT-B1 protein due to single nucleotide polymorphism compared to wild-type protein which may result in impaired protein function. Furthermore, several non-coding SNPs were estimated to affect transcription factor binding and regulation of SLC14A1 gene expression. Additionally, two ncSNPs were found to affect miRNA-based post-transcriptional regulation by creating new seed regions for miRNA binding. This comprehensive in-silico study of SLC14A1 gene variants may serve as a springboard for future large-scale investigations examining SLC14A1 polymorphisms.

Prediction and assessment of deleterious and disease causing nonsynonymous single nucleotide polymorphisms (nsSNPs) in human FOXP4 gene: An in - silico study

In humans, FOXP gene family is involved in embryonic development and cancer progression. The FOXP4 (Forkhead box protein P4) gene belongs to this FOXP gene family. FOXP4 gene plays a crucial role in oncogenesis. Single nucleotide polymorphisms are biological markers and common determinants of human diseases. Mutations can largely affect the function of the corresponding protein. Therefore, the molecular mechanism of nsSNPs in the FOXP4 gene needs to be elucidated. Initially, the SNPs of the FOXP4 gene were extracted from the dbSNP database and a total of 23124 SNPs was found, where 555 nonsynonymous, 20525 intronic, 1114 noncoding transcript, 334 synonymous were obtained and the rest were unspecified. Then, a series of bioinformatics tools (SIFT, PolyPhen2, SNAP2, PhD SNP, PANTHER, I-Mutant2.0, MUpro, GOR IV, ConSurf, NetSurfP 2.0, HOPE, DynaMut2, GeneMANIA, STRING and Schrodinger) were used to explore the effect of nsSNPs on FOXP4 protein function and structural stability. First, 555 nsSNPs were analyzed using SIFT, of which 57 were found as deleterious. Following, PolyPhen2, SNAP2, PhD SNP and PANTHER analyses, 10 nsSNPs (rs372762294, rs141899153, rs142575732, rs376938850, rs367607523, rs112517943, rs140387832, rs373949416, rs373949416 and rs376160648) were common and observed as deleterious, damaging and diseases associated. Following that, using I-Mutant2.0 and MUpro servers, 7 nsSNPs were found to be the most unstable. GOR IV predicted that these seven nsSNPs affect protein structure by altering the protein contents of alpha helixes, extended strands, and random coils. Following DynaMut2, 5 nsSNPs showed a decrease in the ΔΔG value compared with the wild-type and were found to be responsible for destabilizing the corresponding protein. GeneMANIA and STRING network revealed interaction of FOXP4 with other genes. Finally, molecular dynamics simulation analysis revealed consistent fluctuation in RMSD and RMSF values, Rg and hydrogen bonds in the mutant proteins compared with WT, which might alter the functional and structural stability of the corresponding protein. As a result, the aforementioned integrated comprehensive bioinformatic analyses provide insight into how various nsSNPs of the FOXP4 gene change the structural and functional properties of the corresponding protein, potentially proceeding with the pathophysiology of human diseases.

Insights into stereoselective ring formation in canonical strigolactone: Identification of a dirigent domain-containing enzyme catalyzing orobanchol synthesis

Strigolactones (SLs) are plant apocarotenoids with diverse roles and structures. Canonical SLs, widespread and characterized by structural variations in their tricyclic lactone (ABC-ring), are classified into two types based on C-ring configurations. The steric C-ring configuration emerges during the BC-ring closure, downstream of the biosynthetic intermediate, carlactonoic acid (CLA). Most plants produce either type of canonical SLs stereoselectively, e.g., tomato (Solanum lycopersicum) yields orobanchol with an α-oriented C-ring. The mechanisms driving SL structural diversification are partially understood, with limited insight into functional implications. Furthermore, the exact molecular mechanism for the stereoselective BC-ring closure reaction is yet to be known. We identified an enzyme, the stereoselective BC-ring-forming factor (SRF), from the dirigent protein (DIR) family, specifically the DIR-f subfamily, whose biochemical function had not been characterized, making it a key enzyme in stereoselective canonical SL biosynthesis with the α-oriented C-ring. We first confirm the precise catalytic function of the tomato cytochrome P450 SlCYP722C, previously shown to be involved in orobanchol biosynthesis [T. Wakabayashi et al., Sci. Adv. 5, eaax9067 (2019)], to convert CLA to 18-oxocarlactonoic acid. We then show that SRF catalyzes the stereoselective BC-ring closure reaction of 18-oxocarlactonoic acid, forming orobanchol. Our methodology combines experimental and computational techniques, including SRF structure prediction and conducting molecular dynamics simulations, suggesting a catalytic mechanism based on the conrotatory 4π-electrocyclic reaction for the stereoselective BC-ring formation in orobanchol. This study sheds light on the molecular basis of how plants produce SLs with specific stereochemistry in a controlled manner.

A programmable seekRNA guides target selection by IS1111 and IS110 type insertion sequences

IS1111 and IS110 insertion sequence (IS) family members encode an unusual DEDD transposase type and exhibit specific target site selection. The IS1111 group include identifiable subterminal inverted repeats (sTIR) not found in the IS110 type1. IS in both families include a noncoding region (NCR) of significant length and, as each individual IS or group of closely related IS selects a different site, we had previously proposed that an NCR-derived RNA was involved in target selection2. Here, we find that the NCR is usually downstream of the transposase gene in IS1111 family IS and upstream in the IS110 type. Four IS1111 and one IS110 family members that target different sequences are used to demonstrate that the NCR determines a short seeker RNA (seekRNA) that co-purified with the transposase. The seekRNA is essential for transposition of the IS or a cargo flanked by IS ends from and to the preferred target. Short sequences matching both top and bottom strands of the target are present in the seekRNA but their order in IS1111 and IS110 family IS is reversed. Reprogramming the seekRNA and donor flank to target a different site is demonstrated, indicating future biotechnological potential for these systems.

Genetic characterization, structural analysis, and detection of positive selection in small heat shock proteins of Cypriniformes and Clupeiformes

In the current investigation, a total of 42 full-length, non-redundant small heat shock proteins (sHsp) were detected in Cyprinus carpio, Labeo rohita, Danio rerio, Salmo salar, Oncorhynchus mykiss, and Clupea harengus. The sHsp genes were classified into three groups based on phylogenetic analysis. All the sHsps were shown to have higher aliphatic index values, which is an indication that these proteins are more thermally stable. The hydrophilic nature of sHsps was deduced from the fact that all fish species had negative GRAVY scores. In all of the representative fish species, sHsp genes were assigned to distinct chromosomes in an inconsistent and unequal manner. Segmental duplications are the main events that have contributed to the expansion of the sHsp genes in all species. We were also able to determine the selective pressure that was placed on particular codons and discovered several significant coding sites within the coding region of sHsps. Eventually, diversifying positive selection was found to be connected with evolutionary changes in sHsp proteins, which showed that gene evolution controlled the fish adaption event in response to environmental conditions. Clarification of the links between sHsps and environmental stress in fish will be achieved through rigorous genomic comparison, which will also yield substantial new insights.

Identification and analysis of oncogenic non-synonymous single nucleotide polymorphisms in the human NRAS gene: An exclusive in silico study

BACKGROUND

N-ras protein is encoded by the NRAS gene and operates as GDP-GTP-controlled on/off switching. N-ras interacts with cellular signaling networks that regulate various cellular activities including cell proliferation and survival. The nonsynonymous single nucleotide polymorphism (nsSNPs)-mediated alteration can substantially disrupt the structure and activity of the corresponding protein. N-ras has been reported to be associated with numerous diseases including cancers due to the nsSNPs. A comprehensive study on the NRAS gene to unveil the potentially damaging and oncogenic nsSNPs is yet to be accomplished. Hence, this extensive in silico study is intended to identify the disease-associated, specifically oncogenic nsSNPs of the NRAS gene.

RESULTS

Out of 140 missense variants, 7 nsSNPs (I55R, G60E, G60R, Y64D, L79F, D119G, and V152F) were identified to be damaging utilizing 10 computational tools that works based on different algorithms with high accuracy. Among those, G60E, G60R, and D119G variants were further filtered considering their location in the highly conserved region and later identified as oncogenic variants. Interestingly, G60E and G60R variants were revealed to be particularly associated with lung adenocarcinoma, rhabdomyosarcoma, and prostate adenocarcinoma. Therefore, D119G could be subjected to detailed investigation for identifying its association with specific cancer.

CONCLUSION

This in silico study identified the deleterious and oncogenic missense variants of the human NRAS gene that could be utilized for designing further experimental investigation. The outcomes of this study would be worthwhile in future research for developing personalized medicine.

Untangling the Evolution of the Receptor-Binding Motif of SARS-CoV-2

The spike protein determines the host-range specificity of coronaviruses. In particular, the Receptor-Binding Motif in the spike protein from SARS-CoV-2 contains the amino acids involved in molecular recognition of the host Angiotensin Converting Enzyme 2. Therefore, to understand how SARS-CoV-2 acquired its capacity to infect humans it is necessary to reconstruct the evolution of this important motif. Early during the pandemic, it was proposed that the SARS-CoV-2 Receptor-Binding Domain was acquired via recombination with a pangolin infecting coronavirus. This proposal was challenged by an alternative explanation that suggested that the Receptor-Binding Domain from SARS-CoV-2 did not originated via recombination with a coronavirus from a pangolin. Instead, this alternative hypothesis proposed that the Receptor-Binding Motif from the bat coronavirus RaTG13, was acquired via recombination with an unidentified coronavirus. And as a consequence of this event, the Receptor-Binding Domain from the pangolin coronavirus appeared as phylogenetically closer to SARS-CoV-2. Recently, the genomes from coronaviruses from Cambodia (bat_RShST182/200) and Laos (BANAL-20-52/103/247) which are closely related to SARS-CoV-2 were reported. However, no detailed analysis of the evolution of the Receptor-Binding Motif from these coronaviruses was reported. Here we revisit the evolution of the Receptor-Binding Domain and Motif in the light of the novel coronavirus genome sequences. Specifically, we wanted to test whether the above coronaviruses from Cambodia and Laos were the source of the Receptor-Binding Domain from RaTG13. We found that the Receptor-Binding Motif from these coronaviruses is phylogenetically closer to SARS-CoV-2 than to RaTG13. Therefore, the source of the Receptor-Binding Domain from RaTG13 is still unidentified. In accordance with previous studies, our results are consistent with the hypothesis that the Receptor-Binding Motif from SARS-CoV-2 evolved by vertical inheritance from a bat-infecting population of coronaviruses.

P4HA2 activates mTOR via hydroxylation and targeting P4HA2-mTOR inhibits lung adenocarcinoma cell growth

Mammalian target of rapamycin (mTOR) kinase functions as a central regulator of cell growth and metabolism, and its complexes mTORC1 and mTORC2 phosphorylate distinct substrates. Dysregulation of mTOR signaling is commonly implicated in human diseases, including cancer. Despite three decades of active research in mTOR, much remains to be determined. Here, we demonstrate that prolyl 4-hydroxylase alpha-2 (P4HA2) binds directly to mTOR and hydroxylates one highly conserved proline 2341 (P2341) within a kinase domain of mTOR, thereby activating mTOR kinase and downstream effector proteins (e.g. S6K and AKT). Moreover, the hydroxylation of P2341 strengthens mTOR stability and allows mTOR to accurately recognize its substrates such as S6K and AKT. The growth of lung adenocarcinoma cells overexpressing mTORP2341A is significantly reduced when compared with that of cells overexpressing mTORWT. Interestingly, in vivo cell growth assays show that targeting P4HA2-mTOR significantly suppresses lung adenocarcinoma cell growth. In summary, our study reveals an undiscovered hydroxylation-regulatory mechanism by which P4HA2 directly activates mTOR kinase, providing insights for therapeutically targeting mTOR kinase-driven cancers.

In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene

Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.

Three-component systems represent a common pathway for extracytoplasmic addition of pentofuranose sugars into bacterial glycans

Cell surface glycans are major drivers of antigenic diversity in bacteria. The biochemistry and molecular biology underpinning their synthesis are important in understanding host-pathogen interactions and for vaccine development with emerging chemoenzymatic and glycoengineering approaches. Structural diversity in glycostructures arises from the action of glycosyltransferases (GTs) that use an immense catalog of activated sugar donors to build the repeating unit and modifying enzymes that add further heterogeneity. Classical Leloir GTs incorporate α- or β-linked sugars by inverting or retaining mechanisms, depending on the nucleotide sugar donor. In contrast, the mechanism of known ribofuranosyltransferases is confined to β-linkages, so the existence of α-linked ribofuranose in some glycans dictates an alternative strategy. Here, we use Citrobacter youngae O1 and O2 lipopolysaccharide O antigens as prototypes to describe a widespread, versatile pathway for incorporating side-chain α-linked pentofuranoses by extracytoplasmic postpolymerization glycosylation. The pathway requires a polyprenyl phosphoribose synthase to generate a lipid-linked donor, a MATE-family flippase to transport the donor to the periplasm, and a GT-C type GT (founding the GT136 family) that performs the final glycosylation reaction. The characterized system shares similarities, but also fundamental differences, with both cell wall arabinan biosynthesis in mycobacteria, and periplasmic glucosylation of O antigens first discovered in Salmonella and Shigella. The participation of auxiliary epimerases allows the diversification of incorporated pentofuranoses. The results offer insight into a broad concept in microbial glycobiology and provide prototype systems and bioinformatic guides that facilitate discovery of further examples from diverse species, some in currently unknown glycans.

Predicting hotspots for disease-causing single nucleotide variants using sequences-based coevolution, network analysis, and machine learning

To enable personalized medicine, it is important yet highly challenging to accurately predict disease-causing mutations in target proteins at high throughput. Previous computational methods have been developed using evolutionary information in combination with various biochemical and structural features of protein residues to discriminate neutral vs. deleterious mutations. However, the power of these methods is often limited because they either assume known protein structures or treat residues independently without fully considering their interactions. To address the above limitations, we build upon recent progress in machine learning, network analysis, and protein language models, and develop a sequences-based variant site prediction workflow based on the protein residue contact networks: 1. We employ and integrate various methods of building protein residue networks using state-of-the-art coevolution analysis tools (RaptorX, DeepMetaPSICOV, and SPOT-Contact) powered by deep learning. 2. We use machine learning algorithms (Random Forest, Gradient Boosting, and Extreme Gradient Boosting) to optimally combine 20 network centrality scores to jointly predict key residues as hot spots for disease mutations. 3. Using a dataset of 107 proteins rich in disease mutations, we rigorously evaluate the network scores individually and collectively (via machine learning). This work supports a promising strategy of combining an ensemble of network scores based on different coevolution analysis methods (and optionally predictive scores from other methods) via machine learning to predict hotspot sites of disease mutations, which will inform downstream applications of disease diagnosis and targeted drug design.

Diversity of the Rysto gene conferring resistance to potato virus Y in wild relatives of potato

BACKGROUND

Potato virus Y (PVY) is among the economically most damaging viral pathogen in production of potato (Solanum tuberosum) worldwide. The gene Rysto derived from the wild potato relative Solanum stoloniferum confers extreme resistance to PVY.

RESULTS

The presence and diversity of Rysto were investigated in wild relatives of potato (298 genotypes representing 29 accessions of 26 tuber-bearing Solanum species) using PacBio amplicon sequencing. A total of 55 unique Rysto-like sequences were identified in 72 genotypes representing 12 accessions of 10 Solanum species and six resistant controls (potato cultivars Alicja, Bzura, Hinga, Nimfy, White Lady and breeding line PW363). The 55 Rysto-like sequences showed 89.87 to 99.98% nucleotide identity to the Rysto reference gene, and these encoded in total 45 unique protein sequences. While Rysto-like26 identified in Alicja, Bzura, White Lady and Rysto-like16 in PW363 encode a protein identical to the Rysto reference, the remaining 44 predicted Rysto-like proteins were 65.93 to 99.92% identical to the reference. Higher levels of diversity of the Rysto-like sequences were found in the wild relatives of potato than in the resistant control cultivars. The TIR and NB-ARC domains were the most conserved within the Rysto-like proteins, while the LRR and C-JID domains were more variable. Several Solanum species, including S. antipoviczii and S. hougasii, showed resistance to PVY. This study demonstrated Hyoscyamus niger, a Solanaceae species distantly related to Solanum, as a host of PVY.

CONCLUSIONS

The new Rysto-like variants and the identified PVY resistant potato genotypes are potential resistance sources against PVY in potato breeding. Identification of H. niger as a host for PVY is important for cultivation of this plant, studies on the PVY management, its ecology, and migrations. The amplicon sequencing based on PacBio SMRT and the following data analysis pipeline described in our work may be applied to obtain the nucleotide sequences and analyze any full-length genes from any, even polyploid, organisms.

Characterization of missense nonsynonymous single-nucleotide polymorphism of runt-related transcription factor-2 gene - An in silico approach

OBJECTIVES

Single-nucleotide polymorphism (SNP) codes for multiple amino acids, impacting protein functions and disease prognosis. Runt-related transcription factor-2 (RUNX2), a transcription factor linked to osteoblast differentiation, regulates cell proliferation in endothelium and osteoblastic cells. Understanding Runx2's role in nonosseous tissues is rapidly advancing. This study aims to identify harmful SNPs of the RUNX2 gene that may alter disease susceptibility using computational techniques.

METHODS

The study uses various in silico methods to identify nonsynonymous SNPs (nsSNPs) of the RUNX2 gene, which could potentially alter protein structure and functions, with further analyses by I-Mutant, ConSurf, Netsurf 3.0, GeneMANIA, and Have (y)Our Protein Explained.

RESULTS

Six missense nsSNPs were identified as potentially harmful, disease-causing, and damaging. Four were found to be unstable, while five were conserved. All six nsSNPs had a coiled secondary structure. Five nsSNPs were found to be destabilized.

CONCLUSION

The RUNX2 gene's deleterious missense nsSNPs were identified by this study, and they may be exploited in future experimental studies. These high-risk nsSNPs might be considered target molecules in therapeutic and diagnostic therapies in teeth and bone development.

Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR geneComputational analysis of bovine PKLR gene

Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins' stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle.Communicated by Ramaswamy H. Sarma.

Molecular and computational characterization of ABCB11 and ABCG5 variants in Tunisian patients with neonatal/infantile low-GGT intrahepatic cholestasis: Genetic diagnosis and genotype-phenotype correlation assessment

Many inherited conditions cause hepatocellular cholestasis in infancy, including progressive familial intrahepatic cholestasis (PFIC), a heterogeneous group of diseases with highly overlapping symptoms. In our study, six unrelated Tunisian infants with PFIC suspicion were the subject of a panel-target sequencing followed by an exhaustive bioinformatic and modeling investigations. Results revealed five disease-causative variants including known ones: (the p.Asp482Gly and p.Tyr354 * in the ABCB11 gene and the p.Arg446 * in the ABCC2 gene), a novel p.Ala98Cys variant in the ATP-binding cassette subfamily G member 5 (ABCG5) gene and a first homozygous description of the p.Gln312His in the ABCB11 gene. The p.Gln312His disrupts the interaction pattern of the bile salt export pump as well as the flexibility of the second intracellular loop domain harboring this residue. As for the p.Ala98Cys, it modulates both the interactions within the first nucleotide-binding domain of the bile transporter and its accessibility. Two additional potentially modifier variants in cholestasis-associated genes were retained based on their pathogenicity (p.Gly758Val in the ABCC2 gene) and functionality (p.Asp19His in the ABCG8 gene). Molecular findings allowed a PFIC2 diagnosis in five patients and an unexpected diagnosis of sisterolemia in one case. The absence of genotype/phenotype correlation suggests the implication of environmental and epigenetic factors as well as modifier variants involved directly or indirectly in the bile composition, which could explain the cholestasis phenotypic variability.

Cyclodipeptide oxidase is an enzyme filament

Modified cyclic dipeptides represent a widespread class of secondary metabolites with diverse pharmacological activities, including antibacterial, antifungal, and antitumor. Here, we report the structural characterization of the Streptomyces noursei enzyme AlbAB, a cyclodipeptide oxidase (CDO) carrying out α,β-dehydrogenations during the biosynthesis of the antibiotic albonoursin. We show that AlbAB is a megadalton heterooligomeric enzyme filament containing covalently bound flavin mononucleotide cofactors. We highlight that AlbAB filaments consist of alternating dimers of AlbA and AlbB and that enzyme activity is crucially dependent on filament formation. We show that AlbA-AlbB interactions are highly conserved suggesting that other CDO-like enzymes are likely enzyme filaments. As CDOs have been employed in the structural diversification of cyclic dipeptides, our results will be useful for future applications of CDOs in biocatalysis and chemoenzymatic synthesis.

Comprehensive in silico analysis of prolactin receptor (PRLR) gene nonsynonymous single nucleotide polymorphisms (nsSNPs) reveals multifaceted impact on protein structure, function, and interactions

This study delves into the functional and structural implications of non-synonymous single nucleotide polymorphisms (nsSNPs) within the Prolactin Receptor (PRLR) gene. Thirteen deleterious nsSNPs were identified through bioinformatics tools, with SIFT predicting 168 out of 395 nsSNPs as detrimental, exhibiting tolerance index (TI) scores ranging from 0 to 0.05. Polyphen2 assigned likelihood scores >0.99 to all 13 nsSNPs, indicating high probability of harm, while Panther scores classified most nsSNPs as 'probably damaging', with specific mutations like W218R scoring 0.74, suggesting a higher impact. Stability analysis using DDG I-Mutant and DDG Mupro consistently predicted decreased stability for all mutations, with CUPSAT indicating mutations like V125G and W218R significantly decreasing stability. Structural analysis through DynaMut predicted destabilization for all mutations except L196I and L292H. MutPred2 highlighted structural alterations for all nsSNPs except L196I, L293V, R315W, and S353N. Domain analysis revealed key mutations within essential functional domains, with five nsSNPs located within Fibronectin type-III domains. Bayesian analysis through ConSurf identified 9 critical residues, with 11 nsSNPs exhibiting notably high conservation. STRING analysis unveiled a complex interaction network, indicating involvement in vital biological processes like lactation. Molecular dynamics (MD) simulations, spanning 100 nanoseconds, elucidated structural dynamics induced by detrimental missense SNPs. Post-translational modification (PTM) analysis identified specific mutations, such as R351, involved in methylation, while S353 was implicated in phosphorylation and glycosylation. These findings offer comprehensive insights into the molecular and phenotypic effects of deleterious nsSNPs in the PRLR gene, crucial for selective breeding.Communicated by Ramaswamy H. Sarma.

The biosynthesis of the odorant 2-methylisoborneol is compartmentalized inside a protein shell

Terpenoids are the largest class of natural products, found across all domains of life. One of the most abundant bacterial terpenoids is the volatile odorant 2-methylisoborneol (2-MIB), partially responsible for the earthy smell of soil and musty taste of contaminated water. Many bacterial 2-MIB biosynthetic gene clusters were thought to encode a conserved transcription factor, named EshA in the model soil bacterium Streptomyces griseus. Here, we revise the function of EshA, now referred to as Sg Enc, and show that it is a Family 2B encapsulin shell protein. Using cryo-electron microscopy, we find that Sg Enc forms an icosahedral protein shell and encapsulates 2-methylisoborneol synthase (2-MIBS) as a cargo protein. Sg Enc contains a cyclic adenosine monophosphate (cAMP) binding domain (CBD)-fold insertion and a unique metal-binding domain, both displayed on the shell exterior. We show that Sg Enc CBDs do not bind cAMP. We find that 2-MIBS cargo loading is mediated by an N-terminal disordered cargo-loading domain and that 2-MIBS activity and Sg Enc shell structure are not modulated by cAMP. Our work redefines the function of EshA and establishes Family 2B encapsulins as cargo-loaded protein nanocompartments involved in secondary metabolite biosynthesis.

PySupercharge: a python algorithm for enabling ABC transporter bacterial secretion of all proteins through amino acid mutation

BACKGROUND

The process of producing proteins in bacterial systems and secreting them through ATP-binding cassette (ABC) transporters is an area that has been actively researched and used due to its high protein production capacity and efficiency. However, some proteins are unable to pass through the ABC transporter after synthesis, a phenomenon we previously determined to be caused by an excessive positive charge in certain regions of their amino acid sequence. If such an excessive charge is removed, the secretion of any protein through ABC transporters becomes possible.

RESULTS

In this study, we introduce 'linear charge density' as the criteria for possibility of protein secretion through ABC transporters and confirm that this criterion can be applied to various non-secretable proteins, such as SARS-CoV-2 spike proteins, botulinum toxin light chain, and human growth factors. Additionally, we develop a new algorithm, PySupercharge, that enables the secretion of proteins containing regions with high linear charge density. It selectively converts positively charged amino acids into negatively charged or neutral amino acids after linear charge density analysis to enable protein secretion through ABC transporters.

CONCLUSIONS

PySupercharge, which also minimizes functional/structural stability loss of the pre-mutation proteins through the use of sequence conservation data, is currently being operated on an accessible web server. We verified the efficacy of PySupercharge-driven protein supercharging by secreting various previously non-secretable proteins commonly used in research, and so suggest this tool for use in future research requiring effective protein production.

The CD4 transmembrane GGXXG and juxtamembrane (C/F)CV+C motifs mediate pMHCII-specific signaling independently of CD4-LCK interactions

CD4+ T cell activation is driven by five-module receptor complexes. The T cell receptor (TCR) is the receptor module that binds composite surfaces of peptide antigens embedded within MHCII molecules (pMHCII). It associates with three signaling modules (CD3γε, CD3δε, and CD3ζζ) to form TCR-CD3 complexes. CD4 is the coreceptor module. It reciprocally associates with TCR-CD3-pMHCII assemblies on the outside of a CD4+ T cells and with the Src kinase, LCK, on the inside. Previously, we reported that the CD4 transmembrane GGXXG and cytoplasmic juxtamembrane (C/F)CV+C motifs found in eutherian (placental mammal) CD4 have constituent residues that evolved under purifying selection (Lee et al., 2022). Expressing mutants of these motifs together in T cell hybridomas increased CD4-LCK association but reduced CD3ζ, ZAP70, and PLCγ1 phosphorylation levels, as well as IL-2 production, in response to agonist pMHCII. Because these mutants preferentially localized CD4-LCK pairs to non-raft membrane fractions, one explanation for our results was that they impaired proximal signaling by sequestering LCK away from TCR-CD3. An alternative hypothesis is that the mutations directly impacted signaling because the motifs normally play an LCK-independent role in signaling. The goal of this study was to discriminate between these possibilities. Using T cell hybridomas, our results indicate that: intracellular CD4-LCK interactions are not necessary for pMHCII-specific signal initiation; the GGXXG and (C/F)CV+C motifs are key determinants of CD4-mediated pMHCII-specific signal amplification; the GGXXG and (C/F)CV+C motifs exert their functions independently of direct CD4-LCK association. These data provide a mechanistic explanation for why residues within these motifs are under purifying selection in jawed vertebrates. The results are also important to consider for biomimetic engineering of synthetic receptors.

Modulation of peroxisomal import by the PEX13 SH3 domain and a proximal FxxxF binding motif

Import of proteins into peroxisomes depends on PEX5, PEX13 and PEX14. By combining biochemical methods and structural biology, we show that the C-terminal SH3 domain of PEX13 mediates intramolecular interactions with a proximal FxxxF motif. The SH3 domain also binds WxxxF peptide motifs in the import receptor PEX5, demonstrating evolutionary conservation of such interactions from yeast to human. Strikingly, intramolecular interaction of the PEX13 FxxxF motif regulates binding of PEX5 WxxxF/Y motifs to the PEX13 SH3 domain. Crystal structures reveal how FxxxF and WxxxF/Y motifs are recognized by a non-canonical surface on the SH3 domain. The PEX13 FxxxF motif also mediates binding to PEX14. Surprisingly, the potential PxxP binding surface of the SH3 domain does not recognize PEX14 PxxP motifs, distinct from its yeast ortholog. Our data show that the dynamic network of PEX13 interactions with PEX5 and PEX14, mediated by diaromatic peptide motifs, modulates peroxisomal matrix import.

Structural bioinformatics studies of glutamate transporters and their AlphaFold2 predicted water-soluble QTY variants and uncovering the natural mutations of L->Q, I->T, F->Y and Q->L, T->I and Y->F

Glutamate transporters play key roles in nervous physiology by modulating excitatory neurotransmitter levels, when malfunctioning, involving in a wide range of neurological and physiological disorders. However, integral transmembrane proteins including the glutamate transporters remain notoriously difficult to study, due to their localization within the cell membrane. Here we present the structural bioinformatics studies of glutamate transporters and their water-soluble variants generated through QTY-code, a protein design strategy based on systematic amino acid substitutions. These include 2 structures determined by X-ray crystallography, cryo-EM, and 6 predicted by AlphaFold2, and their predicted water-soluble QTY variants. In the native structures of glutamate transporters, transmembrane helices contain hydrophobic amino acids such as leucine (L), isoleucine (I), and phenylalanine (F). To design water-soluble variants, these hydrophobic amino acids are systematically replaced by hydrophilic amino acids, namely glutamine (Q), threonine (T) and tyrosine (Y). The QTY variants exhibited water-solubility, with four having identical isoelectric focusing points (pI) and the other four having very similar pI. We present the superposed structures of the native glutamate transporters and their water-soluble QTY variants. The superposed structures displayed remarkable similarity with RMSD 0.528Å-2.456Å, despite significant protein transmembrane sequence differences (41.1%->53.8%). Additionally, we examined the differences of hydrophobicity patches between the native glutamate transporters and their QTY variants. Upon closer inspection, we discovered multiple natural variations of L->Q, I->T, F->Y and Q->L, T->I, Y->F in these transporters. Some of these natural variations were benign and the remaining were reported in specific neurological disorders. We further investigated the characteristics of hydrophobic to hydrophilic substitutions in glutamate transporters, utilizing variant analysis and evolutionary profiling. Our structural bioinformatics studies not only provided insight into the differences between the hydrophobic helices and hydrophilic helices in the glutamate transporters, but they are also expected to stimulate further study of other water-soluble transmembrane proteins.

Probing the CRL4DCAF12 interactions with MAGEA3 and CCT5 di-Glu C-terminal degrons

Damaged DNA-binding protein-1 (DDB1)- and CUL4-associated factor 12 (DCAF12) serves as the substrate recognition component within the Cullin4-RING E3 ligase (CRL4) complex, capable of identifying C-terminal double-glutamic acid degrons to promote the degradation of specific substrates through the ubiquitin proteasome system. Melanoma-associated antigen 3 (MAGEA3) and T-complex protein 1 subunit epsilon (CCT5) proteins have been identified as cellular targets of DCAF12. To further characterize the interactions between DCAF12 and both MAGEA3 and CCT5, we developed a suite of biophysical and proximity-based cellular NanoBRET assays showing that the C-terminal degron peptides of both MAGEA3 and CCT5 form nanomolar affinity interactions with DCAF12 in vitro and in cells. Furthermore, we report here the 3.17 Å cryo-EM structure of DDB1-DCAF12-MAGEA3 complex revealing the key DCAF12 residues responsible for C-terminal degron recognition and binding. Our study provides new insights and tools to enable the discovery of small molecule handles targeting the WD40-repeat domain of DCAF12 for future proteolysis targeting chimera design and development.

Molecular insights into capsular polysaccharide secretion

Capsular polysaccharides (CPSs) fortify the cell boundaries of many commensal and pathogenic bacteria1. Through the ABC-transporter-dependent biosynthesis pathway, CPSs are synthesized intracellularly on a lipid anchor and secreted across the cell envelope by the KpsMT ABC transporter associated with the KpsE and KpsD subunits1,2. Here we use structural and functional studies to uncover crucial steps of CPS secretion in Gram-negative bacteria. We show that KpsMT has broad substrate specificity and is sufficient for the translocation of CPSs across the inner bacterial membrane, and we determine the cell surface organization and localization of CPSs using super-resolution fluorescence microscopy. Cryo-electron microscopy analyses of the KpsMT-KpsE complex in six different states reveal a KpsE-encaged ABC transporter, rigid-body conformational rearrangements of KpsMT during ATP hydrolysis and recognition of a glycolipid inside a membrane-exposed electropositive canyon. In vivo CPS secretion assays underscore the functional importance of canyon-lining basic residues. Combined, our analyses suggest a molecular model of CPS secretion by ABC transporters.

Conserved proline residues prevent dimerization and aggregation in the β-lactamase BlaC

Evolution leads to conservation of amino acid residues in protein families. Conserved proline residues are usually considered to ensure the correct folding and to stabilize the three-dimensional structure. Surprisingly, proline residues that are highly conserved in class A β-lactamases were found to tolerate various substitutions without large losses in enzyme activity. We investigated the roles of three conserved prolines at positions 107, 226, and 258 in the β-lactamase BlaC from Mycobacterium tuberculosis and found that mutations can lead to dimerization of the enzyme and an overall less stable protein that is prone to aggregate over time. For the variant Pro107Thr, the crystal structure shows dimer formation resembling domain swapping. It is concluded that the proline substitutions loosen the structure, enhancing multimerization. Even though the enzyme does not lose its properties without the conserved proline residues, the prolines ensure the long-term structural integrity of the enzyme.

Modulating the pH profile of vanillin dehydrogenase enzyme from extremophile Bacillus ligniniphilus L1 through computational guided site-directed mutagenesis

Vanillin dehydrogenase (VDH) has recently come forward as an important enzyme for the commercial production of vanillic acid from vanillin in a one-step enzymatic process. However, VDH with high alkaline tolerance and efficiency is desirable to meet the biorefinery requirements. In this study, computationally guided site-directed mutagenesis was performed by increasing the positive and negative charges on the surface and near the active site of the VDH from the alkaliphilic marine bacterium Bacillus ligniniphilus L1, respectively. In total, 20 residues including 15 from surface amino acids and 5 near active sites were selected based on computational analysis and were subjected to site-directed mutations. The optimum pH of the two screened mutants including I132R, and T235E from surface residue and near active site mutant was shifted to 9, and 8.6, with a 2.82- and 2.95-fold increase in their activity compared to wild enzyme at pH 9, respectively. A double mutant containing both these mutations i.e., I132R/T235E was produced which showed a shift in optimum pH of VDH from 7.4 to 9, with an increase of 74.91 % in enzyme activity. Therefore, the double mutant of VDH from the L1 strain (I132R/T235E) produced in this study represents a potential candidate for industrial applications.

An in-silico analysis of OGT gene association with diabetes mellitus

O-GlcNAcylation is a nutrient-sensing post-translational modification process. This cycling process involves two primary proteins: the O-linked N-acetylglucosamine transferase (OGT) catalysing the addition, and the glycoside hydrolase OGA (O-GlcNAcase) catalysing the removal of the O-GlCNAc moiety on nucleocytoplasmic proteins. This process is necessary for various critical cellular functions. The O-linked N-acetylglucosamine transferase (OGT) gene produces the OGT protein. Several studies have shown the overexpression of this protein to have biological implications in metabolic diseases like cancer and diabetes mellitus (DM). This study retrieved 159 SNPs with clinical significance from the SNPs database. We probed the functional effects, stability profile, and evolutionary conservation of these to determine their fit for this research. We then identified 7 SNPs (G103R, N196K, Y228H, R250C, G341V, L367F, and C845S) with predicted deleterious effects across the four tools used (PhD-SNPs, SNPs&Go, PROVEAN, and PolyPhen2). Proceeding with this, we used ROBETTA, a homology modelling tool, to model the proteins with these point mutations and carried out a structural bioinformatics method- molecular docking- using the Glide model of the Schrodinger Maestro suite. We used a previously reported inhibitor of OGT, OSMI-1, as the ligand for these mutated protein models. As a result, very good binding affinities and interactions were observed between this ligand and the active site residues within 4Å of OGT. We conclude that these mutation points may be used for further downstream analysis as drug targets for treating diabetes mellitus.