Facilitated Diffusion Mechanisms in DNA Base Excision Repair and Transcriptional Activation

Quantitative morphological analysis of Deinococcus radiodurans elucidates complex dose-dependent nucleoid condensation during recovery from ionizing radiation

The extremophile Deinococcus radiodurans maintains a highly organized and condensed nucleoid as its default state, possibly contributing to its high tolerance to ionizing radiation (IR). Previous studies of the D. radiodurans nucleoid were limited by reliance on manual image annotation and qualitative metrics. Here, we introduce a high-throughput approach to quantify the geometric properties of cells and nucleoids using confocal microscopy, digital reconstructions of cells, and computational modeling. We utilize this novel approach to investigate the dynamic process of nucleoid condensation in response to IR stress. Our quantitative analysis reveals that at the population level, exposure to IR induced nucleoid compaction and decreased the size of D. radiodurans cells. Morphological analysis and clustering identified six distinct sub-populations across all tested experimental conditions. Results indicate that exposure to IR induced fractional redistributions of cells across sub-populations to exhibit morphologies associated with greater nucleoid condensation and decreased the abundance of sub-populations associated with cell division. Nucleoid-associated proteins (NAPs) may link nucleoid compaction and stress tolerance, but their roles in regulating compaction in D. radiodurans are unknown. Imaging of genomic mutants of known and suspected NAPs that contribute to nucleoid condensation found that deletion of nucleic acid-binding proteins, not previously described as NAPs, can remodel the nucleoid by driving condensation or decondensation in the absence of stress and that IR increased the abundance of these morphological states. Thus, our integrated analysis introduces a new methodology for studying environmental influences on bacterial nucleoids and provides an opportunity to further investigate potential regulators of nucleoid condensation.IMPORTANCEDeinococcus radiodurans, an extremophile known for its stress tolerance, constitutively maintains a highly condensed nucleoid. Qualitative studies have described nucleoid behavior under a variety of conditions. However, a lack of quantitative data regarding nucleoid organization and dynamics has limited our understanding of the regulatory mechanisms controlling nucleoid organization in D. radiodurans. Here, we introduce a quantitative approach that enables high-throughput quantitative measurements of subcellular spatial characteristics in bacterial cells. Applying this to wild-type or single-protein-deficient populations of D. radiodurans subjected to ionizing radiation, we identified significant stress-responsive changes in cell shape, nucleoid organization, and morphology. These findings highlight this methodology's adaptability and capacity for quantitatively analyzing the cellular response to stressors for screening cellular proteins involved in bacterial nucleoid organization.

Dynamic Structures and Fast Transition Kinetics of Oxidized G-Quadruplexes

8-oxoguanines (8-oxoG) in cells form compromised G-quadruplexes (GQs), which may vary GQ mediated gene regulations. By mimicking molecularly crowded cellular environment using 40% DMSO or sucrose, here it is found that oxidized human telomeric GQs have stabilities close to the wild-type (WT) GQs. Surprisingly, while WT GQs show negative formation cooperativity between a Pt(II) binder and molecularly crowded environment, positive cooperativity is observed for oxidized GQ formation. Single-molecule mechanical unfolding reveals that 8-oxoG sequence formed more diverse and flexible structures with faster folding/unfolding transition kinetics, which facilitates the Pt(II) ligand to bind the best-fit structures with positive cooperativity. These findings offer new understanding on structures and properties of oxidized G-rich species in crowded environments. They also provide insights into the design of better ligands to target oxidized G-rich structures formed under oxidative cell stress.

Three-State Diffusion Model of DNA Glycosylase Translocation along Stretched DNA as Revealed by Free Energy Landscapes at the All-Atom Level

DNA glycosylases play key roles in the maintenance of genomic integrity. These enzymes effectively find rare damaged sites in DNA and participate in subsequent base excision repair. Single-molecule and ensemble experiments have revealed key aspects of this damage-site searching mechanism and the involvement of facilitated diffusion. In this study, we describe free energy landscapes of enzyme translocation along nonspecific DNA obtained using a fully atomistic molecular dynamics (MD) simulation of a well-known DNA glycosylase, human 8-oxoguanine DNA glycosylase 1 (hOGG1). Based on an analysis of simulated free energy profiles, we propose a three-state model for the damage-site searching mechanism. In the three states, named the L1, L2, and L3 states, the L1 state is a helical sliding mode in close contact with DNA, whereas the L2 state is a major- or minor-groove tracking mode in loose contact with DNA and the L3 state is a two-dimensional freely diffusing mode during which hOGG1 is somewhat removed from the DNA surface (∼24 Å away from the surface). This three-state model well describes key experimental findings obtained from single-molecule and ensemble experiments and provides a unified molecular picture of the DNA lesion-searching mechanism of hOGG1.

Analyzing paramagnetic NMR data on target DNA search by proteins using a discrete-state kinetic model for translocation

Before reaching their targets, sequence-specific DNA-binding proteins nonspecifically bind to DNA through electrostatic interactions and stochastically change their locations on DNA. Investigations into the dynamics of DNA-scanning by proteins are nontrivial due to the simultaneous presence of multiple translocation mechanisms and many sites for the protein to nonspecifically bind to DNA. Nuclear magnetic resonance (NMR) spectroscopy can provide information about the target DNA search processes at an atomic level. Paramagnetic relaxation enhancement (PRE) is particularly useful to study how the proteins scan DNA in the search process. Previously, relatively simple two-state or three-state exchange models were used to explain PRE data reflecting the target search process. In this work, using more realistic discrete-state stochastic kinetics models embedded into an NMR master equation, we analyzed the PRE data for the HoxD9 homeodomain interacting with DNA. The kinetic models that incorporate sliding, dissociation, association, and intersegment transfer can reproduce the PRE profiles observed at some different ionic strengths. The analysis confirms the previous interpretation of the PRE data and shows that the protein's probability distribution among nonspecific sites is nonuniform during the target DNA search process.

Chromatin organization drives the search mechanism of nuclear factors

Nuclear factors rapidly scan the genome for their targets, but the role of nuclear organization in such search is uncharted. Here we analyzed how multiple factors explore chromatin, combining live-cell single-molecule tracking with multifocal structured illumination of DNA density. We find that factors displaying higher bound fractions sample DNA-dense regions more exhaustively. Focusing on the tumor-suppressor p53, we demonstrate that it searches for targets by alternating between rapid diffusion in the interchromatin compartment and compact sampling of chromatin dense regions. Efficient targeting requires balanced interactions with chromatin: fusing p53 with an exogenous intrinsically disordered region potentiates p53-mediated target gene activation at low concentrations, but leads to condensates at higher levels, derailing its search and downregulating transcription. Our findings highlight the role of disordered regions on factors search and showcase a powerful method to generate traffic maps of the eukaryotic nucleus to dissect how its organization guides nuclear factors action.

Backbone Conformational Equilibrium in Mismatched DNA Correlates with Enzyme Activity

T:G mismatches in mammals arise primarily from the deamination of methylated CpG sites or the incorporation of improper nucleotides. The process by which repair enzymes such as thymine DNA glycosylase (TDG) identify a canonical DNA base in the incorrect pairing context remains a mystery. However, the abundant contacts of the repair enzymes with the DNA backbone suggest a role for protein-phosphate interaction in the recognition and repair processes, where conformational properties may facilitate the proper interactions. We have previously used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion and the effect of a mismatch or lesion compared to canonical DNA and found stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA. We have currently compared our results to substrate dependence for TDG, MBD4, M. HhaI, and CEBPβ, testing for correlations to sequence and base-pair dependence. We found strong correlations of our DNA phosphate backbone equilibrium (Keq) to different enzyme kinetics or binding parameters of these varied enzymes, suggesting that the backbone equilibrium may play an important role in mismatch recognition and/or conformational rearrangement and energetics during nucleotide flipping or other aspects of enzyme interrogation of the DNA substrate.

Single-molecule tracking to determine the abundances and stoichiometries of freely-diffusing protein complexes in living cells: Past applications and future prospects

Most biological processes in living cells rely on interactions between proteins. Live-cell compatible approaches that can quantify to what extent a given protein participates in homo- and hetero-oligomeric complexes of different size and subunit composition are therefore critical to advance our understanding of how cellular physiology is governed by these molecular interactions. Biomolecular complex formation changes the diffusion coefficient of constituent proteins, and these changes can be measured using fluorescence microscopy-based approaches, such as single-molecule tracking, fluorescence correlation spectroscopy, and fluorescence recovery after photobleaching. In this review, we focus on the use of single-molecule tracking to identify, resolve, and quantify the presence of freely-diffusing proteins and protein complexes in living cells. We compare and contrast different data analysis methods that are currently employed in the field and discuss experimental designs that can aid the interpretation of the obtained results. Comparisons of diffusion rates for different proteins and protein complexes in intracellular aqueous environments reported in the recent literature reveal a clear and systematic deviation from the Stokes-Einstein diffusion theory. While a complete and quantitative theoretical explanation of why such deviations manifest is missing, the available data suggest the possibility of weighing freely-diffusing proteins and protein complexes in living cells by measuring their diffusion coefficients. Mapping individual diffusive states to protein complexes of defined molecular weight, subunit stoichiometry, and structure promises to provide key new insights into how protein-protein interactions regulate protein conformational, translational, and rotational dynamics, and ultimately protein function.

Dynamic interplay between target search and recognition for a Type I CRISPR-Cas system

CRISPR-Cas effector complexes enable the defense against foreign nucleic acids and have recently been exploited as molecular tools for precise genome editing at a target locus. To bind and cleave their target, the CRISPR-Cas effectors have to interrogate the entire genome for the presence of a matching sequence. Here we dissect the target search and recognition process of the Type I CRISPR-Cas complex Cascade by simultaneously monitoring DNA binding and R-loop formation by the complex. We directly quantify the effect of DNA supercoiling on the target recognition probability and demonstrate that Cascade uses facilitated diffusion for its target search. We show that target search and target recognition are tightly linked and that DNA supercoiling and limited 1D diffusion need to be considered when understanding target recognition and target search by CRISPR-Cas enzymes and engineering more efficient and precise variants.

Targeting the p53 signaling pathway in cancers: Molecular mechanisms and clinical studies

Tumor suppressor p53 can transcriptionally activate downstream genes in response to stress, and then regulate the cell cycle, DNA repair, metabolism, angiogenesis, apoptosis, and other biological responses. p53 has seven functional domains and 12 splice isoforms, and different domains and subtypes play different roles. The activation and inactivation of p53 are finely regulated and are associated with phosphorylation/acetylation modification and ubiquitination modification, respectively. Abnormal activation of p53 is closely related to the occurrence and development of cancer. While targeted therapy of the p53 signaling pathway is still in its early stages and only a few drugs or treatments have entered clinical trials, the development of new drugs and ongoing clinical trials are expected to lead to the widespread use of p53 signaling-targeted therapy in cancer treatment in the future. TRIAP1 is a novel p53 downstream inhibitor of apoptosis. TRIAP1 is the homolog of yeast mitochondrial intermembrane protein MDM35, which can play a tumor-promoting role by blocking the mitochondria-dependent apoptosis pathway. This work provides a systematic overview of recent basic research and clinical progress in the p53 signaling pathway and proposes that TRIAP1 is an important therapeutic target downstream of p53 signaling.

Correlated Target Search by Vaccinia Virus Uracil-DNA Glycosylase, a DNA Repair Enzyme and a Processivity Factor of Viral Replication Machinery

The protein encoded by the vaccinia virus D4R gene has base excision repair uracil-DNA N-glycosylase (vvUNG) activity and also acts as a processivity factor in the viral replication complex. The use of a protein unlike PolN/PCNA sliding clamps is a unique feature of orthopoxviral replication, providing an attractive target for drug design. However, the intrinsic processivity of vvUNG has never been estimated, leaving open the question whether it is sufficient to impart processivity to the viral polymerase. Here, we use the correlated cleavage assay to characterize the translocation of vvUNG along DNA between two uracil residues. The salt dependence of the correlated cleavage, together with the similar affinity of vvUNG for damaged and undamaged DNA, support the one-dimensional diffusion mechanism of lesion search. Unlike short gaps, covalent adducts partly block vvUNG translocation. Kinetic experiments show that once a lesion is found it is excised with a probability ~0.76. Varying the distance between two uracils, we use a random walk model to estimate the mean number of steps per association with DNA at ~4200, which is consistent with vvUNG playing a role as a processivity factor. Finally, we show that inhibitors carrying a tetrahydro-2,4,6-trioxopyrimidinylidene moiety can suppress the processivity of vvUNG.

Theoretical insights into the full description of DNA target search by subdiffusing proteins

DNA binding proteins (DBPs) diffuse in the cytoplasm to recognise and bind with their respective target sites on the DNA to initiate several biologically important processes. The first passage time distributions (FPTDs) of DBPs are useful in quantifying the timescales of the most-probable search paths in addition to the mean value of the distribution which, strikingly, are decades of order apart in time. However, extremely crowded in vivo conditions or the viscoelasticity of the cellular medium among other factors causes biomolecules to exhibit anomalous diffusion which is usually overlooked in most theoretical studies. We have obtained approximate analytical expressions of a general FPTD and the two characteristic timescales that are valid for any single subdiffusing protein searching for its target in vivo. Our results can be applied to single-particle tracking experiments of target search.

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.

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.

Photoreaction of DNA Containing 5-Halouracil and its Products

5-Halouracil, which is a DNA base analog in which the methyl group at the C5 position of thymine is replaced with a halogen atom, has been used in studies of DNA damage. In DNA strands, the uracil radical generated from 5-halouracil causes DNA damage via a hydrogen-abstraction reaction. We analyzed the photoreaction of 5-halouracil in various DNA structures and revealed that the reaction is DNA structure-dependent. In this review, we summarize the results of the analysis of the reactivity of 5-halouracil in various DNA local structures. Among the 5-halouracil molecules, 5-bromouracil has been used as a probe in the analysis of photoinduced electron transfer through DNA. The analysis of groove-binder/DNA and protein/DNA complexes using a 5-bromouracil-based electron transfer system is also described.

Influence of Nonspecific Interactions between Proteins and In Vivo Cytoplasmic Crowders in Facilitated Diffusion of Proteins: Theoretical Insights

The binding of proteins to their respective specific sites on the DNA through facilitated diffusion serves as the initial step of various important biological processes. While this search process has been thoroughly investigated via in vitro studies, the cellular environment is complex and may interfere with the protein's search dynamics. The cytosol is heavily crowded, which can potentially modify the search by nonspecifically interacting with the protein that has been mostly overlooked. In this work, we probe the target search dynamics in the presence of explicit crowding agents that have an affinity toward the protein. We theoretically investigate the role of such protein-crowder associations in the target search process using a discrete-state stochastic framework that allows for the analytical description of dynamic properties. It is found that stronger nonspecific associations between the crowder and proteins can accelerate the facilitated diffusion of proteins in comparison with a purely inert, rather weakly interacting cellular environment. This effect depends on how strong these associations are, the spatial positions of the target with respect to the crowders, and the size of the crowded region. Our theoretical results are also tested with Monte Carlo computer simulations. Our predictions are in qualitative agreement with existing experimental observations and computational studies.

Cognate DNA Recognition by Engrailed Homeodomain Involves a Conformational Change Controlled via an Electrostatic-Spring-Loaded Latch

Transcription factors must scan genomic DNA, recognize the cognate sequence of their control element(s), and bind tightly to them. The DNA recognition process is primarily carried out by their DNA binding domains (DBD), which interact with the cognate site with high affinity and more weakly with any other DNA sequence. DBDs are generally thought to bind to their cognate DNA without changing conformation (lock-and-key). Here, we used nuclear magnetic resonance and circular dichroism to investigate the interplay between DNA recognition and DBD conformation in the engrailed homeodomain (enHD), as a model case for the homeodomain family of eukaryotic DBDs. We found that the conformational ensemble of enHD is rather flexible and becomes gradually more disordered as ionic strength decreases following a Debye–Hückel’s dependence. Our analysis indicates that enHD’s response to ionic strength is mediated by a built-in electrostatic spring-loaded latch that operates as a conformational transducer. We also found that, at moderate ionic strengths, enHD changes conformation upon binding to cognate DNA. This change is of larger amplitude and somewhat orthogonal to the response to ionic strength. As a consequence, very high ionic strengths (e.g., 700 mM) block the electrostatic-spring-loaded latch and binding to cognate DNA becomes lock-and-key. However, the interplay between enHD conformation and cognate DNA binding is robust across a range of ionic strengths (i.e., 45 to 300 mM) that covers the physiologically-relevant conditions. Therefore, our results demonstrate the presence of a mechanism for the conformational control of cognate DNA recognition on a eukaryotic DBD. This mechanism can function as a signal transducer that locks the DBD in place upon encountering the cognate site during active DNA scanning. The electrostatic-spring-loaded latch of enHD can also enable the fine control of DNA recognition in response to transient changes in local ionic strength induced by variate physiological processes.

Pseudouridine synthase 7 is an opportunistic enzyme that binds and modifies substrates with diverse sequences and structures

Pseudouridine (Ψ) is a ubiquitous RNA modification incorporated by pseudouridine synthase (Pus) enzymes into hundreds of noncoding and protein-coding RNA substrates. Here, we determined the contributions of substrate structure and protein sequence to binding and catalysis by pseudouridine synthase 7 (Pus7), one of the principal messenger RNA (mRNA) modifying enzymes. Pus7 is distinct among the eukaryotic Pus proteins because it modifies a wider variety of substrates and shares limited homology with other Pus family members. We solved the crystal structure of Saccharomyces cerevisiae Pus7, detailing the architecture of the eukaryotic-specific insertions thought to be responsible for the expanded substrate scope of Pus7. Additionally, we identified an insertion domain in the protein that fine-tunes Pus7 activity both in vitro and in cells. These data demonstrate that Pus7 preferentially binds substrates possessing the previously identified UGUAR (R = purine) consensus sequence and that RNA secondary structure is not a strong requirement for Pus7-binding. In contrast, the rate constants and extent of Ψ incorporation are more influenced by RNA structure, with Pus7 modifying UGUAR sequences in less-structured contexts more efficiently both in vitro and in cells. Although less-structured substrates were preferred, Pus7 fully modified every transfer RNA, mRNA, and nonnatural RNA containing the consensus recognition sequence that we tested. Our findings suggest that Pus7 is a promiscuous enzyme and lead us to propose that factors beyond inherent enzyme properties (e.g., enzyme localization, RNA structure, and competition with other RNA-binding proteins) largely dictate Pus7 substrate selection.

Complete Mapping of DNA‐Protein Interactions at the Single‐Molecule Level

DNA–protein interaction plays an essential role in the storage, expression, and regulation of genetic information. A 1D/3D facilitated diffusion mechanism has been proposed to explain the extraordinarily rapid rate of DNA‐binding protein (DBP) searching for cognate sequence along DNA and further studied by single‐molecule experiments. However, direct observation of the detailed chronological protein searching image is still a formidable challenge. Here, for the first time, a single‐molecule electrical monitoring technique is utilized to realize label‐free detection of the DBP–DNA interaction process based on high‐gain silicon nanowire field‐effect transistors (SiNW FETs). The whole binding process of WRKY domain and DNA has been visualized with high sensitivity and single‐base resolution. Impressively, the swinging of hydrogen bonds between amino acid residues and bases in DNA induce the dynamic collective motion of DBP–DNA. This in situ, label‐free electrical detection platform provides a practical experimental methodology for dynamic studies of various biomolecules.

A DNAzyme-driven random biped DNA walking nanomachine for sensitive detection of uracil-DNA glycosylase activity

Highly specific and ultrasensitive detection of uracil-DNA glycosylase (UDG) activity is of great significance for maintaining genomic integrity and medical research of related diseases. Here, we constructed a random DNA walking nanomachine based on a DNAzyme for UDG activity detection on the AuNP (Au nanoparticle) surface. When UDG is present, the U bases in the Y structure are removed, resulting in AP sites, which will be cleaved by Endo-IV to generate a 3' concave end for Exo-III, causing the locking strand of the DNAzyme to be completely hydrolyzed by the Exo-III and release the walking strand to randomly pair with the substrate strand on the AuNP surface; then, the walking strand exerts its cleavage activity with the assistance of Mg2+ to cleave the substrate strand and keep the fluorophore 6-carboxyfluorescein (FAM) away from the surface of the AuNP, which restores the fluorescence signal of this system. In this way, sensitive detection of UDG can be realized, and the detection limit is as low as 3.69 × 10-6 U mL-1. In addition, we found that this method is highly specific to UDG and can be used to detect UDG specifically in complex samples, which has certain application prospects in biomedical research and clinical diagnosis related to UDG.

The needle and the haystack: single molecule tracking to probe the transcription factor search in eukaryotes

Transcription factors (TFs) regulate transcription of their target genes by identifying and binding to regulatory regions of the genome among billions of potential non-specific decoy sites, a task that is often presented as a 'needle in the haystack' challenge. The TF search process is now well understood in bacteria, but its characterization in eukaryotes needs to account for the complex organization of the nuclear environment. Here we review how live-cell single molecule tracking is starting to shed light on the TF search mechanism in the eukaryotic cell and we outline the future challenges to tackle in order to understand how nuclear organization modulates the TF search process in physiological and pathological conditions.

Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field

The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.

Recognition and repair of oxidatively generated DNA lesions in plasmid DNA by a facilitated diffusion mechanism

The oxidatively generated genotoxic spiroiminodihydantoin (Sp) lesions are well-known substrates of the base excision repair (BER) pathway initiated by the bifunctional DNA glycosylase NEIL1. In this work, we reported that the excision kinetics of the single Sp lesions site-specifically embedded in the covalently closed circular DNA plasmids (contour length 2686 base pairs) by NEIL1 are biphasic under single-turnover conditions ([NEIL1] ≫ [SpDNApl]) in contrast with monophasic excision kinetics of the same lesions embedded in147-mer Sp-modified DNA duplexes. Under conditions of a large excess of plasmid DNA base pairs over NEIL1 molecules, the kinetics of excision of Sp lesions are biphasic in nature, exhibiting an initial burst phase, followed by a slower rate of formation of excision products The burst phase is associated with NEIL1-DNA plasmid complexes, while the slow kinetic phase is attributed to the dissociation of non-specific NEIL1-DNA complexes. The amplitude of the burst phase is limited because of the competing non-specific binding of NEIL1 to unmodified DNA sequences flanking the lesion. A numerical analysis of the incision kinetics yielded a value of φ ≍ 0.03 for the fraction of NEIL1 encounters with plasmid molecules that result in the excision of the Sp lesion, and a characteristic dissociation time of non-specific NEIL1-DNA complexes (τ-ns ≍ 8 s). The estimated average DNA translocation distance of NEIL1 is ∼80 base pairs. This estimate suggests that facilitated diffusion enhances the probability that NEIL1 can locate its substrate embedded in an excess of unmodified plasmid DNA nucleotides by a factor of ∼10.

How to assess the structural dynamics of transcription factors by integrating sparse NMR and EPR constraints with molecular dynamics simulations

We review recent advances in modeling structural ensembles of transcription factors from nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopic data, integrated with molecular dynamics (MD) simulations. We focus on approaches that confirm computed conformational ensembles by sparse constraints obtained from magnetic resonance. This combination enables the deduction of functional and structural protein models even if nuclear Overhauser effects (NOEs) are too scarce for conventional structure determination. We highlight recent insights into the folding-upon-DNA binding transitions of intrinsically disordered transcription factors that could be assessed using such integrative approaches.

Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA

DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N⁶-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.

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.

The N-terminal domain of uracil-DNA glycosylase: Roles for disordered regions

The presence of uracil in DNA calls for rapid removal facilitated by the uracil-DNA glycosylase superfamily of enzymes, which initiates the base excision repair (BER) pathway. In humans, uracil excision is accomplished primarily by the human uracil-DNA glycosylase (hUNG) enzymes. In addition to BER, hUNG enzymes play a key role in somatic hypermutation to generate antibody diversity. hUNG has several isoforms, with hUNG1 and hUNG2 being the two major isoforms. Both isoforms contain disordered N-terminal domains, which are responsible for a wide range of functions, with minimal direct impact on catalytic efficiency. Subcellular localization of hUNG enzymes is directed by differing N-terminal sequences, with hUNG1 dedicated to mitochondria and hUNG2 dedicated to the nucleus. An alternative isoform of hUNG1 has also been identified to localize to the nucleus in mouse and human cell models. Furthermore, hUNG2 has been observed at replication forks performing both pre- and post-replicative uracil excision to maintain genomic integrity. Replication protein A (RPA) and proliferating cell nuclear antigen (PCNA) are responsible for recruitment to replication forks via protein-protein interactions with the N-terminus of hUNG2. These interactions, along with protein degradation, are regulated by various post-translational modifications within the N-terminal tail, which are primarily cell-cycle dependent. Finally, translocation on DNA is also mediated by interactions between the N-terminus and DNA, which is enhanced under molecular crowding conditions by preventing diffusion events and compacting tail residues. This review summarizes recent research supporting the emerging roles of the N-terminal domain of hUNG.

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

Lost in the Crowd: How Does Human 8-Oxoguanine DNA Glycosylase 1 (OGG1) Find 8-Oxoguanine in the Genome?

The most frequent DNA lesion resulting from an oxidative stress is 7,8-dihydro-8-oxoguanine (8-oxoG). 8-oxoG is a premutagenic base modification due to its capacity to pair with adenine. Thus, the repair of 8-oxoG is critical for the preservation of the genetic information. Nowadays, 8-oxoG is also considered as an oxidative stress-sensor with a putative role in transcription regulation. In mammalian cells, the modified base is excised by the 8-oxoguanine DNA glycosylase (OGG1), initiating the base excision repair (BER) pathway. OGG1 confronts the massive challenge that is finding rare occurrences of 8-oxoG among a million-fold excess of normal guanines. Here, we review the current knowledge on the search and discrimination mechanisms employed by OGG1 to find its substrate in the genome. While there is considerable data from in vitro experiments, much less is known on how OGG1 is recruited to chromatin and scans the genome within the cellular nucleus. Based on what is known of the strategies used by proteins searching for rare genomic targets, we discuss the possible scenarios allowing the efficient detection of 8-oxoG by OGG1.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

Target search and recognition mechanisms of glycosylase AlkD revealed by scanning FRET-FCS and Markov state models

DNA glycosylase repairs DNA damage to maintain the genome integrity, and thus it is essential for the survival of all organisms. However, it remains a long-standing puzzle how glycosylase diffuses along the genomic DNA to locate the sparse and aberrant lesion sites efficiently and accurately in the genome containing numerous base pairs. Previously, only the high-speed–low-accuracy search mode has been characterized experimentally, while the low-speed–high-accuracy mode is undetectable. Here, we observed the low-speed mode of glycosylase AlkD translocating, and further dissected its molecular mechanisms. To achieve this, we developed an integrated platform by combining scanning FRET-FCS with Markov state model. We expect that this platform can be widely applied to investigate other glycosylases and DNA-binding proteins.

DNA glycosylase is responsible for repairing DNA damage to maintain the genome stability and integrity. However, how glycosylase can efficiently and accurately recognize DNA lesions across the enormous DNA genome remains elusive. It has been hypothesized that glycosylase translocates along the DNA by alternating between a fast but low-accuracy diffusion mode and a slow but high-accuracy mode when searching for DNA lesions. However, the slow mode has not been successfully characterized due to the limitation in the spatial and temporal resolutions of current experimental techniques. Using a newly developed scanning fluorescence resonance energy transfer (FRET)–fluorescence correlation spectroscopy (FCS) platform, we were able to observe both slow and fast modes of glycosylase AlkD translocating on double-stranded DNA (dsDNA), reaching the temporal resolution of microsecond and spatial resolution of subnanometer. The underlying molecular mechanism of the slow mode was further elucidated by Markov state model built from extensive all-atom molecular dynamics simulations. We found that in the slow mode, AlkD follows an asymmetric diffusion pathway, i.e., rotation followed by translation. Furthermore, the essential role of Y27 in AlkD diffusion dynamics was identified both experimentally and computationally. Our results provided mechanistic insights on how conformational dynamics of AlkD–dsDNA complex coordinate different diffusion modes to accomplish the search for DNA lesions with high efficiency and accuracy. We anticipate that the mechanism adopted by AlkD to search for DNA lesions could be a general one utilized by other glycosylases and DNA binding proteins.

Remarkable Enhancement of Nucleotide Excision Repair of a Bulky Guanine Lesion in a Covalently Closed Circular DNA Plasmid Relative to the Same Linearized Plasmid

The excision of DNA lesions by human nucleotide excision repair (NER) has been extensively studied in human cell extracts. Employing DNA duplexes with fewer than 200 bp containing a single bulky, benzo[a]pyrene-derived guanine lesion (B[a]P-dG), the NER yields are typically on the order of ∼5-10%, or less. Remarkably, the NER yield is enhanced by a factor of ∼6 when the B[a]P-dG lesion is embedded in a covalently closed circular pUC19NN plasmid (contour length of 2686 bp) rather than in the same plasmid linearized by a restriction enzyme with the B[a]P-dG adduct positioned at the 945th nucleotide counted from the 5'-end of the linearized DNA molecules. Furthermore, the NER yield in the circular pUC19NN plasmid is ∼9 times greater than in a short 147-mer DNA duplex with the B[a]P-dG adduct positioned in the middle. Although the NER factors responsible for these differences were not explicitly identified here, we hypothesize that the initial DNA damage sensor XPC-RAD23B is a likely candidate; it is known to search for DNA lesions by a constrained one-dimensional search mechanism [Cheon, N. Y., et al. (2019) Nucleic Acids Res. 47, 8337-8347], and our results are consistent with the notion that it dissociates more readily from the blunt ends than from the inner regions of linear DNA duplexes, thus accounting for the remarkable enhancement in NER yields associated with the single B[a]P-dG adduct embedded in covalently closed circular plasmids.

Displacement of Slow-Turnover DNA Glycosylases by Molecular Traffic on DNA

In the base excision repair pathway, the initiating enzymes, DNA glycosylases, remove damaged bases and form long-living complexes with the abasic DNA product, but can be displaced by AP endonucleases. However, many nuclear proteins can move along DNA, either actively (such as DNA or RNA polymerases) or by passive one-dimensional diffusion. In most cases, it is not clear whether this movement is disturbed by other bound proteins or how collisions with moving proteins affect the bound proteins, including DNA glycosylases. We have used a two-substrate system to study the displacement of human OGG1 and NEIL1 DNA glycosylases by DNA polymerases in both elongation and diffusion mode and by D4, a passively diffusing subunit of a viral DNA polymerase. The OGG1–DNA product complex was disrupted by DNA polymerase β (POLβ) in both elongation and diffusion mode, Klenow fragment (KF) in the elongation mode and by D4. NEIL1, which has a shorter half-life on DNA, was displaced more efficiently. Hence, both possibly specific interactions with POLβ and nonspecific collisions (KF, D4) can displace DNA glycosylases from DNA. The protein movement along DNA was blocked by very tightly bound Cas9 RNA-targeted nuclease, providing an upper limit on the efficiency of obstacle clearance.

Sensitive determination of formamidopyrimidine DNA glucosylase based on phosphate group-modulated multi-enzyme catalysis and fluorescent copper nanoclusters

In this work, a method for quantifying the activity of formamidopyrimidine DNA glucosylase (Fpg) was designed based on phosphate group (P)-modulated multi-enzyme catalysis and fluorescent copper nanoclusters (CuNCs). By eliminating 8-oxoguanine from double-stranded DNA, Fpg generates a nick with P at both 3' and 5' termini. Subsequently, part of the DNA is digested by 5'P-activated lambda exonuclease (λ Exo), and the generated 3'P disables exonuclease I (Exo I), resulting in the generation of single-stranded DNA containing poly(thymine) (poly(T)). Using poly(T) as templates, CuNCs were prepared to emit intense fluorescence as the readout of this method. However, in the absence of Fpg, the originally modified 5'P triggers the digestion of λ Exo. In this case, fluorescence emission is not obtained because CuNCs cannot be formed without DNA templates. Therefore, the catalysis of λ Exo and Exo I can be tuned by 5'P and 3'P, which can be further used to determine the activity of Fpg. The fluorescent Fpg biosensor works in a "signal-on" manner with the feature of "zero" background noise, and thus shows desirable analytical features and good performance. Besides, Fpg in serum samples and cell lysate could be accurately detected with the biosensor, indicating the great value of the proposed system in practical and clinical analysis.

Alkyltransferase-like protein clusters scan DNA rapidly over long distances and recruit NER to alkyl-DNA lesions

Alkylation of guanine bases in DNA is detrimental to cells due to its high mutagenic and cytotoxic potential and is repaired by the alkyltransferase AGT. Additionally, alkyltransferase-like proteins (ATLs), which are structurally similar to AGTs, have been identified in many organisms. While ATLs are per se catalytically inactive, strong evidence has suggested that ATLs target alkyl lesions to the nucleotide excision repair system (NER). Using a combination of single-molecule and ensemble approaches, we show here recruitment of UvrA, the initiating enzyme of prokaryotic NER, to an alkyl lesion by ATL. We further characterize lesion recognition by ATL and directly visualize DNA lesion search by highly motile ATL and ATL-UvrA complexes on DNA at the molecular level. Based on the high similarity of ATLs and the DNA-interacting domain of AGTs, our results provide important insight in the lesion search mechanism, not only by ATL but also by AGT, thus opening opportunities for controlling the action of AGT for therapeutic benefit during chemotherapy.

Eukaryotic Base Excision Repair: New Approaches Shine Light on Mechanism

Genomic DNA is susceptible to endogenous and environmental stresses that modify DNA structure and its coding potential. Correspondingly, cells have evolved intricate DNA repair systems to deter changes to their genetic material. Base excision DNA repair involves a number of enzymes and protein cofactors that hasten repair of damaged DNA bases. Recent advances have identified macromolecular complexes that assemble at the DNA lesion and mediate repair. The repair of base lesions generally requires five enzymatic activities: glycosylase, endonuclease, lyase, polymerase, and ligase. The protein cofactors and mechanisms for coordinating the sequential enzymatic steps of repair are being revealed through a range of experimental approaches. We discuss the enzymes and protein cofactors involved in eukaryotic base excision repair, emphasizing the challenge of integrating findings from multiple methodologies. The results provide an opportunity to assimilate biochemical findings with cell-based assays to uncover new insights into this deceptively complex repair pathway.

Macromolecular crowding induces compaction and DNA binding in the disordered N-terminal domain of hUNG2

Many human DNA repair proteins have disordered domains at their N- or C-termini with poorly defined biological functions. We recently reported that the partially structured N-terminal domain (NTD) of human uracil DNA glycosylase 2 (hUNG2), functions to enhance DNA translocation in crowded environments and also targets the enzyme to single-stranded/double-stranded DNA junctions. To understand the structural basis for these effects we now report high-resolution heteronuclear NMR studies of the isolated NTD in the presence and absence of an inert macromolecular crowding agent (PEG8K). Compared to dilute buffer, we find that crowding reduces the degrees of freedom for the structural ensemble, increases the order of a PCNA binding motif and dramatically promotes binding of the NTD for DNA through a conformational selection mechanism. These findings shed new light on the function of this disordered domain in the context of the crowded nuclear environment.

Understanding biochemical processes in the presence of sub-diffusive behavior of biomolecules in solution and living cells

In order to maintain cellular function, biomolecules like protein, DNA, and RNAs have to diffuse to the target spaces within the cell. Changes in the cytosolic microenvironment or in the nucleus during the fulfillment of these cellular processes affect their mobility, folding, and stability thereby impacting the transient or stable interactions with their adjacent neighbors in the organized and dynamic cellular interior. Using classical Brownian motion to elucidate the diffusion behavior of these biomolecules is hard considering their complex nature. The understanding of biomolecular diffusion inside cells still remains elusive due to the lack of a proper model that can be extrapolated to these cases. In this review, we have comprehensively addressed the progresses in this field, laying emphasis on the different aspects of anomalous diffusion in the different biochemical reactions in cell interior. These experiment-based models help to explain the diffusion behavior of biomolecules in the cytosolic and nuclear microenvironment. Moreover, since understanding of biochemical reactions within living cellular system is our main focus, we coupled the experimental observations with the concept of sub-diffusion from in vitro to in vivo condition. We believe that the pairing between the understanding of complex behavior and structure-function paradigm of biological molecules would take us forward by one step in order to solve the puzzle around diseases caused by cellular dysfunction.

Enhancement of RNA/Ligand Association Kinetics via an Electrostatic Anchor

The diverse biological processes mediated by RNA rest upon its recognition of various ligands, including small molecules and nucleic acids. Nevertheless, a recent literature survey suggests that RNA molecular recognition of these ligands is slow, with association rate constants orders of magnitude below the diffusional limit. Thus, we were prompted to consider strategies for increasing RNA association kinetics. Proteins can accelerate ligand association via electrostatic forces, and here, using the Tetrahymena group I ribozyme, we provide evidence that electrostatic forces can accelerate RNA/ligand association. This RNA enzyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. The G 2′- and 3′-OH groups interact with an active site metal ion, termed M_C, within E·S·G, and we perturbed each of these contacts via −NH₃⁺ substitution. New and prior data indicate that G(2′NH₃⁺) and G(3′NH₃⁺) bind as strongly as G, suggesting that the −NH₃⁺ substituents of these analogues avoid repulsive interactions with M_C and make alternative interactions. Unexpectedly, removal of the adjacent −OH via −H substitution to give G(2′H,3′NH₃⁺) and G(2′NH₃⁺_,3′H) enhanced binding, in stark contrast to the deleterious effect of these substitutions on G binding. Pulse–chase experiments indicate that the −NH₃⁺ moiety of G(2′H,3′NH₃⁺) increases the rate of G association. These results suggest that the positively charged −NH₃⁺ group can act as a molecular “anchor” to increase the residence time of the encounter complex and thereby enhance productive binding. Electrostatic anchors may provide a broadly applicable strategy for the development of fast binding RNA ligands and RNA-targeted therapeutics.

Direct Observation and Analysis of the Dynamics of the Photoresponsive Transcription Factor GAL4

Herein, the direct visualization of the dynamic interaction between a photoresponsive transcription factor fusion, GAL4-VVD, and DNA using high-speed atomic force microscopy (HS-AFM) is reported. A series of different GAL4-VVD movements, such as binding, sliding, stalling, and dissociation, was observed. Inter-strand jumping on two double-stranded (ds) DNAs was also observed. Detailed analysis using a long substrate DNA strand containing five GAL4-binding sites revealed that GAL4-VVD randomly moved on the dsDNA using sliding and hopping to rapidly find specific binding sites, and then stalled to the specific sites to form a stable complex formation. These results suggest the existence of different conformations of the protein to enable sliding and stalling. This single-molecule imaging system with nanoscale resolution provides an insight into the searching mechanism used by DNA-binding proteins.