Structure and Flexibility Adaptation in Nonspecific and Specific Protein-DNA Complexes

Accuracy and Reproducibility of Lipari-Szabo Order Parameters From Molecular Dynamics

The Lipari-Szabo generalized order parameter probes the picosecond to nanosecond time scale motions of a protein and is useful for rationalizing a multitude of biological processes such as protein recognition and ligand binding. Although these fast motions are an important and intrinsic property of proteins, it remains unclear what simulation conditions are most suitable to reproduce methyl symmetry axis side chain order parameter data (Saxis2) from molecular dynamics simulations. In this study, we show that, while Saxis2 tends to converge within tens of nanoseconds, it is essential to run 10 to 20 replicas starting from configurations close to the experimental structure to obtain the best agreement with experimental Saxis2 values. Additionally, in a comparison of force fields, AMBER ff14SB outperforms CHARMM36m in accurately capturing these fast time scale motions, and we suggest that the origin of this performance gap is likely attributed to differences in side chain torsional parametrization and not due to differences in the global protein conformations sampled by the force fields. This study provides insight into obtaining accurate and reproducible Saxis2 values from molecular simulations and underscores the necessity of using replica simulations to compute equilibrium properties.

How to scan naked DNA using promiscuous recognition and no clamping: a model for pioneer transcription factors

Most DNA scanning proteins uniquely recognize their cognate sequence motif and slide on DNA assisted by some sort of clamping interface. The pioneer transcription factors that control cell fate in eukaryotes must forgo both elements to gain access to DNA in naked and chromatin forms; thus, whether or how these factors scan naked DNA is unknown. Here, we use single-molecule techniques to investigate naked DNA scanning by the Engrailed homeodomain (enHD) as paradigm of highly promiscuous recognition and open DNA binding interface. We find that enHD scans naked DNA quite effectively, and about 200000-fold faster than expected for a continuous promiscuous slide. To do so, enHD scans about 675 bp of DNA in 100 ms and then redeploys stochastically to another location 530 bp afar in just 10 ms. During the scanning phase enHD alternates between slow- and medium-paced modes every 3 and 40 ms, respectively. We also find that enHD binds nucleosomes and does so with enhanced affinity relative to naked DNA. Our results demonstrate that pioneer-like transcription factors can in principle do both, target nucleosomes and scan active DNA efficiently. The hybrid scanning mechanism used by enHD appears particularly well suited for the highly complex genomic signals of eukaryotic cells.

Viral DNA polymerase structures reveal mechanisms of antiviral drug resistance

DNA polymerases are important drug targets, and many structural studies have captured them in distinct conformations. However, a detailed understanding of the impact of polymerase conformational dynamics on drug resistance is lacking. We determined cryoelectron microscopy (cryo-EM) structures of DNA-bound herpes simplex virus polymerase holoenzyme in multiple conformations and interacting with antivirals in clinical use. These structures reveal how the catalytic subunit Pol and the processivity factor UL42 bind DNA to promote processive DNA synthesis. Unexpectedly, in the absence of an incoming nucleotide, we observed Pol in multiple conformations with the closed state sampled by the fingers domain. Drug-bound structures reveal how antivirals may selectively bind enzymes that more readily adopt the closed conformation. Molecular dynamics simulations and the cryo-EM structure of a drug-resistant mutant indicate that some resistance mutations modulate conformational dynamics rather than directly impacting drug binding, thus clarifying mechanisms that drive drug selectivity.

Role of Shape Deformation of DNA-Binding Sites in Regulating the Efficiency and Specificity in Their Recognition by DNA-Binding Proteins

Accurate transcription of genetic information is crucial, involving precise recognition of the binding motifs by DNA-binding proteins. While some proteins rely on short-range hydrophobic and hydrogen bonding interactions at binding sites, others employ a DNA shape readout mechanism for specific recognition. In this mechanism, variations in DNA shape at the binding motif resulted from either inherent flexibility or binding of proteins at adjacent sites are sensed and capitalized by the searching proteins to locate them specifically. Through extensive computer simulations, we investigate both scenarios to uncover the underlying mechanism and origin of specificity in the DNA shape readout mechanism. Our findings reveal that deformation in shape at the binding motif creates an entropy funnel, allowing information about altered shapes to manifest as fluctuations in minor groove widths. This signal enhances the efficiency of nonspecific search of nearby proteins by directing their movement toward the binding site, primarily driven by a gain in entropy. We propose this as a generic mechanism for DNA shape readout, where specificity arises from the alignment between the molecular frustration of the searching protein and the ruggedness of the entropic funnel governed by molecular features of the protein and arrangement of the DNA bases at the binding site, respectively.

Analyzing paramagnetic NMR data on target DNA search by proteins using a discrete-state kinetic model for translocation

Before reaching their targets, sequence-specific DNA-binding proteins nonspecifically bind to DNA through electrostatic interactions and stochastically change their locations on DNA. Investigations into the dynamics of DNA-scanning by proteins are nontrivial due to the simultaneous presence of multiple translocation mechanisms and many sites for the protein to nonspecifically bind to DNA. Nuclear magnetic resonance (NMR) spectroscopy can provide information about the target DNA search processes at an atomic level. Paramagnetic relaxation enhancement (PRE) is particularly useful to study how the proteins scan DNA in the search process. Previously, relatively simple two-state or three-state exchange models were used to explain PRE data reflecting the target search process. In this work, using more realistic discrete-state stochastic kinetics models embedded into an NMR master equation, we analyzed the PRE data for the HoxD9 homeodomain interacting with DNA. The kinetic models that incorporate sliding, dissociation, association, and intersegment transfer can reproduce the PRE profiles observed at some different ionic strengths. The analysis confirms the previous interpretation of the PRE data and shows that the protein's probability distribution among nonspecific sites is nonuniform during the target DNA search process.

DNA-Nanostructure-Guided Assembly of Proteins into Programmable Shapes

The development of methods to synthesize artificial protein complexes with precisely controlled configurations will enable diverse biological and medical applications. Using DNA to link proteins provides programmability that can be difficult to achieve with other methods. Here, we use DNA origami as an "assembler" to guide the linking of protein-DNA conjugates using a series of oligonucleotide hybridization and displacement operations. We constructed several isomeric protein nanostructures, including a dimer, two types of trimer structures, and three types of tetramer assemblies, on a DNA origami platform by using a C3-symmetric building block composed of a protein trimer modified with DNA handles. Our approach expands the scope for the precise assembly of protein-based nanostructures and will enable the formulation of functional protein complexes with stoichiometric and geometric control.

Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes

The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.

Lambda CI Binding to Related Phage Operator Sequences Validates Alignment Algorithm and Highlights the Importance of Overlooked Bonds

Bacteriophage λ's CI repressor protein controls a genetic switch between the virus's lysogenic and lytic lifecycles, in part, by selectively binding to six different DNA sequences within the phage genome-collectively referred to as operator sites. However, the minimal level of information needed for CI to recognize and specifically bind these six unique-but-related sequences is unclear. In a previous study, we introduced an algorithm that extracts the minimal direct readout information needed for λ-CI to recognize and bind its six binding sites. We further revealed direct readout information shared among three evolutionarily related lambdoid phages: λ-phage, Enterobacteria phage VT2-Sakai, and Stx2 converting phage I, suggesting that the λ-CI protein could bind to the operator sites of these other phages. In this study, we show that λ-CI can indeed bind the other two phages' cognate binding sites as predicted using our algorithm, validating the hypotheses from that paper. We go on to demonstrate the importance of specific hydrogen bond donors and acceptors that are maintained despite changes to the nucleobase itself, and another that has an important role in recognition and binding. This in vitro validation of our algorithm supports its use as a tool to predict alternative binding sites for DNA-binding proteins.

Studying Nucleoid-Associated Protein-DNA Interactions Using Polymer Microgels as Synthetic Mimics

Microfluidically fabricated polymer microgels are used as an experimental platform to analyze protein-DNA interactions regulating bacterial cell division. In particular, we focused on the nucleoid-associated protein SlmA, which forms a nucleoprotein complex with short DNA binding sequences (SBS) that acts as a negative regulator of the division ring stability in Escherichia coli. To mimic the bacterial nucleoid as a dense DNA region of a bacterial cell and investigate the influence of charge and permeability on protein binding and diffusion in there, we have chosen nonionic polyethylene glycol and anionic hyaluronic acid as precursor materials for hydrogel formation, previously functionalized with SBS. SlmA binds specifically to the coupled SBS for both types of microgels while preferentially accumulating at the microgels' surface. We could control the binding specificity by adjusting the buffer composition of the DNA-functionalized microgels. The microgel charge did not impact protein binding; however, hyaluronic acid-based microgels exhibit a higher permeability, promoting protein diffusion; thus, they were the preferred choice for preparing nucleoid mimics. The approaches described here provide attractive tools for bottom-up reconstitution of essential cellular processes in media that more faithfully reproduce intracellular environments.

Multiple modes of DNA compaction by protamine

In sperm cells, protamine replaces histones to compact DNA 10-20 times more than in somatic cells. To characterize the extreme compaction, we employed confocal microscopy and optical tweezers to determine the conformations and stability of protamine-bound λ-DNA. Confocal images show increasing compaction of λ-DNA at increasing protamine concentration. In the presence of protamine, single λ-DNA molecules form bends and loops that unravel at 10-40 pN forces as well as coils that shorten the contour length by up to 40% and withstand forces strong enough (~55 pN) for strand separation. Strand separation nucleates coils, indicating protamine insertion into DNA bases. Protamine may participate in both local and higher-order chromatin organization, leading to extreme compaction and global transcription silencing.

Genetic switching by the Lac repressor is based on two-state Monod-Wyman-Changeux allostery

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod-Wyman-Changeux model.

Nonspecific vs. specific DNA binding free energetics of a transcription factor domain protein

Transcription factor (TF) proteins regulate gene expression by binding to specific sites on the genome. In the facilitated diffusion model, an optimized search process is achieved by the TF alternating between 3D diffusion in the bulk and 1D diffusion along DNA. While undergoing 1D diffusion, the protein can switch from a search mode for fast diffusion along nonspecific DNA to a recognition mode for stable binding to specific DNA. It was recently noticed that, for a small TF domain protein, reorientations on DNA happen between the nonspecific and specific DNA binding. We here conducted all-atom molecular dynamics simulations with steering forces to reveal the protein-DNA binding free energetics, confirming that the search and recognition modes are distinguished primarily by protein orientations on the DNA. As the binding free energy difference between the specific and nonspecific DNA system slightly deviates from that being estimated directly from dissociation constants on 15-bp DNA constructs, we hypothesize that the discrepancy can come from DNA sequences flanking the 6-bp central binding sites that impact on the dissociation kinetics measurements. The hypothesis is supported by a simplified spherical protein-DNA model along with stochastic simulations and kinetic modeling.

Short tandem repeats bind transcription factors to tune eukaryotic gene expression

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

Can molecular dynamics be used to simulate biomolecular recognition?

There are many problems in biochemistry that are difficult to study experimentally. Simulation methods are appealing due to direct availability of atomic coordinates as a function of time. However, direct molecular simulations are challenged by the size of systems and the time scales needed to describe relevant motions. In theory, enhanced sampling algorithms can help to overcome some of the limitations of molecular simulations. Here, we discuss a problem in biochemistry that offers a significant challenge for enhanced sampling methods and that could, therefore, serve as a benchmark for comparing approaches that use machine learning to find suitable collective variables. In particular, we study the transitions LacI undergoes upon moving between being non-specifically and specifically bound to DNA. Many degrees of freedom change during this transition and that the transition does not occur reversibly in simulations if only a subset of these degrees of freedom are biased. We also explain why this problem is so important to biologists and the transformative impact that a simulation of it would have on the understanding of DNA regulation.

Targeting MYC with modular synthetic transcriptional repressors derived from bHLH DNA-binding domains

Despite unequivocal roles in disease, transcription factors (TFs) remain largely untapped as pharmacologic targets due to the challenges in targeting protein-protein and protein-DNA interactions. Here we report a chemical strategy to generate modular synthetic transcriptional repressors (STRs) derived from the bHLH domain of MAX. Our synthetic approach yields chemically stabilized tertiary domain mimetics that cooperatively bind the MYC/MAX consensus E-box motif with nanomolar affinity, exhibit specificity that is equivalent to or beyond that of full-length TFs and directly compete with MYC/MAX protein for DNA binding. A lead STR directly inhibits MYC binding in cells, downregulates MYC-dependent expression programs at the proteome level and inhibits MYC-dependent cell proliferation. Co-crystallization and structure determination of a STR:E-box DNA complex confirms retention of DNA recognition in a near identical manner as full-length bHLH TFs. We additionally demonstrate structure-blind design of STRs derived from alternative bHLH-TFs, confirming that STRs can be used to develop highly specific mimetics of TFs targeting other gene regulatory elements.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Ligand-specific changes in conformational flexibility mediate long-range allostery in the lac repressor

Biological regulation ubiquitously depends on protein allostery, but the regulatory mechanisms are incompletely understood, especially in proteins that undergo ligand-induced allostery with few structural changes. Here we used hydrogen-deuterium exchange with mass spectrometry (HDX/MS) to map allosteric effects in a paradigm ligand-responsive transcription factor, the lac repressor (LacI), in different functional states (apo, or bound to inducer, anti-inducer, and/or DNA). Although X-ray crystal structures of the LacI core domain in these states are nearly indistinguishable, HDX/MS experiments reveal widespread differences in flexibility. We integrate these results with modeling of protein-ligand-solvent interactions to propose a revised model for allostery in LacI, where ligand binding allosterically shifts the conformational ensemble as a result of distinct changes in the rigidity of secondary structures in the different states. Our model provides a mechanistic basis for the altered function of distal mutations. More generally, our approach provides a platform for characterizing and engineering protein allostery.

Evolutionary conservation of protein dynamics: insights from all-atom molecular dynamics simulations of 'peptidase' domain of Spt16

Protein function is encoded in its sequence, manifested in its three-dimensional structure, and facilitated by its dynamics. Studies have suggested that protein structures with higher sequence similarity could have more similar patterns of dynamics. However, such studies of protein dynamics within and across protein families typically rely on coarse-grained models, or approximate metrics like crystallographic B-factors. This study uses µs scale molecular dynamics (MD) simulations to explore the conservation of dynamics among homologs of ∼50 kDa N-terminal module of Spt16 (Spt16N). Spt16N from Saccharomyces cerevisiae (Sc-Spt16N) and three of its homologs with 30-40% sequence identities were available in the PDB. To make our data-set more comprehensive, the crystal structure of an additional homolog (62% sequence identity with Sc-Spt16N) was solved at 1.7 Å resolution. Cumulative MD simulations of 6 µs were carried out on these Spt16N structures and on two additional protein structures with varying degrees of similarity to it. The simulations revealed that correlation in patterns of backbone fluctuations vary linearly with sequence identity. This trend could not be inferred using crystallographic B-factors. Further, normal mode analysis suggested a similar pattern of inter-domain (inter-lobe) motions not only among Spt16N homologs, but also in the M24 peptidase structure. On the other hand, MD simulation results highlighted conserved motions that were found unique for Spt16N protein, this along with electrostatics trends shed light on functional aspects of Spt16N.Communicated by Ramaswamy H. Sarma.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

Alignment of major-groove hydrogen bond arrays uncovers shared information between different DNA sequences that bind the same protein

Protein-DNA binding is of a great interest due to its importance in many biological processes. Previous studies have presented many factors responsible for the recognition and specificity, but understanding the minimal informational requirements for proteins that bind to multiple DNA-sites is still an understudied area of bioinformatics. Here we focus on the hydrogen bonds displayed by the target DNA in the major groove that take part in protein-binding. We show that analyses focused on the base pair identity may overlook key hydrogen bonds. We have developed an algorithm that converts a nucleotide sequence into an array of hydrogen bond donors and acceptors and methyl groups. It then aligns these non-covalent interaction arrays to identify what information is being maintained among multiple DNA sequences. For three different DNA-binding proteins, Lactose repressor, controller protein and λ-CI repressor, we uncovered the minimal pattern of hydrogen bonds that are common amongst all the binding sequences. Notably in the three proteins, key interacting hydrogen bonds are maintained despite nucleobase mutations in the corresponding binding sites. We believe this work will be useful for developing new DNA binding proteins and shed new light on evolutionary relationships.

Flexibility of flanking DNA is a key determinant of transcription factor affinity for the core motif

Selective gene regulation is mediated by recognition of specific DNA sequences by transcription factors (TFs). The extremely challenging task of searching out specific cognate DNA binding sites among several million putative sites within the eukaryotic genome is achieved by complex molecular recognition mechanisms. Elements of this recognition code include the core binding sequence, the flanking sequence context, and the shape and conformational flexibility of the composite binding site. To unravel the extent to which DNA flexibility modulates TF binding, in this study, we employed experimentally guided molecular dynamics simulations of ternary complex of closely related Hox heterodimers Exd-Ubx and Exd-Scr with DNA. Results demonstrate that flexibility signatures embedded in the flanking sequences impact TF binding at the cognate binding site. A DNA sequence has intrinsic shape and flexibility features. While shape features are localized, our analyses reveal that flexibility features of the flanking sequences percolate several basepairs and allosterically modulate TF binding at the core. We also show that lack of flexibility in the motif context can render the cognate site resistant to protein-induced shape changes and subsequently lower TF binding affinity. Overall, this study suggests that flexibility-guided DNA shape, and not merely the static shape, is a key unexplored component of the complex DNA-TF recognition code.

"Flexible hinge" dynamics in mismatched DNA revealed by fluorescence correlation spectroscopy

Altered unwinding/bending fluctuations at DNA lesion sites are implicated as plausible mechanisms for damage sensing by DNA-repair proteins. These dynamics are expected to occur on similar timescales as one-dimensional (1D) diffusion of proteins on DNA if effective in stalling these proteins as they scan DNA. We examined the flexibility and dynamics of DNA oligomers containing 3 base pair (bp) mismatched sites specifically recognized in vitro by nucleotide excision repair protein Rad4 (yeast ortholog of mammalian XPC). A previous Forster resonance energy transfer (FRET) study mapped DNA conformational distributions with cytosine analog FRET pair primarily sensitive to DNA twisting/unwinding deformations (Chakraborty et al. Nucleic Acids Res. 46: 1240-1255 (2018)). These studies revealed B-DNA conformations for nonspecific (matched) constructs but significant unwinding for mismatched constructs specifically recognized by Rad4, even in the absence of Rad4. The timescales of these unwinding fluctuations, however, remained elusive. Here, we labeled DNA with Atto550/Atto647N FRET dyes suitable for fluorescence correlation spectroscopy (FCS). With these probes, we detected higher FRET in specific, mismatched DNA compared with matched DNA, reaffirming unwinding/bending deformations in mismatched DNA. FCS unveiled the dynamics of these spontaneous deformations at ~ 300 µs with no fluctuations detected for matched DNA within the ~ 600 ns-10 ms FCS time window. These studies are the first to visualize anomalous unwinding/bending fluctuations in mismatched DNA on timescales that overlap with the < 500 µs "stepping" times of repair proteins on DNA. Such "flexible hinge" dynamics at lesion sites could arrest a diffusing protein to facilitate damage interrogation and recognition.

Repetitive DNA symmetry elements negatively regulate gene expression in embryonic stem cells

Transcription factor (TF) binding to genomic DNA elements constitutes one of the key mechanisms that regulates gene expression program in cells. Both consensus and nonconsensus DNA sequence elements influence the recognition specificity of TFs. Based on the analysis of experimentally determined c-Myc binding preferences to genomic DNA, here we statistically predict that certain repetitive, nonconsensus DNA symmetry elements can relatively reduce TF-DNA binding preferences. This is in contrast to a different set of repetitive, nonconsensus symmetry elements that can increase the strength of TF-DNA binding. Using c-Myc enhancer reporter system containing consensus motif flanked by nonconsensus sequences in embryonic stem cells, we directly demonstrate that the enrichment in such negatively regulating repetitive symmetry elements is sufficient to reduce the gene expression level compared with native genomic sequences. Negatively regulating repetitive symmetry elements around consensus c-Myc motif and DNA sequences containing consensus c-Myc motif flanked by entirely randomized sequences show similar expression baseline. A possible explanation for this observation is that rather than complete repression, negatively regulating repetitive symmetry elements play a regulatory role in fine-tuning the reduction of gene expression, most probably by binding TFs other than c-Myc.

Thermodynamic analysis of Zα domain-nucleic acid interactions

DNA/RNA molecules adopting the left-handed conformation (Z-form) have been attributed with immunogenic properties. However, their biological role and importance has been a topic of debate for many years. The discovery of Z-DNA/RNA binding domains (Zα domains) in varied proteins that are involved in the innate immune response, such as the interferon inducible form of the RNA editing enzyme ADAR1 (p150), Z-DNA binding protein 1 (ZBP1), the fish kinase PKZ and the poxvirus inhibitor of interferon response E3L, indicates important roles of Z-DNA/RNA in immunity and self/non-self-discrimination. Such Zα domain-containing proteins recognise left-handed Z-DNA/RNA in a conformation-specific manner. Recent studies have implicated these domains in virus recognition. Given these important emerging roles for the Zα domains, it is pivotal to understand the mechanism of recognition of the Z-DNA/Z-RNA by these domains. To this end, we assessed the binding thermodynamics of Zα domain from ORF112 and ADAR1 on T(CG)3 and T(CG)6 oligonucleotides which have high propensity to adopt the Z-conformation. Our study highlights important differences in the mode of oligonucleotide binding by the two Zα domains originating from different proteins. Site-directed mutagenesis was employed together with isothermal titration calorimetry to tease apart finer details of the binding thermodynamics. Our work advances the understanding on binding thermodynamics of Zα domains to their cognate nucleic acid substrates and paves the ground for future efforts to gain a complete appreciation of this process.

Lighting up Nobel Prize-winning studies with protein intrinsic disorder

Intrinsically disordered proteins and regions (IDPs and IDRs) and their importance in biology are becoming increasingly recognized in biology, biochemistry, molecular biology and chemistry textbooks, as well as in current protein science and structural biology curricula. We argue that the sequence → dynamic conformational ensemble → function principle is of equal importance as the classical sequence → structure → function paradigm. To highlight this point, we describe the IDPs and/or IDRs behind the discoveries associated with 17 Nobel Prizes, 11 in Physiology or Medicine and 6 in Chemistry. The Nobel Laureates themselves did not always mention that the proteins underlying the phenomena investigated in their award-winning studies are in fact IDPs or contain IDRs. In several cases, IDP- or IDR-based molecular functions have been elucidated, while in other instances, it is recognized that the respective protein(s) contain IDRs, but the specific IDR-based molecular functions have yet to be determined. To highlight the importance of IDPs and IDRs as general principle in biology, we present here illustrative examples of IDPs/IDRs in Nobel Prize-winning mechanisms and processes.

Structural and dynamic studies of DNA recognition by NF-κB p50 RHR homodimer using methyl NMR spectroscopy

Protein dynamics involving higher-energy sparsely populated conformational substates are frequently critical for protein function. This study describes the dynamics of the homodimer (p50)2 of the p50 Rel homology region (RHR) of the transcription factor NF-κB, using 13C relaxation dispersion experiments with specifically (13C, 1H)-labeled methyl groups of Ile (δ), Leu and Val. Free (p50)2 is highly dynamic in solution, showing μs-ms relaxation dispersion consistent with exchange between the ground state and higher energy substates. These fluctuations propagate from the DNA-binding loops through the core of the domain. The motions are damped in the presence of κB DNA, but the NMR spectra of the DNA complexes reveal multiple local conformations of the p50 RHR homodimer bound to certain κB DNA sequences. Varying the length and sequence of κB DNA revealed two factors that promote a single bound conformation for the complex: the length of the κB site in the duplex and a symmetrical sequence of guanine nucleotides at both ends of the recognition motif. The dynamic nature of the DNA-binding loops, together with the multiple bound conformations of p50 RHR with certain κB sites, is consistent with variations in the transcriptional activity of the p50 homodimer with different κB sequences.

Function and regulation of transcription factors during mitosis-to-G1 transition

During cell division, drastic cellular changes characteristic of mitosis result in the inactivation of the transcriptional machinery, and global downregulation of transcription. Sequence-specific transcription factors (TFs) have thus been considered mere bystanders, devoid of any regulatory function during mitosis. This view changed significantly in recent years, upon the conclusion that many TFs associate with condensed chromosomes during cell division, even occupying a fraction of their genomic target sites in mitotic chromatin. This finding was at the origin of the concept of mitotic bookmarking by TFs, proposed as a mechanism to propagate gene regulatory information across cell divisions, by facilitating the reactivation of specific bookmarked genes. While the underlying mechanisms and biological significance of this model remain elusive, recent developments in this fast-moving field have cast new light into TF activity during mitosis, beyond a bookmarking role. Here, we start by reviewing the most recent findings on the complex nature of TF-chromatin interactions during mitosis, and on mechanisms that may regulate them. Next, and in light of recent reports describing how transcription is reinitiated in temporally distinct waves during mitosis-to-G1 transition, we explore how TFs may contribute to defining this hierarchical gene expression process. Finally, we discuss how TF activity during mitotic exit may impact the acquisition of cell identity upon cell division, and propose a model that integrates dynamic changes in TF-chromatin interactions during this cell-cycle period, with the execution of cell-fate decisions.

Cognate DNA Recognition by Engrailed Homeodomain Involves a Conformational Change Controlled via an Electrostatic-Spring-Loaded Latch

Transcription factors must scan genomic DNA, recognize the cognate sequence of their control element(s), and bind tightly to them. The DNA recognition process is primarily carried out by their DNA binding domains (DBD), which interact with the cognate site with high affinity and more weakly with any other DNA sequence. DBDs are generally thought to bind to their cognate DNA without changing conformation (lock-and-key). Here, we used nuclear magnetic resonance and circular dichroism to investigate the interplay between DNA recognition and DBD conformation in the engrailed homeodomain (enHD), as a model case for the homeodomain family of eukaryotic DBDs. We found that the conformational ensemble of enHD is rather flexible and becomes gradually more disordered as ionic strength decreases following a Debye–Hückel’s dependence. Our analysis indicates that enHD’s response to ionic strength is mediated by a built-in electrostatic spring-loaded latch that operates as a conformational transducer. We also found that, at moderate ionic strengths, enHD changes conformation upon binding to cognate DNA. This change is of larger amplitude and somewhat orthogonal to the response to ionic strength. As a consequence, very high ionic strengths (e.g., 700 mM) block the electrostatic-spring-loaded latch and binding to cognate DNA becomes lock-and-key. However, the interplay between enHD conformation and cognate DNA binding is robust across a range of ionic strengths (i.e., 45 to 300 mM) that covers the physiologically-relevant conditions. Therefore, our results demonstrate the presence of a mechanism for the conformational control of cognate DNA recognition on a eukaryotic DBD. This mechanism can function as a signal transducer that locks the DBD in place upon encountering the cognate site during active DNA scanning. The electrostatic-spring-loaded latch of enHD can also enable the fine control of DNA recognition in response to transient changes in local ionic strength induced by variate physiological processes.

The Non-continuum Nature of Eukaryotic Transcriptional Regulation

Eukaryotic transcription factors are versatile mediators of specificity in gene regulation. This versatility is achieved through mutual specification by context-specific DNA binding on the one hand, and identity-specific protein-protein partnerships on the other. This interactivity, known as combinatorial control, enables a repertoire of complex transcriptional outputs that are qualitatively disjoint, or non-continuum, with respect to binding affinity. This feature contrasts starkly with prokaryotic gene regulators, whose activities in general vary quantitatively in step with binding affinity. Biophysical studies on prokaryotic model systems and more recent investigations on transcription factors highlight an important role for folded state dynamics and molecular hydration in protein/DNA recognition. Analysis of molecular models of combinatorial control and recent literature in low-affinity gene regulation suggest that transcription factors harbor unique conformational dynamics that are inaccessible or unused by prokaryotic DNA-binding proteins. Thus, understanding the intrinsic dynamics involved in DNA binding and co-regulator recruitment appears to be a key to understanding how transcription factors mediate non-continuum outcomes in eukaryotic gene expression, and how such capability might have evolved from ancient, structurally conserved counterparts.

Herpesvirus DNA polymerase: Structures, functions, and mechanisms

Herpesviruses comprise a family of DNA viruses that cause a variety of human and veterinary diseases. During productive infection, mammalian, avian, and reptilian herpesviruses replicate their genomes using a set of conserved viral proteins that include a two subunit DNA polymerase. This enzyme is both a model system for family B DNA polymerases and a target for inhibition by antiviral drugs. This chapter reviews the structure, function, and mechanisms of the polymerase of herpes simplex viruses 1 and 2 (HSV), with only occasional mention of polymerases of other herpesviruses such as human cytomegalovirus (HCMV). Antiviral polymerase inhibitors have had the most success against HSV and HCMV. Detailed structural information regarding HSV DNA polymerase is available, as is much functional information regarding the activities of the catalytic subunit (Pol), which include a DNA polymerization activity that can utilize both DNA and RNA primers, a 3'-5' exonuclease activity, and other activities in DNA synthesis and repair and in pathogenesis, including some remaining to be biochemically defined. Similarly, much is known regarding the accessory subunit, which both resembles and differs from sliding clamp processivity factors such as PCNA, and the interactions of this subunit with Pol and DNA. Both subunits contribute to replication fidelity (or lack thereof). The availability of both pharmacologic and genetic tools not only enabled the initial identification of Pol and the pol gene, but has also helped dissect their functions. Nevertheless, important questions remain for this long-studied enzyme, which is still an attractive target for new drug discovery.

How DNA affects the hyperthermophilic protein Ape10b2 for oligomerization: an investigation using multiple short molecular dynamics simulations

Alba2 is a hyperthermophilic DNA-binding protein, and DNA plays a crucial role in the Alba2 oligomerization process. It is a pity that there is limited research in terms of how DNA affects the conformational change of Alba2 in oligomerization. Herein, we complement the crystal structure of the Ape10b2 (belongs to Alba2)-dsDNA complex (PDB ID: 3U6Y) and employ multiple short molecular dynamics (MSMD) simulations to illuminate the influence of DNA on Ape10b2 at four temperatures (300, 343, 363, and 373 K). Our results indicate that DNA could cause the conformational changes of two important regions (loop1 and loop5), which may be beneficial for protein oligomerization. The results of hydrogen bond analysis show that the increasing number of hydrogen bonds between two monomers of Ape10b2 may also be a favorable factor for oligomerization. In addition, Ape10b2 can stabilize DNA by electrostatic interactions with an increase in temperature, and five residues (Arg40, Arg42, Asn43, Asn45, and Arg46) play a stabilizing role during protein binding to DNA. Our findings could help in understanding the favorable factors leading to protein oligomerization, which contributes to enzyme engineering research from an industrial perspective.

Probing allosteric regulations with coevolution-driven molecular simulations

Protein-mediated allosteric regulations are essential in biology, but their quantitative characterization continues to posit formidable challenges for both experiments and computations. Here, we combine coevolutionary information, multiscale molecular simulations, and free-energy methods to interrogate and quantify the allosteric regulation of functional changes in protein complexes. We apply this approach to investigate the regulation of adenylyl cyclase (AC) by stimulatory and inhibitory G proteins—a prototypical allosteric system that has long escaped from in-depth molecular characterization. We reveal a surprisingly simple ON/OFF regulation of AC functional dynamics through multiple pathways of information transfer. The binding of G proteins reshapes the free-energy landscape of AC following the classical population-shift paradigm. The model agrees with structural and biochemical data and reveals previously unknown experimentally consistent intermediates. Our approach showcases a general strategy to explore uncharted functional space in complex biomolecular regulations.

Exploring the molecular basis of UG-rich RNA recognition by the human splicing factor TDP-43 using molecular dynamics simulation and free energy calculation

Transactivation response element RNA/DNA-binding protein 43 (TDP-43) is involved in the regulation of alternative splicing of human neurodegenerative disease-related genes through binding to long UG-rich RNA sequences. Mutations in TDP-43, most in the homeodomain, cause neurological disorders such as amyotrophic lateral sclerosis and fronto temporal lobar degeneration. Several mutants destabilize the structure and disrupt RNA-binding activity. The biological functions of these mutants have been characterized, but the structural basis behind the loss of RNA-binding activity is unclear. Focused on the specific TDP-43-ssRNA complex (PDB code 4BS2), we applied molecular dynamics simulations and the molecular mechanics Poisson-Boltzmann surface area free energy calculation to characterize and explore the structural and dynamic effects between ssRNA and TDP-43. The energetic analysis indicated that the intermolecular van der Waals interaction and nonpolar solvation energy play an important role in the binding process of TDP-43 and ssRNA. Compared with the wild-type TDP-43, the reduction of the polar or non-polar interaction between all the mutants F149A, D105A/S254A, R171A/D174A, F147L/F149L/F229L/F231L and ssRNA is the main reason for the reduction of its binding free energy. Decomposing energies suggested that the extensive interactions between TDP-43 and the nitrogenous bases of ssRNA are responsible for the specific ssRNA recognition by TDP-43. These results elucidated the TDP-43-ssRNA interaction comprehensively and further extended our understanding of the previous experimental data. The uncovering of TDP-43-ssRNA recognition mechanism will provide us useful insights and new chances for the development of anti-neurodegenerative drugs.

Structural variation affecting DNA backbone interactions underlies adaptation of B3 DNA binding domains to constraints imposed by protein architecture

Functional and architectural diversification of transcription factor families has played a central role in the independent evolution of complex development in plants and animals. Here, we investigate the role of architectural constraints on evolution of B3 DNA binding domains that regulate plant embryogenesis. B3 domains of ABI3, FUS3, LEC2 and VAL1 proteins recognize the same cis-element. Complex architectures of ABI3 and VAL1 integrate cis-element recognition with other signals, whereas LEC2 and FUS3 have reduced architectures conducive to roles as pioneer activators. In yeast and plant in vivo assays, B3 domain functions correlate with architectural complexity of the parent transcription factor rather than phylogenetic relatedness. In a complex architecture, attenuated ABI3-B3 and VAL1-B3 activities enable integration of cis-element recognition with hormone signaling, whereas hyper-active LEC2-B3 and FUS3-B3 over-ride hormonal control. Three clade-specific amino acid substitutions (β4-triad) implicated in interactions with the DNA backbone account for divergence of LEC2-B3 and ABI3-B3. We find a striking correlation between differences in in vitro DNA binding affinity and in vivo activities of B3 domains in plants and yeast. Our results highlight the role of DNA backbone interactions that preserve DNA sequence specificity in adaptation of B3 domains to functional constraints associated with domain architecture.

Hierarchical reactivation of transcription during mitosis-to-G1 transition by Brn2 and Ascl1 in neural stem cells

During mitosis, chromatin condensation is accompanied by a global arrest of transcription. Recent studies suggest transcriptional reactivation upon mitotic exit occurs in temporally coordinated waves, but the underlying regulatory principles have yet to be elucidated. In particular, the contribution of sequence-specific transcription factors (TFs) remains poorly understood. Here we report that Brn2, an important regulator of neural stem cell identity, associates with condensed chromatin throughout cell division, as assessed by live-cell imaging of proliferating neural stem cells. In contrast, the neuronal fate determinant Ascl1 dissociates from mitotic chromosomes. ChIP-seq analysis reveals that Brn2 mitotic chromosome binding does not result in sequence-specific interactions prior to mitotic exit, relying mostly on electrostatic forces. Nevertheless, surveying active transcription using single-molecule RNA-FISH against immature transcripts reveals differential reactivation kinetics for key targets of Brn2 and Ascl1, with transcription onset detected in early (anaphase) versus late (early G1) phases, respectively. Moreover, by using a mitotic-specific dominant-negative approach, we show that competing with Brn2 binding during mitotic exit reduces the transcription of its target gene Nestin Our study shows an important role for differential binding of TFs to mitotic chromosomes, governed by their electrostatic properties, in defining the temporal order of transcriptional reactivation during mitosis-to-G1 transition.

Using delayed decoupling to attenuate residual signals in editing filters

Isotope filtering methods are instrumental in biomolecular nuclear magnetic resonance (NMR) studies as they isolate signals of chemical moieties of interest within complex molecular assemblies. However, isotope filters suppress undesired signals of isotopically enriched molecules through scalar couplings, and variations in scalar couplings lead to imperfect suppressions, as occurs for aliphatic and aromatic moieties in proteins. Here, we show that signals that have escaped traditional filters can be attenuated with mitigated sensitivity losses for the desired signals of unlabeled moieties. The method uses a shared evolution between the detection and preceding preparation period to establish non-observable antiphase coherences and eliminates them through composite pulse decoupling. We demonstrate the method by isolating signals of an unlabeled post-translational modification tethered to an isotopically enriched protein.

Substitutions at Nonconserved Rheostat Positions Modulate Function by Rewiring Long-Range, Dynamic Interactions

Amino acid substitutions at nonconserved protein positions can have noncanonical and "long-distance" outcomes on protein function. Such outcomes might arise from changes in the internal protein communication network, which is often accompanied by changes in structural flexibility. To test this, we calculated flexibilities and dynamic coupling for positions in the linker region of the lactose repressor protein. This region contains nonconserved positions for which substitutions alter DNA-binding affinity. We first chose to study 11 substitutions at position 52. In computations, substitutions showed long-range effects on flexibilities of DNA-binding positions, and the degree of flexibility change correlated with experimentally measured changes in DNA binding. Substitutions also altered dynamic coupling to DNA-binding positions in a manner that captured other experimentally determined functional changes. Next, we broadened calculations to consider the dynamic coupling between 17 linker positions and the DNA-binding domain. Experimentally, these linker positions exhibited a wide range of substitution outcomes: Four conserved positions tolerated hardly any substitutions ("toggle"), ten nonconserved positions showed progressive changes from a range of substitutions ("rheostat"), and three nonconserved positions tolerated almost all substitutions ("neutral"). In computations with wild-type lactose repressor protein, the dynamic couplings between the DNA-binding domain and these linker positions showed varied degrees of asymmetry that correlated with the observed toggle/rheostat/neutral substitution outcomes. Thus, we propose that long-range and noncanonical substitutions outcomes at nonconserved positions arise from rewiring long-range communication among functionally important positions. Such calculations might enable predictions for substitution outcomes at a range of nonconserved positions.

Introduction to a special issue of Magnetic Resonance in honour of Robert Kaptein at the occasion of his 80th birthday

This publication, in honour of Robert Kaptein's 80th birthday, contains contributions from colleagues, many of whom have worked with him, and others who admire his work and have been stimulated by his research. The contributions show current research in biomolecular NMR, spin hyperpolarisation and spin chemistry, including CIDNP (chemically induced dynamic nuclear polarisation), topics to which he has contributed enormously. His proposal of the radical pair mechanism was the birth of the field of spin chemistry, and the laser CIDNP NMR experiment on a protein was a major breakthrough in hyperpolarisation research. He set milestones for biomolecular NMR by developing computational methods for protein structure determination, including restrained molecular dynamics and 3D NMR methodology. With a lac repressor headpiece, he determined one of the first protein structures determined by NMR. His studies of the lac repressor provided the first examples of detailed studies of protein nucleic acid complexes by NMR. This deepened our understanding of protein DNA recognition and led to a molecular model for protein sliding along the DNA. Furthermore, he played a leading role in establishing the cluster of NMR large-scale facilities in Europe. This editorial gives an introduction to the publication and is followed by a biography describing his contributions to magnetic resonance.

Nonspecific Binding-Fundamental Concepts and Consequences for Biosensing Applications

Nature achieves differentiation of specific and nonspecific binding in molecular interactions through precise control of biomolecules in space and time. Artificial systems such as biosensors that rely on distinguishing specific molecular binding events in a sea of nonspecific interactions have struggled to overcome this issue. Despite the numerous technological advancements in biosensor technologies, nonspecific binding has remained a critical bottleneck due to the lack of a fundamental understanding of the phenomenon. To date, the identity, cause, and influence of nonspecific binding remain topics of debate within the scientific community. In this review, we discuss the evolution of the concept of nonspecific binding over the past five decades based upon the thermodynamic, intermolecular, and structural perspectives to provide classification frameworks for biomolecular interactions. Further, we introduce various theoretical models that predict the expected behavior of biosensors in physiologically relevant environments to calculate the theoretical detection limit and to optimize sensor performance. We conclude by discussing existing practical approaches to tackle the nonspecific binding challenge in vitro for biosensing platforms and how we can both address and harness nonspecific interactions for in vivo systems.

A Bit Stickier, a Bit Slower, a Lot Stiffer: Specific vs. Nonspecific Binding of Gal4 to DNA

Transcription factors regulate gene activity by binding specific regions of genomic DNA thanks to a subtle interplay of specific and nonspecific interactions that is challenging to quantify. Here, we exploit Reflective Phantom Interface (RPI), a label-free biosensor based on optical reflectivity, to investigate the binding of the N-terminal domain of Gal4, a well-known gene regulator, to double-stranded DNA fragments containing or not its consensus sequence. The analysis of RPI-binding curves provides interaction strength and kinetics and their dependence on temperature and ionic strength. We found that the binding of Gal4 to its cognate site is stronger, as expected, but also markedly slower. We performed a combined analysis of specific and nonspecific binding—equilibrium and kinetics—by means of a simple model based on nested potential wells and found that the free energy gap between specific and nonspecific binding is of the order of one kcal/mol only. We investigated the origin of such a small value by performing all-atom molecular dynamics simulations of Gal4–DNA interactions. We found a strong enthalpy–entropy compensation, by which the binding of Gal4 to its cognate sequence entails a DNA bending and a striking conformational freezing, which could be instrumental in the biological function of Gal4.

A perspective on the molecular simulation of DNA from structural and functional aspects

As genetic material, DNA not only carries genetic information by sequence, but also affects biological functions ranging from base modification to replication, transcription and gene regulation through its structural and dynamic properties and variations. The motion and structural properties of DNA involved in related biological processes are also multi-scale, ranging from single base flipping to local DNA deformation, TF binding, G-quadruplex and i-motif formation, TAD establishment, compartmentalization and even chromosome territory formation, just to name a few. The sequence-dependent physical properties of DNA play vital role in all these events, and thus it is interesting to examine how simple sequence information affects DNA and the formation of the chromatin structure in these different hierarchical orders. Accordingly, molecular simulations can provide atomistic details of interactions and conformational dynamics involved in different biological processes of DNA, including those inaccessible by current experimental methods. In this perspective, which is mainly based on our recent studies, we provide a brief overview of the atomistic simulations on how the hierarchical structure and dynamics of DNA can be influenced by its sequences, base modifications, environmental factors and protein binding in the context of the protein-DNA interactions, gene regulation and structural organization of chromatin. We try to connect the DNA sequence, the hierarchical structures of DNA and gene regulation.

Transient non-specific DNA binding dominates the target search of bacterial DNA-binding proteins

Despite their diverse biochemical characteristics and functions, all DNA-binding proteins share the ability to accurately locate their target sites among the vast excess of non-target DNA. Toward identifying universal mechanisms of the target search, we used single-molecule tracking of 11 diverse DNA-binding proteins in living Escherichia coli. The mobility of these proteins during the target search was dictated by DNA interactions rather than by their molecular weights. By generating cells devoid of all chromosomal DNA, we discovered that the nucleoid is not a physical barrier for protein diffusion but significantly slows the motion of DNA-binding proteins through frequent short-lived DNA interactions. The representative DNA-binding proteins (irrespective of their size, concentration, or function) spend the majority (58%–99%) of their search time bound to DNA and occupy as much as ∼30% of the chromosomal DNA at any time. Chromosome crowding likely has important implications for the function of all DNA-binding proteins.

•

Protein motion was compared between unperturbed cells and DNA-free cells

•

Protein mobility was dictated by DNA interactions rather than molecular weight

•

The nucleoid is not a physical barrier for protein diffusion

•

The proteins studied spend most (58%–99%) of their search time bound to DNA

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

Characterization of nucleosome sediments for protein interaction studies by solid-state NMR spectroscopy

Regulation of DNA-templated processes such as gene transcription and DNA repair depend on the interaction of a wide range of proteins with the nucleosome, the fundamental building block of chromatin. Both solution and solid-state NMR spectroscopy have become an attractive approach to study the dynamics and interactions of nucleosomes, despite their high molecular weight of ~ 200 kDa. For solid-state NMR (ssNMR) studies, dilute solutions of nucleosomes are converted to a dense phase by sedimentation or precipitation. Since nucleosomes are known to self-associate, these dense phases may induce extensive interactions between nucleosomes, which could interfere with protein-binding studies. Here, we characterized the packing of nucleosomes in the dense phase created by sedimentation using NMR and small-angle X-ray scattering (SAXS) experiments. We found that nucleosome sediments are gels with variable degrees of solidity, have nucleosome concentration close to that found in crystals, and are stable for weeks under high-speed magic angle spinning (MAS). Furthermore, SAXS data recorded on recovered sediments indicate that there is no pronounced long-range ordering of nucleosomes in the sediment. Finally, we show that the sedimentation approach can also be used to study low-affinity protein interactions with the nucleosome. Together, our results give new insights into the sample characteristics of nucleosome sediments for ssNMR studies and illustrate the broad applicability of sedimentation-based NMR studies.

A lattice model on the rate of in vivo site-specific DNA-protein interactions

We develop a lattice model of site-specific DNA-protein interactions under in vivo conditions where DNA is modelled as a self-avoiding random walk that is embedded in a cubic lattice box resembling the living cell. The protein molecule searches for its cognate site on DNA via a combination of three dimensional (3D) and one dimensional (1D) random walks. Hopping and intersegmental transfers occur depending on the conformational state of DNA. Results show that the search acceleration ratio (= search time in pure 3D route/search time in 3D and 1D routes) asymptotically increases towards a limiting value as the dilution factor of DNA (= volume of the cell/the volume of DNA) tends towards infinity. When the dilution ratio is low, then hopping and intersegmental transfers significantly enhance the search efficiency over pure sliding. At high dilution ratio, hopping does not enhance the search efficiency much since under such situation DNA will be in a relaxed conformation that favors only sliding. In the absence of hopping and intersegmental transfers, there exists an optimum sliding time at which the search acceleration ratio attains a maximum in line with the current theoretical results. However, existence of such optimum sliding length disappears in the presence of hopping. When the DNA is confined in a small volume inside the cell resembling a natural cell system, then there exists an optimum dilution and compression ratios (= total cell volume/volume in which DNA is confined) at which the search acceleration factor attains a maximum especially in the presence of hopping and intersegmental transfers. These optimum values are consistent with the values observed in the Escherichia coli cell system. In the absence of confinement of DNA, position of the specific binding site on the genomic DNA significantly influences the search acceleration. However, such position dependent changes in the search acceleration ratio will be nullified in the presence of hopping and intersegmental transfers especially when the DNA is confined in a small volume that is embedded in an outer cell.

Exact solution to the first-passage problem for a particle with a dichotomous diffusion coefficient

We consider the problem of first-passage time for reaching a boundary of a particle which diffuses in one dimension and is confined to the region x∈(0,L), with a diffusion coefficient that switches randomly between two states, having diffusivities that are different. Exact analytical expressions are found for the survival probability of the particle as a function of time. The survival probability has a multiexponential decay, and to characterize it, we use the average rate constant k, as well as the instantaneous rate r(t). Our approach can easily be extended to the case where the diffusion coefficient takes n different values. The model should be of interest to biological processes, in which a reactant searches for a target in a heterogeneous environment, making the diffusion coefficient a random function of time. The best example for this is a protein searching for a target site on the DNA.

Structural basis of a novel repressor, SghR, controlling Agrobacterium infection by cross-talking to plants

Agrobacterium tumefaciens infects various plants and causes crown gall diseases involving temporal expression of virulence factors. SghA is a newly identified virulence factor enzymatically releasing salicylic acid from its glucoside conjugate and controlling plant tumor development. Here, we report the structural basis of SghR, a LacI-type transcription factor highly conserved in Rhizobiaceae family, regulating the expression of SghA and involved in tumorigenesis. We identified and characterized the binding site of SghR on the promoter region of sghA and then determined the crystal structures of apo-SghR, SghR complexed with its operator DNA, and ligand sucrose, respectively. These results provide detailed insights into how SghR recognizes its cognate DNA and shed a mechanistic light on how sucrose attenuates the affinity of SghR with DNA to modulate the expression of SghA. Given the important role of SghR in mediating the signaling cross-talk during Agrobacterium infection, our results pave the way for structure-based inducer analog design, which has potential applications for agricultural industry.