Boian S. Alexandrov

DNA breathing dynamics in the presence of a terahertz field

We consider the influence of a terahertz field on the breathing dynamics of double-stranded DNA. We model the spontaneous formation of spatially localized openings of a damped and driven DNA chain, and find that linear instabilities lead to dynamic dimerization, while true local strand separations require a threshold amplitude mechanism. Based on our results we argue that a specific terahertz radiation exposure may significantly affect the natural dynamics of DNA, and thereby influence intricate molecular processes involved in gene expression and DNA replication.

Mammalian Stem Cells Reprogramming in Response to Terahertz Radiation

We report that extended exposure to broad-spectrum terahertz radiation results in specific changes in cellular functions that are closely related to DNA-directed gene transcription. Our gene chip survey of gene expression shows that whereas 89% of the protein coding genes in mouse stem cells do not respond to the applied terahertz radiation, certain genes are activated, while other are repressed. RT-PCR experiments with selected gene probes corresponding to transcripts in the three groups of genes detail the gene specific effect. The response was not only gene specific but also irradiation conditions dependent. Our findings suggest that the applied terahertz irradiation accelerates cell differentiation toward adipose phenotype by activating the transcription factor peroxisome proliferator-activated receptor gamma (PPARG). Finally, our molecular dynamics computer simulations indicate that the local breathing dynamics of the PPARG promoter DNA coincides with the gene specific response to the THz radiation. We propose that THz radiation is a potential tool for cellular reprogramming.

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.

Specificity and Heterogeneity of Terahertz Radiation Effect on Gene Expression in Mouse Mesenchymal Stem Cells

We report that terahertz (THz) irradiation of mouse mesenchymal stem cells (mMSCs) with a single-frequency (SF) 2.52 THz laser or pulsed broadband (centered at 10 THz) source results in irradiation specific heterogenic changes in gene expression. The THz effect depends on irradiation parameters such as the duration and type of THz source, and on the degree of stem cell differentiation. Our microarray survey and RT-PCR experiments demonstrate that prolonged broadband THz irradiation drives mMSCs toward differentiation, while 2-hour irradiation (regardless of THz sources) affects genes transcriptionally active in pluripotent stem cells. The strictly controlled experimental environment indicates minimal temperature changes and the absence of any discernable response to heat shock and cellular stress genes imply a non-thermal response. Computer simulations of the core promoters of two pluripotency markers reveal association between gene upregulation and propensity for DNA breathing. We propose that THz radiation has potential for non-contact control of cellular gene expression.

Non-thermal effects of terahertz radiation on gene expression in mouse stem cells

Abstract In recent years, terahertz radiation sources are increasingly being exploited in military and civil applications. However, only a few studies have so far been conducted to examine the biological effects associated with terahertz radiation. In this study, we evaluated the cellular response of mesenchymal mouse stem cells exposed to THz radiation. We apply low-power radiation from both a pulsed broad-band (centered at 10 THz) source and from a CW laser (2.52 THz) source. Modeling, empirical characterization, and monitoring techniques were applied to minimize the impact of radiation-induced increases in temperature. qRT-PCR was used to evaluate changes in the transcriptional activity of selected hyperthermic genes. We found that temperature increases were minimal, and that the differential expression of the investigated heat shock proteins (HSP105, HSP90, and CPR) was unaffected, while the expression of certain other genes (Adiponectin, GLUT4, and PPARG) showed clear effects of the THz irradiation after prolonged, broad-band exposure.

DNA dynamics play a role as a basal transcription factor in the positioning and regulation of gene transcription initiation

We assess the role of DNA breathing dynamics as a determinant of promoter strength and transcription start site (TSS) location. We compare DNA Langevin dynamic profiles of representative gene promoters, calculated with the extended non-linear PBD model of DNA with experimental data on transcription factor binding and transcriptional activity. Our results demonstrate that DNA dynamic activity at the TSS can be suppressed by mutations that do not affect basal transcription factor binding–DNA contacts. We use this effect to establish the separate contributions of transcription factor binding and DNA dynamics to transcriptional activity. Our results argue against a purely ‘transcription factor-centric’ view of transcription initiation, suggesting that both DNA dynamics and transcription factor binding are necessary conditions for transcription initiation.

Pre-melting dynamics of DNA and its relation to specific functions

We discuss connections between the nonlinear dynamics of double-stranded DNA, experimental findings, and specific DNA functions. We begin by discussing how thermally induced localized openings (bubbles) of the DNA double strand are important for interpreting dynamic force spectroscopy data. Then we demonstrate a correlation between a sequence-dependent propensity for pre-melting bubble formation and transcription initiation and other regulatory effects in viral DNA. Finally, we discuss the possibility of a connection between DNA dynamics and the ability of repair proteins to recognize ultraviolet (UV) radiation damage sites.

Toward a Detailed Description of the Thermally Induced Dynamics of the Core Promoter

Establishing the general and promoter-specific mechanistic features of gene transcription initiation requires improved understanding of the sequence-dependent structural/dynamic features of promoter DNA. Experimental data suggest that a spontaneous dsDNA strand separation at the transcriptional start site is likely to be a requirement for transcription initiation in several promoters. Here, we use Langevin molecular dynamic simulations based on the Peyrard-Bishop-Dauxois nonlinear model of DNA (PBD LMD) to analyze the strand separation (bubble) dynamics of 80-bp-long promoter DNA sequences. We derive three dynamic criteria, bubble probability, bubble lifetime, and average strand separation, to characterize bubble formation at the transcriptional start sites of eight mammalian gene promoters. We observe that the most stable dsDNA openings do not necessarily coincide with the most probable openings and the highest average strand displacement, underscoring the advantages of proper molecular dynamic simulations. The dynamic profiles of the tested mammalian promoters differ significantly in overall profile and bubble probability, but the transcriptional start site is often distinguished by large (longer than 10 bp) and long-lived transient openings in the double helix. In support of these results are our experimental transcription data demonstrating that an artificial bubble-containing DNA template is transcribed bidirectionally by human RNA polymerase alone in the absence of any other transcription factors.

DNA breathing dynamics distinguish binding from nonbinding consensus sites for transcription factor YY1 in cells

The genome-wide mapping of the major gene expression regulators, the transcription factors (TFs) and their DNA binding sites, is of great importance for describing cellular behavior and phenotypic diversity. Presently, the methods for prediction of genomic TF binding produce a large number of false positives, most likely due to insufficient description of the physiochemical mechanisms of protein–DNA binding. Growing evidence suggests that, in the cell, the double-stranded DNA (dsDNA) is subject to local transient strands separations (breathing) that contribute to genomic functions. By using site-specific chromatin immunopecipitations, gel shifts, BIOBASE data, and our model that accurately describes the melting behavior and breathing dynamics of dsDNA we report a specific DNA breathing profile found at YY1 binding sites in cells. We find that the genomic flanking sequence variations and SNPs, may exert long-range effects on DNA dynamics and predetermine YY1 binding. The ubiquitous TF YY1 has a fundamental role in essential biological processes by activating, initiating or repressing transcription depending upon the sequence context it binds. We anticipate that consensus binding sequences together with the related DNA dynamics profile may significantly improve the accuracy of genomic TF binding sites and TF binding-related functional SNPs.

Binding of Nucleoid-Associated Protein Fis to DNA Is Regulated by DNA Breathing Dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.