Bending of DNA by transcription factors

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

COL1A1 proximal promoter topology regulates its transcriptional response to transforming growth factor β

OBJECTIVE

The COL1A1 proximal promoter contains two GC-rich regions and two inverted CCAAT boxes. The transcription factors Sp1 and CBF bind to the GC sequence at -122 to -115 bp and the inverted CCAAT box at -101 to -96 bp, respectively, and stimulate COL1A1 transcriptional activity.

METHODS

To further define the regulatory mechanisms controlling COL1A1 expression by Sp1 and CBF, we introduced 2, 4, 6, or 8 thymidine nucleotides (T-tracts) at position -111 bp of the COL1A1 gene promoter to increase the physical distance between these two binding sites and examined in vitro the transcriptional activities of the resulting constructs and their response to TGF-β1.`.

RESULTS

Insertion of 2 or 4 nucleotides decreased COL1A1 promoter activity by up to 70%. Furthermore, the expected increase in COL1A1 transcription in response to TGF-β1 was abolished. Computer modeling of the modified DNA structure indicated that increasing the physical distance between the Sp1 and CBF binding sites introduces a rotational change in the DNA topology that disrupts the alignment of Sp1 and CBF binding sites and likely alters protein-protein interactions among these transcription factors or their associated co-activators.

CONCLUSION

The topology of the COL1A1 proximal promoter is crucial in determining the transcriptional activity of the gene and its response to the stimulatory effects of TGF-β1.

Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research

The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.

-

The interaction of nucleic acids with proteins plays an important role in many fundamental biological processes in living cells, including replication, transcription, and translation. Therefore, understanding nucleic acid-protein interaction is of high relevance in many areas of biology, medicine and technology. During almost four decades of its existence atomic force microscopy (AFM) accumulated a significant experience in investigation of biological molecules at a single-molecule level. AFM has become a powerful tool of molecular biology and biophysics providing unique information about properties, structure, and functioning of biomolecules. Despite a great variety of nucleic acid-protein systems under AFM investigations, there are a number of typical approaches for such studies. This review is devoted to the analysis of the typical AFM-based approaches of investigation of DNA (RNA)-protein complexes with a major focus on transcription studies. The basic strategies of AFM analysis of nucleic acid-protein complexes including investigation of the products of DNA-protein reactions and real-time dynamics of DNA-protein interaction are categorized and described by the example of the most relevant research studies. The described approaches and protocols have many universal features and, therefore, are applicable for future AFM studies of various nucleic acid-protein systems.

Structure of a 28.5 kDa duplex-embedded G-quadruplex system resolved to 7.4 Å resolution with cryo-EM

Genomic regions with high guanine content can fold into non-B form DNA four-stranded structures known as G-quadruplexes (G4s). Extensive in vivo investigations have revealed that promoter G4s are transcriptional regulators. Little structural information exists for these G4s embedded within duplexes, their presumed genomic environment. Here, we report the 7.4 Å resolution structure and dynamics of a 28.5 kDa duplex-G4-duplex (DGD) model system using cryo-EM, molecular dynamics, and small-angle X-ray scattering (SAXS) studies. The DGD cryo-EM refined model features a 53° bend induced by a stacked duplex-G4 interaction at the 5' G-tetrad interface with a persistently unstacked 3' duplex. The surrogate complement poly dT loop preferably stacks onto the 3' G-tetrad interface resulting in occlusion of both 5' and 3' tetrad interfaces. Structural analysis shows that the DGD model is quantifiably more druggable than the monomeric G4 structure alone and represents a new structural drug target. Our results illustrate how the integration of cryo-EM, MD, and SAXS can reveal complementary detailed static and dynamic structural information on DNA G4 systems.

A DNA origami-based device for investigating DNA bending proteins by transmission electron microscopy

The DNA origami technique offers precise positioning of nanoscale objects with high accuracy. This has facilitated the development of DNA origami-based functional nanomechanical devices that enable the investigation of DNA-protein interactions at the single particle level. Herein, we used the DNA origami technique to fabricate a nanoscale device for studying DNA bending proteins. For a proof of concept, we used TATA-box binding protein (TBP) to evaluate our approach. Upon binding to the TATA box, TBP causes a bend to DNA of ∼90°. Our device translates this bending into an angular change that is readily observable with a conventional transmission electron microscope (TEM). Furthermore, we investigated the roles of transcription factor II A (TF(II)A) and transcription factor II B (TF(II)B). Our results indicate that TF(II)A introduces additional bending, whereas TF(II)B does not significantly alter the TBP-DNA structure. Our approach can be readily adopted to a wide range of DNA-bending proteins and will aid the development of DNA-origami-based devices tailored for the investigation of DNA-protein interactions.

Single-Molecule Methods for Investigating the Double-Stranded DNA Bendability

The various DNA-protein interactions associated with the expression of genetic information involve double-stranded DNA (dsDNA) bending. Due to the importance of the formation of the dsDNA bending structure, dsDNA bending properties have long been investigated in the biophysics field. Conventionally, DNA bendability is characterized by innate averaging data from bulk experiments. The advent of single-molecule methods, such as atomic force microscopy, optical and magnetic tweezers, tethered particle motion, and single-molecule fluorescence resonance energy transfer measurement, has provided valuable tools to investigate not only the static structures but also the dynamic properties of bent dsDNA. Here, we reviewed the single-molecule methods that have been used for investigating dsDNA bendability and new findings related to dsDNA bending. Single-molecule approaches are promising tools for revealing the unknown properties of dsDNA related to its bending, particularly in cells.

OUP accepted manuscript

DNA mechanical properties play a critical role in every aspect of DNA-dependent biological processes. Recently a high throughput assay named loop-seq has been developed to quantify the intrinsic bendability of a massive number of DNA fragments simultaneously. Using the loop-seq data, we develop a software tool, DNAcycP, based on a deep-learning approach for intrinsic DNA cyclizability prediction. We demonstrate DNAcycP predicts intrinsic DNA cyclizability with high fidelity compared to the experimental data. Using an independent dataset from in vitro selection for enrichment of loopable sequences, we further verified the predicted cyclizability score, termed C-score, can well distinguish DNA fragments with different loopability. We applied DNAcycP to multiple species and compared the C-scores with available high-resolution chemical nucleosome maps. Our analyses showed that both yeast and mouse genomes share a conserved feature of high DNA bendability spanning nucleosome dyads. Additionally, we extended our analysis to transcription factor binding sites and surprisingly found that the cyclizability is substantially elevated at CTCF binding sites in the mouse genome. We further demonstrate this distinct mechanical property is conserved across mammalian species and is inherent to CTCF binding DNA motif.

The molecular basis for the development of Adult T-cell leukemia/lymphoma in patients with an IRF4K59R mutation

Interferon regulatory factor 4 (IRF4) is an essential regulator in the development of many immune cells, including B and T cells and has been implicated directly in numerous haematological malignancies, including Adult T-cell leukemia/lymphoma (ATLL). Recently, an activating mutation in the DNA binding domain of IRF4 (IRF4^K59R ), was found as a recurrent somatic mutation in ATLL patients. However, it remains unknown how this mutation gives rise to the observed oncogenic effect. To understand the mode of IRF4^K59R mediated gain of function in ATLL pathogenesis, we have determined the structural and affinity basis of IRF4^K59R /DNA homodimer complex using X-ray crystallography and surface plasmon resonance. Our study shows that arginine substitution (R59) results in the reorientation of the side chain, enabling the guanidium group to interact with the phosphate backbone of the DNA helix. This markedly contrasts with the IRF4^WT wherein the K59 interacts exclusively with DNA bases. Further, the arginine mutation causes enhanced DNA bending, enabling the IRF4^K59R to interact more robustly with known DNA targets, as evidenced by increased binding affinity of the protein-DNA complex. Together, we demonstrate how key structural features underpin the basis for this activating mutation, thereby providing a molecular rationale for IRF4^K59R -mediated ATLL development. This article is protected by copyright. All rights reserved.

Structure-based virtual screening identifies an 8-hydroxyquinoline as a small molecule GLI1 inhibitor

The glioma-associated family of transcription factors (GLI) have emerged as a promising therapeutic target for a variety of human cancers. In particular, GLI1 plays a central role as a transcriptional regulator for multiple oncogenic signaling pathways, including the hedgehog (Hh) signaling pathway. We undertook a computational screening approach to identify small molecules that directly bind GLI1 for potential development as inhibitors of GLI-mediated transcription. Through these studies, we identified compound 1, which is an 8-hydroxyquinoline, as a high-affinity binder of GLI1. Compound 1 inhibits GLI1-mediated transcriptional activity in several Hh-dependent cellular models, including a primary model of murine medulloblastoma. We also performed a series of computational analyses to define more clearly the mechanism(s) through which 1 inhibits GLI1 function after binding. Our results strongly suggest that binding of 1 to GLI1 does not prevent GLI1/DNA binding nor disrupt the GLI1/DNA complex, but rather, it induces specific conformational changes in the overall complex that prevent proper GLI function. These results highlight the potential of this compound for further development as an anti-cancer agent that targets GLI1.

HSP27 role in cardioprotection by modulating chemotherapeutic doxorubicin-induced cell death

The common phenomenon expected from any anti-cancer drug in use is to kill the cancer cells without any side effects to non-malignant cells. Doxorubicin is an anthracycline derivative anti-cancer drug active over different types of cancers with anti-cancer activity but attributed to unintended cytotoxicity and genotoxicity triggering mitogenic signals inducing apoptosis. Administration of doxorubicin tends to both acute and chronic toxicity resulting in cardiomyopathy (left ventricular dysfunction) and congestive heart failure (CHF). Cardiotoxicity is prevented through administration of different cardioprotectants along with the drug. This review elaborates on mechanism of drug-mediated cardiotoxicity and attenuation principle by different cardioprotectants, with a focus on Hsp27 as cardioprotectant by prevention of drug-induced oxidative stress, cell survival pathways with suppression of intrinsic cell death. In conclusion, Hsp27 may offer an exciting/alternating cardioprotectant, with a wider study being need of the hour, specifically on primary cell line and animal models in conforming its cardioprotectant behaviour.

The ChiS-Family DNA-Binding Domain Contains a Cryptic Helix-Turn-Helix Variant

Sequence-specific DNA-binding domains (DBDs) are conserved in all domains of life. These proteins carry out a variety of cellular functions, and there are a number of distinct structural domains already described that allow for sequence-specific DNA binding, including the ubiquitous helix-turn-helix (HTH) domain. In the facultative pathogen Vibrio cholerae, the chitin sensor ChiS is a transcriptional regulator that is critical for the survival of this organism in its marine reservoir. We recently showed that ChiS contains a cryptic DBD in its C terminus. This domain is not homologous to any known DBD, but it is a conserved domain present in other bacterial proteins. Here, we present the crystal structure of the ChiS DBD at a resolution of 1.28 Å. We find that the ChiS DBD contains an HTH domain that is structurally similar to those found in other DNA-binding proteins, like the LacI repressor. However, one striking difference observed in the ChiS DBD is that the canonical tight turn of the HTH is replaced with an insertion containing a β-sheet, a variant which we term the helix-sheet-helix. Through systematic mutagenesis of all positively charged residues within the ChiS DBD, we show that residues within and proximal to the ChiS helix-sheet-helix are critical for DNA binding. Finally, through phylogenetic analyses we show that the ChiS DBD is found in diverse proteobacterial proteins that exhibit distinct domain architectures. Together, these results suggest that the structure described here represents the prototypical member of the ChiS-family of DBDs.IMPORTANCE Regulating gene expression is essential in all domains of life. This process is commonly facilitated by the activity of DNA-binding transcription factors. There are diverse structural domains that allow proteins to bind to specific DNA sequences. The structural basis underlying how some proteins bind to DNA, however, remains unclear. Previously, we showed that in the major human pathogen Vibrio cholerae, the transcription factor ChiS directly regulates gene expression through a cryptic DNA-binding domain. This domain lacked homology to any known DNA-binding protein. In the current study, we determined the structure of the ChiS DNA-binding domain (DBD) and found that the ChiS-family DBD is a cryptic variant of the ubiquitous helix-turn-helix (HTH) domain. We further demonstrate that this domain is conserved in diverse proteins that may represent a novel group of transcriptional regulators.

Mechanical Flexibility of DNA: A Quintessential Tool for DNA Nanotechnology

The mechanical properties of DNA have enabled it to be a structural and sensory element in many nanotechnology applications. While specific base-pairing interactions and secondary structure formation have been the most widely utilized mechanism in designing DNA nanodevices and biosensors, the intrinsic mechanical rigidity and flexibility are often overlooked. In this article, we will discuss the biochemical and biophysical origin of double-stranded DNA rigidity and how environmental and intrinsic factors such as salt, temperature, sequence, and small molecules influence it. We will then take a critical look at three areas of applications of DNA bending rigidity. First, we will discuss how DNA's bending rigidity has been utilized to create molecular springs that regulate the activities of biomolecules and cellular processes. Second, we will discuss how the nanomechanical response induced by DNA rigidity has been used to create conformational changes as sensors for molecular force, pH, metal ions, small molecules, and protein interactions. Lastly, we will discuss how DNA's rigidity enabled its application in creating DNA-based nanostructures from DNA origami to nanomachines.

Ecdysteroid-mimicking compounds act as both agonists and antagonists to the crustacean ecdysone receptor

To characterize the potential endocrine-disrupting chemicals (EDCs) in the environment that interact with the crustacean ecdysone receptor (EcR), we established a method involving in silico modeling/molecular docking and in vitro reporter gene assay. Cherry shrimp (Neocaridina davidi) EcR (NdEcR) and retinoid X receptor (NdRxR) were identified and cloned for use in this method. A theoretical 3D model of NdEcR ligand-binding domain (LBD) was built in silico based on sequence homology with the established X-ray structure of insect EcR. The interaction of the NdEcR LBD with ecdysteroids, diacylhydrazine (DAH) pesticides, and other potential EDCs was evaluated using molecular docking programs. The results revealed that the ligand-binding pocket in the NdEcR LBD was flexible and adaptive for accommodating ligands of different shapes. The agonistic and antagonistic activities of the candidate compounds were further assessed by in vitro reporter gene assay using human cell lines transiently transfected with NdEcR and NdRxR expression plasmids and a reporter plasmid containing synthesized ecdysone response element. The assay was validated by the dose-dependent responses of EcR-mediated gene transcription after treating the transfected cell lines with ecdysteroids, 20-hydroxyecdysone, and ponasterone A. Examination of the candidate compounds using the reporter gene assay revealed restricted functional specificity to ecdysteroids and DAHs. Three of the tested DAH pesticides originally targeting the insect EcR were found to be weak agonists and strong antagonists of NdEcR. These results suggest that DAHs are potential EDCs for crustaceans that disrupt their ecdysteroid signals by functioning as EcR agonists or antagonists.

Parallelized DNA tethered bead measurements to scrutinize DNA mechanical structure

Tethering beads to DNA offers a panel of single molecule techniques for the refined analysis of the conformational dynamics of DNA and the elucidation of the mechanisms of enzyme activity. Recent developments include the massive parallelization of these techniques achieved by the fabrication of dedicated nanoarrays by soft nanolithography. We focus here on two of these techniques: the Tethered Particle motion and Magnetic Tweezers allowing analysis of the behavior of individual DNA molecules in the absence of force and under the application of a force and/or a torque, respectively. We introduce the experimental protocols for the parallelization and discuss the benefits already gained, and to come, for these single molecule investigations.

Genetic and Epigenetic Mechanisms of β-Globin Gene Switching

Vertebrates have multiple forms of hemoglobin that differ in the composition of their polypeptide chains. During ontogenesis, the composition of these subunits changes. Genes encoding different α- and β-polypeptide chains are located in two multigene clusters on different chromosomes. Each cluster contains several genes that are expressed at different stages of ontogenesis. The phenomenon of stage-specific transcription of globin genes is referred to as globin gene switching. Mechanisms of expression switching, stage-specific activation, and repression of transcription of α- and β-globin genes are of interest from both theoretical and practical points of view. Alteration of balanced expression of globin genes, which usually occurs due to damage to adult β-globin genes, leads to development of severe diseases - hemoglobinopathies. In most cases, reactivation of the fetal hemoglobin gene in patients with β-thalassemia and sickle cell disease can reduce negative consequences of irreversible alterations of expression of the β-globin genes. This review focuses on the current state of research on genetic and epigenetic mechanisms underlying stage-specific switching of β-globin genes.

Chromatin looping as a target for altering erythroid gene expression

The β-hemoglobinopathies are the most common monogenic disorders in humans, with symptoms arising after birth when the fetal γ-globin genes are silenced and the adult β-globin gene is activated. There is a growing appreciation that genome organization and the folding of chromosomes are key determinants of gene transcription. Underlying this function is the activity of transcriptional enhancers that increase the transcription of target genes over long linear distances. To accomplish this, enhancers engage in close physical contact with target promoters through chromosome folding or looping that is orchestrated by protein complexes that bind to both sites and stabilize their interaction. We find that enhancer activity can be redirected with concomitant changes in gene transcription. Both targeting the β-globin locus control region (LCR) to the γ-globin gene in adult erythroid cells by tethering and epigenetic unmasking of a silenced γ-globin gene lead to increased frequency of LCR/γ-globin contacts and reduced LCR/β-globin contacts. The outcome of these manipulations is robust, pancellular γ-globin transcription activation with a concomitant reduction in β-globin transcription. These examples show that chromosome looping may be considered a therapeutic target for gene activation in β-thalassemia and sickle cell disease.

DNA bending by proteins: utilizing plasmid pBendAT as a tool

Protein-induced DNA bending plays a key role in many essential biological processes, such as DNA replication, recombination, and transcription, and can be analyzed by a variety of biochemical and biophysical methods, such as electrophoresis mobility shift assay (EMSA), X-ray crystallography, nuclear magnetic resonance (NMR), and DNA ring closure assay. In this chapter, I will provide a detailed protocol for studying protein-induced DNA bending by utilizing the plasmid pBendAT. pBendAT carries a 230 bp DNA segment containing five pairs of restriction-endonuclease recognition sites and can be used to produce a set of five DNA fragments of identical length, with each fragment having a different positioning of a protein-binding site. Binding of a protein to this site will divide a DNA fragment into two DNA segments. If protein binding leads to DNA bending, then the two DNA segments will be at an angle to each other and the distance between the DNA ends will shorten. As a result, the gel mobility of the protein-DNA complex will be affected as the mobility of a rigid DNA fragment is inversely proportional to the end-to-end distance. Therefore, by analyzing how the position of protein binding within a fragment affects the gel mobility of the protein-DNA complexes, we are able to determine the DNA bending angle and the location of the bend. The DNA fragments of identical length can also be conveniently generated by PCR amplification using pBendAT as the DNA template. Since the 230 bp DNA fragment of pBendAT does not contain more than two consecutive AT base pairs, pBendAT is particularly suitable for the assessment of DNA bending induced by proteins recognizing AT-rich DNA sequences cloned in the 230 bp DNA fragment.

Communication of genome regulatory elements in a folded chromosome

The most popular model of gene activation by remote enhancers postulates that the enhancers interact directly with target promoters via the looping of intervening DNA fragments. This interaction is thought to be necessary for the stabilization of the Pol II pre-initiation complex and/or for the transfer of transcription factors and Pol II, which are initially accumulated at the enhancer, to the promoter. The direct interaction of enhancer(s) and promoter(s) is only possible when these elements are located in close proximity within the nuclear space. Here, we discuss the molecular mechanisms for maintaining the close proximity of the remote regulatory elements of the eukaryotic genome. The models of an active chromatin hub (ACH) and an active nuclear compartment are considered, focusing on the role of chromatin folding in juxtaposing remote DNA sequences. The interconnection between the functionally dependent architecture of the interphase chromosome and nuclear compartmentalization is also discussed.

Attractive hydration forces in DNA-dendrimer interactions on the nanometer scale

The energetic contribution of attractive hydration forces arising from water ordering is an interesting but often neglected aspect of macromolecular interactions. Ordering effects of water can bring about cooperativity in many intermolecular transactions, in both the short and long range. Given its high charge density, this is of particular importance for DNA. For instance, in nanotechnology, highly charged dendrimers are used for DNA compaction and transfection. Hypothesizing that water ordering and hydration forces should be maximal for DNA complexes that show charge complementarity (positive-negative), we present here an analysis of water ordering from molecular dynamics simulations and free energy calculations of the interaction between DNA and a nanoparticle with a high positive charge density. Our results indicate not only that complexation of the dendrimer with DNA affects the local water structure but also that ordered water molecules facilitate long-range interactions between the molecules. This contributes significantly to the free energy of binding of dendrimers to DNA and extends the interaction well beyond the electrostatic range of the DNA. Such water effects are of potentially substantial importance in cases when molecules appear to recognize each other across sizable distances, or for which kinetic rates are too fast to be due to pure diffusion. Our results are in good agreement with experiments on the role of solvent in DNA condensation by multivalent cations and exemplify a microscopic realization of mean-field phenomenological theories for hydration forces between mesoscopic surfaces.

Multiscale modeling of double-helical DNA and RNA: a unification through Lie groups

Several different mechanical models of double-helical nucleic-acid structures that have been presented in the literature are reviewed here together with a new analysis method that provides a reconciliation between these disparate models. In all cases, terminology and basic results from the theory of Lie groups are used to describe rigid-body motions in a coordinate-free way, and when necessary, coordinates are introduced in a way in which simple equations result. We consider double-helical DNAs and RNAs which, in their unstressed referential state, have backbones that are either straight, slightly precurved, or bent by the action of a protein or other bound molecule. At the coarsest level, we consider worm-like chains with anisotropic bending stiffness. Then, we show how bi-rod models converge to this for sufficiently long filament lengths. At a finer level, we examine elastic networks of rigid bases and show how these relate to the coarser models. Finally, we show how results from molecular dynamics simulation at full atomic resolution (which is the finest scale considered here) and AFM experimental measurements (which is at the coarsest scale) relate to these models.

Single-molecule spectroscopic study of dynamic nanoscale DNA bending behavior of HIV-1 nucleocapsid protein

We have studied the conformational dynamics associated with the nanoscale DNA bending induced by human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein using single-molecule Förster resonance energy transfer (SM-FRET). To gain molecular-level insights into how the HIV-1 NC locally distorts the structures of duplexed DNA segments, the dynamics, reversibility, and sequence specificity of the DNA bending behavior of NC have been systematically studied. We have performed SM-FRET measurements on a series of duplexed DNA segments with varying sequences, lengths, and local structures in the presence of the wide-type HIV-1 NC and NC mutants lacking either the basic N-terminal domain or the zinc fingers. On the basis of the SM-FRET results, we have proposed a possible mechanism for the NC-induced DNA bending in which both NC's zinc fingers and N-terminal domain are found to play crucial roles. The SM-FRET results reported here add new mechanistic insights into the biological behaviors and functions of HIV-1 NC as a retroviral DNA-architectural protein which may play critical roles in the compaction, nuclear import, and integration of the proviral DNA during the retroviral life cycle.

Microscopic basis for the mesoscopic extensibility of dendrimer-compacted DNA

The mechanism of DNA compaction by dendrimers is key to the design of nanotechnologies that can deliver genetic material into cells. We present atomistic simulations, mesoscopic modeling and single-molecule pulling experiments describing DNA dendrimer interactions. All-atom molecular dynamics were used to characterize pulling-force-dependent interactions between DNA and generation-3 PAMAM amine-terminated dendrimers, and a free energy profile and mean forces along the interaction coordinate are calculated. The energy, force, and geometry parameters computed at the atomic level are input for a Monte Carlo model yielding mesoscopic force-extension curves. Actual experimental single-molecule curves obtained with optical tweezers are also presented, and they show remarkable agreement with the virtual curves from our model. The calculations reveal the microscopic origin of the hysteresis observed in the phase transition underlying compaction. A broad range of ionic and pulling parameters is sampled, and suggestions for windows of conditions to probe new single-molecule behavior are made.

DNA linking number change induced by sequence-specific DNA-binding proteins

Sequence-specific DNA-binding proteins play a key role in many fundamental biological processes, such as transcription, DNA replication and recombination. Very often, these DNA-binding proteins introduce structural changes to the target DNA-binding sites including DNA bending, twisting or untwisting and wrapping, which in many cases induce a linking number change (ΔLk) to the DNA-binding site. Due to the lack of a feasible approach, ΔLk induced by sequence-specific DNA-binding proteins has not been fully explored. In this paper we successfully constructed a series of DNA plasmids that carry many tandem copies of a DNA-binding site for one sequence-specific DNA-binding protein, such as λ O, LacI, GalR, CRP and AraC. In this case, the protein-induced ΔLk was greatly amplified and can be measured experimentally. Indeed, not only were we able to simultaneously determine the protein-induced ΔLk and the DNA-binding constant for λ O and GalR, but also we demonstrated that the protein-induced ΔLk is an intrinsic property for these sequence-specific DNA-binding proteins. Our results also showed that protein-mediated DNA looping by AraC and LacI can induce a ΔLk to the plasmid DNA templates. Furthermore, we demonstrated that the protein-induced ΔLk does not correlate with the protein-induced DNA bending by the DNA-binding proteins.

Regulation of porcine pituitary glycoprotein hormone alpha subunit gene with LIM-homeobox transcription factor Lhx3

The aim of this study was to characterize the promoter activity of the porcine pituitary glycoprotein hormone common alpha gene (Cga) promoter (-1059/+12) and the role of LIM homeodomain transcription factor Lhx3. A transfection assay using Chinese hamster ovary (CHO) cells showed that the -1059/-101 region of the Cga promoter definitely responds to Lhx3 and that the -1059/-240 region exhibits a high basal transcriptional level in a pituitary-derived cell line, LbetaT2. A DNA binding and DNase I footprinting assay demonstrated that Lhx3 has seven binding sites in the -1059/+12 region of Cga, including a pituitary glycoprotein hormone basal element (PGBE) known as a LIM homeodomain factor-binding site. A transfection assay of the sequence of Lhx3-binding sites fused with minimal promoter vector confirmed their Lhx3-dependent stimulations in LbetaT2 cells. RT-PCR analysis of porcine pituitary ontogeny demonstrated that porcine Lhx3 showed striking changes of expression in both sexes during the fetal period but a stable high level of expression after birth. Thus, the porcine Cga promoter is regulated by Lhx3 through seven sites in the distal and proximal regions.

Positively charged C-terminal subdomains of EcoRV endonuclease: contributions to DNA binding, bending, and cleavage

The carboxy-terminal subdomains of the homodimeric EcoRV restriction endonuclease each bear a net charge of +4 and are positioned on the inner concave surface of the 50 degree DNA bend that is induced by the enzyme. A complete kinetic and structural analysis of a truncated EcoRV mutant lacking these domains was performed to assess the importance of this diffuse charge in facilitating DNA binding, bending, and cleavage. At the level of formation of an enzyme-DNA complex, the association rate for the dimeric mutant enzyme was sharply decreased by 10(3)-fold, while the equilibrium dissociation constant was weakened by nearly 10(6)-fold compared with that of wild-type EcoRV. Thus, the C-terminal subdomains strongly stabilize the enzyme-DNA ground-state complex in which the DNA is known to be bent. Further, the extent of DNA bending as observed by fluorescence resonance energy transfer was also significantly decreased. The crystal structure of the truncated enzyme bound to DNA and calcium ions at 2.4 A resolution reveals that the global fold is preserved and suggests that a divalent metal ion crucial to catalysis is destabilized in the active site. This may explain the 100-fold decrease in the rate of metal-dependent phosphoryl transfer observed for the mutant. These results show that diffuse positive charge associated with the C-terminal subdomains of EcoRV plays a key role in DNA association, bending, and cleavage.

Conformational Statistics of Semi-Flexible Macromolecular Chains with Internal Joints

Fluctuations in the bending angles at internal irregularities of DNA and RNA (such as symmetric loops, bulges, and nicks/gaps) have been observed from various experiments. However, little effort has been made to computationally predict and explain the statistical behavior of semi-flexible chains with internal defects. In this paper, we describe the general structure of these macromolecular chains as inextensible elastic chains with one or more internal joints which have limited ranges of rotation, and propose a method to compute the probability density functions of the end-to-end pose of these macromolecular chains. Our method takes advantage of the operational properties of the non-commutative Fourier transform for the group of rigid-body motions in three-dimensional space, SE(3). Two representative types of joints, the hinge for planar rotation and the ball joint for spatial rotation, are discussed in detail. The proposed method applies to various stiffness models of semi-flexible chain-like macromolecules. Examples are calculated using the Kratky-Porod model with specified stiffness, angular fluctuation, and joint locations. Entropic effects associated with internal angular fluctuations of semi-flexible macromolecular chains with internal joints can be computed using this formulation. Our method also provides a potential tool to detect the existence of internal irregularities.

Analysis of scanning force microscopy images of protein-induced DNA bending using simulations

Bending of DNA is a feature essential to the function of many DNA-binding proteins. Bending angles can be estimated with a variety of techniques, but most directly from images obtained using scanning force microscopy (SFM). Direct measurement of the bending angle using a tangent method often produces angles that deviate significantly from values obtained using other techniques. Here, we describe the application of SFM in combination with simulations of DNA as a means to estimate protein-induced bending angles in a reliable and unbiased fashion. In this manner, we were able to obtain accurate estimates for the bending angles induced by nuclear factor I, octamer-binding transcription factor 1, the human XPC-Rad23B complex and integration host factor [correction]

A precise DNA bend angle is essential for the function of the phage phi29 transcriptional regulator

Bacteriophage φ29 protein p4 is essential for the regulation of the switch from early to late phage transcription. The protein binds to two regions of the phage genome located between the regulated promoters. Each region contains two inverted repeats separated by 1 bp. We used circular permutation assays to study the topology of the DNA upon binding of the protein and found that p4 induced the same extent of bending independent of the topology of the binding region. In addition, the results revealed that the p4-induced bending is not dependent on the affinity to the binding site but is intrinsic to p4 binding. Independent binding sites were identified through the characterization of the minimal sequence required for p4 binding. The protein has different affinity for each of its binding sites, with those overlapping the A2c and A2b promoter cores (sites 1 and 3), having the highest affinity. The functionality of the p4 binding sites and the contribution of p4-mediated promoter restructuring in transcription regulation is discussed.

Core promoter elements of eukaryotic genes have a highly distinctive mechanical property

In spite of the abundant data on DNA sequence, the mechanical aspects of promoter DNA remain poorly understood. We classified 1871 human and 196 mouse RNA polymerase II promoters and investigated average flexibility profiles of the human promoters containing either a TATA box or an initiator (Inr) sequence only. Here, we show that TATA boxes and Inr sequences have a common anomalous mechanical property: they are comprised of distinctively flexible and rigid sequences, compared with the other parts of the promoter region. The +2 position in the Inr consensus sequence does not favor adenine to keep the high flexibility and thus this position is more accurately represented as 'T, G, C>>A'. Additionally, it was also found that DNA region upstream of TATA box or Inr sequence is more rigid than region downstream of each element. These properties may function as a marker for recognition by TATA-binding protein and Inr-binding protein.

Bending of adenovirus origin DNA by nuclear factor I as shown by scanning force microscopy is required for optimal DNA replication

Nuclear factor I (NFI) is a transcription factor that binds to the adenovirus type 5 (Ad5) origin of replication and recruits the adenovirus DNA polymerase, thereby stimulating initiation of DNA replication in vitro. Using scanning force microscopy, we demonstrate that NFI induces a 60 degrees bend upon binding to the origin. The A/T-rich region preceding the core recognition sequence of NFI influences the DNA bend angle, since substitution of A/T base pairs by G/C base pairs severely decreases bending. Mutations in the A/T-rich region do not affect binding of NFI to DNA. However, mutations that reduce the protein-induced bend lead to a loss of NFI-stimulated replication, indicating that DNA bending is functionally important. In contrast, basal initiation or DNA binding of the polymerase is not impaired by these origin mutations. We conclude that binding of NFI to the Ad5 origin causes structural changes in DNA that are essential for the stimulatory function of NFI in replication. We propose that NFI-induced origin bending facilitates the assembly of a functional initiation complex.

Statistical mechanics of sequence-dependent circular DNA and its application for DNA cyclization

DNA cyclization is potentially the most powerful approach for systematic quantitation of sequence-dependent DNA bending and flexibility. We extend the statistical mechanics of the homogeneous DNA circle to a model that considers discrete basepairs, thus allowing for inhomogeneity, and apply the model to analysis of DNA cyclization. The theory starts from an iterative search for the minimum energy configuration of circular DNA. Thermodynamic quantities such as the J factor, which is essentially the ratio of the partition functions of circular and linear forms, are evaluated by integrating the thermal fluctuations around the configuration under harmonic approximation. Accurate analytic expressions are obtained for equilibrium configurations of homogeneous circular DNA with and without bending anisotropy. J factors for both homogeneous and inhomogeneous DNA are evaluated. Effects of curvature, helical repeat, and bending and torsional flexibility in DNA cyclization are analyzed in detail, revealing that DNA cyclization can detect as little as one degree of curvature and a few percent change in flexibility. J factors calculated by our new approach are well consistent with Monte Carlo simulations, whereas the new theory has much greater efficiency in computations. Simulation of experimental results has been demonstrated.

The curved DNA structure in the 5'-upstream region of the light-responsive genes: its universality, binding factor and function for cyanobacterial psbA transcription

A unique DNA curvature, the CIT, has been found in the 5'-upstream region of the psbA2 gene, which exhibits basal, light-responsive and circadian rhythmic transcription, in a unicellular photosynthetic cyanobacterium, Microcystis aeruginosa K-81. In this study, we report the universality of curvatures found in 5'-upstream regions in the psbA family and the function of the curvature in gene expression. Intrinsic curvatures were identified within 1000 bp upstream from the psbA genes in another cyanobacterium, a red alga and in plants (monocot and dicot). Mutagenized curvatures were constructed and confirmed to have disrupted architecture by gel electrophoresis and atomic force microscopy. Relatively small amounts but light-responsive transcripts of psbA2 were observed in cyanobacterial transformants harboring the mutagenized curvature under light/dark and light/high-light conditions. This shows that the curvature is important for basal transcription. In vitro primer extension and DNA mobility shift assay revealed that factors which might bind to the region upstream from the bending center contribute to the effective basal transcription of psbA2.

Multimerization of restriction fragments by magnesium-mediated stable base pairing between overhangs: a cause of electrophoretic mobility shift

The electrophoretic mobility shift assay (EMSA) or simply the "gel shift assay" is one of the most sensitive methods for studying the ability of a protein to bind to DNA. EMSAs are also widely used to investigate protein- or sequence-dependent DNA bending. Here we report that electrophoresis using physiological concentrations of Mg(2+) can cause a mobility shift of restriction fragments in nondenaturing polyacrylamide gels as the overhangs form stable base pairs. This phenomenon was observed at even 37 degrees C. The retardation was, however, more pronounced at low temperatures, where a three-nucleotide overhang 5'-GAC also caused a mobility shift. The stability of the pairing was generally high when the overhangs of four nucleotides display high GC content, while the mobility shift caused by 5'-AATT was greater than those caused by 5'-GATC, 5'-TCGA, and 5'-CTAG. This observation should be taken into account to avoid misinterpretation of the data when the EMSA, especially the circular permutation gel mobility shift assay, is performed using a running buffer that contains Mg(2+) ions. The stable adhesion between short overhangs may present an important basis for genome stability and many genetic processes occurring in living cells.

Differential modulation of DNA conformation by estrogen receptors alpha and beta

The human estrogen receptor (ER) induces transcription of estrogen-responsive genes upon binding to estrogen and the estrogen response element (ERE). To determine whether receptor-induced changes in DNA structure are related to transactivation, we compared the abilities of ER alpha and ER beta to activate transcription and induce distortion and bending in DNA. ER alpha induced higher levels of transcription than ER beta in the presence of 17 beta-estradiol. In circular permutation experiments ER alpha induced greater distortion in DNA fragments containing the consensus ERE sequence than ER beta. Phasing analysis indicated that ER alpha induced a bend directed toward the major groove of the DNA helix but that ER beta failed to induce a directed DNA bend. Likewise, the ER alpha DNA binding domain (DBD) and hinge region induced a bend directed toward the major groove of the DNA helix, but the ER beta DBD and hinge region failed to bend ERE-containing DNA fragments. Using receptor chimeras we demonstrated that the ER alpha DBD C-terminal extension is required for directed DNA bending. Transient transfection assays revealed that appropriately oriented DNA bending enhances receptor-mediated transactivation. The different abilities of ER alpha and ER beta to induce change in DNA structure could foster or inhibit the interaction of regulatory proteins with the receptor and other transcription factors and help to explain how estrogen-responsive genes are differentially regulated by these two receptors.

The smooth muscle gamma-actin gene promoter is a molecular target for the mouse bagpipe homologue, mNkx3-1, and serum response factor

An evolutionarily conserved vertebrate homologue of the Drosophila NK-3 homeodomain gene bagpipe, Nkx3-1, is expressed in vascular and visceral mesoderm-derived muscle tissues and may influence smooth muscle cell differentiation. Nkx3-1 was evaluated for mediating smooth muscle gamma-actin (SMGA) gene activity, a specific marker of smooth muscle differentiation. Expression of mNkx3-1 in heterologous CV-1 fibroblasts was unable to elicit SMGA promoter activity but required the coexpression of serum response factor (SRF) to activate robust SMGA transcription. A novel complex element containing a juxtaposed Nkx-binding site (NKE) and an SRF-binding element (SRE) in the proximal promoter region was found to be necessary for the Nkx3-1/SRF coactivation of SMGA transcription. Furthermore, Nkx3-1 and SRF associate through protein-protein interactions and the homeodomain region of Nkx3-1 facilitated SRF binding to the complex NKE.SRE. Mutagenesis of Nkx3-1 revealed an inhibitory domain within its C-terminal segment. In addition, mNkx3-1/SRF cooperative activity required an intact Nkx3-1 homeodomain along with the MADS box of SRF, which contains DNA binding and dimerization structural domains, and the contiguous C-terminal SRF activation domain. Thus, SMGA is a novel target for Nkx3-1, and the activity of Nkx3-1 on the SMGA promoter is dependent upon SRF.

A portion of the nucleotide sequence corresponding to the N-terminal coding region of livJ is essential for its transcriptional regulation

We investigated the regulation of the livJ and livKHMGF operons, which are involved in branched-chain amino acid high-affinity transport in Salmonella typhimurium. When livJ was fused to lacZ at the second codon of livJ to make a livJ-lacZ protein fusion, expression from the livJ promoter was not repressed even under repressing growth conditions; however, expression of an analogous construct of livK-lacZ was repressed. When livJ was fused to lacZ at the twelfth codon of livJ, the expression level under unrepressing growth conditions was elevated, resulting in apparent repressibility of the livJ-lacZ protein fusion. Expression from the livJ-lacZ operon fusion, in which livJ was fused to lacZ 159 bp downstream from the A of the start codon of livJ, was relatively normal under unrepressing growth conditions. Deletion analysis and site-directed base-substitution analysis strongly suggested that cis-acting element for regulation of livJ transcription, 5'-GGCAGGATGTATCG-3', starting at +21 and ending at +34 downstream from the A of the start codon of livJ, was present in the N-terminal coding region of livJ.

Preferential platination of an activated cellular promoter by cis-diamminedichloroplatinum

This study examines how accessibility to cisplatin on various genomic regions in T47D breast cancer cells, including the retinoic acid receptor beta gene promoter and coding region and the dihydrofolate reductase gene promoter and coding region, is affected by treatment of the cells with 9-cis retinoic acid, a treatment that activates the retinoic acid receptor beta gene promoter in these cells. A PCR-based assay was used to measure cisplatin adduct density based on the inhibition of PCR amplification of templates from cisplatin treated versus untreated cells. Treatment of cells with 9-cis retinoic acid enhanced accessibility to cisplatin on the retinoic acid receptor beta gene promoter region, but not on the coding regions of that gene nor on the dihydrofolate reductase gene promoter or coding regions, where accessibilities to cisplatin remained 2-4 times lower than on the activated retinoic acid receptor beta gene promoter. Examination of smaller regions within this promoter region showed a repression of platination in the 500 bp region surrounding the TATA box in cells prior to 9-cis retinoic acid treatment, which was abolished following promoter activation. Differences in sequence composition between the various regions could not fully account for differences in platination, suggesting that structural features such as bends in retinoic acid receptor beta gene promoter DNA following gene activation, create energetically favorable sites for platination, and contribute to the cytotoxicity of the drug.

DNA damage by nitrogen mustard in a gene containing multiple Sp1-binding sites

The human cytochrome c(1) gene TATA-less promoter contains 10 Sp1-binding elements that regulate the activation of transcription of this gene. Quantitative PCR was used to show that nitrogen mustard induces DNA lesions within this Sp1-binding region following exposure of HeLa cells to clinical levels of the drug. Alkylation of the cytochrome c(1) gene in HeLa cells increased with reaction time up to 4 h following exposure to nitrogen mustard, with 50% of the lesions (approximately 0.8/kb) forming within 1 h. An Sp1 competition assay showed that nitrogen mustard inhibited the binding of Sp1 to the promoter region of the cytochrome c(1) gene in HeLa cells. These results show that nitrogen mustard-induced damage to Sp1-binding sites may contribute to the toxicity of this compound by interfering with the activation of specific genes.

DNA bending induced by the archaebacterial histone-like protein MC1

The conformational changes induced by the binding of the histone-like protein MC1 to DNA duplexes have been analyzed by dark-field electron microscopy and polyacrylamide gel electrophoresis. Visualisation of the DNA molecules by electron microscopy reveals that the binding of MC1 induces sharp kinks. Linear DNA duplexes (176 bp) which contained a preferential site located at the center were used for quantitative analysis. Measurements of the angle at the center of all duplexes, at a fixed DNA concentration, as a function of the MC1 concentration, were very well fitted by a simple model of an isotropic flexible junction and an equilibrium between the two conformations of DNA with bound or unbound MC1. This model amounts to double-folded Gaussian distributions and yields an equilibrium deflection angle of theta0=116 degrees for the DNA with bound MC1. It allowed measurements of the fraction of DNA with bound MC1 to be taken as a function of MC1 concentrations and yields an equilibrium dissociation constant of Kd=100 nM. It shows that the flexibility of DNA is reduced by the binding of MC1 and the formation of a kink. The equilibrium dissociation constant value was corroborated by gel electrophoresis. Control of the model by the computation of the reduced chi2 shows that the measurements are consistent and that electron microscopy can be used to quantify precisely the DNA deformations induced by the binding of a protein to a preferential site.

DNA structure in human RNA polymerase II promoters

The fact that DNA three-dimensional structure is important for transcriptional regulation begs the question of whether eukaryotic promoters contain general structural features independently of what genes they control. We present an analysis of a large set of human RNA polymerase II promoters with a very low level of sequence similarity. The sequences, which include both TATA-containing and TATA-less promoters, are aligned by hidden Markov models. Using three different models of sequence-derived DNA bendability, the aligned promoters display a common structural profile with bendability being low in a region upstream of the transcriptional start point and significantly higher downstream. Investigation of the sequence composition in the two regions shows that the bendability profile originates from the sequential structure of the DNA, rather than the general nucleotide composition. Several trinucleotides known to have high propensity for major groove compression are found much more frequently in the regions downstream of the transcriptional start point, while the upstream regions contain more low-bendability triplets. Within the region downstream of the start point, we observe a periodic pattern in sequence and bendability, which is in phase with the DNA helical pitch. The periodic bendability profile shows bending peaks roughly at every 10 bp with stronger bending at 20 bp intervals. These observations suggest that DNA in the region downstream of the transcriptional start point is able to wrap around protein in a manner reminiscent of DNA in a nucleosome. This notion is further supported by the finding that the periodic bendability is caused mainly by the complementary triplet pairs CAG/CTG and GGC/GCC, which previously have been found to correlate with nucleosome positioning. We present models where the high-bendability regions position nucleosomes at the downstream end of the transcriptional start point, and consider the possibility of interaction between histone-like TAFs and this area. We also propose the use of this structural signature in computational promoter-finding algorithms.

Polymer chain statistics and conformational analysis of DNA molecules with bends or sections of different flexibility

The worm-like chain model has often been employed to describe the average conformation of long, intrinsically straight polymer molecules, including DNA. The present study extends the applicability of the worm-like chain model to polymers containing bends or sections of different flexibility. Several cases have been explicitly considered: (i) polymers with a single bend; (ii) polymers with multiple coplanar bends; (iii) polymers with two non-coplanar bends; and (iv) polymers comprised of sections with different persistence lengths. Expressions describing the average conformation of such polymers in terms of the mean-square end-to-end distance have been derived for each case. For cases (i) and (iv), expressions for the projection of the end-to-end vector onto the initial orientation of the chain are presented. The expressions derived here have been used to investigate DNA molecules with sequence-induced bending (A-tracts). Mean-square end-to-end distance values determined from a large number of A-tract containing DNA molecules visualized by scanning force microscopy resulted in an average bend angle of 13.5 degrees per A-tract. A similar study was performed to characterize the flexibility of double-strandedDNA molecules containing a single-stranded region. Analysis of their mean-square end-to-end distance yielded a persistence length of 1.3 nm for single-stranded DNA.

PTH-responsive osteoblast nuclear matrix architectural transcription factor binds to the rat type I collagen promoter

In connective tissue, cell structure contributes to type I collagen expression. Differences in osteoblast microarchitecture may account for the two distinct cis elements regulating basal expression, in vivo and in vitro, of the rat type I collagen alpha1(I) polypeptide chain (COL1A1). The COL1A1 promoter conformation may be the penultimate culmination of osteoblast structure. Architectural transcription factors bind to the minor groove of AT-rich DNA and bend it, altering interactions between other trans-acting proteins. Similarly, nuclear matrix (NM) proteins bind to the minor groove of AT-rich matrix-attachment regions, regulating transcription by altering DNA structure. We propose that osteoblast NM architectural transcription factors link cell structure to promoter geometry and COL1A1 transcription. Our objective was to identify potential osteoblast NM architectural transcription factors near the in vitro and in vivo regulatory regions of the rat COL1A1 promoter. Nuclear protein-promoter interactions were analyzed by gel shift analysis and related techniques. NM extracts were derived from rat osteosarcoma cells and from rat bone. The NM protein, NMP4, and a soluble nuclear protein, NP, both bound to two homologous poly(dT) elements within the COL1A1 in vitro regulatory region and proximal to the in vivo regulatory element. These proteins bound within the minor groove and bent the DNA. Parathyroid hormone increased NP/NMP4 binding to both poly(dT) elements and decreased COL1A1 mRNA in the osteosarcoma cells. NP/NMP4-COL1A1 promoter interactions may represent a molecular pathway by which osteoblast structure is coupled to COL1A1 expression.