Advancing Transcription Factor Binding Site Prediction Using DNA Breathing Dynamics and Sequence Transformers via Cross Attention

Understanding the impact of genomic variants on transcription factor binding and gene regulation remains a key area of research, with implications for unraveling the complex mechanisms underlying various functional effects. Our study delves into the role of DNA's biophysical properties, including thermodynamic stability, shape, and flexibility in transcription factor (TF) binding. We developed a multi-modal deep learning model integrating these properties with DNA sequence data. Trained on ChIP-Seq (chromatin immunoprecipitation sequencing) data in vivo involving 690 TF-DNA binding events in human genome, our model significantly improves prediction performance in over 660 binding events, with up to 9.6% increase in AUROC metric compared to the baseline model when using no DNA biophysical properties explicitly. Further, we expanded our analysis to in vitro high-throughput Systematic Evolution of Ligands by Exponential enrichment (SELEX) and Protein Binding Microarray (PBM) datasets, comparing our model with established frameworks. The inclusion of DNA breathing features consistently improved TF binding predictions across different cell lines in these datasets. Notably, for complex ChIP-Seq datasets, integrating DNABERT2 with a cross-attention mechanism provided greater predictive capabilities and insights into the mechanisms of disease-related non-coding variants found in genome-wide association studies. This work highlights the importance of DNA biophysical characteristics in TF binding and the effectiveness of multi-modal deep learning models in gene regulation studies.

Examining DNA breathing with pyDNA-EPBD

MOTIVATION

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

RESULTS

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm.

AVAILABILITY AND IMPLEMENTATION

pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.

Examining DNA Breathing with pyDNA-EPBD

Motivation

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion.This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

Results

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop- Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo (MCMC) algorithm.

Describing nonequilibrium soft matter with mean field game theory

We demonstrate that combining an emerging approach to game theory with self-consistent mean field theory provides realistic treatments of diblock copolymer phase evolution. We especially examine order-order phase transformations upon quenched temperature change involving hexagonal cylinders, lamellae, and the gyroid. Our findings demonstrate that (i) the game theoretical dynamics produce realistic trajectories for the evolution of the local compositions, (ii) the predicted small-angle scattering follows experimentally observed trends, (iii) nucleation and growth is active when the system is quenched far from the critical point, and (iv) epitaxial growth is manifest. To our knowledge, the methodology presented provides the first merger of mean field game theory and statistical mechanics for soft matter systems, giving a new inroad to studying polymer dynamics.

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

The role of structural parameters in DNA cyclization

BACKGROUND

The intrinsic bendability of DNA plays an important role with relevance for myriad of essential cellular mechanisms. The flexibility of a DNA fragment can be experimentally and computationally examined by its propensity for cyclization, quantified by the Jacobson-Stockmayer J factor. In this study, we use a well-established coarse-grained three-dimensional model of DNA and seven distinct sets of experimentally and computationally derived conformational parameters of the double helix to evaluate the role of structural parameters in calculating DNA cyclization.

RESULTS

We calculate the cyclization rates of 86 DNA sequences with previously measured J factors and lengths between 57 and 325 bp as well as of 20,000 randomly generated DNA sequences with lengths between 350 and 4000 bp. Our comparison with experimental data is complemented with analysis of simulated data.

CONCLUSIONS

Our data demonstrate that all sets of parameters yield very similar results for longer DNA fragments, regardless of the nucleotide sequence, which are in agreement with experimental measurements. However, for DNA fragments shorter than 100 bp, all sets of parameters performed poorly yielding results with several orders of magnitude difference from the experimental measurements. Our data show that DNA cyclization rates calculated using conformational parameters based on nucleosome packaging data are most similar to the experimental measurements. Overall, our study provides a comprehensive large-scale assessment of the role of structural parameters in calculating DNA cyclization rates.

Allostery through protein-induced DNA bubbles

Allostery through DNA is increasingly recognized as an important modulator of DNA functions. Here, we show that the coalescence of protein-induced DNA bubbles can mediate allosteric interactions that drive protein aggregation. We propose that such allostery may regulate DNA's flexibility and the assembly of the transcription machinery. Mitochondrial transcription factor A (TFAM), a dual-function protein involved in mitochondrial DNA (mtDNA) packaging and transcription initiation, is an ideal candidate to test such a hypothesis owing to its ability to locally unwind the double helix. Numerical simulations demonstrate that the coalescence of TFAM-induced bubbles can explain experimentally observed TFAM oligomerization. The resulting melted DNA segment, approximately 10 base pairs long, around the joints of the oligomers act as flexible hinges, which explains the efficiency of TFAM in compacting DNA. Since mitochondrial polymerase (mitoRNAP) is involved in melting the transcription bubble, TFAM may use the same allosteric interaction to both recruit mitoRNAP and initiate transcription.

Determination of the shear viscosity of the one-component plasma

The shear viscosity coefficient of the one-component plasma is calculated with unprecedented accuracy using equilibrium molecular dynamics simulations and the Green-Kubo relation. Numerical and statistical uncertainties and their mitigation for improving accuracy are analyzed. In the weakly coupled regime, our results agree with the Landau-Spitzer prediction. In the moderately and strongly coupled regimes, our results are found in good agreement with recent results obtained for the Yukawa one-component plasma using nonequilibrium molecular dynamics. A practical formula is provided for evaluating the viscosity coefficient across coupling regimes, from the weakly coupled regime to solidification threshold. The results are used to test theoretical predictions of the viscosity coefficients found in the literature.

Statistical mechanics of a discrete Schrödinger equation with saturable nonlinearity

We study the statistical mechanics of the one-dimensional discrete nonlinear Schrödinger (DNLS) equation with saturable nonlinearity. Our study represents an extension of earlier work [Phys. Rev. Lett. 84, 3740 (2000)] regarding the statistical mechanics of the one-dimensional DNLS equation with a cubic nonlinearity. As in this earlier study, we identify the spontaneous creation of localized excitations with a discontinuity in the partition function. The fact that this phenomenon is retained in the saturable DNLS is nontrivial, since in contrast to the cubic DNLS whose nonlinear character is enhanced as the excitation amplitude increases, the saturable DNLS, in fact, becomes increasingly linear as the excitation amplitude increases. We explore the nonlinear dynamics of this phenomenon by direct numerical simulations.

Models for energy and charge transport and storage in biomolecules

Two models for energy and charge transport and storage in biomolecules are considered. A model based on the discrete nonlinear Schrödinger equation with long-range dispersive interactions (LRI's) between base pairs of DNA is offered for the description of nonlinear dynamics of the DNA molecule. We show that LRI's are responsible for the existence of an interval of bistability where two stable stationary states, a narrow, pinned state and a broad, mobile state, coexist at each value of the total energy. The possibility of controlled switching between pinned and mobile states is demonstrated. The mechanism could be important for controlling energy storage and transport in DNA molecules. Another model is offered for the description of nonlinear excitations in proteins and other anharmonic biomolecules. We show that in the highly anharmonic systems a bound state of Davydov and Boussinesq solitons can exist.

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.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

Anharmonic dynamics of intramolecular hydrogen bonds driven by DNA breathing

We study the effects of the anharmonic strand-separation dynamics of double-stranded DNA on the infrared spectra of the intramolecular base-pairing hydrogen bonds. Using the extended Peyrard-Bishop-Dauxois model for the DNA breathing dynamics coupled with the Lippincott-Schroeder potential for N-H· · ·N and N-H· · ·O hydrogen bonding, we identify a high-frequency (~96 THz) feature in the infrared spectra. We show that this sharp peak arises as a result of the anharmonic base-pair breathing dynamics of DNA. In addition, we study the effects of friction on the infrared spectra. For higher temperatures (~300 K), where the anharmonicity of DNA dynamics is pronounced, the high-frequency peak is always present irrespective of the friction strength.

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.

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

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 is likely to be a factor in the genomic nucleotide repeats expansions related to diseases

Trinucleotide repeats sequences (TRS) represent a common type of genomic DNA motif whose expansion is associated with a large number of human diseases. The driving molecular mechanisms of the TRS ongoing dynamic expansion across generations and within tissues and its influence on genomic DNA functions are not well understood. Here we report results for a novel and notable collective breathing behavior of genomic DNA of tandem TRS, leading to propensity for large local DNA transient openings at physiological temperature. Our Langevin molecular dynamics (LMD) and Markov Chain Monte Carlo (MCMC) simulations demonstrate that the patterns of openings of various TRSs depend specifically on their length. The collective propensity for DNA strand separation of repeated sequences serves as a precursor for outsized intermediate bubble states independently of the G/C-content. We report that repeats have the potential to interfere with the binding of transcription factors to their consensus sequence by altered DNA breathing dynamics in proximity of the binding sites. These observations might influence ongoing attempts to use LMD and MCMC simulations for TRS–related modeling of genomic DNA functionality in elucidating the common denominators of the dynamic TRS expansion mutation with potential therapeutic applications.

Feigenbaum cascade of discrete breathers in a model of DNA

We demonstrate that period-doubled discrete breathers appear from the anticontinuum limit of the driven Peyrard-Bishop-Dauxois model of DNA. These novel breathers result from a stability overlap between subharmonic solutions of the driven Morse oscillator. Subharmonic breathers exist whenever a stability overlap is present within the Feigenbaum cascade to chaos and therefore an entire cascade of such breathers exists. This phenomenon is present in any driven lattice where the on-site potential admits subharmonic solutions. In DNA these breathers may have ramifications for cellular gene expression.

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.

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.

Nonlinearity in DNA and its relation to specific functions

In memory of Alwyn Scott, we discuss the connection 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 the sequence-dependent propensity for 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-radiation damage.

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.

Accessing the information encoded in DNA requires that RNA polymerases recognize the core promoter, a sequence that marks the start of a gene. Statistical analysis of known promoter sequences has failed to reveal a simple code for identifying promoters, leading to the suggestion that promoter DNA is distinguished by certain structural/dynamic properties encoded in nonobvious ways by the literal sequence. Because the DNA strands at the promoter need to be separated for transcription to begin, we previously proposed that promoter sequences exhibit a propensity for spontaneous strand separation. Here, we conduct simulations of the ultrafast, small-scale strand separation motions of eight mammalian promoters and show that start sites tend to form larger and more stable openings in the double helix compared to other sequences. Experimentally, we show that an artificial permanent opening in the double helix is sufficient for transcription in the absence of sequence-specific protein–DNA contacts. These findings support a view of DNA as a structurally active participant in gene expression, rather than the commonly envisioned passive digital storage device. Our analysis suggests that functionally relevant structural variation in genomic DNA occurs at the level of fast motions not readily observed by traditional molecular structure analysis.

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.

Ordering and reverse ordering mechanisms of triblock copolymers in the presence of solvent

Self-consistent field theory is used to study the self-assembly of a triblock copolymer melt. Two different external factors (temperature and solvent) are shown to affect the self-assembly. Either one or two-step self-assembly can be found as a function of temperature in the case of a neat triblock melt, or as a function of increasing solvent content (for non-selective solvents) in the case of a triblock-solvent mixture. For selective solvents, it is shown that increasing the solvent content leads to more complicated self-assembly mechanisms, including a reversed transition where order is found to increase instead of decreasing as expected, and re-entrant behavior where order is found to increase at first, and then decrease to a previous state of disorder.

Stress distributions in diblock copolymers

We demonstrate how a generalized self-consistent field theory for polymer melts that includes elastic stress and strain fields can be applied to the study of AB diblock copolymers melts. By obtaining the stress distributions for volume conserving strain loadings where lamellar and hexagonal morphologies are stable, we show that the local stress is reduced at the domain interface but slightly enhanced in the immediate vicinity of the interface. The overall stress profile is the result of the combined effects of chain connectivity across the interface, which yields a positive contribution, and the immiscible nature of the monomers, which leads to a stress reduction because of interfacial tension.

Opening rates of DNA hairpins: experiment and model

We present single-molecule measurements of the opening rate of DNA hairpins under mechanical tension and compare with the results obtained from a reduced-degrees-of-freedom statistical mechanics model. We extract the apparent position of the transition state s and find that the model, with no fitting parameters, reproduces the experimental measurements surprisingly well. Our values for s are different from the ones obtained in previous experiments, where, however, the experimental conditions were different (different force fields, different salt concentrations). Thus it appears that the values of s measured for relatively short hairpins are strongly affected by these experimental conditions.

Bubble statistics and dynamics in double-stranded DNA

The dynamical properties of double-stranded DNA are studied in the framework of the Peyrard-Bishop-Dauxois model using Langevin dynamics. Our simulations are analyzed in terms of two distribution functions describing localized separations ("bubbles") of the double strand. The result that the bubble distributions are more sharply peaked at the active sites than thermodynamically obtained distributions is ascribed to the fact that the bubble lifetimes affect the distributions. Certain base-pair sequences are found to promote long-lived bubbles, and we argue that this is a result of length scale competition between the nonlinearity and disorder present in the system.

Selectivity- and size-induced segregation of molecular and nanoscale species in microphase-ordered triblock copolymers

Microphase-ordered block copolymers serve as model systems to elucidate the potential of molecular self-assembly and organic templates to fabricate functionalized polymeric materials. Both aspects are related to the incorporation of secondary species such as low-molar-mass compounds or nanoparticles within the copolymer matrices. Since the resulting properties of such functionalized copolymers critically depend on the morphology of the blend or composite, the nonrandom distribution of such inclusions within the copolymer matrix must be understood. Using a self-consistent field theoretical approach, we quantitatively evaluate the segregation and interfacial excess of low-molar-mass and nanoscale species in ordered triblock copolymers as functions of block selectivity and inclusion size. The predictions are found to agree with the morphology observed in a model triblock copolymer/nanoparticle composite, thereby demonstrating the generality of this approach. Our results suggest a wide correspondence in the structure-forming effect of molecular and nanoscale inclusions that will have implications in the design and processing of functional nanostructured polymers.

Molecularly asymmetric triblock copolymers as a single-molecule route to ordered bidisperse polymer brushes

The conditions signaling the formation of bidisperse brushes in ordered block copolymers are investigated as an A(2) block is progressively grown onto an A(1)B diblock copolymer to form a series of molecularly asymmetric, isomorphic A(1)BA(2) triblock copolymers. Small-angle scattering and self-consistent field theory confirm that the microphase-ordered period decreases when the A(2) block is short relative to the A(1) block, but then increases as A(1)+A(2) bidisperse brushes develop. The mechanical properties systematically follow the spatial distribution of the A(2) block.

Discrete nonlinear Schrödinger equations with arbitrarily high-order nonlinearities

A class of discrete nonlinear Schrödinger equations with arbitrarily high-order nonlinearities is introduced. These equations are derived from the same Hamiltonian using different Poisson brackets and include as particular cases the saturable discrete nonlinear Schrödinger equation and the Ablowitz-Ladik equation. As a common property, these equations possess three kinds of exact analytical stationary solutions for which the Peierls-Nabarro barrier is zero. Several properties of these solutions, including stability, discrete breathers, and moving solutions, are investigated.

Probing the mechanical unzipping of DNA

A study of the micromechanical unzipping of DNA in the framework of the Peyrard-Bishop-Dauxois model is presented. We introduce a Monte Carlo technique that allows accurate determination of the dependence of the unzipping forces on unzipping speed and temperature. Our findings agree quantitatively with experimental results for homogeneous DNA, and for lamda-phage DNA we reproduce the recently obtained experimental force-temperature phase diagram. Finally, we argue that there may be fundamental differences between in vivo and in vitro DNA unzipping.

Statistical mechanics of general discrete nonlinear Schrödinger models: localization transition and its relevance for Klein-Gordon lattices

We extend earlier work [Phys. Rev. Lett. 84, 3740 (2000)]] on the statistical mechanics of the cubic one-dimensional discrete nonlinear Schrödinger (DNLS) equation to a more general class of models, including higher dimensionalities and nonlinearities of arbitrary degree. These extensions are physically motivated by the desire to describe situations with an excitation threshold for creation of localized excitations, as well as by recent work suggesting noncubic DNLS models to describe Bose-Einstein condensates in deep optical lattices, taking into account the effective condensate dimensionality. Considering ensembles of initial conditions with given values of the two conserved quantities, norm and Hamiltonian, we calculate analytically the boundary of the "normal" Gibbsian regime corresponding to infinite temperature, and perform numerical simulations to illuminate the nature of the localization dynamics outside this regime for various cases. Furthermore, we show quantitatively how this DNLS localization transition manifests itself for small-amplitude oscillations in generic Klein-Gordon lattices of weakly coupled anharmonic oscillators (in which energy is the only conserved quantity), and determine conditions for the existence of persistent energy localization over large time scales.

Multipeaked polarons in soft potentials

We consider a minimal coupled charge / excitation-lattice model capturing a competition between linear polaronic self-trapping and the self-focusing effects of a soft nonlinear on-site potential. The standard single-humped polaron ceases to exist above a critical value of the coupling strength, closely related to the inflection point in the nonlinear potential. For couplings beyond this critical value, we find that successive multihumped polaronic solutions correspond to the lowest-energy stationary states of the system, which may admit interesting quantum resonance behavior.

Elastic moduli of multiblock copolymers in the lamellar phase

We study the linear elastic response of multiblock copolymer melts in the lamellar phase, where the molecules are composed of tethered symmetric AB diblock copolymers. We use a self-consistent field theory method, and introduce a real space approach to calculate the tensile and shear moduli as a function of block number. The former is found to be in qualitative agreement with experiment. We find that the increase in bridging fraction with block number, that follows the increase in modulus, is not responsible for the increase in modulus. It is demonstrated that the change in modulus is due to an increase in mixing of repulsive A and B monomers. Under extension, this increase originates from a widening of the interface, and more molecules pulled free of the interface. Under compression, only the second of these two processes acts to increase the modulus.

Improved convergence in block copolymer self-consistent field theory by Anderson mixing

A modification to real space polymeric self-consistent field theory algorithms that greatly improves the convergence properties is presented. The method is based on Anderson mixing [D. G. Anderson, J. Assoc. Comput. Mach. 12, 547 (1965)], and each iteration computed takes negligibly longer to perform than with other methods, but the number of iterations required to reach a high accuracy solution is greatly reduced. No a priori knowledge of possible phases is required to apply this method. We apply our approach to a standard diblock copolymer melt, and demonstrate iteration reductions of more than a factor of 5 in some cases.

Hysteresis and metastability in the quenched turbulent dynamics of the complex Ginzburg-Landau equation

We consider the quenched dynamics of the two-dimensional complex Ginzburg-Landau equation in its turbulent regime. We initialize the system in a frustrated state and observe how frustration affects the evolution towards the turbulent state. This process is performed for parameter values where, for random initial conditions, the system evolves into the turbulent state. We observe that the glassiness of the initial condition can inhibit the occurrence of the absolute instability close to the critical point for that instability in parameter space. Sufficiently far from the critical point, the turbulent state will develop, but only after spending considerable time in a transient metastable state of fixed vortex density. The parameter distance from the critical point is found to scale as an exponential of a power of the lifetime of the metastable state, and with a power exponent depending on the "depth" of the original quench. The limiting regimes of shallow and deep quench are identified by their respective values of the exponent, and the distinct mechanisms leading to the relaxation to turbulence in each case are highlighted.

Three-dimensional elastic compatibility and varieties of twins in martensites

We model a cubic-to-tetragonal martensitic transition by a Ginzburg-Landau free energy in the symmetric strain tensor. We show in three dimensions (3D) that solving the St. Venant compatibility relations for strain, treated as independent field equations, generates three anisotropic long-range potentials between the two order parameter components. These potentials encode 3D discrete symmetries, express the energetics of lattice integrity, and determine 3D textures. Simulation predictions include twins with temperature-varying orientation, helical twins, competing metastable states, and compatibility-induced elastic frustration. Our approach also applies to improper ferroelastics.

Twisted localized modes

In a number of recent papers the so-called twisted localized mode of the discrete nonlinear Schrödinger equation has been proposed. Herein, we study the existence and stability properties of such modes. We analyze the persistence of quasiperiodic modes and study the domains of existence and numerical stability of the exact form of such solutions. We identify the bifurcations through which they lose their stability and follow the behavior of the intrinsic localized modes and their eigenmodes even in the unstable regime.