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

Molecular mechanism of cellulose depolymerization by the two-domain BlCel9A enzyme from the glycoside hydrolase family 9

Carbohydrate-active enzymes from the glycoside hydrolase family 9 (GH9) play a key role in processing lignocellulosic biomass. Although the structural features of some GH9 enzymes are known, the molecular mechanisms that drive their interactions with cellulosic substrates remain unclear. To investigate the molecular mechanisms that the two-domain Bacillus licheniformis BlCel9A enzyme utilizes to depolymerize cellulosic substrates, we used a combination of biochemical assays, X-ray crystallography, small-angle X-ray scattering, and molecular dynamics simulations. The results reveal that BlCel9A breaks down cellulosic substrates, releasing cellobiose and glucose as the major products, but is highly inefficient in cleaving oligosaccharides shorter than cellotetraose. In addition, fungal lytic polysaccharide oxygenase (LPMO) TtLPMO9H enhances depolymerization of crystalline cellulose by BlCel9A, while exhibiting minimal impact on amorphous cellulose. The crystal structures of BlCel9A in both apo form and bound to cellotriose and cellohexaose were elucidated, unveiling the interactions of BlCel9A with the ligands and their contribution to substrate binding and products release. MD simulation analysis reveals that BlCel9A exhibits higher interdomain flexibility under acidic conditions, and SAXS experiments indicate that the enzyme flexibility is induced by pH and/or temperature. Our findings provide new insights into BlCel9A substrate specificity and binding, and synergy with the LPMOs.

Oligomerization mediated by the D2 domain of DTX3L is critical for DTX3L-PARP9 reading function of mono-ADP-ribosylated androgen receptor

Deltex proteins are a family of E3 ubiquitin ligases that encode C-terminal RING and DTC domains that mediate interactions with E2 ubiquitin-conjugating enzymes and recognize ubiquitination substrates. DTX3L is unique among the Deltex proteins based on its N-terminal domain architecture. The N-terminal D1 and D2 domains of DTX3L mediate homo-oligomerization, and the D3 domain interacts with PARP9, a protein that contains tandem macrodomains with ADP-ribose reader function. While DTX3L and PARP9 are known to heterodimerize, and assemble into a high molecular weight oligomeric complex, the nature of the oligomeric structure, including whether this contributes to the ADP-ribose reader function is unknown. Here, we report a crystal structure of the DTX3L N-terminal D2 domain and show that it forms a tetramer with, conveniently, D2 symmetry. We identified two interfaces in the structure: a major, conserved interface with a surface of 973 Å2 and a smaller one of 415 Å2. Using native mass spectrometry, we observed molecular species that correspond to monomers, dimers and tetramers of the D2 domain. Reconstitution of DTX3L knockout cells with a D1-D2 deletion mutant showed the domain is dispensable for DTX3L-PARP9 heterodimer formation, but necessary to assemble an oligomeric complex with efficient reader function for ADP-ribosylated androgen receptor. Our results suggest that homo-oligomerization of DTX3L is important for the DTX3L-PARP9 complex to read mono-ADP-ribosylation on a ligand-regulated transcription factor.

Genetic variation in the BLM gene and its expression in the ovaries is closely related to kidding number in goats

Bloom (BLM) helicase plays an important role in DNA replication and the maintenance of genome integrity. BLM protein deficiency, which plays a vital role in the sperm-egg union and germ-cell development during reproduction, can lead to severe DNA damage in goats. However, the effect of BLM protein deficiency on goat litter size has not been reported. Herein, we studied the association between the genetic variation in the BLM gene and the number of kids per litter in Guizhou white goats. We explored differences in the expression of the BLM protein in the follicles of single and multi-kid nanny goats. We also analyzed the effects of dysregulated BLM gene expression on the proliferation and apoptosis of ovarian granulosa cells and the expression of genes related to follicle development in goats. Five single nucleotide polymorphism (SNP) loci, including the non-synonymous mutations g.38179 A > G, g.40626 G > C and g.89621 T > G; the intron synonymous mutation g.56961 G > A and the exon synonymous mutation g.65796 C > T were found in the BLM gene. All SNPs loci were in Hardy-Weinberg equilibrium, and correlation analysis showed that the g.65796 C > T and g.89621 T > G loci polymorphism was strongly associated with litter size in the first three litters (P < 0.05). The diplogenotype Hap 2/2 (AAGGAACCTT) showed no significant difference in litter size between different births, indicating that the diploid genotype is stable in different litter sizes. Bioinformatics analysis showed that three non-synonymous mutation loci (p.T488A, p.A662S, and p.S1373A) could affect BLM protein stability, and mutations in p.T488A and p.S1373A led to changes in amino acid polarity and associated interactions. qPCR results showed that the expression level of the BLM gene in the uterus and ovaries of TT genotype nanny goats was significantly higher than that of GG genotype nanny goats. Indirect immunofluorescence assay (IF) showed that the BLM protein was significantly overexpressed in both the primordial and growing follicles of nanny goats with multiple kids (P < 0.01). Disrupting BLM gene expression in the ovarian granulosa cells down-regulated the expression of the Cyp19A1 gene. It also significantly inhibited the proliferation of follicles and induces early apoptosis of the granulosa cells. These findings confirm that polymorphism in the BLM gene is closely related to the littering traits of Guizhou white goats, and it affects the reproductive performance of nanny goats by regulating the development of the oocytes and granulosa cells. This work provides new evidence on the regulatory effect of the BLM gene on the litter size of nanny goats.

Benchmarking the accuracy of structure-based binding affinity predictors on Spike-ACE2 deep mutational interaction set

Since the start of COVID-19 pandemic, a huge effort has been devoted to understanding the Spike (SARS-CoV-2)-ACE2 recognition mechanism. To this end, two deep mutational scanning studies traced the impact of all possible mutations across receptor binding domain (RBD) of Spike and catalytic domain of human ACE2. By concentrating on the interface mutations of these experimental data, we benchmarked six commonly used structure-based binding affinity predictors (FoldX, EvoEF1, MutaBind2, SSIPe, HADDOCK, and UEP). These predictors were selected based on their user-friendliness, accessibility, and speed. As a result of our benchmarking efforts, we observed that none of the methods could generate a meaningful correlation with the experimental binding data. The best correlation is achieved by FoldX (R = -0.51). When we simplified the prediction problem to a binary classification, that is, whether a mutation is enriching or depleting the binding, we showed that the highest accuracy is achieved by FoldX with a 64% success rate. Surprisingly, on this set, simple energetic scoring functions performed significantly better than the ones using extra evolutionary-based terms, as in Mutabind and SSIPe. Furthermore, we demonstrated that recent AI approaches, mmCSM-PPI and TopNetTree, yielded comparable performances to the force field-based techniques. These observations suggest plenty of room to improve the binding affinity predictors in guessing the variant-induced binding profile changes of a host-pathogen system, such as Spike-ACE2. To aid such improvements we provide our benchmarking data at https://github.com/CSB-KaracaLab/RBD-ACE2-MutBench with the option to visualize our mutant models at https://rbd-ace2-mutbench.github.io/.

Preclinical Targeting of the PGRMC1-CK2 Axis with Silmitasertib: A Potential Strategy for Lung Adenocarcinoma Therapy

Progesterone receptor membrane component 1 (PGRMC1) is a pleiotropic protein over-expressed in lung adenocarcinoma (LUAD). The precise molecular mechanisms underlying the signature motif of Casein kinase (CK2) presence in PGRMC1 and their role in LUAD remain unclear. X-ray crystallographic structure for CK2 and PGRMC1 from the PubChem database was obtained and subjected to protein-protein interaction (PPI) analysis to identify their interactions. In addition, the CK2 inhibitor - Silmitasertib was also utilised to understand the interaction between PGRMC1-CK2. The PPI complex (PGRMC1-CK2) and the PPI-ligand interaction analysis and their Molecular Dynamics (MD) studies revealed the stability of their interactions and critical amino acid contacts within the 5Ǻ vicinity of the CK2 signature motif "T/S-x-x-E/D". Moreover, in-vitro colony formation assay, migration assay, and gene expression analysis using quantitative Real-time PCR revealed that Silmitasertib (IC50-2.5 μM) was highly influential in suppressing the PGRMC1-CK2 expression axis. In conclusion, our study infers that PGRMC1-CK-2 axis inhibition could be a potential therapeutic option to limit the promotion and progression of lung cancer.

Growth-Coupled Evolutionary Pressure Improving Epimerases for D-Allulose Biosynthesis Using a Biosensor-Assisted In Vivo Selection Platform

Fast screening strategies that enable high-throughput evaluation and identification of desired variants from diversified enzyme libraries are crucial to tailoring biocatalysts for the synthesis of D-allulose, which is currently limited by the poor catalytic performance of ketose 3-epimerases (KEases). Here, the study designs a minimally equipment-dependent, high-throughput, and growth-coupled in vivo screening platform founded on a redesigned D-allulose-dependent biosensor system. The genetic elements modulating regulator PsiR expression levels undergo systematic optimization to improve the growth-responsive dynamic range of the biosensor, which presents ≈30-fold facilitated growth optical density with a high signal-to-noise ratio (1.52 to 0.05) toward D-allulose concentrations from 0 to 100 mm. Structural analysis and evolutionary conservation analysis of Agrobacterium sp. SUL3 D-allulose 3-epimerase (ADAE) reveal a highly conserved catalytic active site and variable hydrophobic pocket, which together regulate substrate recognition. Structure-guided rational design and directed evolution are implemented using the growth-coupled in vivo screening platform to reprogram ADAE, in which a mutant M42 (P38N/V102A/Y201L/S207N/I251R) is identified with a 6.28-fold enhancement of catalytic activity and significantly improved thermostability with a 2.5-fold increase of the half-life at 60 °C. The research demonstrates that biosensor-assisted growth-coupled evolutionary pressure combined with structure-guided rational design provides a universal route for engineering KEases.

The deleterious variants of N-acetylgalactosamine-6-sulfatase (GalN6S) enzyme trigger Morquio a syndrome by disrupting protein foldings

Lysosomal enzymes degrade cellular macromolecules, while their inactivation causes human hereditary metabolic disorders. Mucopolysaccharidosis IVA (MPS IVA; Moquio A syndrome) is one of the lysosomal storage disorders caused by a defective Galactosamine-6-sulfatase (GalN6S) enzyme. In several populations, disease incidence is elevated due to missense mutations brought on by non-synonymous allelic variation in the GalN6S enzyme. Here, we studied the effect of non-synonymous single nucleotide polymorphism (nsSNPs) on the structural dynamics of the GalN6S enzyme and its binding with N-acetylgalactosamine (GalNAc) using all-atom molecular dynamics simulation and an essential dynamics approach. Consequently, in this study, we have identified three functionally disruptive mutations in domain-I and domain-II, that is, S80L, R90W, and S162F, which presumably contribute to post-translational modifications. The study delineated that both domains work cooperatively, and alteration in domain II (S80L, R90W) leads to conformational changes in the catalytic site in domain-I, while mutation S162F mainly provokes higher residual flexibility of domain II. These results show that these mutations impair the hydrophobic core, implying that Morquio A syndrome is caused by misfolding of the GalN6S enzyme. The results also show the instability of the GalN6S-GalNAc complex upon substitution. Overall, the structural dynamics resulting from point mutations give the molecular rationale for Moquio A syndrome and, more importantly, the Mucopolysaccharidoses (MPS) family of diseases, re-establishing MPS IVA as a protein-folding disease.Communicated by Ramaswamy H. Sarma.

A computational approach to analyzing the functional and structural impacts of Tripeptidyl-Peptidase 1 missense mutations in neuronal ceroid lipofuscinosis

Neuronal ceroid-lipofuscinosis (NCLs) are a group of severe neurodegenerative conditions, most likely present in infantile, late infantile, juvenile, and adult-onset forms. Their phenotypic characteristics comprise eyesight damage, reduced motor activity and cognitive function, and sometimes tend to die in the initial stage. In recent studies, NCLs have been categorized into at least 14 genetic collections (CLN1-14). CLN2 gene encodes Tripeptidyl peptidase 1 (TPP1), which affects late infantile-onset form. In this study, we retrieved a mutational dataset screening for TPP1 protein from various databases (ClinVar, UniProt, HGMD). Fifty-six missense mutants were enumerated with computational methods to perceive the significant mutants (G475R and G501C) and correlated with clinical and literature data. A structure-based screening method was initiated to understand protein-ligand interaction and dynamic simulation. The docking procedure was performed for the native (3EDY) and mutant (G473R and G501C) structures with Gemfibrozil (gem), which lowers the lipid level, decreases the triglycerides amount in the blood circulation, and controls hyperlipidemia. The Native had an interaction score of -5.57 kcal/mol, and the mutants had respective average binding scores of -6.24 (G473R) and - 5.17 (G501C) kcal/mol. Finally, molecular dynamics simulation showed that G473R and G501C mutants had better flexible and stable orientation in all trajectory analyses. Therefore, this work gives an extended understanding of both functional and structural levels of influence for the mutant form that leads to NCL disorder.

AlphaFold-assisted structure determination of a bacterial protein of unknown function using X-ray and electron crystallography

Macromolecular crystallography generally requires the recovery of missing phase information from diffraction data to reconstruct an electron-density map of the crystallized molecule. Most recent structures have been solved using molecular replacement as a phasing method, requiring an a priori structure that is closely related to the target protein to serve as a search model; when no such search model exists, molecular replacement is not possible. New advances in computational machine-learning methods, however, have resulted in major advances in protein structure predictions from sequence information. Methods that generate predicted structural models of sufficient accuracy provide a powerful approach to molecular replacement. Taking advantage of these advances, AlphaFold predictions were applied to enable structure determination of a bacterial protein of unknown function (UniProtKB Q63NT7, NCBI locus BPSS0212) based on diffraction data that had evaded phasing attempts using MIR and anomalous scattering methods. Using both X-ray and micro-electron (microED) diffraction data, it was possible to solve the structure of the main fragment of the protein using a predicted model of that domain as a starting point. The use of predicted structural models importantly expands the promise of electron diffraction, where structure determination relies critically on molecular replacement.

Heterozygous gain of function variant in GUCY1A2 may cause autonomous ovarian hyperfunction

PURPOSE

The purpose of this study was to characterize the phenotype associated with a de novo gain-of-function variant in the GUCY1A2 gene.

METHODS

An individual carrying the de novo heterozygous variant c.1458G>T p.(E486D) in GUCY1A2 was identified by exome sequencing. The effect of the corresponding enzyme variant α2E486D/β1 was evaluated using concentration-response measurements with wild-type enzyme and the variant in cytosolic fractions of HEK293 cells, UV-vis absorbance spectra of the corresponding purified enzymes, and examination of overexpressed fluorescent protein-tagged constructs by confocal laser scanning microscopy.

RESULTS

The patient presented with precocious peripheral puberty resembling the autonomous ovarian puberty seen in McCune-Albright syndrome. Additionally, the patient displayed severe intellectual disability. In vitro activity assays revealed an increased nitric oxide affinity for the mutant enzyme. The response to carbon monoxide was unchanged, while thermostability was decreased compared to wild type. Heme content, susceptibility to oxidation, and subcellular localization upon overexpression were unchanged.

CONCLUSION

Our data define a syndromic autonomous ovarian puberty likely due to the activating allele p.(E486D) in GUCY1A2 leading to an increase in cGMP. The overlap with the ovarian symptoms of McCune-Albright syndrome suggests an impact of this cGMP increase on the cAMP pathway in the ovary. Additional cases will be needed to ensure a causal link.

New insights into GATOR2-dependent interactions and its conformational changes in amino acid sensing

Eukaryotic cells coordinate growth under different environmental conditions via mechanistic target of rapamycin complex 1 (mTORC1). In the amino-acid-sensing signalling pathway, the GATOR2 complex, containing five evolutionarily conserved subunits (WDR59, Mios, WDR24, Seh1L and Sec13), is required to regulate mTORC1 activity by interacting with upstream CASTOR1 (arginine sensor) and Sestrin2 (leucine sensor and downstream GATOR1 complex). GATOR2 complex utilizes β-propellers to engage with CASTOR1, Sestrin2 and GATOR1, removal of these β-propellers results in substantial loss of mTORC1 capacity. However, structural information regarding the interface between amino acid sensors and GATOR2 remains elusive. With the recent progress of the AI-based tool AlphaFold2 (AF2) for protein structure prediction, structural models were predicted for Sentrin2-WDR24-Seh1L and CASTOR1-Mios β-propeller. Furthermore, the effectiveness of relevant residues within the interface was examined using biochemical experiments combined with molecular dynamics (MD) simulations. Notably, fluorescence resonance energy transfer (FRET) analysis detected the structural transition of GATOR2 in response to amino acid signals, and the deletion of Mios β-propeller severely impeded that change at distinct arginine levels. These findings provide structural perspectives on the association between GATOR2 and amino acid sensors and can facilitate future research on structure determination and function.

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

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

Ultralow Background Membrane Editors for Spatiotemporal Control of Phosphatidic Acid Metabolism and Signaling

Phosphatidic acid (PA) is a multifunctional lipid with important metabolic and signaling functions, and efforts to dissect its pleiotropy demand strategies for perturbing its levels with spatiotemporal precision. Previous membrane editing approaches for generating local PA pools used light-mediated induced proximity to recruit a PA-synthesizing enzyme, phospholipase D (PLD), from the cytosol to the target organelle membrane. Whereas these optogenetic PLDs exhibited high activity, their residual activity in the dark led to undesired chronic lipid production. Here, we report ultralow background membrane editors for PA wherein light directly controls PLD catalytic activity, as opposed to localization and access to substrates, exploiting a light-oxygen-voltage (LOV) domain-based conformational photoswitch inserted into the PLD sequence and enabling their stable and nonperturbative targeting to multiple organelle membranes. By coupling organelle-targeted LOVPLD activation to lipidomics analysis, we discovered different rates of metabolism for PA and its downstream products depending on the subcellular location of PA production. We also elucidated signaling roles for PA pools on different membranes in conferring local activation of AMP-activated protein kinase signaling. This work illustrates how membrane editors featuring acute, optogenetic conformational switches can provide new insights into organelle-selective lipid metabolic and signaling pathways.

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

Monoamine transporters including transporters for serotonin, dopamine, and norepinephrine play key roles in monoaminergic synaptic signaling, involving in the molecular etiology of a wide range of neurological and physiological disorders. Despite being crucial drug targets, the study of transmembrane proteins remains challenging due to their localization within the cell membrane. To address this, we present the structural bioinformatics studies of 7 monoamine transporters and their water-soluble variants designed using the QTY code, by systematically replacing the hydrophobic amino acids leucine (L), valine (V), isoleucine (I) and phenylalanine (F) with hydrophilic amino acids (glutamine (Q), threonine (T) and tyrosine (Y). The resulting QTY variants, despite significant protein transmembrane sequence differences (44.27%-51.85%), showed similar isoelectric points (pI) and molecular weights. While their hydrophobic surfaces significantly reduced, this change resulted in a minimal structural alteration. Quantitatively, Alphafold2 predicted QTY variant structures displayed remarkable similarity with RMSD 0.492Å-1.619Å. Accompanied by the structural similarities of substituted amino acids in the context of 1.5Å electron density maps, our study revealed multiple QTY and reverse QTY variations in genomic databases. We further analyzed their phenotypical and topological characteristics. By extending evolutionary game theory to the molecular foundations of biology, we provided insights into the evolutionary dynamics of chemically distinct alpha-helices, their usage in different chemotherapeutic applications, and open possibilities of diagnostic medicine. Our study rationalizes that QTY variants of monoamine transporters may not only become distinct tools for medical, structural, and evolutionary research, but these transporters may also emerge as contemporary therapeutic targets, providing a new approach to treatment for several conditions.

Structural basis for peroxidase encapsulation inside the encapsulin from the Gram-negative pathogen Klebsiella pneumoniae

Encapsulins are self-assembling protein nanocompartments capable of selectively encapsulating dedicated cargo proteins, including enzymes involved in iron storage, sulfur metabolism, and stress resistance. They represent a unique compartmentalization strategy used by many pathogens to facilitate specialized metabolic capabilities. Encapsulation is mediated by specific cargo protein motifs known as targeting peptides (TPs), though the structural basis for encapsulation of the largest encapsulin cargo class, dye-decolorizing peroxidases (DyPs), is currently unknown. Here, we characterize a DyP-containing encapsulin from the enterobacterial pathogen Klebsiella pneumoniae. By combining cryo-electron microscopy with TP and TP-binding site mutagenesis, we elucidate the molecular basis for cargo encapsulation. TP binding is mediated by cooperative hydrophobic and ionic interactions as well as shape complementarity. Our results expand the molecular understanding of enzyme encapsulation inside protein nanocompartments and lay the foundation for rationally modulating encapsulin cargo loading for biomedical and biotechnological applications.

Triaging of α-helical proteins to the mitochondrial outer membrane by distinct chaperone machinery based on substrate topology

Mitochondrial outer membrane ⍺-helical proteins play critical roles in mitochondrial-cytoplasmic communication, but the rules governing the targeting and insertion of these biophysically diverse proteins remain unknown. Here, we first defined the complement of required mammalian biogenesis machinery through genome-wide CRISPRi screens using topologically distinct membrane proteins. Systematic analysis of nine identified factors across 21 diverse ⍺-helical substrates reveals that these components are organized into distinct targeting pathways that act on substrates based on their topology. NAC is required for the efficient targeting of polytopic proteins, whereas signal-anchored proteins require TTC1, a cytosolic chaperone that physically engages substrates. Biochemical and mutational studies reveal that TTC1 employs a conserved TPR domain and a hydrophobic groove in its C-terminal domain to support substrate solubilization and insertion into mitochondria. Thus, the targeting of diverse mitochondrial membrane proteins is achieved through topological triaging in the cytosol using principles with similarities to ER membrane protein biogenesis systems.

Molecular dynamics of the human RhD and RhAG blood group proteins

Introduction: Blood group antigens of the RH system (formerly known as "Rhesus") play an important role in transfusion medicine because of the severe haemolytic consequences of antibodies to these antigens. No crystal structure is available for RhD proteins with its partner RhAG, and the precise stoichiometry of the trimer complex remains unknown. Methods: To analyse their structural properties, the trimers formed by RhD and/or RhAG subunits were generated by protein modelling and molecular dynamics simulations were performed. Results: No major differences in structural behaviour were found between trimers of different compositions. The conformation of the subunits is relatively constant during molecular dynamics simulations, except for three large disordered loops. Discussion: This work makes it possible to propose a reasonable stoichiometry and demonstrates the potential of studying the structural behaviour of these proteins to investigate the hundreds of genetic variants relevant to transfusion medicine.

Autosomal Recessive Rod-Cone Dystrophy with Mild Extra-Ocular Manifestations Due to a Splice-Affecting Variant in BBS9

Bardet-Biedl syndrome (BBS), one of the most common forms of syndromic inherited retinal diseases (IRDs), is characterized by the combination of retinal degeneration with additional extra-ocular manifestations, including obesity, intellectual disability, kidney disease, polydactyly and other skeletal abnormalities. We observed an Israeli patient with autosomal recessive apparently non-syndromic rod-cone dystrophy (RCD). Extra-ocular findings were limited to epilepsy and dental problems. Genetic analysis with a single molecule molecular inversion probes-based panel that targets the exons and splice sites of 113 genes associated with retinitis pigmentosa and Leber congenital amaurosis revealed a homozygous rare missense variant in the BBS9 gene (c.263C>T;p.(Ser88Leu)). This variant, which affects a highly conserved amino acid, is also located in the last base of Exon 3, and predicted to be splice-altering. An in vitro minigene splice assay demonstrated that this variant leads to the partial aberrant splicing of Exon 3. Therefore, we suggest that this variant is likely hypomorphic. This is in agreement with the relatively mild phenotype observed in the patient. Hence, the findings in our study expand the phenotypic spectrum associated with BBS9 variants and indicate that variants in this gene should be considered not only in BBS patients but also in individuals with non-syndromic IRD or IRD with very mild extra-ocular manifestations.

Exploring the role of TLK2 mutation in tropical calcific pancreatitis: an in silico and molecular dynamics simulation study

Tropical calcific pancreatitis (TCP) is a juvenile form of non-alcoholic chronic pancreatitis seen exclusively in tropical countries. The disease poses a high risk of complications, including pancreatic diabetes and cancer, leading to significant mortality due to poor diagnosis and ineffective treatments. This study employed whole exome sequencing (WES) of 5 TCP patient samples to identify genetic variants associated with TCP. Advanced computational techniques were used to gain atomic-level insights into disease progression, including microsecond-scale long MD simulations and essential dynamics. In silico virtual screening was performed to identify potential therapeutic compounds targeting the mutant protein using the Asinex and DrugBank compound library. WES analysis predicted several single nucleotide variants (SNVs) associated with TCP, including a novel missense variant (c.T1802A or p.V601E) in the TLK2 gene. Computational analysis revealed that the p.V601E mutation significantly affected the structure of the TLK2 kinase domain and its conformational dynamics, altering the interaction profile between ATP and the binding pocket. These changes could impact TLK2's kinase activity and functions, potentially correlating with TCP progression. Promising lead compounds that selectively bind to the TLK2 mutant protein were identified, offering potential for therapeutic interventions in TCP. These findings hold great potential for future research.Communicated by Ramaswamy H. Sarma.

Structural and mechanistic insights into activation of the human RNA ligase RTCB by Archease

RNA ligases of the RTCB-type play an essential role in tRNA splicing, the unfolded protein response and RNA repair. RTCB is the catalytic subunit of the pentameric human tRNA ligase complex. RNA ligation by the tRNA ligase complex requires GTP-dependent activation of RTCB. This active site guanylylation reaction relies on the activation factor Archease. The mechanistic interplay between both proteins has remained unknown. Here, we report a biochemical and structural analysis of the human RTCB-Archease complex in the pre- and post-activation state. Archease reaches into the active site of RTCB and promotes the formation of a covalent RTCB-GMP intermediate through coordination of GTP and metal ions. During the activation reaction, Archease prevents futile RNA substrate binding to RTCB. Moreover, monomer structures of Archease and RTCB reveal additional states within the RNA ligation mechanism. Taken together, we present structural snapshots along the reaction cycle of the human tRNA ligase.

PTCH1 Gene Variants, mRNA Expression, and Bioinformatics Insights in Mexican Cutaneous Squamous Cell Carcinoma Patients

BACKGROUND

Skin cancer is one of the most frequent types of cancer, and cutaneous squamous cell carcinoma (cSCC) constitutes 20% of non-melanoma skin cancer (NMSC) cases. PTCH1, a tumor suppressor gene involved in the Sonic hedgehog signaling pathway, plays a crucial role in neoplastic processes.

METHODS

An analytical cross-sectional study, encompassing 211 cSCC patients and 290 individuals in a control group (CG), was performed. A subgroup of samples was considered for the relative expression analysis, and the results were obtained using quantitative real-time PCR (qPCR) with TaqMan® probes. The functional, splicing, and disease-causing effects of the proposed variants were explored via bioinformatics.

RESULTS

cSCC was predominant in men, especially in sun-exposed areas such as the head and neck. No statistically significant differences were found regarding the rs357564, rs2236405, rs2297086, and rs41313327 variants of PTCH1, or in the risk of cSCC, nor in the mRNA expression between the cSCC group and CG. A functional effect of rs357564 and a disease-causing relation to rs41313327 was identified.

CONCLUSION

The proposed variants were not associated with cSCC risk in this Mexican population, but we recognize the need for analyzing larger population groups to elucidate the disease-causing role of rare variants.

Targeting imidazole-glycerol phosphate dehydratase in plants: novel approach for structural and functional studies, and inhibitor blueprinting

The histidine biosynthetic pathway (HBP) is targeted for herbicide design with preliminary success only regarding imidazole-glycerol phosphate dehydratase (IGPD, EC 4.2.1.19), or HISN5, as referred to in plants. HISN5 catalyzes the sixth step of the HBP, in which imidazole-glycerol phosphate (IGP) is dehydrated to imidazole-acetol phosphate. In this work, we present high-resolution cryoEM and crystal structures of Medicago truncatula HISN5 (MtHISN5) in complexes with an inactive IGP diastereoisomer and with various other ligands. MtHISN5 can serve as a new model for plant HISN5 structural studies, as it enables resolving protein-ligand interactions at high (2.2 Å) resolution using cryoEM. We identified ligand-binding hotspots and characterized the features of plant HISN5 enzymes in the context of the HISN5-targeted inhibitor design. Virtual screening performed against millions of small molecules not only revealed candidate molecules but also identified linkers for fragments that were experimentally confirmed to bind. Based on experimental and computational approaches, this study provides guidelines for designing symmetric HISN5 inhibitors that can reach two neighboring active sites. Finally, we conducted analyses of sequence similarity networks revealing that plant HISN5 enzymes derive from cyanobacteria. We also adopted a new approach to measure MtHISN5 enzymatic activity using isothermal titration calorimetry and enzymatically synthesized IGP.

Computational Analysis of Non-synonymous SNPs in ATM Kinase: Structural Insights, Functional Implications, and Inhibitor Discovery

Ataxia telangiectasia-mutated (ATM) protein kinase, a key player in cellular integrity regulation, is known for its role in DNA damage response. This study investigates the broader impact of ATM on cellular processes and potential clinical manifestations arising from mutations, aiming to expand our understanding of ATM's diverse functions beyond conventional roles. The research employs a comprehensive set of computational techniques for a thorough analysis of ATM mutations. The mutation data are curated from dbSNP and HuVarBase databases. A meticulous assessment is conducted, considering factors such as deleterious effects, protein stability, oncogenic potential, and biophysical characteristics of the identified mutations. Conservation analysis, utilizing diverse computational tools, provides insights into the evolutionary significance of these mutations. Molecular docking and dynamic simulation analyses are carried out for selected mutations, investigating their interactions with Y2080D, AZD0156, and quercetin inhibitors to gauge potential therapeutic implications. Among the 419 mutations scrutinized, five (V1913C, Y2080D, L2656P, C2770G, and C2930G) are identified as both disease causing and protein destabilizing. The study reveals the oncogenic potential of these mutations, supported by findings from the COSMIC database. Notably, Y2080D is associated with haematopoietic and lymphoid cancers, while C2770G shows a correlation with squamous cell carcinomas. Molecular docking and dynamic simulation analyses highlight strong binding affinities of quercetin for Y2080D and AZD0156 for C2770G, suggesting potential therapeutic options. In summary, this computational analysis provides a comprehensive understanding of ATM mutations, revealing their potential implications in cellular integrity and cancer development. The study underscores the significance of Y2080D and C2770G mutations, offering valuable insights for future precision medicine targeting-specific ATM. Despite informative computational analyses, a significant research gap exists, necessitating essential in vitro and in vivo studies to validate the predicted effects of ATM mutations on protein structure and function.

Role of two modules controlling the interaction between SKAP1 and SRC kinases comparison with SKAP2 architecture and consequences for evolution

SRC kinase associated phosphoprotein 1 (SKAP1), an adaptor for protein assembly, plays an important role in the immune system such as stabilizing immune synapses. Understanding how these functions are controlled at the level of the protein-protein interactions is necessary to describe these processes and to develop therapeutics. Here, we dissected the SKAP1 modular organization to recognize SRC kinases and compared it to that of its paralog SRC kinase associated phosphoprotein 2 (SKAP2). Different conserved motifs common to either both proteins or specific to SKAP2 were found using this comparison. Two modules harboring different binding properties between SKAP1 and SKAP2 were identified: one composed of two conserved motifs located in the second interdomain interacting at least with the SH2 domain of SRC kinases and a second one composed of the DIM domain modulated by the SH3 domain and the activation of SRC kinases. This work suggests a convergent evolution of the binding properties of some SRC kinases interacting specifically with either SKAP1 or SKAP2.

ZP2 cleavage blocks polyspermy by modulating the architecture of the egg coat

Following the fertilization of an egg by a single sperm, the egg coat or zona pellucida (ZP) hardens and polyspermy is irreversibly blocked. These events are associated with the cleavage of the N-terminal region (NTR) of glycoprotein ZP2, a major subunit of ZP filaments. ZP2 processing is thought to inactivate sperm binding to the ZP, but its molecular consequences and connection with ZP hardening are unknown. Biochemical and structural studies show that cleavage of ZP2 triggers its oligomerization. Moreover, the structure of a native vertebrate egg coat filament, combined with AlphaFold predictions of human ZP polymers, reveals that two protofilaments consisting of type I (ZP3) and type II (ZP1/ZP2/ZP4) components interlock into a left-handed double helix from which the NTRs of type II subunits protrude. Together, these data suggest that oligomerization of cleaved ZP2 NTRs extensively cross-links ZP filaments, rigidifying the egg coat and making it physically impenetrable to sperm.

Structural basis for Retriever-SNX17 assembly and endosomal sorting

During endosomal recycling, Sorting Nexin 17 (SNX17) facilitates the transport of numerous membrane cargo proteins by tethering them to the Retriever complex. Despite its importance, the mechanisms underlying this interaction have remained elusive. Here, we report the structure of the Retriever-SNX17 complex determined using cryogenic electron microscopy (cryo-EM). Our structure reveals that the C-terminal tail of SNX17 engages with a highly conserved interface between the VPS35L and VPS26C subunits of Retriever. Through comprehensive biochemical, cellular, and proteomic analyses, we demonstrate that disrupting this interface impairs the Retriever-SNX17 interaction, subsequently affecting the recycling of SNX17-dependent cargos and altering the composition of the plasma membrane proteome. Intriguingly, we find that the SNX17-binding pocket on Retriever can be utilized by other ligands that share a consensus acidic C-terminal tail motif. By showing how SNX17 is linked to Retriever, our findings uncover a fundamental mechanism underlying endosomal trafficking of critical cargo proteins and reveal a mechanism by which Retriever can engage with other regulatory factors.

Distinct characteristics of putative archaeal 5-methylcytosine RNA methyltransferases unveil their substrate specificities and evolutionary ancestries

5-Methylcytosine methyltransferases (m5C MTases) are known to be involved in the modification of RNA. Although these enzymes have been relatively well characterized in bacteria and eukarya, a complete understanding of the archaeal counterparts is lacking. In this study, the identification and characterization of archaeal RNA m5C MTases were performed. As a case study, a hyperthermophilic archaeon, Pyrococcus horikoshii OT3, which possesses five putative RNA m5C MTases, was chosen. Among the five putative RNA m5C MTases, two proteins (PH0851 and PH1991) have been characterized as homologs of a bacterial rRNA MTase (RsmB) and eukaryal tRNA MTase (NSUN6), respectively. The in-depth characterization of the remaining three putative RNA m5C MTases (PH1078, PH1374, and PH1537) in this study suggests the presence of the signature architecture and catalytic residues plausibly involved in the binding of their cognate RNA substrates. Additionally, the results also suggest the existence of two RsmB-like proteins (PH0851 and PH1078) belonging to the same subfamily IV of m5C RNA MTase. However, the proteins PH1374 and PH1537 belong to the same subfamily V but bind to different substrates, rRNA and tRNA, respectively. The findings further indicate that archaeal RNA m5C MTases link those from bacteria and eukarya.Communicated by Ramaswamy H. Sarma.

Sec7 regulatory domains scaffold autoinhibited and active conformations

The late stages of Golgi maturation involve a series of sequential trafficking events in which cargo-laden vesicles are produced and targeted to multiple distinct subcellular destinations. Each of these vesicle biogenesis events requires activation of an Arf GTPase by the Sec7/BIG guanine nucleotide exchange factor (GEF). Sec7 localization and activity is regulated by autoinhibition, positive feedback, and interaction with other GTPases. Although these mechanisms have been characterized biochemically, we lack a clear picture of how GEF localization and activity is modulated by these signals. Here, we report the cryogenic electron microscopy structure of full-length Sec7 in its autoinhibited form, revealing the architecture of its multiple regulatory domains. We use functional experiments to determine the basis for autoinhibition and use structural predictions to produce a model for an active conformation of the GEF that is supported empirically. This study therefore elucidates the conformational transition that Sec7 undergoes to become active on the organelle membrane surface.

Structural Requirements for Activity of Mind bomb1 in Notch Signaling

Mind bomb 1 (MIB1) is a RING E3 ligase that ubiquitinates Notch ligands, a necessary step for induction of Notch signaling. The structural basis for binding of the JAG1 ligand by the N-terminal region of MIB1 is known, yet how the ankyrin (ANK) and RING domains of MIB1 cooperate to catalyze ubiquitin transfer from E2~Ub to Notch ligands remains unclear. Here, we show that the third RING domain and adjacent coiled coil region of MIB1 (ccRING3) drives MIB1 dimerization and that ubiquitin transfer activity of MIB1 relies solely on RING3. We report x-ray crystal structures of a UbcH5B-ccRING3 complex as a fusion protein and of the ANK region. Directly tethering the N-terminal region to ccRING3 forms a minimal MIB1 protein, which is sufficient to induce a Notch response in receiver cells. Together, these studies define the functional elements of an E3 ligase needed for ligands to induce a Notch signaling response.

SMARCAL1 ubiquitylation controls its association with RPA-coated ssDNA and promotes replication fork stability

Impediments in replication fork progression cause genomic instability, mutagenesis, and severe pathologies. At stalled forks, RPA-coated single-stranded DNA (ssDNA) activates the ATR kinase and directs fork remodeling, 2 key early events of the replication stress response. RFWD3, a recently described Fanconi anemia (FA) ubiquitin ligase, associates with RPA and promotes its ubiquitylation, facilitating late steps of homologous recombination (HR). Intriguingly, RFWD3 also regulates fork progression, restart and stability via poorly understood mechanisms. Here, we used proteomics to identify putative RFWD3 substrates during replication stress in human cells. We show that RFWD3 interacts with and ubiquitylates the SMARCAL1 DNA translocase directly in vitro and following DNA damage in vivo. SMARCAL1 ubiquitylation does not trigger its subsequent proteasomal degradation but instead disengages it from RPA thereby regulating its function at replication forks. Proper regulation of SMARCAL1 by RFWD3 at stalled forks protects them from excessive MUS81-mediated cleavage in response to UV irradiation, thereby limiting DNA replication stress. Collectively, our results identify RFWD3-mediated SMARCAL1 ubiquitylation as a novel mechanism that modulates fork remodeling to avoid genome instability triggered by aberrant fork processing.

The crystal structure of mycothiol disulfide reductase (Mtr) provides mechanistic insight into the specific low-molecular-weight thiol reductase activity of Actinobacteria

Low-molecular-weight (LMW) thiols are involved in many processes in all organisms, playing a protective role against reactive species, heavy metals, toxins and antibiotics. Actinobacteria, such as Mycobacterium tuberculosis, use the LMW thiol mycothiol (MSH) to buffer the intracellular redox environment. The NADPH-dependent FAD-containing oxidoreductase mycothiol disulfide reductase (Mtr) is known to reduce oxidized mycothiol disulfide (MSSM) to MSH, which is crucial to maintain the cellular redox balance. In this work, the first crystal structures of Mtr are presented, expanding the structural knowledge and understanding of LMW thiol reductases. The structural analyses and docking calculations provide insight into the nature of Mtrs, with regard to the binding and reduction of the MSSM substrate, in the context of related oxidoreductases. The putative binding site for MSSM suggests a similar binding to that described for the homologous glutathione reductase and its respective substrate glutathione disulfide, but with distinct structural differences shaped to fit the bulkier MSSM substrate, assigning Mtrs as uniquely functioning reductases. As MSH has been acknowledged as an attractive antitubercular target, the structural findings presented in this work may contribute towards future antituberculosis drug development.

FAM81A is a postsynaptic protein that regulates the condensation of postsynaptic proteins via liquid-liquid phase separation

Proteome analyses of the postsynaptic density (PSD), a proteinaceous specialization beneath the postsynaptic membrane of excitatory synapses, have identified several thousands of proteins. While proteins with predictable functions have been well studied, functionally uncharacterized proteins are mostly overlooked. In this study, we conducted a comprehensive meta-analysis of 35 PSD proteome datasets, encompassing a total of 5,869 proteins. Employing a ranking methodology, we identified 97 proteins that remain inadequately characterized. From this selection, we focused our detailed analysis on the highest-ranked protein, FAM81A. FAM81A interacts with PSD proteins, including PSD-95, SynGAP, and NMDA receptors, and promotes liquid-liquid phase separation of those proteins in cultured cells or in vitro. Down-regulation of FAM81A in cultured neurons causes a decrease in the size of PSD-95 puncta and the frequency of neuronal firing. Our findings suggest that FAM81A plays a crucial role in facilitating the interaction and assembly of proteins within the PSD, and its presence is important for maintaining normal synaptic function. Additionally, our methodology underscores the necessity for further characterization of numerous synaptic proteins that still lack comprehensive understanding.

Gossypium hirsutum gene of unknown function Gohir.A02G161000 encodes a potential transmembrane Root UVB Sensitive 4 Protein with a putative protein-protein interaction interface

A gene of unknown function, Gohir.A02G161000.1, identified in Gossypium hirsutum was studied using computational sequence and structure bioinformatics tools. The associated protein GhRUS4-A0A1U8JPV7 (UniProt A0A1U8JPV7) is predicted to be a plastid-localized, transmembrane root UVB-sensitive 4 (RUS4) protein with a newly identified potential dimerization surface. Evidence from homology and sequence conservation suggest involvement in auxin transport and pollen maturation.

The structural basis of eukaryotic chaperonin TRiC/CCT: Action and folding

Accurate folding of proteins in living cells often requires the cooperative support of molecular chaperones. Eukaryotic group II chaperonin Tailless complex polypeptide 1-Ring Complex (TRiC) accomplishes this task by providing a folding chamber for the substrate that is regulated by an Adenosine triphosphate (ATP) hydrolysis-dependent cycle. Once delivered to and recognized by TRiC, the nascent substrate enters the folding chamber and undergoes folding and release in a stepwise manner. During the process, TRiC subunits and cochaperones such as prefoldin and phosducin-like proteins interact with the substrate to assist the overall folding process in a substrate-specific manner. Coevolution between the components is supposed to consult the binding specificity and ultimately expand the substrate repertoire assisted by the chaperone network. This review describes the TRiC chaperonin and the substrate folding process guided by the TRiC network in cooperation with cochaperones, specifically focusing on recent progress in structural analyses.

DNMT3B PWWP mutations cause hypermethylation of heterochromatin

The correct establishment of DNA methylation patterns is vital for mammalian development and is achieved by the de novo DNA methyltransferases DNMT3A and DNMT3B. DNMT3B localises to H3K36me3 at actively transcribing gene bodies via its PWWP domain. It also functions at heterochromatin through an unknown recruitment mechanism. Here, we find that knockout of DNMT3B causes loss of methylation predominantly at H3K9me3-marked heterochromatin and that DNMT3B PWWP domain mutations or deletion result in striking increases of methylation in H3K9me3-marked heterochromatin. Removal of the N-terminal region of DNMT3B affects its ability to methylate H3K9me3-marked regions. This region of DNMT3B directly interacts with HP1α and facilitates the bridging of DNMT3B with H3K9me3-marked nucleosomes in vitro. Our results suggest that DNMT3B is recruited to H3K9me3-marked heterochromatin in a PWWP-independent manner that is facilitated by the protein's N-terminal region through an interaction with a key heterochromatin protein. More generally, we suggest that DNMT3B plays a role in DNA methylation homeostasis at heterochromatin, a process which is disrupted in cancer, aging and Immunodeficiency, Centromeric Instability and Facial Anomalies (ICF) syndrome.

Regulatory mechanisms of dopamine metabolism in a marine Meyerozyma guilliermondii GXDK6 under NaCl stress as revealed by integrative multi-omics analysis

Dopamine can be used to treat depression, myocardial infarction, and other diseases. However, few reports are available on the de novo microbial synthesis of dopamine from low-cost substrate. In this study, integrated omics technology was used to explore the dopamine metabolism of a novel marine multi-stress-tolerant aromatic yeast Meyerozyma guilliermondii GXDK6. GXDK6 was found to have the ability to biosynthesize dopamine when using glucose as the substrate. 14 key genes for the biosynthesis of dopamine were identified by whole genome-wide analysis. Transcriptomic and proteomic data showed that the expression levels of gene AAT2 encoding aspartate aminotransferase (regulating dopamine anabolism) were upregulated, while gene AO-I encoding copper amine oxidase (involved in dopamine catabolism) were downregulated under 10 % NaCl stress compared with non-NaCl stress, thereby contributing to biosynthesis of dopamine. Further, the amount of dopamine under 10 % NaCl stress was 2.51-fold higher than that of zero NaCl, which was consistent with the multi-omics results. Real-time fluorescence quantitative PCR (RT-qPCR) and high-performance liquid chromatography (HPLC) results confirmed the metabolic model of dopamine. Furthermore, by overexpressing AAT2, AST enzyme activity was increased by 24.89 %, the expression of genes related to dopamine metabolism was enhanced, and dopamine production was increased by 56.36 % in recombinant GXDK6AAT2. In conclusion, Meyerozyma guilliermondii GXDK6 could utilize low-cost carbon source to synthesize dopamine, and NaCl stress promoted the biosynthesis of dopamine.

Structural studies with crotoxin B from Crotalus durissus collilineatus venom suggest a heterodimeric assembly formed by two new isoforms

In accidents involving Crotalus snakes, the crotoxin complex (CTX) plays lethal action due to its neurotoxic activity. On the other hand, CTX have potential biotechnological application due to its anti-tumoral, anti-inflammatory, antimicrobial, analgesic and immunomodulatory properties. CTX is a heterodimer composed of Crotoxin A (CA or crotapotin), the acidic nontoxic and non-enzymatic component and; Crotoxin B (CB), a basic, toxic and catalytic PLA2. Currently, there are two classes of CTX isoforms, whose differences in their biological activities have been attributed to features presented in CB isoforms. Here, we present the crystal structure of CB isolated from the Crotalus durissus collilineatus venom. It amino acid sequence was assigned using the SEQUENCE SLIDER software, which revealed that the crystal structure is a heterodimer composed of two new CB isoforms (colCB-A and colCB-B). Bioinformatic and biophysical analyses showed that the toxin forms a tetrameric assembly in solution similar to CB from Crotalus durissus terrificus venom, despite some differences observed at the dimeric interface. By the previously proposed classification, the colCB-B presents features of the class I isoforms while colCB-A cannot be classified into classes I and II based on its amino acid sequence. Due to similar features observed for other CB isoforms found in the NCBI database and the results obtained for colCB-A, we suggest that there are more than two classes of CTX and CB isoforms in crotalic venoms.

Regulatory sites of CaM-sensitive adenylyl cyclase AC8 revealed by cryo-EM and structural proteomics

Membrane adenylyl cyclase AC8 is regulated by G proteins and calmodulin (CaM), mediating the crosstalk between the cAMP pathway and Ca2+ signalling. Despite the importance of AC8 in physiology, the structural basis of its regulation by G proteins and CaM is not well defined. Here, we report the 3.5 Å resolution cryo-EM structure of the bovine AC8 bound to the stimulatory Gαs protein in the presence of Ca2+/CaM. The structure reveals the architecture of the ordered AC8 domains bound to Gαs and the small molecule activator forskolin. The extracellular surface of AC8 features a negatively charged pocket, a potential site for unknown interactors. Despite the well-resolved forskolin density, the captured state of AC8 does not favour tight nucleotide binding. The structural proteomics approaches, limited proteolysis and crosslinking mass spectrometry (LiP-MS and XL-MS), allowed us to identify the contact sites between AC8 and its regulators, CaM, Gαs, and Gβγ, as well as to infer the conformational changes induced by these interactions. Our results provide a framework for understanding the role of flexible regions in the mechanism of AC regulation.

Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Bioinformatic analysis of human ZPR1 gene pathogenic exome mutations

Advanced sequencing technologies enable rapid detection of sequence variants, aiming to uncover the molecular foundations of human genetic disorders. The challenge lies in interpreting the influence of new exome variants that lead to diverse phenotypes. Our study introduces a detailed, multi-tiered method for assessing the impact of novel variants, particularly focusing on the zinc finger protein 1 (ZPR1) gene. Herein, we employed a combination of variant effect predictors, protein stability analyses, and the American College of Medical Genetics and Association of Molecular Pathology (ACMG/AMP) guidelines. Our structural analysis pinpoints specific amino acid residues in the ZPR1 zinc finger domains that are sensitive to changes, distinguishing between benign and disease-causing coding variants using rigorous in silico tools. We examined 223 germline ZPR1 exome variants, uncovering significant ethnic disparities in the frequency of heterozygous harmful ZPR1 variants, ranging from 0.04% in the Ashkenazi Jewish population to 0.34% in African/African Americans. Additionally, the discovery of three homozygous carriers in European and South Asian groups suggests a higher occurrence of ZPR1 variants in these demographics, meriting further exploration. This research provides insights into the prevalence and implications of amino acid substitutions in the ZPR1 protein.

Structural-functional analysis of drug target aspartate semialdehyde dehydrogenase

Aspartate β-semialdehyde dehydrogenase (ASADH) is a key enzyme in the biosynthesis of essential amino acids in microorganisms and some plants. Inhibition of ASADHs can be a potential drug target for developing novel antimicrobial and herbicidal compounds. This review covers up-to-date information about sequence diversity, ligand/inhibitor-bound 3D structures, potential inhibitors, and key pharmacophoric features of ASADH useful in designing novel and target-specific inhibitors of ASADH. Most reported ASADH inhibitors have two highly electronegative functional groups that interact with two key arginyl residues present in the active site of ASADHs. The structural information, active site binding modes, and key interactions between the enzyme and inhibitors serve as the basis for designing new and potent inhibitors against the ASADH family.

ATG9B is a tissue-specific homotrimeric lipid scramblase that can compensate for ATG9A

Macroautophagy/autophagy is a fundamental aspect of eukaryotic biology, and the autophagy-related protein ATG9A is part of the core machinery facilitating this process. In addition to ATG9A vertebrates encode ATG9B, a poorly characterized paralog expressed in a subset of tissues. Herein, we characterize the structure of human ATG9B revealing the conserved homotrimeric quaternary structure and explore the conformational dynamics of the protein. Consistent with the experimental structure and computational chemistry, we establish that ATG9B is a functional lipid scramblase. We show that ATG9B can compensate for the absence of ATG9A in starvation-induced autophagy displaying similar subcellular trafficking and steady-state localization. Finally, we demonstrate that ATG9B can form a heteromeric complex with ATG2A. By establishing the molecular structure and function of ATG9B, our results inform the exploration of niche roles for autophagy machinery in more complex eukaryotes and reveal insights relevant across species.Abbreviation: ATG: autophagy related; CHS: cholesteryl hemisuccinate; cryo-EM: single-particle cryogenic electron microscopy; CTF: contrast transfer function: CTH: C- terminal α helix; FSC: fourier shell correlation; HDIR: HORMA domain interacting region; LMNG: lauryl maltose neopentyl glycol; MD: molecular dynamics simulations; MSA: multiple sequence alignment; NBD-PE: 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(7-nitro-2-1,3-benzoxadiazol-4-yl ammonium salt); POPC: palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine; RBG: repeating beta groove domain; RMSD: root mean square deviation; SEC: size-exclusion chromatography; TMH: transmembrane helix.

In silico analysis of alpha-synuclein protein variants and posttranslational modifications related to Parkinson's disease

Parkinson's disease (PD) is among the most prevalent neurodegenerative disorders, affecting over 10 million people worldwide. The protein encoded by the SNCA gene, alpha-synuclein (ASYN), is the major component of Lewy body (LB) aggregates, a histopathological hallmark of PD. Mutations and posttranslational modifications (PTMs) in ASYN are known to influence protein aggregation and LB formation, possibly playing a crucial role in PD pathogenesis. In this work, we applied computational methods to characterize the effects of missense mutations and PTMs on the structure and function of ASYN. Missense mutations in ASYN were compiled from the literature/databases and underwent a comprehensive predictive analysis. Phosphorylation and SUMOylation sites of ASYN were retrieved from databases and predicted by algorithms. ConSurf was used to estimate the evolutionary conservation of ASYN amino acids. Molecular dynamics (MD) simulations of ASYN wild-type and variants A30G, A30P, A53T, and G51D were performed using the GROMACS package. Seventy-seven missense mutations in ASYN were compiled. Although most mutations were not predicted to affect ASYN stability, aggregation propensity, amyloid formation, and chaperone binding, the analyzed mutations received relatively high rates of deleterious predictions and predominantly occurred at evolutionarily conserved sites within the protein. Moreover, our predictive analyses suggested that the following mutations may be possibly harmful to ASYN and, consequently, potential targets for future investigation: K6N, T22I, K34E, G36R, G36S, V37F, L38P, G41D, and K102E. The MD analyses pointed to remarkable flexibility and essential dynamics alterations at nearly all domains of the studied variants, which could lead to impaired contact between NAC and the C-terminal domain triggering protein aggregation. These alterations may have functional implications for ASYN and provide important insight into the molecular mechanism of PD, supporting the design of future biomedical research and improvements in existing therapies for the disease.

Structural dissection of two redox proteins from the shipworm symbiont Teredinibacter turnerae

The discovery of lytic polysaccharide monooxygenases (LPMOs), a family of copper-dependent enzymes that play a major role in polysaccharide degradation, has revealed the importance of oxidoreductases in the biological utilization of biomass. In fungi, a range of redox proteins have been implicated as working in harness with LPMOs to bring about polysaccharide oxidation. In bacteria, less is known about the interplay between redox proteins and LPMOs, or how the interaction between the two contributes to polysaccharide degradation. We therefore set out to characterize two previously unstudied proteins from the shipworm symbiont Teredinibacter turnerae that were initially identified by the presence of carbohydrate binding domains appended to uncharacterized domains with probable redox functions. Here, X-ray crystal structures of several domains from these proteins are presented together with initial efforts to characterize their functions. The analysis suggests that the target proteins are unlikely to function as LPMO electron donors, raising new questions as to the potential redox functions that these large extracellular multi-haem-containing c-type cytochromes may perform in these bacteria.

Structure of glyceraldehyde-3-phosphate dehydrogenase from Paracoccidioides lutzii in complex with an aldonic sugar acid

The pathogen Paracoccidioides lutzii (Pb01) is found in South America countries Colombia, Ecuador, Venezuela and Brazil, especially in the central, west, and north regions of the latter. It belongs to the Ajellomycetaceae family, Onygenales order, and is typically thermodimorphic, presenting yeast cells when it grows in animal tissues, but mycelia when in the environment, where it produces the infectious propagule. This fungus is one of the etiologic agents of Paracoccidioidomycosis (PCM), the most important endemic fungal infection in Latin America. Investigations on its genome have contributed to a better understanding about its metabolism and revealed the complexity of several metabolic glycolytic pathways. Glyceraldehyde-3-Phosphate Dehydrogenase from Paracoccidioides lutzii (PlGAPDH) is considered a moonlighting protein and participates in several biological processes of this pathogen. The enzyme was expressed and purified, as seen in SDS-PAGE gel, crystallized and had its three dimensional structure (3D) determined in complex with NAD+, a sulphate ion and d-galactonic acid, therefore, a type of 'GAA site'. It is the first GAPDH structure to show this chemical type in this site and how this protein can bind an acid derived from oxidation of a linear hexose.

The binding of Mint/X11 PDZ domains to Ca V 2 calcium channels predates bilaterian animals

PDZ domain mediated interactions with voltage-gated calcium (Ca V ) channel C-termini play important roles in localizing and compartmentalizing membrane Ca 2+ signaling. The first such interaction discovered was between the neuronal multi-domain protein Mint-1, and the presynaptc calcium channel Ca V 2.2 in mammals. Although the physiological significance of this interaction is unclear, its occurrence in vertebrates and bilaterian invertebrates suggests important and conserved functions. In this study, we explore the evolutionary origins of Mint and its interaction with Ca V 2 channels. Phylogenetic and structural in silico analyses revealed that Mint is an animal-specific gene, like Ca V 2 channels, which bears a highly divergent N-terminus but strongly conserved C-terminus comprised of a phosphotyrosine binding domain, two tandem PDZ domains (PDZ-1 and PDZ-2), and a C-terminal auto-inhibitory element that binds and inhibits PDZ-1. Also deeply conserved are other Mint interacting proteins, namely amyloid precursor and related proteins, presenilins, neurexin, as well as CASK and Veli which form a tripartite complex with Mint in bilaterians. Through yeast 2-hybrid and bacterial 2-hybrid experiments, we show that Mint and Ca V 2 channels from cnidarians and placozoans interact in vitro , and in situ hybridization revealed co-expression of corresponding transcripts in dissociated neurons from the cnidarian Nematostella vectensis . Unexpectedly, the Mint orthologue from the ctenophore Hormiphora californiensis was able to strongly bind the divergent C-terminal ligands of cnidarian and placozoan Ca V 2 channels, despite neither the ctenophore Mint, nor the placozoan and cnidarian orthologues, binding the ctenophore Ca V 2 channel C-terminus. Altogether, our analyses provide a model for the emergence of this interaction in early animals first via adoption of a PDZ ligand by Ca V 2 channels, followed by sequence changes in the ligand that caused a modality switch for binding to Mint.

Flexible structural arrangement and DNA-binding properties of protein p6 from Bacillus subtillis phage φ29

The genome-organizing protein p6 of Bacillus subtilis bacteriophage φ29 plays an essential role in viral development by activating the initiation of DNA replication and participating in the early-to-late transcriptional switch. These activities require the formation of a nucleoprotein complex in which the DNA adopts a right-handed superhelix wrapping around a multimeric p6 scaffold, restraining positive supercoiling and compacting the viral genome. Due to the absence of homologous structures, prior attempts to unveil p6's structural architecture failed. Here, we employed AlphaFold2 to engineer rational p6 constructs yielding crystals for three-dimensional structure determination. Our findings reveal a novel fold adopted by p6 that sheds light on its self-association mechanism and its interaction with DNA. By means of protein-DNA docking and molecular dynamic simulations, we have generated a comprehensive structural model for the nucleoprotein complex that consistently aligns with its established biochemical and thermodynamic parameters. Besides, through analytical ultracentrifugation, we have confirmed the hydrodynamic properties of the nucleocomplex, further validating in solution our proposed model. Importantly, the disclosed structure not only provides a highly accurate explanation for previously experimental data accumulated over decades, but also enhances our holistic understanding of the structural and functional attributes of protein p6 during φ29 infection.

Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease

The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.Communicated by Ramaswamy H. Sarma.

Cancer-Related Mutations in the Sam Domains of EphA2 Receptor and Ship2 Lipid Phosphatase: A Computational Study

The lipid phosphatase Ship2 interacts with the EphA2 receptor by forming a heterotypic Sam (sterile alpha motif)-Sam complex. Ship2 works as a negative regulator of receptor endocytosis and consequent degradation, and anti-oncogenic effects in cancer cells should be induced by hindering its association with EphA2. Herein, a computational approach is presented to investigate the relationship between Ship2-Sam/EphA2-Sam interaction and cancer onset and further progression. A search was first conducted through the COSMIC (Catalogue of Somatic Mutations in Cancer) database to identify cancer-related missense mutations positioned inside or close to the EphA2-Sam and Ship2-Sam reciprocal binding interfaces. Next, potential differences in the chemical-physical properties of mutant and wild-type Sam domains were evaluated by bioinformatics tools based on analyses of primary sequences. Three-dimensional (3D) structural models of mutated EphA2-Sam and Ship2-Sam domains were built as well and deeply analysed with diverse computational instruments, including molecular dynamics, to classify potentially stabilizing and destabilizing mutations. In the end, the influence of mutations on the EphA2-Sam/Ship2-Sam interaction was studied through docking techniques. This in silico approach contributes to understanding, at the molecular level, the mutation/cancer relationship by predicting if amino acid substitutions could modulate EphA2 receptor endocytosis.

Structural and Functional Characterization of Obesumbacterium proteus Phytase: A Comprehensive In-Silico Study

Phytate, also known as myoinositol hexakisphosphate, exhibits anti-nutritional properties and possesses a negative environmental impact. Phytase enzymes break down phytate, showing potential in various industries, necessitating thorough biochemical and computational characterizations. The present study focuses on Obesumbacterium proteus phytase (OPP), indicating its similarities with known phytases and its potential through computational analyses. Structure, functional, and docking results shed light on OPP's features, structural stability, strong and stable interaction, and dynamic conformation, with flexible sidechains that could adapt to different temperatures or specific functions. Root Mean Square fluctuation (RMSF) highlighted fluctuating regions in OPP, indicating potential sites for stability enhancement through mutagenesis. The systematic approach developed here could aid in enhancing enzyme properties via a rational engineering approach. Computational analysis expedites enzyme discovery and engineering, complementing the traditional biochemical methods to accelerate the quest for superior enzymes for industrial applications.