Engineering and application of LacI mutants with stringent expressions

Optimal transcriptional regulatory circuits are expected to exhibit stringent control, maintaining silence in the absence of inducers while exhibiting a broad induction dynamic range upon the addition of effectors. In the Plac /LacI pair, the promoter of the lac operon in Escherichia coli is characterized by its leakiness, attributed to the moderate affinity of LacI for its operator target. In response to this limitation, the LacI regulatory protein underwent engineering to enhance its regulatory properties. The M7 mutant, carrying I79T and N246S mutations, resulted in the lac promoter displaying approximately 95% less leaky expression and a broader induction dynamic range compared to the wild-type LacI. An in-depth analysis of each mutation revealed distinct regulatory profiles. In contrast to the wild-type LacI, the M7 mutant exhibited a tighter binding to the operator sequence, as evidenced by surface plasmon resonance studies. Leveraging the capabilities of the M7 mutant, a high-value sugar biosensor was constructed. This biosensor facilitated the selection of mutant galactosidases with approximately a seven-fold improvement in specific activity for transgalactosylation. Consequently, this advancement enabled enhanced biosynthesis of galacto-oligosaccharides (GOS).

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

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

A tuneable genetic switch for tight control of tac promoters in Escherichia coli boosts expression of synthetic injectisomes

Biosafety of engineered bacteria as living therapeutics requires a tight regulation to control the specific delivery of protein effectors, maintaining minimum leakiness in the uninduced (OFF) state and efficient expression in the induced (ON) state. Here, we report a three repressors (3R) genetic circuit that tightly regulates the expression of multiple tac promoters (Ptac) integrated in the chromosome of E. coli and drives the expression of a complex type III secretion system injectisome for therapeutic protein delivery. The 3R genetic switch is based on the tetracycline repressor (TetR), the non-inducible lambda repressor cI (ind-) and a mutant lac repressor (LacIW220F ) with higher activity. The 3R switch was optimized with different protein translation and degradation signals that control the levels of LacIW220F . We demonstrate the ability of an optimized switch to fully repress the strong leakiness of the Ptac promoters in the OFF state while triggering their efficient activation in the ON state with anhydrotetracycline (aTc), an inducer suitable for in vivo use. The implementation of the optimized 3R switch in the engineered synthetic injector E. coli (SIEC) strain boosts expression of injectisomes upon aTc induction, while maintaining a silent OFF state that preserves normal growth in the absence of the inducer. Since Ptac is a commonly used promoter, the 3R switch may have multiple applications for tight control of protein expression in E. coli. In addition, the modularity of the 3R switch may enable its tuning for the control of Ptac promoters with different inducers.

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

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

Deep representation learning improves prediction of LacI-mediated transcriptional repression

Recent progress in DNA synthesis and sequencing technology has enabled systematic studies of protein function at a massive scale. We explore a deep mutational scanning study that measured the transcriptional repression function of 43,669 variants of the Escherichia coli LacI protein. We analyze structural and evolutionary aspects that relate to how the function of this protein is maintained, including an in-depth look at the C-terminal domain. We develop a deep neural network to predict transcriptional repression mediated by the lac repressor of Escherichia coli using experimental measurements of variant function. When measured across 10 separate training and validation splits using 5,009 single mutations of the lac repressor, our best-performing model achieved a median Pearson correlation of 0.79, exceeding any previous model. We demonstrate that deep representation learning approaches, first trained in an unsupervised manner across millions of diverse proteins, can be fine-tuned in a supervised fashion using lac repressor experimental datasets to more effectively predict a variant's effect on repression. These findings suggest a deep representation learning model may improve the prediction of other important properties of proteins.

The genotype-phenotype landscape of an allosteric protein

Allostery is a fundamental biophysical mechanism that underlies cellular sensing, signaling, and metabolism. Yet a quantitative understanding of allosteric genotype-phenotype relationships remains elusive. Here, we report the large-scale measurement of the genotype-phenotype landscape for an allosteric protein: the lac repressor from Escherichia coli, LacI. Using a method that combines long-read and short-read DNA sequencing, we quantitatively measure the dose-response curves for nearly 105 variants of the LacI genetic sensor. The resulting data provide a quantitative map of the effect of amino acid substitutions on LacI allostery and reveal systematic sequence-structure-function relationships. We find that in many cases, allosteric phenotypes can be quantitatively predicted with additive or neural-network models, but unpredictable changes also occur. For example, we were surprised to discover a new band-stop phenotype that challenges conventional models of allostery and that emerges from combinations of nearly silent amino acid substitutions.

Characterization of Gene Repression by Designed Transcription Activator-like Effector Dimer Proteins

Gene regulation by control of transcription initiation is a fundamental property of living cells. Much of our understanding of gene repression originated from studies of the Escherichia coli lac operon switch, in which DNA looping plays an essential role. To validate and generalize principles from lac for practical applications, we previously described artificial DNA looping driven by designed transcription activator-like effector dimer (TALED) proteins. Because TALE monomers bind the idealized symmetrical lac operator sequence in two orientations, our prior studies detected repression due to multiple DNA loops. We now quantitatively characterize gene repression in living E. coli by a collection of individual TALED loops with systematic loop length variation. Fitting of a thermodynamic model allows unequivocal demonstration of looping and comparison of the engineered TALED repression system with the natural lac repressor system.

Artificial Cells Capable of Long-Lived Protein Synthesis by Using Aptamer Grafted Polymer Hydrogel

Herein we report a new type of artificial cells capable of long-term protein expression and regulation. We constructed the artificial cells by grafting anti-His-tag aptamer into the polymer backbone of the hydrogel particles, and then immobilizing the His-tagged proteinaceous factors of the transcription and translation system into the hydrogel particles. Long-term protein expression for at least 16 days was achieved by continuously flowing feeding buffer through the artificial cells. The effect of various metal ions on the protein expression in the artificial cells was investigated. Utilizing the lac operator-repressor system, we could regulate the expression level of eGFP in the artificial cells by controlling the β-D-1-thiogalatopyranoside (IPTG) concentration in the feeding buffer. The artificial cells based on the aptamer grafted hydrogel provide a useful platform for gene circuit engineering, metabolic engineering, drug delivery, and biosensors.

Predictive shifts in free energy couple mutations to their phenotypic consequences

Mutation is a critical mechanism by which evolution explores the functional landscape of proteins. Despite our ability to experimentally inflict mutations at will, it remains difficult to link sequence-level perturbations to systems-level responses. Here, we present a framework centered on measuring changes in the free energy of the system to link individual mutations in an allosteric transcriptional repressor to the parameters which govern its response. We find that the energetic effects of the mutations can be categorized into several classes which have characteristic curves as a function of the inducer concentration. We experimentally test these diagnostic predictions using the well-characterized LacI repressor of Escherichia coli, probing several mutations in the DNA binding and inducer binding domains. We find that the change in gene expression due to a point mutation can be captured by modifying only the model parameters that describe the respective domain of the wild-type protein. These parameters appear to be insulated, with mutations in the DNA binding domain altering only the DNA affinity and those in the inducer binding domain altering only the allosteric parameters. Changing these subsets of parameters tunes the free energy of the system in a way that is concordant with theoretical expectations. Finally, we show that the induction profiles and resulting free energies associated with pairwise double mutants can be predicted with quantitative accuracy given knowledge of the single mutants, providing an avenue for identifying and quantifying epistatic interactions.

Sliding of a single lac repressor protein along DNA is tuned by DNA sequence and molecular switching

In any living cell, genome maintenance is carried out by DNA-binding proteins that recognize specific sequences among a vast amount of DNA. This includes fundamental processes such as DNA replication, DNA repair, and gene expression and regulation. Here, we study the mechanism of DNA target search by a single lac repressor protein (LacI) with ultrafast force-clamp spectroscopy, a sub-millisecond and few base-pair resolution technique based on laser tweezers. We measure 1D-diffusion of proteins on DNA at physiological salt concentrations with 20 bp resolution and find that sliding of LacI along DNA is sequence dependent. We show that only allosterically activated LacI slides along non-specific DNA sequences during target search, whereas the inhibited conformation does not support sliding and weakly interacts with DNA. Moreover, we find that LacI undergoes a load-dependent conformational change when it switches between sliding and strong binding to the target sequence. Our data reveal how DNA sequence and molecular switching regulate LacI target search process and provide a comprehensive model of facilitated diffusion for LacI.

Structure-guided approach to site-specific fluorophore labeling of the lac repressor LacI

The lactose operon repressor protein LacI has long served as a paradigm of the bacterial transcription factors. However, the mechanisms whereby LacI rapidly locates its cognate binding site on the bacterial chromosome are still elusive. Single-molecule fluorescence imaging approaches are well suited for the study of these mechanisms but rely on a functionally compatible fluorescence labeling of LacI. Particularly attractive for protein fluorescence labeling are synthetic fluorophores due to their small size and favorable photophysical characteristics. Synthetic fluorophores are often conjugated to natively occurring cysteine residues using maleimide chemistry. For a site-specific and functionally compatible labeling with maleimide fluorophores, the target protein often needs to be redesigned to remove unwanted native cysteines and to introduce cysteines at locations better suited for fluorophore attachment. Biochemical screens can then be employed to probe for the functional activity of the redesigned protein both before and after dye labeling. Here, we report a mutagenesis-based redesign of LacI to enable a functionally compatible labeling with maleimide fluorophores. To provide an easily accessible labeling site in LacI, we introduced a single cysteine residue at position 28 in the DNA-binding headpiece of LacI and replaced two native cysteines with alanines where derivatization with bulky substituents is known to compromise the protein’s activity. We find that the redesigned LacI retains a robust activity in vitro and in vivo, provided that the third native cysteine at position 281 is retained in LacI. In a total internal reflection microscopy assay, we observed individual Cy3-labeled LacI molecules bound to immobilized DNA harboring the cognate O₁ operator sequence, indicating that the dye-labeled LacI is functionally active. We have thus been able to generate a functional fluorescently labeled LacI that can be used to unravel mechanistic details of LacI target search at the single molecule level.

Bacterial gene control by DNA looping using engineered dimeric transcription activator like effector (TALE) proteins

Genetic switches must alternate between states whose probabilities are dependent on regulatory signals. Classical examples of transcriptional control in bacteria depend on repressive DNA loops anchored by proteins whose structures are sensitive to small molecule inducers or co-repressors. We are interested in exploiting these natural principles to engineer artificial switches for transcriptional control of bacterial genes. Here, we implement designed homodimeric DNA looping proteins ('Transcription Activator-Like Effector Dimers'; TALEDs) for this purpose in living bacteria. Using well-studied FKBP dimerization domains, we build switches that mimic regulatory characteristics of classical Escherichia coli lactose, galactose and tryptophan operon promoters, including induction or co-repression by small molecules. Engineered DNA looping using TALEDs is thus a new approach to tuning gene expression in bacteria. Similar principles may also be applicable for gene control in eukaryotes.

Seeing red; the development of pON.mCherry, a broad-host range constitutive expression plasmid for Gram-negative bacteria

The development of plasmid-mediated gene expression control in bacteria revolutionized the field of bacteriology. Many of these expression control systems rely on the addition of small molecules, generally metabolites or non-metabolized analogs thereof, to the growth medium to induce expression of the genes of interest. The paradigmatic example of an expression control system is the lac system from Escherichia coli, which typically relies on the Ptac promoter and the Lac repressor, LacI. In many cases, however, constitutive gene expression is desired, and other experimental approaches require the coordinated control of multiple genes. While multiple systems have been developed for use in E. coli and its close relatives, the utility and/or functionality of these tools does not always translate to other species. For example, for the Gram-negative pathogen, Legionella pneumophila, a causative agent of Legionnaires' Disease, the aforementioned Ptac system represents the only well-established expression control system. In order to enhance the tools available to study bacterial gene expression in L. pneumophila, we developed a plasmid, pON.mCherry, which confers constitutive gene expression from a mutagenized LacI binding site. We demonstrate that pON.mCherry neither interferes with other plasmids harboring an intact LacI-Ptac expression system nor alters the growth of Legionella species during intracellular growth. Furthermore, the broad-host range plasmid backbone of pON.mCherry allows constitutive gene expression in a wide variety of Gram-negative bacterial species, making pON.mCherry a useful tool for the greater research community.

Solution NMR investigation of the response of the lactose repressor core domain dimer to hydrostatic pressure

Previous investigations of the sensitivity of the lac repressor to high-hydrostatic pressure have led to varying conclusions. Here high-pressure solution NMR spectroscopy is used to provide an atomic level view of the pressure induced structural transition of the lactose repressor regulatory domain (LacI* RD) bound to the ligand IPTG. As the pressure is raised from ambient to 3kbar the native state of the protein is converted to a partially unfolded form. Estimates of rotational correlation times using transverse optimized relaxation indicates that a monomeric state is never reached and that the predominate form of the LacI* RD is dimeric throughout this pressure change. Spectral analysis suggests that the pressure-induced transition is localized and is associated with a volume change of approximately -115mlmol^-1 and an average pressure dependent change in compressibility of approximately 30mlmol^-1kbar^-1. In addition, a subset of resonances emerge at high-pressures indicating the presence of a non-native but folded alternate state.

Anomalous Dynamics of Water Confined in Protein-Protein and Protein-DNA Interfaces

Confined water often exhibits anomalous properties not observable in the bulk phase. Although water in hydrophobic confinement has been the focus of intense investigation, the behavior of water confined between hydrophilic surfaces, which are more frequently found in biological systems, has not been fully explored. Here, we investigate using molecular dynamics simulations dynamical properties of the water confined in hydrophilic protein-protein and protein-DNA interfaces. We find that the interfacial water exhibits glassy slow relaxations even at 300 K. In particular, the rotational dynamics show a logarithmic decay that was observed in glass-forming liquids at deeply supercooled states. We argue that such slow water dynamics are indeed induced by the hydrophilic binding surfaces, which is in opposition to the picture that the hydration water slaves protein motions. Our results will significantly impact the view on the role of water in biomolecular interactions.

Dissecting the stochastic transcription initiation process in live Escherichia coli

We investigate the hypothesis that, in Escherichia coli, while the concentration of RNA polymerases differs in different growth conditions, the fraction of RNA polymerases free for transcription remains approximately constant within a certain range of these conditions. After establishing this, we apply a standard model-fitting procedure to fully characterize the in vivo kinetics of the rate-limiting steps in transcription initiation of the Plac/ara-1 promoter from distributions of intervals between transcription events in cells with different RNA polymerase concentrations. We find that, under full induction, the closed complex lasts ∼788 s while subsequent steps last ∼193 s, on average. We then establish that the closed complex formation usually occurs multiple times prior to each successful initiation event. Furthermore, the promoter intermittently switches to an inactive state that, on average, lasts ∼87 s. This is shown to arise from the intermittent repression of the promoter by LacI. The methods employed here should be of use to resolve the rate-limiting steps governing the in vivo dynamics of initiation of prokaryotic promoters, similar to established steady-state assays to resolve the in vitro dynamics.

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

Spec-seq: determining protein-DNA-binding specificity by sequencing

The specificity of protein-DNA interactions can be determined directly by sequencing the bound and unbound fractions in a standard binding reaction. The procedure is easy and inexpensive, and the accuracy can be high for thousands of sequences assayed in parallel. From the measurements, simple models of specificity, such as position weight matrices, can be assessed for their accuracy and more complex models developed if useful. Those may provide more accurate predictions of in vivo binding sites and can help us to understand the details of recognition. As an example, we demonstrate new information gained about the binding of lac repressor. One can apply the same method to combinations of factors that bind simultaneously to a single DNA and determine both the specificity of the individual factors and the cooperativity between them.

What controls DNA looping?

The looping of DNA provides a means of communication between sequentially distant genomic sites that operate in tandem to express, copy, and repair the information encoded in the DNA base sequence. The short loops implicated in the expression of bacterial genes suggest that molecular factors other than the naturally stiff double helix are involved in bringing the interacting sites into close spatial proximity. New computational techniques that take direct account of the three-dimensional structures and fluctuations of protein and DNA allow us to examine the likely means of enhancing such communication. Here, we describe the application of these approaches to the looping of a 92 base-pair DNA segment between the headpieces of the tetrameric Escherichia coli Lac repressor protein. The distortions of the double helix induced by a second protein--the nonspecific nucleoid protein HU--increase the computed likelihood of looping by several orders of magnitude over that of DNA alone. Large-scale deformations of the repressor, sequence-dependent features in the DNA loop, and deformability of the DNA operators also enhance looping, although to lesser degrees. The correspondence between the predicted looping propensities and the ease of looping derived from gene-expression and single-molecule measurements lends credence to the derived structural picture.

Design and characterization of a dual-mode promoter with activation and repression capability for tuning gene expression in yeast

Modularity in controlling gene expression artificially is becoming an essential aspect of synthetic biology. Artificial transcriptional control of gene expression is one of the most well-developed methods for the design of novel synthetic regulatory networks. Such networks are intended to help understand natural cellular phenomena and to enable new biotechnological applications. Promoter sequence manipulation with cis-regulatory elements is a key approach to control gene expression transcriptionally. Here, we have designed a promoter that can be both activated and repressed, as a contribution to the library of synthetic biological 'parts'. Starting with the minimal cytochrome C (minCYC) promoter in yeast, we incorporated five steroid hormone responsive elements (SHREs) and one lac operator site, respectively, upstream and downstream of the TATA box. This allows activation through the testosterone-responsive androgen receptor, and repression through the LacI repressor. Exposure to varying concentrations of testosterone (to vary activation) and IPTG (to vary repression) demonstrated the ability to tune the promoter's output curve over a wide range. By integrating activating and repressing signals, the promoter permits a useful form of signal integration, and we are optimistic that it will serve as a component in future regulatory networks, including feedback controllers.

Protein distributions from a stochastic model of the lac operon of E. coli with DNA looping: analytical solution and comparison with experiments

Although noisy gene expression is widely accepted, its mechanisms are subjects of debate, stimulated largely by single-molecule experiments. This work is concerned with one such study, in which Choi et al., 2008, obtained real-time data and distributions of Lac permease in E. coli. They observed small and large protein bursts in strains with and without auxiliary operators. They also estimated the size and frequency of these bursts, but these were based on a stochastic model of a constitutive promoter. Here, we formulate and solve a stochastic model accounting for the existence of auxiliary operators and DNA loops. We find that DNA loop formation is so fast that small bursts are averaged out, making it impossible to extract their size and frequency from the data. In contrast, we can extract not only the size and frequency of the large bursts, but also the fraction of proteins derived from them. Finally, the proteins follow not the negative binomial distribution, but a mixture of two distributions, which reflect the existence of proteins derived from small and large bursts.

Lac repressor mediated DNA looping: Monte Carlo simulation of constrained DNA molecules complemented with current experimental results

Tethered particle motion (TPM) experiments can be used to detect time-resolved loop formation in a single DNA molecule by measuring changes in the length of a DNA tether. Interpretation of such experiments is greatly aided by computer simulations of DNA looping which allow one to analyze the structure of the looped DNA and estimate DNA-protein binding constants specific for the loop formation process. We here present a new Monte Carlo scheme for accurate simulation of DNA configurations subject to geometric constraints and apply this method to Lac repressor mediated DNA looping, comparing the simulation results with new experimental data obtained by the TPM technique. Our simulations, taking into account the details of attachment of DNA ends and fluctuations of the looped subsegment of the DNA, reveal the origin of the double-peaked distribution of RMS values observed by TPM experiments by showing that the average RMS value for anti-parallel loop types is smaller than that of parallel loop types. The simulations also reveal that the looping probabilities for the anti-parallel loop types are significantly higher than those of the parallel loop types, even for loops of length 600 and 900 base pairs, and that the correct proportion between the heights of the peaks in the distribution can only be attained when loops with flexible Lac repressor conformation are taken into account. Comparison of the in silico and in vitro results yields estimates for the dissociation constants characterizing the binding affinity between O1 and Oid DNA operators and the dimeric arms of the Lac repressor.

Direct measurement of transcription factor dissociation excludes a simple operator occupancy model for gene regulation

Transcription factors mediate gene regulation by site-specific binding to chromosomal operators. It is commonly assumed that the level of repression is determined solely by the equilibrium binding of a repressor to its operator. However, this assumption has not been possible to test in living cells. Here we have developed a single-molecule chase assay to measure how long an individual transcription factor molecule remains bound at a specific chromosomal operator site. We find that the lac repressor dimer stays bound on average 5 min at the native lac operator in Escherichia coli and that a stronger operator results in a slower dissociation rate but a similar association rate. Our findings do not support the simple equilibrium model. The discrepancy with this model can, for example, be accounted for by considering that transcription initiation drives the system out of equilibrium. Such effects need to be considered when predicting gene activity from transcription factor binding strengths.

Rheostats and toggle switches for modulating protein function

The millions of protein sequences generated by genomics are expected to transform protein engineering and personalized medicine. To achieve these goals, tools for predicting outcomes of amino acid changes must be improved. Currently, advances are hampered by insufficient experimental data about nonconserved amino acid positions. Since the property “nonconserved” is identified using a sequence alignment, we designed experiments to recapitulate that context: Mutagenesis and functional characterization was carried out in 15 LacI/GalR homologs (rows) at 12 nonconserved positions (columns). Multiple substitutions were made at each position, to reveal how various amino acids of a nonconserved column were tolerated in each protein row. Results showed that amino acid preferences of nonconserved positions were highly context-dependent, had few correlations with physico-chemical similarities, and were not predictable from their occurrence in natural LacI/GalR sequences. Further, unlike the “toggle switch” behaviors of conserved positions, substitutions at nonconserved positions could be rank-ordered to show a “rheostatic”, progressive effect on function that spanned several orders of magnitude. Comparisons to various sequence analyses suggested that conserved and strongly co-evolving positions act as functional toggles, whereas other important, nonconserved positions serve as rheostats for modifying protein function. Both the presence of rheostat positions and the sequence analysis strategy appear to be generalizable to other protein families and should be considered when engineering protein modifications or predicting the impact of protein polymorphisms.

Transcription-factor binding and sliding on DNA studied using micro- and macroscopic models

Transcription factors search for specific operator sequences by alternating rounds of 3D diffusion with rounds of 1D diffusion (sliding) along the DNA. The details of such sliding have largely been beyond direct experimental observation. For this purpose we devised an analytical formulation of umbrella sampling along a helical coordinate, and from extensive and fully atomistic simulations we quantified the free-energy landscapes that underlie the sliding dynamics and dissociation kinetics for the LacI dimer. The resulting potential of mean force distributions show a fine structure with an amplitude of 1 k(B)T for sliding and 12 k(B)T for dissociation. Based on the free-energy calculations the repressor slides in close contact with DNA for 8 bp on average before making a microscopic dissociation. By combining the microscopic molecular-dynamics calculations with Brownian simulation including rotational diffusion from the microscopically dissociated state we estimate a macroscopic residence time of 48 ms at the same DNA segment and an in vitro sliding distance of 240 bp. The sliding distance is in agreement with previous in vitro sliding-length estimates. The in vitro prediction for the macroscopic residence time also compares favorably to what we measure by single-molecule imaging of nonspecifically bound fluorescently labeled LacI in living cells. The investigation adds to our understanding of transcription-factor search kinetics and connects the macro-/mesoscopic rate constants to the microscopic dynamics.

Using synthetic biology to distinguish and overcome regulatory and functional barriers related to nitrogen fixation

Biological nitrogen fixation is a complex process requiring multiple genes working in concert. To date, the Klebsiella pneumoniae nif gene cluster, divided into seven operons, is one of the most studied systems. Its nitrogen fixation capacity is subject to complex cascade regulation and physiological limitations. In this report, the entire K. pneumoniae nif gene cluster was reassembled as operon-based BioBrick parts in Escherichia coli. It provided ~100% activity of native K. pneumoniae system. Based on the expression levels of these BioBrick parts, a T7 RNA polymerase-LacI expression system was used to replace the σ(54)-dependent promoters located upstream of nif operons. Expression patterns of nif operons were critical for the maximum activity of the recombinant system. By mimicking these expression levels with variable-strength T7-dependent promoters, ~42% of the nitrogenase activity of the σ(54)-dependent nif system was achieved in E. coli. When the newly constructed T7-dependent nif system was challenged with different genetic and physiological conditions, it bypassed the original complex regulatory circuits, with minor physiological limitations. Therefore, we have successfully replaced the nif regulatory elements with a simple expression system that may provide the first step for further research of introducing nif genes into eukaryotic organelles, which has considerable potentials in agro-biotechnology.

A single mutation in the core domain of the lac repressor reduces leakiness

BACKGROUND

The lac operon provides cells with the ability to switch from glucose to lactose metabolism precisely when necessary. This metabolic switch is mediated by the lac repressor (LacI), which in the absence of lactose binds to the operator DNA sequence to inhibit transcription. Allosteric rearrangements triggered by binding of the lactose isomer allolactose to the core domain of the repressor impede DNA binding and lift repression. In Nature, the ability to detect and respond to environmental conditions comes at the cost of the encoded enzymes being constitutively expressed at low levels. The readily-switched regulation provided by LacI has resulted in its widespread use for protein overexpression, and its applications in molecular biology represent early examples of synthetic biology. However, the leakiness of LacI that is essential for the natural function of the lac operon leads to an increased energetic burden, and potentially toxicity, in heterologous protein production.

RESULTS

Analysis of the features that confer promiscuity to the inducer-binding site of LacI identified tryptophan 220 as a target for saturation mutagenesis. We found that phenylalanine (similarly to tryptophan) affords a functional repressor that is still responsive to IPTG. Characterisation of the W220F mutant, LacIWF, by measuring the time dependence of GFP production at different IPTG concentrations and at various incubation temperatures showed a 10-fold reduction in leakiness and no decrease in GFP production. Cells harbouring a cytotoxic protein under regulatory control of LacIWF showed no decrease in viability in the early phases of cell growth. Changes in responsiveness to IPTG observed in vivo are supported by the thermal shift assay behaviour of purified LacIWF with IPTG and operator DNA.

CONCLUSIONS

In LacI, long-range communications are responsible for the transmission of the signal from the inducer binding site to the DNA binding domain and our results are consistent with the involvement of position 220 in modulating these. The mutation of this single tryptophan residue to phenylalanine generated an enhanced repressor with a 10-fold decrease in leakiness. By minimising the energetic burden and cytotoxicity caused by leakiness, LacIWF constitutes a useful switch for protein overproduction and synthetic biology.

Theoretical and experimental analysis of the forced LacI-AraC oscillator with a minimal gene regulatory model

Oscillatory dynamics have been observed in multiple cellular functions and synthetic constructs; and here, we study the behavior of a synthetic oscillator under temporal perturbations. We use a minimal model, involving a single transcription factor with delayed self-repression and enzymatic degradation, together with a first-order perturbative approach, to derive an analytical expression for the power spectrum of the system, which characterizes its response to external forces and molecular noise. Experimentally, we force and monitor the dynamics of the LacI-AraC oscillator in single cells during long time intervals by constructing a microfluidics device. Pulse dynamics of IPTG with different periods serve to perturb this system. Due to the resonance of the system, we predict theoretically and confirm experimentally the dependence on the forcing frequency of the variability in gene expression with time and the synchronization of the population to the input signal. The reported results show that the engineering of gene circuits can provide test cases for dynamical models, which could be further exploited in synthetic biology.

The effect of LacI autoregulation on the performance of the lactose utilization system in Escherichia coli

The lactose operon of Escherichia coli is a paradigm system for quantitative understanding of gene regulation in prokaryotes. Yet, none of the many mathematical models built so far to study the dynamics of this system considered the fact that the Lac repressor regulates its own transcription by forming a transcriptional roadblock at the O3 operator site. Here we study the effect of autoregulation on intracellular LacI levels and also show that cAMP-CRP binding does not affect the efficiency of autoregulation. We built a mathematical model to study the role of LacI autoregulation in the lactose utilization system. Previously, it has been argued that negative autoregulation can significantly reduce noise as well as increase the speed of response. We show that the particular molecular mechanism, a transcriptional roadblock, used to achieve self-repression in the lac system does neither. Instead, LacI autoregulation balances two opposing states, one that allows quicker response to smaller pulses of external lactose, and the other that minimizes production costs in the absence of lactose.

Sliding and target location of DNA-binding proteins: an NMR view of the lac repressor system

In non-specific lac headpiece-DNA complexes selective NMR line broadening is observed that strongly depends on length and composition of the DNA fragments. This broadening involves amide protons found in the non-specific lac-DNA structure to be interacting with the DNA phosphate backbone, and can be ascribed to DNA sliding of the protein along the DNA. This NMR exchange broadening has been used to estimate the 1D diffusion constant for sliding along non-specific DNA. The observed 1D diffusion constant of 4×10(-12) cm(2)/s is two orders of magnitude smaller than derived from previous kinetic experiments, but falls in the range of values determined more recently using single molecule methods. This strongly supports the notion that sliding could play at most a minor role in the association kinetics of binding of lac repressor to lac operator and that other processes such as hopping and intersegment transfer contribute to facilitate the DNA recognition process.

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.

Comparative genomics of CytR, an unusual member of the LacI family of transcription factors

CytR is a transcription regulator from the LacI family, present in some gamma-proteobacteria including Escherichia coli and known not only for its cellular role, control of transport and utilization of nucleosides, but for a number of unusual structural properties. The present study addressed three related problems: structure of CytR-binding sites and motifs, their evolutionary conservation, and identification of new members of the CytR regulon. While the majority of CytR-binding sites are imperfect inverted repeats situated between binding sites for another transcription factor, CRP, other architectures were observed, in particular, direct repeats. While the similarity between sites for different genes in one genome is rather low, and hence the consensus motif is weak, there is high conservation of orthologous sites in different genomes (mainly in the Enterobacteriales) arguing for the presence of specific CytR-DNA contacts. On larger evolutionary distances candidate CytR sites may migrate but the approximate distance between flanking CRP sites tends to be conserved, which demonstrates that the overall structure of the CRP-CytR-DNA complex is gene-specific. The analysis yielded candidate CytR-binding sites for orthologs of known regulon members in less studied genomes of the Enterobacteriales and Vibrionales and identified a new candidate member of the CytR regulon, encoding a transporter named NupT (YcdZ).

Adaptive evolution of the lactose utilization network in experimentally evolved populations of Escherichia coli

Adaptation to novel environments is often associated with changes in gene regulation. Nevertheless, few studies have been able both to identify the genetic basis of changes in regulation and to demonstrate why these changes are beneficial. To this end, we have focused on understanding both how and why the lactose utilization network has evolved in replicate populations of Escherichia coli. We found that lac operon regulation became strikingly variable, including changes in the mode of environmental response (bimodal, graded, and constitutive), sensitivity to inducer concentration, and maximum expression level. In addition, some classes of regulatory change were enriched in specific selective environments. Sequencing of evolved clones, combined with reconstruction of individual mutations in the ancestral background, identified mutations within the lac operon that recapitulate many of the evolved regulatory changes. These mutations conferred fitness benefits in environments containing lactose, indicating that the regulatory changes are adaptive. The same mutations conferred different fitness effects when present in an evolved clone, indicating that interactions between the lac operon and other evolved mutations also contribute to fitness. Similarly, changes in lac regulation not explained by lac operon mutations also point to important interactions with other evolved mutations. Together these results underline how dynamic regulatory interactions can be, in this case evolving through mutations both within and external to the canonical lactose utilization network.

[Machine learning study of DNA binding by transcription factors from the LacI family]

We studied 1372 LacI-family transcription factors and their 4484 DNA binding sites using machine learning algorithms and feature selection techniques. The Naive Bayes classifier and Logistic Regression were used to predict binding sites given transcription factor sequences and to classify factor-site pairs on binding and non-binding ones. Prediction accuracy was estimated using 10-fold cross-validation. Experiments showed that the best prediction of nucleotide densities at selected site positions is obtained using only a few key protein sequence positions. These positions are stably selected by the forward feature selection based on the mutual information of factor-site position pairs.

A tale of two repressors

Few proteins have had such a strong impact on a field, as the lac repressor and λ repressor have had in Molecular Biology in bacteria. The genes required for lactose utilization are negatively regulated; the lac repressor binds to an upstream operator blocking the transcription of the enzymes necessary for lactose utilization. A similar switch regulates the virus life cycle; λ repressor binds to an operator site and blocks transcription of the phage genes necessary for lytic development. It is now 50 years since Jacob and Monod first proposed a model for gene regulation, which survives essentially unchanged in contemporary textbooks. Jacob, F. & Monod, J. (1961). Genetic regulatory mechanisms in the synthesis of proteins. J. Mol. Biol. 3, 318-356. This model provides a cogent depiction of how a set of genes can be coordinately transcribed in response to environmental conditions and regulates metabolic events in the cell. A historical perspective that illustrates the role these two repressor molecules played and their contribution to our understanding of gene regulation is presented.

Monitoring DNA binding to Escherichia coli lactose repressor using quartz crystal microbalance with dissipation

Lactose repressor protein (LacI) functions as a negative transcription regulator in Escherichia coli by binding to the operator DNA sequence. Our understanding of the immobilized LacI function and the effect of ligand binding on the conformation of LacI-DNA complexes remains poorly understood. Here, we have examined the difference in functionality of wild-type and mutant LacI binding to the target DNA using quartz crystal microbalance with dissipation (QCM-D). To direct the orientation of LacI binding to the gold surface, residue 334 was substituted with cysteine (T334C) to generate a sulfur-gold linkage. Position 334 is located on the surface opposite the DNA-binding domain and remote from the site for inducer binding. With T334C immobilized on the gold surface, our sensors successfully detect operator binding as well as the release of the operator DNA from the repressor in the presence of inducer isopropyl-β-D-thiogalactoside (IPTG). Besides the natural operator DNA sequence (O(1)), a symmetric high-affinity DNA sequence (O(sym)), and a non-specific DNA (O(ns)) sequence with low affinity were also used. In addition, the impact of anti-inducer o-nitrophenyl-beta-d-fucoside (ONPF), which stabilizes LacI operator binding, was examined. The results from immobilized mutant LacI are in good agreement with known solution parameters for LacI-ligand interactions, demonstrating that QCM-D provides a rapid and efficient measurement of DNA binding and impact of ligands upon binding for this complex oligomeric protein.

The in vivo kinetics of RNA polymerase II elongation during co-transcriptional splicing

Kinetic analysis shows that RNA polymerase elongation kinetics are not modulated by co-transcriptional splicing and that post-transcriptional splicing can proceed at the site of transcription without the presence of the polymerase.

RNA processing events that take place on the transcribed pre-mRNA include capping, splicing, editing, 3′ processing, and polyadenylation. Most of these processes occur co-transcriptionally while the RNA polymerase II (Pol II) enzyme is engaged in transcriptional elongation. How Pol II elongation rates are influenced by splicing is not well understood. We generated a family of inducible gene constructs containing increasing numbers of introns and exons, which were stably integrated in human cells to serve as actively transcribing gene loci. By monitoring the association of the transcription and splicing machineries on these genes in vivo, we showed that only U1 snRNP localized to the intronless gene, consistent with a splicing-independent role for U1 snRNP in transcription. In contrast, all snRNPs accumulated on intron-containing genes, and increasing the number of introns increased the amount of spliceosome components recruited. This indicates that nascent RNA can assemble multiple spliceosomes simultaneously. Kinetic measurements of Pol II elongation in vivo, Pol II ChIP, as well as use of Spliceostatin and Meayamycin splicing inhibitors showed that polymerase elongation rates were uncoupled from ongoing splicing. This study shows that transcription elongation kinetics proceed independently of splicing at the model genes studied here. Surprisingly, retention of polyadenylated mRNA was detected at the transcription site after transcription termination. This suggests that the polymerase is released from chromatin prior to the completion of splicing, and the pre-mRNA is post-transcriptionally processed while still tethered to chromatin near the gene end.

The pre-mRNA emerging from RNA polymerase II during eukaryotic transcription undergoes a series of processing events. These include 5′-capping, intron excision and exon ligation during splicing, 3′-end processing, and polyadenylation. Processing events occur co-transcriptionally, meaning that a variety of enzymes assemble on the pre-mRNA while the polymerase is still engaged in transcription. The concept of co-transcriptional mRNA processing raises questions about the possible coupling between the transcribing polymerase and the processing machineries. Here we examine how the co-transcriptional assembly of the splicing machinery (the spliceosome) might affect the elongation kinetics of the RNA polymerase. Using live-cell microscopy, we followed the kinetics of transcription of genes containing increasing numbers of introns and measured the recruitment of transcription and splicing factors. Surprisingly, a sub-set of splicing factors was recruited to an intronless gene, implying that there is a polymerase-coupled scanning mechanism for intronic sequences. There was no difference in polymerase elongation rates on genes with or without introns, suggesting that the spliceosome does not modulate elongation kinetics. Experiments including inhibition of splicing or transcription, together with stochastic computational simulation, demonstrated that pre-mRNAs can be retained on the gene when polymerase termination precedes completion of splicing. Altogether we show that polymerase elongation kinetics are not affected by splicing events on the emerging pre-mRNA, that increased splicing leads to more splicing factors being recruited to the mRNA, and that post-transcriptional splicing can proceed at the site of transcription in the absence of the polymerase.

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.

Binding of lac repressor-GFP fusion protein to lac operator sites inserted in the tobacco chloroplast genome examined by chromatin immunoprecipitation

Chromatin immunoprecipitation (ChIP) has been used to detect binding of DNA-binding proteins to sites in nuclear and mitochondrial genomes. Here, we describe a method for detecting protein-binding sites on chloroplast DNA, using modifications to the nuclear ChIP procedures. The method was developed using the lac operator (lacO)/lac repressor (LacI) system from Escherichia coli. The lacO sequences were integrated into a single site between the rbcL and accD genes in tobacco plastid DNA and homoplasmic transplastomic plants were crossed with transgenic tobacco plants expressing a nuclear-encoded plastid-targeted GFP-LacI fusion protein. In the progeny, the GFP-LacI fusion protein could be visualized in living tissues using confocal microscopy, and was found to co-localize with plastid nucleoids. Isolated chloroplasts from the lacO/GFP-LacI plants were lysed, treated with micrococcal nuclease to digest the DNA to fragments of approximately 600 bp and incubated with antibodies to GFP and protein A-Sepharose. PCR analysis on DNA extracted from the immunoprecipitate demonstrated IPTG (isopropylthiogalactoside)-sensitive binding of GFP-LacI to lacO. Binding of GFP-LacI to endogenous sites in plastid DNA showing sequence similarity to lacO was also detected, but required reversible cross-linking with formaldehyde. This may provide a general method for the detection of binding sites on plastid DNA for specific proteins.

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.

Why Hofmeister effects of many salts favor protein folding but not DNA helix formation

The majority ( approximately 70%) of surface buried in protein folding is hydrocarbon, whereas in DNA helix formation, the majority ( approximately 65%) of surface buried is relatively polar nitrogen and oxygen. Our previous quantification of salt exclusion from hydrocarbon (C) accessible surface area (ASA) and accumulation at amide nitrogen (N) and oxygen (O) ASA leads to a prediction of very different Hofmeister effects on processes that bury mostly polar (N, O) surface compared to the range of effects commonly observed for processes that bury mainly nonpolar (C) surface, e.g., micelle formation and protein folding. Here we quantify the effects of salts on folding of the monomeric DNA binding domain (DBD) of lac repressor (lac DBD) and on formation of an oligomeric DNA duplex. In accord with this prediction, no salt investigated has a stabilizing Hofmeister effect on DNA helix formation. Our ASA-based analyses of model compound data and estimates of the surface area buried in protein folding and DNA helix formation allow us to predict Hofmeister effects on these processes. We observe semiquantitative to quantitative agreement between these predictions and the experimental values, obtained from a novel separation of coulombic and Hofmeister effects. Possible explanations of deviations, including salt-dependent unfolded ensembles and interactions with other types of surface, are discussed.

Femtonewton entropic forces can control the formation of protein-mediated DNA loops

We show that minuscule entropic forces, on the order of 100 fN, can prevent the formation of DNA loops-a ubiquitous means of regulating the expression of genes. We observe a tenfold decrease in the rate of LacI-mediated DNA loop formation when a tension of 200 fN is applied to the substrate DNA, biasing the thermal fluctuations that drive loop formation and breakdown events. Conversely, once looped, the DNA-protein complex is insensitive to applied force. Our measurements are in excellent agreement with a simple polymer model of loop formation in DNA, and show that an antiparallel topology is the preferred LacI-DNA loop conformation for a generic loop-forming construct.

Mathematical biology modules based on modern molecular biology and modern discrete mathematics

We describe an ongoing collaborative curriculum materials development project between Sweet Briar College and Western Michigan University, with support from the National Science Foundation. We present a collection of modules under development that can be used in existing mathematics and biology courses, and we address a critical national need to introduce students to mathematical methods beyond the interface of biology with calculus. Based on ongoing research, and designed to use the project-based-learning approach, the modules highlight applications of modern discrete mathematics and algebraic statistics to pressing problems in molecular biology. For the majority of projects, calculus is not a required prerequisite and, due to the modest amount of mathematical background needed for some of the modules, the materials can be used for an early introduction to mathematical modeling. At the same time, most modules are connected with topics in linear and abstract algebra, algebraic geometry, and probability, and they can be used as meaningful applied introductions into the relevant advanced-level mathematics courses. Open-source software is used to facilitate the relevant computations. As a detailed example, we outline a module that focuses on Boolean models of the lac operon network.

Towards evolving a better repressor

Transcriptional regulation is an essential component of all metabolic pathways. At the most basic level, a protein binds to a particular DNA sequence (operator) on the genome and either positively or negatively alters the level of transcription. Together, the protein and its operator form an epigenetic switch that regulates gene expression. In an effort to produce a 'better' switch, we have discovered novel facets of the lac operon that are responsible for optimal functionality. We have uncovered a relationship between operator binding affinity and inducibility and demonstrated that the operator DNA is not a passive component of a genetic switch; it is responsible for establishing binding affinity, specificity as well as translational efficiency. In addition, an operator's directionality can indirectly affect gene expression. Unraveling the basic properties of this classical epigenetic switch demonstrates that multiple factors must be optimized in designing a better switch.

Structure-based predictive models for allosteric hot spots

In allostery, a binding event at one site in a protein modulates the behavior of a distant site. Identifying residues that relay the signal between sites remains a challenge. We have developed predictive models using support-vector machines, a widely used machine-learning method. The training data set consisted of residues classified as either hotspots or non-hotspots based on experimental characterization of point mutations from a diverse set of allosteric proteins. Each residue had an associated set of calculated features. Two sets of features were used, one consisting of dynamical, structural, network, and informatic measures, and another of structural measures defined by Daily and Gray [1]. The resulting models performed well on an independent data set consisting of hotspots and non-hotspots from five allosteric proteins. For the independent data set, our top 10 models using Feature Set 1 recalled 68–81% of known hotspots, and among total hotspot predictions, 58–67% were actual hotspots. Hence, these models have precision P = 58–67% and recall R = 68–81%. The corresponding models for Feature Set 2 had P = 55–59% and R = 81–92%. We combined the features from each set that produced models with optimal predictive performance. The top 10 models using this hybrid feature set had R = 73–81% and P = 64–71%, the best overall performance of any of the sets of models. Our methods identified hotspots in structural regions of known allosteric significance. Moreover, our predicted hotspots form a network of contiguous residues in the interior of the structures, in agreement with previous work. In conclusion, we have developed models that discriminate between known allosteric hotspots and non-hotspots with high accuracy and sensitivity. Moreover, the pattern of predicted hotspots corresponds to known functional motifs implicated in allostery, and is consistent with previous work describing sparse networks of allosterically important residues.

Allostery is the process whereby a molecule binds to one site in a protein and alters the function of a distant site. This phenomenon is ubiquitous, as proteins frequently must adapt their behavior to changes in the cellular milieu. The mechanism(s) underlying allostery remains incompletely understood. In particular, predictive models are needed that distinguish amino-acid residues that are critical to allostery, or “hotspots”, from non-hotspots. Here we have used data-mining approaches to infer rules that distinguish hotspots from non-hotspots. Starting with a data set of known hotspot and non-hotspot residues from a diverse set of allosteric proteins, the training data set, we applied machine learning to this data to “learn” models, or sets of rules, for distinguishing hotspots and non-hotspots by inferring associations between the classification (hotspot or non-hotspot) and an associated set of calculated attributes. Many models that showed the highest predictive power on the training data also exhibited high accuracy and sensitivity when applied to an independent data set. Moreover, the pattern of predicted hotspots in the proteins we studied was consistent with known structure/function relationships and previous work suggesting that a network of essential residues mediates the allosteric transition.

Analysis of in-vivo LacR-mediated gene repression based on the mechanics of DNA looping

Interactions of E. coli lac repressor (LacR) with a pair of operator sites on the same DNA molecule can lead to the formation of looped nucleoprotein complexes both in vitro and in vivo. As a major paradigm for loop-mediated gene regulation, parameters such as operator affinity and spacing, repressor concentration, and DNA bending induced by specific or non-specific DNA-binding proteins (e.g., HU), have been examined extensively. However, a complete and rigorous model that integrates all of these aspects in a systematic and quantitative treatment of experimental data has not been available. Applying our recent statistical-mechanical theory for DNA looping, we calculated repression as a function of operator spacing (58-156 bp) from first principles and obtained excellent agreement with independent sets of in-vivo data. The results suggest that a linear extended, as opposed to a closed v-shaped, LacR conformation is the dominant form of the tetramer in vivo. Moreover, loop-mediated repression in wild-type E. coli strains is facilitated by decreased DNA rigidity and high levels of flexibility in the LacR tetramer. In contrast, repression data for strains lacking HU gave a near-normal value of the DNA persistence length. These findings underscore the importance of both protein conformation and elasticity in the formation of small DNA loops widely observed in vivo, and demonstrate the utility of quantitatively analyzing gene regulation based on the mechanics of nucleoprotein complexes.