Asymmetrical roles of zinc fingers in dynamic DNA-scanning process by the inducible transcription factor Egr-1

Role of Shape Deformation of DNA-Binding Sites in Regulating the Efficiency and Specificity in Their Recognition by DNA-Binding Proteins

Accurate transcription of genetic information is crucial, involving precise recognition of the binding motifs by DNA-binding proteins. While some proteins rely on short-range hydrophobic and hydrogen bonding interactions at binding sites, others employ a DNA shape readout mechanism for specific recognition. In this mechanism, variations in DNA shape at the binding motif resulted from either inherent flexibility or binding of proteins at adjacent sites are sensed and capitalized by the searching proteins to locate them specifically. Through extensive computer simulations, we investigate both scenarios to uncover the underlying mechanism and origin of specificity in the DNA shape readout mechanism. Our findings reveal that deformation in shape at the binding motif creates an entropy funnel, allowing information about altered shapes to manifest as fluctuations in minor groove widths. This signal enhances the efficiency of nonspecific search of nearby proteins by directing their movement toward the binding site, primarily driven by a gain in entropy. We propose this as a generic mechanism for DNA shape readout, where specificity arises from the alignment between the molecular frustration of the searching protein and the ruggedness of the entropic funnel governed by molecular features of the protein and arrangement of the DNA bases at the binding site, respectively.

Design, synthesis, and biological evaluation of 2-(2-oxoindolin-3-ylidene)hydrazinecarbothioamides as a potential EGR-1 inhibitor for targeted therapy of atopic dermatitis

Atopic dermatitis is a chronic inflammatory skin disease characterized by intense itching and frequent skin barrier dysfunctions. EGR-1 is a transcription factor that aggravates the pathogenesis of atopic dermatitis by promoting the production of various inflammatory cytokines. Three 2-(2-oxoindolin-3-ylidene)hydrazinecarbothioamides (IT21, IT23, and IT25) were identified as novel inhibitors of EGR-1 DNA-binding activity. In silico docking experiments were performed to elucidate the binding conditions of the EGR-1 zinc-finger (ZnF) DNA-binding domain. Electrophoretic mobility shift assays confirmed the targeted binding effect on the EGR-1 ZnF DNA-binding domain, leading to dose-dependent dissociation of the EGR-1-DNA complex. At the functional cellular level, IT21, IT23, and IT25 effectively reduced mRNA expression of TNFα-induced EGR-1-regulated inflammatory genes, particularly in HaCaT keratinocytes inflamed by TNFα. In the in vivo efficacy study, IT21, IT23, and IT25 demonstrated the potential to alleviate atopic dermatitis-like skin lesions in the ear skin of BALB/c mice. These findings suggest that targeting the EGR-1 ZnF DNA-binding domain with 2-(2-oxoindolin-3-ylidene)hydrazinecarbothioamide derivatives (IT21, IT23, and IT25) could serve as lead compounds for the development of potential therapeutic agents against inflammatory skin disorders, including atopic dermatitis.

A mechanistic insight on how Compromised Hydrolysis of Triacylglycerol 7 (CHT7) restrains the involvement of it's CXC domain from quiescence repression

CHT7 is a regulator of quiescence repression in Chlamydomonas reinhardtii. Initially, CHT7's repression activity was thought to be managed by its DNA-binding CXC domain. Later, it was found that the CHT7-CXC domain is dispensable for CHT7's activities. Rather, CHT7's predicted protein domains were proposed to be involved in regulation activities by binding to other repressors in the cell. Yet, it remains unclear why and how CHT7 refrains its CXC domain from participating in any transcriptional activities. The question becomes more intriguing, since CXC binding regions are available in promoter regions of some of the misregulated genes in CHT7 mutant (cht7). Through biophysical experiments and molecular dynamics approaches, we studied the DNA recognition behavior of CHT7-CXC. The results indicate that this domain possesses sequence selectivity due to the differential binding abilities of its subdomains. Further, to understand if the case is that CXC loses its DNA binding capabilities in the vicinity of other repressors, we examined CHT7-CXC's DNA binding stability under the spatial constraint conditions created through fusing CHT7-CXC with AsLOV2. The results show limited ability of CHT7-CXC to withstand steric forces and provide insights to why and how algal cells may hold back CHT7-CXC's indulgence in quiescence repression. CLASSIFICATIONS: Biological Sciences, Biophysics and Computational Biology.

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.

Negatively charged, intrinsically disordered regions can accelerate target search by DNA-binding proteins

In eukaryotes, many DNA/RNA-binding proteins possess intrinsically disordered regions (IDRs) with large negative charge, some of which involve a consecutive sequence of aspartate (D) or glutamate (E) residues. We refer to them as D/E repeats. The functional role of D/E repeats is not well understood, though some of them are known to cause autoinhibition through intramolecular electrostatic interaction with functional domains. In this work, we investigated the impacts of D/E repeats on the target DNA search kinetics for the high-mobility group box 1 (HMGB1) protein and the artificial protein constructs of the Antp homeodomain fused with D/E repeats of varied lengths. Our experimental data showed that D/E repeats of particular lengths can accelerate the target association in the overwhelming presence of non-functional high-affinity ligands ('decoys'). Our coarse-grained molecular dynamics (CGMD) simulations showed that the autoinhibited proteins can bind to DNA and transition into the uninhibited complex with DNA through an electrostatically driven induced-fit process. In conjunction with the CGMD simulations, our kinetic model can explain how D/E repeats can accelerate the target association process in the presence of decoys. This study illuminates an unprecedented role of the negatively charged IDRs in the target search process.

Engineering and Design of Programmable Genome Editors

Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.

Zinc Fingers 10 and 11 of Miz-1 undergo conformational exchange to achieve specific DNA binding

Miz-1 (ZBTB17) is a poly-zinc finger BTB/POZ transcription factor with 12 consecutive C2H2 zinc fingers (ZFs) that binds transcriptional start sites (TSSs) to regulate the expression of genes involved in cell development and proliferation. As of now, it is not known which of the 12 consecutive ZFs are responsible for the recognition of the 24 base pair consensus sequence found at these TSSs. Evidence suggests ZFs 7-12 plays this role. We provide validation for this and describe the structural and dynamical characterization of unprecedented conformational exchange in the linker between ZFs 10 and 11. This conformational exchange uncouples ZFs 7-10 from 11 and 12 and promotes a scanning-recognition mechanism through which the two segments cooperate to bind two sub-sites at both ends of the consensus. We further show that this can result in the coiling of TSSs as part of Miz-1's mechanism of transcriptional transactivation.

Single molecule characterization of the binding kinetics of a transcription factor and its modulation by DNA sequence and methylation

The interaction of transcription factors with their response elements in DNA is emerging as a highly complex process, whose characterization requires measuring the full distribution of binding and dissociation times in a well-controlled assay. Here, we present a single-molecule assay that exploits the thermal fluctuations of a DNA hairpin to detect the association and dissociation of individual, unlabeled transcription factors. We demonstrate this new approach by following the binding of Egr1 to its consensus motif and the three binding sites found in the promoter of the Lhb gene, and find that both association and dissociation are modulated by the 9 bp core motif and the sequences around it. In addition, CpG methylation modulates the dissociation kinetics in a sequence and position-dependent manner, which can both stabilize or destabilize the complex. Together, our findings show how variations in sequence and methylation patterns synergistically extend the spectrum of a protein's binding properties, and demonstrate how the proposed approach can provide new insights on the function of transcription factors.

The role of zinc finger linkers in zinc finger protein binding to DNA

Zinc finger proteins (ZFP) play important roles in cellular processes. The DNA binding region of ZFP consists of 3 zinc finger DNA binding domains connected by amino acid linkers, the sequence TGQKP connects ZF1 and ZF2, and TGEKP connects ZF2 with ZF3. Linkers act to tune the zinc finger protein in the right position to bind its DNA target, the type of amino acid residues and length of linkers reflect on ZF1-ZF2-ZF3 interactions and contribute to the search and recognition process of ZF protein to its DNA target. Linker mutations and the affinity of the resulting mutants to specific and nonspecific DNA targets were studied by MD simulations and MM_GB(PB)SA. The affinity of mutants to DNA varied with type and position of amino acid residue. Mutation of K in TGQKP resulted in loss in affinity due to the loss of positive K interaction with phosphates, mutation of G showed loss in affinity to DNA, WT protein and all linker mutants showed loss in affinity to a nonspecific DNA target, this finding confirms previous reports which interpreted this loss in affinity as due to ZF1 having an anchoring role, and ZF3 playing an explorer role in the binding mechanism. The change in ZFP-DNA affinity with linker mutations is discussed in view of protein structure and role of linker residues in binding.

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.

Conformational Plasticity and DNA-Binding Specificity of the Eukaryotic Transcription Factor Pax5

The eukaryotic transcription factor Pax5 has a DNA-binding Paired domain composed of two independent helical bundle subdomains joined by a flexible linker. Previously, we showed distinct biophysical properties of the N-terminal (NTD) and C-terminal (CTD) subdomains, with implications for how these two regions cooperate to distinguish nonspecific and cognate DNA sites [Perez-Borrajero, C., et al. (2016) J. Mol. Biol. 428, 2372-2391]. In this study, we combined experimental methods and molecular dynamics (MD) simulations to dissect the mechanisms underlying the functional differences between the Pax5 subdomains. Both subdomains showed a similar dependence of DNA-binding affinity on ionic strength. However, due to a greater contribution of non-ionic interactions, the NTD bound its cognate DNA half-site with an affinity approximately 10-fold higher than that of the CTD with its half-site. These interactions involve base-mediated contacts as evidenced by nuclear magnetic resonance spectroscopy-monitored chemical shift perturbations. Isothermal titration calorimetry revealed that favorable enthalpic and compensating unfavorable entropic changes were substantially larger for DNA binding by the NTD than by the CTD. Complementary MD simulations indicated that the DNA recognition helix H3 of the NTD is particularly flexible in the absence of DNA and undergoes the largest changes in conformational dynamics upon binding. Overall, these data suggest that the differences observed for the subdomains of Pax5 are due to the coupling of DNA binding with dampening of motions in the NTD required for specific base contacts. Thus, the conformational plasticity of the Pax5 Paired domain underpins the differing roles of its subdomains in association with nonspecific versus cognate DNA sites.

Prerecognition Diffusion Mechanism of Human DNA Mismatch Repair Proteins along DNA: Msh2-Msh3 versus Msh2-Msh6

DNA mismatch repair (MMR) is an important postreplication process that eliminates mispaired or unpaired nucleotides to ensure genomic replication fidelity. In humans, Msh2-Msh6 and Msh2-Msh3 are the two mismatch repair initiation factors that recognize DNA lesions. While X-ray crystal structures exist for these proteins in complex with DNA lesions, little is known about their structures during the initial search along nonspecific double-stranded DNA, because they are short-lived and difficult to determine experimentally. In this study, various computational approaches were used to sidestep these difficulties. All-atom and coarse-grained simulations based on the crystal structures of Msh2-Msh3 and Msh2-Msh6 showed no translation along the DNA, suggesting that the initial search conformation differs from the lesion-bound crystal structure. We modeled probable search-mode structures of MSH proteins and showed, using coarse-grained molecular dynamics simulations, that they can perform rotation-coupled diffusion on DNA, which is a suitable and efficient search mechanism for their function and one predicted earlier by fluorescence resonance energy transfer and fluorescence microscopy studies. This search mechanism is implemented by electrostatic interactions among the mismatch-binding domain (MBD), the clamp domains, and the DNA backbone. During simulations, their diffusion rate did not change significantly with an increasing salt concentration, which is consistent with observations from experimental studies. When the gap between their DNA-binding clamps was increased, Msh2-Msh3 diffused mostly via the clamp domains while Msh2-Msh6 still diffused using the MBD, reproducing the experimentally measured lower diffusion coefficient of Msh2-Msh6. Interestingly, Msh2-Msh3 was capable of dissociating from the DNA, whereas Msh2-Msh6 always diffused on the DNA duplex. This is consistent with the experimental observation that Msh2-Msh3, unlike Msh2-Msh6, can overcome obstacles such as nucleosomes. Our models provide a molecular picture of the different mismatch search mechanisms undertaken by Msh2-Msh6 and Msh2-Msh3, despite the similarity of their structures.

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.

Balance between asymmetry and abundance in multi-domain DNA-binding proteins may regulate the kinetics of their binding to DNA

DNA sequences are often recognized by multi-domain proteins that may have higher affinity and specificity than single-domain proteins. However, the higher affinity to DNA might be coupled with slower recognition kinetics. In this study, we address this balance between stability and kinetics for multi-domain Cys₂His₂- (C₂H₂-) type zinc-finger (ZF) proteins. These proteins are the most prevalent DNA-binding domain in eukaryotes and C₂H₂ type zinc-finger proteins (C₂H₂-ZFPs) constitute nearly one-half of all known and predicted transcription factors in human. Extensive contact with DNA via tandem ZF domains confers high stability on the sequence-specific complexes. However, this can limit target search efficiency, especially for low abundance ZFPs. Earlier, we found that asymmetrical distribution of electrostatic charge among the three ZF domains of the low abundance transcription factor Egr-1 facilitates its DNA search process. Here, on a diverse set of 273 human C₂H₂-ZFP comprised of 3–15 tandem ZF domains, we find that, in many cases, electrostatic charge and binding specificity are asymmetrically distributed among the ZF domains so that neighbouring domains have different DNA-binding properties. For proteins containing 3–6 ZF domains, we show that the low abundance proteins possess a higher degree of non-specific asymmetry and vice versa. Our findings suggest that where the electrostatics of tandem ZF domains are similar (i.e., symmetrical), the ZFPs are more abundant to optimize their DNA search efficiency. This study reveals new insights into the fundamental determinants of recognition by C₂H₂-ZFPs of their DNA binding sites in the cellular landscape. The importance of electrostatic asymmetry with respect to binding site recognition by C₂H₂-ZFPs suggests the possibility that it may also be important in other ZFP systems and reveals a new design feature for zinc finger engineering.

Optimal recognition of proteins to DNA is governed by various factors among them the thermodynamics, kinetics and specificity of the protein-DNA complex. Multi-domain DNA-binding proteins are expected to have higher affinity and specificity due to the extensive interface they form with DNA. However, larger interface may result with higher friction when these proteins scan the DNA for the target site via the sliding mechanism. A way to overcome this drawback is to have asymmetry in the protein so that the interface with DNA is smaller. Alternatively, higher abundance can also increase the search speed. Here, using computational analysis of large data set of multi-domain zinc finger DNA-binding proteins, we report a trade-off between asymmetry and abundance.

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other's substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.

Diffusion of ring-shaped proteins along DNA: case study of sliding clamps

Several DNA-binding proteins, such as topoisomerases, helicases and sliding clamps, have a toroidal (i.e. ring) shape that topologically traps DNA, with this quality being essential to their function. Many DNA-binding proteins that function, for example, as transcription factors or enzymes were shown to be able to diffuse linearly (i.e. slide) along DNA during the search for their target binding sites. The protein's sliding properties and ability to search DNA, which often also involves hopping and dissociation, are expected to be different when it encircles the DNA. In this study, we explored the linear diffusion of four ring-shaped proteins of very similar structure: three sliding clamps (PCNA, β-clamp, and the gp45) and the 9-1-1 protein, with a particular focus on PCNA. Coarse-grained molecular dynamics simulations were performed to decipher the sliding mechanism adopted by these ring-shaped proteins and to determine how the molecular properties of the inner and outer ring govern its search speed. We designed in silico variants to dissect the contributions of ring geometry and electrostatics to the sliding speed of ring-shaped proteins along DNA. We found that the toroidal proteins diffuse when they are tilted relative to the DNA axis and able to rotate during translocation, but that coupling between rotation and translocation is quite weak. Their diffusion speed is affected by the shape of the inner ring and, to a lesser extent, by its electrostatic properties. However, breaking the symmetry of the electrostatic potential can result in deviation of the DNA from the center of the ring and cause slower linear diffusion. The findings are discussed in light of earlier computational and experimental studies on the sliding of clamps.

Single-molecule DNA unzipping reveals asymmetric modulation of a transcription factor by its binding site sequence and context

Most functional transcription factor (TF) binding sites deviate from their ‘consensus’ recognition motif, although their sites and flanking sequences are often conserved across species. Here, we used single-molecule DNA unzipping with optical tweezers to study how Egr-1, a TF harboring three zinc fingers (ZF1, ZF2 and ZF3), is modulated by the sequence and context of its functional sites in the Lhb gene promoter. We find that both the core 9 bp bound to Egr-1 in each of the sites, and the base pairs flanking them, modulate the affinity and structure of the protein–DNA complex. The effect of the flanking sequences is asymmetric, with a stronger effect for the sequence flanking ZF3. Characterization of the dissociation time of Egr-1 revealed that a local, mechanical perturbation of the interactions of ZF3 destabilizes the complex more effectively than a perturbation of the ZF1 interactions. Our results reveal a novel role for ZF3 in the interaction of Egr-1 with other proteins and the DNA, providing insight on the regulation of Lhb and other genes by Egr-1. Moreover, our findings reveal the potential of small changes in DNA sequence to alter transcriptional regulation, and may shed light on the organization of regulatory elements at promoters.

Tunable order-disorder continuum in protein-DNA interactions

DNA-binding protein domains (DBDs) sample diverse conformations in equilibrium facilitating the search and recognition of specific sites on DNA over millions of energetically degenerate competing sites. We hypothesize that DBDs have co-evolved to sense and exploit the strong electric potential from the array of negatively charged phosphate groups on DNA. We test our hypothesis by employing the intrinsically disordered DBD of cytidine repressor (CytR) as a model system. CytR displays a graded increase in structure, stability and folding rate on increasing the osmolarity of the solution that mimics the non-specific screening by DNA phosphates. Electrostatic calculations and an Ising-like statistical mechanical model predict that CytR exhibits features of an electric potential sensor modulating its dimensions and landscape in a unique distance-dependent manner, while DNA plays the role of a non-specific macromolecular chaperone. Accordingly, CytR binds its natural half-site faster than the diffusion-controlled limit and even random DNA conforming to an electrostatic-steering binding mechanism. Our work unravels for the first time the synergistic features of a natural electrostatic potential sensor, a novel binding mechanism driven by electrostatic frustration and disorder, and the role of DNA in promoting distance-dependent protein structural transitions critical for switching between specific and non-specific DNA-binding modes.

Quantifying the two-state facilitated diffusion model of protein-DNA interactions

The current report extends the facilitated diffusion model to account for conflict between the search and recognition binding modes adopted by DNA-binding proteins (DBPs) as they search DNA and subsequently recognize and bind to their specific binding site. The speed of the search dynamics is governed by the energetic ruggedness of the protein-DNA landscape, whereas the rate for the recognition process is mostly dictated by the free energy barrier for the transition between the DBP's search and recognition binding modes. We show that these two modes are negatively coupled, such that fast 1D sliding and rapid target site recognition probabilities are unlikely to coexist. Thus, a tradeoff occurs between optimizing the timescales for finding and binding the target site. We find that these two kinetic properties can be balanced to produce a fast timescale for the total target search and recognition process by optimizing frustration. Quantification of the facilitated diffusion model by including a frustration term enables it to explain several experimental observations concerning search and recognition speeds. The extended model captures experimental estimate of the energetic ruggedness of the protein-DNA landscape and predicts how various molecular properties of protein-DNA binding affect recognition kinetics. Particularly, point mutations may change the frustration and so affect protein association with DNA, thus providing a means to modulate protein-DNA affinity by manipulating the protein's association or dissociation reactions.

NMR-based investigations into target DNA search processes of proteins

To perform their function, transcription factors and DNA-repair/modifying enzymes must first locate their targets in the vast presence of nonspecific, but structurally similar sites on genomic DNA. Before reaching their targets, these proteins stochastically scan DNA and dynamically move from one site to another on DNA. Solution NMR spectroscopy provides unique atomic-level insights into the dynamic DNA-scanning processes, which are difficult to gain by any other experimental means. In this review, we provide an introductory overview on the NMR methods for the structural, dynamic, and kinetic investigations of target DNA search by proteins. We also discuss advantages and disadvantages of these NMR methods over other methods such as single-molecule techniques and biochemical approaches.

Role of protein structure and the role of individual fingers in zinc finger protein-DNA recognition: a molecular dynamics simulation study and free energy calculations

Molecular dynamics and MM_GBSA energy calculations on various zinc finger proteins containing three and four fingers bound to their target DNA gave insights into the role of each finger in the DNA binding process as part of the protein structure. The wild type Zif 268 (PDB code: 1AAY) gave a ΔG value of - 76.1 (14) kcal/mol. Zinc fingers ZF1, ZF2 and ZF3 were mutated in one experiment and in another experiment one finger was cut and the rest of the protein was studied for binding. The ΔΔG values for the Zinc Finger protein with both ZF1 and ZF2 mutated was + 80 kcal/mol, while mutating only ZF1 the ΔΔG value was + 52 kcal/mol (relative to the wild type). Cutting ZF3 and studying the protein consisting only of ZF1 linked to ZF2 gave a ΔΔG value of + 68 kcal/mol. Upon cutting ZF1, the resulting ZF2 linked to ZF3 protein gave a ΔΔG value of + 41 kcal/mol. The above results shed light on the importance of each finger in the binding process, especially the role of ZF1 as the anchoring finger followed in importance by ZF2 and ZF3. The energy difference between the binding of the wild type protein Zif268 (1AAY) and that for individual finger binding to DNA according to the formula: ΔΔG_{linkers, otherstructuralfactors} = ΔG_zif268 - (ΔG_F1+F2+F3) gave a value = - 44.5 kcal/mol. This stabilization can be attributed to the contribution of linkers and other structural factors in the intact protein in the DNA binding process. DNA binding energies of variant proteins of the wild type Zif268 which differ in their ZF1 amino acid sequence gave evidence of a good relationship between binding energy and recognition and specificity, this finding confirms the reported vital role of ZF1 in the ZF protein scanning and anchoring to the target DNA sequence. The role of hydrogen bonds in both specific and nonspecific amino acid-DNA contacts is discussed in relation to mutations. The binding energies of variant Zinc Finger proteins confirmed the role of ZF1 in the recognition, specificity and anchoring of the zinc finger protein to DNA.

Roles of conformational disorder and downhill folding in modulating protein-DNA recognition

Transcription factors are thought to efficiently search for their target DNA site via a combination of conventional 3D diffusion and 1D diffusion along the DNA molecule mediated by non-specific electrostatic interactions. This process requires the DNA-binding protein to quickly exchange between a search competent and a target recognition mode, but little is known as to how these two binding modes are encoded in the conformational properties of the protein. Here, we investigate this issue on the engrailed homeodomain (EngHD), a DNA-binding domain that folds ultrafast and exhibits a complex conformational behavior consistent with the downhill folding scenario. We explore the interplay between folding and DNA recognition using a coarse-grained computational model that allows us to manipulate the folding properties of the protein and monitor its non-specific and specific binding to DNA. We find that conformational disorder increases the search efficiency of EngHD by promoting a fast gliding search mode in addition to sliding. When gliding, EngHD remains loosely bound to DNA moving linearly along its length. A partially disordered EngHD also binds more dynamically to the target site, reducing the half-life of the specific complex via a spring-loaded mechanism. These findings apply to all conditions leading to partial disorder. However, we also find that at physiologically relevant temperatures EngHD is well folded and can only obtain the conformational flexibility required to accelerate 1D diffusion when it folds/unfolds within the downhill scenario (crossing a marginal free energy barrier). In addition, the conformational flexibility of native downhill EngHD enables its fast reconfiguration to lock into the specific binding site upon arrival, thereby affording finer control of the on- and off-rates of the specific complex. Our results provide key mechanistic insights into how DNA-binding domains optimize specific DNA recognition through the control of their conformational dynamics and folding mechanism.

Intersegmental transfer of proteins between DNA regions in the presence of crowding

Intersegmental transfer that involves direct relocation of a DNA-binding protein from one nonspecific DNA site to another was previously shown to contribute to speeding up the identification of the DNA target site. This mechanism is promoted when the protein is composed of at least two domains that have different DNA binding affinities and thus show a degree of mobility. In this study, we investigate the effect of particle crowding on the ability of a multi-domain protein to perform intersegmental transfer. We show that although crowding conditions often favor 1D diffusion of proteins along DNA over 3D diffusion, relocation of one of the tethered domains to initiate intersegmental transfer is possible even under crowding conditions. The tendency to perform intersegmental transfer by a multi-domain protein under crowding conditions is much higher for larger crowding particles than smaller ones and can be even greater than under no-crowding conditions. We report that the asymmetry of the two domains is even magnified by the crowders. The observations that crowding supports intersegmental transfer serve as another example that in vivo complexity does not necessarily slow down DNA search kinetics by proteins.

Potential role of DNA methylation as a facilitator of target search processes for transcription factors through interplay with methyl-CpG-binding proteins

Eukaryotic genomes contain numerous non-functional high-affinity sequences for transcription factors. These sequences potentially serve as natural decoys that sequester transcription factors. We have previously shown that the presence of sequences similar to the target sequence could substantially impede association of the transcription factor Egr-1 with its targets. In this study, using a stopped-flow fluorescence method, we examined the kinetic impact of DNA methylation of decoys on the search process of the Egr-1 zinc-finger protein. We analyzed its association with an unmethylated target site on fluorescence-labeled DNA in the presence of competitor DNA duplexes, including Egr-1 decoys. DNA methylation of decoys alone did not affect target search kinetics. In the presence of the MeCP2 methyl-CpG-binding domain (MBD), however, DNA methylation of decoys substantially (∼10-30-fold) accelerated the target search process of the Egr-1 zinc-finger protein. This acceleration did not occur when the target was also methylated. These results suggest that when decoys are methylated, MBD proteins can block them and thereby allow Egr-1 to avoid sequestration in non-functional locations. This effect may occur in vivo for DNA methylation outside CpG islands (CGIs) and could facilitate localization of some transcription factors within regulatory CGIs, where DNA methylation is rare.

Discrete-State Kinetics Model for NMR-Based Analysis of Protein Translocation on DNA at Equilibrium

In the target DNA search process, sequence-specific DNA-binding proteins first nonspecifically bind to DNA and stochastically move from one site to another before reaching their targets. To rigorously assess how the translocation process influences NMR signals from proteins interacting with nonspecific DNA, we incorporated a discrete-state kinetic model for protein translocation on DNA into the McConnell equation. Using this equation, we simulated line shapes of NMR signals from proteins undergoing translocations on DNA through sliding, dissociation/reassociation, and intersegment transfer. Through this analysis, we validated an existing NMR approach for kinetic investigations of protein translocation on DNA, which utilizes NMR line shapes of two nonspecific DNA-protein complexes and their mixture. We found that, despite its use of simplistic two-state approximation neglecting the presence of many microscopic states, the previously proposed NMR approach provides accurate kinetic information on the intermolecular translocations of proteins between two DNA molecules. Interestingly, our results suggest that the same NMR approach can also provide qualitative information about the one-dimensional diffusion coefficient for proteins sliding on DNA.

Optimal Length of Conformational Transition Region in Protein Search for Targets on DNA

The starting point of many fundamental biological processes is associated with protein molecules finding and recognizing specific sites on DNA. However, despite a large number of experimental and theoretical studies on protein search for targets on DNA, many molecular aspects of underlying mechanisms are still not well understood. Experiments show that proteins bound to DNA can switch between slow recognition and fast search conformations. However, from a theoretical point of view, such conformational transitions should slow down the protein search for specific sites on DNA, in contrast to available experimental observations. In addition, experiments indicate that the nucleotide composition near the target site is more symmetrically homogeneous, leading to stronger effective interactions between proteins and DNA at these locations. However, as has been shown theoretically, this should also make the search less efficient, which is not observed. We propose a possible resolution of these problems by suggesting that conformational transitions occur only within a segment around the target where stronger interactions between proteins and DNA are observed. Two theoretical methods, based on continuum and discrete-state stochastic calculations, are developed, allowing us to obtain a comprehensive dynamic description for the protein search process in this system. The existence of an optimal length of the conformational transition zone with the shortest mean search time is predicted.

In-depth study of DNA binding of Cys2His2 finger domains in testis zinc-finger protein

Previously, we identified that both fingers 1 and 2 in the three Cys₂His₂ zinc-finger domains (TZD) of testis zinc-finger protein specifically bind to its cognate DNA; however, finger 3 is non-sequence–specific. To gain insights into the interaction mechanism, here we further investigated the DNA-binding characteristics of TZD bound to non-specific DNAs and its finger segments bound to cognate DNA. TZD in non-specific DNA binding showed smaller chemical shift perturbations, as expected. However, the direction of shift perturbation, change of DNA imino-proton NMR signal, and dynamics on the ¹⁵N backbone atom significantly differed between specific and non-specific binding. Using these unique characteristics, we confirmed that the three single-finger segments (TZD₁, TZD₂ and TZD₃) and the two-finger segment (TZD₂₃) non-specifically bind to the cognate DNA. In comparison, the other two-finger segment (TZD₁₂) binding to the cognate DNA features simultaneous non-specific and semi-specific binding, both slowly exchanged in terms of NMR timescale. The process of TZD binding to the cognate DNA is likely stepwise: initially TZD non-specifically binds to DNA, then fingers 1 and 2 insert cooperatively into the major groove of DNA by semi-specific binding, and finally finger 3 non-specifically binds to DNA, which promotes the specific binding on fingers 1 and 2 and stabilizes the formation of a specific TZD–DNA complex.

The long zinc finger domain of PRDM9 forms a highly stable and long-lived complex with its DNA recognition sequence

PR domain containing protein 9 (PRDM9) is a meiosis-specific, multi-domain protein that regulates the location of recombination hotspots by targeting its DNA recognition sequence for double-strand breaks (DSBs). PRDM9 specifically recognizes DNA via its tandem array of zinc fingers (ZnFs), epigenetically marks the local chromatin by its histone methyltransferase activity, and is an important tether that brings the DNA into contact with the recombination initiation machinery. A strong correlation between PRDM9-ZnF variants and specific DNA motifs at recombination hotspots has been reported; however, the binding specificity and kinetics of the ZnF domain are still obscure. Using two in vitro methods, gel mobility shift assays and switchSENSE, a quantitative biophysical approach that measures binding rates in real time, we determined that the PRDM9-ZnF domain forms a highly stable and long-lived complex with its recognition sequence, with a dissociation halftime of many hours. The ZnF domain exhibits an equilibrium dissociation constant (K _D) in the nanomolar (nM) range, with polymorphisms in the recognition sequence directly affecting the binding affinity. We also determined that alternative sequences (15-16 nucleotides in length) can be specifically bound by different subsets of the ZnF domain, explaining the binding plasticity of PRDM9 for different sequences. Finally, longer binding targets are preferred than predicted from the numbers of ZnFs contacting the DNA. Functionally, a long-lived complex translates into an enzymatically active PRDM9 at specific DNA-binding sites throughout meiotic prophase I that might be relevant in stabilizing the components of the recombination machinery to a specific DNA target until DSBs are initiated by Spo11.

Structural Insights into c-Myc-interacting Zinc Finger Protein-1 (Miz-1) Delineate Domains Required for DNA Scanning and Sequence-specific Binding

c-Myc-interacting zinc finger protein-1 (Miz-1) is a poly-Cys2His2 zinc finger (ZF) transcriptional regulator of many cell cycle genes. A Miz-1 DNA sequence consensus has recently been identified and has also unveiled Miz-1 functions in other cellular processes, underscoring its importance in the cell. Miz-1 contains 13 ZFs, but it is unknown why Miz-1 has so many ZFs and whether they recognize and bind DNA sequences in a typical fashion. Here, we used NMR to deduce the role of Miz-1 ZFs 1-4 in detecting the Miz-1 consensus sequence and preventing nonspecific DNA binding. In the construct containing the first 4 ZFs, we observed that ZFs 3 and 4 form an unusual compact and stable structure that restricts their motions. Disruption of this compact structure by an electrostatically mismatched A86K mutation profoundly affected the DNA binding properties of the WT construct. On the one hand, Miz1-4^WT was found to bind the Miz-1 DNA consensus sequence weakly and through ZFs 1-3 only. On the other hand, the four ZFs in the structurally destabilized Miz1-4^A86K mutant bound to the DNA consensus with a 30-fold increase in affinity (100 nm). The formation of such a thermodynamically stable but nonspecific complex is expected to slow down the rate of DNA scanning by Miz-1 during the search for its consensus sequence. Interestingly, we found that the motif stabilizing the compact structure between ZFs 3 and 4 is conserved and enriched in other long poly-ZF proteins. As discussed in detail, our findings support a general role of compact inter-ZF structures in minimizing the formation of off-target DNA complexes.

The Role of DNA Looping in the Search for Specific Targets on DNA by Multisite Proteins

Many cellular processes involve simultaneous interactions between DNA and protein molecules at several locations. They are regulated and controlled by special protein-DNA complexes, which are known as synaptic complexes or synaptosomes. Because of the multisite nature of involved proteins, it was suggested that during the formation of synaptic complexes DNA loops might appear, but their role is unclear. We developed a theoretical model that allowed us to evaluate the effect of transient DNA loop formation. It is based on a discrete-state stochastic method that explicitly takes into account the free-energy contributions due to the appearance of DNA loops. The formation of the synaptic complexes is viewed as a search for a specific binding site on DNA by the protein molecule already bound to DNA at another location. It was found that the search might be optimized by varying the position of the target and the total length of DNA. Furthermore, the formation of transient DNA loops leads to faster dynamics if it is associated with favorable enthalpic contributions to nonspecific protein-DNA interactions. It is also shown that DNA looping might reduce stochastic noise in the system.

Thermodynamic Additivity for Impacts of Base-Pair Substitutions on Association of the Egr-1 Zinc-Finger Protein with DNA

The transcription factor Egr-1 specifically binds as a monomer to its 9 bp target DNA sequence, GCGTGGGCG, via three zinc fingers and plays important roles in the brain and cardiovascular systems. Using fluorescence-based competitive binding assays, we systematically analyzed the impacts of all possible single-nucleotide substitutions in the target DNA sequence and determined the change in binding free energy for each. Then, we measured the changes in binding free energy for sequences with multiple substitutions and compared them with the sum of the changes in binding free energy for each constituent single substitution. For the DNA variants with two or three nucleotide substitutions in the target sequence, we found excellent agreement between the measured and predicted changes in binding free energy. Interestingly, however, we found that this thermodynamic additivity broke down with a larger number of substitutions. For DNA sequences with four or more substitutions, the measured changes in binding free energy were significantly larger than predicted. On the basis of these results, we analyzed the occurrences of high-affinity sequences in the genome and found that the genome contains millions of such sequences that might functionally sequester Egr-1.

Single-Molecule Characterization of DNA-Protein Interactions Using Nanopore Biosensors

Detection and characterization of nucleic acid-protein interactions, particularly those involving DNA and proteins such as transcription factors, enzymes, and DNA packaging proteins, remain significant barriers to our understanding of genetic regulation. Nanopores are an extremely sensitive and versatile sensing platform for label-free detection of single biomolecules. Analyte molecules are drawn to and through a nanoscale aperture by an electrophoretic force, which acts upon their native charge while in the sensing region of the pore. When the nanopore's diameter is only slightly larger than the biopolymer's cross section (typically a few nm); the latter must translocate through the pore in a linear fashion due to the constricted geometry in this region. These features allow nanopores to interrogate protein-nucleic acids in multiple sensing modes: first, by scanning and mapping the locations of binding sites along an analyte molecule, and second, by probing the strength of the bond between a protein and nucleic acid, using the native charge of the nucleic acid to apply an electrophoretic force to the complex while the protein is geometrically prevented from passing through the nanopore. In this chapter, we describe progress toward nanopore sensing of protein-nucleic acid complexes in the context of both mapping binding sites and performing force spectroscopy to determine the strength of interactions. We conclude by reviewing the strengths and challenges of the nanopore technique in the context of studying DNA-protein interactions.

NMR Scalar Couplings across Intermolecular Hydrogen Bonds between Zinc-Finger Histidine Side Chains and DNA Phosphate Groups

NMR scalar couplings across hydrogen bonds represent direct evidence for the partial covalent nature of hydrogen bonds and provide structural and dynamic information on hydrogen bonding. In this article, we report heteronuclear 15N-31P and 1H-31P scalar couplings across the intermolecular hydrogen bonds between protein histidine (His) imidazole and DNA phosphate groups. These hydrogen-bond scalar couplings were observed for the Egr-1 zinc-finger-DNA complex. Although His side-chain NH protons are typically undetectable in heteronuclear 1H-15N correlation spectra due to rapid hydrogen exchange, this complex exhibited two His side-chain NH signals around 1H 14.3 ppm and 15N 178 ppm at 35 °C. Through various heteronuclear multidimensional NMR experiments, these signals were assigned to two zinc-coordinating His side chains in contact with DNA phosphate groups. The data show that the Nδ1 atoms of these His side chains are protonated and exhibit the 1H-15N cross-peaks. Using heteronuclear 1H, 15N, and 31P NMR experiments, we observed the hydrogen-bond scalar couplings between the His 15Nδ1/1Hδ1 and DNA phosphate 31P nuclei. These results demonstrate the direct involvement of the zinc-coordinating His side chains in the recognition of DNA by the Cys2His2-class zinc fingers in solution.

Sequence heterogeneity accelerates protein search for targets on DNA

The process of protein search for specific binding sites on DNA is fundamentally important since it marks the beginning of all major biological processes. We present a theoretical investigation that probes the role of DNA sequence symmetry, heterogeneity, and chemical composition in the protein search dynamics. Using a discrete-state stochastic approach with a first-passage events analysis, which takes into account the most relevant physical-chemical processes, a full analytical description of the search dynamics is obtained. It is found that, contrary to existing views, the protein search is generally faster on DNA with more heterogeneous sequences. In addition, the search dynamics might be affected by the chemical composition near the target site. The physical origins of these phenomena are discussed. Our results suggest that biological processes might be effectively regulated by modifying chemical composition, symmetry, and heterogeneity of a genome.

Conformational Dynamics and the Binding of Specific and Nonspecific DNA by the Autoinhibited Transcription Factor Ets-1

The affinity of the Ets-1 transcription factor for DNA is autoinhibited by an intrinsically disordered serine-rich region (SRR) and a helical inhibitory module (IM) appended to its winged helix-turn-helix ETS domain. Using NMR spectroscopy, we investigated how Ets-1 recognizes specific versus nonspecific DNA, with a focus on the roles of protein dynamics and autoinhibition in these processes. Upon binding either DNA, the two marginally stable N-terminal helices of the IM predominantly unfold, but still sample partially ordered conformations. Also, on the basis of amide chemical shift perturbation mapping, Ets-1 associates with both specific and nonspecific DNA through the same canonical ETS domain interface. These interactions are structurally independent of the SRR, and thus autoinhibition does not impart DNA-binding specificity. However, relative to the pronounced NMR spectroscopic changes in Ets-1 resulting from specific DNA binding, the spectra of the nonspecific DNA complexes showed conformational exchange broadening and lacked several diagnostic amide and indole signals attributable to hydrogen bonding interactions seen in reported X-ray crystallographic structures of this transcription factor with its cognate DNA sequences. Such differences are highlighted by the chemical shift and relaxation properties of several interfacial lysine and arginine side chains. Collectively, these data support a general model in which Ets-1 interacts with nonspecific DNA via dynamic electrostatic interactions, whereas hydrogen bonding drives the formation of well-ordered complexes with specific DNA.

Crowding on DNA in Protein Search for Targets

Proteins searching and recognizing specific sites on DNA is required for initiating all major biological processes. While the details of the protein search for targets on DNA in purified in vitro systems are reasonably well understood, the situation in real cells is much less clear. The presence of other types of molecules on DNA should prevent reaching the targets, but experiments show that, surprisingly, the molecular crowding on DNA influences the search dynamics much less than expected. We develop a theoretical method that allowed us to clarify the mechanisms of the protein search on DNA in the presence of crowding. It is found that the dimensionality of the search trajectories specifies whether the crowding will affect the target finding. For 3D search pathways it is minimal, while the strongest effect is for 1D search pathways when the crowding particle can block the search. In addition, for 1D search we determined that the critical parameter is a mobility of crowding agents: highly mobile molecules do not affect the search dynamics, while the slow particles can significantly slow down the process. Physical-chemical explanations of the observed phenomena are presented. Our theoretical predictions thus explain the experimental observations, and they are also supported by extensive numerical simulations.

Protein search for multiple targets on DNA

Protein-DNA interactions are crucial for all biological processes. One of the most important fundamental aspects of these interactions is the process of protein searching and recognizing specific binding sites on DNA. A large number of experimental and theoretical investigations have been devoted to uncovering the molecular description of these phenomena, but many aspects of the mechanisms of protein search for the targets on DNA remain not well understood. One of the most intriguing problems is the role of multiple targets in protein search dynamics. Using a recently developed theoretical framework we analyze this question in detail. Our method is based on a discrete-state stochastic approach that takes into account most relevant physical-chemical processes and leads to fully analytical description of all dynamic properties. Specifically, systems with two and three targets have been explicitly investigated. It is found that multiple targets in most cases accelerate the search in comparison with a single target situation. However, the acceleration is not always proportional to the number of targets. Surprisingly, there are even situations when it takes longer to find one of the multiple targets in comparison with the single target. It depends on the spatial position of the targets, distances between them, average scanning lengths of protein molecules on DNA, and the total DNA lengths. Physical-chemical explanations of observed results are presented. Our predictions are compared with experimental observations as well as with results from a continuum theory for the protein search. Extensive Monte Carlo computer simulations fully support our theoretical calculations.

Changes in conformational dynamics of basic side chains upon protein-DNA association

Basic side chains play major roles in recognition of nucleic acids by proteins. However, dynamic properties of these positively charged side chains are not well understood. In this work, we studied changes in conformational dynamics of basic side chains upon protein–DNA association for the zinc-finger protein Egr-1. By nuclear magnetic resonance (NMR) spectroscopy, we characterized the dynamics of all side-chain cationic groups in the free protein and in the complex with target DNA. Our NMR order parameters indicate that the arginine guanidino groups interacting with DNA bases are strongly immobilized, forming rigid interfaces. Despite the strong short-range electrostatic interactions, the majority of the basic side chains interacting with the DNA phosphates exhibited high mobility, forming dynamic interfaces. In particular, the lysine side-chain amino groups exhibited only small changes in the order parameters upon DNA-binding. We found a similar trend in the molecular dynamics (MD) simulations for the free Egr-1 and the Egr-1–DNA complex. Using the MD trajectories, we also analyzed side-chain conformational entropy. The interfacial arginine side chains exhibited substantial entropic loss upon binding to DNA, whereas the interfacial lysine side chains showed relatively small changes in conformational entropy. These data illustrate different dynamic characteristics of the interfacial arginine and lysine side chains.

Solution structure of the 13th C2H2 Zinc Finger of Miz-1

Miz-1 is a BTB/POZ transcription factor that contains 13C2H2 Zinc Finger domains (ZF). Miz-1 transactivates and represses the transcription of a myriad of genes involved in many aspects of the biology of the cell. The detailed molecular interactions through which Miz-1 controls transcription, including its specific DNA binding via its ZF domains, remain to be understood and documented. In our effort to shed light into the structural biology of Miz-1, we have undertaken the determination of the structure of all its ZF and the characterization of their interactions with cognate DNA. The structure of ZF 1 to 10 have already been solved and characterized. Here, we present the structure of the synthetic Miz-1 ZF13 determined by 2D (1)H-(1)H NMR.

Residence Times of Molecular Complexes in Solution from NMR Data of Intermolecular Hydrogen-Bond Scalar Coupling

The residence times of molecular complexes in solution are important for understanding biomolecular functions and drug actions. We show that NMR data of intermolecular hydrogen-bond scalar couplings can yield information on the residence times of molecular complexes in solution. The molecular exchange of binding partners via the breakage and reformation of a complex causes self-decoupling of intermolecular hydrogen-bond scalar couplings, and this self-decoupling effect depends on the residence time of the complex. For protein-DNA complexes, we investigated the salt concentration dependence of intermolecular hydrogen-bond scalar couplings between the protein side-chain (15)N and DNA phosphate (31)P nuclei, from which the residence times were analyzed. The results were consistent with those obtained by (15)Nz-exchange spectroscopy. This self-decoupling-based kinetic analysis is unique in that it does not require any different signatures for the states involved in the exchange, whereas such conditions are crucial for kinetic analyses by typical NMR and other methods.

Influence of quasi-specific sites on kinetics of target DNA search by a sequence-specific DNA-binding protein

Functions of transcription factors require formation of specific complexes at particular sites in cis-regulatory elements of genes. However, chromosomal DNA contains numerous sites that are similar to the target sequences recognized by transcription factors. The influence of such “quasi-specific” sites on functions of the transcription factors is not well understood at present by experimental means. In this work, using fluorescence methods, we have investigated the influence of quasi-specific DNA sites on the efficiency of target location by the zinc finger DNA-binding domain of the inducible transcription factor Egr-1, which recognizes a 9 bp sequence. By stopped-flow assays, we measured the kinetics of Egr-1’s association with a target site on 143 bp DNA in the presence of various competitor DNAs, including nonspecific and quasi-specific sites. The presence of quasi-specific sites on competitor DNA significantly decelerated the target association by the Egr-1 protein. The impact of the quasi-specific sites depended strongly on their affinity, their concentration, and the degree of their binding to the protein. To quantitatively describe the kinetic impact of the quasi-specific sites, we derived an analytical form of the apparent kinetic rate constant for the target association and used it for fitting to the experimental data. Our kinetic data with calf thymus DNA as a competitor suggested that there are millions of high-affinity quasi-specific sites for Egr-1 among the 3 billion bp of genomic DNA. This study quantitatively demonstrates that naturally abundant quasi-specific sites on DNA can considerably impede the target search processes of sequence-specific DNA-binding proteins.

Balancing between affinity and speed in target DNA search by zinc-finger proteins via modulation of dynamic conformational ensemble

Although engineering of transcription factors and DNA-modifying enzymes has drawn substantial attention for artificial gene regulation and genome editing, most efforts focus on affinity and specificity of the DNA-binding proteins, typically overlooking the kinetic properties of these proteins. However, a simplistic pursuit of high affinity can lead to kinetically deficient proteins that spend too much time at nonspecific sites before reaching their targets on DNA. We demonstrate that structural dynamic knowledge of the DNA-scanning process allows for kinetically and thermodynamically balanced engineering of DNA-binding proteins. Our current study of the zinc-finger protein Egr-1 (also known as Zif268) and its nuclease derivatives reveals kinetic and thermodynamic roles of the dynamic conformational equilibrium between two modes during the DNA-scanning process: one mode suitable for search and the other for recognition. By mutagenesis, we were able to shift this equilibrium, as confirmed by NMR spectroscopy. Using fluorescence and biochemical assays as well as computational simulations, we analyzed how the shifts of the conformational equilibrium influence binding affinity, target search kinetics, and efficiency in displacing other proteins from the target sites. A shift toward the recognition mode caused an increase in affinity for DNA and a decrease in search efficiency. In contrast, a shift toward the search mode caused a decrease in affinity and an increase in search efficiency. This accelerated site-specific DNA cleavage by the zinc-finger nuclease, without enhancing off-target cleavage. Our study shows that appropriate modulation of the dynamic conformational ensemble can greatly improve zinc-finger technology, which has used Egr-1 (Zif268) as a major scaffold for engineering.

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Dynamic Equilibria of Short-Range Electrostatic Interactions at Molecular Interfaces of Protein-DNA Complexes

Intermolecular ion pairs (salt bridges) are crucial for protein-DNA association. For two protein-DNA complexes, we demonstrate that the ion pairs of protein side-chain NH3+ and DNA phosphate groups undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. While the crystal structures of the complexes show only the solvent-separated ion pair (SIP) state for some interfacial lysine side chains, our NMR hydrogen-bond scalar coupling data clearly indicate the presence of the contact ion pair (CIP) state for the same residues. The 0.6-μs molecular dynamics (MD) simulations confirm dynamic transitions between the CIP and SIP states. This behavior is consistent with our NMR order parameters and scalar coupling data for the lysine side chains. Using the MD trajectories, we also analyze the free energies of the CIP-SIP equilibria. This work illustrates the dynamic nature of short-range electrostatic interactions in DNA recognition by proteins.

A multiscale approach to simulating the conformational properties of unbound multi-C₂H₂ zinc finger proteins

The conformational properties of unbound multi-Cys2 His2 (mC2H2) zinc finger proteins, in which zinc finger domains are connected by flexible linkers, are studied by a multiscale approach. Three methods on different length scales are utilized. First, atomic detail molecular dynamics simulations of one zinc finger and its adjacent flexible linker confirmed that the zinc finger is more rigid than the flexible linker. Second, the end-to-end distance distributions of mC2H2 zinc finger proteins are computed using an efficient atomistic pivoting algorithm, which only takes excluded volume interactions into consideration. The end-to-end distance distribution gradually changes its profile, from left-tailed to right-tailed, as the number of zinc fingers increases. This is explained by using a worm-like chain model. For proteins of a few zinc fingers, an effective bending constraint favors an extended conformation. Only for proteins containing more than nine zinc fingers, is a somewhat compacted conformation preferred. Third, a mesoscale model is modified to study both the local and the global conformational properties of multi-C2H2 zinc finger proteins. Simulations of the CCCTC-binding factor (CTCF), an important mC2H2 zinc finger protein for genome spatial organization, are presented.