Structural basis for cooperative DNA binding by CAP and lac repressor

Minicircle Biopharmaceuticals–An Overview of Purification Strategies

Minicircles are non-viral delivery vectors with promising features for biopharmaceutical applications. These vectors are plasmid-derived circular DNA molecules that are obtained in vivo in Escherichia coli by the intramolecular recombination of a parental plasmid, which generates a minicircle containing the eukaryotic therapeutic cassette of interest and a miniplasmid containing the prokaryotic backbone. The production process results thus in a complex mixture, which hinders the isolation of minicircle molecules from other DNA molecules. Several strategies have been proposed over the years to meet the challenge of purifying and obtaining high quality minicircles in compliance with the regulatory guidelines for therapeutic use. In minicircle purification, the characteristics of the strain and parental plasmid used have a high impact and strongly affect the purification strategy that can be applied. This review summarizes the different methods developed so far, focusing not only on the purification method itself but also on its dependence on the upstream production strategy used.

Ageing and sources of transcriptional heterogeneity

Cellular heterogeneity is an important contributor to biological function and is employed by cells, tissues and organisms to adapt, compensate, respond, defend and/or regulate specific processes. Research over the last decades has revealed that transcriptional noise is a major driver for cell-to-cell variability. In this review we will discuss sources of transcriptional variability, in particular bursting of gene expression and how it could contribute to cellular states and fate decisions. We will highlight recent developments in single cell sequencing technologies that make it possible to address cellular heterogeneity in unprecedented detail. Finally, we will review recent literature, in which these new technologies are harnessed to address pressing questions in the field of ageing research, such as transcriptional noise and cellular heterogeneity in the course of ageing.

Long Time-Scale Atomistic Simulations of the Structure and Dynamics of Transcription Factor-DNA Recognition

Recent years have witnessed an explosion of interest in computational studies of DNA binding proteins, including both coarse-grained and atomistic simulations of transcription factor-DNA recognition, to understand how these transcription factors recognize their binding sites on the DNA with such exquisite specificity. The present study performs microsecond time scale all-atom simulations of the dimeric form of the lactose repressor (LacI), both in the absence of any DNA and in the presence of both specific and nonspecific complexes, considering three different DNA sequences. We examine, specifically, the conformational differences between specific and nonspecific protein-DNA interactions, as well as the behavior of the helix-turn-helix motif of LacI when interacting with the DNA. Our simulations suggest that stable LacI binding occurs primarily to bent A-form DNA, with a loss of LacI conformational entropy and optimization of correlated conformational equilibria across the protein. In addition, binding to the specific operator sequence involves a slightly larger number of stabilizing DNA-protein hydrogen bonds (in comparison to nonspecific complexes), which may account for the experimentally observed specificity for this operator. In doing so, our simulations provide a detailed atomistic description of potential structural drivers for LacI selectivity.

The mechanism and high-free-energy transition state of lac repressor-lac operator interaction

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)–lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (k_a) and dissociation (k_d) rate constants. RO-destabilizing solutes (urea, KCl) reduce k_a comparably (urea) or more than (KCl) they increase k_d, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce k_d without affecting k_a significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Direct observation of a 91 bp LacI-mediated, negatively supercoiled DNA loop by atomic force microscope

Escherichia coli lactose repressor (LacI), a tetrameric protein, is a paradigmatic transcriptional factor that controls the expression of lacZYA in the lac operon. It specifically binds to the O1, O2, and O3 operators of the lac promoter, forms DNA loops, and regulates transcription of the lac operon. In this article, utilizing combined techniques of DNA-nicking assay and AFM imaging, we directly observed a 91 bp LacI-mediated, negatively supercoiled DNA loop mimicking the DNA loop between the O1 and O3 operators in the lac promoter.

DNA supercoiling, a critical signal regulating the basal expression of the lac operon in Escherichia coli

Escherichia coli lac repressor (LacI) is a paradigmatic transcriptional factor that controls the expression of lacZYA in the lac operon. This tetrameric protein specifically binds to the O1, O2 and O3 operators of the lac operon and forms a DNA loop to repress transcription from the adjacent lac promoter. In this article, we demonstrate that upon binding to the O1 and O2 operators at their native positions LacI constrains three (−) supercoils within the 401-bp DNA loop of the lac promoter and forms a topological barrier. The stability of LacI-mediated DNA topological barriers is directly proportional to its DNA binding affinity. However, we find that DNA supercoiling modulates the basal expression from the lac operon in E. coli. Our results are consistent with the hypothesis that LacI functions as a topological barrier to constrain free, unconstrained (−) supercoils within the 401-bp DNA loop of the lac promoter. These constrained (−) supercoils enhance LacI’s DNA-binding affinity and thereby the repression of the promoter. Thus, LacI binding is superhelically modulated to control the expression of lacZYA in the lac operon under varying growth conditions.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

Conformational Modeling of Continuum Structures in Robotics and Structural Biology: A Review

Hyper-redundant (or snakelike) manipulators have many more degrees of freedom than are required to position and orient an object in space. They have been employed in a variety of applications ranging from search-and-rescue to minimally invasive surgical procedures, and recently they even have been proposed as solutions to problems in maintaining civil infrastructure and the repair of satellites. The kinematic and dynamic properties of snakelike robots are captured naturally using a continuum backbone curve equipped with a naturally evolving set of reference frames, stiffness properties, and mass density. When the snakelike robot has a continuum architecture, the backbone curve corresponds with the physical device itself. Interestingly, these same modeling ideas can be used to describe conformational shapes of DNA molecules and filamentous protein structures in solution and in cells. This paper reviews several classes of snakelike robots: (1) hyper-redundant manipulators guided by backbone curves; (2) flexible steerable needles; and (3) concentric tube continuum robots. It is then shown how the same mathematical modeling methods used in these robotics contexts can be used to model molecules such as DNA. All of these problems are treated in the context of a common mathematical framework based on the differential geometry of curves, continuum mechanics, and variational calculus. Both coordinate-dependent Euler-Lagrange formulations and coordinate-free Euler-Poincaré approaches are reviewed.

Kinetic approaches to lactose operon induction and bimodality

The quasi-equilibrium approximation is acceptable when molecular interactions are fast enough compared to circuit dynamics, but is no longer allowed when cellular activities are governed by rare events. A typical example is the lactose operon (lac), one of the most famous paradigms of transcription regulation, for which several theories still coexist to describe its behaviors. The lac system is generally analyzed by using equilibrium constants, contradicting single-event hypotheses long suggested by Novick and Weiner (1957). Enzyme induction as an all-or-none phenomenon. Proc. Natl. Acad. Sci. USA 43, 553-566) and recently refined in the study of (Choi et al., 2008. A stochastic single-molecule event triggers phenotype switching of a bacterial cell. Science 322, 442-446). In the present report, a lac repressor (LacI)-mediated DNA immunoprecipitation experiment reveals that the natural LacI-lac DNA complex built in vivo is extremely tight and long-lived compared to the time scale of lac expression dynamics, which could functionally disconnect the abortive expression bursts and forbid using the standard modes of lac bistability. As alternatives, purely kinetic mechanisms are examined for their capacity to restrict induction through: (i) widely scattered derepression related to the arrival time variance of a predominantly backward asymmetric random walk and (ii) an induction threshold arising in a single window of derepression without recourse to nonlinear multimeric binding and Hill functions. Considering the complete disengagement of the lac repressor from the lac promoter as the probabilistic consequence of a transient stepwise mechanism, is sufficient to explain the sigmoidal lac responses as functions of time and of inducer concentration. This sigmoidal shape can be misleadingly interpreted as a phenomenon of equilibrium cooperativity classically used to explain bistability, but which has been reported to be weak in this system.

Interplay of protein and DNA structure revealed in simulations of the lac operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.

Gene repression by minimal lac loops in vivo

The inflexibility of double-stranded DNA with respect to bending and twisting is well established in vitro. Understanding apparent DNA physical properties in vivo is a greater challenge. Here, we exploit repression looping with components of the Escherichia coli lac operon to monitor DNA flexibility in living cells. We create a minimal system for testing the shortest possible DNA repression loops that contain an E. coli promoter, and compare the results to prior experiments. Our data reveal that loop-independent repression occurs for certain tight operator/promoter spacings. When only loop-dependent repression is considered, fits to a thermodynamic model show that DNA twisting limits looping in vivo, although the apparent DNA twist flexibility is 2- to 4-fold higher than in vitro. In contrast, length-dependent resistance to DNA bending is not observed in these experiments, even for the shortest loops constraining <0.4 persistence lengths of DNA. As observed previously for other looping configurations, loss of the nucleoid protein heat unstable (HU) markedly disables DNA looping in vivo. Length-independent DNA bending energy may reflect the activities of architectural proteins and the structure of the DNA topological domain. We suggest that the shortest loops are formed in apical loops rather than along the DNA plectonemic superhelix.

A loop-controlled rrnB P1 promoter for high-level expression of heterologous proteins in Escherichia coli

An effective protein expression system was constructed in Escherichia coli utilizing the rRNA rrnB P1 promoter and the regulatory element of the lac operon (lacO). To regulate the transcriptional activity of the rrnB P1 promoter, we designed two lacO sites with an intervening loop structure; expression was verified by measuring the levels of the β-1,4-glucanase gene, cel5G. Basal expression from the looped promoter construct was reduced by 92% when compared to expression from the T7 promoter. We also found that the host cell type had a significant effect on the regulation of the rrnB P1 promoter: E. coli DH5α and DH10B had high expression levels, whereas the expression in BL21(DE3) was more stringent.

Group theory and biomolecular conformation: I. Mathematical and computational models

Biological macromolecules, and the complexes that they form, can be described in a variety of ways ranging from quantum mechanical and atomic chemical models, to coarser grained models of secondary structure and domains, to continuum models. At each of these levels, group theory can be used to describe both geometric symmetries and conformational motion. In this survey, a detailed account is provided of how group theory has been applied across computational structural biology to analyze the conformational shape and motion of macromolecules and complexes.

Wrapped-around models for the lac operon complex

The protein-DNA complex, involved in the lac operon of enteric bacteria, is paradigmatic in understanding the extent of DNA bending and plasticity due to interactions with protein assemblies acting as DNA regulators. For the lac operon, two classes of structures have been proposed: 1), with the protein tetramer lying away from the DNA loop (wrapped-away model); and 2), with the protein tetramer lying inside the DNA loop (wrapped-around model). A recently developed electrostatic analytical model shows that the size and net charge of the Lac protein tetramer allow the bending of DNA, which is consistent with another wrapped-around model from the literature. Coarse-grained models, designed based on this observation, are extensively investigated and show three kinds of wrapped-around arrangements of DNA and a lower propensity for wrapped-away configurations. Molecular dynamics simulations of an all-atom model, built on the basis of the most tightly collapsed coarse-grained model, show that most of the DNA double-helical architecture is maintained in the region between O3 and O1 DNA operators, that the DNA distortion is concentrated in the chain beyond the O1 operator, and that the protein tetramer can adapt the N-terminal domains to the DNA tension.

DNA curvature and flexibility in vitro and in vivo

It has been more than 50 years since the elucidation of the structure of double-helical DNA. Despite active research and progress in DNA biology and biochemistry, much remains to be learned in the field of DNA biophysics. Predicting the sequence-dependent curvature and flexibility of DNA is difficult. Applicability of the conventional worm-like chain polymer model of DNA has been challenged. The fundamental forces responsible for the remarkable resistance of DNA to bending and twisting remain controversial. The apparent 'softening' of DNA measured in vivo in the presence of kinking proteins and superhelical strain is incompletely understood. New methods and insights are being applied to these problems. This review places current work on DNA biophysics in historical context and illustrates the ongoing interplay between theory and experiment in this exciting field.

Tetramer opening in LacI-mediated DNA looping

Lactose repressor protein (LacI) controls transcription of the genes involved in lactose metabolism in bacteria. Essential to optimal LacI-mediated regulation is its ability to bind simultaneously to two operators, forming a loop on the intervening DNA. Recently, several lines of evidence (both theoretical and experimental) have suggested various possible loop structures associated with different DNA binding topologies and LacI tetramer structural conformations (adopted by flexing about the C-terminal tetramerization domain). We address, specifically, the role of protein opening in loop formation by employing the single-molecule tethered particle motion method on LacI protein mutants chemically cross-linked at different positions along the cleft between the two dimers. Measurements on the wild-type and uncross-linked LacI mutants led to the observation of two distinct levels of short tether length, associated with two different DNA looping structures. Restricting conformational flexibility of the protein by chemical cross-linking induces pronounced effects. Crosslinking the dimers at the level of the N-terminal DNA binding head (E36C) completely suppresses looping, whereas cross-linking near the C-terminal tetramerization domain (Q231C) results in changes of looping geometry detected by the measured tether length distributions. These observations lead to the conclusion that tetramer opening plays a definite role in at least a subset of LacI/DNA loop conformations.

Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.

The Stochastic Elastica and Excluded-Volume Perturbations of DNA Conformational Ensembles

A coordinate-free Lie-group formulation for generating ensembles of DNA conformations in solution is presented. In this formulation, stochastic differential equations define sample paths on the Euclidean motion group. The ensemble of these paths exhibits the same behavior as solutions of the Fokker-Planck equation for the stochastically forced elastica. Longer chains for which the effects of excluded volume become important are handled by piecing together shorter chains and modeling their interactions. It is assumed that the final chain lengths of interest are long enough for excluded volume effects to become important, but not so long that the semi-flexible nature of the chain is lost. The effect of excluded volume is then taken into account by grouping short self-avoiding conformations into 'bundles' with common end constraints and computing average interaction effects between bundles. The accuracy of this approximation is shown to be good when using a numerically generated ensemble of self-avoiding sample paths as the baseline for comparison.

Mesoscale modeling of multi-protein-DNA assemblies: the role of the catabolic activator protein in Lac-repressor-mediated looping

DNA looping plays a key role in the regulation of the lac operon in Escherichia coli. The presence of a tightly bent loop (between sequentially distant sites of Lac repressor protein binding) purportedly hinders the binding of RNA polymerase and subsequent transcription of the genetic message. The unexpectedly favorable binding interaction of this protein-DNA assembly with the catabolic activator protein (CAP), a protein that also bends DNA and paradoxically facilitates the binding of RNA polymerase, stimulated extension of our base-pair level theory of DNA elasticity to the treatment of DNA loops formed in the presence of several proteins. Here we describe in detail a procedure to determine the structures and free energies of multi-protein-DNA assemblies and illustrate the predicted effects of CAP binding on the configurations of the wild-type 92-bp Lac repressor-mediated O3-O1 DNA loop. We show that the DNA loop adopts an antiparallel orientation in the most likely structure and that this loop accounts for the published experimental observation that, when CAP is bound to the loop, one of the arms of LacR binds to an alternative site that is displaced from the original site by 5 bp.

Interconvertible lac repressor-DNA loops revealed by single-molecule experiments

At many promoters, transcription is regulated by simultaneous binding of a protein to multiple sites on DNA, but the structures and dynamics of such transcription factor-mediated DNA loops are poorly understood. We directly examined in vitro loop formation mediated by Escherichia coli lactose repressor using single-molecule structural and kinetics methods. Small ( approximately 150 bp) loops form quickly and stably, even with out-of-phase operator spacings. Unexpectedly, repeated spontaneous transitions between two distinct loop structures were observed in individual protein-DNA complexes. The results imply a dynamic equilibrium between a novel loop structure with the repressor in its crystallographic "V" conformation and a second structure with a more extended linear repressor conformation that substantially lessens the DNA bending strain. The ability to switch between different loop structures may help to explain how robust transcription regulation is maintained even though the mechanical work required to form a loop may change substantially with metabolic conditions.

Single-molecule manipulation reveals supercoiling-dependent modulation of lac repressor-mediated DNA looping

Gene expression regulation is a fundamental biological process which deploys specific sets of genomic information depending on physiological or environmental conditions. Several transcription factors (including lac repressor, LacI) are present in the cell at very low copy number and increase their local concentration by binding to multiple sites on DNA and looping the intervening sequence. In this work, we employ single-molecule manipulation to experimentally address the role of DNA supercoiling in the dynamics and stability of LacI-mediated DNA looping. We performed measurements over a range of degrees of supercoiling between -0.026 and +0.026, in the absence of axial stretching forces. A supercoiling-dependent modulation of the lifetimes of both the looped and unlooped states was observed. Our experiments also provide evidence for multiple structural conformations of the LacI-DNA complex, depending on torsional constraints. The supercoiling-dependent modulation demonstrated here adds an important element to the model of the lac operon. In fact, the complex network of proteins acting on the DNA in a living cell constantly modifies its topological and mechanical properties: our observations demonstrate the possibility of establishing a signaling pathway from factors affecting DNA supercoiling to transcription factors responsible for the regulation of specific sets of genes.

Intrinsic curvature of DNA influences LacR-mediated looping

Protein-mediated DNA looping is a common mechanism for regulating gene expression. Loops occur when a protein binds to two operators on the same DNA molecule. The probability of looping is controlled, in part, by the basepair sequence of inter-operator DNA, which influences its structural properties. One structural property is the intrinsic or stress-free curvature. In this article, we explore the influence of sequence-dependent intrinsic curvature by exercising a computational rod model for the inter-operator DNA as applied to looping of the LacR-DNA complex. Starting with known sequences for the inter-operator DNA, we first compute the intrinsic curvature of the helical axis as input to the rod model. The crystal structure of the LacR (with bound operators) then defines the requisite boundary conditions needed for the dynamic rod model that predicts the energetics and topology of the intervening DNA loop. A major contribution of this model is its ability to predict a broad range of published experimental data for highly bent (designed) sequences. The model successfully predicts the loop topologies known from fluorescence resonance energy transfer measurements, the linking number distribution known from cyclization assays with the LacR-DNA complex, the relative loop stability known from competition assays, and the relative loop size known from gel mobility assays. In addition, the computations reveal that highly curved sequences tend to lower the energetic cost of loop formation, widen the energy distribution among stable and meta-stable looped states, and substantially alter loop topology. The inclusion of sequence-dependent intrinsic curvature also leads to nonuniform twist and necessitates consideration of eight distinct binding topologies from the known crystal structure of the LacR-DNA complex.

Multilevel deconstruction of the In vivo behavior of looped DNA-protein complexes

Protein-DNA complexes with loops play a fundamental role in a wide variety of cellular processes, ranging from the regulation of DNA transcription to telomere maintenance. As ubiquitous as they are, their precise in vivo properties and their integration into the cellular function still remain largely unexplored. Here, we present a multilevel approach that efficiently connects in both directions molecular properties with cell physiology and use it to characterize the molecular properties of the looped DNA-lac repressor complex while functioning in vivo. The properties we uncover include the presence of two representative conformations of the complex, the stabilization of one conformation by DNA architectural proteins, and precise values of the underlying twisting elastic constants and bending free energies. Incorporation of all this molecular information into gene-regulation models reveals an unprecedented versatility of looped DNA-protein complexes at shaping the properties of gene expression.

Combinatorial transcriptional control of the lactose operon of Escherichia coli

The goal of systems biology is to understand the behavior of the whole in terms of knowledge of the parts. This is hard to achieve in many cases due to the difficulty of characterizing the many constituents involved in a biological system and their complex web of interactions. The lac promoter of Escherichia coli offers the possibility of confronting "system-level" properties of transcriptional regulation with the known biochemistry of the molecular constituents and their mutual interactions. Such confrontations can reveal previously unknown constituents and interactions, as well as offer insight into how the components work together as a whole. Here we study the combinatorial control of the lac promoter by the regulators Lac repressor (LacR) and cAMP-receptor protein (CRP). A previous in vivo study [Setty Y, Mayo AE, Surette MG, Alon U (2003) Proc Natl Acad Sci USA 100:7702-7707] found gross disagreement between the observed promoter activities and the expected behavior based on the known molecular mechanisms. We repeated the study by identifying and removing several extraneous factors that significantly modulated the expression of the lac promoter. Through quantitative, systematic characterization of promoter activity for a number of key mutants and guided by the thermodynamic model of transcriptional regulation, we were able to account for the combinatorial control of the lac promoter quantitatively, in terms of a cooperative interaction between CRP and LacR-mediated DNA looping. Specifically, our analysis indicates that the sensitivity of the inducer response results from LacR-mediated DNA looping, which is significantly enhanced by CRP.

How phosphotransferase system-related protein phosphorylation regulates carbohydrate metabolism in bacteria

The phosphoenolpyruvate(PEP):carbohydrate phosphotransferase system (PTS) is found only in bacteria, where it catalyzes the transport and phosphorylation of numerous monosaccharides, disaccharides, amino sugars, polyols, and other sugar derivatives. To carry out its catalytic function in sugar transport and phosphorylation, the PTS uses PEP as an energy source and phosphoryl donor. The phosphoryl group of PEP is usually transferred via four distinct proteins (domains) to the transported sugar bound to the respective membrane component(s) (EIIC and EIID) of the PTS. The organization of the PTS as a four-step phosphoryl transfer system, in which all P derivatives exhibit similar energy (phosphorylation occurs at histidyl or cysteyl residues), is surprising, as a single protein (or domain) coupling energy transfer and sugar phosphorylation would be sufficient for PTS function. A possible explanation for the complexity of the PTS was provided by the discovery that the PTS also carries out numerous regulatory functions. Depending on their phosphorylation state, the four proteins (domains) forming the PTS phosphorylation cascade (EI, HPr, EIIA, and EIIB) can phosphorylate or interact with numerous non-PTS proteins and thereby regulate their activity. In addition, in certain bacteria, one of the PTS components (HPr) is phosphorylated by ATP at a seryl residue, which increases the complexity of PTS-mediated regulation. In this review, we try to summarize the known protein phosphorylation-related regulatory functions of the PTS. As we shall see, the PTS regulation network not only controls carbohydrate uptake and metabolism but also interferes with the utilization of nitrogen and phosphorus and the virulence of certain pathogens.

Effect of supercoiling on formation of protein-mediated DNA loops

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short (50-100 basepair), protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay.

Molecular mechanisms of cellular mechanics

Mechanical forces play an essential role in cellular processes as input, output, and signals. Various protein complexes in the cell are designed to handle, transform and use such forces. For instance, proteins of muscle and the extracellular matrix can withstand considerable stretching forces, hearing-related and mechanosensory proteins can transform weak mechanical stimuli into electrical signals, and regulatory proteins are suited to forcing DNA into loops to control gene expression. Here we review the structure-function relationship of four protein complexes with well defined and representative mechanical functions. The first example is titin, a protein that confers passive elasticity on muscle. The second system is the elastic extracellular matrix protein, fibronectin, and its cellular receptor integrin. The third protein system is the transduction apparatus in hearing and other mechanical senses, likely containing cadherin and ankyrin repeats. The last system is the lac repressor protein, which regulates gene expression by looping DNA. This review focuses on atomic level descriptions of the physical mechanisms underlying the various mechanical functions of the stated proteins.

Do femtonewton forces affect genetic function? A review

Protein-Mediated DNA looping is intricately related to gene expression. Therefore any mechanical constraint that disrupts loop formation can play a significant role in gene regulation. Polymer physics models predict that less than a piconewton of force may be sufficient to prevent the formation of DNA loops. Thus, it appears that tension can act as a molecular switch that controls the much larger forces associated with the processive motion of RNA polymerase. Since RNAP can exert forces over 20 pN before it stalls, a 'substrate tension switch' could offer a force advantage of two orders of magnitude. Evidence for such a mechanism is seen in recent in vitro micromanipulation experiments. In this article we provide new perspective on existing theory and experimental data on DNA looping in vitro and in vivo. We elaborate on the connection between tension and a variety of other intracellular mechanical constraints including sequence specific curvature and supercoiling. In the process, we emphasize that the richness and versatility of DNA mechanics opens up a whole new paradigm of gene regulation to explore.

Modeling DNA loops using the theory of elasticity

An elastic rod model of a protein-bound DNA loop is adapted for application in multi-scale simulations of protein-DNA complexes. The classical Kirchhoff system of equations which describes the equilibrium structure of the elastic loop is modified to account for the intrinsic twist and curvature, anisotropic bending properties, and electrostatic charge of DNA. The effects of bending anisotropy and electrostatics are studied for the DNA loop clamped by the lac repressor protein. For two possible lengths of the loop, several topologically different conformations are predicted and extensively analyzed over the broad range of model parameters describing DNA bending and electrostatic properties. The scope and applications of the model in already accomplished and in future multi-scale studies of protein-DNA complexes are discussed.

Conformational Analysis of Stiff Chiral Polymers with End-Constraints

We present a Lie-group-theoretic method for the kinematic and dynamic analysis of chiral semi-flexible polymers with end constraints. The first is to determine the minimum energy conformations of semi-flexible polymers with end constraints, and the second is to perform normal mode analysis based on the determined minimum energy conformations. In this paper, we use concepts from the theory of Lie groups and principles of variational calculus to model such polymers as inextensible or extensible chiral elastic rods with coupling between twisting and bending stiffnesses, and/or between twisting and extension stiffnesses. This method is general enough to include any stiffness and chirality parameters in the context of elastic filament models with the quadratic elastic potential energy function. As an application of this formulation, the analysis of DNA conformations is discussed. We demonstrate our method with examples of DNA conformations in which topological properties such as writhe, twist, and linking number are calculated from the results of the proposed method. Given these minimum energy conformations, we describe how to perform the normal mode analysis. The results presented here build both on recent experimental work in which DNA mechanical properties have been measured, and theoretical work in which the mechanics of non-chiral elastic rods has been studied.

Scalable molecular dynamics with NAMD

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

A gamut of loops: meandering DNA

Nucleoprotein complexes comprising short DNA loops (150 base pairs or less) are involved in a wide variety of DNA transactions (e.g. transcription regulation, replication and recombination) in both prokaryotes and eukaryotes, and also can be useful in designing nanostructures. In these higher-order nucleoprotein complexes, proteins bound to spatially separated sites on a DNA interact with each other by looping out the relatively stiff intervening DNA. Recent technological developments have enabled determination of DNA trajectories in a few DNA-loop-containing regulatory complexes. Results show that, in a given system, a specific DNA trajectory is preferred over others.

Structural dynamics of the lac repressor-DNA complex revealed by a multiscale simulation

A multiscale simulation of a complex between the lac repressor protein (LacI) and a 107-bp-long DNA segment is reported. The complex between the repressor and two operator DNA segments is described by all-atom molecular dynamics; the size of the simulated system comprises either 226,000 or 314,000 atoms. The DNA loop connecting the operators is modeled as a continuous elastic ribbon, described mathematically by the nonlinear Kirchhoff differential equations with boundary conditions obtained from the coordinates of the terminal base pairs of each operator. The forces stemming from the looped DNA are included in the molecular dynamics simulations; the loop structure and the forces are continuously recomputed because the protein motions during the simulations shift the operators and the presumed termini of the loop. The simulations reveal the structural dynamics of the LacI-DNA complex in unprecedented detail. The multiple domains of LacI exhibit remarkable structural stability during the simulation, moving much like rigid bodies. LacI is shown to absorb the strain from the looped DNA mainly through its mobile DNA-binding head groups. Even with large fluctuating forces applied, the head groups tilt strongly and keep their grip on the operator DNA, while the remainder of the protein retains its V-shaped structure. A simulated opening of the cleft of LacI by 500-pN forces revealed the interactions responsible for locking LacI in the V-conformation.

Disruption of protein-mediated DNA looping by tension in the substrate DNA

Protein-mediated DNA looping is important in a variety of biological processes, including gene regulation and genetic transformation. Although the biochemistry of loop formation is well established, the mechanics of loop closure in a constrained cellular environment has received less attention. Recent single molecule measurements show that mechanical constraints have a significant impact on DNA looping and motivate the need for a more comprehensive characterization of the effects of tension. By modeling DNA as a wormlike chain, we calculate how continuous stretching of the substrate DNA affects the loop formation probability. We find that when the loop size is >100 bp, a tension of 500 fN can increase the time required for loop closure by two orders of magnitude. This force is small compared to the piconewton forces that are associated with RNA polymerases and other molecular motors, indicating that intracellular mechanical forces might affect transcriptional regulation. In contrast to existing theory, we find that for loops <200 bp, the effect of tension is partly dependent on the relative orientation of the DNA-binding domains in the linker protein. Our results provide perspective on recent DNA looping experiments and suggestions for future micromechanical studies.

Modelling DNA loops using continuum and statistical mechanics

The classical Kirchhoff elastic-rod model applied to DNA is extended to account for sequence-dependent intrinsic twist and curvature, anisotropic bending rigidity, electrostatic force interactions, and overdamped Brownian motion in a solvent. The zero-temperature equilibrium rod model is then applied to study the structural basis of the function of the lac repressor protein in the lac operon of Escherichia coli. The structure of a DNA loop induced by the clamping of two distant DNA operator sites by lac repressor is investigated and the optimal geometries for the loop of length 76 bp are predicted. Further, the mimicked binding of catabolite gene activator protein (CAP) inside the loop provides solutions that might explain the experimentally observed synergy in DNA binding between the two proteins. Finally, a combined Monte Carlo and Brownian dynamics solver for a worm-like chain model is described and a preliminary analysis of DNA loop-formation kinetics is presented.