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

Revealing the conformational dynamics of UDP-GlcNAc recognition by O-GlcNAc transferase via Markov state model

The O-linked N-acetylglucosamine (O-GlcNAc) glycosylation is a critical post-translational modification and closely linked to various physiological and pathological conditions. The O-GlcNAc transferase (OGT) functions as the only glycosyltransferase of O-GlcNAc glycosylation by transferring GlcNAc from UDP-GlcNAc to serine or threonine residues on protein substrates. The interaction mode of UDP-GlcNAc against OGT has been preliminarily revealed by the crystal structures, yet an atomic-level comprehension for the conformational dynamics of the recognition process remains elusive. Here, we construct the Markov state model based on extensive all-atom molecular dynamics (MD) simulations with an aggregated simulation time of ∼9 μs, and reveal that the UDP-GlcNAc recognition process by OGT encompasses four key metastable states, occurring within an estimated timescale of ∼10 μs. During UDP-GlcNAc recognition process, we find the pyrophosphate moiety (P2O52-) initially anchors to the active pocket via salt bridge and hydrogen bonds, facilitating subsequent binding of the uridine and GlcNAc moieties. Furthermore, the functional roles of K842 involved in the salt bridge with P2O52- were evaluated through extra mutant MD simulations. Overall, our study provides valuable insights into the UDP-GlcNAc recognition mechanism by OGT, which could further aid in mechanistic studies of O-GlcNAc glycosylation and drug development targeting on OGT.

All-or-None Selectivity in Probing Polarity-Determined Trinucleotide Repeat Foldings with a Parity Resolution by a Beyond-Size-Matching Ligand

Abnormal amplification of trinucleotide repeats (TNRs) is associated with neurodegenerative diseases by forming a particular hairpin bulge. It is well known that the polarity and parity of TNRs can regulate the formed hairpin structures. Therefore, there is a great challenge to efficiently discriminate the hairpin structures of TNRs with substantial selectivity. Herein, we developed a fluorescent ligand of pseudohypericin (Pse) with a beyond-size-matching (BSM) geometry to selectively sense hairpin structures of GTC and CTG TNRs. The GTC hairpin structures can bind with Pse dominantly at extreme T-T mismatches by the virtue of their most extrahelical conformations, while there is no binding event to occur with the polarity-inverted counterpart CTG hairpin structures because of the limited space provided by their intrahelical T-T mismatches. In addition, this all-or-none response with the polarity-dependent folding (PoDF) is independent of the length of these TNRs. Interestingly, the parity-dependent folding (PaDF) of GTC hairpin structures can also be resolved. Besides pure TNRs, the competency of this BSM ligand to sense the PoDF and PaDF effects was also generalized to DNAs with TNRs occurring at loop and stem end regions. To our knowledge, this is the first experimental observation with the state-of-the-art performance over the fluorescence measurement of PoDF and PaDF in TNRs. Our work provides an expedient way to elucidate the TNR folding by designing ligands having BSM features.

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.

Automated Path Searching Reveals the Mechanism of Hydrolysis Enhancement by T4 Lysozyme Mutants

Bacteriophage T4 lysozyme (T4L) is a glycosidase that is widely applied as a natural antimicrobial agent in the food industry. Due to its wide applications and small size, T4L has been regarded as a model system for understanding protein dynamics and for large-scale protein engineering. Through structural insights from the single conformation of T4L, a series of mutations (L99A,G113A,R119P) have been introduced, which have successfully raised the fractional population of its only hydrolysis-competent excited state to 96%. However, the actual impact of these substitutions on its dynamics remains unclear, largely due to the lack of highly efficient sampling algorithms. Here, using our recently developed travelling-salesman-based automated path searching (TAPS), we located the minimum-free-energy path (MFEP) for the transition of three T4L mutants from their ground states to their excited states. All three mutants share a three-step transition: the flipping of F114, the rearrangement of α0/α1 helices, and final refinement. Remarkably, the MFEP revealed that the effects of the mutations are drastically beyond the expectations of their original design: (a) the G113A substitution not only enhances helicity but also fills the hydrophobic Cavity I and reduces the free energy barrier for flipping F114; (b) R119P barely changes the stability of the ground state but stabilizes the excited state through rarely reported polar contacts S117_OG:N132_ND2, E11_OE1:R145_NH1, and E11_OE2:Q105_NE2; (c) the residue W138 flips into Cavity I and further stabilizes the excited state for the triple mutant L99A,G113A,R119P. These novel insights that were unexpected in the original mutant design indicated the necessity of incorporating path searching into the workflow of rational protein engineering.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

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.

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.

Kinetic Analysis of the Effect of N-Terminal Acetylation on Thymine DNA Glycosylase

Thymine DNA glycosylase (TDG) is tasked with initiating DNA base excision repair by recognizing and removing T, U, the chemotherapeutic 5-fluorouracil (5-FU), and many other oxidized and halogenated pyrimidine bases. TDG contains a long, unstructured N-terminus that contains four known sites of acetylation: lysine (K) residues 59, 83, 84, and 87. Here, K to glutamine (Q) mutants are used as acetyl-lysine (AcK) analogues to probe the effect of N-terminal acetylation on the kinetics of TDG. We find that mimicking acetylation affects neither the maximal single-turnover rate k_max nor the turnover rate k_TO, indicating that the steps after initial binding, through chemistry and product release, are not affected. Under subsaturating conditions, however, acetylation changes the processing of U substrates. Subtle differences among AcK analogues are revealed with 5-FU in single-stranded DNA. We propose that the subtleties observed among the AcK analogues may be amplified on the genomic scale, leading to regulation of TDG activity. N-terminal acetylation, though, may also play a structural, rather than kinetic role in vivo.

Kinetics and thermodynamics of BI-BII interconversion altered by T:G mismatches in DNA

T:G mismatches in DNA result in humans primarily from deamination of methylated CpG sites. They are repaired by redundant systems, such as thymine DNA glycosylase (TDG) and methyl-binding domain enzyme (MBD4), and maintenance of these sites has been implicated in epigenetic processes. The process by which these enzymes identify a canonical DNA base in the incorrect basepairing context remains a mystery. However, the conserved contacts of the repair enzymes with the DNA backbone suggests a role for protein-phosphate interaction in the recognition and repair processes. We have used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion, and for this work have focused on alterations to the activation barriers to interconversion and the effect of a mismatch compared with canonical DNA. We have found that alterations to the ΔG of interconversion for T:G basepairs are remarkably similar to U:G basepairs in the form of stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA, suggesting a universality of this result for TDG substrates. Likewise, we see perturbations to the free energy (∼1 kcal/mol) and enthalpy (2-5 kcal/mol) of activation for the BI-BII interconversion localized to the phosphates flanking the mismatch. Overall our results strongly suggest that the perturbed backbone energetics in T:G basepairs play a significant role in the recognition process of DNA repair enzymes.

The rate of formation and stability of abasic site interstrand crosslinks in the DNA duplex

DNA interstrand crosslinks (ICLs) strands pose an impenetrable barrier for DNA replication. Different ICLs are known to recruit distinct DNA repair pathways. NEIL3 glycosylase has been known to remove an abasic (Ap) site derived DNA crosslink (Ap-ICL). An Ap-ICL forms spontaneously from the Ap site with an adjacent adenine in the opposite strand. Lack of genetic models and a poor understanding of the fate of these lesions leads to many questions about the occurrence and the toxicity of Ap-ICL in cells. Here, we investigate the circumstances of Ap-ICL formation. With an array of different oligos, we have investigated the rates of formation, the yields, and the stability of Ap-ICL. Our findings point out how different bases in the vicinity of the Ap site change crosslink formation in vitro. We reveal that AT-rich rather than GC-rich regions in the surrounding Ap site lead to higher rates of Ap-ICL formation. Overall, our data reveal that Ap-ICL can be formed in virtually any DNA sequence context surrounding a hot spot of a 5'-Ap-dT pair, albeit with significantly different rates and yields. Based on Ap-ICL formation in vitro, we attempt to predict the number of Ap-ICLs in the cell.

Oxidation Enhances Binding of Extrahelical 5-Methyl-Cytosines by Thymine DNA Glycosylase

The DNA repair protein thymine DNA glycosylase (TDG) removes mispaired or damaged bases, such as oxidized methyl-cytosine, from DNA by cleavage of the glycosidic bond between the sugar and the target base flipped into the enzyme's active site. The enzyme is active against formyl-cytosine and carboxyl-cytosine, whereas the lower oxidized hydroxymethyl-cytosine and methyl-cytosine itself are not processed by the enzyme. Molecular dynamics simulations with thermodynamic integration of TDG complexed to DNA carrying one of four different (oxidized) methyl-cytosine bases in extrahelcial conformation, methyl-cytosine (mC), hydroxymethyl-cytosine (hmC), formyl-cytosine (fC), or carboxyl-cytosine (caC), show a more favorable binding affinity of the higher oxidized forms, fC and caC, than the nonsubstrate bases hmC and mC. Despite rather comparable, reaction-competent conformations of the flipped bases in the active site of the enzyme, more and stronger interactions with active site residues account for the preferred binding of the higher oxidized bases. Binding of the negatively charged caC and the neutral fC are strengthened by interactions with positively charged His151. Our calculated proton affinities find this protonation state of His151 the preferred one in the presence of caC and conceivable in the presence of fC as well as increasing the binding affinity toward the two bases. Discrimination of the substrate bases is further achieved by the backbone of Tyr152 that forms a strong hydrogen bond to the carboxyl and formyl oxygen atoms of caC and fC, respectively, a contact that is completely lacking in mC and much weaker in hmC. Overall, our computational results indicate that the enzyme discriminates the different oxidation forms of methyl-cytosine already at the formation of the extrahelical complexes.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

Markov State Models to Study the Functional Dynamics of Proteins in the Wake of Machine Learning

Markov state models (MSMs) based on molecular dynamics (MD) simulations are routinely employed to study protein folding, however, their application to functional conformational changes of biomolecules is still limited. In the past few years, the field of computational chemistry has experienced a surge of advancements stemming from machine learning algorithms, and MSMs have not been left out. Unlike global processes, such as protein folding, the application of MSMs to functional conformational changes is challenging because they mostly consist of localized structural transitions. Therefore, it is critical to properly select a subset of structural features that can describe the slowest dynamics of these functional conformational changes. To address this challenge, we recommend several automatic feature selection methods such as Spectral-OASIS. To identify states in MSMs, the chosen features can be subject to dimensionality reduction methods such as TICA or deep learning based VAMPNets to project MD conformations onto a few collective variables for subsequent clustering. Another challenge for the application of MSMs to the study of functional conformational changes is the ability to comprehend their biophysical mechanisms, as MSMs built for these processes often require a large number of states. We recommend the recently developed quasi-MSMs (qMSMs) to address this issue. Compared to MSMs, qMSMs encode the non-Markovian dynamics via the generalized master equation and can significantly reduce the number of states. As a result, qMSMs can be built with a handful of states to facilitate the interpretation of functional conformational changes. In the wake of machine learning, we believe that the rapid advancement in the MSM methodology will lead to their wider application in studying functional conformational changes of biomolecules.

Interaction of Thymine DNA Glycosylase with Oxidised 5-Methyl-cytosines in Their Amino- and Imino-Forms

Thymine DNA Glycosylase (TDG) is an enzyme of the base excision repair mechanism and removes damaged or mispaired bases from DNA via hydrolysis of the glycosidic bond. Specificity is of high importance for such a glycosylase, so as to avoid the damage of intact DNA. Among the substrates reported for TDG are mispaired uracil and thymine but also formyl-cytosine and carboxyl-cytosine. Methyl-cytosine and hydroxylmethyl-cytosine are, in contrast, not processed by the TDG enzyme. We have in this work employed molecular dynamics simulations to explore the conformational dynamics of DNA carrying a formyl-cytosine or carboxyl-cytosine and compared those to DNA with the non-cognate bases methyl-cytosine and hydroxylmethyl-cytosine, as amino and imino tautomers. Whereas for the mispairs a wobble conformation is likely decisive for recognition, all amino tautomers of formyl-cytosine and carboxyl-cytosine exhibit the same Watson–Crick conformation, but all imino tautomers indeed form wobble pairs. The conformational dynamics of the amino tautomers in free DNA do not exhibit differences that could be exploited for recognition, and also complexation to the TDG enzyme does not induce any alteration that would indicate preferable binding to one or the other oxidised methyl-cytosine. The imino tautomers, in contrast, undergo a shift in the equilibrium between a closed and a more open, partially flipped state, towards the more open form upon complexation to the TDG enzyme. This stabilisation of the more open conformation is most pronounced for the non-cognate bases methyl-cytosine and hydroxyl-cytosine and is thus not a likely mode for recognition. Moreover, calculated binding affinities for the different forms indicate the imino forms to be less likely in the complexed DNA. These findings, together with the low probability of imino tautomers in free DNA and the indifference of the complexed amino tautomers, suggest that discrimination of the oxidised methyl-cytosines does not take place in the initial complex formation.

The Research of G-Motif Construction and Chirality in Deoxyguanosine Monophosphate Nucleotide Complexes

A detailed understanding of the mismatched base-pairing interactions in DNA will help reveal genetic diseases and provide a theoretical basis for the development of targeted drugs. Here, we utilized mononucleotide fragment to simulate mismatch DNA interactions in a local hydrophobic microenvironment. The bipyridyl-type bridging ligands were employed as a mild stabilizer to stabilize the GG mismatch containing complexes, allowing mismatch to be visualized based on X-ray crystallography. Five single crystals of 2′-deoxyguanosine–5′–monophosphate (dGMP) metal complexes were designed and obtained via the process of self-assembly. Crystallographic studies clearly reveal the details of the supramolecular interaction between mononucleotides and guest intercalators. A novel guanine–guanine base mismatch pattern with unusual (high anti)–(high anti) type of arrangement around the glycosidic angle conformations was successfully constructed. The solution state ¹H–NMR, ESI–MS spectrum studies, and UV titration experiments emphasize the robustness of this g–motif in solution. Additionally, we combined the methods of single-crystal and solution-, solid-state CD spectrum together to discuss the chirality of the complexes. The complexes containing the g–motif structure, which reduces the energy of the system, following the solid-state CD signals, generally move in the long-wave direction. These results provided a new mismatched base pairing, that is g–motif. The interaction mode and full characterizations of g–motif will contribute to the study of the mismatched DNA interaction.

Hg(II) Binding to Thymine Bases in DNA

The compounds of mercury can be highly toxic and can interfere with a range of biological processes, although many aspects of the mechanism of toxicity are still obscure or unknown. One especially intriguing property of Hg(II) is its ability to bind DNA directly, making interstrand cross-links between thymine nucleobases in AT-rich sequences. We have used a combination of small molecule X-ray diffraction, X-ray spectroscopies, and computational chemistry to study the interactions of Hg(II) with thymine. We find that the energetically preferred mode of thymine binding in DNA is to the N3 and predict only minor distortions of the DNA structure on binding one Hg(II) to two cross-adjacent thymine nucleotides. The preferred geometry is predicted to be twisted away from coplanar through a torsion angle of between 32 and 43°. Using 1-methylthymine as a model, the bis-thymine coordination of Hg(II) is found to give a highly characteristic X-ray spectroscopic signature that is quite distinct from other previously described biological modes of binding of Hg(II). This work enlarges and deepens our view of significant biological targets of Hg(II) and demonstrates tools that can provide a characteristic signature for the binding of Hg(II) to DNA in more complex matrices including intact cells and tissues, laying the foundation for future studies of mechanisms of mercury toxicity.

Assessing the sequence dependence of pyrimidine-pyrimidone (6-4) photoproduct in a duplex double-stranded DNA: A pitfall for microsecond range simulation

Sequence dependence of the (6-4) photoproduct conformational landscape when embedded in six 25-bp duplexes is evaluated along extensive unbiased and enhanced (replica exchange with solute tempering, REST2) molecular dynamics simulations. The structural reorganization as the central pyrimidines become covalently tethered is traced back in terms of non-covalent interactions, DNA bending, and extrusion of adenines of the opposite strands. The close sequence pattern impacts the conformational landscape around the lesion, inducing different upstream and downstream flexibilities. Moreover, REST2 simulations allow us to probe structures possibly important for damaged DNA recognition.

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.

Reaction Coordinate and Thermodynamics of Base Flipping in RNA

Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.

Mapping the DNA-Binding Motif of Scabin Toxin, a Guanine Modifying Enzyme from Streptomyces scabies

Scabin is a mono-ADP-ribosyltransferase toxin/enzyme and possible virulence factor produced by the agriculture pathogen, Streptomyces scabies. Recently, molecular dynamic approaches and MD simulations revealed its interaction with both NAD⁺ and DNA substrates. An Essential Dynamics Analysis identified a crab-claw-like mechanism, including coupled changes in the exposed motifs, and the R_β1-R_La-N_Lc-STT_β2-W_PN-W_ARTT-(QxE)_ARTT sequence motif was proposed as a catalytic signature of the Pierisin family of DNA-acting toxins. A new fluorescence assay was devised to measure the kinetics for both RNA and DNA substrates. Several protein variants were prepared to probe the Scabin-NAD-DNA molecular model and to reveal the reaction mechanism for the transfer of ADP-ribose to the guanine base in the DNA substrate. The results revealed that there are several lysine and arginine residues in Scabin that are important for binding the DNA substrate; also, key residues such as Asn110 in the mechanism of ADP-ribose transfer to the guanine base were identified. The DNA-binding residues are shared with ScARP from Streptomyces coelicolor but are not conserved with Pierisin-1, suggesting that the modification of guanine bases by ADP-ribosyltransferases is divergent even in the Pierisin family.

Elucidating molecular mechanisms of functional conformational changes of proteins via Markov state models

Functional conformational changes of proteins can facilitate numerous biological events in cells. The Markov state model (MSM) built from molecular dynamics simulations provide a powerful approach to study them. We here introduce a protocol that is tailor-made for constructing MSMs to study the functional conformational changes of proteins. In this protocol, one of the important steps is to select proper molecular features that can collectively describe the slowest timescales of conformational changes of interest. We recommend spectral oASIS, the modified version of oASIS, as a promising approach for automatic feature selection. Recently developed deep learning methods could also serve efficient approaches for selecting features and finding collective variables. Using DNA repair enzymes and RNA polymerases as examples, we review recent applications of MSMs to elucidate molecular mechanisms of functional conformational changes. Finally, we discuss remaining challenges and future perspectives for constructing MSMs to study functional conformational changes of proteins.

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.

Opening dynamics of HIV-1 gp120 upon receptor binding is dictated by a key hydrophobic core

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.

Polymorphic G:G mismatches act as hotspots for inducing right-handed Z DNA by DNA intercalation

DNA mismatches are highly polymorphic and dynamic in nature, albeit poorly characterized structurally. We utilized the antitumour antibiotic Co^II(Chro)₂ (Chro = chromomycin A3) to stabilize the palindromic duplex d(TTGGCGAA) DNA with two G:G mismatches, allowing X-ray crystallography-based monitoring of mismatch polymorphism. For the first time, the unusual geometry of several G:G mismatches including syn–syn, water mediated anti–syn and syn–syn-like conformations can be simultaneously observed in the crystal structure. The G:G mismatch sites of the d(TTGGCGAA) duplex can also act as a hotspot for the formation of alternative DNA structures with a GC/GA-5′ intercalation site for binding by the GC-selective intercalator actinomycin D (ActiD). Direct intercalation of two ActiD molecules to G:G mismatch sites causes DNA rearrangements, resulting in backbone distortion to form right-handed Z-DNA structures with a single-step sharp kink. Our study provides insights on intercalators-mismatch DNA interactions and a rationale for mismatch interrogation and detection via DNA intercalation.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Effect of the thymine-DNA glycosylase rs4135050 variant on Saudi smoker population

Background

Thymine‐DNA glycosylase (TDG) is an essential DNA‐repair enzyme which works in both epigenetic regulation and genome maintenance. It is also responsible for efficient correction of multiple endogenous DNA lesions which occur commonly in mammalian genomes. Research of genetic variants such as SNPs, resulting in disease, is predicted to yield clinical advancements through the identification of sensitive genetic markers and the development of disease prevention and therapy. To that end, the main objective of the present study is to identify the possible interactions between cigarette smoking and the rs4135050 variant of the TDG gene, situated in the intron position, among Saudi individuals.

Methods

TDG rs4135050 (A/T) was investigated by genotyping 239, and 235 blood specimens were obtained from nonsmokers and smokers of cigarette respectively.

Results

T allele frequency was found which showed a significant protective effect on Saudi male smokers (OR = 0.64, p = 0.0187) compared to nonsmoking subjects, but not in female smokers. Furthermore, smokers aged less than 29 years, the AT and AT+TT genotypes decreased more than four times the risk of initiation of smoking related‐diseases compare to the ancestral AA homozygous genotype. Paradoxically, the AT (OR = 3.88, p = 0.0169) and AT+TT (OR = 2.86, p = 0.0420) genotypes were present at a higher frequency in smoking patients aged more than 29 years as compared to nonsmokers at the same ages.

Conclusion

Depending on the gender and age of patients, TDG rs4135050 may provide a novel biomarker for the early diagnosis and prevention of several diseases caused by cigarette smoking.

Interactions of the DNA Repair Enzyme Human Thymine DNA Glycosylase with Cognate and Noncognate DNA

Glycosylases specifically recognize and flip their target base out of the DNA helix into the enzyme's active site. Our simulations show that a partially flipped state, already present in free DNA carrying a T:G mispair, becomes the more probable state compared to the closed state after binding of thymine DNA glycosylase (TDG). Paired thymine (T:A) or methyl-cytosine (mC:G) does not exhibit a partially flipped state in free or complexed DNA. Important enzyme-DNA interactions exhibit significant strength in the intrahelical and extrahelical TDG-DNA complexes. The computed binding free energy differences suggest these interactions account for the stabilization of the partially flipped state, thereby driving the T:G mispair toward base flip. In the fully flipped state, the cognate base thymine is significantly better accommodated in the enzyme's active site than noncognate bases are, suggesting the hydrolysis step as the last of several stages at which base recognition can be achieved.

Nucleosomes and the three glycosylases: High, medium, and low levels of excision by the uracil DNA glycosylase superfamily

Human cells express the UDG superfamily of glycosylases, which excise uracil (U) from the genome. The three members of this structural superfamily are uracil DNA glycosylase (UNG/UDG), single-strand selective monofunctional uracil DNA glycosylase (SMUG1), and thymine DNA glycosylase (TDG). We previously reported that UDG is efficient at removing U from DNA packaged into nucleosome core particles (NCP) and is minimally affected by the histone proteins when acting on an outward-facing U in the dyad region. In an effort to determine whether this high activity is a general property of the UDG superfamily of glycosylases, we compare the activity of UDG, SMUG1, and TDG on a U:G wobble base pair using NCP assembled from Xenopus laevis histones and the Widom 601 positioning sequence. We found that while UDG is highly active, SMUG1 is severely inhibited on NCP and this inhibition is independent of sequence context. Here we also provide the first report of TDG activity on an NCP, and found that TDG has an intermediate level of activity in excision of U and is severely inhibited in its excision of T. These results are discussed in the context of cellular roles for each of these enzymes.