ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules

De novo TLK1 and MDM1 mutations in a patient with a neurodevelopmental disorder and immunodeficiency

The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants underlie the neurodevelopmental disorder (NDD) 'Intellectual Disability, Autosomal Dominant 57' (MRD57), characterized by intellectual disability and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown. A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a heterozygous TLK1 variant (c.1435C>G, p.Q479E), as well as a mutation in MDM1 (c.1197dupT, p.K400∗). Cells expressing TLK1 p.Q479E exhibited reduced cytokine responses and elevated DNA damage, but not increased radiation sensitivity or DNA repair defects. The TLK1 p.Q479E variant impaired kinase activity but not proximal protein interactions. Our study provides the first functional characterization of NDD-associated TLK1 variants and suggests that, such as TLK2, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.

A non-symmetrical p97 conformation initiates a multistep recruitment of Ufd1/Npl4

In vitro experiments and cryo-EM structures of p97 and its cofactor, Ufd1/Npl4 (UN), elucidated substrate processing. Yet, the structural transitions and the related ATPase cycle upon UN binding remain unresolved. We captured two discrete conformations: One in which D1 protomers are ATP bound, while the D2 subunits are in the ADP state, presumably required for substrate engagement with the D2 pore; and a heterologous nucleotide state within the D1 ring in which only two NTDs are in the "up" ATP state that favors UN binding. Further analysis suggests that initially, UN binds p97's non-symmetrical conformation, this association promotes a structural transition upon which five NTDs shift to an "up" state and are poised to bind ATP. The UBXL domain of Npl4 was captured bound to an NTD in the ADP state, demonstrating a conformation that may provide directionality to incoming substrate and introduce the flexibility needed for substrate processing.

Polymorphism in the leucine-rich repeats of TLR7 in different breeds of chicken and in silico analysis of its effect on TLR7 structure and function

Chicken toll-like receptor 7 (chTLR7) is a viral sensing pattern recognition receptor and detects ssRNA. The ligand binding site comprises leucine-rich repeats (LRRs) located in the ectodomain of chTLR7. Hence, any polymorphism in the binding site would modify its functional interaction with the ligand, resulting in varied strength of immune response. This study first aimed to compare the single nucleotide polymorphisms (SNPs) associated with the ligand binding site of TLR7 in three indigenous chicken breeds namely Aseel, Kadaknath, Nicobari along with an exotic breed White Leghorn. Four synonymous SNPs (P123P, I171I, N339N and L421L) and four non-synonymous SNPs (I121V, S135T, F356S and S447G) were identified among various breeds. We employed in silico tools to screen the pathogenic nsSNPs and one nsSNP was identified as having potential impact on chTLR7 protein. Moreover, sequence and structure-based methods were used to determine the effect of nsSNPs on protein stability. It revealed I121V, F356S, and S447G as decreasing the stability while S135T increasing the stability of chTLR7. Additionally, docking analysis confirmed that I121V and F356S reduced the binding affinity of ligands (R-848 and polyU) to chTLR7 protein. The results suggest that the nsSNPs found in this study could alter the ligand binding of chTLR7 and modify the immune response between different breeds further contributing to disease susceptibility or resistance. Further, in vitro and in vivo studies are needed to analyze the effect of these SNPs on susceptibility or resistance against various viral diseases in poultry.

Mutations found in cancer patients compromise DNA binding of the winged helix protein STK19

Serine/threonine protein kinase 19 (STK19) has been reported to phosphorylate and activate oncogenic NRAS to promote melanomagenesis. However, concerns have been raised about whether STK19 is a kinase. STK19 has also been identified as a putative factor involved in the transcription-coupled nucleotide excision repair (TC-NER) pathway. In this study, we determined the 1.32 Å crystal structure of human STK19. The structure reveals that STK19 is a winged helix (WH) protein consisting of three tandem WH domains. STK19 binds more strongly to double-stranded DNA and RNA (dsDNA/dsRNA) than to ssDNA. A positively charged patch centered on helix WH3-H1 contributes to dsDNA binding, which is unusual because the WH domain typically uses helix H3 as the recognition helix. Importantly, mutations of the conserved residues in the basic patch, K186N, R200W, and R215W, are found in cancer patients, and these mutations compromise STK19 DNA binding. Other mutations have been predicted to produce a similar effect, including two mutations that disrupt the nuclear localization signal (NLS) motif. These mutations may indirectly impact the DNA binding capacity of STK19 by interfering with its nuclear localization.

SCAPER-Related Autosomal Recessive Retinitis Pigmentosa with Intellectual Disability: Confirming and Extending the Phenotypic Spectrum and Bioinformatics Analyses

Mutations in the gene SCAPER (S phase Cyclin A-Associated Protein residing in the Endoplasmic Reticulum) have recently been associated with retinitis pigmentosa (RP) and intellectual disability (ID). In 2011, a possible involvement of SCAPER in human diseases was discovered for the first time due to the identification of a homozygous mutation causing ID in an Iranian family. Later, five studies were published in 2019 that described patients with autosomal recessive syndromic retinitis pigmentosa (arRP) accompanied by ID and attention-deficit/hyperactivity disorder (ADHD). This present study describes three patients from an Arab consanguineous family in Israel with similar clinical features of the SCAPER syndrome. In addition, new manifestations of ocular symptoms, nystagmus, glaucoma, and elevator palsy, were observed. Genetic testing of the patients and both parents via whole-exome sequencing revealed the homozygous mutation c.2023-2A>G in SCAPER. Phenotypic and genotypic descriptions for all available cases described in the literature including our current three cases (37 cases) were carried out, in addition to a bioinformatics analysis for all the genetic variants that was undertaken. Our study confirms and extends the clinical manifestations of SCAPER-related disorders.

Structure, regulation, and mechanisms of nonmuscle myosin-2

Members of the myosin superfamily of molecular motors are large mechanochemical ATPases that are implicated in an ever-expanding array of cellular functions. This review focuses on mammalian nonmuscle myosin-2 (NM2) paralogs, ubiquitous members of the myosin-2 family of filament-forming motors. Through the conversion of chemical energy into mechanical work, NM2 paralogs remodel and shape cells and tissues. This process is tightly controlled in time and space by numerous synergetic regulation mechanisms to meet cellular demands. We review how recent advances in structural biology together with elegant biophysical and cell biological approaches have contributed to our understanding of the shared and unique mechanisms of NM2 paralogs as they relate to their kinetics, regulation, assembly, and cellular function.

Identification of high-risk non-synonymous SNPs (nsSNPs) in DNAH1 and DNAH17 genes associated with male infertility: a bioinformatics analysis

Male infertility is a significant reproductive issue affecting a considerable number of couples worldwide. While there are various causes of male infertility, genetic factors play a crucial role in its development. We focused on identifying and analyzing the high-risk nsSNPs in DNAH1 and DNAH17 genes, which encode proteins involved in sperm motility. A total of 20 nsSNPs for DNAH1 and 10 nsSNPs for DNAH17 were analyzed using various bioinformatics tools including SIFT, PolyPhen-2, CADD, PhD-SNPg, VEST-4, and MutPred2. As a result, V1287G, L2071R, R2356W, R3169C, R3229C, E3284K, R4096L, R4133C, and A4174T in DNAH1 gene and C1803Y, C1829Y, R1903C, and L3595P in DNAH17 gene were identified as high-risk nsSNPs. These nsSNPs were predicted to decrease protein stability, and almost all were found in highly conserved amino acid positions. Additionally, 4 nsSNPs were observed to alter post-translational modification status. Furthermore, the interaction network analysis revealed that DNAH1 and DNAH17 interact with DNAH2, DNAH3, DNAH5, DNAH7, DNAH8, DNAI2, DNAL1, CFAP70, DNAI3, DNAI4, ODAD1, and DNAI7, demonstrating the importance of DNAH1 and DNAH17 proteins in the overall functioning of the sperm motility machinery. Taken together, these findings revealed the detrimental effects of identified high-risk nsSNPs on protein structure and function and highlighted their potential relevance to male infertility. Further studies are warranted to validate these findings and to elucidate the underlying mechanisms.

CryoEM reveals the structure of an archaeal pilus involved in twitching motility

Amongst the major types of archaeal filaments, several have been shown to closely resemble bacterial homologues of the Type IV pili (T4P). Within Sulfolobales, member species encode for three types of T4P, namely the archaellum, the UV-inducible pilus system (Ups) and the archaeal adhesive pilus (Aap). Whereas the archaellum functions primarily in swimming motility, and the Ups in UV-induced cell aggregation and DNA-exchange, the Aap plays an important role in adhesion and twitching motility. Here, we present a cryoEM structure of the Aap of the archaeal model organism Sulfolobus acidocaldarius. We identify the component subunit as AapB and find that while its structure follows the canonical T4P blueprint, it adopts three distinct conformations within the pilus. The tri-conformer Aap structure that we describe challenges our current understanding of pilus structure and sheds new light on the principles of twitching motility.

Characterization and Structural Analysis of a Novel Carbohydrate-Binding Module from Family 96 with Chondroitin Sulfate-Specific Binding Capacity

Chondroitin sulfate (CS) is the predominant glycosaminoglycan within the human body and is widely applied in various industries. Carbohydrate-binding modules (CBMs) possessing the capacity for carbohydrate recognition are verified to be important tools for polysaccharide investigation. Only one CS-specific CBM, PhCBM100, has hitherto been characterized. In the present study, two CBM96 domains present in the same putative PL8_3 chondroitin AC lyase were discovered and recombinantly expressed. The results of microtiter plate assays and affinity gel electrophoresis assays showed that the two corresponding proteins, DmCBM96-1 and DmCBM96-2, bind specifically to CSs. The crystal structure of DmCBM96-1 was determined at a 2.20 Å resolution. It adopts a β-sandwich fold comprising two antiparallel β-sheets, showing structural similarities to TM6-N4, which is the founding member of the CBM96 family. Site mutagenesis analysis revealed that the residues of Arg27, Lys45, Tyr51, Arg53, and Arg157 are critical for CS binding. The characterization of the two CBM96 proteins demonstrates the diverse ligand specificity of the CBM96 family and provides promising tools for CS investigation.

Computer-driven Evolution of Myrosinase from the Cabbage Aphid for Efficient Production of (R)-Sulforaphane

Myrosinase (Myr) catalyzes the hydrolysis of glucosinolates, yielding biologically active metabolites. In this study, glucoraphanin (GRA) extracted from broccoli seeds was effectively hydrolyzed using a Myr-obtained cabbage aphid (Brevicoryne brassicae) (BbMyr) to produce (R)-sulforaphane (SFN). The gene encoding BbMyr was successfully heterologously expressed in Escherichia coli, resulting in the production of 1.6 g/L (R)-SFN, with a remarkable yield of 20.8 mg/gbroccoli seeds, achieved using recombination E. coli whole-cell catalysis under optimal conditions (pH 4.5, 45 °C). Subsequently, BbMyr underwent combinatorial simulation-driven mutagenesis, yielding a mutant, DE9 (N321D/Y426S), showing a remarkable 2.91-fold increase in the catalytic efficiency (kcat/KM) compared with the original enzyme. Molecular dynamics simulations demonstrated that the N321D mutation in loopA of mutant DE9 enhanced loopA stability by inducing favorable alterations in hydrogen bonds, while the Y426S mutation in loopB decreased spatial resistance. This research lays a foundation for the environmentally sustainable enzymatic (R)-SFN synthesis.

Structure of the human heparan-α-glucosaminide N-acetyltransferase (HGSNAT)

Degradation of heparan sulfate (HS), a glycosaminoglycan (GAG) comprised of repeating units of N-acetylglucosamine and glucuronic acid, begins in the cytosol and is completed in the lysosomes. Acetylation of the terminal non-reducing amino group of a-D-glucosamine of HS is essential for its complete breakdown into monosaccharides and free sulfate. Heparan-a-glucosaminide N-acetyltransferase (HGSNAT), a resident of the lysosomal membrane, catalyzes this essential acetylation reaction by accepting and transferring the acetyl group from cytosolic acetyl-CoA to terminal a-D-glucosamine of HS in the lysosomal lumen. Mutation-induced dysfunction in HGSNAT causes abnormal accumulation of HS within the lysosomes and leads to an autosomal recessive neurodegenerative lysosomal storage disorder called mucopolysaccharidosis IIIC (MPS IIIC). There are no approved drugs or treatment strategies to cure or manage the symptoms of, MPS IIIC. Here, we use cryo-electron microscopy (cryo-EM) to determine a high-resolution structure of the HGSNAT-acetyl-CoA complex, the first step in HGSNAT catalyzed acetyltransferase reaction. In addition, we map the known MPS IIIC mutations onto the structure and elucidate the molecular basis for mutation-induced HGSNAT dysfunction.

Enhancing the Thermal Stability and Enzyme Activity of Ketopantoate Hydroxymethyltransferase through Interface Modification Engineering

Ketopantoate hydroxymethyltransferase (KPHMT) plays a pivotal role in d-pantothenic acid biosynthesis. Most KPHMTs are homodecamers with low thermal stability, posing challenges for protein engineering and limiting output enhancement. Previously, a high-enzyme activity KPHMT mutant (K25A/E189S) from Corynebacterium glutamicum was screened as mother strain (M0). Building upon this strain, our study focused on interface engineering modifications, employing a multifaceted approach including integrating folding-free energy calculation, B-factor analysis, and conserved site analysis. Preliminary screening led to the selection of five mutants in the interface─E106S, E98T, E98N, S247I, and S247D─showing improved thermal stability, culminating in the double-site mutant M8 (M0-E98N/S247D). M8 exhibited a T1/2 value of 288.79 min at 50 °C, showing a 3.29-fold increase compared to M0. Meanwhile, the Tm value of M8 was elevated from 53.2 to 59.6 °C. Investigations of structural and molecular dynamics simulations revealed alterations in surface electrostatic charge distribution and the formation of increased hydrogen bonds between subunits, contributing to enhanced thermal stability. This investigation corroborates the efficacy of interface engineering modifications in bolstering KPHMT stability while showing its potential for positively impacting industrial d-pantothenic acid synthesis.

Prediction of the most deleterious non-synonymous SNPs in the human IL1B gene: evidence from bioinformatics analyses

BACKGROUND

Polymorphisms in IL1B play a significant role in depression, multiple inflammatory-associated disorders, and susceptibility to infection. Functional non-synonymous SNPs (nsSNPs) result in changes in the encoded amino acids, potentially leading to structural and functional alterations in the mutant proteins. So far, most genetic studies have concentrated on SNPs located in the IL1B promoter region, without addressing nsSNPs and their association with multifactorial diseases. Therefore, this study aimed to explore the impact of deleterious nsSNPs retrieved from the dbSNP database on the structure and functions of the IL1B protein.

RESULTS

Six web servers (SIFT, PolyPhen-2, PROVEAN, SNPs&GO, PHD-SNP, PANTHER) were used to analyze the impact of 222 missense SNPs on the function and structure of IL1B protein. Five novel nsSNPs (E100K, T240I, S53Y, D128Y, and F228S) were found to be deleterious and had a mutational impact on the structure and function of the IL1B protein. The I-mutant v2.0 and MUPro servers predicted that these mutations decreased the stability of the IL1B protein. Additionally, these five mutations were found to be conserved, underscoring their significance in protein structure and function. Three of them (T240I, D128Y, and F228S) were predicted to be cancer-causing nsSNPs. To analyze the behavior of the mutant structures under physiological conditions, we conducted a 50 ns molecular dynamics simulation using the WebGro online tool. Our findings indicate that the mutant values differ from those of the IL1B wild type in terms of RMSD, RMSF, Rg, SASA, and the number of hydrogen bonds.

CONCLUSIONS

This study provides valuable insights into nsSNPs located in the coding regions of IL1B, which lead to direct deleterious effects on the functional and structural aspects of the IL1B protein. Thus, these nsSNPs could be considered significant candidates in the pathogenesis of disorders caused by IL1B dysfunction, contributing to effective drug discovery and the development of precision medications. Thorough research and wet lab experiments are required to verify our findings. Moreover, bioinformatic tools were found valuable in the prediction of deleterious nsSNPs.

Developing deprotectase biocatalysts for synthesis

Organic synthesis often requires multiple steps where a functional group (FG) is concealed from reaction by a protecting group (PG). Common PGs include N-carbobenzyloxy (Cbz or Z) of amines and tert-butyloxycarbonyl (OtBu) of acids. An essential step is the removal of the PG, but this often requires excess reagents, extensive time and can have low % yield. An overarching goal of biocatalysis is to use "green" or "enzymatic" methods to catalyse chemical transformations. One under-utilised approach is the use of "deprotectase" biocatalysts to selectively remove PGs from various organic substrates. The advantage of this methodology is the exquisite selectivity of the biocatalyst to only act on its target, leaving other FGs and PGs untouched. A number of deprotectase biocatalysts have been reported but they are not commonly used in mainstream synthetic routes. This study describes the construction of a cascade to deprotect doubly-protected amino acids. The well known Bacillus BS2 esterase was used to remove the OtBu PG from various amino acid substrates. The more obscure Sphingomonas Cbz-ase (amidohydrolase) was screened with a range of N-Cbz-modified amino acid substrates. We then combined both the BS2 and Cbz-ase together for a 1 pot, 2 step deprotection of the model substrate CBz-L-Phe OtBu to produce the free L-Phe. We also provide some insight into the residues involved in substrate recognition and catalysis using docked ligands in the crystal structure of BS2. Similarly, a structural model of the Cbz-ase identifies a potential di-metal binding site and reveals conserved active site residues. This new biocatalytic cascade should be further explored for its application in chemical synthesis.

Retron-Eco1 assembles NAD+-hydrolyzing filaments that provide immunity against bacteriophages

Retrons are toxin-antitoxin systems protecting bacteria against bacteriophages via abortive infection. The Retron-Eco1 antitoxin is formed by a reverse transcriptase (RT) and a non-coding RNA (ncRNA)/multi-copy single-stranded DNA (msDNA) hybrid that neutralizes an uncharacterized toxic effector. Yet, the molecular mechanisms underlying phage defense remain unknown. Here, we show that the N-glycosidase effector, which belongs to the STIR superfamily, hydrolyzes NAD+ during infection. Cryoelectron microscopy (cryo-EM) analysis shows that the msDNA stabilizes a filament that cages the effector in a low-activity state in which ADPr, a NAD+ hydrolysis product, is covalently linked to the catalytic E106 residue. Mutations shortening the msDNA induce filament disassembly and the effector's toxicity, underscoring the msDNA role in immunity. Furthermore, we discovered a phage-encoded Retron-Eco1 inhibitor (U56) that binds ADPr, highlighting the intricate interplay between retron systems and phage evolution. Our work outlines the structural basis of Retron-Eco1 defense, uncovering ADPr's pivotal role in immunity.

Novel molecular, structural and clinical findings in an Italian cohort of congenital cataract

The current genetic diagnostic workup of congenital cataract (CC) is mainly based on NGS panels, whereas exome sequencing (ES) has occasionally been employed. In this multicentre study, we investigated by ES the detection yield, mutational spectrum and genotype-phenotype correlations in a CC cohort recruited between 2020 and mid-2022. The cohort consisted of 67 affected individuals from 51 unrelated families and included both non-syndromic (75%) and syndromic (25%) phenotypes, with extra-CC ocular/visual features present in both groups (48% and 76%, respectively). The functional effect of variants was predicted by 3D modelling and hydropathy properties changes. Variant clustering was used for the in-depth assessment of genotype-phenotype correlations. A diagnostic (pathogenic or likely pathogenic) variant was identified in 19 out of 51 probands/families (~37%). In a further 14 probands/families a candidate variant was identified: in 12 families a VUS was detected, of which 9 were considered plausibly pathogenic (i.e., 4 or 5 points according to ACMG criteria), while in 2 probands ES identified a single variant in an autosomal recessive gene associated with CC. Eighteen probands/families, manifesting primarily non-syndromic CC (15/18, 83%), remained unsolved. The identified variants (8 P, 12 LP, 10 VUS-PP, and 5 VUS), half of which were unreported in the literature, affected five functional categories of genes involved in transcription/splicing, lens formation/homeostasis (i.e., crystallin genes), membrane signalling, cell-cell interaction, and immune response. A phenotype-specific variant clustering was observed in four genes (KIF1A, MAF, PAX6, SPTAN1), whereas variable expressivity and potential phenotypic expansion in two (BCOR, NHS) and five genes (CWC27, KIF1A, IFIH1, PAX6, SPTAN1), respectively. Finally, ES allowed to detect variants in six genes not commonly included in commercial CC panels. These findings broaden the genotype-phenotype correlations in one of the largest CC cohorts tested by ES, providing novel insights into the underlying pathogenetic mechanisms and emphasising the power of ES as first-tier test.

Multi-omics reveals new links between Fructosamine-3-Kinase (FN3K) and core metabolic pathways

Fructosamine-3-kinases (FN3Ks) are a conserved family of repair enzymes that phosphorylate reactive sugars attached to lysine residues in peptides and proteins. Although FN3Ks are present across the Tree of Life and share detectable sequence similarity to eukaryotic protein kinases, the biological processes regulated by these kinases are largely unknown. To address this knowledge gap, we leveraged the FN3K CRISPR Knock-Out (KO) HepG2 cell line alongside an integrative multi-omics study combining transcriptomics, metabolomics, and interactomics to place these enzymes in a pathway context. The integrative analyses revealed the enrichment of pathways related to oxidative stress response, lipid biosynthesis (cholesterol and fatty acids), and carbon and co-factor metabolism. Moreover, enrichment of nicotinamide adenine dinucleotide (NAD) binding proteins and localization of human FN3K (HsFN3K) to mitochondria suggests potential links between FN3K and NAD-mediated energy metabolism and redox balance. We report specific binding of HsFN3K to NAD compounds in a metal and concentration-dependent manner and provide insight into their binding mode using modeling and experimental site-directed mutagenesis. Our studies provide a framework for targeting these understudied kinases in diabetic complications and metabolic disorders where redox balance and NAD-dependent metabolic processes are altered.

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.

Decoding the Structure-Function Relationship of the Muramidase Domain in E. coli O157.H7 Bacteriophage Endolysin: A Potential Building Block for Chimeric Enzybiotics

Bacteriophage endolysins are potential alternatives to conventional antibiotics for treating multidrug-resistant gram-negative bacterial infections. However, their structure-function relationships are poorly understood, hindering their optimization and application. In this study, we focused on the individual functionality of the C-terminal muramidase domain of Gp127, a modular endolysin from E. coli O157:H7 bacteriophage PhaxI. This domain is responsible for the enzymatic activity, whereas the N-terminal domain binds to the bacterial cell wall. Through protein modeling, docking experiments, and molecular dynamics simulations, we investigated the activity, stability, and interactions of the isolated C-terminal domain with its ligand. We also assessed its expression, solubility, toxicity, and lytic activity using the experimental data. Our results revealed that the C-terminal domain exhibits high activity and toxicity when tested individually, and its expression is regulated in different hosts to prevent self-destruction. Furthermore, we validated the muralytic activity of the purified refolded protein by zymography and standardized assays. These findings challenge the need for the N-terminal binding domain to arrange the active site and adjust the gap between crucial residues for peptidoglycan cleavage. Our study shed light on the three-dimensional structure and functionality of muramidase endolysins, thereby enriching the existing knowledge pool and laying a foundation for accurate in silico modeling and the informed design of next-generation enzybiotic treatments.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

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.

A comprehensive examination of ACE2 receptor and prediction of spike glycoprotein and ACE2 interaction based on in silico analysis of ACE2 receptor

The ACE2 receptor plays a vital role not only in the SARS-CoV-induced epidemic but also in various other diseases, including cardiovascular diseases and ARDS. While studies have explored the interactions between ACE2 and SARS-CoV proteins, comprehensive research utilizing bioinformatic tools on the ACE2 protein has been lacking. The one aim of present study was to extensively analyze the regions of the ACE2 protein. After utilizing all bioinformatics tools especially G104 and L108 regions on ACE2 were come forward. The results of our analysis revealed that possible mutations or deletions in the G104 and L108 regions play a critical role in both the biological functioning and the determination of the chemical-physical properties of ACE2. Additionally, these regions were found to be more susceptible to mutations or deletions compared to other regions of the ACE2 protein. Notably, the randomly selected peptide, LQQNGSSVLS (100-109), which includes G104 and L108, exhibited a crucial role in binding the RBD of the spike protein, as supported by docking scores. Furthermore, both MDs and iMODs results provided evidence that G104 and L108 influence the dynamics of ACE2-spike complexes. This study is expected to offer a new perspective on the ACE2-SARS-CoV interaction and other research areas where ACE2 plays a significant role, such as biotechnology (protein engineering, enzyme optimization), medicine (RAS, pulmonary and cardiac diseases), and basic research (structural motifs, stabilizing protein folds, or facilitating important inter molecular contacts, protein's proper structure and function).Communicated by Ramaswamy H. Sarma.

Structural organization of the retriever-CCC endosomal recycling complex

The recycling of membrane proteins from endosomes to the cell surface is vital for cell signaling and survival. Retriever, a trimeric complex of vacuolar protein-sorting-associated protein (VPS)35L, VPS26C and VPS29, together with the CCC complex comprising coiled-coil domain-containing (CCDC)22, CCDC93 and copper metabolism domain-containing (COMMD) proteins, plays a crucial role in this process. The precise mechanisms underlying retriever assembly and its interaction with CCC have remained elusive. Here, we present a high-resolution structure of retriever in humans determined using cryogenic electron microscopy. The structure reveals a unique assembly mechanism, distinguishing it from its remotely related paralog retromer. By combining AlphaFold predictions and biochemical, cellular and proteomic analyses, we further elucidate the structural organization of the entire retriever-CCC complex across evolution and uncover how cancer-associated mutations in humans disrupt complex formation and impair membrane protein homeostasis. These findings provide a fundamental framework for understanding the biological and pathological implications associated with retriever-CCC-mediated endosomal recycling.

Fructose-1-kinase has pleiotropic roles in Escherichia coli

In Escherichia coli, the master transcription regulator catabolite repressor activator (Cra) regulates >100 genes in central metabolism. Cra binding to DNA is allosterically regulated by binding to fructose-1-phosphate (F-1-P), but the only documented source of F-1-P is from the concurrent import and phosphorylation of exogenous fructose. Thus, many have proposed that fructose-1,6-bisphosphate (F-1,6-BP) is also a physiological regulatory ligand. However, the role of F-1,6-BP has been widely debated. Here, we report that the E. coli enzyme fructose-1-kinase (FruK) can carry out its "reverse" reaction under physiological substrate concentrations to generate F-1-P from F-1,6-BP. We further show that FruK directly binds Cra with nanomolar affinity and forms higher order, heterocomplexes. Growth assays with a ΔfruK strain and fruK complementation show that FruK has a broader role in metabolism than fructose catabolism. Since fruK itself is repressed by Cra, these newly-reported events add layers to the dynamic regulation of E. coli's central metabolism that occur in response to changing nutrients. These findings might have wide-spread relevance to other γ-proteobacteria, which conserve both Cra and FruK.

Mechanistic insights into the conformational changes and alterations in residual communications due to the mutations in the pncA Gene of Mycobacterium tuberculosis: A computational perspective for effective therapeutic solutions

Due to its emerging resistance to first-line anti-TB medications, tuberculosis (TB) is one of the most contagious illness in the world. According to reports, the effectiveness of treating TB is severely impacted by drug resistance, notably resistance caused by mutations in the pncA gene-encoded pyrazinamidase (PZase) to the front-line drug pyrazinamide (PZA). The present study focused on investigating the resistance mechanism caused by the mutations D12N, T47A, and H137R to better understand the structural and molecular events responsible for the resistance acquired by the pncA gene of Mycobacterium tuberculosis (MTB) at the structural level. Bioinformatics analysis predicted that all three mutations were deleterious and located near the active centre of the pncA, affecting its functional activity. Furthermore, molecular dynamics simulation (MDS) results established that mutations significantly reduced the structural stability and caused the rearrangement of FE2+ in the active centre of pncA. Moreover, essential dynamics analysis, including principal component analysis (PCA) and free energy landscape (FEL), concluded variations in the protein motion and decreased conformational space in the mutants. Additionally, the mutations potentially impacted the network topologies and altered the residual communications in the network. The complex simulation study results established the significant movement of the flap region from the active centre of mutant complexes, further supporting the flap region's significance in developing resistance to the PZA drug. This study advances our knowledge of the primary cause of the mechanism of PZA resistance and the structural dynamics of pncA mutants, which will help us to design new and potent chemical scaffolds to treat drug-resistant TB (DR-TB).

Structural Basis for the Interaction of Redβ Single-Strand Annealing Protein with Escherichia coli Single-Stranded DNA-Binding Protein

Redβ is a protein from bacteriophage λ that binds to single-stranded DNA (ssDNA) to promote the annealing of complementary strands. Together with λ-exonuclease (λ-exo), Redβ is part of a two-component DNA recombination system involved in multiple aspects of genome maintenance. The proteins have been exploited in powerful methods for bacterial genome engineering in which Redβ can anneal an electroporated oligonucleotide to a complementary target site at the lagging strand of a replication fork. Successful annealing in vivo requires the interaction of Redβ with E. coli single-stranded DNA-binding protein (SSB), which coats the ssDNA at the lagging strand to coordinate access of numerous replication proteins. Previous mutational analysis revealed that the interaction between Redβ and SSB involves the C-terminal domain (CTD) of Redβ and the C-terminal tail of SSB (SSB-Ct), the site for binding of numerous host proteins. Here, we have determined the x-ray crystal structure of Redβ CTD in complex with a peptide corresponding to the last nine residues of SSB (MDFDDDIPF). Formation of the complex is predominantly mediated by hydrophobic interactions between two phenylalanine side chains of SSB (Phe-171 and Phe-177) and an apolar groove on the CTD, combined with electrostatic interactions between the C-terminal carboxylate of SSB and Lys-214 of the CTD. Mutation of any of these residues to alanine significantly disrupts the interaction of full-length Redβ and SSB proteins. Structural knowledge of this interaction will help to expand the utility of Redβ-mediated recombination to a wider range of bacterial hosts for applications in synthetic biology.

Influence of point mutations on PR65 conformational adaptability: Insights from molecular simulations and nanoaperture optical tweezers

PR65 is the HEAT repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem repeat protein. Its conformational mechanics plays a crucial role in PP2A function by opening/closing substrate binding/catalysis interface. Using in silico saturation mutagenesis, we identified PR65 "hinge" residues whose substitutions could alter its conformational adaptability and thereby PP2A function, and selected six mutations that were verified to be expressed and soluble. Molecular simulations and nanoaperture optical tweezers revealed consistent results on the specific effects of the mutations on the structure and dynamics of PR65. Two mutants observed in simulations to stabilize extended/open conformations exhibited higher corner frequencies and lower translational scattering in experiments, indicating a shift toward extended conformations, whereas another displayed the opposite features, confirmed by both simulations and experiments. The study highlights the power of single-molecule nanoaperture-based tweezers integrated with in silico approaches for exploring the effect of mutations on protein structure and dynamics.

Structural conversion of the spidroin C-terminal domain during assembly of spider silk fibers

The major ampullate Spidroin 1 (MaSp1) is the main protein of the dragline spider silk. The C-terminal (CT) domain of MaSp1 is crucial for the self-assembly into fibers but the details of how it contributes to the fiber formation remain unsolved. Here we exploit the fact that the CT domain can form silk-like fibers by itself to gain knowledge about this transition. Structural investigations of fibers from recombinantly produced CT domain from E. australis MaSp1 reveal an α-helix to β-sheet transition upon fiber formation and highlight the helix No4 segment as most likely to initiate the structural conversion. This prediction is corroborated by the finding that a peptide corresponding to helix No4 has the ability of pH-induced conversion into β-sheets and self-assembly into nanofibrils. Our results provide structural information about the CT domain in fiber form and clues about its role in triggering the structural conversion of spidroins during fiber assembly.

Protein target similarity is positive predictor of in vitro antipathogenic activity: a drug repurposing strategy for Plasmodium falciparum

Drug discovery is an intricate and costly process. Repurposing existing drugs and active compounds offers a viable pathway to develop new therapies for various diseases. By leveraging publicly available biomedical information, it is possible to predict compounds' activity and identify their potential targets across diverse organisms. In this study, we aimed to assess the antiplasmodial activity of compounds from the Repurposing, Focused Rescue, and Accelerated Medchem (ReFRAME) library using in vitro and bioinformatics approaches. We assessed the in vitro antiplasmodial activity of the compounds using blood-stage and liver-stage drug susceptibility assays. We used protein sequences of known targets of the ReFRAME compounds with high antiplasmodial activity (EC50 < 10 uM) to conduct a protein-pairwise search to identify similar Plasmodium falciparum 3D7 proteins (from PlasmoDB) using NCBI protein BLAST. We further assessed the association between the compounds' in vitro antiplasmodial activity and level of similarity between their known and predicted P. falciparum target proteins using simple linear regression analyses. BLAST analyses revealed 735 P. falciparum proteins that were similar to the 226 known protein targets associated with the ReFRAME compounds. Antiplasmodial activity of the compounds was positively associated with the degree of similarity between the compounds' known targets and predicted P. falciparum protein targets (percentage identity, E value, and bit score), the number of the predicted P. falciparum targets, and their respective mutagenesis index and fitness scores (R2 between 0.066 and 0.92, P < 0.05). Compounds predicted to target essential P. falciparum proteins or those with a druggability index of 1 showed the highest antiplasmodial activity.

In Silico Characterization of Pathogenic Homeodomain Missense Mutations in the PITX2 Gene

Paired homologous domain transcription factor 2 (PITX2) is critically involved in ocular and cardiac development. Mutations in PITX2 are consistently reported in association with Axenfeld-Rieger syndrome, an autosomal dominant genetic disorder and atrial fibrillation, a common cardiac arrhythmia. In this study, we have mined missense mutations in PITX2 gene from NCBI-dbSNP and Ensembl databases, evaluated the pathogenicity of the missense variants in the homeodomain and C-terminal region using five in silico prediction tools SIFT, PolyPhen2, GERP, Mutation Assessor and CADD. Fifteen homeodomain mutations G42V, G42R, R45W, S49Y, R53W, E53D, E55V, R62H, P65S, R69H, G75R, R84G, R86K, R87W, R91P were found to be highly pathogenic by both SIFT, PolyPhen2 were further functionally characterized using I-Mutant 2.0, Consurf, MutPred and Project Hope. The findings of the study can be used for prioritizing mutations in the context of genetic studies.

An In Silico Analysis of Genetic Variants and Structural Modeling of the Human Frataxin Protein in Friedreich's Ataxia

Friedreich's Ataxia (FRDA) stands out as the most prevalent form of hereditary ataxias, marked by progressive movement ataxia, loss of vibratory sensitivity, and skeletal deformities, severely affecting daily functioning. To date, the only medication available for treating FRDA is Omaveloxolone (Skyclarys®), recently approved by the FDA. Missense mutations within the human frataxin (FXN) gene, responsible for intracellular iron homeostasis regulation, are linked to FRDA development. These mutations induce FXN dysfunction, fostering mitochondrial iron accumulation and heightened oxidative stress, ultimately triggering neuronal cell death pathways. This study amalgamated 226 FXN genetic variants from the literature and database searches, with only 18 previously characterized. Predictive analyses revealed a notable prevalence of detrimental and destabilizing predictions for FXN mutations, predominantly impacting conserved residues crucial for protein function. Additionally, an accurate, comprehensive three-dimensional model of human FXN was constructed, serving as the basis for generating genetic variants I154F and W155R. These variants, selected for their severe clinical implications, underwent molecular dynamics (MD) simulations, unveiling flexibility and essential dynamic alterations in their N-terminal segments, encompassing FXN42, FXN56, and FXN78 domains pivotal for protein maturation. Thus, our findings indicate potential interaction profile disturbances in the FXN42, FXN56, and FXN78 domains induced by I154F and W155R mutations, aligning with the existing literature.

Novel β-galactosidase activity and first crystal structure of Glycoside Hydrolase family 154

Polysaccharide Utilization Loci (PULs) are physically linked gene clusters conserved in the Gram-negative phylum of Bacteroidota and are valuable sources for Carbohydrate Active enZyme (CAZyme) discovery. This study focuses on BD-β-Gal, an enzyme encoded in a metagenomic PUL and member of the Glycoside Hydrolase family 154 (GH154). BD-β-Gal showed exo-β-galactosidase activity with regiopreference for hydrolyzing β-d-(1,6) glycosidic linkages. Notably, it exhibited a preference for d-glucopyranosyl (d-Glcp) over d-galactopyranosyl (d-Galp) and d-fructofuranosyl (d-Fruf) at the reducing end of the investigated disaccharides. In addition, we determined the high resolution crystal structure of BD-β-Gal, thus providing the first structural characterization of a GH154 enzyme. Surprisingly, this revealed an (α/α)6 topology, which has not been observed before for β-galactosidases. BD-β-Gal displayed low structural homology with characterized CAZymes, but conservation analysis suggested that the active site was located in a central cavity, with conserved E73, R252, and D253 as putative catalytic residues. Interestingly, BD-β-Gal has a tetrameric structure and a flexible loop from a neighboring protomer may contribute to its reaction specificity. Finally, we showed that the founding member of GH154, BT3677 from Bacteroides thetaiotaomicron, described as β-glucuronidase, displayed exo-β-galactosidase activity like BD-β-Gal but lacked a tetrameric structure.

Darobactin Substrate Engineering and Computation Show Radical Stability Governs Ether versus C-C Bond Formation

The Gram-negative selective antibiotic darobactin A has attracted interest owing to its intriguing fused bicyclic structure and unique targeting of the outer membrane protein BamA. Darobactin, a ribosomally synthesized and post-translationally modified peptide (RiPP), is produced by a radical S-adenosyl methionine (rSAM)-dependent enzyme (DarE) and contains one ether and one C-C cross-link. Herein, we analyze the substrate tolerance of DarE and describe an underlying catalytic principle of the enzyme. These efforts produced 51 enzymatically modified darobactin variants, revealing that DarE can install the ether and C-C cross-links independently and in different locations on the substrate. Notable variants with fused bicyclic structures were characterized, including darobactin W3Y, with a non-Trp residue at the twice-modified central position, and darobactin K5F, which displays a fused diether ring pattern. While lacking antibiotic activity, quantum mechanical modeling of darobactins W3Y and K5F aided in the elucidation of the requisite features for high-affinity BamA engagement. We also provide experimental evidence for β-oxo modification, which adds support for a proposed DarE mechanism. Based on these results, ether and C-C cross-link formation was investigated computationally, and it was determined that more stable and longer-lived aromatic Cβ radicals correlated with ether formation. Further, molecular docking and transition state structures based on high-level quantum mechanical calculations support the different indole connectivity observed for ether (Trp-C7) and C-C (Trp-C6) cross-links. Finally, mutational analysis and protein structural predictions identified substrate residues that govern engagement to DarE. Our work informs on darobactin scaffold engineering and further unveils the underlying principles of rSAM catalysis.

Point mutation in a virus-like capsid drives symmetry reduction to form tetrahedral cages

Protein capsids are a widespread form of compartmentalization in nature. Icosahedral symmetry is ubiquitous in capsids derived from spherical viruses, as this geometry maximizes the internal volume that can be enclosed within. Despite the strong preference for icosahedral symmetry, we show that simple point mutations in a virus-like capsid can drive the assembly of unique symmetry-reduced structures. Starting with the encapsulin from Myxococcus xanthus, a 180-mer bacterial capsid that adopts the well-studied viral HK97 fold, we use mass photometry and native charge detection mass spectrometry to identify a triple histidine point mutant that forms smaller dimorphic assemblies. Using cryoelectron microscopy, we determine the structures of a precedented 60-mer icosahedral assembly and an unexpected 36-mer tetrahedron that features significant geometric rearrangements around a new interaction surface between capsid protomers. We subsequently find that the tetrahedral assembly can be generated by triple-point mutation to various amino acids and that even a single histidine point mutation is sufficient to form tetrahedra. These findings represent a unique example of tetrahedral geometry when surveying all characterized encapsulins, HK97-like capsids, or indeed any virus-derived capsids reported in the Protein Data Bank, revealing the surprising plasticity of capsid self-assembly that can be accessed through minimal changes in the protein sequence.

The crystal structure of bacteriophage λ RexA provides novel insights into the DNA binding properties of Rex-like phage exclusion proteins

RexA and RexB function as an exclusion system that prevents bacteriophage T4rII mutants from growing on Escherichia coli λ phage lysogens. Recent data established that RexA is a non-specific DNA binding protein that can act independently of RexB to bias the λ bistable switch toward the lytic state, preventing conversion back to lysogeny. The molecular interactions underlying these activities are unknown, owing in part to a dearth of structural information. Here, we present the 2.05-Å crystal structure of the λ RexA dimer, which reveals a two-domain architecture with unexpected structural homology to the recombination-associated protein RdgC. Modelling suggests that our structure adopts a closed conformation and would require significant domain rearrangements to facilitate DNA binding. Mutagenesis coupled with electromobility shift assays, limited proteolysis, and double electron-electron spin resonance spectroscopy support a DNA-dependent conformational change. In vivo phenotypes of RexA mutants suggest that DNA binding is not a strict requirement for phage exclusion but may directly contribute to modulation of the bistable switch. We further demonstrate that RexA homologs from other temperate phages also dimerize and bind DNA in vitro. Collectively, these findings advance our mechanistic understanding of Rex functions and provide new evolutionary insights into different aspects of phage biology.

Biomolecular investigations into BBI reveal an enzymatic mechanism for PUFA isomerisation in bifidobacterium CFA bioconversion strains

Multiple species of Bifidobacterium exhibit the ability to bioconvert conjugated fatty acids (CFAs), which is considered an important pathway for these strains to promote host health. However, there has been limited progress in understanding the enzymatic mechanism of CFA bioconversion by bifidobacteria, despite the increasing number of studies identifying CFA-producing strains. The protein responsible for polyunsaturated fatty acid (PUFA) isomerization in B. breve CCFM683 has recently been discovered and named BBI, providing a starting point for exploring Bifidobacterium isomerases (BIs). This study presents the sequence classification of membrane-bound isomerases from four common Bifidobacterium species that produce CFA. Heterologous expression, purification, and enzymatic studies of the typical sequences revealed that all possess a single c9, t11 isomer as the product and share common features in terms of enzymatic properties and catalytic kinetics. Using molecular docking and alanine scanning, Lys84, Tyr198, Asn202, and Leu245 located in the binding pocket were identified as critical to the catalytic activity, a finding further confirmed by site-directed mutagenesis-based screening assays. Overall, these findings provide insightful knowledge concerning the molecular mechanisms of BIs. This will open up additional opportunities for the use of bifidobacteria and CFAs in probiotic foods and precision nutrition.

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.

Structural flexibility and heterogeneity of recombinant human glial fibrillary acidic protein (GFAP)

Glial fibrillary acidic protein (GFAP) is a promising biomarker for brain and spinal cord disorders. Recent studies have highlighted the differences in the reliability of GFAP measurements in different biological matrices. The reason for these discrepancies is poorly understood as our knowledge of the protein's 3-dimensional conformation, proteoforms, and aggregation remains limited. Here, we investigate the structural properties of GFAP under different conditions. For this, we characterized recombinant GFAP proteins from various suppliers and applied hydrogen-deuterium exchange mass spectrometry (HDX-MS) to provide a snapshot of the conformational dynamics of GFAP in artificial cerebrospinal fluid (aCSF) compared to the phosphate buffer. Our findings indicate that recombinant GFAP exists in various conformational species. Furthermore, we show that GFAP dimers remained intact under denaturing conditions. HDX-MS experiments show an overall decrease in H-bonding and an increase in solvent accessibility of GFAP in aCSF compared to the phosphate buffer, with clear indications of mixed EX2 and EX1 kinetics. To understand possible structural interface regions and the evolutionary conservation profiles, we combined HDX-MS results with the predicted GFAP-dimer structure by AlphaFold-Multimer. We found that deprotected regions with high structural flexibility in aCSF overlap with predicted conserved dimeric 1B and 2B domain interfaces. Structural property predictions combined with the HDX data show an overall deprotection and signatures of aggregation in aCSF. We anticipate that the outcomes of this research will contribute to a deeper understanding of the structural flexibility of GFAP and ultimately shed light on its behavior in different biological matrices.

Evolutionary dynamics of the cytoskeletal profilin gene family in Brassica juncea L. reveal its roles in silique development and stress resilience

Essential to plant adaptation, cell wall (CW) integrity is maintained by CW-biosynthesis genes. Cytoskeletal actin-(de)polymerizing, phospholipid-binding profilin (PRF) proteins play important roles in maintaining cellular homeostasis across kingdoms. However, evolutionary selection of PRF genes and their systematic characterization in family Brassicaceae, especially in Brassica juncea remain unexplored. Here, a comprehensive analysis of genome-wide identification of BjPRFs, their phylogenetic association, genomic localization, gene structure, and transcriptional profiling were performed in an evolutionary framework. Identification of 23 BjPRFs in B. juncea indicated an evolutionary conservation within Brassicaceae. The BjPRFs evolved through paralogous and orthologous gene formation in Brassica genomes. Evolutionary divergence of BjPRFs indicated purifying selection, with nonsynonymous (dN)/synonymous (dS) value of 0.090 for orthologous gene-pairs. Hybrid homology-modeling identified evolutionary distinct and conserved domains in BjPRFs which suggested that these proteins evolved following the divergence of monocot and eudicot plants. RNA-seq profiles of BjPRFs revealed their functional evolution in spatiotemporal manner during plant-development and stress-conditions in diploid/amphidiploid Brassica species. Real-Time PCR experiments in seedling, vegetative, floral and silique tissues of B. juncea suggested their essential roles in systematic plant development. These observations underscore the expansion of BjPRFs in B. juncea, and offer valuable evolutionary insights for exploring cellular mechanisms, and stress resilience.

Computational analysis of sodium-dependent phosphate transporter SLC20A1/PiT1 gene identifies missense variations C573F, and T58A as high-risk deleterious SNPs

SLC20A1/PiT1 is a sodium-dependent inorganic phosphate transporter, initially recognized as the retroviral receptor for Gibbon Ape Leukemia Virus in humans. SNPs in SLC20A1 is associated with Combined Pituitary Hormone Deficiency and Sodium Lithium Counter transport. Using in silico techniques, we have screened the nsSNPs for their deleterious effect on the structure and function of SLC20A1. Screening with sequence and structure-based tools on 430 nsSNPs, filtered 17 nsSNPs which are deleterious. To evaluate the role of these SNPs, protein modeling and MD simulations were performed. A comparative analysis of model generated with SWISS-MODEL and AlphaFold shows that many residues are in the disallowed region of Ramachandran plot. Since SWISS-MODEL structure has a 25-residue deletion, the AlphaFold structure was used to perform MD simulation for equilibration and structure refinement. Further, to understand perturbation of energetics, we performed in silico mutagenesis and ΔΔG calculation using FoldX on MD refined structures, which yielded SNPs that are neutral (3), destabilizing (12) and stabilizing (2) on protein structure. Furthermore, to elucidate the impact of SNPs on structure, we performed MD simulations to discern the changes in RMSD, Rg, RMSF and LigPlot of interacting residues. RMSF profiles of representative SNPs revealed that A114V (neutral) and T58A (positive) were more flexible & C573F (negative) was more rigid compared to wild type, which is also reflected in the changes in number of local interacting residues in LigPlot and ΔΔG. Taken together, our results show that SNPs can lead to structural perturbations and impact the function of SLC20A1 with potential implications for disease.Communicated by Ramaswamy H. Sarma.

Genetic analysis of PALB2 gene WD40 domain in canine mammary tumour patients

BACKGROUND

DNA repair mechanisms are essential for tumorigenesis and disruption of HR mechanism is an important predisposing factor of human breast cancers (BC). PALB2 is an important part of the HR. There are similarities between canine mammary tumours (CMT) and BCs. As its human counterpart, PALB2 mutations could be a predisposing factor of CMT.

OBJECTIVES

In this study, we aimed to investigate the impacts of PALB2 variants on tumorigenesis and canine mammary tumor (CMT) malignancy.

METHODS

We performed Sanger sequencing to detect germline mutations in the WD40 domain of the canine PALB2 gene in CMT patients. We conducted in silico analysis to investigate the variants, and compared the germline PALB2 mutations in humans that cause breast cancer (BC) with the variants detected in dogs with CMT.

RESULTS

We identified an intronic (c.3096+8C>G) variant, two exonic (p.A1050V and p.R1354R) variants, and a 3' UTR variant (c.4071T>C). Of these, p.R1354R and c.4071T>C novel variants were identified for the first time in this study. We found that the p.A1050V mutation had a significant effect. However, we could not determine sufficient similarity due to the differences in nucleotide/amino acid sequences between two species. Nonetheless, possible variants of human sequences in the exact location as their dog counterparts are associated with several cancer types, implying that the variants could be crucial for tumorigenesis in dogs. Our results did not show any effect of the variants on tumor malignancy.

CONCLUSIONS

The current project is the first study investigating the relationship between the PALB2 gene WD40 domain and CMTs. Our findings will contribute to a better understanding of the pathogenic mechanism of the PALB2 gene in CMTs. In humans, variant positions in canines have been linked to cancer-related phenotypes such as familial BC, endometrial tumor, and hereditary cancer predisposition syndrome. The results of bioinformatics analyses should be investigated through functional tests or case-control studies.

Simulation and Computational Study of RING Domain Mutants of BRCA1 and Ube2k in AD/PD Pathophysiology

Lysine-based post-translational modification (PTM) such as acylation, acetylation, deamination, methylation, SUMOylation, and ubiquitination has proven to be a major regulator of gene expression, chromatin structure, protein stability, protein-protein interaction, protein degradation, and cellular localization. However, besides all the PTMs, ubiquitination stands as the second most common PTM after phosphorylation that is involved in the etiology of neurodegenerative diseases (NDDs) namely, Alzheimer's disease (AD) and Parkinson's disease (PD). NDDs are characterized by the accumulation of misfolded protein aggregates in the brain that lead to disease-related gene mutation and irregular protein homeostasis. The ubiquitin-proteasome system (UPS) is in charge of degrading these misfolded proteins, which involve an interplay of E1, E2, E3, and deubiquitinase enzymes. Impaired UPS has been commonly observed in NDDs and E3 ligases are the key members of the UPS, thus, dysfunction of the same can accelerate the neurodegeneration process. Therefore, the aim of this study is firstly, to find E3 ligases that are common in both AD and PD through data mining. Secondly, to study the impact of mutation on its structure and function. The study deciphered 74 E3 ligases that were common in both AD and PD. Later, 10 hub genes were calculated of which protein-protein interaction, pathway enrichment, lysine site prediction, domain, and motif analysis were performed. The results predicted BRCA1, PML, and TRIM33 as the top three putative lysine-modified E3 ligases involved in AD and PD pathogenesis. However, based on structural characterization, BRCA1 was taken further to study RING domain mutation that inferred K32Y, K32L, K32C, K45V, K45Y, and K45G as potential mutants that alter the structural and functional ability of BRCA1 to interact with Ube2k, E2-conjugating enzyme. The most probable mutant observed after molecular dynamics simulation of 50 ns is K32L. Therefore, our study concludes BRCA1, a potential E3 ligase common in AD and PD, and RING domain mutation at sites K32 and K45 possibly disturbs its interaction with its E2, Ube2k.

Unravelling the Relacatib activity against the CTSK proteins causing pycnodysostosis: a molecular docking and dynamics approach

Pycnodysostosis is an atypical autosomal recessive condition of Lysosomal storage disorder that originated due to the deficit of the enzyme Cathepsin K which is vital for normal osteoclast action in bone resorption. Abnormal degradation of type 1 collagen and accumulation of toxic undigested collagen fibers in lysosomes of the osteoclast cells resulting in high bone density, brittle bones, and a short stature is caused in CTSK protein-carrying individuals. The broad aim of this study is to identify the most significant variant through various computational pipelines. This study was initiated by retrieving a total number of thirty-six variants from NCBI, HGMD, and UniProt databases, and the Y283C variant was found to be more significant by various standard computational tools. A structural investigation was performed to understand and gain a better knowledge about the interaction profile for the native (1BY8) and variant (Y283C) with Relacatib (a small-molecule drug that blocks the function of Cathepsin K, an enzyme that has been linked to osteoporosis, osteoarthritis, and other bone-degrading diseases). The interaction profile was analyzed using molecular docking. Relacatib (ligand) had an average binding affinity for both native (-7.16 kcal/mol) and Y283C (-6.76 kcal/mol). Finally, Molecular dynamics simulations were done in duplicates to recognize the variant (Y283C) activity of the protein structure against Relacatib for 100 ns. This study assists in comprehending the most pathogenic amino-acid variant, the ligand interaction with the protein structure, and paves the way for understanding the steadiness of the ligand with the native and selected significant amino-acid variant.Communicated by Ramaswamy H. Sarma.

M. Ahmed S, Alsafar HS, Halwani R. Multiple inborn errors of type I IFN immunity in a 33-year-old male with a fatal case of COVID-19

The host genetic inborn errors of immunity (IEIs) have been shown to contribute to susceptibility to life-threatening coronavirus disease 2019 (COVID-19), as it had been associated previously with other viral infections. Most genetic association studies have described IEIs as a monogenic defect, while there have been no reports of patients with multiple inherited immune deficiencies. This is a complex case of IEIs predisposing to severe viral infections in an unvaccinated 33-year-old male patient. The patient was admitted with no respiratory symptoms, showed a SARS-CoV-2 PCR positive test on the second day of admission, started developing progressive lung consolidation within three days of hospitalization, and was moved from non-invasive to mechanical ventilation within 12 days of hospitalization. Impaired production of type I IFN was detected in patient PBMCs treated with poly(I:C), at both mRNA and protein levels. Whole exome sequencing revealed three mutations across type I IFN production pathway, which were predicted to be loss-of-function (pLOF). The three mutations were predicted to predispose to severe viral infections: monoallelic R488X TLR3, monoallelic His684Arg TLR3, and biallelic Val363Met IRF3. Functional analysis confirmed that all these mutations dysregulated the type I IFN pathway. Evaluation of TLR3 and IRF3 IFN-β1 luciferase reporter activity showed a hypomorphic suppression of function. TOPO TA cloning was used to ascertain the positioning of both TLR3 variants, indicating that both variants were on the same allele. We have described a unique complex IEI patient with multiple mutations, particularly along type I IFN production pathway.

High-throughput crystallographic fragment screening of Zika virus NS3 Helicase

The Zika virus (ZIKV), discovered in Africa in 1947, swiftly spread across continents, causing significant concern due to its recent association with microcephaly in newborns and Guillain-Barré syndrome in adults. Despite a decrease in prevalence, the potential for a resurgence remains, necessitating urgent therapeutic interventions. Like other flaviviruses, ZIKV presents promising drug targets within its replication machinery, notably the NS3 helicase (NS3Hel) protein, which plays critical roles in viral replication. However, a lack of structural information impedes the development of specific inhibitors targeting NS3Hel. Here we applied high-throughput crystallographic fragment screening on ZIKV NS3Hel, which yielded structures that reveal 3D binding poses of 46 fragments at multiple sites of the protein, including 11 unique fragments in the RNA-cleft site. These fragment structures provide templates for direct design of hit compounds and should thus assist the development of novel direct-acting antivirals against ZIKV and related flaviviruses, thus opening a promising avenue for combating future outbreaks.

The D Gene in CDR H3 Determines a Public Class of Human Antibodies to SARS-CoV-2

Public antibody responses have been found against many infectious agents. Structural convergence of public antibodies is usually determined by immunoglobulin V genes. Recently, a human antibody public class against SARS-CoV-2 was reported, where the D gene (IGHD3-22) encodes a common YYDxxG motif in heavy-chain complementarity-determining region 3 (CDR H3), which determines specificity for the receptor-binding domain (RBD). In this review, we discuss the isolation, structural characterization, and genetic analyses of this class of antibodies, which have been isolated from various cohorts of COVID-19 convalescents and vaccinees. All eleven YYDxxG antibodies with available structures target the SARS-CoV-2 RBD in a similar binding mode, where the CDR H3 dominates the interaction with antigen. The antibodies target a conserved site on the RBD that does not overlap with the receptor-binding site, but their particular angle of approach results in direct steric hindrance to receptor binding, which enables both neutralization potency and breadth. We also review the properties of CDR H3-dominant antibodies that target other human viruses. Overall, unlike most public antibodies, which are identified by their V gene usage, this newly discovered public class of YYDxxG antibodies is dominated by a D-gene-encoded motif and uncovers further opportunities for germline-targeting vaccine design.

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.

Evolutionary trajectory of TRPM2 channel activation by adenosine diphosphate ribose and calcium

Ion channel activation upon ligand gating triggers a myriad of biological events and, therefore, evolution of ligand gating mechanism is of fundamental importance. TRPM2, a typical ancient ion channel, is activated by adenosine diphosphate ribose (ADPR) and calcium and its activation has evolved from a simple mode in invertebrates to a more complex one in vertebrates, but the evolutionary process is still unknown. Molecular evolutionary analysis of TRPM2s from more than 280 different animal species has revealed that, the C-terminal NUDT9-H domain has evolved from an enzyme to a ligand binding site for activation, while the N-terminal MHR domain maintains a conserved ligand binding site. Calcium gating pattern has also evolved, from one Ca2+-binding site as in sea anemones to three sites as in human. Importantly, we identified a new group represented by olTRPM2, which has a novel gating mode and fills the missing link of the channel gating evolution. We conclude that the TRPM2 ligand binding or activation mode evolved through at least three identifiable stages in the past billion years from simple to complicated and coordinated. Such findings benefit the evolutionary investigations of other channels and proteins.

Structural basis of ligand specificity and channel activation in an insect gustatory receptor

Gustatory receptors (GRs) are critical for insect chemosensation and are potential targets for controlling pests and disease vectors, making their structural investigation a vital step toward such applications. We present structures of Bombyx mori Gr9 (BmGr9), a fructose-gated cation channel, in agonist-free and fructose-bound states. BmGr9 forms a tetramer similar to distantly related insect odorant receptors (ORs). Upon fructose binding, BmGr9's channel gate opens through helix S7b movements. In contrast to ORs, BmGr9's ligand-binding pocket, shaped by a kinked helix S4 and a shorter extracellular S3-S4 loop, is larger and solvent accessible in both agonist-free and fructose-bound states. Also, unlike ORs, fructose binding by BmGr9 involves helix S5 and a pocket lined with aromatic and polar residues. Structure-based sequence alignments reveal distinct patterns of ligand-binding pocket residue conservation in GR subfamilies associated with different ligand classes. These data provide insight into the molecular basis of GR ligand specificity and function.

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.