Promoter melting triggered by bacterial RNA polymerase occurs in three steps

Resistance mechanism of Escherichia coli strains with different ampicillin resistance levels

Antibiotic resistance is an important problem that threatens medical treatment. Differences in the resistance levels of microorganisms cause great difficulties in understanding the mechanisms of antibiotic resistance. Therefore, the molecular reasons underlying the differences in the level of antibiotic resistance need to be clarified. For this purpose, genomic and transcriptomic analyses were performed on three Escherichia coli strains with varying degrees of adaptive resistance to ampicillin. Whole-genome sequencing of strains with different levels of resistance detected five mutations in strains with 10-fold resistance and two additional mutations in strains with 95-fold resistance. Overall, three of the seven mutations occurred as a single base change, while the other four occurred as insertions or deletions. While it was thought that 10-fold resistance was achieved by the effect of mutations in the ftsI, marAR, and rpoC genes, it was found that 95-fold resistance was achieved by the synergistic effect of five mutations and the ampC mutation. In addition, when the general transcriptomic profiles were examined, it was found that similar transcriptomic responses were elicited in strains with different levels of resistance. This study will improve our view of resistance mechanisms in bacteria with different levels of resistance and provide the basis for our understanding of the molecular mechanism of antibiotic resistance in ampicillin-resistant E. coli strains. KEY POINTS: •The mutation of the ampC promoter may act synergistically with other mutations and lead to higher resistance. •Similar transcriptomic responses to ampicillin are induced in strains with different levels of resistance. •Low antibiotic concentrations are the steps that allow rapid achievement of high antibiotic resistance.

Coarse-Grained Models to Study Protein-DNA Interactions and Liquid-Liquid Phase Separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define nonbonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), which is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond time scale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micrometer-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Brewing COFFEE: A Sequence-Specific Coarse-Grained Energy Function for Simulations of DNA-Protein Complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription and gene expression to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse-grained Force Field for Energy Estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the self-organized polymer model with side-chains for proteins and the three interaction site model for DNA in a modular fashion, without recalibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a data set of high-resolution crystal structures. The only parameter in COFFEE is the strength (λDNAPRO) of the DNA-protein contact potential. For an optimal choice of λDNAPRO, the crystallographic B-factors for DNA-protein complexes with varying sizes and topologies are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts scattering profiles that are in quantitative agreement with small-angle X-ray scattering experiments, as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which do not alter the balance of electrostatic interactions but affect chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

Regulatory fine-tuning of mcr-1 increases bacterial fitness and stabilises antibiotic resistance in agricultural settings

Antibiotic resistance tends to carry fitness costs, making it difficult to understand how resistance can be maintained in the absence of continual antibiotic exposure. Here we investigate this problem in the context of mcr-1, a globally disseminated gene that confers resistance to colistin, an agricultural antibiotic that is used as a last resort for the treatment of multi-drug resistant infections. Here we show that regulatory evolution has fine-tuned the expression of mcr-1, allowing E. coli to reduce the fitness cost of mcr-1 while simultaneously increasing colistin resistance. Conjugative plasmids have transferred low-cost/high-resistance mcr-1 alleles across an incredible diversity of E. coli strains, further stabilising mcr-1 at the species level. Regulatory mutations were associated with increased mcr-1 stability in pig farms following a ban on the use of colistin as a growth promoter that decreased colistin consumption by 90%. Our study shows how regulatory evolution and plasmid transfer can combine to stabilise resistance and limit the impact of reducing antibiotic consumption.

Brewing COFFEE: A sequence-specific coarse-grained energy function for simulations of DNA-protein complexes

DNA-protein interactions are pervasive in a number of biophysical processes ranging from transcription, gene expression, to chromosome folding. To describe the structural and dynamic properties underlying these processes accurately, it is important to create transferable computational models. Toward this end, we introduce Coarse grained force field for energy estimation, COFFEE, a robust framework for simulating DNA-protein complexes. To brew COFFEE, we integrated the energy function in the Self-Organized Polymer model with Side Chains for proteins and the Three Interaction Site model for DNA in a modular fashion, without re-calibrating any of the parameters in the original force-fields. A unique feature of COFFEE is that it describes sequence-specific DNA-protein interactions using a statistical potential (SP) derived from a dataset of high-resolution crystal structures. The only parameter in COFFEE is the strength (λ D N A P R O ) of the DNA-protein contact potential. For an optimal choice of λ D N A P R O , the crystallographic B-factors for DNA-protein complexes, with varying sizes and topologies, are quantitatively reproduced. Without any further readjustments to the force-field parameters, COFFEE predicts the scattering profiles that are in quantitative agreement with SAXS experiments as well as chemical shifts that are consistent with NMR. We also show that COFFEE accurately describes the salt-induced unraveling of nucleosomes. Strikingly, our nucleosome simulations explain the destabilization effect of ARG to LYS mutations, which does not alter the balance of electrostatic interactions, but affects chemical interactions in subtle ways. The range of applications attests to the transferability of COFFEE, and we anticipate that it would be a promising framework for simulating DNA-protein complexes at the molecular length-scale.

A coarse-grained DNA model to study protein-DNA interactions and liquid-liquid phase separation

Recent advances in coarse-grained (CG) computational models for DNA have enabled molecular-level insights into the behavior of DNA in complex multiscale systems. However, most existing CG DNA models are not compatible with CG protein models, limiting their applications for emerging topics such as protein-nucleic acid assemblies. Here, we present a new computationally efficient CG DNA model. We first use experimental data to establish the model's ability to predict various aspects of DNA behavior, including melting thermodynamics and relevant local structural properties such as the major and minor grooves. We then employ an all-atom hydropathy scale to define non-bonded interactions between protein and DNA sites, to make our DNA model compatible with an existing CG protein model (HPS-Urry), that is extensively used to study protein phase separation, and show that our new model reasonably reproduces the experimental binding affinity for a prototypical protein-DNA system. To further demonstrate the capabilities of this new model, we simulate a full nucleosome with and without histone tails, on a microsecond timescale, generating conformational ensembles and provide molecular insights into the role of histone tails in influencing the liquid-liquid phase separation (LLPS) of HP1α proteins. We find that histone tails interact favorably with DNA, influencing the conformational ensemble of the DNA and antagonizing the contacts between HP1α and DNA, thus affecting the ability of DNA to promote LLPS of HP1α. These findings shed light on the complex molecular framework that fine-tunes the phase transition properties of heterochromatin proteins and contributes to heterochromatin regulation and function. Overall, the CG DNA model presented here is suitable to facilitate micron-scale studies with sub-nm resolution in many biological and engineering applications and can be used to investigate protein-DNA complexes, such as nucleosomes, or LLPS of proteins with DNA, enabling a mechanistic understanding of how molecular information may be propagated at the genome level.

Modeling DNA Opening in the Eukaryotic Transcription Initiation Complexes via Coarse-Grained Models

Recently, the molecular mechanisms of transcription initiation have been intensively studied. Especially, the cryo-electron microscopy revealed atomic structure details in key states in the eukaryotic transcription initiation. Yet, the dynamic processes of the promoter DNA opening in the pre-initiation complex remain obscured. In this study, based on the three cryo-electron microscopic yeast structures for the closed, open, and initially transcribing complexes, we performed multiscale molecular dynamics (MD) simulations to model structures and dynamic processes of DNA opening. Combining coarse-grained and all-atom MD simulations, we first obtained the atomic model for the DNA bubble in the open complexes. Then, in the MD simulation from the open to the initially transcribing complexes, we found a previously unidentified intermediate state which is formed by the bottleneck in the fork loop 1 of Pol II: The loop opening triggered the escape from the intermediate, serving as a gatekeeper of the promoter DNA opening. In the initially transcribing complex, the non-template DNA strand passes a groove made of the protrusion, the lobe, and the fork of Rpb2 subunit of Pol II, in which several positively charged and highly conserved residues exhibit key interactions to the non-template DNA strand. The back-mapped all-atom models provided further insights on atomistic interactions such as hydrogen bonding and can be used for future simulations.

Theory and simulations of condensin mediated loop extrusion in DNA

Condensation of hundreds of mega-base-pair-long human chromosomes in a small nuclear volume is a spectacular biological phenomenon. This process is driven by the formation of chromosome loops. The ATP consuming motor, condensin, interacts with chromatin segments to actively extrude loops. Motivated by real-time imaging of loop extrusion (LE), we created an analytically solvable model, predicting the LE velocity and step size distribution as a function of external load. The theory fits the available experimental data quantitatively, and suggests that condensin must undergo a large conformational change, induced by ATP binding, bringing distant parts of the motor to proximity. Simulations using a simple model confirm that the motor transitions between an open and a closed state in order to extrude loops by a scrunching mechanism, similar to that proposed in DNA bubble formation during bacterial transcription. Changes in the orientation of the motor domains are transmitted over ~50 nm, connecting the motor head and the hinge, thus providing an allosteric basis for LE.

Effects of Gold Nanoparticles on the Stepping Trajectories of Kinesin

A substantial increase in the temporal resolution of the stepping of dimeric molecular motors is possible by tracking the position of a large gold nanoparticle (GNP) attached to a labeled site on one of the heads. This technique was employed to measure the stepping trajectories of conventional kinesin (Kin1) using the time-dependent position of the GNP as a proxy. The trajectories revealed that the detached head always passes to the right of the head that is tightly bound to the microtubule (MT) during a step. In interpreting the results of such experiments, it is assumed that the GNP does not significantly alter the diffusive motion of the detached head. We used coarse-grained simulations of a system consisting of the MT-Kin1 complex with and without attached GNP to investigate how the stepping trajectories are affected. The two significant findings are: (1) The GNP does not faithfully track the position of the stepping head, and (2) the rightward bias is typically exaggerated by the GNP. Both these findings depend on the precise residue position to which the GNP is attached. Surprisingly, the stepping trajectories of kinesin are not significantly affected if, in addition to the GNP, a 1 μm diameter cargo is attached to the coiled coil. Our simulations suggest the effects of the large probe have to be considered when inferring the stepping mechanisms using GNP tracking experiments.

Asymmetry in histone rotation in forced unwrapping and force quench rewrapping in a nucleosome

Single molecule pulling experiments have shown that DNA in the nucleosomes unwraps in two stages from the histone protein core (HPC). The first stage, attributed to the rupture of the outer DNA turn, occurs between 3 and 5 pNs, and is reversible. The inner DNA turn ruptures irreversibly at forces between 9 and 15 pNs (or higher) in the second stage. Molecular simulations using the Self-Organized Polymer model capture the experimental findings. The unwrapping of the outer DNA turn is independent of the pulling direction. The rupture of the DNA inner turn depends on the pulling direction and involves overcoming substantial energetic (most likely electrostatic in origin) and kinetic barriers. They arise because the mechanical force has to generate sufficient torque to rotate the HPC by 180°. On the other hand, during the rewrapping process, HPC rotation is stochastic, with force playing no role. The assembly of the outer DNA wrap upon force quench nearly coincides with the unwrapping process, confirming the reversibility of the outer turn rupture. The asymmetry in HPC rotation during unwrapping and rewrapping explains the observed hysteresis in the stretch-release cycles in experiments. We propose experiments to test the prediction that HPC rotation produces kinetic barriers in the unwrapping process.

Learning the Regulatory Code of Gene Expression

Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.

Thermodynamic uncertainty relation to assess biological processes

We review the trade-offs between speed, fluctuations, and thermodynamic cost involved with biological processes in nonequilibrium states and discuss how optimal these processes are in light of the universal bound set by the thermodynamic uncertainty relation (TUR). The values of the uncertainty product Q of TUR, which can be used as a measure of the precision of enzymatic processes realized for a given thermodynamic cost, are suboptimal when the substrate concentration is at the Michaelis constant, and some of the key biological processes are found to work around this condition. We illustrate the utility of Q in assessing how close the molecular motors and biomass producing machineries are to the TUR bound, and for the cases of biomass production (or biological copying processes), we discuss how their optimality quantified in terms of Q is balanced with the error rate in the information transfer process. We also touch upon the trade-offs in other error-minimizing processes in biology, such as gene regulation and chaperone-assisted protein folding. A spectrum of Q recapitulating the biological processes surveyed here provides glimpses into how biological systems are evolved to optimize and balance the conflicting functional requirements.

Vitamin B and Vitamin C Affect DNA Methylation and Amino Acid Metabolism in Mycobacterium bovis BCG

Vitamins are essential nutrients and key cofactors of enzymes that regulate cellular metabolism, and also activate the immune system. Recent studies have shown that vitamin B1 (V_{B 1}) and vitamin C (Vc) can inhibit Mycobacterium tuberculosis growth, but the precise mechanism is still not well understood. In the present study, we have used RNA-sequencing (RNA-seq), liquid chromatography coupled to mass spectrometry (LC-MS) and single-molecule real-time (SMRT) sequencing to analyze the transcriptional, metabolic and methylation profiles of Mycobacterium bovis BCG when treated with V_{B 1} and Vc. Our results show that, after vitamin treatment, variant metabolites were mainly clustered in pathways related to amino acid metabolism. Treatment with both vitamins significantly up-regulated the gene encoding cysteine synthase A. Additionally, only BCG that was treated with V_C showed m4c modifications. Genes harboring this methylation were up-regulated, suggesting that m4c methylation can promote gene transcription to some extent. Overall, this study contributes to the understanding of the effects of V_{B 1} and V_C, and suggests that these vitamins constitute potential anti-tuberculosis drugs.

Stepwise Promoter Melting by Bacterial RNA Polymerase

Transcription initiation requires formation of the open promoter complex (RPo). To generate RPo, RNA polymerase (RNAP) unwinds the DNA duplex to form the transcription bubble and loads the DNA into the RNAP active site. RPo formation is a multi-step process with transient intermediates of unknown structure. We use single-particle cryoelectron microscopy to visualize seven intermediates containing Escherichia coli RNAP with the transcription factor TraR en route to forming RPo. The structures span the RPo formation pathway from initial recognition of the duplex promoter in a closed complex to the final RPo. The structures and supporting biochemical data define RNAP and promoter DNA conformational changes that delineate steps on the pathway, including previously undetected transient promoter-RNAP interactions that contribute to populating the intermediates but do not occur in RPo. Our work provides a structural basis for understanding RPo formation and its regulation, a major checkpoint in gene expression throughout evolution.

Gene Silencing Mechanisms Revealed by Dynamics of Guide, Target, and Duplex Binding to Argonaute

Argonaute (Ago) protein plays a central role in silencing gene expression by binding a "guide" strand to the base-pair with a complementary mRNA and degrading the mRNA. The current understanding of how Ago-guide and Ago-guide-mRNA complexes assemble is based mainly on static crystal structures; the associated kinetic pathways remain unknown/unclear. By simulating the successive binding of guide/target strand to Thermus thermophilus Ago (TtAgo) and computing the respective free energy landscapes, we directly visualize how TtAgo silencing complexes form and function. We show that the guide binding rate depends on its initial loading position onto TtAgo. Subsequent target recognition beyond the scissile 10-11 nucleotides must overcome a substantial energy barrier for TtAgo's nucleotide-binding groove to expand widely. This work reveals novel roles for the core TtAgo domains and shows how the kinetic barriers that must be overcome for critical structural changes to occur lead to target repression/cleavage.

Recent Advances in Understanding σ70-Dependent Transcription Initiation Mechanisms

Prokaryotic transcription is one of the most studied biological systems, with relevance to many fields including the development and use of antibiotics, the construction of synthetic gene networks, and the development of many cutting-edge methodologies. Here, we discuss recent structural, biochemical, and single-molecule biophysical studies targeting the mechanisms of transcription initiation in bacteria, including the formation of the open complex, the reaction of initial transcription, and the promoter escape step that leads to elongation. We specifically focus on the mechanisms employed by the RNA polymerase holoenzyme with the housekeeping sigma factor σ⁷⁰. The recent progress provides answers to long-held questions, identifies intriguing new behaviours, and opens up fresh questions for the field of transcription.

Forced-rupture of cell-adhesion complexes reveals abrupt switch between two brittle states

Cell adhesion complexes (CACs), which are activated by ligand binding, play key roles in many cellular functions ranging from cell cycle regulation to mediation of cell extracellular matrix adhesion. Inspired by single molecule pulling experiments using atomic force spectroscopy on leukocyte function-associated antigen-1 (LFA-1), expressed in T-cells, bound to intercellular adhesion molecules (ICAM), we performed constant loading rate (r_f) and constant force (F) simulations using the self-organized polymer model to describe the mechanism of ligand rupture from CACs. The simulations reproduce the major experimental finding on the kinetics of the rupture process, namely, the dependence of the most probable rupture forces (f*s) on ln r_f (r_f is the loading rate) exhibits two distinct linear regimes. The first, at low r_f, has a shallow slope, whereas the slope at high r_f is much larger, especially for a LFA-1/ICAM-1 complex with the transition between the two occurring over a narrow r_f range. Locations of the two transition states (TSs) extracted from the simulations show an abrupt change from a high value at low r_f or constant force, F, to a low value at high r_f or F. This unusual behavior in which the CACs switch from one brittle (TS position is a constant over a range of forces) state to another brittle state is not found in forced-rupture in other protein complexes. We explain this novel behavior by constructing the free energy profiles, F(Λ)s, as a function of a collective reaction coordinate (Λ), involving many key charged residues and a critical metal ion (Mg²⁺). The TS positions in F(Λ), which quantitatively agree with the parameters extracted using the Bell-Evans model, change abruptly at a critical force, demonstrating that it, rather than the molecular extension, is a good reaction coordinate. Our combined analyses using simulations performed in both the pulling modes (constant r_f and F) reveal a new mechanism for the two loading regimes observed in the rupture kinetics in CACs.

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

Determining the Transcription Rates Yielding Steady-State Production of mRNA in the Lac Genetic Switch of Escherichia coli

To elucidate the regulatory dynamics of the gene expression activation and inactivation, an in silico biochemical model of the lac circuit in Escherichia coli was used to evaluate the transcription rates that yield the steady-state mRNA production in active and inactive states of the lac circuit. This result can be used in synthetic biology applications to understand the limits of the genetic synthesis. Since most genetic networks involve many interconnected components with positive and negative feedback control, intuitive understanding of their dynamics is often difficult to obtain. Although the kinetic model of the lac circuit considered involves only a single positive feedback, the developed computational framework can be used to evaluate supported ranges of other reaction rates in genetic circuits with more complex regulatory networks. More specifically, the inducible lac gene switch in E. coli is regulated by unbinding and binding of the inducer-repressor complexes to or from the DNA operator to switch the gene expression on and off. The dependency of mRNA production at steady state on different transcription rates and the repressor complexes has been studied by computer simulations in the Lattice Microbe software. Provided that the lac circuit is in active state, the transcription rate is independent of the inducer-repressor complexes present in the cell. In inactive state, the transcription rate is dependent on the specific inducer-repressor complex bound to the operator that inactivates the gene expression. We found that the repressor complex with the largest affinity to the operator yields the smallest range of the feasible transcription rates to yield the steady state while the lac circuit is in inactive state. In contrast, the steady state in active state can be obtained for any value of the transcription rate.

Studying transcription initiation by RNA polymerase with diffusion-based single-molecule fluorescence

Over the past decade, fluorescence-based single-molecule studies significantly contributed to characterizing the mechanism of RNA polymerase at different steps in transcription, especially in transcription initiation. Transcription by bacterial DNA-dependent RNA polymerase is a multistep process that uses genomic DNA to synthesize complementary RNA molecules. Transcription initiation is a highly regulated step in E. coli, but it has been challenging to study its mechanism because of its stochasticity and complexity. In this review, we describe how single-molecule approaches have contributed to our understanding of transcription and have uncovered mechanistic details that were not observed in conventional assays because of ensemble averaging.

A model for genesis of transcription systems

Repeating sequences generated from RNA gene fusions/ligations dominate ancient life, indicating central importance of building structural complexity in evolving biological systems. A simple and coherent story of life on earth is told from tracking repeating motifs that generate α/β proteins, 2-double-Ψ-β-barrel (DPBB) type RNA polymerases (RNAPs), general transcription factors (GTFs), and promoters. A general rule that emerges is that biological complexity that arises through generation of repeats is often bounded by solubility and closure (i.e., to form a pseudo-dimer or a barrel). Because the first DNA genomes were replicated by DNA template-dependent RNA synthesis followed by RNA template-dependent DNA synthesis via reverse transcriptase, the first DNA replication origins were initially 2-DPBB type RNAP promoters. A simplifying model for evolution of promoters/replication origins via repetition of core promoter elements is proposed. The model can explain why Pribnow boxes in bacterial transcription (i.e., (-12)TATAATG(-6)) so closely resemble TATA boxes (i.e., (-31)TATAAAAG(-24)) in archaeal/eukaryotic transcription. The evolution of anchor DNA sequences in bacterial (i.e., (-35)TTGACA(-30)) and archaeal (BRE(up); BRE for TFB recognition element) promoters is potentially explained. The evolution of BRE(down) elements of archaeal promoters is potentially explained.

Initial events in bacterial transcription initiation

Transcription initiation is a highly regulated step of gene expression. Here, we discuss the series of large conformational changes set in motion by initial specific binding of bacterial RNA polymerase (RNAP) to promoter DNA and their relevance for regulation. Bending and wrapping of the upstream duplex facilitates bending of the downstream duplex into the active site cleft, nucleating opening of 13 bp in the cleft. The rate-determining opening step, driven by binding free energy, forms an unstable open complex, probably with the template strand in the active site. At some promoters, this initial open complex is greatly stabilized by rearrangements of the discriminator region between the -10 element and +1 base of the nontemplate strand and of mobile in-cleft and downstream elements of RNAP. The rate of open complex formation is regulated by effects on the rapidly-reversible steps preceding DNA opening, while open complex lifetime is regulated by effects on the stabilization of the initial open complex. Intrinsic DNA opening-closing appears less regulated. This noncovalent mechanism and its regulation exhibit many analogies to mechanisms of enzyme catalysis.

Base flipping in open complex formation at bacterial promoters

In the process of transcription initiation, the bacterial RNA polymerase binds double-stranded (ds) promoter DNA and subsequently effects strand separation of 12 to 14 base pairs (bp), including the start site of transcription, to form the so-called "open complex" (also referred to as RP(o)). This complex is competent to initiate RNA synthesis. Here we will review the role of σ70 and its homologs in the strand separation process, and evidence that strand separation is initiated at the -11A (the A of the non-template strand that is 11 bp upstream from the transcription start site) of the promoter. By using the fluorescent adenine analog, 2-aminopurine, it was demonstrated that the -11A on the non-template strand flips out of the DNA helix and into a hydrophobic pocket where it stacks with tyrosine 430 of σ70. Open complexes are remarkably stable, even though in vivo, and under most experimental conditions in vitro, dsDNA is much more stable than its strand-separated form. Subsequent structural studies of other researchers have confirmed that in the open complex the -11A has flipped into a hydrophobic pocket of σ70. It was also revealed that RPo was stabilized by three additional bases of the non-template strand being flipped out of the helix and into hydrophobic pockets, further preventing re-annealing of the two complementary DNA strands.

Antibiotics trapping transcription initiation intermediates: To melt or to bend, what's first?

Promoter DNA melting, culminating in the loading of the single-stranded DNA template into the RNA polymerase active site, is a key step in transcription initiation. Recently, the first transcription inhibitors found to block distinct steps of promoter melting were characterized. Here, the impact of these studies is discussed with respect to the current models of transcription initiation.

Bacterial sigma factors: a historical, structural, and genomic perspective

Transcription initiation is the crucial focal point of gene expression in prokaryotes. The key players in this process, sigma factors (σs), associate with the catalytic core RNA polymerase to guide it through the essential steps of initiation: promoter recognition and opening, and synthesis of the first few nucleotides of the transcript. Here we recount the key advances in σ biology, from their discovery 45 years ago to the most recent progress in understanding their structure and function at the atomic level. Recent data provide important structural insights into the mechanisms whereby σs initiate promoter opening. We discuss both the housekeeping σs, which govern transcription of the majority of cellular genes, and the alternative σs, which direct RNA polymerase to specialized operons in response to environmental and physiological cues. The review concludes with a genome-scale view of the extracytoplasmic function σs, the most abundant group of alternative σs.

Perspective: Reaches of chemical physics in biology

Chemical physics as a discipline contributes many experimental tools, algorithms, and fundamental theoretical models that can be applied to biological problems. This is especially true now as the molecular level and the systems level descriptions begin to connect, and multi-scale approaches are being developed to solve cutting edge problems in biology. In some cases, the concepts and tools got their start in non-biological fields, and migrated over, such as the idea of glassy landscapes, fluorescence spectroscopy, or master equation approaches. In other cases, the tools were specifically developed with biological physics applications in mind, such as modeling of single molecule trajectories or super-resolution laser techniques. In this introduction to the special topic section on chemical physics of biological systems, we consider a wide range of contributions, all the way from the molecular level, to molecular assemblies, chemical physics of the cell, and finally systems-level approaches, based on the contributions to this special issue. Chemical physicists can look forward to an exciting future where computational tools, analytical models, and new instrumentation will push the boundaries of biological inquiry.

Structure based design of protein linkers for zinc finger nuclease

Zinc finger nucleases are a promising tool to edit DNA in many biological applications, in particular for gene knockout. Despite many efforts the number of genes that can be effectively targeted with ZFNs remains severely limited, as available constructs cannot address arbitrary gene sequences. Here, we develop a novel concept to significantly enhance the number of DNA sequences that can be targeted by ZFN. Using an efficient computational model, we provide an extensive library of possible linker molecules between individual zinc finger motifs in the construct that can skip up to 10 base pairs between adjacent zinc finger recognition sites in the DNA sequence, which increases the number of genes that can be efficiently targeted by more than an order of magnitude.

Cooperativity and interaction energy threshold effects in recognition of the -10 promoter element by bacterial RNA polymerase

RNA polymerase (RNAP) melts promoter DNA to form transcription-competent open promoter complex (RPo). Interaction of the RNAP σ subunit with non-template strand bases of a conserved -10 element (consensus sequence T-12A-11T-10A-9A-8T-7) is an important source of energy-driving localized promoter melting. Here, we used an RNAP molecular beacon assay to investigate interdependencies of RNAP interactions with -10 element nucleotides. The results reveal a strong cooperation between RNAP interactions with individual -10 element non-template strand nucleotides and indicate that recognition of the -10 element bases occurs only when free energy of the overall RNAP -10 element binding reaches a certain threshold level. The threshold-like mode of the -10 element recognition may be related to the energetic cost of attaining a conformation of the -10 element that is recognizable by RNAP. The RNAP interaction with T/A-12 base pair was found to be strongly stimulated by RNAP interactions with other -10 element bases and with promoter spacer between the -10 and -35 promoter elements. The data also indicate that unmelted -10 promoter element can impair RNAP interactions with promoter DNA upstream of the -11 position. We suggest that cooperativity and threshold effects are important factors guiding the dynamics and selectivity of RPo formation.

Tunable thermal switching via DNA-based nano-devices

DNA has a well-defined structural transition--the denaturation of its double-stranded form into two single strands--that strongly affects its thermal transport properties. We show that, according to a widely implemented model for DNA denaturation, one can engineer DNA 'heattronic' devices that have a rapidly increasing thermal conductance over a narrow temperature range across the denaturation transition (∼350 K). The origin of this rapid increase of conductance, or 'switching', is the softening of the lattice and suppression of nonlinear effects as the temperature crosses the transition temperature and DNA denatures. Most importantly, we demonstrate that DNA nano-junctions have a broad range of thermal tunability by varying the sequence and length, and exploiting the underlying nonlinear behavior. We discuss the role of disorder in the base sequence, as well as the relation to genomic DNA. These results set the basis for developing thermal devices out of materials with nonlinear structural dynamics, as well as understanding the underlying mechanisms of DNA denaturation.

Clostridium difficile TcdC protein binds four-stranded G-quadruplex structures

Clostridium difficile infections are increasing worldwide due to emergence of virulent strains. Infections can result in diarrhea and potentially fatal pseudomembranous colitis. The main virulence factors of C. difficile are clostridial toxins TcdA and TcdB. Transcription of the toxins is positively regulated by the sigma factor TcdR. Negative regulation is believed to occur through TcdC, a proposed anti-sigma factor. Here, we describe the biochemical properties of TcdC to understand the mechanism of TcdC action. Bioinformatic analysis of the TcdC protein sequence predicted the presence of a hydrophobic stretch [amino acids (aa) 30–50], a potential dimerization domain (aa 90–130) and a C-terminal oligonucleotide-binding fold. Gel filtration chromatography of two truncated recombinant TcdC proteins (TcdCΔ1-89 and TcdCΔ1-130) showed that the domain between aa 90 and 130 is involved in dimerization. Binding of recombinant TcdC to single-stranded DNA was studied using a single-stranded Systematic Evolution of Ligands by Exponential enrichment approach. This involved specific binding of single-stranded DNA sequences from a pool of random oligonucleotides, as monitored by electrophoretic-mobility shift assay. Analysis of the oligonucleotides bound showed that the oligonucleotide-binding fold domain of TcdC can bind specifically to DNA folded into G-quadruplex structures containing repetitive guanine nucleotides forming a four-stranded structure. In summary, we provide evidence for DNA binding of TcdC, which suggests an alternative function for this proposed anti-sigma factor.

Single-molecule nanometry for biological physics

Precision measurement is a hallmark of physics but the small length scale (∼nanometer) of elementary biological components and thermal fluctuations surrounding them challenge our ability to visualize their action. Here, we highlight the recent developments in single-molecule nanometry where the position of a single fluorescent molecule can be determined with nanometer precision, reaching the limit imposed by the shot noise, and the relative motion between two molecules can be determined with ∼0.3 nm precision at ∼1 ms time resolution, as well as how these new tools are providing fundamental insights into how motor proteins move on cellular highways. We will also discuss how interactions between three and four fluorescent molecules can be used to measure three and six coordinates, respectively, allowing us to correlate the movements of multiple components. Finally, we will discuss recent progress in combining angstrom-precision optical tweezers with single-molecule fluorescent detection, opening new windows for multi-dimensional single-molecule nanometry for biological physics.

E. coli 6S RNA release from RNA polymerase requires σ70 ejection by scrunching and is orchestrated by a conserved RNA hairpin

The 6S RNA in Escherichia coli suppresses housekeeping transcription by binding to RNA polymerase holoenzyme (core polymerase + σ⁷⁰) under low nutrient conditions and rescues σ⁷⁰-dependent transcription in high nutrient conditions by the synthesis of a short product RNA (pRNA) using itself as a template. Here we characterize a kinetic intermediate that arises during 6S RNA release. This state, consisting of 6S RNA and core polymerase, is related to the formation of a top-strand "release" hairpin that is conserved across the γ-proteobacteria. Deliberately slowing the intrinsic 6S RNA release rate by nucleotide feeding experiments reveals that σ⁷⁰ ejection occurs abruptly once a pRNA length of 9 nucleotides (nt) is reached. After σ⁷⁰ ejection, an additional 4 nt of pRNA synthesis is required before the 6S:pRNA complex is finally released from core polymerase. Changing the E. coli 6S RNA sequence to preclude formation of the release hairpin dramatically slows the speed of 6S RNA release but, surprisingly, does not alter the abruptness of σ⁷⁰ ejection. Rather, the pRNA size required to trigger σ⁷⁰ release increases from 9 nt to 14 nt. That a precise pRNA length is required to trigger σ⁷⁰ release either with or without a hairpin implicates an intrinsic "scrunching"-type release mechanism. We speculate that the release hairpin serves two primary functions in the γ-proteobacteria: First, its formation strips single-stranded "-10" 6S RNA interactions away from σ⁷⁰. Second, the formation of the hairpin accumulates RNA into a region of the polymerase complex previously associated with DNA scrunching, further destabilizing the 6S:pRNA:polymerase complex.

The origin of minus-end directionality and mechanochemistry of Ncd motors

Adaptation of molecular structure to the ligand chemistry and interaction with the cytoskeletal filament are key to understanding the mechanochemistry of molecular motors. Despite the striking structural similarity with kinesin-1, which moves towards plus-end, Ncd motors exhibit minus-end directionality on microtubules (MTs). Here, by employing a structure-based model of protein folding, we show that a simple repositioning of the neck-helix makes the dynamics of Ncd non-processive and minus-end directed as opposed to kinesin-1. Our computational model shows that Ncd in solution can have both symmetric and asymmetric conformations with disparate ADP binding affinity, also revealing that there is a strong correlation between distortion of motor head and decrease in ADP binding affinity in the asymmetric state. The nucleotide (NT) free-ADP (φ-ADP) state bound to MTs favors the symmetric conformation whose coiled-coil stalk points to the plus-end. Upon ATP binding, an enhanced flexibility near the head-neck junction region, which we have identified as the important structural element for directional motility, leads to reorienting the coiled-coil stalk towards the minus-end by stabilizing the asymmetric conformation. The minus-end directionality of the Ncd motor is a remarkable example that demonstrates how motor proteins in the kinesin superfamily diversify their functions by simply rearranging the structural elements peripheral to the catalytic motor head domain.

Proteins belonging to the kinesin superfamily are responsible for vesicle or organelle transport, spindle morphogenesis, and chromosome sorting during cell division. Interestingly, while most proteins in kinesin superfamily that share the common catalytic motor head domain have plus-end directionality along microtubules, kinesin-14 (Ncd) exhibits minus-end directionality. Despite the several circumstantial evidences on the determining factors for the motor directionality in the last decade, a comprehensive understanding of the mechanism governing the Ncd minus-end directionality is still missing. Our studies provide a clear explanation for this minus-end directionality and the associated mechanochemical cycle. Here, we modeled an Ncd motor by employing structural details available in the literature to simulate its conformational dynamics. Simulations using our structure-based model of Ncd assert that the dynamics due to a simple rearrangement of structural elements, peripheral to the catalytic motor domain, is the key player in determining both the directionality and mechanochemistry unique to Ncd motors. Although Ncd has a different directionality, it uses a similar strategy to kinesin-1 of structural adaptation of the catalytic motor domain. Therefore using the same physical principle of protein folding and very similar structural elements, motors in the kinesin superfamily are able to achieve a variety of biological function.

A common feature from different subunits of a homomeric AAA+ protein contacts three spatially distinct transcription elements

Initiation of σ(54)-dependent transcription requires assistance to melt DNA at the promoter site but is impeded by numerous protein-protein and nucleo-protein interactions. To alleviate these inhibitory interactions, hexameric bacterial enhancer binding proteins (bEBP), a subset of the ATPases associated with various cellular activities (AAA+) protein family, are required to remodel the transcription complex using energy derived from ATP hydrolysis. However, neither the process of energy conversion nor the internal architecture of the closed promoter complex is well understood. Escherichia coli Phage shock protein F (PspF), a well-studied bEBP, contains a surface-exposed loop 1 (L1). L1 is key to the energy coupling process by interacting with Region I of σ(54) (σ(54)(RI)) in a nucleotide dependent manner. Our analyses uncover new levels of complexity in the engagement of a multimeric bEBP with a basal transcription complex via several L1s. The mechanistic implications for these multivalent L1 interactions are elaborated in the light of available structures for the bEBP and its target complexes.

Influence of DNA template choice on transcription and inhibition of Escherichia coli RNA polymerase

In recent decades, quantitative transcription assays using bacterial RNA polymerase (RNAP) have been performed under widely diverse experimental conditions. We demonstrate that the template choice can influence the inhibitory potency of RNAP inhibitors. Furthermore, we illustrate that the sigma factor (σ(70)) surprisingly increases the transcription efficiency of templates with nonphysiological nonprokaryotic promoters. Our results might be a useful guideline in the early stages of using RNAP for drug discovery.

A structural perspective on the dynamics of kinesin motors

Despite significant fluctuation under thermal noise, biological machines in cells perform their tasks with exquisite precision. Using molecular simulation of a coarse-grained model and theoretical arguments, we envisaged how kinesin, a prototype of biological machines, generates force and regulates its dynamics to sustain persistent motor action. A structure-based model, which can be versatile in adapting its structure to external stresses while maintaining its native fold, was employed to account for several features of kinesin dynamics along the biochemical cycle. This analysis complements our current understandings of kinesin dynamics and connections to experiments. We propose a thermodynamic cycle for kinesin that emphasizes the mechanical and regulatory role of the neck linker and clarify issues related to the motor directionality, and the difference between the external stalling force and the internal tension responsible for the head-head coordination. The comparison between the thermodynamic cycle of kinesin and macroscopic heat engines highlights the importance of structural change as the source of work production in biomolecular machines.

Beyond DNA origami: the unfolding prospects of nucleic acid nanotechnology

Nucleic acid nanotechnology exploits the programmable molecular recognition properties of natural and synthetic nucleic acids to assemble structures with nanometer-scale precision. In 2006, DNA origami transformed the field by providing a versatile platform for self-assembly of arbitrary shapes from one long DNA strand held in place by hundreds of short, site-specific (spatially addressable) DNA 'staples'. This revolutionary approach has led to the creation of a multitude of two-dimensional and three-dimensional scaffolds that form the basis for functional nanodevices. Not limited to nucleic acids, these nanodevices can incorporate other structural and functional materials, such as proteins and nanoparticles, making them broadly useful for current and future applications in emerging fields such as nanomedicine, nanoelectronics, and alternative energy.

Coarse-grained molecular simulations of large biomolecules

Recently, we saw a dramatic increase in the number of researches that rely on coarse-grained (CG) simulations for large biomolecules. Here, first, we briefly describe recently developed and used CG models for proteins and nucleic acids. Balance between structure-based and physico-chemical terms is a key issue. We also discuss the multiscale algorithms used to derive CG parameters. Next, we comment on the dynamics used in CG simulations with an emphasis on the importance of hydrodynamic interactions. We then discuss the pros and cons of CG simulations. Finally, we overview recent exciting applications of CG simulations. Publicly available tools and software for CG simulations are also summarized.

Structural basis for promoter-10 element recognition by the bacterial RNA polymerase σ subunit

The key step in bacterial promoter opening is recognition of the -10 promoter element (T(-12)A(-11)T(-10)A(-9)A(-8)T(-7) consensus sequence) by the RNA polymerase σ subunit. We determined crystal structures of σ domain 2 bound to single-stranded DNA bearing-10 element sequences. Extensive interactions occur between the protein and the DNA backbone of every -10 element nucleotide. Base-specific interactions occur primarily with A(-11) and T(-7), which are flipped out of the single-stranded DNA base stack and buried deep in protein pockets. The structures, along with biochemical data, support a model where the recognition of the -10 element sequence drives initial promoter opening as the bases of the nontemplate strand are extruded from the DNA double-helix and captured by σ. These results provide a detailed structural basis for the critical roles of A(-11) and T(-7) in promoter melting and reveal important insights into the initiation of transcription bubble formation.

Link between allosteric signal transduction and functional dynamics in a multisubunit enzyme: S-adenosylhomocysteine hydrolase

S-adenosylhomocysteine hydrolase (SAHH), a cellular enzyme that plays a key role in methylation reactions including those required for maturation of viral mRNA, is an important drug target in the discovery of antiviral agents. While targeting the active site is a straightforward strategy of enzyme inhibition, evidence of allosteric modulation of active site in many enzymes underscores the molecular origin of signal transduction. Information of co-evolving sequences in SAHH family and the key residues for functional dynamics that can be identified using native topology of the enzyme provide glimpses into how the allosteric signaling network, dispersed over the molecular structure, coordinates intra- and intersubunit conformational dynamics. To study the link between the allosteric communication and functional dynamics of SAHHs, we performed Brownian dynamics simulations by building a coarse-grained model based on the holo and ligand-bound structures. The simulations of ligand-induced transition revealed that the signal of intrasubunit closure dynamics is transmitted to form intersubunit contacts, which in turn invoke a precise alignment of active site, followed by the dimer-dimer rotation that compacts the whole tetrameric structure. Further analyses of SAHH dynamics associated with ligand binding provided evidence of both induced fit and population shift mechanisms and also showed that the transition-state ensemble is akin to the ligand-bound state. Besides the formation of enzyme-ligand contacts at the active site, the allosteric couplings from the residues distal to the active site are vital to the enzymatic function.

Convergent transcription in the butyrolactone regulon in Streptomyces coelicolor confers a bistable genetic switch for antibiotic biosynthesis

cis-encoded antisense RNAs (cis asRNA) have been reported to participate in gene expression regulation in both eukaryotic and prokaryotic organisms. Its presence in Streptomyces coelicolor has also been reported recently; however, its role has yet to be fully investigated. Using mathematical modeling we explore the role of cis asRNA produced as a result of convergent transcription in scbA-scbR genetic switch. scbA and scbR gene pair, encoding repressor–amplifier proteins respectively, mediates the synthesis of a signaling molecule, the γ-butyrolactone SCB1 and controls the onset of antibiotic production. Our model considers that transcriptional interference caused by convergent transcription of two opposing RNA polymerases results in fatal collision and transcriptional termination, which suppresses transcription efficiency. Additionally, convergent transcription causes sense and antisense interactions between complementary sequences from opposing strands, rendering the full length transcript inaccessible for translation. We evaluated the role of transcriptional interference and the antisense effect conferred by convergent transcription on the behavior of scbA-scbR system. Stability analysis showed that while transcriptional interference affects the system, it is asRNA that confers scbA-scbR system the characteristics of a bistable switch in response to the signaling molecule SCB1. With its critical role of regulating the onset of antibiotic synthesis the bistable behavior offers this two gene system the needed robustness to be a genetic switch. The convergent two gene system with potential of transcriptional interference is a frequent feature in various genomes. The possibility of asRNA regulation in other such gene-pairs is yet to be examined.

Promoter spacer DNA plays an active role in integrating the functional consequences of RNA polymerase contacts with -10 and -35 promoter elements

Bacterial RNA polymerase (RNAP) interacts with conserved -10 and -35 promoter elements to recognize the promoter and to form an open complex in which DNA duplex around transcription start site melts. Using model DNA constructs (fork junction DNA) that mimic DNA structure found in the open complex we observed that the consequences of mutations in -10 promoter element for RNAP binding exhibited a striking dependence on the presence or absence of a functional -35 promoter element. A role of spacer DNA (a non-conserved DNA sequence connecting -10 and -35 promoter elements) in this phenomenon was probed with a series of fork junction DNA constructs containing perturbations to the spacer DNA. In the absence of a physical connection between the -10 and -35 DNA elements, or when -10 and -35 DNA elements were connected by a long flexible non-DNA linker, the dependence of RNAP interactions with -10 element on the strength of -35 element was lost. When these DNA elements were linked by a rigid DNA duplex or by a DNA duplex containing a short single-stranded gap, the coupling between the -10 and -35 binding activities was observed. These results indicated that promoter spacer DNA played an active role in integrating the functional consequences of RNA polymerase contacts with -10 and -35 promoter element. This role likely involves physical deformation of the spacer occurring in parallel with promoter melting as shown by Fluorescence Resonance Energy Transfer (FRET) experiments with the probes incorporated into spacer DNA.

Driving denaturation: nanoscale thermal transport as a probe of DNA melting

DNA denaturation has been a subject of intense study due to its relationship to DNA transcription and its fundamental importance as a nonlinear structural transition. Many aspects of this phenomenon, however, remain poorly understood. Existing models fit quite well with experimental results on the fraction of unbound base pairs versus temperature, but yield incorrect results for other essential quantities such as the base pair fluctuation time scales. Here we demonstrate that nanoscale thermal transport can serve as a sensitive probe of the underlying microscopic physics responsible for the dynamics of DNA denaturation. Specifically, we show that the heat transport properties of DNA are altered significantly as it denatures, and this alteration encodes detailed information on the dynamics of thermal fluctuations and their interaction along the strand. This finding allows for the discrimination between models of DNA denaturation and will help shed new light on the nonlinear vibrational dynamics of this important molecule.