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

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

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

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

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

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

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

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

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

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

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

Structural organization of the retriever-CCC endosomal recycling complex

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

Fructose-1-kinase has pleiotropic roles in Escherichia coli

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

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

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

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

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

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

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

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

BACKGROUND

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

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

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

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

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

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

High-throughput crystallographic fragment screening of Zika virus NS3 Helicase

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

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

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

Cyclodipeptide oxidase is an enzyme filament

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

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

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

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

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

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

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

Molecular Dynamics Simulation of Kir6.2 Variants Reveals Potential Association with Diabetes Mellitus

Diabetes mellitus (DM) represents a problem for the healthcare system worldwide. DM has very serious complications such as blindness, kidney failure, and cardiovascular disease. In addition to the very bad socioeconomic impacts, it influences patients and their families and communities. The global costs of DM and its complications are huge and expected to rise by the year 2030. DM is caused by genetic and environmental risk factors. Genetic testing will aid in early diagnosis and identification of susceptible individuals or populations using ATP-sensitive potassium (KATP) channels present in different tissues such as the pancreas, myocardium, myocytes, and nervous tissues. The channels respond to different concentrations of blood sugar, stimulation by hormones, or ischemic conditions. In pancreatic cells, they regulate the secretion of insulin and glucagon. Mutations in the KCNJ11 gene that encodes the Kir6.2 protein (a major constituent of KATP channels) were reported to be associated with Type 2 DM, neonatal diabetes mellitus (NDM), and maturity-onset diabetes of the young (MODY). Kir6.2 harbors binding sites for ATP and phosphatidylinositol 4,5-diphosphate (PIP2). The ATP inhibits the KATP channel, while the (PIP2) activates it. A Kir6.2 mutation at tyrosine330 (Y330) was demonstrated to reduce ATP inhibition and predisposes to NDM. In this study, we examined the effect of mutations on the Kir6.2 structure using bioinformatics tools and molecular dynamic simulations (SIFT, PolyPhen, SNAP2, PANTHER, PhD&SNP, SNP&Go, I-Mutant, MuPro, MutPred, ConSurf, HOPE, and GROMACS). Our results indicated that M199R, R201H, R206H, and Y330H mutations influence Kir6.2 structure and function and therefore may cause DM. We conclude that MD simulations are useful techniques to predict the effects of mutations on protein structure. In addition, the M199R, R201H, R206H, and Y330H variant in the Kir6.2 protein may be associated with DM. These results require further verification in protein-protein interactions, Kir6.2 function, and case-control studies.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Choose Your Own Adventure: A Comprehensive Database of Reactions Catalyzed by Cytochrome P450 BM3 Variants

Cytochrome P450 BM3 monooxygenase is the topic of extensive research as many researchers have evolved this enzyme to generate a variety of products. However, the abundance of information on increasingly diversified variants of P450 BM3 that catalyze a broad array of chemistry is not in a format that enables easy extraction and interpretation. We present a database that categorizes variants by their catalyzed reactions and includes details about substrates to provide reaction context. This database of >1500 P450 BM3 variants is downloadable and machine-readable and includes instructions to maximize ease of gathering information. The database allows rapid identification of commonly reported substitutions, aiding researchers who are unfamiliar with the enzyme in identifying starting points for enzyme engineering. For those actively engaged in engineering P450 BM3, the database, along with this review, provides a powerful and user-friendly platform to understand, predict, and identify the attributes of P450 BM3 variants, encouraging the further engineering of this enzyme.

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

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

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

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

ANN based prediction of ligand binding sites outside deep cavities to facilitate drug designing

The ever-changing environmental conditions and pollution are the prime reasons for the onset of several emerging and re-merging diseases. This demands the faster designing of new drugs to curb the deadly diseases in less waiting time to cure the animals and humans. Drug molecules interact with only protein surface on specific locations termed as ligand binding sites (LBS). Therefore, the knowledge of LBS is required for rational drug designing. Existing geometrical LBS prediction methods rely on search of cavities based on the fact that 83% of the LBS found in deep cavities, however, these methods usually fail where LBS localize outside deep cavities. To overcome this challenge, the present work provides an artificial neural network (ANN) based method to predict LBS outside deep cavities in animal proteins including human to facilitate drug designing. In the present work a feed-forward backpropagation neural network was trained by utilizing 38 structural, atomic, physiochemical, and evolutionary discriminant features of LBS and non-LBS residues localized in the extracted roughest patch on protein surface. The performance of this ANN based prediction method was found 76% better for those proteins where cavity subspace (extracted by MetaPocket 2.0, a consensus method) failed to predict LBS due to their localization outside the deep cavities. The prediction of LBS outside deep cavities will facilitate in drug designing for the proteins where it is not possible due to lack of LBS information as the geometrical LBS prediction methods rely on extraction of deep cavities.

Proteins à la carte: riboproteogenomic exploration of bacterial N-terminal proteoform expression

In recent years, it has become evident that the true complexity of bacterial proteomes remains underestimated. Gene annotation tools are known to propagate biases and overlook certain classes of truly expressed proteins, particularly proteoforms-protein isoforms arising from a single gene. Recent (re-)annotation efforts heavily rely on ribosome profiling by providing a direct readout of translation to fully describe bacterial proteomes. In this study, we employ a robust riboproteogenomic pipeline to conduct a systematic census of expressed N-terminal proteoform pairs, representing two isoforms encoded by a single gene raised by annotated and alternative translation initiation, in Salmonella. Intriguingly, conditional-dependent changes in relative utilization of annotated and alternative translation initiation sites (TIS) were observed in several cases. This suggests that TIS selection is subject to regulatory control, adding yet another layer of complexity to our understanding of bacterial proteomes.

IMPORTANCE

With the emerging theme of genes within genes comprising the existence of alternative open reading frames (ORFs) generated by translation initiation at in-frame start codons, mechanisms that control the relative utilization of annotated and alternative TIS need to be unraveled and our molecular understanding of resulting proteoforms broadened. Utilizing complementary ribosome profiling strategies to map ORF boundaries, we uncovered dual-encoding ORFs generated by in-frame TIS usage in Salmonella. Besides demonstrating that alternative TIS usage may generate proteoforms with different characteristics, such as differential localization and specialized function, quantitative aspects of conditional retapamulin-assisted ribosome profiling (Ribo-RET) translation initiation maps offer unprecedented insights into the relative utilization of annotated and alternative TIS, enabling the exploration of gene regulatory mechanisms that control TIS usage and, consequently, the translation of N-terminal proteoform pairs.

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

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

Deciphering the Impact of Rare Missense Variants in EGFR-TKI-Resistant Non-Small-Cell Lung Cancer through Whole Exome Sequencing: A Computational Approach

Targeted therapy revolutionizes the treatment of non-small-cell lung cancer (NSCLC), harboring molecular change. Epidermal growth factor receptor(EGFR) mutations play a crucial role in the development of NSCLC, serving as a pivotal factor in its pathogenesis. We elucidated the mechanisms of resistance and potential therapeutic strategies in NSCLC resistant to the EGFR-tyrosine kinase inhibitor (EGFR-TKI). This is achieved by identifying rare missense variants through whole exome sequencing (WES). The goal is to enhance our understanding, identify biomarkers, and lay the groundwork for targeted interventions, thereby offering hope for an improved NSCLC treatment landscape. We conducted WES analysis on 16 NSCLC samples with EGFR-TKI-resistant NSCLC obtained from SRA-NCBI (PRJEB50602) to reveal genomic profiles within the EGFR-TKI. Our findings showed that 48% of the variants were missense, and after filtering with the Ensembl variant effect predictor, 53 rare missense variants in 23 genes were identified as highly deleterious. Further examination using pathogenic tools like PredictSNP revealed 12 deleterious rare missense variants in 7 genes: ZNF717, PSPH, ESRRA, SEMA3G, PTPN7, CAVIN4, and MYBBP1A. Molecular dynamics simulation (MDS) suggested that the L385P variant alters the structural flexibility of ESRRA, potentially leading to unfolding of ERRα proteins. This could impact their function and alter ERRα expression. These insights from MDS enhance our understanding of the structural and dynamic consequences of the L385P ESRRA variant and provide valuable implications for subsequent therapeutic considerations and targeted interventions.

Emergence of an Orphan Nitrogenase Protein Following Atmospheric Oxygenation

Molecular innovations within key metabolisms can have profound impacts on element cycling and ecological distribution. Yet, much of the molecular foundations of early evolved enzymes and metabolisms are unknown. Here, we bring one such mystery to relief by probing the birth and evolution of the G-subunit protein, an integral component of certain members of the nitrogenase family, the only enzymes capable of biological nitrogen fixation. The G-subunit is a Paleoproterozoic-age orphan protein that appears more than 1 billion years after the origin of nitrogenases. We show that the G-subunit arose with novel nitrogenase metal dependence and the ecological expansion of nitrogen-fixing microbes following the transition in environmental metal availabilities and atmospheric oxygenation that began ∼2.5 billion years ago. We identify molecular features that suggest early G-subunit proteins mediated cofactor or protein interactions required for novel metal dependency, priming ancient nitrogenases and their hosts to exploit these newly diversified geochemical environments. We further examined the degree of functional specialization in G-subunit evolution with extant and ancestral homologs using laboratory reconstruction experiments. Our results indicate that permanent recruitment of the orphan protein depended on the prior establishment of conserved molecular features and showcase how contingent evolutionary novelties might shape ecologically important microbial innovations.

Lipoxygenases at the Intersection of Infection and Carcinogenesis

The persisting presence of opportunistic pathogens like Pseudomonas aeruginosa poses a significant threat to many immunocompromised cancer patients with pulmonary infections. This review highlights the complexity of interactions in the host's defensive eicosanoid signaling network and its hijacking by pathogenic bacteria to their own advantage. Human lipoxygenases (ALOXs) and their mouse counterparts are integral elements of the innate immune system, mostly operating in the pro-inflammatory mode. Taking into account the indispensable role of inflammation in carcinogenesis, lipoxygenases have counteracting roles in this process. In addition to describing the structure-function of lipoxygenases in this review, we discuss their roles in such critical processes as cancer cell signaling, metastases, death of cancer and immune cells through ferroptosis, as well as the roles of ALOXs in carcinogenesis promoted by pathogenic infections. Finally, we discuss perspectives of novel oncotherapeutic approaches to harness lipoxygenase signaling in tumors.

Refolding, Crystallization, and Crystal Structure Analysis of a Scavenger Receptor Cysteine-Rich Domain of Human Salivary Agglutinin Expressed in Escherichia coli

Scavenger receptors are a protein superfamily that typically consists of one or more repeats of the scavenger receptor cysteine-rich structural domain (SRCRD), which is an ancient and highly conserved protein module. The expression and purification of eukaryotic proteins containing multiple disulfide bonds has always been challenging. The expression systems that are commonly used to express SRCRD proteins mainly consist of eukaryotic protein expression systems. Herein, we established a high-level expression strategy of a Type B SRCRD unit from human salivary agglutinin using the Escherichia coli expression system, followed by a refolding and purification process. The untagged recombinant SRCRD was expressed in E. coli using the pET-32a vector, which was followed by a refolding process using the GSH/GSSG redox system. The SRCRD expressed in E. coli SHuffle T7 showed better solubility after refolding than that expressed in E. coli BL21(DE3), suggesting the importance of the disulfide bond content prior to refolding. The quality of the refolded protein was finally assessed using crystallization and crystal structure analysis. As proteins refolded from inclusion bodies exhibit a high crystal quality and reproducibility, this method is considered a reliable strategy for SRCRD protein expression and purification. To further confirm the structural integrity of the refolded SRCRD protein, the purified protein was subjected to crystallization using sitting-drop vapor diffusion method. The obtained crystals of SRCRD diffracted X-rays to a resolution of 1.47 Å. The solved crystal structure appeared to be highly conserved, with four disulfide bonds appropriately formed. The surface charge distribution of homologous SRCRD proteins indicates that the negatively charged region at the surface is associated with their calcium-dependent ligand recognition. These results suggest that a high-quality SRCRD protein expressed by E. coli SHuffle T7 can be successfully folded and purified, providing new options for the expression of members of the scavenger receptor superfamily.

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

Structural insights into the mechanism of protein transport by the Type 9 Secretion System translocon

Secretion systems are protein export machines that enable bacteria to exploit their environment through the release of protein effectors. The Type 9 Secretion System (T9SS) is responsible for protein export across the outer membrane (OM) of bacteria of the phylum Bacteroidota. Here we trap the T9SS of Flavobacterium johnsoniae in the process of substrate transport by disrupting the T9SS motor complex. Cryo-EM analysis of purified substrate-bound T9SS translocons reveals an extended translocon structure in which the previously described translocon core is augmented by a periplasmic structure incorporating the proteins SprE, PorD and a homologue of the canonical periplasmic chaperone Skp. Substrate proteins bind to the extracellular loops of a carrier protein within the translocon pore. As transport intermediates accumulate on the translocon when energetic input is removed, we deduce that release of the substrate-carrier protein complex from the translocon is the energy-requiring step in T9SS transport.

Rational evolution for altering the ligand preference of estrogen receptor alpha

Estrogen receptor α is commonly used in synthetic biology to control the activity of genome editing tools. The activating ligands, estrogens, however, interfere with various cellular processes, thereby limiting the applicability of this receptor. Altering its ligand preference to chemicals of choice solves this hurdle but requires adaptation of unspecified ligand-interacting residues. Here, we provide a solution by combining rational protein design with multi-site-directed mutagenesis and directed evolution of stably integrated variants in Saccharomyces cerevisiae. This method yielded an estrogen receptor variant, named TERRA, that lost its estrogen responsiveness and became activated by tamoxifen, an anti-estrogenic drug used for breast cancer treatment. This tamoxifen preference of TERRA was maintained in mammalian cells and mice, even when fused to Cre recombinase, expanding the mammalian synthetic biology toolbox. Not only is our platform transferable to engineer ligand preference of any steroid receptor, it can also profile drug-resistance landscapes for steroid receptor-targeted therapies.