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

Prediction of deleterious non-synonymous SNPs of human MDC1 gene: an in silico approach

MDC1 (Mediator of DNA damage Checkpoint protein 1) functions to facilitate the localization of numerous DNA damage response (DDR) components to DNA double-strand break sites. MDC1 is an integral component in preserving genomic stability and appropriate DDR regulation. There haven't been systematic investigations of MDC1 mutations that induce cancer and genomic instability. Variations in nsSNPs have the potential to modify the protein chemistry and their function. Describing functional SNPs in disease-associated genes presents a significant conundrum for investigators, it is possible to assess potential functional SNPs before conducting larger population examinations. Multiple sequences and structure-based bioinformatics strategies were implemented in the current in-silico investigation to discern potential nsSNPs of the MDC1 genes. The nsSNPs were identified with SIFT, SNAP2, Align GVGD, PolyPhen-2, and PANTHER, and their stability was determined with MUpro. The conservation, solvent accessibility, and structural effects of the mutations were identified with ConSurf, NetSurfP-2.0, and SAAFEC-SEQ respectively. Cancer-related analysis of the nsSNPs was conducted using cBioPortal and TCGA web servers. The present study appraised five nsSNPs (P1426T, P69S, P194R, P203L, and H131Y) as probably mutilating due to their existence in highly conserved regions and propensity to deplete protein stability. The nsSNPs P194R, P203L, and H131Y were concluded as deleterious and possibly damaging from the 5 prediction tools. The functional nsSNP P194R mutation is associated with skin cutaneous melanoma while no significant records were found for other nsSNPs. The present study concludes that the highly deleterious P194R mutations can potentially induce genomic instability and contribute to various cancers' pathogenesis. Developing drugs targeting these mutations can undoubtedly be advantageous in large population-based studies, particularly in the development of precision medicine.

Structure of the human systemic RNAi defective transmembrane protein 1 (hSIDT1) reveals the conformational flexibility of its lipid binding domain

In Caenorhabditis elegans, inter-cellular transport of the small non-coding RNA causing systemic RNAi is mediated by the transmembrane protein SID1, encoded by the sid1 gene in the systemic RNAi defective (sid) loci. SID1 shares structural and sequence similarity with cholesterol uptake protein 1 (CHUP1) and is classified as a member of the ChUP family. Although systemic RNAi is not an evolutionarily conserved process, the sid gene products are found across the animal kingdom, suggesting the existence of other novel gene regulatory mechanisms mediated by small non-coding RNAs. Human homologs of sid gene products-hSIDT1 and hSIDT2-mediate contact-dependent lipophilic small non-coding dsRNA transport. Here, we report the structure of recombinant human SIDT1. We find that the extra-cytosolic domain of hSIDT1 adopts a double jelly roll fold, and the transmembrane domain exists as two modules-a flexible lipid binding domain and a rigid transmembrane domain core. Our structural analyses provide insights into the inherent conformational dynamics within the lipid binding domain in ChUP family members.

Positive selection analysis of Cyt proteins from Bacillus thuringiensis: A conservative trend driven by negative (purifying) selection

Bacillus thuringiensis is a Gram-positive entomopathogenic bacterium that produces different pesticidal proteins: vegetative insecticidal proteins (Vpb1/Vpa2, Vip3, and Vpb4) during vegetative growth, which are secreted to the culture medium, and δ-endotoxins (Cry and Cyt) during sporulation, which accumulate into parasporal crystals. Cyt proteins are the smaller subset of δ-endotoxins targeting Diptera species. While Cry and Vip3 proteins undergo positive selection, our analysis suggests that Cyt proteins evolve following a conservative trend driven negative (purifying) selection.

Structure-guided Evolutionary Analysis of Interactome Network Rewiring at Single Residue Resolution in Yeasts

Protein-protein interactions (PPIs) are known to rewire extensively during evolution leading to lineage-specific and species-specific changes in molecular processes. However, the detailed molecular evolutionary mechanisms underlying interactome network rewiring are not well-understood. Here, we combine high-confidence PPI data, high-resolution three-dimensional structures of protein complexes, and homology-based structural annotation transfer to construct structurally-resolved interactome networks for the two yeasts S. cerevisiae and S. pombe. We then classify PPIs according to whether they are preserved or different between the two yeast species and compare site-specific evolutionary rates of interfacial versus non-interfacial residues for these different categories of PPIs. We find that residues in PPI interfaces evolve significantly more slowly than non-interfacial residues when using lineage-specific measures of evolutionary rate, but not when using non-lineage-specific measures. Furthermore, both lineage-specific and non-lineage-specific evolutionary rate measures can distinguish interfacial residues from non-interfacial residues for preserved PPIs between the two yeasts, but only the lineage-specific measure is appropriate for rewired PPIs. Finally, both lineage-specific and non-lineage-specific evolutionary rate measures are appropriate for elucidating structural determinants of protein evolution for residues outside of PPI interfaces. Overall, our results demonstrate that unlike tertiary structures of single proteins, PPIs and PPI interfaces can be highly volatile in their evolution, thus requiring the use of lineage-specific measures when studying their evolution. These results yield insight into the evolutionary design principles of PPIs and the mechanisms by which interactions are preserved or rewired between species, improving our understanding of the molecular evolution of PPIs and PPI interfaces at the residue level.

PCAF-mediated acetylation of METTL3 impairs mRNA translation efficiency in response to oxidative stress

METTL3 methylates RNA and regulates the fate of mRNA through its methyltransferase activity. METTL3 enhances RNA translation independently of its catalytic activity. However, the underlying mechanism is still elusive. Here, we report that METTL3 is both interacted with and acetylated at lysine 177 by the acetyltransferase PCAF and deacetylated by SIRT3. Neither the methyltransferase activity nor the stability of METTL3 is affected by its acetylation at K177. Importantly, acetylation of METTL3 blocks its interaction with EIF3H, a subunit of the translation initiation factor, thereby reducing mRNA translation efficiency. Interestingly, acetylation of METTL3 responds to oxidative stress. Mechanistically, oxidative stress enhances the interaction of PCAF with METTL3, increases METTL3 acetylation, and suppresses the interaction of METTL3 with EIF3H, thereby decreasing the translation efficiency of ribosomes and inhibiting cell proliferation. Altogether, we suggest a mechanism by which oxidative stress regulates RNA translation efficiency by the modulation of METTL3 acetylation mediated by PCAF.

MAGPIE: An interactive tool for visualizing and analyzing protein-ligand interactions

Quantitative tools to compile and analyze biomolecular interactions among chemically diverse binding partners would improve therapeutic design and aid in studying molecular evolution. Here we present Mapping Areas of Genetic Parsimony In Epitopes (MAGPIE), a publicly available software package for simultaneously visualizing and analyzing thousands of interactions between a single protein or small molecule ligand (the "target") and all of its protein binding partners ("binders"). MAGPIE generates an interactive three-dimensional visualization from a set of protein complex structures that share the target ligand, as well as sequence logo-style amino acid frequency graphs that show all the amino acids from the set of protein binders that interact with user-defined target ligand positions or chemical groups. MAGPIE highlights all the salt bridge and hydrogen bond interactions made by the target in the visualization and as separate amino acid frequency graphs. Finally, MAGPIE collates the most common target-binder interactions as a list of "hotspots," which can be used to analyze trends or guide the de novo design of protein binders. As an example of the utility of the program, we used MAGPIE to probe how different antibody fragments bind a viral antigen; how a common metabolite binds diverse protein partners; and how two ligands bind orthologs of a well-conserved glycolytic enzyme for a detailed understanding of evolutionarily conserved interactions involved in its activation and inhibition. MAGPIE is implemented in Python 3 and freely available at https://github.com/glasgowlab/MAGPIE, along with sample datasets, usage examples, and helper scripts to prepare input structures.

MTR3D-AF2: Expanding the coverage of spatially derived missense tolerance scores across the human proteome using AlphaFold2

The missense tolerance ratio (MTR) was developed as a novel approach to assess the deleteriousness of variants. Its three-dimensional successor, MTR3D, was demonstrated powerful at discriminating pathogenic from benign variants. However, its reliance on experimental structures and homologs limited its coverage of the proteome. We have now utilized AlphaFold2 models to develop MTR3D-AF2, which covers 89.31% of proteins and 85.39% of residues across the human proteome. This work has improved MTR3D's ability to distinguish clinically established pathogenic from benign variants. MTR3D-AF2 is freely available as an interactive web server at https://biosig.lab.uq.edu.au/mtr3daf2/.

Structures of BlEst2 from Bacillus licheniformis in its propeptide and mature forms reveal autoinhibitory effects of the C-terminal domain

Carboxylesterases comprise a major class of α/β-fold hydrolases responsible for the cleavage and formation of ester bonds. Found ubiquitously in nature, these enzymes are crucial for the metabolism of both endogenous and exogenous carboxyl esters in animals, plants and microorganisms. Beyond their essential physiological roles, carboxylesterases stand out as one of the important classes of biocatalysts for biotechnology. BlEst2, an enzyme previously classified as Bacillus licheniformis esterase, remains largely uncharacterized. In the present study, we elucidate the structural biology, molecular dynamics and biochemical features of BlEst2. Our findings reveal a canonical α/β-hydrolase fold similar to the ESTHER block L of lipases, further augmented by two additional accessory C-terminal domains. Notably, the catalytic domain demonstrates two insertions, which occupy conserved locations in α/β-hydrolase proteins and commonly form the lid domain in lipase structures. Intriguingly, our in vitro cleavage of C-terminal domains revealed the structure of the active form of BlEst2. Upon activation, BlEst2 showed a markedly elevated hydrolytic activity. This observation implies that the intramolecular C-terminal domain serves as a regulatory intramolecular inhibitor. Interestingly, despite exhibiting esterase-like activity, BlEst2 structural characteristics align more closely with lipases. This suggests that BlEst2 could potentially represent a previously unrecognized subgroup within the realm of carboxyl ester hydrolases.

Evolution and stress response potential of the plant splicing factor U1C

Alternative splicing is a crucial process in multicellular eukaryote, facilitated by the assembly of spliceosomal complexes comprising numerous small ribonucleoproteins. At an early stage, U1C is thought to be required for 5' splice site recognition and base pairing. However, a systematic analysis of the U1C gene family in response to developmental cues and stress conditions has not yet been conducted in plants. This study identified 114 U1C genes in 72 plant species using basic bioinformatics analyses. Phylogenetic analysis was used to compare gene and protein structures, promoter motifs, and tissue- and stress-specific expression levels, revealing their functional commonalities or diversity in response to developmental cues, such as embryonic expression, or stress treatments, including drought and heat. Fluorescence quantitative expression analysis showed that U1C gene expression changed under salt, low temperature, drought, and Cd stress in rice seedlings. However, gene expression in shoots and roots was not consistent under different stress conditions, suggesting a complex regulatory mechanism. This research provides foundational insights into the U1C gene family's role in plant development and stress responses, highlighting potential targets for future studies.

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.

Secondary interactions in ubiquitin-binding domains achieve linkage or substrate specificity

Small ubiquitin-binding domains (UBDs) recognize small surface patches on ubiquitin with weak affinity, and it remains a conundrum how specific cellular responses may be achieved. Npl4-type zinc-finger (NZF) domains are ∼30 amino acid, compact UBDs that can provide two ubiquitin-binding interfaces, imposing linkage specificity to explain signaling outcomes. We here comprehensively characterize the linkage preference of human NZF domains. TAB2 prefers Lys6 and Lys63 linkages phosphorylated on Ser65, explaining why TAB2 recognizes depolarized mitochondria. Surprisingly, most NZF domains do not display chain linkage preference, despite conserved, secondary interaction surfaces. This suggests that some NZF domains may specifically bind ubiquitinated substrates by simultaneously recognizing substrate and an attached ubiquitin. We show biochemically and structurally that the NZF1 domain of the E3 ligase HOIPbinds preferentially to site-specifically ubiquitinated forms of NEMO and optineurin. Thus, despite their small size, UBDs may impose signaling specificity via multivalent interactions with ubiquitinated substrates.

YjgA plays dual roles in enhancing PTC maturation

Ribosome biogenesis is a highly regulated cellular process that involves the control of numerous assembly factors. The small protein YjgA has been reported to play a role in the late stages of 50S assembly. However, the precise molecular mechanism underlying its function remains unclear. In this study, cryo-electron microscopy (cryo-EM) structures revealed that depletion of YjgA or its N-terminal loop in Escherichia coli both lead to the accumulation of immature 50S particles with structural abnormalities mainly in peptidyl transferase center (PTC) and H68/69 region. CryoDRGN analysis uncovered 8 and 6 distinct conformations of pre50S for ΔyjgA and YjgA-ΔNloop, respectively. These conformations highlighted the role of the N-terminal loop of YjgA in integrating uL16 and stabilizing H89 in PTC, which was further verified by the pull-down assays of YjgA and its mutants with uL16. Together with the function of undocking H68 through the binding of its C-terminal CTLH-like domain to the base of the L1 stalk, YjgA facilitates the maturation of PTC. This study identified critical domains of YjgA contributing to 50S assembly efficiency, providing a comprehensive understanding of the dual roles of YjgA in accelerating ribosome biogenesis and expanding our knowledge of the intricate processes governing cellular protein synthesis.

Understanding P-Rex regulation: structural breakthroughs and emerging perspectives

Rho GTPases are a family of highly conserved G proteins that regulate numerous cellular processes, including cytoskeleton organisation, migration, and proliferation. The 20 canonical Rho GTPases are regulated by ∼85 guanine nucleotide exchange factors (GEFs), with the largest family being the 71 Diffuse B-cell Lymphoma (Dbl) GEFs. Dbl GEFs promote GTPase activity through the highly conserved Dbl homology domain. The specificity of GEF activity, and consequently GTPase activity, lies in the regulation and structures of the GEFs themselves. Dbl GEFs contain various accessory domains that regulate GEF activity by controlling subcellular localisation, protein interactions, and often autoinhibition. This review focuses on the two phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P3)-dependent Rac exchangers (P-Rex), particularly the structural basis of P-Rex1 autoinhibition and synergistic activation. First, we discuss structures that highlight the conservation of P-Rex catalytic and phosphoinositide binding activities. We then explore recent breakthroughs in uncovering the structural basis for P-Rex1 autoinhibition and detail the proposed minimal two-step model of how PI(3,4,5)P3 and Gβγ synergistically activate P-Rex1 at the membrane. Additionally, we discuss the further layers of P-Rex regulation provided by phosphorylation and P-Rex2-PTEN coinhibitory complex formation, although these mechanisms remain incompletely understood. Finally, we leverage the available data to infer how cancer-associated mutations in P-Rex2 destabilise autoinhibition and evade PTEN coinhibitory complex formation, leading to increased P-Rex2 GEF activity and driving cancer progression and metastasis.

Systematic computational strategies for identifying protein targets and lead discovery

Computational algorithms and tools have retrenched the drug discovery and development timeline. The applicability of computational approaches has gained immense relevance owing to the dramatic surge in the structural information of biomacromolecules and their heteromolecular complexes. Computational methods are now extensively used in identifying new protein targets, druggability assessment, pharmacophore mapping, molecular docking, the virtual screening of lead molecules, bioactivity prediction, molecular dynamics of protein-ligand complexes, affinity prediction, and for designing better ligands. Herein, we provide an overview of salient components of recently reported computational drug-discovery workflows that includes algorithms, tools, and databases for protein target identification and optimized ligand selection.

Molecular characterization of Rft1, an ER membrane protein associated with congenital disorder of glycosylation RFT1-CDG

The oligosaccharide needed for protein N-glycosylation is assembled on a lipid carrier via a multi-step pathway. Synthesis is initiated on the cytoplasmic face of the endoplasmic reticulum (ER) and completed on the luminal side after transbilayer translocation of a heptasaccharide lipid intermediate. More than 30 Congenital Disorders of Glycosylation (CDGs) are associated with this pathway, including RFT1-CDG which results from defects in the membrane protein Rft1. Rft1 is essential for the viability of yeast and mammalian cells and was proposed as the transporter needed to flip the heptasaccharide lipid intermediate across the ER membrane. However, other studies indicated that Rft1 is not required for heptasaccharide lipid flipping in microsomes or unilamellar vesicles reconstituted with ER membrane proteins, nor is it required for the viability of at least one eukaryote. It is therefore not known what essential role Rft1 plays in N-glycosylation. Here, we present a molecular characterization of human Rft1, using yeast cells as a reporter system. We show that it is a multi-spanning membrane protein located in the ER, with its N and C-termini facing the cytoplasm. It is not N-glycosylated. The majority of RFT1-CDG mutations map to highly conserved regions of the protein. We identify key residues that are important for Rft1's ability to support N-glycosylation and cell viability. Our results provide a necessary platform for future work on this enigmatic protein.

Revisiting the Roles of Catalytic Residues in Human Ornithine Transcarbamylase

Human ornithine transcarbamylase (hOTC) is a mitochondrial transferase protein involved in the urea cycle and is crucial for the conversion of toxic ammonia to urea. Structural analysis coupled with kinetic studies of Escherichia coli, rat, bovine, and other transferase proteins has identified residues that play key roles in substrate recognition and conformational changes but has not provided direct evidence for all of the active residues involved in OTC function. Here, computational methods were used to predict the likely active residues of hOTC; the function of these residues was then probed with site-directed mutagenesis and biochemical characterization. This process identified previously reported active residues, as well as distal residues that contribute to activity. Mutation of active site residue D263 resulted in a substantial loss of activity without a decrease in protein stability, suggesting a key catalytic role for this residue. Mutation of predicted second-layer residues H302, K307, and E310 resulted in significant decreases in enzymatic activity relative to that of wild-type (WT) hOTC with respect to l-ornithine. The mutation of fourth-layer residue H107 to produce the hOTC H107N variant resulted in a 66-fold decrease in catalytic efficiency relative to that of WT hOTC with respect to carbamoyl phosphate and a substantial loss of thermal stability. Further investigation identified H107 and to a lesser extent E98Q as key residues involved in maintaining the hOTC quaternary structure. This work biochemically demonstrates the importance of D263 in hOTC catalytic activity and shows that residues remote from the active site also play key roles in activity.

Structure analysis of the telomere resolvase from the Lyme disease spirochete Borrelia garinii reveals functional divergence of its C-terminal domain

Borrelia spirochetes are the causative agents of Lyme disease and relapsing fever, two of the most common tick-borne illnesses. A characteristic feature of these spirochetes is their highly segmented genomes which consists of a linear chromosome and a mixture of up to approximately 24 linear and circular extrachromosomal plasmids. The complexity of this genomic arrangement requires multiple strategies for efficient replication and partitioning during cell division, including the generation of hairpin ends found on linear replicons mediated by the essential enzyme ResT, a telomere resolvase. Using an integrative structural biology approach employing advanced modelling, circular dichroism, X-ray crystallography and small-angle X-ray scattering, we have generated high resolution structural data on ResT from B. garinii. Our data provides the first high-resolution structures of ResT from Borrelia spirochetes and revealed active site positioning in the catalytic domain. We also demonstrate that the C-terminal domain of ResT is required for both transesterification steps of telomere resolution, and is a requirement for DNA binding, distinguishing ResT from other telomere resolvases from phage and bacteria. These results advance our understanding of the molecular function of this essential enzyme involved in genome maintenance in Borrelia pathogens.

Exploring the antiviral inhibitory activity of Niloticin against the NS2B/NS3 protease of Dengue virus (DENV2)

Dengue virus (DENV2) is the cause of dengue disease and a worldwide health problem. DENV2 replicates in the host cell using polyproteins such as NS3 protease in conjugation with NS2B cofactor, making NS3 protease a promising antiviral drug-target. This study investigated the efficacy of 'Niloticin' against NS2B/NS3-protease. In silico and in vitro analyses were performed which included interaction of niloticin with NS2B/NS3-protease, protein stability and flexibility, mutation effect, betweenness centrality of residues and analysis of cytotoxicity, protein expression and WNV NS3-protease activity. Similar like acyclovir, niloticin forms strong H-bonds and hydrophobic interactions with residues LEU149, ASN152, LYS74, GLY148 and ALA164. The stability of the niloticin-NS2B/NS3-protease complex was found to be stable compared to the apo NS2B/NS3-protease in structural deviation, PCA, compactness and FEL analysis. The IC50 value of niloticin was 0.14 μM in BHK cells based on in vitro cytotoxicity analysis and showed significant activity at 2.5 μM in a concentration-dependent manner. Western blotting and qRT-PCR analyses showed that niloticin reduced DENV2 protein transcription in a dose-dependent manner. Besides, niloticin confirmed the inhibition of NS3-protease by the SensoLyte 440 WNV protease detection kit. These promising results suggest that niloticin could be an effective antiviral drug against DENV2 and other flaviviruses.

Defects in Glabrous 3 (GL3) functionality underlie the absence of trichomes in Brassica napus

Previously, expression of the Arabidopsis thaliana GLABRA3 (GL3) induced trichome formation in Brassica napus. GL3 orthologues were examined from glabrous (B. oleracea), semi-glabrous (B. napus), moderately hirsute (B. rapa), and very hirsute (B. villosa) Brassica species. Ectopic expression of BnGL3, BrGL3 alleles, or BvGL3 induced trichome formation in glabrous B. napus with the effect on trichome number commensurate with density in the original accessions. Chimeric GL3 proteins in which the B. napus amino terminal region, which interacts with MYB proteins, or the middle region, which interacts with the WD40 protein TTG1, was exchanged with corresponding regions from A. thaliana were as stimulatory to trichome production as AtGL3. Exchange of the carboxy-terminal region containing a bHLH domain and an ACT domain did not alter the trichome stimulatory activity, although modeling of the ACT domain identified differences that could affect GL3 dimerization. B. napus A- and C-genomes orthologues differed in their abilities to form homo- and heterodimers. Modeling of the amino-terminal region revealed a conserved domain that may represent the MYB factor binding pocket. This region interacted with the MYB factors GL1, CPC, and TRY, as well as with JAZ8, which is involved in jasmonic acid-mediated regulation of MYC-like transcription factors. Protein interaction studies indicated that GL1 interaction with GL3 from B. napus and A. thaliana may underlie the difference in their respective abilities to induce trichome formation.

Cefiderocol susceptibility endows hope in treating carbapenem-resistant Pseudomonas aeruginosa: insights from in vitro and in silico evidence

'High-risk' hypermutable clones of Pseudomonas aeruginosa disseminating extensive drug-resistance (XDR) have raised global health concerns with escalating mortality rates in immunocompromised patients. Mutations in conventional drug-targets under antibiotic stress necessitate structural understanding to formulate sustainable therapeutics. In the present study, the major β-lactam antibiotic target, penicillin-binding protein-3 (PBP3) with mutations F533L and T91A, were identified in carbapenemase-positive P. aeruginosa isolates (n = 6) using whole genome sequencing. Antibiotic susceptibility tests showed susceptibility to cefiderocol (MIC ≤ 4 μg ml-1) despite pan-β-lactam resistance in the isolates. Both the mutations reduced local intra-chain interactions in PBP3 that marginally increased the local flexibility (∼1%) in the structures to affect antibiotic-interactions. Molecular dynamics simulations confirmed the overall stability of the PBP3 mutants through root-mean square deviations, radius of gyration, solvent-accessibility and density curves, which favored their selection. Docking studies unveiled that the mutations in PBP3 elicited unfavorable stereochemical clashes with the conventional antibiotics thereby increasing their inhibition constants (IC) up to ∼50 fold. It was deciphered that cefiderocol retained its susceptibility despite mutations in PBP3, due to its higher average binding affinity (ΔG: -8.2 ± 0.4 kcal mol-1) towards multiple PBP-targets and lower average binding affinity (ΔG: -6.7 ± 0.7 kcal mol-1) to β-lactamases than the other β-lactam antibiotics. The molecular dynamics simulations and molecular mechanics Poisson Boltzmann surface area calculations further indicated energetically favorable binding for cefiderocol with PBP3 proteins. The study gave structural insight into emerging non-polar amino acid substitutions in PBP3 causing XDR and recommends prioritizing available antibiotics based on multi-target affinities to overcome challenges imposed by target-protein mutations.

Evolution of orphan and atypical histidine kinases and response regulators for microbial signaling diversity

Two-component signaling systems (TCS) are the predominant means of microbes for sensing and responding to environmental stimuli. Typically, TCS is comprised of a sensor histidine kinase (HK) and a cognate response regulator (RR), which might have coevolved together. They usually involve the phosphoryl transfer signaling mechanism. However, there are also some orphan and atypical HK and RR homologs, and their evolutionary origins are still not very clear. They are not associated with cognate pairs or lack the conserved residues for phosphoryl transfer, but they could receive or respond to signals from other regulators. The objective of this study is to reveal the evolutionary history of these orphan and atypical HK and RR homologs. Structural, domain, sequence, and phylogenetic analyses indicated that their evolution process might undergo gene duplication, divergence, and domain shuffling. Meanwhile, lateral gene transfer might also be involved for their gene distribution. Evolution of orphan and atypical HK and RR homologs have increased their signaling diversity, which could be helpful for microbial adaption in complex environments.

SigmaR1 shapes rough endoplasmic reticulum membrane sheets

Rough endoplasmic reticulum (ER) sheets are a fundamental domain of the ER and the gateway into the secretory pathway. Although reticulon proteins stabilize high-curvature ER tubules, it is unclear whether other proteins scaffold the flat membranes of rough ER sheets. Through a proteomics screen using ER sheet-localized RNA-binding proteins as bait, we identify the sigma-1 receptor (SigmaR1) as an ER sheet-shaping factor. High-resolution live cell imaging and electron tomography assign SigmaR1 as an ER sheet-localized factor whose levels determine the amount of rough ER sheets in cells. Structure-guided mutagenesis and in vitro reconstitution on giant unilamellar vesicles further support a mechanism whereby SigmaR1 oligomers use their extended arrays of amphipathic helices to bind and flatten the lumenal leaflet of ER membranes to oppose membrane curvature and stabilize rough ER sheets.

Polysaccharide utilization loci encoded DUF1735 likely functions as membrane-bound spacer for carbohydrate active enzymes

Proteins featuring the Domain of Unknown Function 1735 are frequently found in Polysaccharide Utilization Loci, yet their role remains unknown. The domain and vicinity analyzer programs we developed mine the Kyoto Encyclopedia of Genes and Genomes and UniProt to enhance the functional prediction of DUF1735. Our datasets confirmed the exclusive presence of DUF1735 in Bacteroidota genomes, with Bacteroidetes thetaiotaomicron harboring 46 copies. Notably, 97.8% of DUF1735 are encoded in PULs, and 89% are N-termini of multimodular proteins featuring C-termini like Laminin_G_3, F5/8-typeC, and GH18 domains. Predominantly possessing a predicted lipoprotein signal peptide and sharing an immunoglobulin-like β-sandwich fold with the BACON domain and the N-termini of SusE/F, DUF1735 likely functions as N-terminal, membrane-bound spacer for diverse C-termini involved in PUL-mediated carbohydrate utilization.

The palisade layer of the poxvirus core is composed of flexible A10 trimers

Due to its asymmetric shape, size and compactness, the structure of the infectious mature virus (MV) of vaccinia virus (VACV), the best-studied poxvirus, remains poorly understood. Instead, subviral particles, in particular membrane-free viral cores, have been studied with cryo-electron microscopy. Here, we compared viral cores obtained by detergent stripping of MVs with cores in the cellular cytoplasm, early in infection. We focused on the prominent palisade layer on the core surface, combining cryo-electron tomography, subtomogram averaging and AlphaFold2 structure prediction. We showed that the palisade is composed of densely packed trimers of the major core protein A10. Trimers display a random order and their classification indicates structural flexibility. A10 on cytoplasmic cores is organized in a similar manner, indicating that the structures obtained in vitro are physiologically relevant. We discuss our results in the context of the VACV replicative cycle, and the assembly and disassembly of the infectious MV.

A protective and broadly binding antibody class engages the influenza virus hemagglutinin head at its stem interface

Influenza infection and vaccination impart strain-specific immunity that protects against neither seasonal antigenic variants nor the next pandemic. However, antibodies directed to conserved sites can confer broad protection. Here we identify and characterize a class of human antibodies that engage a previously undescribed, conserved epitope on the influenza hemagglutinin (HA) protein. Prototype antibody S8V1-157 binds at the normally occluded interface between the HA head and stem. Antibodies to this HA head-stem interface epitope are non-neutralizing in vitro but protect against lethal influenza infection in mice. Antibody isotypes that direct clearance of infected cells enhance this protection. Head-stem interface antibodies bind to most influenza A serotypes and seasonal human variants, and are present at low frequencies in the memory B cell populations of multiple human donors. Vaccines designed to elicit these antibodies might contribute to "universal" influenza immunity.

Computational screening of pathogenic missense nsSNPs in heme oxygenase 1 (HMOX1) gene and their structural and functional consequences

Heme Oxygenase 1 (HMOX1) is a cytoprotective enzyme, exhibiting the highest activity in the spleen, catalyzing the heme ring breakdown into products of biological significance- biliverdin, CO, and Fe2+. In vascular cells, HMOX1 possesses strong anti-apoptotic, antioxidant, anti-proliferative, anti-inflammatory, and immunomodulatory actions. The majority of these activities are crucial for the prevention of atherogenesis. Single amino acid substitutions in proteins generated by missense non-synonymous single nucleotide polymorphism (nsSNPs) in the protein-encoding regions of genes are potent enough to cause significant medical challenges due to the alteration of protein structure and function. The current study aimed at characterizing and analyzing high-risk nsSNPs associated with the human HMOX1 gene. Preliminary screening of the total available 288 missense SNPs was performed through the lens of deleteriousness and stability prediction tools. Finally, a total of seven nsSNPs (Y58D, A131T, Y134H, F166S, F167S, R183S and M186V) were found to be most deleterious by all tools that are present at highly conserved positions. Molecular dynamics simulations (MDS) analysis explained the mutational effects on the dynamic action of the wild-type and mutant proteins. In a nutshell, R183S (rs749644285) was identified as a highly detrimental mutation that could significantly render the enzymatic activity of HMOX1. The finding of this computational analysis might help subject the experimental confirmatory analysis to characterize the role of nsSNPs in HMOX1.Communicated by Ramaswamy H. Sarma.

Improved prediction of site-rates from structure with averaging across homologs

Variation in mutation rates at sites in proteins can largely be understood by the constraint that proteins must fold into stable structures. Models that calculate site-specific rates based on protein structure and a thermodynamic stability model have shown a significant but modest ability to predict empirical site-specific rates calculated from sequence. Models that use detailed atomistic models of protein energetics do not outperform simpler approaches using packing density. We demonstrate that a fundamental reason for this is that empirical site-specific rates are the result of the average effect of many different microenvironments in a phylogeny. By analyzing the results of evolutionary dynamics simulations, we show how averaging site-specific rates across many extant protein structures can lead to correct recovery of site-rate prediction. This result is also demonstrated in natural protein sequences and experimental structures. Using predicted structures, we demonstrate that atomistic models can improve upon contact density metrics in predicting site-specific rates from a structure. The results give fundamental insights into the factors governing the distribution of site-specific rates in protein families.

Deep mutational scanning quantifies DNA binding and predicts clinical outcomes of PAX6 variants

Nonsense and missense mutations in the transcription factor PAX6 cause a wide range of eye development defects, including aniridia, microphthalmia and coloboma. To understand how changes of PAX6:DNA binding cause these phenotypes, we combined saturation mutagenesis of the paired domain of PAX6 with a yeast one-hybrid (Y1H) assay in which expression of a PAX6-GAL4 fusion gene drives antibiotic resistance. We quantified binding of more than 2700 single amino-acid variants to two DNA sequence elements. Mutations in DNA-facing residues of the N-terminal subdomain and linker region were most detrimental, as were mutations to prolines and to negatively charged residues. Many variants caused sequence-specific molecular gain-of-function effects, including variants in position 71 that increased binding to the LE9 enhancer but decreased binding to a SELEX-derived binding site. In the absence of antibiotic selection, variants that retained DNA binding slowed yeast growth, likely because such variants perturbed the yeast transcriptome. Benchmarking against known patient variants and applying ACMG/AMP guidelines to variant classification, we obtained supporting-to-moderate evidence that 977 variants are likely pathogenic and 1306 are likely benign. Our analysis shows that most pathogenic mutations in the paired domain of PAX6 can be explained simply by the effects of these mutations on PAX6:DNA association, and establishes Y1H as a generalisable assay for the interpretation of variant effects in transcription factors.

Enzymatic characterization and thermostability improvement of an acidophilic endoxylanase PphXyn11 from Paenibacillus physcomitrellae XB

GH11 enzyme is known to be specific and efficient for the hydrolysis of xylan. It has been isolated from many microorganisms, and its enzymatic characteristics and thermostability vary between species. In this study, a GH11 enzyme PphXyn11 from a novel xylan-degrading strain of Paenibacillus physcomitrellae XB was characterized, and five mutants were constructed to try to improve the enzyme's thermostability. The results showed that PphXyn11 was an acidophilic endo-β-1,4-xylanase with the optimal reaction pH of 3.0-4.0, and it could deconstruct different kinds of xylan substrates efficiently, such as beechwood xylan, wheat arabinoxylan and xylo-oligosaccharides, to produce xylobiose and xylotriose as the main products at the optimal reaction temperature of 40 °C. Improvement of the thermal stability of PphXyn11 using site-directed mutagenesis revealed that three mutants, W33C/N47C, S127C/N174C and S49E, designed by adding the disulfide bonds at the N-terminal, C-terminal and increasing the charged residues on the surface of PphXyn11 respectively, could increase the enzymatic activity and thermal stablility significantly and make the optimal reaction temperature reach 50 °C. Molecular dynamics simulations as well as computed the numbers of salt bridges and hydrogen bonds indicated that the protein structures of these three mutants were more stable than the wild type, which provided theoretical support for their improved thermal stability. Certainly, further research is necessary to improve the enzymatic characteristics of PphXyn11 to achieve the bioconversion of hemicellulosic biomass on an applicable scale.

Biochemical, kinetic, and structural characterization of a Bacillus tequilensis nitroreductase

Nitroreductases (NRs) are NAD(P)H-dependent flavoenzymes that reduce nitro aromatic compounds to their corresponding arylamines via the nitroso and hydroxylamine intermediates. Because of their broad substrate scope and versatility, NRs have found application in multiple fields such as biocatalysis, bioremediation, cell-imaging and prodrug activation. However, only a limited number of members of the broad NR superfamily (> 24 000 sequences) have been experimentally characterized. Within this group of enzymes, only few are capable of amine synthesis, which is a fundamental chemical transformation for the pharmaceutical, agricultural, and textile industries. Herein, we provide a comprehensive description of a recently discovered NR from Bacillus tequilensis, named BtNR. This enzyme has previously been demonstrated to have the capability to fully convert nitro aromatic and heterocyclic compounds to their respective primary amines. In this study, we determined its biochemical, kinetic and structural properties, including its apparent melting temperature (Tm) of 59 °C, broad pH activity range (from pH 3 to 10) and a notably low redox potential (-236 ± 1 mV) in comparison to other well-known NRs. We also determined its steady-state and pre-steady-state kinetic parameters, which are consistent with other NRs. Additionally, we elucidated the crystal structure of BtNR, which resembles the well-characterized Escherichia coli oxygen-insensitive NAD(P)H nitroreductase (NfsB), and investigated the substrate binding in its active site through docking and molecular dynamics studies with four nitro aromatic substrates. Guided by these structural analyses, we probed the functional roles of active site residues by site-directed mutagenesis. Our findings provide valuable insights into the biochemical and structural properties of BtNR, as well as its potential applications in biotechnology.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

The eRF1 degrader SRI-41315 acts as a molecular glue at the ribosomal decoding center

Translation termination is an essential cellular process, which is also of therapeutic interest for diseases that manifest from premature stop codons. In eukaryotes, translation termination requires eRF1, which recognizes stop codons, catalyzes the release of nascent proteins from ribosomes and facilitates ribosome recycling. The small molecule SRI-41315 triggers eRF1 degradation and enhances translational readthrough of premature stop codons. However, the mechanism of action of SRI-41315 on eRF1 and translation is not known. Here we report cryo-EM structures showing that SRI-41315 acts as a metal-dependent molecular glue between the N domain of eRF1 responsible for stop codon recognition and the ribosomal subunit interface near the decoding center. Retention of eRF1 on ribosomes by SRI-41315 leads to ribosome collisions, eRF1 ubiquitylation and a higher frequency of translation termination at near-cognate stop codons. Our findings reveal a new mechanism of release factor inhibition and additional implications for pharmacologically targeting eRF1.

Identification of a Novel Homozygous GLS Gene Variant Associated with Developmental and Epileptic Encephalopathy (DEE) Type 71

Developmental and epileptic encephalopathy (DEEs) (OMIM#618,328) is characterized by seizures, hypotonia, and brain abnormalities, often arising from mutations in genes crucial for brain function. Among these genes, GLS stands out due to its vital role in the central nervous system (CNS), with homozygous variants potentially causing DEE type 71. Using Whole Exome Sequencing (WES) on a patient exhibiting symptoms of epileptic encephalopathy, we identified a novel homozygous variant, NM_014905.5:c.1849G > T; p.(Asp617Tyr), in the GLS gene. The 5-year-old patient, born to consanguineous parents, presented with developmental delay, encephalopathy, frequent seizures, and hypotonia. Sanger sequencing further validated the GLS gene variant in both the patient and his family. Furthermore, our bioinformatics analysis indicated that this missense variant could lead to alteration of splicing, resulting in the activation of a cryptic donor site and potentially causing loss of protein function. Our finding highlights the pathogenic significance of the GLS gene, particularly in the context of brain disorders, specifically DEE71.

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.

Late-onset Krabbe disease presenting as spastic paraplegia - implications of GCase and CTSB/D

OBJECTIVE

Krabbe disease (KD) is a multisystem neurodegenerative disorder with severe disability and premature death, mostly with an infancy/childhood onset. In rare cases of late-onset phenotypes, symptoms are often milder and difficult to diagnose. We here present a translational approach combining diagnostic and biochemical analyses of a male patient with a progressive gait disorder starting at the age of 44 years, with a final diagnosis of late-onset KD (LOKD).

METHODS

Additionally to cerebral MRI, protein structural analyses of the β-galactocerebrosidase protein (GALC) were performed. Moreover, expression, lysosomal localization, and activities of β-glucocerebrosidase (GCase), cathepsin B (CTSB), and cathepsin D (CTSD) were analyzed in leukocytes, fibroblasts, and lysosomes of fibroblasts.

RESULTS

Exome sequencing revealed biallelic likely pathogenic variants: GALC exons 11-17: 33 kb deletion; exon 4: missense variant (c.334A>G, p.Thr112Ala). We detected a reduced GALC activity in leukocytes and fibroblasts. While histological KD phenotypes were absent in fibroblasts, they showed a significantly decreased activities of GCase, CTSB, and CTSD in lysosomal fractions, while expression levels were unaffected.

INTERPRETATION

The presented LOKD case underlines the age-dependent appearance of a mildly pathogenic GALC variant and its interplay with other lysosomal proteins. As GALC malfunction results in reduced ceramide levels, we assume this to be causative for the here described decrease in CTSB and CTSD activity, potentially leading to diminished GCase activity. Hence, we emphasize the importance of a functional interplay between the lysosomal enzymes GALC, CTSB, CTSD, and GCase, as well as between their substrates, and propose their conjoined contribution in KD pathology.

Comprehensive bioinformatics analysis of structural and functional consequences of deleterious missense mutations in the human QDPR gene

Quinonoid dihydropteridine reductase (QDPR) is an enzyme that regulates tetrahydrobiopterin (BH4), a cofactor for enzymes involved in neurotransmitter synthesis and blood pressure regulation. Reduced QDPR activity can cause dihydrobiopterin (BH2) accumulation and BH4 depletion, leading to impaired neurotransmitter synthesis, oxidative stress, and increased risk of Parkinson's disease. A total of 10,236 SNPs were identified in the QDPR gene, with 217 being missense SNPs. Over 18 different sequence-based and structure-based tools were employed to assess the protein's biological activity, with several computational tools identifying deleterious SNPs. Additionally, the article provides detailed information about the QDPR gene and protein structure and conservation analysis. The results showed that 10 mutations were harmful and linked to brain and central nervous system disorders, and were predicted to be oncogenic by Dr. Cancer and CScape. Following conservation analysis, the HOPE server was used to analyse the effect of six selected mutations (L14P, V15G, G23S, V54G, M107K, G151S) on the protein structure. Overall, the study provides insights into the biological and functional impact of nsSNPs on QDPR activity and the potential induced pathogenicity and oncogenicity. In the future, research can be conducted to systematically evaluate QDPR gene variation through clinical studies, investigate mutation prevalence across different geographical regions, and validate computational results with conclusive experiments.Communicated by Ramaswamy H. Sarma.

Comparison of Xrn1 and Rat1 5' → 3' exoribonucleases in budding yeast supports the specific role of Xrn1 in cotranslational mRNA decay

The yeast Saccharomyces cerevisiae and most eukaryotes carry two 5' → 3' exoribonuclease paralogs. In yeast, they are called Xrn1, which shuttles between the nucleus and the cytoplasm, and executes major cytoplasmic messenger RNA (mRNA) decay, and Rat1, which carries a strong nuclear localization sequence (NLS) and localizes to the nucleus. Xrn1 is 30% identical to Rat1 but has an extra ~500 amino acids C-terminal extension. In the cytoplasm, Xrn1 can degrade decapped mRNAs during the last round of translation by ribosomes, a process referred to as "cotranslational mRNA decay." The division of labor between the two enzymes is still enigmatic and serves as a paradigm for the subfunctionalization of many other paralogs. Here we show that Rat1 is capable of functioning in cytoplasmic mRNA decay, provided that Rat1 remains cytoplasmic due to its NLS disruption (cRat1). This indicates that the physical segregation of the two paralogs plays roles in their specific functions. However, reversing segregation is not sufficient to fully complement the Xrn1 function. Specifically, cRat1 can partially restore the cell volume, mRNA stability, the proliferation rate, and 5' → 3' decay alterations that characterize xrn1Δ cells. Nevertheless, cotranslational decay is only slightly complemented by cRat1. The use of the AlphaFold prediction for cRat1 and its subsequent docking with the ribosome complex and the sequence conservation between cRat1 and Xrn1 suggest that the tight interaction with the ribosome observed for Xrn1 is not maintained in cRat1. Adding the Xrn1 C-terminal domain to Rat1 does not improve phenotypes, which indicates that lack of the C-terminal is not responsible for partial complementation. Overall, during evolution, it appears that the two paralogs have acquired specific characteristics to make functional partitioning beneficial.

Investigating the structural and functional consequences of germline single nucleotide polymorphisms located in the genes of the alternative lengthening of telomere (ALT) pathway

Background

The Alternative Lengthening of Telomeres (ALT) pathway represents a non-canonical mechanism of telomere maintenance that operates independently of the conventional telomerase activity. The three biologically significant proteins, designated as SMARCAL1 (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A-like protein 1), DAXX (Death domain-associated protein 6) and ATRX (alpha-thalassemia/mental retardation, X-linked) are associated with ALT in certain cancer types. The purpose of this study was to identify the most high-risk nsSNPs (non-synonymous Single Nucleotide Polymorphisms) within these three genes and assess their impacts on the structure and function of the proteins they encode.

Methods

The reported genetic polymorphisms of SMARCAL1, DAXX and ATRX genes were retrieved from the Ensembl database. Later, various computational tools like PROVEAN, PolyPhen2, SNPs and GO, SNAP2, Predict-SNP, Panther and PMut were used to predict the most deleterious nsSNPs. MutPred was used to understand the underlying molecular reasons of those nsSNPs being deleterious, followed by prediction of Post Translational Modification Sites (PTMs) using ModPred. I-Mutant and MUpro were used to predict the effect of SNP on energy stability. Later, 3D clustering analysis was done using Mutation 3D server. Moreover, ConSurf was utilized to identify the conservation scores of wild-type amino acids. Additionally, the NCBI conserved domain search tool was employed to pinpoint conserved domains within these three proteins. Project-Hope helped for biophysical validation, followed by prediction of these genes' interaction and function by using GeneMANIA.

Result

Analysis on SMARCAL1 protein revealed that among 665 nsSNPs, four were identified as the most deleterious: L578S, T581S, P582A, and P582S. Similarly, within the DAXX protein, among a pool of 480 nsSNPs, P284S, R230C, and R230S were found out to be the most deleterious variants. In case of ATRX protein, V178D, R246C, and V277G, from the total of 1009 nsSNPs, were predicted to be the most deleterious. All these nsSNPs were found to occur at residue positions that are 100 % conserved within protein domains and were predicted to be most damaging from both structural and functional perspectives and highly destabilizing to their corresponding proteins.

Conclusion

Computational investigation on the 3 proteins-SMARCAL1, DAXX and ATRX through different bioinformatics analysis tools concludes that the identified high risk nsSNPs of these proteins are pathogenic SNPs. These variants potentially exert functional and structural influences, thus making them valuable candidates for future genetic studies.

Structural characterization of macro domain-containing Thoeris antiphage defense systems

Thoeris defense systems protect bacteria from infection by phages via abortive infection. In these systems, ThsB proteins serve as sensors of infection and generate signaling nucleotides that activate ThsA effectors. Silent information regulator and SMF/DprA-LOG (SIR2-SLOG) containing ThsA effectors are activated by cyclic ADP-ribose (ADPR) isomers 2'cADPR and 3'cADPR, triggering abortive infection via nicotinamide adenine dinucleotide (NAD+) depletion. Here, we characterize Thoeris systems with transmembrane and macro domain (TM-macro)-containing ThsA effectors. We demonstrate that ThsA macro domains bind ADPR and imidazole adenine dinucleotide (IAD), but not 2'cADPR or 3'cADPR. Combining crystallography, in silico predictions, and site-directed mutagenesis, we show that ThsA macro domains form nucleotide-induced higher-order oligomers, enabling TM domain clustering. We demonstrate that ThsB can produce both ADPR and IAD, and we identify a ThsA TM-macro-specific ThsB subfamily with an active site resembling deoxy-nucleotide and deoxy-nucleoside processing enzymes. Collectively, our study demonstrates that Thoeris systems with SIR2-SLOG and TM-macro ThsA effectors trigger abortive infection via distinct mechanisms.

Major-groove sequence-specific RNA recognition by LoaP, a paralog of transcription elongation factor NusG

LoaP is a member of the universal NusG protein family. Previously, we reported that unlike other characterized homologs, LoaP binds RNA sequence-specifically, recognizing a stem-loop in the 5'-untranslated region of operons it regulates. To elucidate how this NusG homolog acquired this ability, we now determined the co-crystal structure of Thermoanaerobacter pseudethanolicus LoaP bound to its cognate 26-nucleotide dfn RNA element. Our structure reveals that the LoaP C-terminal KOW domain recognizes the helical portion of the RNA by docking into a broadened major groove, while a protruding β-hairpin of the N-terminal NusG-like domain binds the UNCG tetraloop capping the stem-loop. Major-groove RNA recognition is unusual and is made possible by conserved features of the dfn hairpin. Superposition with structures of other NusG proteins implies that LoaP can bind concurrently to the dfn RNA and the transcription elongation complex, suggesting a new level of co-transcriptional regulation by proteins of this conserved family.

Deafness DFNB128 Associated with a Recessive Variant of Human MAP3K1 Recapitulates Hearing Loss of Map3k1-Deficient Mice

Deafness in vertebrates is associated with variants of hundreds of genes. Yet, many mutant genes causing rare forms of deafness remain to be discovered. A consanguineous Pakistani family segregating nonsyndromic deafness in two sibships were studied using microarrays and exome sequencing. A 1.2 Mb locus (DFNB128) on chromosome 5q11.2 encompassing six genes was identified. In one of the two sibships of this family, a novel homozygous recessive variant NM_005921.2:c.4460G>A p.(Arg1487His) in the kinase domain of MAP3K1 co-segregated with nonsyndromic deafness. There are two previously reported Map3k1-kinase-deficient mouse models that are associated with recessively inherited syndromic deafness. MAP3K1 phosphorylates serine and threonine and functions in a signaling pathway where pathogenic variants of HGF, MET, and GAB1 were previously reported to be associated with human deafness DFNB39, DFNB97, and DFNB26, respectively. Our single-cell transcriptome data of mouse cochlea mRNA show expression of Map3k1 and its signaling partners in several inner ear cell types suggesting a requirement of wild-type MAP3K1 for normal hearing. In contrast to dominant variants of MAP3K1 associated with Disorders of Sex Development 46,XY sex-reversal, our computational modeling of the recessive substitution p.(Arg1487His) predicts a subtle structural alteration in MAP3K1, consistent with the limited phenotype of nonsyndromic deafness.

Stress-dependent condensate formation regulated by the ubiquitin-related modifier Urm1

The ability of proteins and RNA to coalesce into phase-separated assemblies, such as the nucleolus and stress granules, is a basic principle in organizing membraneless cellular compartments. While the constituents of biomolecular condensates are generally well documented, the mechanisms underlying their formation under stress are only partially understood. Here, we show in yeast that covalent modification with the ubiquitin-like modifier Urm1 promotes the phase separation of a wide range of proteins. We find that the drop in cellular pH induced by stress triggers Urm1 self-association and its interaction with both target proteins and the Urm1-conjugating enzyme Uba4. Urmylation of stress-sensitive proteins promotes their deposition into stress granules and nuclear condensates. Yeast cells lacking Urm1 exhibit condensate defects that manifest in reduced stress resilience. We propose that Urm1 acts as a reversible molecular "adhesive" to drive protective phase separation of functionally critical proteins under cellular stress.

Borg extrachromosomal elements of methane-oxidizing archaea have conserved and expressed genetic repertoires

Borgs are huge extrachromosomal elements (ECE) of anaerobic methane-consuming "Candidatus Methanoperedens" archaea. Here, we used nanopore sequencing to validate published complete genomes curated from short reads and to reconstruct new genomes. 13 complete and four near-complete linear genomes share 40 genes that define a largely syntenous genome backbone. We use these conserved genes to identify new Borgs from peatland soil and to delineate Borg phylogeny, revealing two major clades. Remarkably, Borg genes encoding nanowire-like electron-transferring cytochromes and cell surface proteins are more highly expressed than those of host Methanoperedens, indicating that Borgs augment the Methanoperedens activity in situ. We reconstructed the first complete 4.00 Mbp genome for a Methanoperedens that is inferred to be a Borg host and predicted its methylation motifs, which differ from pervasive TC and CC methylation motifs of the Borgs. Thus, methylation may enable Methanoperedens to distinguish their genomes from those of Borgs. Very high Borg to Methanoperedens ratios and structural predictions suggest that Borgs may be capable of encapsulation. The findings clearly define Borgs as a distinct class of ECE with shared genomic signatures, establish their diversification from a common ancestor with genetic inheritance, and raise the possibility of periodic existence outside of host cells.

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.

Exploring the effect of disease causing mutations in metal binding sites of human ARSA in metachromatic leukodystrophy

The arylsulfatase A (ARSA) gene is observed to be deficient in patients with metachromatic leukodystrophy (MLD), a type of lysosomal storage disease. MLD is a severe neurodegenerative disorder characterized by an autosomal recessive inheritance pattern. This study aimed to map the most deleterious mutations at the metal binding sites of ARSA and the amino acids in proximity to the mutated positions. We utilized an array of computational tools, including PredictSNP, MAPP, PhD-SNP, PolyPhen-1, PolyPhen-2, SIFT, SNAP, and ConSurf, to identify the most detrimental mutations potentially implicated in MLD collected from UniProt, ClinVar, and HGMD. Two mutations, D29N and D30H, as being extremely deleterious based on assessments of pathogenicity, conservation, biophysical characteristics, and stability analysis. The D29 and D30 are located at the metal-interacting regions of ARSA and found to undergo post-translational modification, specifically phosphorylation. Henceforth, the in-depth effect of metal binding upon mutation was examined using molecular dynamics simulations (MDS) before and after phosphorylation. The MDS results exhibited high deviation for the D29N and D30H mutations in comparison to the native, and the same was confirmed by significant residue fluctuation and reduced compactness. These structural alterations suggest that such mutations may influence protein functionality, offering potential avenues for personalized therapeutic and providing a basis for potential mutation-specific treatments for severe MLD patients.

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.

Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

Molecular characterization of Rft1, an ER membrane protein associated with congenital disorder of glycosylation RFT1-CDG

The oligosaccharide needed for protein N-glycosylation is assembled on a lipid carrier via a multi-step pathway. Synthesis is initiated on the cytoplasmic face of the endoplasmic reticulum (ER) and completed on the luminal side after transbilayer translocation of a heptasaccharide lipid intermediate. More than 30 Congenital Disorders of Glycosylation (CDGs) are associated with this pathway, including RFT1-CDG which results from defects in the membrane protein Rft1. Rft1 is essential for the viability of yeast and mammalian cells and was proposed as the transporter needed to flip the heptasaccharide lipid intermediate across the ER membrane. However, other studies indicated that Rft1 is not required for heptasaccharide lipid flipping in microsomes or unilamellar vesicles reconstituted with ER membrane proteins, nor is it required for the viability of at least one eukaryote. It is therefore not known what essential role Rft1 plays in N-glycosylation. Here, we present a molecular characterization of human Rft1, using yeast cells as a reporter system. We show that it is a multi-spanning membrane protein located in the ER, with its N and C-termini facing the cytoplasm. It is not N-glycosylated. The majority of RFT1-CDG mutations map to highly conserved regions of the protein. We identify key residues that are important for Rft1's ability to support N-glycosylation and cell viability. Our results provide a necessary platform for future work on this enigmatic protein.

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.