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.

TMEFF1 is a neuron-specific restriction factor for herpes simplex virus.

The brain is highly sensitive to damage caused by infection and inflammation1,2. Herpes simplex virus 1 (HSV-1) is a neurotropic virus and the cause of herpes simplex encephalitis3. It is unknown whether neuron-specific antiviral factors control virus replication to prevent infection and excessive inflammatory responses, hence protecting the brain. Here we identify TMEFF1 as an HSV-1 restriction factor using genome-wide CRISPR screening. TMEFF1 is expressed specifically in neurons of the central nervous system and is not regulated by type I interferon, the best-known innate antiviral system controlling virus infections. Depletion of TMEFF1 in stem-cell-derived human neurons led to elevated viral replication and neuronal death following HSV-1 infection. TMEFF1 blocked the HSV-1 replication cycle at the level of viral entry through interactions with nectin-1 and non-muscle myosin heavy chains IIA and IIB, which are core proteins in virus-cell binding and virus-cell fusion, respectively4-6. Notably, Tmeff1-/- mice exhibited increased susceptibility to HSV-1 infection in the brain but not in the periphery. Within the brain, elevated viral load was observed specifically in neurons. Our study identifies TMEFF1 as a neuron-specific restriction factor essential for prevention of HSV-1 replication in the central nervous system.

Cryo-EM structures of a mycobacterial ABC transporter that mediates rifampicin resistance.

Drug-resistant Tuberculosis (TB) is a global public health problem. Resistance to rifampicin, the most effective drug for TB treatment, is a major growing concern. The etiological agent, Mycobacterium tuberculosis (Mtb), has a cluster of ATP-binding cassette (ABC) transporters which are responsible for drug resistance through active export. Here, we describe studies characterizing Mtb Rv1217c-1218c as an ABC transporter that can mediate mycobacterial resistance to rifampicin and have determined the cryo-electron microscopy structures of Rv1217c-1218c. The structures show Rv1217c-1218c has a type V exporter fold. In the absence of ATP, Rv1217c-1218c forms a periplasmic gate by two juxtaposed-membrane helices from each transmembrane domain (TMD), while the nucleotide-binding domains (NBDs) form a partially closed dimer which is held together by four salt-bridges. Adenylyl-imidodiphosphate (AMPPNP) binding induces a structural change where the NBDs become further closed to each other, which downstream translates to a closed conformation for the TMDs. AMPPNP binding results in the collapse of the outer leaflet cavity and the opening of the periplasmic gate, which was proposed to play a role in substrate export. The rifampicin-bound structure shows a hydrophobic and periplasm-facing cavity is involved in rifampicin binding. Phospholipid molecules are observed in all determined structures and form an integral part of the Rv1217c-1218c transporter system. Our results provide a structural basis for a mycobacterial ABC exporter that mediates rifampicin resistance, which can lead to different insights into combating rifampicin resistance.

Structure of the priming arabinosyltransferase AftA required for AG biosynthesis of Mycobacterium tuberculosis.

Arabinogalactan (AG) is an essential cell wall component in mycobacterial species, including the deadly human pathogen Mycobacterium tuberculosis. It plays a pivotal role in forming the rigid mycolyl-AG-peptidoglycan core for in vitro growth. AftA is a membrane-bound arabinosyltransferase and a key enzyme involved in AG biosynthesis which bridges the assembly of the arabinan chain to the galactan chain. It is known that AftA catalyzes the transfer of the first arabinofuranosyl residue from the donor decaprenyl-monophosphoryl-arabinose to the mature galactan chain (i.e., priming); however, the priming mechanism remains elusive. Herein, we report the cryo-EM structure of Mtb AftA. The detergent-embedded AftA assembles as a dimer with an interface maintained by both the transmembrane domain (TMD) and the soluble C-terminal domain (CTD) in the periplasm. The structure shows a conserved glycosyltransferase-C fold and two cavities converging at the active site. A metal ion participates in the interaction of TMD and CTD of each AftA molecule. Structural analyses combined with functional mutagenesis suggests a priming mechanism catalyzed by AftA in Mtb AG biosynthesis. Our data further provide a unique perspective into anti-TB drug discovery.

Allele-Specific Suppression of Variant MHC With High-Precision RNA Nuclease CRISPR-Cas13d Prevents Hypertrophic Cardiomyopathy.

BACKGROUND: Familial hypertrophic cardiomyopathy has severe clinical complications of heart failure, arrhythmia, and sudden cardiac death. Heterozygous single nucleotide variants (SNVs) of sarcomere genes such as MYH7 are the leading cause of this type of disease. CRISPR-Cas13 (clustered regularly interspaced short palindromic repeats and their associated protein 13) is an emerging gene therapy approach for treating genetic disorders, but its therapeutic potential in genetic cardiomyopathy remains unexplored. METHODS: We developed a sensitive allelic point mutation reporter system to screen the mutagenic variants of Cas13d. On the basis of Cas13d homology structure, we rationally designed a series of Cas13d variants and obtained a high-precision Cas13d variant (hpCas13d) that specifically cleaves the MYH7 variant RNAs containing 1 allelic SNV. We validated the high precision and low collateral cleavage activity of hpCas13d through various in vitro assays. We generated 2 HCM mouse models bearing distinct MYH7 SNVs and used adenovirus-associated virus serotype 9 to deliver hpCas13d specifically to the cardiomyocytes. We performed a large-scale library screening to assess the potency of hpCas13d in resolving 45 human MYH7 allelic pathogenic SNVs. RESULTS: Wild-type Cas13d cannot distinguish and specifically cleave the heterozygous MYH7 allele with SNV. hpCas13d, with 3 amino acid substitutions, had minimized collateral RNase activity and was able to resolve various human MYH7 pathological sequence variations that cause hypertrophic cardiomyopathy. In vivo application of hpCas13d to 2 hypertrophic cardiomyopathy models caused by distinct human MYH7 analogous sequence variations specifically suppressed the altered allele and prevented cardiac hypertrophy. CONCLUSIONS: Our study unveils the great potential of CRISPR-Cas nucleases with high precision in treating inheritable cardiomyopathy and opens a new avenue for therapeutic management of inherited cardiac diseases.

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.

Antifouling Asymmetric Block Copolymer Nanofilms via Freestanding Interfacial Polymerization for Efficient and Sustainable Water Purification.

Membrane materials that resist nonspecific or specific adsorption are urgently required in widespread practical applications, such as water purification, food processing, and life sciences. In water purification, inevitable membrane fouling not only limits membrane separation performance, leading to a decline in both permeance and selectivity, but also remarkably increases operation requirements, and augments extra maintenance costs and higher energy consumption. In this work, we report a freestanding interfacial polymerization (IP) fabrication strategy for in situ creation of asymmetric block copolymer (BCP) nanofilms with antifouling properties, greatly outperforming the conventional surface post-modification approaches. The resultant free-standing asymmetric BCP nanofilms with highly-dense, highly-hydrophilic polyethylene glycol (PEG) brushes on one side, can be readily formed via a typical IP process of a well-defined double-hydrophilic BCP composed of a highly-efficient antifouling PEG block and a membrane-forming multiamine block. The asymmetric BCP nanofilms have been applied for efficient and sustainable natural water purification, demonstrating extraordinary antifouling capabilities accompanied with superior separation performance far beyond commercial polyamide nanofiltration membranes. The antifouling behaviors of asymmetric BCP nanofilms derived from the combined effect of the hydration layer, electrostatic repulsion and steric hindrance were further elucidated by water flux and fouling resistance in combination with all-atom molecular dynamics (MD) simulation. This work opens up a new avenue for the large-scale and low-cost creation of broad-spectrum, asymmetric membrane materials with diverse functional "defect-free" surfaces in real-world applications.

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.