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.

Structural and mechanistic investigations on CC bond forming α-oxoamine synthase allowing L-glutamate as substrate.

Carbon­carbon (C-C) bonds serve as the fundamental structural backbone of organic molecules. As a critical CC bond forming enzyme, α-oxoamine synthase is responsible for the synthesis of α-amino ketones by performing the condensation reaction between amino acids and acyl-CoAs. We previously identified an α-oxoamine synthase (AOS), named as Alb29, involved in albogrisin biosynthesis in Streptomyces albogriseolus MGR072. This enzyme belongs to the α-oxoamine synthase family, a subfamily under the pyridoxal 5'-phosphate (PLP) dependent enzyme superfamily. In this study, we report the crystal structures of Alb29 bound to PLP and L-Glu, which provide the atomic-level structural insights into the substrate recognition by Alb29. We discover that Alb29 can catalyze the amino transformation from L-Gln to L-Glu, besides the condensation of L-Glu with ß-methylcrotonyl coenzyme A. Subsequent structural analysis has revealed that one flexible loop in Alb29 plays an important role in both amino transformation and condensation. Based on the crystal structure of the S87G mutant in the loop region, we capture two distinct conformations of the flexible loop in the active site, compared with the wild-type Alb29. Our study offers valuable insights into the catalytic mechanism underlying substrate recognition of Alb29.

Investigation of the Catalytic Mechanism of a Soluble N-glycosyltransferase Allows Synthesis of N-glycans at Noncanonical Sequons.

The soluble N-glycosyltransferase from Actinobacillus pleuropneumoniae (ApNGT) can establish an N-glycosidic bond at the asparagine residue in the Asn-Xaa-Ser/Thr consensus sequon and is one of the most promising tools for N-glycoprotein production. Here, by integrating computational and experimental strategies, we revealed the molecular mechanism of the substrate recognition and following catalysis of ApNGT. These findings allowed us to pinpoint a key structural motif (215DVYM218) in ApNGT responsible for the peptide substrate recognition. Moreover, Y222 and H371 of ApNGT were found to participate in activating the acceptor Asn. The constructed models were supported by further crystallographic studies and the functional roles of the identified residues were validated by measuring the glycosylation activity of various mutants against a library of synthetic peptides. Intriguingly, with particular mutants, site-selective N-glycosylation of canonical or noncanonical sequons within natural polypeptides from the SARS-CoV-2 spike protein could be achieved, which were used to investigate the biological roles of the N-glycosylation in membrane fusion during virus entry. Our study thus provides in-depth molecular mechanisms underlying the substrate recognition and catalysis for ApNGT, leading to the synthesis of previously unknown chemically defined N-glycoproteins for exploring the biological importance of the N-glycosylation at a specific site.

Dynamic Phosphorylation of G9a Regulates its Repressive Activity on Chromatin Accessibility and Mitotic Progression.

Phosphorylation of Ser10 of histone H3 (H3S10p), together with the adjacent methylation of Lys9 (H3K9me), has been proposed to function as a 'phospho-methyl switch' to regulate mitotic chromatin architecture. Despite of immense understanding of the roles of H3S10 phosphorylation, how H3K9me2 are dynamically regulated during mitosis is poorly understood. Here, it is identified that Plk1 kinase phosphorylates the H3K9me1/2 methyltransferase G9a/EHMT2 at Thr1045 (pT1045) during early mitosis, which attenuates its catalytic activity toward H3K9me2. Cells bearing Thr1045 phosphomimic mutant of G9a (T1045E) show decreased H3K9me2 levels, increased chromatin accessibility, and delayed mitotic progression. By contrast, dephosphorylation of pT1045 during late mitosis by the protein phosphatase PPP2CB reactivates G9a activity and upregulates H3K9me2 levels, correlated with decreased levels of H3S10p. Therefore, the results provide a mechanistic explanation of the essential of a 'phospho-methyl switch' and highlight the importance of Plk1 and PPP2CB-mediated dynamic regulation of G9a activity in chromatin organization and mitotic progression.

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.

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.

Legumain inhibitor prevents breast cancer bone metastasis by attenuating osteoclast differentiation and function.

Breast cancer is the main lethal disease among females, and metastasis to lung and bone poses a serious threat to patients' life. Therefore, identification of novel molecular mediators that can potentially be exploited as therapeutic targets for treating osteolytic bone metastases is needed. A murine model of breast cancer bone metastasis was developed by injection of 4 T1.2 cells into the left ventricle and hence directly into the arterial system leading to bone. AEP (Asparagine endopeptidase) inhibitor combined with epirubicin or epirubicin alone was administered by intraperitoneal injection into animal model. The presence of bone metastatic and osteolytic lesions in bone were assessed by bioluminescent imaging and X-rays analysis. The expression of EMT (Epithelial-Mesenchymal Transition) relevant genes were examined by Western blotting. Cell migration and invasion were investigated with a transwell assay. Compound BIC-113, small molecule inhibitors of AEP, inhibited AEP enzymatic activity in breast cancer cell lines, and affected invasion and migration of cancer cells, but had no effect on cell growth. In animal model of breast cancer bone metastasis, compound BIC-113 combined with epirubicin inhibited breast cancer bone metastasis and attenuated breast cancer osteolytic lesions in bone by inhibiting osteoclast differentiation and EMT. These results indicate that compound BIC-113 combined with epirubicin has the potential to be used in breast cancer therapy by preventing bone metastasis via improving E-cadherin expression and inhibition of osteoclast formation.

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.

Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion.

One-dimensional (1-D) sliding of transcription factor (TF) protein along DNA is essential for facilitated diffusion of the TF to locate target DNA site for genetic regulation. Detecting base-pair (bp) resolution of the TF sliding or stepping on the DNA is still experimentally challenging. We have recently performed all-atom molecular dynamics (MD) simulations capturing spontaneous 1-bp stepping of a small WRKY domain TF protein along DNA. Based on the 10 µs WRKY stepping path obtained from such simulations, the protocol here shows how to conduct more extensive conformational samplings of the TF-DNA systems, by constructing the Markov state model (MSM) for the 1-bp protein stepping, with various numbers of micro- and macro-states tested for the MSM construction. In order to examine processive 1-D diffusional search of the TF protein along DNA with structural basis, the protocol further shows how to conduct coarse-grained (CG) MD simulations to sample long-time scale dynamics of the system. Such CG modeling and simulations are particularly useful to reveal the protein-DNA electrostatic impacts on the processive diffusional motions of the TF protein above tens of microseconds, in comparison with sub-microseconds to microseconds protein stepping motions revealed from the all-atom simulations.

Early aggregation mechanism of Aß16-22 revealed by Markov state models.

Aß16-22 is believed to have critical role in early aggregation of full length amyloids that are associated with the Alzheimer's disease and can aggregate to form amyloid fibrils. However, the early aggregation mechanism is still unsolved. Here, multiple long-term molecular dynamics simulations combining with Markov state model were used to probe the early oligomerization mechanism of Aß16-22 peptides. The identified dimeric form adopted either globular random-coil or extended ß-strand like conformations. The observed dimers of these variants shared many overall conformational characteristics but differed in several aspects at detailed level. In all cases, the most common type of secondary structure was intermolecular antiparallel ß-sheets. The inter-state transitions were very frequent ranges from few to hundred nanoseconds. More strikingly, those states which contain fraction of ß secondary structure and significant amount of extended coiled structures, therefore exposed to the solvent, were majorly participated in aggregation. The assembly of low-energy dimers, in which the peptides form antiparallel ß sheets, occurred by multiple pathways with the formation of an obligatory intermediates. We proposed that these states might facilitate the Aß16-22 aggregation through a significant component of the conformational selection mechanism, because they might increase the aggregates population by promoting the inter-chain hydrophobic and the hydrogen bond contacts. The formation of early stage antiparallel ß sheet structures is critical for oligomerization, and at the same time provided a flat geometry to seed the ordered ß-strand packing of the fibrils. Our findings hint at reorganization of this part of the molecule as a potentially critical step in Aß aggregation and will insight into early oligomerization for large ß amyloids.

RQC helical hairpin in Bloom's syndrome helicase regulates DNA unwinding by dynamically intercepting nascent nucleotides.

The RecQ family of helicases are important for maintenance of genomic integrity. Although functions of constructive subdomains of this family of helicases have been extensively studied, the helical hairpin (HH) in the RecQ-C-terminal domain (RQC) has been underappreciated and remains poorly understood. Here by using single-molecule fluorescence resonance energy transfer, we found that HH in the human BLM transiently intercepts different numbers of nucleotides when it is unwinding a double-stranded DNA. Single-site mutations in HH that disrupt hydrogen bonds and/or salt bridges between DNA and HH change the DNA binding conformations and the unwinding features significantly. Our results, together with recent clinical tests that correlate single-site mutations in HH of human BLM with the phenotype of cancer-predisposing syndrome or Bloom's syndrome, implicate pivotal roles of HH in BLM's DNA unwinding activity. Similar mechanisms might also apply to other RecQ family helicases, calling for more attention to the RQC helical hairpin.

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.

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.

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.

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.

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.

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.