Stopped-flow fluorescence kinetic study of protein sliding and intersegment transfer in the target DNA search process

Pharmacokinetic-Pharmacodynamic Analysis of pH-Responsive Doxorubicin-Releasing Micelles with Anticancer Activity

This study aimed to investigate the effect of in vivo pH-responsive doxorubicin (DOX) release and the targetability of pilot molecules in folic acid (FA)-modified micelles using a pharmacokinetic-pharmacodynamic (PK-PD) model. The time profiles of intratumoral DOX concentrations in Walker256 tumor-bearing rats were monitored using a microdialysis probe, followed by compartmental analysis, to evaluate intratumoral tissue pharmacokinetics. Maximal DOX was released from micelles 350 min after the administration of pH-responsive DOX-releasing micelles. However, FA modification of the micelles shortened the time to peak drug concentration to 150 min. Additionally, FA modification resulted in a 27-fold increase in the tumor inflow rate constant. Walker256 tumor-bearing rats were subsequently treated with DOX, pH-responsive DOX-releasing micelles, and pH-responsive DOX-releasing FA-modified micelles to monitor the tumor growth-time profiles. An intratumoral threshold concentration of DOX (55-64 ng/g tumor) was introduced into the drug efficacy compartment to construct a PD model, followed by PK-PD analysis of the tumor growth-time profiles. Similar results of threshold concentration and drug potency of DOX were obtained across all three formulations. Cell proliferation was delayed as the drug delivery ability of DOX was improved. The PK model, which was developed using the microdialysis method, revealed the intratumoral pH-responsive DOX distribution profiles. This facilitated the estimation of intratumoral PK parameters. The PD model with threshold concentrations contributed to the estimation of PD parameters in the three formulations, with consistent mechanisms observed. We believe that our PK-PD model can objectively assess the contributions of pH-responsive release ability and pilot molecule targetability to pharmacological effects.

Influence of an Intrinsically Disordered Region on Protein Domains Revealed by NMR-Based Electrostatic Potential Measurements

Many human proteins possess intrinsically disordered regions containing consecutive aspartate or glutamate residues ("D/E repeats"). Approximately half of them are DNA/RNA-binding proteins. In this study, using nuclear magnetic resonance (NMR) spectroscopy, we investigated the electrostatic properties of D/E repeats and their influence on folded domains within the same protein. Local electrostatic potentials were directly measured for the HMGB1 protein, its isolated D/E repeats, and DNA-binding domains by NMR. The data provide quantitative information about the electrostatic interactions between distinct segments of HMGB1. Due to the interactions between the D/E repeats and the DNA-binding domains, local electrostatic potentials of the DNA-binding domains within the full-length HMGB1 protein were largely negative despite the presence of many positively charged residues. Our NMR data on counterions and electrostatic potentials show that the D/E repeats and DNA have similar electrostatic properties and compete for the DNA-binding domains. The competition promotes dissociation of the protein-DNA complex and influences the molecular behavior of the HMGB1 protein. These effects may be general among the DNA/RNA-binding proteins with D/E repeats.

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.

Structural basis for transcription factor ZBTB7A recognition of DNA and effects of ZBTB7A somatic mutations that occur in human acute myeloid leukemia

ZBTB7A belongs to a small family of transcription factors having three members in humans (7A, 7B, and 7C). They share a BTB/POZ protein interaction domain at the amino end and a zinc-finger DNA-binding domain at the carboxyl end. They control the transcription of a wide range of genes, having varied functions in hematopoiesis, oncogenesis, and metabolism (in particular glycolysis). ZBTB7A-binding profiles at gene promoters contain a consensus G(a/c)CCC motif, followed by a CCCC sequence in some instances. Structural and mutational investigations suggest that DNA-specific contacts with the four-finger tandem array of ZBTB7A are formed sequentially, initiated from ZF1-ZF2 binding to G(a/c)CCC before spreading to ZF3-ZF4, which bind the DNA backbone and the 3' CCCC sequence, respectively. Here, we studied some mutations found in t(8;21)-positive acute myeloid leukemia patients that occur within the ZBTB7A DNA-binding domain. We determined that these mutations generally impair ZBTB7A DNA binding, with the most severe disruptions resulting from mutations in ZF1 and ZF2, and the least from a frameshift mutation in ZF3 that results in partial mislocalization. Information provided here on ZBTB7A-DNA interactions is likely applicable to ZBTB7B/C, which have overlapping functions with ZBTB7A in controlling primary metabolism.

Understanding protein diffusion on force-induced stretched DNA conformation

DNA morphology is subjected to environmental conditions and is closely coupled with its function. For example, DNA experiences stretching forces during several biological processes, including transcription and genome transactions, that significantly alter its conformation from that of B-DNA. Indeed, a well-defined 1.5 times extended conformation of dsDNA, known as Σ-DNA, has been reported in DNA complexes with proteins such as Rad51 and RecA. A striking feature in Σ-DNA is that the nucleobases are partitioned into triplets of three locally stacked bases separated by an empty rise gap of ∼ 5 Å. The functional role of such a DNA base triplet was hypothesized to be coupled with the ease of recognition of DNA bases by DNA-binding proteins (DBPs) and the physical origin of three letters (codon/anti-codon) in the genetic code. However, the underlying mechanism of base-triplet formation and the ease of DNA base-pair recognition by DBPs remain elusive. To investigate, here, we study the diffusion of a protein on a force-induced stretched DNA using coarse-grained molecular dynamics simulations. Upon pulling at the 3' end of DNA by constant forces, DNA exhibits a conformational transition from B-DNA to a ladder-like S-DNA conformation via Σ-DNA intermediate. The resulting stretched DNA conformations exhibit non-uniform base-pair clusters such as doublets, triplets, and quadruplets, of which triplets are energetically more stable than others. We find that protein favors the triplet formation compared to its unbound form while interacting non-specifically along DNA, and the relative population of it governs the ruggedness of the protein-DNA binding energy landscape and enhances the efficiency of DNA base recognition. Furthermore, we analyze the translocation mechanism of a DBP under different force regimes and underscore the significance of triplet formation in regulating the facilitated diffusion of protein on DNA. Our study, thus, provides a plausible framework for understanding the structure-function relationship between triplet formation and base recognition by a DBP and helps to understand gene regulation in complex regulatory processes.

Principles of Amino Acid and Nucleotide Revealed by Binding Affinities between Homogeneous Oligopeptides and Single-stranded DNA Molecule s

We have determined the binding strengths between nucleotides of adenine, thymine, guanine and cytosine in homogeneous single stranded DNAs and homo-octapeptides consisting of 20 common amino acids. We use a bead-based fluorescence assay for these measurements in which octapeptides are immobilized on the bead surface and ssDNAs are in solutions. The results provide a molecular basis for analyzing selectivity, specificity and polymorphisms of amino-acid-nucleotide interactions. Comparative analyses of the distribution of the binding energies reveal unique binding strengths patterns assignable to each pair of DNA nucleotide and amino acid originating from the chemical structures. Pronounced favorable (such as Arg-G , etc.) and unfavorable (such as Ile-T , etc.) binding interactions can be identified in selected groups of amino acid and nucleotide pairs that could provide basis to elucidate energetics of amino-acid-nucleotide interactions. Such interaction selectivity, specificity and polymorphism manifest the contributions from DNA backbone, DNA bases, as well as main chain and side chain of the amino acids.

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.

Heterogeneous Dynamical Environment at the Interface of a Protein-DNA Complex

Binding between protein and DNA is an essential process to regulate different biological activities. Two puzzling questions in protein-DNA recognition are (i) how the protein's binding domain identifies the DNA sequence in an aqueous solution and (ii) how the formation of the complex alters the dynamical environment around it. In this work, we present results obtained from molecular dynamics simulations of the N-terminal α-helical domain of the λ-repressor protein (in dimeric form) bound to the corresponding operator DNA. Effects of formation of the complex in modifying the microscopic dynamics of water as well as the kinetics of hydrogen bonds at the interface have been explored. Locally heterogeneous restricted water motions at the complex interface have been observed, the extent of restriction being more significant around the directly bound residues of the protein and the DNA. In particular, the calculation revealed the existence of significantly constrained motionally restricted water layer that can form either bridges around the directly bound residues of the protein and DNA or are engaged in forming water-mediated contacts between a fraction of the unbound residues. More importantly, it is observed that the restricted water motion around the complex is correlated with the hydrogen bond relaxation time scale at the interface. It is further demonstrated that the kinetics of water-water hydrogen bonds involving the bridged water are influenced more due to complex formation.

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

Interactions between DNA and DNA-binding proteins play an important role in many essential cellular processes. A key function of the DNA-binding protein p53 is to search for and bind to target sites incorporated in genomic DNA, which triggers transcriptional regulation. How do p53 molecules achieve "rapid" and "accurate" target search in living cells? The search dynamics of p53 were expected to include 3D diffusion in solution, 1D diffusion along DNA, and intersegmental transfer between two different DNA strands. Single-molecule fluorescence microscopy enabled the tracking of p53 molecules on DNA and the characterization of these dynamics quantitatively. Recent intensive single-molecule studies of p53 succeeded in revealing each of these search dynamics. Here, we review these studies and discuss the target search mechanisms of p53.

Structural Basis of Enhanced Facilitated Diffusion of DNA-Binding Protein in Crowded Cellular Milieu

Although the fast association between DNA-binding proteins (DBPs) and DNA is explained by a facilitated diffusion mechanism, in which DBPs adopt a weighted combination of three-dimensional diffusion and one-dimensional (1D) sliding and hopping modes of transportation, the role of cellular environment that contains many nonspecifically interacting proteins and other biomolecules is mostly overlooked. By performing large-scale computational simulations with an appropriately tuned model of protein and DNA in the presence of nonspecifically interacting bulk and DNA-bound crowders (genomic crowders), we demonstrate the structural basis of the enhanced facilitated diffusion of DBPs inside a crowded cellular milieu through, to our knowledge, novel 1D scanning mechanisms. In this one-dimensional scanning mode, the protein can float along the DNA under the influence of nonspecific interactions of bulk crowder molecules. The search mode is distinctly different compared to usual 1D sliding and hopping dynamics in which protein diffusion is regulated by the DNA electrostatics. In contrast, the presence of genomic crowders expedites the target search process by transporting the protein over DNA segments through the formation of a transient protein-crowder bridged complex. By analyzing the ruggedness of the associated potential energy landscape, we underpin the molecular origin of the kinetic advantages of these search modes and show that they successfully explain the experimentally observed acceleration of facilitated diffusion of DBPs by molecular crowding agents and crowder-concentration-dependent enzymatic activity of transcription factors. Our findings provide crucial insights into gene regulation kinetics inside the crowded cellular milieu.

How do DNA-bound proteins leave their binding sites? The role of facilitated dissociation

Dissociation of a protein from DNA is often assumed to be described by an off rate that is independent of other molecules in solution. Recent experiments and computational analyses have challenged this view by showing that unbinding rates (residence times) of DNA-bound proteins can depend on concentrations of nearby molecules that are competing for binding. This 'facilitated dissociation' (FD) process can occur at the single-binding site level via formation of a ternary complex, and can dominate over 'spontaneous dissociation' at low (submicromolar) concentrations. In the crowded intracellular environment FD introduces new regulatory possibilities at the level of individual biomolecule interactions.

Rational design of a DNA sequence-specific modular protein tag by tuning the alkylation kinetics

Sequence-selective chemical modification of DNA by synthetic ligands has been a long-standing challenge in the field of chemistry. Even when the ligand consists of a sequence-specific DNA binding domain and reactive group, sequence-selective reactions by these ligands are often accompanied by off-target reactions. A basic principle to design DNA modifiers that react at specific sites exclusively governed by DNA sequence recognition remains to be established. We have previously reported selective DNA modification by a self-ligating protein tag conjugated with a DNA-binding domain, termed as a modular adaptor, and orthogonal application of modular adaptors by relying on the chemoselectivity of the protein tag. The sequence-specific crosslinking reaction by the modular adaptor is thought to proceed in two steps: the first step involves the formation of a DNA-protein complex, while in the second step, a proximity-driven intermolecular crosslinking occurs. According to this scheme, the specific crosslinking reaction of a modular adaptor would be driven by the DNA recognition process only when the dissociation rate of the DNA complex is much higher than the rate constant for the alkylation reaction. In this study, as a proof of principle, a set of combinations for modular adaptors and their substrates were utilized to evaluate the reactions. Three types of modular adaptors consisting of a single type of self-ligating tag and three types of DNA binding proteins fulfill the kinetic requirements for the reaction of the self-ligating tag with a substrate and the dissociation of the DNA-protein complex. These modular adaptors actually undergo sequence-specific crosslinking reactions exclusively driven by the recognition of a specific DNA sequence. The design principle of sequence-specific modular adaptors based on the kinetic aspects of complex formation and chemical modification is applicable for developing recognition-driven selective modifiers for proteins and other biological macromolecules.

Flexibility of the Binding Regions of a Protein-DNA Complex and the Structure and Ordering of Interfacial Water

Noncovalent interactions between protein and DNA are important to comprehend different biological activities in living organisms. One important issue is how the protein identifies the target DNA and the influence of the resulting protein-DNA complex on the hydration environment around it. In this study, we have carried out atomistic molecular dynamics simulations of the protein-DNA complex formed by the dimeric form of the α-helical N-terminal domain of the λ-repressor protein with its operator DNA. Local heterogeneous flexibilities of the residues of the protein and the DNA components that are involved in binding and the microscopic structure and ordering of water around those have been investigated in detail. The calculations revealed concurrent existence of highly ordered as well as disordered water molecules at the interface. It is found that a fraction of doubly coordinated water molecules exhibit high degree of ordering at the interface, while the randomly oriented ones are coordinated with three water molecules. The effect has been found to be more around the protein and DNA residues that are in contact in the complexed state. We believe that such highly ordered two-coordinated water molecules are likely to act as an adhesive to facilitate the formation of a protein-DNA complex and maintain its structural stability.

A Plasmonic Approach to Study Protein Interaction Kinetics through the Dimerization of Functionalized Ag Nanoparticles

Understanding the kinetics of protein interactions plays a key role in biology with significant implications for the design of analytical methods for disease monitoring and diagnosis in medical care, research and industrial applications. Herein, we introduce a novel plasmonic approach to study the binding kinetics of protein-ligand interactions following the formation of silver nanoparticles (Ag NPs) dimers by UV-Vis spectroscopy that can be used as probes for antigen detection and quantification. To illustrate and test the method, the kinetics of the prototype biotin-streptavidin (Biot-STV) pair interaction was studied. Controlled aggregates (dimers) of STV functionalized Ag NPs were produced by adding stoichiometric quantities of gliadin-specific biotinylated antibodies (IgG-Biot). The dimerization kinetics was studied in a systematic way as a function of Ag NPs size and at different concentrations of IgG-Biot. The kinetics data have shown to be consistent with a complex reaction mechanism in which only the Ag NPs attached to the IgG-Biot located in a specific STV site are able to form dimers. These results help in elucidating a complex reaction mechanism involved in the dimerization kinetics of functionalized Ag NPs, which can serve as probes in surface plasmon resonance-based bioassays for the detection and quantification of different biomarkers or analytes of interest.

Intrinsically disordered domain of tumor suppressor p53 facilitates target search by ultrafast transfer between different DNA strands

Intersegmental transfer (IST) is an important strategy in the target search used by sequence-specific DNA-binding proteins (DBPs), enabling DBPs to search for targets between multiple DNA strands without dissociation. We examined the IST of the tumor suppressor p53 using ensemble stopped-flow and single-molecule fluorescence measurements. The ensemble measurements demonstrated that p53 exhibits very fast IST, whose rate constant was ∼108 M-1 s-1. To determine the domains of p53 responsible for IST, two mutants with deletions of one of its two DNA binding domains were generated. The mutant lacking the disordered C-terminal (CT) domain (the CoreTet mutant) abolished IST, whereas the mutant lacking the structured core domain (the TetCT mutant) maintained IST, clearly demonstrating the importance of the CT domain. Single-molecule fluorescence measurements further demonstrated the transfer of p53 between two tethered DNA strands. The pseudo-wild-type p53 and the TetCT mutant showed significant transfer efficiencies, whereas the transfer efficiency for the CoreTet mutant was zero. These results suggest that ultrafast IST might be promoted by four copies of the CT domain, by binding to two DNA strands simultaneously. Such ultrafast IST might be important to avoid nearby-bound DBPs during the target search process of p53 in nucleus.

Integrated rate laws for processive and distributive enzymatic turnover

Recently derived steady-state differential rate laws for the catalytic turnover of molecules containing two substrate sites are reformulated as integrated rate laws. The analysis applies to a broad class of Markovian dynamic models, motivated by the varied and often complex mechanisms associated with DNA modifying enzymes. Analysis of experimental data for the methylation kinetics of DNA by Dam (DNA adenine methyltransferase) is drastically improved through the use of integrated rate laws. Data that are too noisy for fitting to differential predictions are reliably interpreted through the integrated rate laws.

Facilitated Diffusion Mechanisms in DNA Base Excision Repair and Transcriptional Activation

Preservation of the coding potential of the genome and highly regulated gene expression over the life span of a human are two fundamental requirements of life. These processes require the action of repair enzymes or transcription factors that efficiently recognize specific sites of DNA damage or transcriptional regulation within a restricted time frame of the cell cycle or metabolism. A failure of these systems to act results in accumulated mutations, metabolic dysfunction, and disease. Despite the multifactorial complexity of cellular DNA repair and transcriptional regulation, both processes share a fundamental physical requirement that the proteins must rapidly diffuse to their specific DNA-binding sites that are embedded within the context of a vastly greater number of nonspecific DNA-binding sites. Superimposed on the needle-in-the-haystack problem is the complex nature of the cellular environment, which contains such high concentrations of macromolecules that the time frame for diffusion is expected to be severely extended as compared to dilute solution. Here we critically review the mechanisms for how these proteins solve the needle-in-the-haystack problem and how the effects of cellular macromolecular crowding can enhance facilitated diffusion processes. We restrict the review to human proteins that use stochastic, thermally driven site-recognition mechanisms, and we specifically exclude systems involving energy cofactors or circular DNA clamps. Our scope includes ensemble and single-molecule studies of the past decade or so, with an emphasis on connecting experimental observations to biological function.

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 Macromolecular Crowding on the Intracellular Diffusion of DNA Binding Proteins

Recent experiments suggest that cellular crowding facilitates the target search dynamics of proteins on DNA, the mechanism of which is not yet known. By using large scale computer simulations, we show that two competing factors, namely the width of the depletion layer that separates the crowder cloud from the DNA molecule and the degree of protein-crowder crosstalk, act in harmony to affect the target search dynamics of proteins. The impacts vary from nonspecific to specific target search regime. During a nonspecific search, dynamics of a protein is only minimally affected, whereas, a significantly different behaviour is observed when the protein starts forming a specific protein-DNA complex. We also find that the severity of impacts largely depends upon physiological crowder concentration and deviation from it leads to attenuation in the binding kinetics. Based on extensive kinetic study and binding energy landscape analysis, we further present a comprehensive molecular description of the search process that allows us to interpret the experimental findings.

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.

Transcription factor clusters regulate genes in eukaryotic cells

Transcription is regulated through binding factors to gene promoters to activate or repress expression, however, the mechanisms by which factors find targets remain unclear. Using single-molecule fluorescence microscopy, we determined in vivo stoichiometry and spatiotemporal dynamics of a GFP tagged repressor, Mig1, from a paradigm signaling pathway of Saccharomyces cerevisiae. We find the repressor operates in clusters, which upon extracellular signal detection, translocate from the cytoplasm, bind to nuclear targets and turnover. Simulations of Mig1 configuration within a 3D yeast genome model combined with a promoter-specific, fluorescent translation reporter confirmed clusters are the functional unit of gene regulation. In vitro and structural analysis on reconstituted Mig1 suggests that clusters are stabilized by depletion forces between intrinsically disordered sequences. We observed similar clusters of a co-regulatory activator from a different pathway, supporting a generalized cluster model for transcription factors that reduces promoter search times through intersegment transfer while stabilizing gene expression.

Microscopic understanding of the conformational features of a protein–DNA complex

Protein–DNA interactions play crucial roles in different stages of genetic activities, such as replication of genome, initiation of transcription,etc.

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.

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.

Time-Resolved Detection of Light-Induced Dimerization of Monomeric Aureochrome-1 and Change in Affinity for DNA

Aureochrome (Aureo) is a recently discovered blue light sensor protein initially from Vaucheria frigida, in which it controls blue light-dependent branch formation and/or development of a sex organ by a light-dependent change in the affinity for DNA. Although photochemical reactions of Aureo-LOV (LOV is a C-terminal light-oxygen-voltage domain) and the N-terminal truncated construct containing a bZIP (N-terminal basic leucine zipper domain) and a LOV domain have previously been reported, the reaction kinetics of the change in affinity for DNA have never been elucidated. The reactions of Aureo where the cysteines are replaced by serines (AureoCS) as well as the kinetics of the change in affinity for a target DNA are investigated in the time-domain. The dimerization rate constant is obtained as 2.8 × 10(4) M(-1) s(-1), which suggests that the photoinduced dimerization occurs in the LOV domain and the bZIP domain dimerizes using the interaction with DNA. Surprisingly, binding with the target DNA is completed very quickly, 7.7 × 10(4) M(-1) s(-1), which is faster than the protein dimerization rate. It is proposed that the nonspecific electrostatic interaction, which is observed as a weak binding with DNA, may play a role in the efficient searching for the target sequence within the DNA.

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.

Single-molecule fluorescence studies on DNA looping

Structure and dynamics of DNA impact how the genetic code is processed and maintained. In addition to its biological importance, DNA has been utilized as building blocks of various nanomachines and nanostructures. Thus, understanding the physical properties of DNA is of fundamental importance to basic sciences and engineering applications. DNA can undergo various physical changes. Among them, DNA looping is unique in that it can bring two distal sites together, and thus can be used to mediate interactions over long distances. In this paper, we introduce a FRET-based experimental tool to study DNA looping at the single molecule level. We explain the connection between experimental measurables and a theoretical concept known as the J factor with the intent of raising awareness of subtle theoretical details that should be considered when drawing conclusions. We also explore DNA looping-assisted protein diffusion mechanism called intersegmental transfer using protein induced fluorescence enhancement (PIFE). We present some preliminary results and future outlooks.

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.

Searching target sites on DNA by proteins: Role of DNA dynamics under confinement

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

Physicochemical Properties of Ion Pairs of Biological Macromolecules

Ion pairs (also known as salt bridges) of electrostatically interacting cationic and anionic moieties are important for proteins and nucleic acids to perform their function. Although numerous three-dimensional structures show ion pairs at functionally important sites of biological macromolecules and their complexes, the physicochemical properties of the ion pairs are not well understood. Crystal structures typically show a single state for each ion pair. However, recent studies have revealed the dynamic nature of the ion pairs of the biological macromolecules. Biomolecular ion pairs undergo dynamic transitions between distinct states in which the charged moieties are either in direct contact or separated by water. This dynamic behavior is reasonable in light of the fundamental concepts that were established for small ions over the last century. In this review, we introduce the physicochemical concepts relevant to the ion pairs and provide an overview of the recent advancement in biophysical research on the ion pairs of biological macromolecules.

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.

DNA-Segment-Facilitated Dissociation of Fis and NHP6A from DNA Detected via Single-Molecule Mechanical Response

The rate of dissociation of a DNA-protein complex is often considered to be a property of that complex, without dependence on other nearby molecules in solution. We study the kinetics of dissociation of the abundant Escherichia coli nucleoid protein Fis from DNA, using a single-molecule mechanics assay. The rate of Fis dissociation from DNA is strongly dependent on the solution concentration of DNA. The off-rate (k(off)) of Fis from DNA shows an initially linear dependence on solution DNA concentration, characterized by an exchange rate of k(ex)≈9×10(-4) (ng/μl)(-1) s(-1) for 100 mM univalent salt buffer, with a very small off-rate at zero DNA concentration. The off-rate saturates at approximately k(off,max)≈8×10(-3) s(-1) for DNA concentrations above ≈20 ng/μl. This exchange reaction depends mainly on DNA concentration with little dependence on the length of the DNA molecules in solution or on binding affinity, but this does increase with increasing salt concentration. We also show data for the yeast HMGB protein NHP6A showing a similar DNA-concentration-dependent dissociation effect, with faster rates suggesting generally weaker DNA binding by NHP6A relative to Fis. Our results are well described by a model with an intermediate partially dissociated state where the protein is susceptible to being captured by a second DNA segment, in the manner of "direct transfer" reactions studied for other DNA-binding proteins. This type of dissociation pathway may be important to protein-DNA binding kinetics in vivo where DNA concentrations are large.

Real sequence effects on the search dynamics of transcription factors on DNA

Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.

Biomolecular Interaction Analysis Using an Optical Surface Plasmon Resonance Biosensor: The Marquardt Algorithm vs Newton Iteration Algorithm

Kinetic analysis of biomolecular interactions are powerfully used to quantify the binding kinetic constants for the determination of a complex formed or dissociated within a given time span. Surface plasmon resonance biosensors provide an essential approach in the analysis of the biomolecular interactions including the interaction process of antigen-antibody and receptors-ligand. The binding affinity of the antibody to the antigen (or the receptor to the ligand) reflects the biological activities of the control antibodies (or receptors) and the corresponding immune signal responses in the pathologic process. Moreover, both the association rate and dissociation rate of the receptor to ligand are the substantial parameters for the study of signal transmission between cells. A number of experimental data may lead to complicated real-time curves that do not fit well to the kinetic model. This paper presented an analysis approach of biomolecular interactions established by utilizing the Marquardt algorithm. This algorithm was intensively considered to implement in the homemade bioanalyzer to perform the nonlinear curve-fitting of the association and disassociation process of the receptor to ligand. Compared with the results from the Newton iteration algorithm, it shows that the Marquardt algorithm does not only reduce the dependence of the initial value to avoid the divergence but also can greatly reduce the iterative regression times. The association and dissociation rate constants, k_a, k_d and the affinity parameters for the biomolecular interaction, K_A, K_D, were experimentally obtained 6.969×10⁵ mL·g^-1·s^-1, 0.00073 s^-1, 9.5466×10⁸ mL·g^-1 and 1.0475×10^-9 g·mL^-1, respectively from the injection of the HBsAg solution with the concentration of 16ng·mL^-1. The kinetic constants were evaluated distinctly by using the obtained data from the curve-fitting results.

Structural impact of complete CpG methylation within target DNA on specific complex formation of the inducible transcription factor Egr-1

The inducible transcription factor Egr-1 binds specifically to 9-bp target sequences containing two CpG sites that can potentially be methylated at four cytosine bases. Although it appears that complete CpG methylation would make an unfavorable steric clash in the previous crystal structures of the complexes with unmethylated or partially methylated DNA, our affinity data suggest that DNA recognition by Egr-1 is insensitive to CpG methylation. We have determined, at a 1.4-Å resolution, the crystal structure of the Egr-1 zinc-finger complex with completely methylated target DNA. Structural comparison of the three different methylation states reveals why Egr-1 can recognize the target sequences regardless of CpG methylation.

Mechanisms of diffusional search for specific targets by DNA-dependent proteins

To perform their functions, many DNA-dependent proteins have to quickly locate specific targets against the vast excess of nonspecific DNA. Although this problem was first formulated over 40 years ago, the mechanism of such search remains one of the unsolved fundamental problems in the field of protein-DNA interactions. Several complementary mechanisms have been suggested: sliding, based on one-dimensional random diffusion along the DNA contour; hopping, in which the protein "jumps" between the closely located DNA fragments; macroscopic association-dissociation of the protein-DNA complex; and intersegmental transfer. This review covers the modern state of the problem of target DNA search, theoretical descriptions, and methods of research at the macroscopic (molecule ensembles) and microscopic (individual molecules) levels. Almost all studied DNA-dependent proteins search for specific targets by combined three-dimensional diffusion and one-dimensional diffusion along the DNA contour.

Search by proteins for their DNA target site: 2. The effect of DNA conformation on the dynamics of multidomain proteins

Multidomain transcription factors, which are especially abundant in eukaryotic genomes, are advantageous to accelerate the search kinetics for target site because they can follow the intersegment transfer via the monkey-bar mechanism in which the protein forms a bridged intermediate between two distant DNA regions. Monkey-bar dynamics highly depends on the properties of the multidomain protein (the affinity of each of the constituent domains to the DNA and the length of the linker) and the DNA molecules (their inter-distance and inter-angle). In this study, we investigate using coarse-grained molecular dynamics simulations how the local conformation of the DNA may affect the DNA search performed by a multidomain protein Pax6 in comparison to that of the isolated domains. Our results suggest that in addition to the common rotation-coupled translation along the DNA major groove, for curved DNA the tethered domains may slide in a rotation-decoupled sliding mode. Furthermore, the multidomain proteins move by longer jumps on curved DNA compared with those performed by the single domain protein. The long jumps originate from the DNA curvature bringing two sequentially distant DNA sites into close proximity with each other and they suggest that multidomain proteins may move on highly curved DNA faster than linear DNA.

Distinct facilitated diffusion mechanisms by E. coli Type II restriction endonucleases

The passive search by proteins for particular DNA sequences involving nonspecific DNA is essential for gene regulation, DNA repair, phage defense, and diverse epigenetic processes. Distinct mechanisms contribute to these searches, and it remains unresolved as to which mechanism or blend of mechanisms best suits a particular protein and, more importantly, its biological role. To address this, we compare the translocation properties of two well-studied bacterial restriction endonucleases (ENases), EcoRI and EcoRV. These dimeric, magnesium-dependent enzymes hydrolyze related sites (EcoRI ENase, 5'-GAATTC-3'; EcoRV ENase, 5'-GATATC-3'), leaving overhangs and blunt DNA segments, respectively. Here, we demonstrate that the extensive sliding by EcoRI ENase, involving sliding up to ∼600 bp prior to dissociating from the DNA, contrasts with a larger reliance on hopping mechanism(s) by EcoRV ENase. The mechanism displayed by EcoRI ENase results in a highly thorough search of DNA, whereas the EcoRV ENase mechanism results in an extended, yet less rigorous, interrogation of DNA sequence space. We describe how these mechanistic distinctions are complemented by other aspects of these endonucleases, such as the 10-fold higher in vivo concentrations of EcoRI ENase compared to that of EcoRV ENase. Further, we hypothesize that the highly diverse enzyme arsenal that bacteria employ against foreign DNA involves seemingly similar enzymes that rely on distinct but complementary search mechanisms. Our comparative approach reveals how different proteins utilize distinct site-locating strategies.

An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding

We present a computational model of transcription factor motion that explains both the observed rapid target finding of transcription factors, and how this motion influences protein and genome structure. Using the Smoldyn software, we modelled transcription factor motion arising from a combination of unrestricted 3D diffusion in the nucleoplasm, sliding along the DNA filament, and transferring directly between filament sections by intersegmental transfer. This presents a fine-grain picture of the way in which transcription factors find their targets two orders of magnitude faster than 3D diffusion alone allows. Eukaryotic genomes contain sections of nucleosome free regions (NFRs) around the promoters; our model shows that the presence and size of these NFRs can be explained as their acting as antennas on which transcription factors slide to reach their targets. Additionally, our model shows that intersegmental transfer may have shaped the quaternary structure of transcription factors: sequence specific DNA binding proteins are unusually enriched in dimers and tetramers, perhaps because these allow intersegmental transfer, which accelerates target site finding. Finally, our model shows that a ‘hopping’ motion can emerge from 3D diffusion on small scales. This explains the apparently long sliding lengths that have been observed for some DNA binding proteins observed in vitro. Together, these results suggest that transcription factor diffusion dynamics help drive the evolution of protein and genome structure.

Effective strategy to assign ¹H- ¹⁵N heteronuclear correlation NMR signals from lysine side-chain NH3₃⁺ groups of proteins at low temperature

Recent studies have shown that lysine side-chain NH3(+) groups are excellent probes for NMR investigations of dynamics involving hydrogen bonds and ion pairs relevant to protein function. However, due to rapid hydrogen exchange, observation of (1)H-(15)N NMR cross peaks from lysine NH3(+) groups often requires use of a relatively low temperature, which renders difficulty in resonance assignment. Here we present an effective strategy to assign (1)H and (15)N resonances of NH3(+) groups at low temperatures. This strategy involves two new (1)H/(13)C/(15)N triple-resonance experiments for lysine side chains. Application to a protein-DNA complex is demonstrated.

Positive and negative impacts of nonspecific sites during target location by a sequence-specific DNA-binding protein: origin of the optimal search at physiological ionic strength

The inducible transcription factor Egr-1, which recognizes a 9-bp target DNA sequence via three zinc-finger domains, rapidly activates particular genes upon cellular stimuli such as neuronal signals and vascular stresses. Here, using the stopped-flow fluorescence method, we measured the target search kinetics of the Egr-1 zinc-finger protein at various ionic strengths between 40 and 400 mM KCl and found the most efficient search at 150 mM KCl. We further investigated the kinetics of intersegment transfer, dissociation, and sliding of this protein on DNA at distinct concentrations of KCl. Our data suggest that Egr-1's kinetic properties are well suited for efficient scanning of chromosomal DNA in vivo. Based on a newly developed theory, we analyzed the origin of the optimal search efficiency at physiological ionic strength. Target association is accelerated by nonspecific binding to nearby sites and subsequent sliding to the target as well as by intersegment transfer. Although these effects are stronger at lower ionic strengths, such conditions also favor trapping of the protein at distant nonspecific sites, decelerating the target association. Our data demonstrate that Egr-1 achieves the optimal search at physiological ionic strength through a compromise between the positive and negative impacts of nonspecific interactions with DNA.