Mesoscale DNA feature in antibody-coding sequence facilitates somatic hypermutation.

Somatic hypermutation (SHM), initiated by activation-induced cytidine deaminase (AID), generates mutations in the antibody-coding sequence to allow affinity maturation. Why these mutations intrinsically focus on the three nonconsecutive complementarity-determining regions (CDRs) remains enigmatic. Here, we found that predisposition mutagenesis depends on the single-strand (ss) DNA substrate flexibility determined by the mesoscale sequence surrounding AID deaminase motifs. Mesoscale DNA sequences containing flexible pyrimidine-pyrimidine bases bind effectively to the positively charged surface patches of AID, resulting in preferential deamination activities. The CDR hypermutability is mimicable in in vitro deaminase assays and is evolutionarily conserved among species using SHM as a major diversification strategy. We demonstrated that mesoscale sequence alterations tune the in vivo mutability and promote mutations in an otherwise cold region in mice. Our results show a non-coding role of antibody-coding sequence in directing hypermutation, paving the way for the synthetic design of humanized animal models for optimal antibody discovery and explaining the AID mutagenesis pattern in lymphoma.

A kinetic model reveals the critical gating motifs for donor-substrate loading into Actinobacillus pleuropneumoniae N-glycosyltransferase.

Soluble N-glycosyltransferase from Actinobacillus pleuropneumoniae (ApNGT) catalyzes the glycosylation of asparagine residues, and represents one of the most encouraging biocatalysts for N-glycoprotein production. Since the sugar tolerance of ApNGT is restricted to limited monosaccharides (e.g., Glc, GlcN, Gal, Xyl, and Man), tremendous efforts are devoted to expanding the substrate scope of ApNGT via enzyme engineering. However, rational design of novel NGT variants suffers from an elusive understanding of the substrate-binding process from a dynamic point of view. Here, by employing extensive all-atom molecular dynamics (MD) simulations integrated with a kinetic model, we reveal, at the atomic level, the complete donor-substrate binding process from the bulk solvent to the ApNGT active-site, and the key intermediate states of UDP-Glc during its loading dynamics. We are able to determine the critical transition event that limits the overall binding rate, which guides us to pinpoint the key ApNGT residues dictating the donor-substrate entry. The functional roles of several identified gating residues were evaluated through site-directed mutagenesis and enzymatic assays. Two single-point mutations, N471A and S496A, could profoundly enhance the catalytic activity of ApNGT. Our work provides deep mechanistic insights into the structural dynamics of the donor-substrate loading process for ApNGT, which sets a rational basis for design of novel NGT variants with desired substrate specificity.

O-glycosylation of SARS-CoV-2 spike protein by host O-glycosyltransferase strengthens its trimeric structure.

Protein O-glycosylation, also known as mucin-type O-glycosylation, is one of the most abundant glycosylation in mammalian cells. It is initially catalyzed by a family of polypeptide GalNAc transferases (ppGalNAc-Ts). The trimeric spike protein (S) of SARS-CoV-2 is highly glycosylated and facilitates the virus's entry into host cells and membrane fusion of the virus. However, the functions and relationship between host ppGalNAc-Ts and O-glycosylation on the S protein remain unclear. Herein, we identify 15 O-glycosites and 10 distinct O-glycan structures on the S protein using an HCD-product-dependent triggered ETD mass spectrometric analysis. We observe that the isoenzyme T6 of ppGalNAc-Ts (ppGalNAc-T6) exhibits high O-glycosylation activity for the S protein, as demonstrated by an on-chip catalytic assay. Overexpression of ppGalNAc-T6 in HEK293 cells significantly enhances the O-glycosylation level of the S protein, not only by adding new O-glycosites but also by increasing O-glycan heterogeneity. Molecular dynamics simulations reveal that O-glycosylation on the protomer-interface regions, modified by ppGalNAc-T6, potentially stabilizes the trimeric S protein structure by establishing hydrogen bonds and non-polar interactions between adjacent protomers. Furthermore, mutation frequency analysis indicates that most O-glycosites of the S protein are conserved during the evolution of SARS-CoV-2 variants. Taken together, our finding demonstrate that host O-glycosyltransferases dynamically regulate the O-glycosylation of the S protein, which may influence the trimeric structural stability of the protein. This work provides structural insights into the functional role of specific host O-glycosyltransferases in regulating the O-glycosylation of viral envelope proteins.

Loading dynamics of one SARS-CoV-2-derived peptide into MHC-II revealed by kinetic models.

Major histocompatibility complex class II (MHC-II) plays an indispensable role in activating CD4+ T cell immune responses by presenting antigenic peptides on the cell surface for recognition by T cell receptors. The assembly of MHC-II and antigenic peptide is therefore a prerequisite for the antigen presentation. To date, however, the atomic-level mechanism underlying the peptide-loading dynamics for MHC-II is still elusive. Here, by constructing Markov state models based on extensive all-atom molecular dynamics simulations, we reveal the complete peptide-loading dynamics into MHC-II for one SARS-CoV-2 S-protein-derived antigenic peptide (235ITRFQTLLALHRSYL249). Our Markov state model identifies six metastable states (S1-S6) during the peptide-loading process and determines two dominant loading pathways. The peptide could potentially approach the antigen-binding groove via either its N- or C-terminus. Then, the consecutive insertion of several anchor residues into the binding pockets profoundly dictates the peptide-loading dynamics. Notably, the MHC-II αA52-E55 motif could guide the peptide loading into the antigen-binding groove via forming ß-sheets conformation with the incoming peptide. The rate-limiting step, namely S5→S6, is mainly attributed to a considerable desolvation penalty triggered by the binding of the peptide C-terminus. Moreover, we further examined the conformational changes associated with the peptide exchange process catalyzed by the chaperon protein HLA-DM. A flipped-out conformation of MHC-II αW43 captured in S1-S3 is considered a critical anchor point for HLA-DM to modulate the structural dynamics. Our work provides deep structural insights into the key regulatory factors in MHC-II responsible for peptide recognition and guides future design for peptide vaccines against SARS-CoV-2.

Phosphorylation Regulation Mechanism of ß2 Integrin for the Binding of Filamin Revealed by Markov State Model.

Leukocyte adhesion deficiency-1 (LAD-1) disorder is a severe immunodeficiency syndrome caused by deficiency or mutation of ß2 integrin. The phosphorylation on threonine 758 of ß2 integrin acts as a molecular switch inhibiting the binding of filamin. However, the switch mechanism of site-specific phosphorylation at the atom level is still poorly understood. To resolve the regulation mechanism, all-atom molecular dynamics simulation and Markov state model were used to study the dynamic regulation pathway of phosphorylation. Wild type system possessed lower binding free energy and fewer number of states than the phosphorylated system. Both systems underwent local disorder-to-order conformation conversion when achieving steady states. To reach steady states, wild type adopted less number of transition paths/shortest path according to the transition path theory than the phosphorylated system. The underlying phosphorylated regulation pathway was from P1 to P0 and then P4 state, and the main driving force should be hydrogen bond and hydrophobic interaction disturbing the secondary structure of phosphorylated states. These studies will shed light on the pathogenesis of LAD-1 disease and lay a foundation for drug development.

Atomic resolution of short-range sliding dynamics of thymine DNA glycosylase along DNA minor-groove for lesion recognition.

Thymine DNA glycosylase (TDG), as a repair enzyme, plays essential roles in maintaining the genome integrity by correcting several mismatched/damaged nucleobases. TDG acquires an efficient strategy to search for the lesions among a vast number of cognate base pairs. Currently, atomic-level details of how TDG translocates along DNA as it approaches the lesion site and the molecular mechanisms of the interplay between TDG and DNA are still elusive. Here, by constructing the Markov state model based on hundreds of molecular dynamics simulations with an integrated simulation time of ∼25 µs, we reveal the rotation-coupled sliding dynamics of TDG along a 9 bp DNA segment containing one G·T mispair. We find that TDG translocates along DNA at a relatively faster rate when distant from the lesion site, but slows down as it approaches the target, accompanied by deeply penetrating into the minor-groove, opening up the mismatched base pair and significantly sculpturing the DNA shape. Moreover, the electrostatic interactions between TDG and DNA are found to be critical for mediating the TDG translocation. Notably, several uncharacterized TDG residues are identified to take part in regulating the conformational switches of TDG occurred in the site-transfer process, which warrants further experimental validations.

Mechanism of REST/NRSF regulation of clustered protocadherin α genes.

Repressor element-1 silencing transcription factor (REST) or neuron-restrictive silencer factor (NRSF) is a zinc-finger (ZF) containing transcriptional repressor that recognizes thousands of neuron-restrictive silencer elements (NRSEs) in mammalian genomes. How REST/NRSF regulates gene expression remains incompletely understood. Here, we investigate the binding pattern and regulation mechanism of REST/NRSF in the clustered protocadherin (PCDH) genes. We find that REST/NRSF directionally forms base-specific interactions with NRSEs via tandem ZFs in an anti-parallel manner but with striking conformational changes. In addition, REST/NRSF recruitment to the HS5-1 enhancer leads to the decrease of long-range enhancer-promoter interactions and downregulation of the clustered PCDHα genes. Thus, REST/NRSF represses PCDHα gene expression through directional binding to a repertoire of NRSEs within the distal enhancer and variable target genes.

DNA Deformation Exerted by Regulatory DNA-Binding Motifs in Human Alkyladenine DNA Glycosylase Promotes Base Flipping.

Human alkyladenine DNA glycosylase (AAG) is a key enzyme that corrects a broad range of alkylated and deaminated nucleobases to maintain genomic integrity. When encountering the lesions, AAG adopts a base-flipping strategy to extrude the target base from the DNA duplex to its active site, thereby cleaving the glycosidic bond. Despite its functional importance, the detailed mechanism of such base extrusion and how AAG distinguishes the lesions from an excess of normal bases both remain elusive. Here, through the Markov state model constructed on extensive all-atom molecular dynamics simulations, we find that the alkylated nucleobase (N3-methyladenine, 3MeA) everts through the DNA major groove. Two key AAG motifs, the intercalation and E131-N146 motifs, play active roles in bending/pressing the DNA backbone and widening the DNA minor groove during 3MeA eversion. In particular, the intercalated residue Y162 is involved in buckling the target site at the early stage of 3MeA eversion. Our traveling-salesman based automated path searching algorithm further revealed that a non-target normal adenine tends to be trapped in an exo site near the active site, which however barely exists for a target base 3MeA. Collectively, these results suggest that the Markov state model combined with traveling-salesman based automated path searching acts as a promising approach for studying complex conformational changes of biomolecules and dissecting the elaborate mechanism of target recognition by this unique enzyme.

Dynamics of peptide loading into major histocompatibility complex class I molecules chaperoned by TAPBPR.

Major histocompatibility complex class I (MHC-I) molecules display antigenic peptides on the cell surface for T cell receptor scanning, thereby activating the immune response. Peptide loading into MHC-I molecules is thus a critical step during the antigen presentation process. Chaperone TAP-binding protein related (TAPBPR) plays a critical role in promoting high-affinity peptide loading into MHC-I, by discriminating against the low-affinity ones. However, the complete peptide loading dynamics into TAPBPR-bound MHC-I is still elusive. Here, we constructed kinetic network models based on hundreds of short-time MD simulations with an aggregated simulation time of ∼21.7 µs, and revealed, at atomic level, four key intermediate states of one antigenic peptide derived from melanoma-associated MART-1/Melan-A protein during its loading process into TAPBPR-bound MHC-I. We find that the TAPBPR binding at the MHC-I pocket-F can substantially reshape the distant pocket-B via allosteric regulations, which in turn promotes the following peptide N-terminal loading. Intriguingly, the partially loaded peptide could profoundly weaken the TAPBPR-MHC stability, promoting the dissociation of the TAPBPR scoop-loop (SL) region from the pocket-F to a more solvent-exposed conformation. Structural inspections further indicate that the peptide loading could remotely affect the SL binding site through both allosteric perturbations and direct contacts. In addition, another structural motif of TAPBPR, the jack hairpin region, was also found to participate in mediating the peptide editing. Our study sheds light on the detailed molecular mechanisms underlying the peptide loading process into TAPBPR-bound MHC-I and pinpoints the key structural factors responsible for dictating the peptide-loading dynamics.

Computational investigations on target-site searching and recognition mechanisms by thymine DNA glycosylase during DNA repair process.

DNA glycosylase, as one member of DNA repair machineries, plays an essential role in correcting mismatched/damaged DNA nucleotides by cleaving the N-glycosidic bond between the sugar and target nucleobase through the base excision repair (BER) pathways. Efficient corrections of these DNA lesions are critical for maintaining genome integrity and preventing premature aging and cancers. The target-site searching/recognition mechanisms and the subsequent conformational dynamics of DNA glycosylase, however, remain challenging to be characterized using experimental techniques. In this review, we summarize our recent studies of sequential structural changes of thymine DNA glycosylase (TDG) during the DNA repair process, achieved mostly by molecular dynamics (MD) simulations. Computational simulations allow us to reveal atomic-level structural dynamics of TDG as it approaches the target-site, and pinpoint the key structural elements responsible for regulating the translocation of TDG along DNA. Subsequently, upon locating the lesions, TDG adopts a base-flipping mechanism to extrude the mispaired nucleobase into the enzyme active-site. The constructed kinetic network model elucidates six metastable states during the base-extrusion process and suggests an active role of TDG in flipping the intrahelical nucleobase. Finally, the molecular mechanism of product release dynamics after catalysis is also summarized. Taken together, we highlight to what extent the computational simulations advance our knowledge and understanding of the molecular mechanism underlying the conformational dynamics of TDG, as well as the limitations of current theoretical work.

Allosteric regulation in CRISPR/Cas1-Cas2 protospacer acquisition mediated by DNA and Cas2.

Cas1 and Cas2 are highly conserved proteins across clustered-regularly-interspaced-short-palindromic-repeat-Cas systems and play a significant role in protospacer acquisition. Based on crystal structure of twofold symmetric Cas1-Cas2 in complex with dual-forked protospacer DNA (psDNA), we conducted all-atom molecular dynamics simulations to study the psDNA binding, recognition, and response to cleavage on the protospacer-adjacent-motif complementary sequence, or PAMc, of Cas1-Cas2. In the simulation, we noticed that two active sites of Cas1 and Cas1' bind asymmetrically to two identical PAMc on the psDNA captured from the crystal structure. For the modified psDNA containing only one PAMc, as that to be recognized by Cas1-Cas2 in general, our simulations show that the non-PAMc association site of Cas1-Cas2 remains destabilized until after the stably bound PAMc being cleaved at the corresponding association site. Thus, long-range correlation appears to exist upon the PAMc cleavage between the two active sites (∼10 nm apart) on Cas1-Cas2, which can be allosterically mediated by psDNA and Cas2 and Cas2' in bridging. To substantiate such findings, we conducted repeated runs and further simulated Cas1-Cas2 in complex with synthesized psDNA sequences psL and psH, which have been measured with low and high frequency in acquisition, respectively. Notably, such intersite correlation becomes even more pronounced for the Cas1-Cas2 in complex with psH but remains low for the Cas1-Cas2 in complex with psL. Hence, our studies demonstrate that PAMc recognition and cleavage at one active site of Cas1-Cas2 may allosterically regulate non-PAMc association or even cleavage at the other site, and such regulation can be mediated by noncatalytic Cas2 and DNA protospacer to possibly support the ensued psDNA acquisition.

DeepAntigen: a novel method for neoantigen prioritization via 3D genome and deep sparse learning.

MOTIVATION: The mutations of cancers can encode the seeds of their own destruction, in the form of T-cell recognizable immunogenic peptides, also known as neoantigens. It is computationally challenging, however, to accurately prioritize the potential neoantigen candidates according to their ability of activating the T-cell immunoresponse, especially when the somatic mutations are abundant. Although a few neoantigen prioritization methods have been proposed to address this issue, advanced machine learning model that is specifically designed to tackle this problem is still lacking. Moreover, none of the existing methods considers the original DNA loci of the neoantigens in the perspective of 3D genome which may provide key information for inferring neoantigens' immunogenicity. RESULTS: In this study, we discovered that DNA loci of the immunopositive and immunonegative MHC-I neoantigens have distinct spatial distribution patterns across the genome. We therefore used the 3D genome information along with an ensemble pMHC-I coding strategy, and developed a group feature selection-based deep sparse neural network model (DNN-GFS) that is optimized for neoantigen prioritization. DNN-GFS demonstrated increased neoantigen prioritization power comparing to existing sequence-based approaches. We also developed a webserver named deepAntigen (http://yishi.sjtu.edu.cn/deepAntigen) that implements the DNN-GFS as well as other machine learning methods. We believe that this work provides a new perspective toward more accurate neoantigen prediction which eventually contribute to personalized cancer immunotherapy. AVAILABILITY AND IMPLEMENTATION: Data and implementation are available on webserver: http://yishi.sjtu.edu.cn/deepAntigen. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Ikarugamycin inhibits pancreatic cancer cell glycolysis by targeting hexokinase 2.

Mangrove-derived actinobacteria strains are well-known for producing novel secondary metabolites. The polycyclic tetramate macrolactam (PTM), ikarugamycin (IKA) isolated from Streptomyces xiamenensis 318, exhibits antiproliferative activities against pancreatic ductal adenocarcinoma (PDAC) in vitro. However, the protein target for bioactive IKA is unclear. In this study, whole transcriptome-based profiling revealed that the glycolysis pathway is significantly affected by IKA. Metabolomic studies demonstrated that IKA treatment induces a significant drop in glucose-6-phosphate and a slight increase in intracellular glucose level. Analysis of glucose consumption, lactate production, and the extracellular acidification rate confirmed the inhibitory role of IKA on the glycolytic flux in PDAC cells. Surface plasmon resonance (SPR) experiments and docking studies identified the key enzyme of glycolysis, hexokinase 2 (HK2), as a molecular target of IKA. Moreover, IKA reduced tumor size without overt cytotoxicity in mice with PDAC xenografts and increased chemotherapy response to gemcitabine in PDAC cells in vitro. Taken together, IKA can block glycolysis in pancreatic cancer by targeting HK2, which may be a potential drug candidate for PDAC treatment.

Determining selection free energetics from nucleotide pre-insertion to insertion in viral T7 RNA polymerase transcription fidelity control.

An elongation cycle of a transcribing RNA polymerase (RNAP) usually consists of multiple kinetics steps, so there exist multiple kinetic checkpoints where non-cognate nucleotides can be selected against. We conducted comprehensive free energy calculations on various nucleotide insertions for viral T7 RNAP employing all-atom molecular dynamics simulations. By comparing insertion free energy profiles between the non-cognate nucleotide species (rGTP and dATP) and a cognate one (rATP), we obtained selection free energetics from the nucleotide pre-insertion to the insertion checkpoints, and further inferred the selection energetics down to the catalytic stage. We find that the insertion of base mismatch rGTP proceeds mainly through an off-path along which both pre-insertion screening and insertion inhibition play significant roles. In comparison, the selection against dATP is found to go through an off-path pre-insertion screening along with an on-path insertion inhibition. Interestingly, we notice that two magnesium ions switch roles of leave and stay during the dATP on-path insertion. Finally, we infer that substantial selection energetic is still required to catalytically inhibit the mismatched rGTP to achieve an elongation error rate ∼10-4 or lower; while no catalytic selection seems to be further needed against dATP to obtain an error rate ∼10-2.

Key structural motifs in Thymine DNA glycosylase responsible for recognizing certain DNA bent conformation revealed by atomic simulations.

Knowledge of how DNA bending facilitates the target-base searching by Thymine DNA glycosylase (TDG) is of major importance for unraveling the recognition mechanism between DNA and TDG in DNA repair process. An atomic-level understanding of the initial encounter between TDG and DNA before base-flipping, however, is still elusive. Here, we employ all-atom molecular dynamics (MD) simulations with an integrated simulation time of ∼3 µs to investigate how TDG responses to different DNA bending conformations. By constructing several TDG-DNA complexes with varied DNA bend angles (ranging from ∼0° to 60°), we pinpoint the key TDG motifs responsible for recognizing certain DNA bending conformations. Particularly, several positively charged residues, i.e., Lys232, Lys240, and Lys246, are critical for the tight binding with DNA backbones. Importantly, the roll-angle patterns, rather than the tilt and twist angles, are found to be strongly correlated with the extent of DNA bending, which in turn, governs the TDG recognition. Further comparisons between the naked and TDG-bound DNA conformations reveal that the TDG binding can impose a substantial DNA deformation, resulting in profound roll-angle alterations. Our studies warrant further experimental validations and provide deep structural insights into the recognition mechanism between TDG and DNA during their initial encounter.

Refolding Dynamics of gp41 from Pre-fusion to Pre-hairpin States during HIV-1 Entry.

The HIV-1 infection is triggered by the binding of the viral envelope glycoprotein (Env) gp120-gp41 trimer to host-cell receptor CD4 and co-receptor CCR5/CXCR4, which leads to substantial conformational changes of Env, that is, structural transition of gp120 from a closed to an open state followed by gp41 refolding from pre-fusion to post-fusion states. The latter finally promotes membrane fusion, likely via visiting a critical pre-hairpin state of gp41. The complete conformational dynamics of the pre-hairpin formation at atomic resolution, however, is still unknown. Here, by constructing a Markov state model based on the all-atom molecular dynamics (MD) with an aggregated simulation time of ∼24 µs, we reveal the gp41 refolding dynamics from pre-fusion to pre-hairpin state and the key metastable states involved. Moreover, we further explored the drug resistance mechanism of two C-terminal heptad repeat-derived gp41 inhibitors, T20 and sifuvirtide, based on the constructed inhibitor-bound gp41 pre-hairpin complexes. The results indicate that these two inhibitors have distinct binding sites on gp41 but share a common drug resistance region that usually exhibits a helical structure in the pre-hairpin state yet adopts various secondary structures in other metastable states. Moreover, we conducted several mutant MD simulations to further investigate the mechanisms of how some drug-resistant mutations might affect the pre-hairpin formation, which in turn prevent the inhibitor recognition. Our findings provide deep structural insights into the molecular mechanisms of the pre-hairpin formation for gp41, which helps to guide future anti-HIV drug design.

Base-flipping dynamics from an intrahelical to an extrahelical state exerted by thymine DNA glycosylase during DNA repair process.

Thymine DNA glycosylase (TDG) is a DNA repair enzyme that excises a variety of mismatched or damaged nucleotides (nts), e.g. dU, dT, 5fC and 5caC. TDG is shown to play essential roles in maintaining genome integrity and correctly programming epigenetic modifications through DNA demethylation. After locating the lesions, TDG employs a base-flipping strategy to recognize the damaged nucleobases, whereby the interrogated nt is extruded from the DNA helical stack and binds into the TDG active site. The dynamic mechanism of the base-flipping process at an atomistic resolution, however, remains elusive. Here, we employ the Markov State Model (MSM) constructed from extensive all-atom molecular dynamics (MD) simulations to reveal the complete base-flipping process for a G.T mispair at a tens of microsecond timescale. Our studies identify critical intermediates of the mispaired dT during its extrusion process and reveal the key TDG residues involved in the inter-state transitions. Notably, we find an active role of TDG in promoting the intrahelical nt eversion, sculpturing the DNA backbone, and penetrating into the DNA minor groove. Three additional TDG substrates, namely dU, 5fC, and 5caC, are further tested to evaluate the substituent effects of various chemical modifications of the pyrimidine ring on base-flipping dynamics.

Dynamics of the excised base release in thymine DNA glycosylase during DNA repair process.

Thymine DNA glycosylase (TDG) initiates base excision repair by cleaving the N-glycosidic bond between the sugar and target base. After catalysis, the release of excised base is a requisite step to terminate the catalytic cycle and liberate the TDG for the following enzymatic reactions. However, an atomistic-level understanding of the dynamics of the product release process in TDG remains unknown. Here, by employing molecular dynamics simulations combined with the Markov State Model, we reveal the dynamics of the thymine release after the excision at microseconds timescale and all-atom resolution. We identify several key metastable states of the thymine and its dominant releasing pathway. Notably, after replacing the TDG residue Gly142 with tyrosine, the thymine release is delayed compared to the wild-type (wt) TDG, as supported by our potential of mean force (PMF) calculations. These findings warrant further experimental tests to potentially trap the excised base in the active site of TDG after the catalysis, which had been unsuccessful by previous attempts. Finally, we extended our studies to other TDG products, including the uracil, 5hmU, 5fC and 5caC bases in order to compare the product release for different targeting bases in the TDG-DNA complex.

Understanding the molecular mechanism of umami recognition by T1R1-T1R3 using molecular dynamics simulations.

Taste receptor T1R1-T1R3 can be activated by binding to several natural ligands, e.g., l-glutamate and 5'-ribonucleotides etc., thereby stimulating the umami taste. The molecular mechanism of umami recognition at atomic details, however, remains elusive. Here, using homology modeling, molecular docking and molecular dynamics (MD) simulations, we investigate the effects of five natural umami ligands on the structural dynamics of T1R1-T1R3. Our work identifies the key residues that are directly involved in recognizing the binding ligands. In addition, two adjacent binding sites in T1R1 are determined for substrate binding, and depending on the molecular size and chemical properties of the incoming ligand, one or both binding sites can be occupied. More interestingly, the ligand binding can modulate the pocket size, which is likely correlated with the closing and opening motions of T1R1. We then classify these five ligands into two groups according to their different binding effects on T1R1, which likely associate with the distinct umami signals stimulated by various ligands. This work warrants new experimental assays to further validate the theoretical model and provides guidance to design more effective umami ligands.

Inhibition of polypeptide N-acetyl-α-galactosaminyltransferases is an underlying mechanism of dietary polyphenols preventing colorectal tumorigenesis.

Ellagitannin-derived ellagic acid (EA) and colonic metabolite urolithins are functional dietary ingredients for cancer prevention, but the underlying mechanism need elucidation. Mucin-type O-glycosylation, initiated by polypeptide N-acetyl-α-galactosaminyltransferases (ppGalNAc-Ts), fine-tunes multiple biological processes and is closely associated with cancer progression. Herein, we aim to explore how specific tannin-based polyphenols affect tumor behavior of colorectal cancer cells (CRC) by modulating O-glycosylation. Utilizing HPLC-based enzyme assay, we find urolithin D (UroD), EA and gallic acid (GA) potently inhibit ppGalNAc-Ts. In particular, UroD inhibits ppGalNAc-T2 through a peptide/protein-competitive manner with nanomolar affinity. Computational simulations combined with site-directed mutagenesis further support the inhibitors' mode of action. Moreover, lectin analysis and metabolic labelling reveal that UroD can reduce cell O-glycans but not N-glycans. Transwell experiments prove that UroD inhibits migration and invasion of CRC cells. Our work proves that specific tannin-based polyphenols can potently inhibit ppGalNAc-Ts activity to reduce cell O-glycosylation and lead to lowering the migration and invasion of CRC cells, suggesting that disturbance of mucin-type O-glycosylation is an important mechanism for the function of dietary polyphenols.