Davit A. Potoyan

Recent successes in coarse-grained modeling of DNA

The growing interest in the DNA‐based mesoscale systems of biological and nonbiological nature has encouraged the computational molecular science community to develop coarse‐grained (CG) representationsof the DNA that will be simple enough to permit exhaustive simulations in a reasonable amount of time, yet complex enough to capture the essential physics at play. In the recent years, there have been some major developments in the DNA coarse‐graining area and several fairly sophisticated models are now available that faithfully reproduce key mechanical and chemical properties of the double‐ and single‐stranded DNA. However, there are still many challenges, which limit the applicability of the present models, and much has to be done yet to develop more reliable schemes which would have a predictive power beyond the target domain of the intrinsic parametrization. A development of robust, controllable, and transferrable CG DNA force fields will provide an invaluable tool for gaining physical insights into the molecular nature of complex DNA‐based nanoscale entities such as the chromatin, virus capsids, and DNA nanocomposites. In the present contribution, we provide an overview of the recent developments in the DNA coarse‐graining field. Our aim is to review the existing CG models of the double‐stranded DNA, where a small selection of models, which we believe provide avenues for promising future development, are discussed in some detail.

Regulation of the H4 tail binding and folding landscapes via Lys-16 acetylation

Intrinsically disordered proteins (IDP) are a broad class of proteins with relatively flat energy landscapes showing a high level of functional promiscuity, which are frequently regulated through posttranslational covalent modifications. Histone tails, which are the terminal segments of the histone proteins, are prominent IDPs that are implicated in a variety of signaling processes, which control chromatin organization and dynamics. Although a large body of work has been done on elucidating the roles of posttranslational modifications in functional regulation of IDPs, molecular mechanisms behind the observed behaviors are not fully understood. Using extensive atomistic molecular dynamics simulations, we found in this work that H4 tail mono-acetylation at LYS-16, which is a key covalent modification, induces a significant reorganization of the tail’s conformational landscape, inducing partial ordering and enhancing the propensity for alpha-helical segments. Furthermore, our calculations of the potentials of mean force between the H4 tail and a DNA fragment indicate that contrary to the expectations based on simple electrostatic reasoning, the Lys-16 mono-acetylated H4 tail binds to DNA stronger than the unacetylated protein. Based on these results, we propose a molecular mechanism for the way Lys-16 acetylation might lead to experimentally observed disruption of compact chromatin fibers.

Molecular stripping in the NF-κB/IκB/DNA genetic regulatory network

Significance The classic models of genetic switches are based on thermodynamic control via binding affinities. Complex nonequilibrium eukaryotic gene networks with many targets such as the NF-κB/IκB/DNA system challenge this picture. Waiting for unbinding from the numerous NF-κB binding sites takes time that would allow wasteful overexpression of many genes. Active molecular stripping of NF-κB from the DNA by IκB solves the timescale problem raised in the classic picture based on thermodynamics alone. Genetic switches based on the NF-κB/IκB/DNA system are master regulators of an array of cellular responses. Recent kinetic experiments have shown that IκB can actively remove NF-κB bound to its genetic sites via a process called “molecular stripping.” This allows the NF-κB/IκB/DNA switch to function under kinetic control rather than the thermodynamic control contemplated in the traditional models of gene switches. Using molecular dynamics simulations of coarse-grained predictive energy landscape models for the constituent proteins by themselves and interacting with the DNA we explore the functional motions of the transcription factor NF-κB and its various binary and ternary complexes with DNA and the inhibitor IκB. These studies show that the function of the NF-κB/IκB/DNA genetic switch is realized via an allosteric mechanism. Molecular stripping occurs through the activation of a domain twist mode by the binding of IκB that occurs through conformational selection. Free energy calculations for DNA binding show that the binding of IκB not only results in a significant decrease of the affinity of the transcription factor for the DNA but also kinetically speeds DNA release. Projections of the free energy onto various reaction coordinates reveal the structural details of the stripping pathways.

The Acetylation Landscape of the H4 Histone Tail: Disentangling the Interplay between the Specific and Cumulative Effects

Histone tails, the intrinsically disordered terminal regions of histone proteins, are key modulators of the structure and dynamics of chromatin and, consequently, are central to many DNA template-directed processes including replication, repair, and transcription. Acetylation of histone tails is a major post-translational modification (PTM) involved in regulating chromatin, yet it remains unclear how acetylation modifies the disordered state of histone tails and affects their function. We investigated the consequences of increasing acetylation on the isolated H4 histone tail by characterizing the conformational ensembles of unacetylated, mono-, di-, tri-, and tetra-acetylated H4 histone tails using Replica Exchange Molecular Dynamics (REMD) simulations. We found that progressive acetylation has a cumulative effect on the H4 tail, decreasing conformational heterogeneity, increasing helical propensity, and increasing hydrogen bond occupancies. The monoacetylation of lysine 16, however, has unique and specific effects: drastically decreasing the conformational heterogeneity of the H4 tail and leading to highly localized helical secondary structure and elongated conformations. We describe how the cumulative effects of acetylation arise from the charge reduction and increased hydrophobicity associated with adding acetyl groups, while the specific effects are a consequence of steric interactions that are sequence specific. Additionally, we found that increasing the level of acetylation results in the formation of spatially clustered lysines that could serve as recognition patches for binding of chromatin regulating proteins. Hence, we explore the mechanisms by which different acetylation patterns may result in specific recognition of the H4 histone tails by protein or DNA binding partners.

Computing free energy of a large-scale allosteric transition in adenylate kinase using all atom explicit solvent simulations.

During allosteric motions proteins navigate rugged energy landscapes. Hence, mapping of these multidimensional landscapes into lower dimensional manifolds is important for gaining deeper insights into allosteric dynamics. Using a recently developed computational technique, we calculated the free energy difference between the open and closed states of adenylate kinase, an allosteric protein which was extensively studied previously using both experimental and theoretical approaches. Two independent simulations indicate reasonable convergence of the computed free energy profiles. The numerical value of the open/closed free energy difference is only 1-2 k(B)T, much smaller than some of the prior estimates. We also found that the conformations structurally close to the open form still retain many LID-NMP contacts, suggesting that the conformational basin of the closed form is larger than expected. The latter suggestion may explain the discrepancy in relative populations of open and closed forms of unligated adenylate kinase, observed in NMR and FRET experiments.

On the dephasing of genetic oscillators

Significance Genetic oscillators are ubiquitous regulatory motifs in the molecular control circuits of living cells. Prominent examples include the cell cycle and cellular signaling. There are two primary sources of noise in these oscillators: the binary character of gene states and the low copy numbers of proteins. This molecular noise induces dephasing in oscillators in analogy to the way imperfections in clocks lead to disparities in time. We can study how these two distinct sources contribute to dephasing and how well it can be approximated by adding noise to a macroscopic model. Using the NFκB/IκB oscillator as an example, we find that gene noise leads to significant deviations from the often-used phenomenological models. The digital nature of genes combined with the associated low copy numbers of proteins regulating them is a significant source of stochasticity, which affects the phase of biochemical oscillations. We show that unlike ordinary chemical oscillators, the dichotomic molecular noise of gene state switching in gene oscillators affects the stochastic dephasing in a way that may not always be captured by phenomenological limit cycle-based models. Through simulations of a realistic model of the NFκB/IκB network, we also illustrate the dephasing phenomena that are important for reconciling single-cell and population-based experiments on gene oscillators.

Molecular stripping, targets and decoys as modulators of oscillations in the NF-κB/IκBα/DNA genetic network

Eukaryotic transcription factors in the NF-κB family are central components of an extensive genetic network that activates cellular responses to inflammation and to a host of other external stressors. This network consists of feedback loops that involve the inhibitor IκBα, numerous downstream functional targets, and still more numerous binding sites that do not appear to be directly functional. Under steady stimulation, the regulatory network of NF-κB becomes oscillatory, and temporal patterns of NF-κB pulses appear to govern the patterns of downstream gene expression needed for immune response. Understanding how the information from external stress passes to oscillatory signals and is then ultimately relayed to gene expression is a general issue in systems biology. Recently, in vitro kinetic experiments as well as molecular simulations suggest that active stripping of NF-κB by IκBα from its binding sites can modify the traditional systems biology view of NF-κB/IκBα gene circuits. In this work, we revise the commonly adopted minimal model of the NF-κB regulatory network to account for the presence of the large number of binding sites for NF-κB along with dissociation from these sites that may proceed either by passive unbinding or by active molecular stripping. We identify regimes where the kinetics of target and decoy unbinding and molecular stripping enter a dynamic tug of war that may either compensate each other or amplify nuclear NF-κB activity, leading to distinct oscillatory patterns. Our finding that decoys and stripping play a key role in shaping the NF-κB oscillations suggests strategies to control NF-κB responses by introducing artificial decoys therapeutically.

PEST Control of Molecular Stripping of NFκB from DNA Transcription Sites

We recently proposed a model for IκBα-mediated molecular stripping of NFκB from transcription sites. IκBα was shown experimentally to form a transient ternary complex with DNA-bound NFκB, but the mechanism by which the IκBα accelerates dissociation of the NFκB from the DNA was unknown. In this paper we construct and compute free energy profiles for the wild-type IκBα-mediated molecular stripping reaction of NFκB from DNA and compare with that for a mutant of IκBα bearing a charge-neutralized PEST. The differences in the free energy profile for stripping originate from the frustrated electrostatic interactions between the negatively charged PEST and the DNA. The PEST occupies two different conformations in the NFκB-IκBα binary complex, one of which occupies the DNA-binding cavity. Specific interactions with positively charged residues in the N-terminal domains of both p50 and p65 apparently draw the domains closer together hindering reassociation of DNA. Comparison with the charge-neutralized mutant reveals that all of these functional consequences result from the negative charges in the PEST sequence of IκBα.

Dichotomous noise models of gene switches

Molecular noise in gene regulatory networks has two intrinsic components, one part being due to fluctuations caused by the birth and death of protein or mRNA molecules which are often present in small numbers and the other part arising from gene state switching, a single molecule event. Stochastic dynamics of gene regulatory circuits appears to be largely responsible for bifurcations into a set of multi-attractor states that encode different cell phenotypes. The interplay of dichotomous single molecule gene noise with the nonlinear architecture of genetic networks generates rich and complex phenomena. In this paper, we elaborate on an approximate framework that leads to simple hybrid multi-scale schemes well suited for the quantitative exploration of the steady state properties of large-scale cellular genetic circuits. Through a path sum based analysis of trajectory statistics, we elucidate the connection of these hybrid schemes to the underlying master equation and provide a rigorous justification for using dichotomous noise based models to study genetic networks. Numerical simulations of circuit models reveal that the contribution of the genetic noise of single molecule origin to the total noise is significant for a wide range of kinetic regimes.

Resolving the NFκB Heterodimer Binding Paradox: Strain and Frustration Guide the Binding of Dimeric Transcription Factors

Many eukaryotic transcription factors function after forming oligomers. The choice of protein partners is a nonrandom event that has distinct functional consequences for gene regulation. In the present work we examine three dimers of transcription factors in the NFκB family: p50p50, p50p65, and p65p65. The NFκB dimers bind to a myriad of genomic sites and switch the targeted genes on or off with precision. The p65p50 heterodimer of NFκB is the strongest DNA binder, and its unbinding is controlled kinetically by molecular stripping from the DNA induced by IκB. In contrast, the homodimeric forms of NFκB, p50p50 and p65p65, bind DNA with significantly less affinity, which places the DNA residence of the homodimers under thermodynamic rather than kinetic control. It seems paradoxical that the heterodimer should bind more strongly than either of the symmetric homodimers since DNA is a nearly symmetric target. Using a variety of energy landscape analysis tools, here we uncover the features in the molecular architecture of NFκB dimers that are responsible for these drastically different binding free energies. We show that frustration in the heterodimer interface gives the heterodimer greater conformational plasticity, allowing the heterodimer to better accommodate the DNA. We also show how the elastic energy and mechanical strain in NFκB dimers can be found by extracting the principal components of the fluctuations in Cartesian coordinates as well as fluctuations in the space of physical contacts, which are sampled via simulations with a predictive energy landscape Hamiltonian. These energetic contributions determine the specific detailed mechanisms of binding and stripping for both homo- and heterodimers.