How p53 Molecules Solve the Target DNA Search Problem: A Review

Skipping events impose repeated binding attempts: profound kinetic implications of protein-DNA conformational changes

The kinetics of protein-DNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. Protein-DNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRY-DNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific protein-DNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance protein-DNA recognition despite its intrinsically slow kinetics.

miR-460b-5p promotes proliferation and differentiation of chicken myoblasts and targets RBM19 gene

The meat production of broilers is crucial to economic benefits of broiler industries, while the slaughter performance of broilers is directly determined by skeletal muscle development. Hence, the broiler breeding for growth traits shows a great importance. As a kind of small noncoding RNA, microRNA (miRNA) can regulate the expression of multiple genes and perform a wide range of regulation in organisms. Currently, more and more studies have confirmed that miRNAs are closely associated with skeletal muscle development of chickens. Based on our previous miR-seq analysis (accession number: PRJNA668199), miR-460b-5p was screened as one of the key miRNAs probably involved in the growth regulation of chickens. However, the regulatory effect of miR-460b-5p on the development of chicken skeletal muscles is still unclear. Therefore, miR-460b-5p was further used for functional validation at the cellular level in this study. The expression pattern of miR-460b-5p was investigated in proliferation and differentiation stages of chicken primary myoblasts. It was showed that the expression level of miR-460b-5p gradually decreased from the proliferation stage (GM 50%) to the lowest at 24 h of differentiation. As differentiation proceeded, miR-460b-5p expression increased significantly, reaching the highest and stabilizing at 72 h and 96 h of differentiation. Through mRNA quantitative analysis of proliferation marker genes, CCK-8 and Edu assays, miR-460b-5p was found to significantly facilitate the transition of myoblasts from G1 to S phase and promote chicken myoblast proliferation. mRNA and protein quantitative analysis of differentiation marker genes, as well as the indirect immunofluorescence results of myotubes, revealed that miR-460b-5p significantly stimulated myotube development and promote chicken myoblast differentiation. In addition, the target relationship was validated for miR-460b-5p according to the dual-luciferase reporter assay and mRNA quantitative analysis, which indicates that miR-460b-5p was able to regulate RBM19 expression by specifically binding to the 3' UTR of RBM19. In summary, miR-460b-5p has positive regulatory effects on the proliferation and differentiation of chicken myoblasts, and RBM19 is a target gene of miR-460b-5p.

DNA Condensation by Peptide-Conjugated PAMAM Dendrimers. Influence of Peptide Charge

Nucleic acid delivery to cells is an important therapeutic strategy that requires the transport of nucleic acids to intracellular compartments and their protection from enzymatic degradation. This can be achieved through the complexation of the nucleic acids with polycations. Poly(amidoamine) (PAMAM) dendrimers and peptide-conjugated dendrimers have been investigated as delivery vectors. Inspired by these studies and the role of flexible peptide domains in protein-DNA interactions, we studied the impact of conjugating two peptides (tails) to generation 2 (G2) PAMAM dendrimers on DNA condensation and polyplex formation. Using gel electrophoresis, dye exclusion assays, atomic force microscopy, and Monte Carlo simulations, it is shown that the steric impact of neutral peptide tails is to hinder the formation of DNA-G2 polyplexes composed of multiple DNA chains. If the tails are negatively charged, which results in overall neutral G2 conjugates, then the interaction of G2 with DNA is hindered. Increasing the net positive charge of the tails resulted in the complexation capacity of G2 with the DNA being restored. While DNA complexation is obtained for a similar net charge balance for G2 and G2 conjugates with positive tails, fewer of the latter are required to achieve a comparable condensation degree. Furthermore, it is shown that about 40% of the DNA remains accessible to binding by small molecules. Overall, this shows that tuning the net charge of peptide tails conjugated to PAMAM dendrimers offers a handle to control the complexation capacity of DNA, which can be explored as a novel route for optimization as gene delivery vehicles.

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.

Single-molecule characterization of target search dynamics of DNA-binding proteins in DNA-condensed droplets

Target search models of DNA-binding proteins in cells typically consider search mechanisms that include 3D diffusion and 1D sliding, which can be characterized by single-molecule tracking on DNA. However, the finding of liquid droplets of DNA and nuclear components in cells cast doubt on extrapolation from the behavior in ideal non-condensed DNA conditions to those in cells. In this study, we investigate the target search behavior of DNA-binding proteins in reconstituted DNA-condensed droplets using single-molecule fluorescence microscopy. To mimic nuclear condensates, we reconstituted DNA-condensed droplets using dextran and PEG polymers. In the DNA-condensed droplets, we measured the translational movement of four DNA-binding proteins (p53, Nhp6A, Fis and Cas9) and p53 mutants possessing different structures, sizes, and oligomeric states. Our results demonstrate the presence of fast and slow mobility modes in DNA-condensed droplets for the four DNA-binding proteins. The slow mobility mode capability is correlated strongly to the molecular size and the number of DNA-binding domains on DNA-binding proteins, but only moderately to the affinity to single DNA segments in non-condensed conditions. The slow mobility mode in DNA-condensed droplets is interpreted as a multivalent interaction mode of the DNA-binding protein to multiple DNA segments.

Coupling of binding and differential subdomain folding of the intrinsically disordered transcription factor CREB

The cyclic AMP response element binding protein (CREB) contains a basic leucine zipper motif (bZIP) that forms a coiled coil structure upon dimerization and specific DNA binding. Although this state is well characterized, key features of CREB bZIP binding and folding are not well understood. We used single-molecule Förster resonance energy transfer (smFRET) to probe conformations of CREB bZIP subdomains. We found differential folding of the basic region and leucine zipper in response to different binding partners; a strong and previously unreported DNA-independent dimerization affinity; folding upon binding to nonspecific DNA; and evidence of long-range interdomain interactions in full-length CREB that modulate DNA binding. These studies provide new insights into DNA binding and dimerization and have implications for CREB function.

Mapping Interactions of the Intrinsically Disordered C-Terminal Regions of Tetrameric p53 by Segmental Isotope Labeling and NMR

The C-terminal region of the tumor suppressor protein p53 contains three domains, nuclear localization signal (NLS), tetramerization domain (TET), and C-terminal regulatory domain (CTD), which are essential for p53 function. Characterization of the structure and interactions of these domains within full-length p53 has been limited by the overall size and flexibility of the p53 tetramer. Using trans-intein splicing, we have generated full-length p53 constructs in which the C-terminal region is isotopically labeled with 15N for NMR analysis, allowing us to obtain atomic-level information on the C-terminal domains in the context of the full-length protein. Resonances of NLS and CTD residues have narrow linewidths, showing that these regions are largely solvent-exposed and dynamically disordered, whereas resonances from the folded TET are broadened beyond detection. Two regions of the CTD, spanning residues 369-374 and 381-388 and with high lysine content, make dynamic and sequence-independent interactions with DNA in regions that flank the p53 recognition element. The population of DNA-bound states increases as the length of the flanking regions is extended up to approximately 20 base pairs on either side of the recognition element. Acetylation of K372, K373, and K382, using a construct of the transcriptional coactivator CBP containing the TAZ2 and acetyltransferase domains, inhibits interaction of the CTD with DNA. This work provides high-resolution insights into the behavior of the intrinsically disordered C-terminal regions of p53 within the full-length tetramer and the molecular basis by which the CTD mediates DNA binding and specificity.

Quantifying Cellular Repair, Misrepair and Apoptosis Induced by Boron Ions, Gamma Rays and PRIMA-1 Using the RHR Formulation

The recent interaction cross-section-based formulation for radiation-induced direct cellular inactivation, mild and severe sublethal damage, DNA-repair and cell survival have been developed to accurately describe cellular repair, misrepair and apoptosis in TP53 wild-type and mutant cells. The principal idea of this new non-homologous repairable-homologous repairable (RHR) damage formulation is to separately describe the mild damage that can be rapidly handled by the most basic repair processes including the non-homologous end joining (NHEJ), and more complex damage requiring longer repair times and high-fidelity homologous recombination (HR) repair. Taking the interaction between these two key mammalian DNA repair processes more accurately into account has significantly improved the method as indicated in the original publication. Based on the principal mechanisms of 7 repair and 8 misrepair processes presently derived, it has been possible to quite accurately describe the probability that some of these repair processes when unsuccessful can induce cellular apoptosis with increasing doses of γrays, boron ions and PRIMA-1. Interestingly, for all LETs studied (≈0.3-160 eV/nm) the increase in apoptosis saturates when the cell survival reaches about 10% and the fraction of un-hit cells is well below the 1% level. It is shown that most of the early cell kill for low-to-medium LETs are due to apoptosis since the cell survival as well as the non-apoptotic cells agree very well at low doses and other death processes dominate beyond D > 1 Gy. The low-dose apoptosis is due to the fact that the full activation of the checkpoint kinases ATM and Chk2 requires >8 and >18 DSBs per cell to phosphorylate p53 at serine 15 and 20. Therefore, DNA repair is not fully activated until well after 1/2 Gy, and the cellular response may be apoptotic by default before the low-dose hyper sensitivity (LDHS) is replaced by an increased radiation tolerance as the DNA repair processes get maximal efficiency. In effect, simultaneously explaining the LDHS and inverse dose rate phenomena. The partial contributions by the eight newly derived misrepair processes was determined so they together accurately described the experimental apoptosis induction data for γ rays and boron ions. Through these partial misrepair contributions it was possible to predict the apoptotic response based solely on carefully analyzed cell survival data, demonstrating the usefulness of an accurate DNA repair-based cell survival approach. The peak relative biological effectiveness (RBE) of the boron ions was 3.5 at 160 eV/nm whereas the analogous peak relative apoptotic effectiveness (RAE) was 3.4 but at 40 eV/nm indicating the clinical value of the lower LET light ion (15 \le {\rm{LET}} \le 55{\rm{\ eV}}/{\rm{nm}},{\rm{\ }}2 \le Z \le 5) in therapeutic applications to maximize tumor apoptosis and senescence. The new survival expressions were also applied on mouse embryonic fibroblasts with key knocked-out repair genes, showing a good agreement between the principal non-homologous and homologous repair terms and also a reasonable prediction of the associated apoptotic induction. Finally, the formulation was used to estimate the increase in DNA repair and apoptotic response in combination with the mutant p53 reactivating compound PRIMA-1 and γ rays, indicating a 10-2 times increase in apoptosis with 5 μM of the compound reaching apoptosis levels not far from peak apoptosis boron ions in a TP53 mutant cell line. To utilize PRIMA-1 induced apoptosis and cellular sensitization for reactive oxygen species (ROS), concomitant biologically optimized radiation therapy is proposed to maximize the complication free tumor cure for the multitude of TP53 mutant tumors seen in the clinic. The experimental data also indicated the clinically very important high-absorbed dose ROS effect of PRIMA-1.

Rapid recruitment of p53 to DNA damage sites directs DNA repair choice and integrity

Our work focuses on the critical longstanding question of the nontranscriptional role of p53 in tumor suppression. We demonstrate here that poly(ADP-ribose) polymerase (PARP)–dependent modification of p53 enables rapid recruitment of p53 to damage sites, where it in turn directs early repair pathway selection. Specifically, p53-mediated recruitment of 53BP1 at early time points promotes nonhomologous end joining over the more error-prone microhomology end-joining. Similarly, p53 directs nucleotide excision repair by mediating DDB1 recruitment. This property of p53 also correlates with tumor suppression in vivo. Our study provides mechanistic insight into how certain transcriptionally deficient p53 mutants may retain tumor-suppressive functions through regulating the DNA damage response.

p53 is primarily known as a downstream transcriptional effector in the DNA damage-response cascade. We report that endogenous p53 rapidly accumulates at DNA damage sites within 2 s of UVA microirradiation. The kinetics of p53 recruitment mimics those of known DNA damage-response proteins, such as Ku70 and poly(ADP-ribose) polymerase (PARP), and precedes recruitment of Nbs1, 53BP1, and DDB1. Mutations in the DNA-binding and C-terminal domains significantly suppress this rapid recruitment. The C-terminal domain of p53 contains key residues for PARP interaction that are required for rapid recruitment of p53 to DNA damage sites, as is PARP-dependent modification. The presence of p53 at damage sites influences the recruitment kinetics of 53BP1 and DDB1 and directs the choice of nonhomologous end joining repair (NHEJ) and nucleotide excision repair. Mutations that suppressed rapid recruitment of p53 promoted error-prone alternative end-joining (alt-NHEJ) and inhibited nucleotide excision repair. Our finding that p53 is a critical early responder to DNA damage stands in contrast with its extensively studied role as a downstream transcriptional regulator in DNA damage repair. We highlight an unrecognized role of p53 in directing DNA repair dynamics and integrity and suggest a parallel mode of p53 tumor suppression apart from its function as a transcription factor.

Dynamic transcription regulation at the single-molecule level

Cell fate changes during development, differentiation, and reprogramming are largely controlled at the transcription level. The DNA-binding transcription factors (TFs) often act in a combinatorial fashion to alter chromatin states and drive cell type-specific gene expression. Recent advances in fluorescent microscopy technologies have enabled direct visualization of biomolecules involved in the process of transcription and its regulatory events at the single-molecule level in living cells. Remarkably, imaging and tracking individual TF molecules at high temporal and spatial resolution revealed that they are highly dynamic in searching and binding cognate targets, rather than static and binding constantly. In combination with investigation using techniques from biochemistry, structure biology, genetics, and genomics, a more well-rounded view of transcription regulation is emerging. In this review, we briefly cover the technical aspects of live-cell single-molecule imaging and focus on the biological relevance and interpretation of the single-molecule dynamic features of transcription regulatory events observed in the native chromatin environment of living eukaryotic cells. We also discuss how these dynamic features might shed light on mechanistic understanding of transcription regulation.

Single-Molecule Microscopy Meets Molecular Dynamics Simulations for Characterizing the Molecular Action of Proteins on DNA and in Liquid Condensates

DNA-binding proteins trigger various cellular functions and determine cellular fate. Before performing functions such as transcription, DNA repair, and DNA recombination, DNA-binding proteins need to search for and bind to their target sites in genomic DNA. Under evolutionary pressure, DNA-binding proteins have gained accurate and rapid target search and binding strategies that combine three-dimensional search in solution, one-dimensional sliding along DNA, hopping and jumping on DNA, and intersegmental transfer between two DNA molecules. These mechanisms can be achieved by the unique structural and dynamic properties of these proteins. Single-molecule fluorescence microscopy and molecular dynamics simulations have characterized the molecular actions of DNA-binding proteins in detail. Furthermore, these methodologies have begun to characterize liquid condensates induced by liquid-liquid phase separation, e.g., molecular principles of uptake and dynamics in droplets. This review discusses the molecular action of DNA-binding proteins on DNA and in liquid condensate based on the latest studies that mainly focused on the model protein p53.

Testing mechanisms of DNA sliding by architectural DNA-binding proteins: dynamics of single wild-type and mutant protein molecules in vitro and in vivo

Architectural DNA-binding proteins (ADBPs) are abundant constituents of eukaryotic or bacterial chromosomes that bind DNA promiscuously and function in diverse DNA reactions. They generate large conformational changes in DNA upon binding yet can slide along DNA when searching for functional binding sites. Here we investigate the mechanism by which ADBPs diffuse on DNA by single-molecule analyses of mutant proteins rationally chosen to distinguish between rotation-coupled diffusion and DNA surface sliding after transient unbinding from the groove(s). The properties of yeast Nhp6A mutant proteins, combined with molecular dynamics simulations, suggest Nhp6A switches between two binding modes: a static state, in which the HMGB domain is bound within the minor groove with the DNA highly bent, and a mobile state, where the protein is traveling along the DNA surface by means of its flexible N-terminal basic arm. The behaviors of Fis mutants, a bacterial nucleoid-associated helix-turn-helix dimer, are best explained by mobile proteins unbinding from the major groove and diffusing along the DNA surface. Nhp6A, Fis, and bacterial HU are all near exclusively associated with the chromosome, as packaged within the bacterial nucleoid, and can be modeled by three diffusion modes where HU exhibits the fastest and Fis the slowest diffusion.

Engineering of the genome editing protein Cas9 to slide along DNA

The genome editing protein Cas9 faces engineering challenges in improving off-target DNA cleavage and low editing efficiency. In this study, we aimed to engineer Cas9 to be able to slide along DNA, which might facilitate genome editing and reduce off-target cleavage. We used two approaches to achieve this: reducing the sliding friction along DNA by removing the interactions of Cas9 residues with DNA and facilitating sliding by introducing the sliding-promoting tail of Nhp6A. Seven engineered mutants of Cas9 were prepared, and their performance was tested using single-molecule fluorescence microscopy. Comparison of the mutations enabled the identification of key residues of Cas9 to enhance the sliding along DNA in the presence and absence of single guide RNA (sgRNA). The attachment of the tail to Cas9 mutants enhanced sliding along DNA, particularly in the presence of sgRNA. Together, using the proposed approaches, the sliding ability of Cas9 was improved up to eightfold in the presence of sgRNA. A sliding model of Cas9 and its engineering action are discussed herein.

What Are the Molecular Requirements for Protein Sliding along DNA?

DNA-binding proteins rely on linear diffusion along the longitudinal DNA axis, supported by their nonspecific electrostatic affinity for DNA, to search for their target recognition sites. One may therefore expect that the ability to engage in linear diffusion along DNA is universal to all DNA-binding proteins, with the detailed biophysical characteristics of that diffusion differing between proteins depending on their structures and functions. One key question is whether the linear diffusion mechanism is defined by translation coupled with rotation, a mechanism that is often termed sliding. We conduct coarse-grained and atomistic molecular dynamics simulations to investigate the minimal requirements for protein sliding along DNA. We show that coupling, while widespread, is not universal. DNA-binding proteins that slide along DNA transition to uncoupled translation–rotation (i.e., hopping) at higher salt concentrations. Furthermore, and consistently with experimental reports, we find that the sliding mechanism is the less dominant mechanism for some DNA-binding proteins, even at low salt concentrations. In particular, the toroidal PCNA protein is shown to follow the hopping rather than the sliding mechanism.

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.

Protein-DNA target search relies on quantum walk

Protein-DNA interactions play a fundamental role in all life systems. A critical issue of such interactions is given by the strategy of protein search for specific targets on DNA. The mechanisms by which the protein are able to find relatively small cognate sequences, typically 15-20 base pairs (bps) for repressors, and 4-6 bps for restriction enzymes among the millions of bp of non-specific chromosomal DNA have hardly engaged researchers for decades. Recent experimental studies have generated new insights on the basic processes of protein-DNA interactions evidencing the underlying complex dynamic phenomena involved, which combine three-dimensional and one-dimensional motion along the DNA chain. It has been demonstrated that protein molecules have an extraordinary ability to find the target very quickly on the DNA chain, in some cases, with two orders of magnitude faster than the diffusion limit. This unique property of protein-DNA search mechanism is known as facilitated diffusion. Several theoretical mechanisms have been suggested to describe the origin of facilitated diffusion. However, none of such models currently has the ability to fully describe the protein search strategy. In this paper, we suggest that the ability of proteins to identify consensus sequences on DNA is based on the entanglement of π-π electrons between DNA nucleotides and protein amino acids. The π-π entanglement is based on Quantum Walk (QW), through Coin-position entanglement (CPE). First, the protein identifies a dimer belonging to the consensus sequence, and localize a π on such dimer, hence, the other π electron scans the DNA chain until the sequence is identified. Focusing on the example of recognition of consensus sequences of EcoRV or EcoRI, we will describe the quantum features of QW on protein-DNA complexes during the search strategy, such as walker quadratic spreading on a coherent superposition of different vertices and environment-supported long-time survival probability of the walker. We will employ both discrete- or continuous-time versions of QW. Biased and unbiased classical Random Walk (CRW) have been used for a long time to describe the Protein-DNA search strategy. QW, the quantum version of CRW, has been widely studied for its applications in quantum information applications. In our biological application, the walker (the protein) resides at a vertex in a graph (the DNA structural topology). Differently to CRW, where the walker moves randomly, the quantum walker can hop along the edges in the graph to reach other vertices entering coherently a superposition across different vertices spreading quadratically faster than CRW analogous evidencing the typical speed up features of the QW. When applied to a protein-DNA target search problem, QW gives the possibility to achieve the experimental diffusional motion of proteins over diffusion classical limits experienced along DNA chains exploiting quantum features such as CPE and long-time survival probability supported by the environment. In turn, we come to the conclusion that, under quantum picture, the protein search strategy does not distinguish between one-dimensional (1D) and three-dimensional (3D) cases.

The disordered DNA-binding domain of p53 is indispensable for forming an encounter complex to and jumping along DNA

The tumor suppressor p53 utilizes a facilitated diffusion mechanism to search for and bind to target DNA sequences. Sub-millisecond single-molecule fluorescence tracking demonstrated that p53 forms a short-lived encounter complex to DNA then converts to the long-lived complex that can move and jump along DNA during the target search. To reveal the role of each DNA-binding domain of p53 in these processes, we investigated two p53 mutants lacking either of two DNA-binding domains; structured core and disordered C-terminal domains, using sub-millisecond single-molecule fluorescence microscopy. We found that the C-terminal domain is required for the encounter complex formation and conversion to the long-lived complex. The long-lived complex is stabilized by the core domain as well as the C-terminal domain. Furthermore, only the C-terminal domain participates in the jump of p53 along DNA at a high salt concentration. We propose that the flexible C-terminal domain of p53 is twined around DNA, which can form the encounter complex, convert to the long-lived complex, and enable p53 to land on DNA after the jump.