Structural basis for ion selectivity in potassium-selective channelrhodopsins.

KCR channelrhodopsins (K+-selective light-gated ion channels) have received attention as potential inhibitory optogenetic tools but more broadly pose a fundamental mystery regarding how their K+ selectivity is achieved. Here, we present 2.5-2.7 Å cryo-electron microscopy structures of HcKCR1 and HcKCR2 and of a structure-guided mutant with enhanced K+ selectivity. Structural, electrophysiological, computational, spectroscopic, and biochemical analyses reveal a distinctive mechanism for K+ selectivity; rather than forming the symmetrical filter of canonical K+ channels achieving both selectivity and dehydration, instead, three extracellular-vestibule residues within each monomer form a flexible asymmetric selectivity gate, while a distinct dehydration pathway extends intracellularly. Structural comparisons reveal a retinal-binding pocket that induces retinal rotation (accounting for HcKCR1/HcKCR2 spectral differences), and design of corresponding KCR variants with increased K+ selectivity (KALI-1/KALI-2) provides key advantages for optogenetic inhibition in vitro and in vivo. Thus, discovery of a mechanism for ion-channel K+ selectivity also provides a framework for next-generation optogenetics.

Mechanism of Substrate Translocation in an Alternating Access Transporter.

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

Cryo-EM reveals mechanistic insights into lipid-facilitated polyamine export by human ATP13A2.

The cytoplasmic polyamine maintains cellular homeostasis by chelating toxic metal cations, regulating transcriptional activity, and protecting DNA. ATP13A2 was identified as a lysosomal polyamine exporter responsible for polyamine release into the cytosol, and its dysfunction is associated with Alzheimer's disease and other neural degradation diseases. ATP13A2 belongs to the P5 subfamily of the P-type ATPase family, but its mechanisms remain unknown. Here, we report the cryoelectron microscopy (cryo-EM) structures of human ATP13A2 under four different conditions, revealing the structural coupling between the polyamine binding and the dephosphorylation. Polyamine is bound at the luminal tunnel and recognized through numerous electrostatic and π-cation interactions, explaining its broad specificity. The unique N-terminal domain is anchored to the lipid membrane to stabilize the E2P conformation, thereby accelerating the E1P-to-E2P transition. These findings reveal the distinct mechanism of P5B ATPases, thereby paving the way for neuroprotective therapy by activating ATP13A2.

A Zn-dependent structural transition of SOD1 modulates its ability to undergo phase separation.

The misfolding and mutation of Cu/Zn superoxide dismutase (SOD1) is commonly associated with amyotrophic lateral sclerosis (ALS). SOD1 can accumulate within stress granules (SGs), a type of membraneless organelle, which is believed to form via liquid-liquid phase separation (LLPS). Using wild-type, metal-deficient, and different ALS disease mutants of SOD1 and computer simulations, we report here that the absence of Zn leads to structural disorder within two loop regions of SOD1, triggering SOD1 LLPS and amyloid formation. The addition of exogenous Zn to either metal-free SOD1 or to the severe ALS mutation I113T leads to the stabilization of the loops and impairs SOD1 LLPS and aggregation. Moreover, partial Zn-mediated inhibition of LLPS was observed for another severe ALS mutant, G85R, which shows perturbed Zn-binding. By contrast, the ALS mutant G37R, which shows reduced Cu-binding, does not undergo LLPS. In addition, SOD1 condensates induced by Zn-depletion exhibit greater cellular toxicity than aggregates formed by prolonged incubation under aggregating conditions. Overall, our work establishes a role for Zn-dependent modulation of SOD1 conformation and LLPS properties that may contribute to amyloid formation.

Folds from fold: Exploring topological isoforms of a single-domain protein.

Expanding the protein fold space beyond linear chains is of fundamental significance, yet remains largely unexplored. Herein, we report the creation of seven topological isoforms (i.e., linear, cyclic, knot, lasso, pseudorotaxane, and catenane) from a single protein fold precursor by rewiring the connectivity of secondary structure elements of the SpyTag-SpyCatcher complex and mutating the reactive residue on SpyTag to abolish the isopeptide bonding. These topological isoforms can be directly expressed in cells. Their topologies were confirmed by combined techniques of proteolytic digestion, fluorescence correlation spectroscopy (FCS), size-exclusion chromatography (SEC), and topological transformation. To study the effects of topology on their structures and properties, their biophysical properties were characterized by differential scanning calorimetry (DSC), heteronuclear single quantum coherence nuclear magnetic resonance spectroscopy (HSQC-NMR), and circular dichroism (CD) spectroscopy. Molecular dynamics (MD) simulations were further performed to reveal the atomic details of structural changes upon unfolding. Both experimental and simulation results suggest that they share a similar, well-folded hydrophobic core but exhibit distinct folding/unfolding dynamic behaviors. These results shed light onto the folding landscape of topological isoforms derived from the same protein fold. As a model system, this work improves our understanding of protein structure and dynamics beyond linear chains and suggests that protein folds are highly amenable to topological variation.

Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1.

De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.

Specific residues and conformational plasticity define the substrate specificity of short-chain dehydrogenases/reductases.

Short-chain dehydrogenases/reductases (SDRs) are one of the most prevalent enzyme families distributed among the sequenced microorganisms. Despite the presence of a conserved catalytic tetrad and high structural similarity, these enzymes exhibit different substrate specificities. The insufficient knowledge regarding the amino acids underlying substrate specificity hinders the understanding of the SDRs' roles in diverse and significant biological processes. Here, we performed bioinformatic analysis, molecular modeling, and mutagenesis studies to identify the key residues that regulate the substrate specificities of two homologous microbial SDRs (i.e., DesE and KduD). Further, we investigated the impact of altering the physicochemical properties of these amino acids on enzyme activity. Interestingly, molecular dynamics simulations also suggest a critical role of enzyme conformational flexibility in substrate recognition and catalysis. Overall, our findings improve the understanding of microbial SDR substrate specificity and shed light on future rational design of more efficient and effective biocatalysts.

High-Strength Polycrystalline Covalent Organic Framework with Abnormal Thermal Transport Insensitive to Grain Boundary.

Grain boundaries (GBs) in two-dimensional (2D) covalent organic frameworks (COFs) unavoidably form during the fabrication process, playing pivotal roles in the physical characteristics of COFs. Herein, molecular dynamics simulations were employed to elucidate the fracture failure and thermal transport mechanisms of polycrystalline COFs (p-COFs). The results revealed that the tilt angle of GBs significantly influences out-of-plane wrinkles and residual stress in monolayer p-COFs. The tensile strength of p-COFs can be enhanced and weakened with the tilt angle, which exhibits an inverse relationship with the defect density. The crack always originates from weaker heptagon rings during uniaxial tension. Notably, the thermal transport in p-COFs is insensitive to the GBs due to the variation of minor polymer chain length at defects, which is abnormal for other 2D crystalline materials. This study contributes insights into the impact of GBs in p-COFs and offers theoretical guidance for structural design and practical applications of advanced COFs.

Different conformational dynamics of SNARE protein Ykt6 among yeast and mammals.

Ykt6 is one of the most conserved SNARE (N-ethylmaleimide-sensitive factor attachment protein receptor) proteins involved in multiple intracellular membrane trafficking processes. The membrane-anchoring function of Ykt6 has been elucidated to result from its conformational transition from a closed state to an open state. Two ways of regulating the conformational transition were proposed: the C-terminal lipidation and the phosphorylation at the SNARE core. Despite many aspects of common properties, Ykt6 displays differential cellular localizations and functional behaviors in different species, such as yeast, mammals, and worms. The structure-function relationship underlying these differences remains elusive. Here, we combined biochemical characterization, single-molecule FRET measurement, and molecular dynamics simulation to compare the conformational dynamics of yeast and rat Ykt6. Compared to rat Ykt6 (rYkt6), yeast Ykt6 (yYkt6) has more open conformations and could not bind dodecylphosphocholine that inhibits rYkt6 in the closed state. A point mutation T46L/Q57A was shown to be able to convert yYkt6 to a more closed and dodecylphosphocholine-bound state, where Leu46 contributes key hydrophobic interactions for the closed state. We also demonstrated that the phospho-mutation S174D could shift the conformation of rYkt6 to a more open state, but the corresponding mutation S176D in yYkt6 leads to a slightly more closed conformation. These observations shed light on the regulatory mechanism underlying the variations of Ykt6 functions across species.

The T1150A cancer mutant of the protein lysine dimethyltransferase NSD2 can introduce H3K36 trimethylation.

Protein lysine methyltransferases (PKMTs) play essential roles in gene expression regulation and cancer development. Somatic mutations in PKMTs are frequently observed in cancer cells. In biochemical experiments, we show here that the NSD1 mutations Y1971C, R2017Q, and R2017L observed mostly in solid cancers are catalytically inactive suggesting that NSD1 acts as a tumor suppressor gene in these tumors. In contrast, the frequently observed T1150A in NSD2 and its T2029A counterpart in NSD1, both observed in leukemia, are hyperactive and introduce up to three methyl groups in H3K36 in biochemical and cellular assays, while wildtype NSD2 and NSD1 only introduce up to two methyl groups. In Molecular Dynamics simulations, we determined key mechanistic and structural features controlling the product specificity of this class of enzymes. Simulations with NSD2 revealed that H3K36me3 formation is possible due to an enlarged active site pocket of T1150A and loss of direct contacts of T1150 to critical residues which regulate the product specificity of NSD2. Bioinformatic analyses of published data suggested that the generation of H3K36me3 by NSD2 T1150A could alter gene regulation by antagonizing H3K27me3 finally leading to the upregulation of oncogenes.

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory.

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

A novel mRNA multi-epitope vaccine of Acinetobacter baumannii based on multi-target protein design in immunoinformatic approach.

Acinetobacter baumannii is a gram-negative bacillus prevalent in nature, capable of thriving under various environmental conditions. As an opportunistic pathogen, it frequently causes nosocomial infections such as urinary tract infections, bacteremia, and pneumonia, contributing to increased morbidity and mortality in clinical settings. Consequently, developing novel vaccines against Acinetobacter baumannii is of utmost importance. In our study, we identified 10 highly conserved antigenic proteins from the NCBI and UniProt databases for epitope mapping. We subsequently screened and selected 8 CTL, HTL, and LBL epitopes, integrating them into three distinct vaccines constructed with adjuvants. Following comprehensive evaluations of immunological and physicochemical parameters, we conducted molecular docking and molecular dynamics simulations to assess the efficacy and stability of these vaccines. Our findings indicate that all three multi-epitope mRNA vaccines designed against Acinetobacter baumannii are promising; however, further animal studies are required to confirm their reliability and effectiveness.

In silico and structure-based evaluation of deleterious mutations identified in human Chk1, Chk2, and Wee1 protein kinase.

Checkpoint kinases Chk1, Chk2, Wee1 are playing a key role in DNA damage response and genomic integrity. Cancer-associated mutations identified in human Chk1, Chk2, and Wee1 were retrieved to understand the function associated with the mutation and also alterations in the folding pattern. Therefore, an attempt has been made to identify deleterious effect of variants using in silico and structure-based approach. Variants of uncertain significance for Chk1, Chk2, and Wee1 were retrieved from different databases and four prediction servers were employed to predict pathogenicity of mutations. Further, Interpro, I-Mutant 3.0, Consurf, TM-align, and have (y)our protein explained were used for comprehensive study of the deleterious effects of variants. The sequences of Chk1, Chk2, and Wee1 were analyzed using Clustal Omega, and the three-dimensional structures of the proteins were aligned using TM-align. The molecular dynamics simulations were performed to explore the differences in folding pattern between Chk1, Chk2, Wee1 wild-type, and mutant protein and also to evaluate the structural integrity. Thirty-six variants in Chk1, 250 Variants in Chk2, and 29 in Wee1 were categorized as pathogenic using in silico prediction tools. Furthermore, 25 mutations in Chk1, 189 in Chk2, and 14 in Wee1 were highly conserved, possessing deleterious effect and also influencing the protein structure and function. These identified mutations may provide underlying genetic intricacies to serve as potential targets for therapeutic inventions and clinical management.

Revisiting Luteolin Against the Mediators of Human Metastatic Colorectal Carcinoma: A Biomolecular Approach.

Metastatic colorectal carcinoma (mCRC) is one of the prevalent subtypes of human cancers and is caused by the alterations of various lifestyle and diet-associated factors. ß-catenin, GSK-3ß, PI3K-α, AKT1, and NF-κB p50 are known to be the critical regulators of tumorigenesis and immunopathogenesis of mCRC. Unfortunately, current drugs have limited efficacy, side effects and can lead to chemoresistance. Therefore, searching for a nontoxic, efficacious anti-mCRC agent is crucial and of utmost interest. The present study demonstrates the identification of a productive and nontoxic anti-mCRC agent through a five-targets (ß-catenin, GSK-3ß, PI3K-α, AKT1, and p50)-based and three-tier (binding affinity, pharmacokinetics, and pharmacophore) screening strategy involving a series of 30 phytocompounds having a background of anti-inflammatory/anti-mCRC efficacy alongside 5-fluorouracil (FU), a reference drug. Luteolin (a phyto-flavonoid) was eventually rendered as the most potent and safe phytocompound. This inference was verified through three rounds of validation. Firstly, luteolin was found to be effective against the different mCRC cell lines (HCT-15, HCT-116, DLD-1, and HT-29) without hampering the viability of non-tumorigenic ones (RWPE-1). Secondly, luteolin was found to curtail the clonogenicity of CRC cells, and finally, it also disrupted the formation of colospheroids, a characteristic of metastasis. While studying the mechanistic insights, luteolin was found to inhibit ß-catenin activity (a key regulator of mCRC) through direct physical interactions, promoting its degradation by activating GSK3-ß and ceasing its activation by inactivating AKT1 and PI3K-α. Luteolin also inhibited p50 activity, which could be useful in mitigating mCRC-associated proinflammatory milieu. In conclusion, our study provides evidence on the efficacy of luteolin against the critical key regulators of immunopathogenesis of mCRC and recommends further studies in animal models to determine the effectiveness efficacy of this natural compound for treating mCRC in the future.

In Silico Prediction of Potential Inhibitors for Targeting RNA CAG Repeats via Molecular Docking and Dynamics Simulation: A Drug Discovery Approach.

Spinocerebellar ataxia (SCA) is a rare neurological illness inherited dominantly that causes severe impairment and premature mortality. While each rare disease may affect individuals infrequently, collectively they pose a significant healthcare challenge. It is mainly carried out due to the expansion of RNA triplet (CAG) repeats, although missense or point mutations can also be induced. Unfortunately, there is no cure; only symptomatic treatments are available. To date, SCA has about 48 subtypes, the most common of these being SCA 1, 2, 3, 6, 7, 12, and 17 having CAG repeats. Using molecular docking and molecular dynamics (MD) simulation, this study seeks to investigate effective natural herbal neuroprotective compounds against CAG repeats, which are therapeutically significant in treating SCA. Initially, virtual screening followed by molecular docking was used to estimate the binding affinity of neuroprotective natural compounds toward CAG repeats. The compound with the highest binding affinity, somniferine, was then chosen for MD simulation. The structural stability, interaction mechanism, and conformational dynamics of CAG repeats and somniferine were investigated via MD simulation. The MD study revealed that during the simulation period, the interaction between CAG repeats and somniferine stabilizes and results in fewer conformational variations. This in silico study suggests that Somniferine can be used as a therapeutic medication against RNA CAG repeats in SCA.

Concurrent inhibition of IR, ITGB1, and CD36 perturbated the interconnected network of energy metabolism and epithelial-to-mesenchymal transition in breast cancer cells.

Altered energy metabolism is an emerging hallmark of cancer and plays a pivotal in cell survival, proliferation, and biosynthesis. In a rapidly proliferating cancer, energy metabolism acts in synergism with epithelial-to-mesenchymal transition (EMT), enabling cancer stemness, dissemination, and metastasis. In this study, an interconnected functional network governing energy metabolism and EMT signaling pathways was targeted through the concurrent inhibition of IR, ITGB1, and CD36 activity. A novel multicomponent MD simulation approach was employed to portray the simultaneous inhibition of IR, ITGB1, and CD36 by a 2:1 combination of Pimozide and Ponatinib. Further, in-vitro studies revealed the synergistic anticancer efficacy of drugs against monolayer as well as tumor spheroids of breast cancer cell lines (MCF-7 and MDA-MB-231). In addition, the combination therapy exerted approximately 40% of the apoptotic population and more than 1.5- to 3-fold reduction in the expression of ITGB1, IR, p-IR, IRS-1, and p-AKT in MCF-7 and MDA-MB-231 cell lines. Moreover, the reduction in fatty acid uptake, lipid droplet accumulation, cancer stemness, and migration properties were also observed. Thus, targeting IR, ITGB1, and CD36 in the interconnected network with the combination of Pimozide and Ponatinib represents a promising therapeutic approach for breast cancer.

Validating Fractalkine receptor as a target and identifying candidates for drug discovery against type 2 diabetes.

Type 2 diabetes mellitus (T2DM) is one of the most common chronic diseases employing abnormal levels of insulin. Enhancing the insulin production is greatly aided by the regulatory mechanisms of the Fractalkine receptor (CX3CR1) system in islet ß-cell function. However, elements including a high-fat diet, obesity, and ageing negatively impact the expression of CX3CR1 in islets. CX3CL1/CX3CR1 receptor-ligand complex is now recognized as a novel therapeutic target. It suggests that T2DM-related ß-cell dysfunction may result from lower amount of these proteins. We analyzed the differential expression of CX3CR1 gene samples taken from persons with T2DM using data obtained from the Gene Expression Omnibus database. Homology modeling enabled us to generate the three-dimensional structure of CX3CR1 and a possible binding pocket. The optimized CX3CR1 structure was subjected to rigorous screening against a massive library of 693 million drug-like molecules from the ZINC15 database. This screening process led to the identification of three compounds with strong binding affinity at the identified binding pocket of CX3CR1. To further evaluate the potential of these compounds, molecular dynamics simulations were conducted over a 50 ns time scale to assess the stability of the protein-ligand complexes. These simulations revealed that ZINC000032506419 emerged as the most promising drug-like compound among the three potent molecules. The discovery of ZINC000032506419 holds exciting promise as a potential therapeutic agent for T2D and other related metabolic disorders. These findings pave the way for the development of effective medications to address the complexities of T2DM and its associated metabolic diseases.

In Silico and In Vitro Study of Isoquercitrin against Kidney Cancer and Inflammation by Triggering Potential Gene Targets.

Kidney cancer has emerged as a major medical problem in recent times. Multiple compounds are used to treat kidney cancer by triggering cancer-causing gene targets. For instance, isoquercitrin (quercetin-3-O-ß-d-glucopyranoside) is frequently present in fruits, vegetables, medicinal herbs, and foods and drinks made from plants. Our previous study predicted using protein-protein interaction (PPI) and molecular docking analysis that the isoquercitrin compound can control kidney cancer and inflammation by triggering potential gene targets of IGF1R, PIK3CA, IL6, and PTGS2. So, the present study is about further in silico and in vitro validation. We performed molecular dynamic (MD) simulation, gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, cytotoxicity assay, and RT-PCR and qRT-PCR validation. According to the MD simulation (250 ns), we found that IGF1R, PIK3CA, and PTGS2, except for IL6 gene targets, show stable binding energy with a stable complex with isoquercitrin. We also performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the final targets to determine their regulatory functions and signaling pathways. Furthermore, we checked the cytotoxicity effect of isoquercitrin (IQ) and found that 5 µg/mL and 10 µg/mL doses showed higher cell viability in a normal kidney cell line (HEK 293) and also inversely showed an inhibition of cell growth at 35% and 45%, respectively, in the kidney cancer cell line (A498). Lastly, the RT-PCR and qRT-PCR findings showed a significant decrease in PTGS2, PIK3CA, and IGF1R gene expression, except for IL6 expression, following dose-dependent treatments with IQ. Thus, we can conclude that isoquercitrin inhibits the expression of PTGS2, PIK3CA, and IGF1R gene targets, which in turn controls kidney cancer and inflammation.

Structural insights in to the atypical type-I ABC Glucose-6-phosphate importer VCA0625-27 of Vibrio cholerae.

Sugar phosphates are potential sources of carbon and phosphate for bacteria. Despite that the process of internalization of Glucose-6-Phosphate (G6P) through plasma membrane remained elusive in several bacteria. VCA0625-27, made of periplasmic ligand binding protein (PLBP) VCA0625, an atypical monomeric permease VCA0626, and a cytosolic ATPase VCA0627, recently emerged as hexose-6-phosphate uptake system of Vibrio cholerae. Here we report high resolution crystal structure of VCA0625 in G6P bound state that largely resembles AfuA of Actinobacillus pleuropneumoniae. MD simulations on VCA0625 in apo and G6P bound states unraveled an 'open to close' and swinging bi-lobal motions, which are diminished upon G6P binding. Mutagenesis followed by biochemical assays on VCA0625 underscored that R34 works as gateway to bind G6P. Although VCA0627 binds ATP, it is ATPase deficient in the absence of VCA0625 and VCA0626, which is a signature phenomenon of type-I ABC importer. Further, modeling, docking and systematic sequence analysis allowed us to envisage the existence of similar atypical type-I G6P importer with fused monomeric permease in 27 other gram-negative bacteria.

Structural basis of pri-let-7 recognition by human pseudouridine synthase TruB1.

Let-7 was one of the first microRNAs (miRNAs) to be discovered and its expression promotes differentiation during development and function as tumor suppressors in various cancers. The maturation process of let-7 miRNA is tightly regulated by multiple RNA-binding proteins. For example, LIN28 binds to the terminal loops of the precursors of let-7 family and block their processing into mature miRNAs. Trim25 promotes the uridylation-mediated degradation of pre-let-7 modified by LIN28/TUT4. Recently, human pseudouridine synthase TruB1 has been reported to facilitate let-7 maturation by directly binding to pri-let-7 and recruiting Drosha-DGCR8 microprocessor. Through biochemical assay and structural investigation, we show that human TruB1 binds specifically the terminal loop of pri-let-7a1 at nucleotides 31-41, which folds as a small stem-loop architecture. Although TruB1 recognizes the terminal loop of pri-let-7a1 in a way similar to how E. coli TruB interacts with tRNA, a conserved KRKK motif in human and other higher eukaryotes adds an extra binding interface and strengthens the recognition of TruB1 for pri-let-7a1 through electrostatic interactions. These findings reveal the structural basis of TruB1-pri-let-7 interaction which may assists the elucidation of precise role of TruB1 in biogenesis of let-7.