Search

Selective serotonin reuptake inhibitors (SSRIs) have shown promise in cancer therapy, particularly for hepatocellular carcinoma (HCC), but their molecular targets and mechanisms remain unclear. Here, we show that SSRIs exhibit significant anti-HCC effects independent of their classical target, the serotonin reuptake transporter (SERT). Using global inverse gene expression profiling, drug affinity responsive target stability assays, and in silico molecular docking, we demonstrate that citalopram targets glucose transporter 1 (GLUT1), resulting in reduced glycolytic flux. A mutant GLUT1 variant at the citalopram binding site (E380) diminishes the drug's inhibitory effects on the Warburg effect and tumor growth. In preclinical models, citalopram dampens the growth of GLUT1 high liver tumors and displays a synergistic effect with anti-PD-1 therapy. Retrospective analysis reveals that SSRI use correlates with a lower risk of metastasis among patients with HCC. Our study describes a role for SSRIs in cancer metabolism, establishing a rationale for their repurposing as potential anti-cancer drugs for HCC., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

The molecular mechanisms underlying epitope recognition by T cell receptors (TCRs) are critical for activating T cell immune responses and rationally designing TCR-based therapeutics. Single-cell sequencing techniques vastly boost the accumulation of TCR sequences, while the limitation of available TCR-pMHC structures hampers further investigations. In this study, we proposed a computational pipeline that incorporates structural information and single-cell sequencing data to investigate the epitope-recognition mechanisms for TCRs without identified epitopes. By antigen specificity clustering, we mapped the epitope sequences between epitope-known and epitope-unknown TCRs from COVID-19 patients. One reported SARS-CoV-2 epitope, NQKLIANQF (S 919 - 927 ), was identified for a TCR expressed by 614 T cells (TCR-614). Epitope screening also identified a potential cross-reactive epitope, KLKTLVATA (NSP3 1790 - 1798 ), for a TCR expressed by 204 T cells (TCR-204). By molecular dynamics (MD) simulations, we revealed the detailed epitope-recognition mechanisms for both TCRs. The structural motifs responsible for epitope recognition revealed by the MD simulations are consistent with the sequential features recognized by the sequence-based clustering method. We hope that this strategy could facilitate the discovery and optimization of TCR-based therapeutics., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

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.

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.

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., Competing Interests: Declaration of competing interest The authors declare no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

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., (© 2023 The Authors. Advanced Science published by Wiley-VCH GmbH.)

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 ( 215 DVYM 218 ) 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., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

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., Competing Interests: Declaration of interests Shanghai Institute of Biochemistry and Cell Biology and Shanghai Jiao Tong University School of Medicine have filed a patent application based on the findings in this article., (Copyright © 2023 Elsevier Inc. All rights reserved.)

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 ( 235 ITRFQTLLALHRSYL 249 ). 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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Declaration of competing interest The authors report no conflicts of interest for this work., (Copyright © 2023. Published by Elsevier Inc.)

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.

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.

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.

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.

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., (Copyright © 2022 Elsevier B.V. All rights reserved.)

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.

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., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

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., (Copyright © 2021 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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., Competing Interests: Declaration of competing interest The authors have declared that no competing interests exist., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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., (© 2020 Federation of American Societies for Experimental Biology.)

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.

XimA is a unique amide synthetase that belongs to the ANL superfamily of adenylating enzymes, but with a special structural fold. In order to improve the enzyme promiscuity, we engineered XimA by site-directed mutagenesis at a specific position based on our theoretical model of XimA. Thus, we were able to produce diverse benzopyran derivatives with up to 15 different l-form and d-form amino acid substitutions, catalyzed by several XimA variants. Molecular docking and molecular dynamics simulations conducted for various XimA systems provide further structural insights into the substitution effects of the phenylalanine-201 as an active site residue on protein dynamics and enzyme catalysis.

HIV-1 entry is mediated firstly by the molecular recognition between the viral glycoprotein gp120 and its receptor CD4 on host T-cells. As a key antigen that can be targeted by neutralizing antibodies, gp120 has been a focus for extensive studies with efforts to understand its structural properties and conformational dynamics upon receptor binding. An atomistic-level revelation of gp120 opening dynamics activated by CD4, however, is still unknown. Here, by constructing a Markov State Model (MSM) based on hundreds of Molecular Dynamics (MD) simulations with an aggregated simulation time of ∼20 microseconds (μs), we identify the key metastable states of gp120 during its opening dynamics upon CD4 binding. The MSM provides a clear dynamic model whereby the identified metastable states coexist and can reach an equilibrium. More importantly, a hydrophobic core flanked by variable loops (V1V2 and V3) and the β20/21 region plays an essential role in triggering the gp120 opening. Any destabilizing effects introduced into the hydrophobic core, therefore, can be expected to promote transition of gp120 to an open state. Moreover, the variable loops demonstrate high flexibilities in fully open gp120. In particular, the V3 region is capable of exploring both closed and open conformations, even with the V1/V2 loops largely adopting an open form. In addition, the bridging sheet formation in gp120 is likely induced by the incoming co-receptor/antibody recognitions, since the V1/V2 structure is highly heterogeneous so that the bridging-sheet formed conformation is not the most populated state. Our studies provide deep insights into the dynamic features of gp120 and its molecular recognitions to the broadly neutralizing antibodies, which guides future attempts to design more effective gp120 immunogens.

The known benzonaphthyridine alkaloid, albogrisin A ( 1 ), and six new compounds, including two pyrazinone stereoisomers, albogrisin B ( 2 )/B' ( 2' ), together with four 4 H -pyrroloquinolinones, two diastereoisomers, albogrisin C ( 3 )/C' ( 3' ), and their methyl esters, albogrisin D ( 4 )/D' ( 4' ), were isolated from mangrove-derived Streptomyces albogriseolus MGR072. 2 and 2' are converted into 1 in acidic aqueous solution but into 3 / 3' and 4 / 4' in 0.05% trifluoroacetic acid acetonitrile. 4 and 4' are new indoleamine 2,3-dioxygenase 1 inhibitors.

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., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The conformational transition of prion protein (PrP) from a native form PrP C to a pathological isoform PrP Sc is the main cause of a number of prion diseases in human and animals. Thus, understanding the molecular basis of conformational transition of PrP will be valuable for unveiling the etiology of PrP-related diseases. Here, to explore the potential misfolding mechanism of PrP under the acidic condition, which is known to promote PrP misfolding and trigger its aggregation, the conventional and accelerated molecular dynamics (MD) simulations combined with the Markov state model (MSM) analysis were performed. The conventional MD simulations reveal that, at an acidic pH, the globular domain of PrP is partially unfolded, particularly for the α2 C-terminus. Structural analysis of the key macrostates obtained by MSM indicates that the α2 C-terminus and the β2-α2 loop may serve as important sites for the pH-induced PrP misfolding. Meanwhile, the α1 may also participate in the pH-induced structural conversion by moving away from the α2-α3 subdomain. Notably, dynamical network analysis of the key metastable states indicates that the protonated H187 weakens the interactions between the α2 C-terminus, α1-β2 loop, and α2-α3 loop, leading these domains, especially the α2 C-terminus, to become unstable and to begin to misfold. Therefore, the α2 C-terminus plays a key role in the PrP misfolding process and serves as a potential site for drug targeting. Overall, our findings can deepen the understanding of the pathogenesis related to PrP and provide useful guidance for the future drug discovery.

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., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA polymerase (RNAP) from bacteriophage T7 is a representative single-subunit viral RNAP that can transcribe with high promoter activities without assistances from transcription factors. We accordingly studied this small transcription machine computationally as a model system to understand underlying mechanisms of mechano-chemical coupling and fidelity control in the RNAP transcription elongation. Here we summarize our computational work from several recent publications to demonstrate first how T7 RNAP translocates via Brownian alike motions along DNA right after the catalytic product release. Then we show how the backward translocation motions are prevented at post-translocation upon successful nucleotide incorporation, which is also subject to stepwise nucleotide selection and acts as a pawl for "selective ratcheting". The structural dynamics and energetics features revealed from our atomistic molecular dynamics (MD) simulations and related analyses on the single-subunit T7 RNAP thus provided detailed and quantitative characterizations on the Brownian-ratchet working scenario of a prototypical transcription machine with sophisticated nucleotide selectivity for fidelity control. The presented mechanisms can be more or less general for structurally similar viral or mitochondrial RNAPs and some of DNA polymerases, or even for the RNAP engine of the more complicated transcription machinery in higher organisms.

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.

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., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Here, we studied the complete process of a viral T7 RNA polymerase (RNAP) translocation on DNA during transcription elongation by implementing extensive all-atom molecular dynamics (MD) simulations to construct a Markov state model (MSM). Our studies show that translocation proceeds in a Brownian motion, and the RNAP thermally transits among multiple metastable states. We observed non-synchronized backbone movements of the nucleic acid (NA) chains with the RNA translocation accomplished first, while the template DNA lagged. Notably, both the O-helix and Y-helix on the fingers domain play key roles in facilitating NA translocation through the helix opening. The helix opening allows a key residue Tyr639 to become inserted into the active site, which pushes the RNA-DNA hybrid forward. Another key residue, Phe644, coordinates the downstream template DNA motions by stacking and un-stacking with a transition nucleotide (TN) and its adjacent nucleotide. Moreover, the O-helix opening at pre-translocation (pre-trans) likely resists backtracking. To test this hypothesis, we computationally designed mutants of T7 RNAP by replacing the amino acids on the O-helix with counterpart residues from a mitochondrial RNAP that is capable of backtracking. The current experimental results support the hypothesis., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The dynamics of the RNA polymerase II (Pol II) backtracking process is poorly understood. We built a Markov State Model from extensive molecular dynamics simulations to identify metastable intermediate states and the dynamics of backtracking at atomistic detail. Our results reveal that Pol II backtracking occurs in a stepwise mode where two intermediate states are involved. We find that the continuous bending motion of the Bridge helix (BH) serves as a critical checkpoint, using the highly conserved BH residue T831 as a sensing probe for the 3'-terminal base paring of RNA:DNA hybrid. If the base pair is mismatched, BH bending can promote the RNA 3'-end nucleotide into a frayed state that further leads to the backtracked state. These computational observations are validated by site-directed mutagenesis and transcript cleavage assays, and provide insights into the key factors that regulate the preferences of the backward translocation.

The aim of this study was to evaluate Eucalyptus badjensis concerning the genetic variation for growth traits and the potential of the species in supporting a breeding programme. The field trial was a provenance/progeny test established in Campina da Alegria, Santa Catarina, Brazil (latitude 26°52'05.1″ S, longitude 51°48'47.5″ W, altitude 1,015 m) in a soil classified as Latossolic Alumino-Ferric Brown Nitosol. The experiment comprised 60 open-pollinated progenies from the provenances Glenbog and Badja State Forest, New South Wales, Australia. Ten replicates and plots with six plants in row were used. At the age of 17 years, 279 trees were assessed for diameter of the bole at breast height (DBH), total tree height (H) and volume of wood with bark (Vol). After submitting the data to statistical genetic analysis, the overall means for DBH, H and Vol were 45.17 cm, 33.30 m and 2.84 m3, and the estimates of additive coefficient of variation [ C V a(%)] were 12.59%, 5.91% and 26.51%, respectively. Heritability coefficients of additive effects ( h a 2) were also estimated and the following values were found: 0.443, 0.312 and 0.358. Thirty-nine trees from 25 different progenies were selected. The expected means of the provenances after improvement were 50.02 cm, 34.35 m and 3.47 m3 for DBH, H and Vol, respectively.

Pyrophosphate ion (PPi) release during transcription elongation is a signature step in each nucleotide addition cycle. The kinetics and energetics of the process as well as how it proceeds with substantial conformational changes of the polymerase complex determine the mechano-chemical coupling mechanism of the transcription elongation. Here we investigated detailed dynamics of the PPi release process in a single-subunit RNA polymerase (RNAP) from bacteriophage T7, implementing all-atom molecular dynamics (MD) simulations. We obtained a jump-from-cavity kinetic model of the PPi release utilizing extensive nanosecond MD simulations. We found that the PPi release in T7 RNAP is initiated by the PPi dissociation from two catalytic aspartic acids, followed by a comparatively slow jump-from-cavity activation process. Combining with a number of microsecond long MD simulations, we also found that the activation process is hindered by charged residue associations as well as by local steric and hydrogen bond interactions. On the other hand, the activation is greatly assisted by a highly flexible lysine residue Lys472 that swings its side chain to pull PPi out. The mechanism can apply in general to single subunit RNA and DNA polymerases with similar molecular structures and conserved key residues. Remarkably, the flexible lysine or arginine residue appears to be a universal module that assists the PPi release even in multi-subunit RNAPs with charge facilitated hopping mechanisms. We also noticed that the PPi release is not tightly coupled to opening motions of an O-helix on the fingers domain of T7 RNAP according to the microsecond MD simulations. Our study thus supports the Brownian ratchet scenario of the mechano-chemical coupling in the transcription elongation of the single-subunit polymerase.

The RNA polymerase II elongation is central in eukaryotic transcription. Although multiple intermediates of the elongation complex have been identified, the dynamical mechanisms remain elusive or controversial. Here we build a structure-based kinetic model of a full elongation cycle of polymerase II, taking into account transition rates and conformational changes characterized from both single molecule experimental studies and computational simulations at atomistic scale. Our model suggests a force-dependent slow transition detected in the single molecule experiments corresponds to an essential conformational change of a trigger loop (TL) opening prior to the polymerase translocation. The analyses on mutant study of E1103G and on potential sequence effects of the translocation substantiate this proposal. Our model also investigates another slow transition detected in the transcription elongation cycle which is independent of mechanical force. If this force-independent slow transition happens as the TL gradually closes upon NTP binding, the analyses indicate that the binding affinity of NTP to the polymerase has to be sufficiently high. Otherwise, one infers that the slow transition happens pre-catalytically but after the TL closing. Accordingly, accurate determination of intrinsic properties of NTP binding is demanded for an improved characterization of the polymerase elongation. Overall, the study provides a working model of the polymerase II elongation under a generic Brownian ratchet mechanism, with most essential structural transition and functional kinetics elucidated.

Nucleotide selection is essential for fidelity control in gene replication and transcription. Recent work on T7 RNA polymerase suggested that a small posttranslocation free energy bias stabilizes Tyr(639) in the active site to aid nucleotide selection. However, it was not clear exactly how Tyr(639) assists the selection. Here we report a molecular-dynamics simulation study revealing atomistic detail of this critical selectivity. The study shows first that Tyr(639) blocks the active site at posttranslocation by marginally stacking to the end basepair of the DNA-RNA hybrid. The study then demonstrates that at the nucleotide preinsertion state, a cognate RNA nucleotide does not affect the local Tyr(639) stabilization, whereas a noncognate nucleotide substantially stabilizes Tyr(639) so that Tyr(639) keeps blocking the active site. As a result, further nucleotide insertion into the active site, which requires moving Tyr(639) out of the site, would be hindered for the noncognate nucleotide, but not for the cognate nucleotide. In particular, we note that water molecules assist the ribose recognition in the RNA nucleotide preinsertion, and help Tyr(639) stacking to the end basepair in the case of a DNA nucleotide. It was also seen that a base-mismatched nucleotide at preinsertion directly grabs Tyr(639) for the active site stabilization. We also find that in a mutant polymerase Y639F the strong stabilization of residue 639 in the active site cannot establish upon the DNA nucleotide preinsertion. The finding explains the reduced differentiation between ribo- and deoxyribonucleotides that has been recorded experimentally for the mutant polymerase.

The aim of this study was to assess the prevalence and risk factors for faecal carriage of extended-spectrum β-lactamase (ESBL) and plasmidic AmpC β-lactamase (pAmpC) Escherichia coli producers in dogs. A three-month cross-sectional study was conducted and 151 rectal swabs were obtained from healthy dogs. ESBL and pAmpC genes were detected by PCR and were sequenced. Logistic regression models were used to investigate risk factors for the carriage of ESBL and pAmpC-producing E. coli. About 15 per cent of the isolates carried ESBL genes (blaCTX-M-32 n=8, blaCTX-M-15 n=5, blaCTX-M-1 n=3, blaCTX-M-9-like n=4) and 20 per cent carried pAmpC genes (blaCMY-2 n=23, blaCMY-2-like n=2). Thirteen dogs carried an E. coli isolate with both an ESBL and a pAmpC gene. One E. coli isolate harboured the human blaDHA-1 pAmpC gene, which has not been previously reported in companion animals in Europe. Dogs with a history of antimicrobial therapy in the past year had a higher risk of being carriers of ESBL-producing (P=0.003, OR =7.85) and pAmpC-producing (P=0.005, OR=6.28) E. coli. Dogs from shelter/breeders were approximately three times more likely to have an ESBL- or a pAmpC-producing E. coli than dogs from private owners. Males have a reduced risk of carrying a pAmpC-producing E. coli than females (P=0.017, OR =0.28). The knowledge of potential risk factors may help to limit the impact of resistance through implementation of effective control measures and judicious antimicrobial therapy., (British Veterinary Association.)

Transcription is a central step in gene expression, in which the DNA template is processively read by RNA polymerase II (Pol II), synthesizing a complementary messenger RNA transcript. At each cycle, Pol II moves exactly one register along the DNA, a process known as translocation. Although X-ray crystal structures have greatly enhanced our understanding of the transcription process, the underlying molecular mechanisms of translocation remain unclear. Here we use sophisticated simulation techniques to observe Pol II translocation on a millisecond timescale and at atomistic resolution. We observe multiple cycles of forward and backward translocation and identify two previously unidentified intermediate states. We show that the bridge helix (BH) plays a key role accelerating the translocation of both the RNA:DNA hybrid and transition nucleotide by directly interacting with them. The conserved BH residues, Thr831 and Tyr836, mediate these interactions. To date, this study delivers the most detailed picture of the mechanism of Pol II translocation at atomic level.

Conformational changes of proteins are an*Author contributed equally with all other contributors. essential part of many biological processes such as: protein folding, ligand binding, signal transduction, allostery, and enzymatic catalysis. Molecular dynamics (MD) simulations can describe the dynamics of molecules at atomic detail, therefore providing a much higher temporal and spatial resolution than most experimental techniques. Although MD simulations have been widely applied to study protein dynamics, the timescales accessible by conventional MD methods are usually limited to timescales that are orders of magnitude shorter than the conformational changes relevant for most biological functions. During the past decades great effort has been devoted to the development of theoretical methods that may enhance the conformational sampling. In recent years, it has been shown that the statistical mechanics framework provided by discrete-state and -time Markov State Models (MSMs) can predict long timescale dynamics from a pool of short MD simulations. In this chapter we provide the readers an account of the basic theory and selected applications of MSMs. We will first introduce the general concepts behind MSMs, and then describe the existing procedures for the construction of MSMs. This will be followed by the discussions of the challenges of constructing and validating MSMs, Finally, we will employ two biologically-relevant systems, the RNA polymerase and the LAO-protein, to illustrate the application of Markov State Models to elucidate the molecular mechanisms of complex conformational changes at biologically relevant timescales.

Aims: HIV-1 expanded in an allogenic system (Al-HIV) represents a cheaper and faster alternative to the autologous virus (Au-HIV) as an antigen in anti-HIV immunotherapy. In this study, chemically inactivated HIV-1 obtained through autologous or allogenic systems were compared., Patients & Methods: Au-HIV and Al-HIV obtained from cultures of peripheral blood mononuclear cells from 11 HIV(+) individuals were tested for virus production, yield and time of culture, and their ability to elicit a specific immune response in vitro., Results: The allogenic system was more efficient than the autologous system. Dendritic cells pulsed with Au-HIV and Al-HIV presented a similar phenotypic profile, but only Al-HIV induced a significant increase in IFN-γ(+) lymphocytes., Conclusion: The use of an allogenic system displays several advantages in terms of cell manipulation, time and cost of culture, and immunogenicity.

The dynamics of the PPi release during the transcription elongation of bacterial RNA polymerase and its effects on the Trigger Loop (TL) opening motion are still elusive. Here, we built a Markov State Model (MSM) from extensive all-atom molecular dynamics (MD) simulations to investigate the mechanism of the PPi release. Our MSM has identified a simple two-state mechanism for the PPi release instead of a more complex four-state mechanism observed in RNA polymerase II (Pol II). We observed that the PPi release in bacterial RNA polymerase occurs at sub-microsecond timescale, which is ∼3-fold faster than that in Pol II. After escaping from the active site, the (Mg-PPi)(2-) group passes through a single elongated metastable region where several positively charged residues on the secondary channel provide favorable interactions. Surprisingly, we found that the PPi release is not coupled with the TL unfolding but correlates tightly with the side-chain rotation of the TL residue R1239. Our work sheds light on the dynamics underlying the transcription elongation of the bacterial RNA polymerase.

The last years of research have been particularly dynamic in establishing the importance of peptide deformylase (PDF), a protein of the N-terminal methionine excision (NME) pathway that removes formyl-methionine from mitochondrial-encoded proteins. The genomic sequence of the human PDF gene is shared with the COG8 gene, which encodes a component of the oligomeric golgi complex, a very unusual case in Eukaryotic genomes. Since PDF is crucial in maintaining mitochondrial function and given the atypical short distance between the end of COG8 coding sequence and the PDF initiation codon, we investigated whether the regulation of the human PDF is affected by the COG8 overlapping partner. Our data reveals that PDF has several transcription start sites, the most important of which only 18 bp from the initiation codon. Furthermore, luciferase-activation assays using differently-sized fragments defined a 97 bp minimal promoter region for human PDF, which is capable of very strong transcriptional activity. This fragment contains a potential Sp1 binding site highly conserved in mammalian species. We show that this binding site, whose mutation significantly reduces transcription activation, is a target for the Sp1 transcription factor, and possibly of other members of the Sp family. Importantly, the entire minimal promoter region is located after the end of COG8's coding region, strongly suggesting that the human PDF preserves an independent regulation from its overlapping partner., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Pyrophosphate ion (PP(i)) release after nucleotide incorporation is a necessary step for RNA polymerase II (pol II) to enter the next nucleotide addition cycle during transcription elongation. However, the role of pol II residues in PP(i) release and the mechanistic relationship between PP(i) release and the conformational change of the trigger loop remain unclear. In this study, we constructed a Markov state model (MSM) from extensive all-atom molecular dynamics (MD) simulations in the explicit solvent to simulate the PP(i) release process along the pol II secondary channel. Our results show that the trigger loop has significantly larger intrinsic motion after catalysis and formation of PP(i), which in turn aids PP(i) release mainly through the hydrogen bonding between the trigger loop residue H1085 and the (Mg-PP(i))(2-) group. Once PP(i) leaves the active site, it adopts a hopping model through several highly conserved positively charged residues such as K752 and K619 to release from the pol II pore region of the secondary channel. These positive hopping sites form favorable interactions with PP(i) and generate four kinetically metastable states as identified by our MSM. Furthermore, our single-mutant simulations suggest that H1085 and K752 aid PP(i) exit from the active site after catalysis, whereas K619 facilitates its passage through the secondary channel. Finally, we suggest that PP(i) release could help the opening motion of the trigger loop, even though PP(i) release precedes full opening of the trigger loop due to faster PP(i) dynamics. Our simulations provide predictions to guide future experimental tests., (© 2011 American Chemical Society)

Amyloid fibrillation of proteins is associated with a great variety of pathologic conditions. Development of new molecules that can monitor amyloidosis kinetics and inhibit fibril formation is of great diagnostic and therapeutic value. In this work, we have developed a biocompatible molecule that functions as an ex situ monitor and an in situ inhibitor for protein fibrillation, using insulin as a model protein. 1,2-Bis[4-(3-sulfonatopropoxyl)phenyl]-1,2-diphenylethene salt (BSPOTPE) is nonemissive when it is dissolved with native insulin in an incubation buffer but starts to fluoresce when it is mixed with preformed insulin fibril, enabling ex situ monitoring of amyloidogenesis kinetics and high-contrast fluorescence imaging of protein fibrils. Premixing BSPOTPE with insulin, on the other hand, inhibits the nucleation process and impedes the protofibril formation. Increasing the dose of BSPOTPE boosts its inhibitory potency. Theoretical modeling using molecular dynamics simulations and docking reveals that BSPOTPE is prone to binding to partially unfolded insulin through hydrophobic interaction of the phenyl rings of BSPOTPE with the exposed hydrophobic residues of insulin. Such binding is assumed to have stabilized the partially unfolded insulin and obstructed the formation of the critical oligomeric species in the protein fibrillogenesis process., (© 2011 American Chemical Society)

A recently discovered small-molecule inhibitor, BMS-488043 (BMS-488), for the invasion of Human immunodeficiency virus Type 1 (HIV-1), shows a high activity against the entry of diversified HIV-1. Docking and molecular dynamic studies have been carried out to understand the binding mode of BMS-488 to gp120 as well as the effect of the small molecule on the conformational change of gp120 induced by CD4 binding. The results indicate that BMS-488 can accommodate in the CD4 binding pocket and interfere the CD4 binding in a noncompetitive mode. The piperazine group of BMS-488 prevents the bridging sheet formation of gp120 induced by the CD4 binding mainly through blocking the rotation of the Trp112 located on the α1 helix of gp120. The bridging sheet formation cannot be blocked for the W112A mutant of gp120 due to the reduced steric hindrance, in agreement with its significant resistance to the BMS inhibitor. The aza-indole ring is likely to interfere the exposure of gp41 by stacking within the β3-β5 and LB loops to disrupt the close packing of Pro212-His66-Phe210. The mode of action of BMS-488 also accommodates many mutagenesis results related to BMS-488 activity., (Copyright © 2011 Wiley-Liss, Inc.)