The loopometer: a quantitative in vivo assay for DNA-looping proteins.

Proteins that can bring together separate DNA sites, either on the same or on different DNA molecules, are critical for a variety of DNA-based processes. However, there are no general and technically simple assays to detect proteins capable of DNA looping in vivo nor to quantitate their in vivo looping efficiency. Here, we develop a quantitative in vivo assay for DNA-looping proteins in Escherichia coli that requires only basic DNA cloning techniques and a LacZ assay. The assay is based on loop assistance, where two binding sites for the candidate looping protein are inserted internally to a pair of operators for the E. coli LacI repressor. DNA looping between the sites shortens the effective distance between the lac operators, increasing LacI looping and strengthening its repression of a lacZ reporter gene. Analysis based on a general model for loop assistance enables quantitation of the strength of looping conferred by the protein and its binding sites. We use this 'loopometer' assay to measure DNA looping for a variety of bacterial and phage proteins.

Effect of Photocaged Isopropyl ß-d-1-thiogalactopyranoside Solubility on the Light Responsiveness of LacI-controlled Expression Systems in Different Bacteria.

Photolabile protecting groups play a significant role in controlling biological functions and cellular processes in living cells and tissues, as light offers high spatiotemporal control, is non-invasive as well as easily tuneable. In the recent past, photo-responsive inducer molecules such as 6-nitropiperonyl-caged IPTG (NP-cIPTG) have been used as optochemical tools for Lac repressor-controlled microbial expression systems. To further expand the applicability of the versatile optochemical on-switch, we have investigated whether the modulation of cIPTG water solubility can improve the light responsiveness of appropriate expression systems in bacteria. To this end, we developed two new cIPTG derivatives with different hydrophobicity and demonstrated both an easy applicability for the light-mediated control of gene expression and a simple transferability of this optochemical toolbox to the biotechnologically relevant bacteria Pseudomonas putida and Bacillus subtilis. Notably, the more water-soluble cIPTG derivative proved to be particularly suitable for light-mediated gene expression in these alternative expression hosts.

Enhancement of LacI binding in vivo.

Transcription factors (TFs) bind to specific sequences in DNA to regulate transcription. Despite extensive measurements of TFs' dissociation constant (Kd) in vitro, their apparent Kdin vivo are usually unknown. LacI, a bacterial TF, is often used to artificially recruit proteins onto eukaryotic genomes. As LacI binds tightly to its recognition site (LacO) in vitro with a Kd about 10 picomolar (pM), it is often assumed that LacI also has high affinity to LacO in vivo. In this work, we measured LacI binding in living yeast cells using a fluorescent repressor operator system and found an apparent Kd of ∼0.6 µM, four orders of magnitude higher than that in vitro. By genetically altering (i) GFP-LacI structure, (ii) GFP-LacI stability, (iii) chromosome accessibility and (iv) LacO sequence, we reduced the apparent Kd to <10 nM. It turns out that the GFP tagging location and the fusion protein stability have a large effect on LacI binding, but surprisingly, chromosome accessibility only plays a mild role. These findings contribute to our quantitative understanding of the features that affect the apparent Kd of TF in cells. They also provide guidance for future design of more specific chromosomal recruitment through high-affinity TFs.

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.

Price of disorder in the lac repressor hinge helix.

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

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 (ka) and dissociation (kd) rate constants. RO-destabilizing solutes (urea, KCl) reduce ka comparably (urea) or more than (KCl) they increase kd, 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 kd without affecting ka 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.

DNA Internal Motion Likely Accelerates Protein Target Search in a Packed Nucleoid.

Transcription factors must diffuse through densely packed and coiled DNA to find their binding sites. Using a coarse-grained model of DNA and lac repressor (LacI) in the Escherichia coli nucleoid, simulations were performed to examine how LacI diffuses in such a space. Despite the canonical picture of LacI diffusing rather freely, in reality the DNA is densely packed, is not rigid but highly mobile, and the dynamics of DNA dictates to a great extent the LacI motion. A possibly better picture of unbound LacI motion is that of gated diffusion, where DNA confines LacI in a cage, but LacI can move between cages when hindering DNA strands move out of the way. Three-dimensional diffusion constants for unbound LacI computed from simulations closely match those for unbound LacI in vivo reported in the literature. The internal motions of DNA appear to be governed by strong internal forces arising from being crowded into the small space of the nucleoid. A consequence of the DNA internal motion is that protein target search may be accelerated.

De novo design of protein-protein interactions through modification of inter-molecular helix-helix interface residues.

For de novo design of protein-protein interactions (PPIs), information on the shape and chemical complementarity of their interfaces is generally required. Recent advances in computational PPI design have allowed for de novo design of protein complexes, and several successful examples have been reported. In addition, a simple and easy-to-use approach has also been reported that arranges leucines on a solvent-accessible region of an α-helix and places charged residues around the leucine patch to induce interactions between the two helical peptides. For this study, we adopted this approach to de novo design a new PPI between the helical bundle proteins sulerythrin and LARFH. A non-polar patch was created on an α-helix of LARFH around which arginine residues were introduced to retain its solubility. The strongest interaction found was for the LARFH variant cysLARFH-IV-3L3R and the sulerythrin mutant 6L6D (KD=0.16 µM). This artificial protein complex is maintained by hydrophobic and ionic interactions formed by the inter-molecular helical bundle structure. Therefore, by the simple and easy-to-use approach to create de novo interfaces on the α-helices, we successfully generated an artificial PPI. We also created a second LARFH variant with the non-polar patch surrounded by positively charged residues at each end. Upon mixing this LARFH variant with 6L6D, mesh-like fibrous nanostructures were observed by atomic force microscopy. Our method may, therefore, also be applicable to the de novo design of protein nanostructures.

Nucleic acid sensing with enzyme-DNA binding protein conjugates cascade and simple DNA nanostructures.

A versatile and universal DNA sensing platform is presented based on enzyme-DNA binding protein tags conjugates and simple DNA nanostructures. Two enzyme conjugates were thus prepared, with horseradish peroxidase linked to the dimeric single-chain bacteriophage Cro repressor protein (HRP-scCro) and glucose oxidase linked to the dimeric headpiece domain of Escherichia coli LacI repressor protein (GOx-dHP), and used in conjunction with a hybrid ssDNA-dsDNA detection probe. This probe served as a simple DNA nanostructure allowing first for target recognition through its target-complementary single-stranded DNA (ssDNA) part and then for signal generation after conjugate binding on the double-stranded DNA (dsDNA) containing the specific binding sites for the dHP and scCro DNA binding proteins. The DNA binding proteins chosen in this work have different sequence specificity, high affinity, and lack of cross-reactivity. The proposed sensing system was validated for the detection of model target ssDNA from high-risk human papillomavirus (HPV16) and the limits of detection of 45, 26, and 21 pM were achieved using the probes with scCro/dHP DNA binding sites ratio of 1:1, 2:1, and 1:2, respectively. The performance of the platform in terms of limit of detection was comparable to direct HRP systems using target-specific oligonucleotide-HRP conjugates. The ratio of the two enzymes can be easily manipulated by changing the number of binding sites on the detection probe, offering further optimization possibilities of the signal generation step. Moreover, since the signal is obtained in the absence of externally added hydrogen peroxide, the described platform is compatible with paper-based assays for molecular diagnostics applications. Finally, just by changing the ssDNA part of the detection probe, this versatile nucleic acid platform can be used for the detection of different ssDNA target sequences or in a multiplex detection configuration without the need to change any of the conjugates. Graphical abstract DNA sensing platform based on an immobilized ssDNA capture probe and a hybrid ssDNA-dsDNA detection probe that specifically hybridize with the ssDNA target. The hybrid ssDNA-dsDNA detection probe also provides the binding sites for the enzyme-DNA binding protein conjugates (HRP-scCro and GOx-dHP) that generate the colorimetric signal.

Flexibility and Disorder in Gene Regulation: LacI/GalR and Hox Proteins.

To modulate transcription, a variety of input signals must be sensed by genetic regulatory proteins. In these proteins, flexibility and disorder are emerging as common themes. Prokaryotic regulators generally have short, flexible segments, whereas eukaryotic regulators have extended regions that lack predicted secondary structure (intrinsic disorder). Two examples illustrate the impact of flexibility and disorder on gene regulation: the prokaryotic LacI/GalR family, with detailed information from studies on LacI, and the eukaryotic family of Hox proteins, with specific insights from investigations of Ultrabithorax (Ubx). The widespread importance of structural disorder in gene regulatory proteins may derive from the need for flexibility in signal response and, particularly in eukaryotes, in protein partner selection.

Digital switching in a biosensor circuit via programmable timing of gene availability.

Transient delivery of gene circuits is required in many potential applications of synthetic biology, yet the pre-steady-state processes that dominate this delivery route pose major challenges for robust circuit deployment. Here we show that site-specific recombinases can rectify undesired effects by programmable timing of gene availability in multigene circuits. We exemplify the concept with a proportional sensor for endogenous microRNA (miRNA) and show a marked reduction in its ground state leakage due to desynchronization of the circuit's repressor components and their repression target. The new sensors display a dynamic range of up to 1,000-fold compared to 20-fold in the standard configuration. We applied the approach to classify cell types on the basis of miRNA expression profile and measured >200-fold output differential between positively and negatively identified cells. We also showed major improvements in specificity with cytotoxic output. Our study opens new venues in gene circuit design via judicious temporal control of circuits' genetic makeup.

Optimized expression and purification of biophysical quantities of Lac repressor and Lac repressor regulatory domain.

The recombinant production of Lac repressor (LacI) in Escherichia coli is complicated by its ubiquitous use as a regulatory element in commercially-available expression vectors and host strains. While LacI-regulated expression systems are often used to produce recombinant LacI, the product can be heterogeneous and unsuitable for some studies. Alternative approaches include using unregulated vectors which typically suffer from low yield or vectors with promoters induced by metabolically active sugars which can dilute isotope labels necessary for certain biophysical studies. Here, an optimized expression system and isolation protocol for producing various constructs of LacI is introduced which eliminates these complications. The expression vector is an adaptation of the pASK backbone wherein expression of the lacI gene is regulated by an anhydrotetracyline inducible tetA promoter and the host strain lacks the lacI gene. Typical yields in highly deuterated minimal medium are nearly 2-fold greater than those previously reported. Notably, the new expression system is also able to produce the isolated regulatory domain of LacI without co-expression of the full-length protein and without any defects in cell viability, eliminating the inconvenient requirement for strict monitoring of cell densities during pre-culturing. Typical yields in highly deuterated minimal medium are significantly greater than those previously reported. Characterization by solution NMR shows that LacI constructs produced using this expression system are highly homogenous and functionally active.

Multiple LacI-mediated loops revealed by Bayesian statistics and tethered particle motion.

The bacterial transcription factor LacI loops DNA by binding to two separate locations on the DNA simultaneously. Despite being one of the best-studied model systems for transcriptional regulation, the number and conformations of loop structures accessible to LacI remain unclear, though the importance of multiple coexisting loops has been implicated in interactions between LacI and other cellular regulators of gene expression. To probe this issue, we have developed a new analysis method for tethered particle motion, a versatile and commonly used in vitro single-molecule technique. Our method, vbTPM, performs variational Bayesian inference in hidden Markov models. It learns the number of distinct states (i.e. DNA-protein conformations) directly from tethered particle motion data with better resolution than existing methods, while easily correcting for common experimental artifacts. Studying short (roughly 100 bp) LacI-mediated loops, we provide evidence for three distinct loop structures, more than previously reported in single-molecule studies. Moreover, our results confirm that changes in LacI conformation and DNA-binding topology both contribute to the repertoire of LacI-mediated loops formed in vitro, and provide qualitatively new input for models of looping and transcriptional regulation. We expect vbTPM to be broadly useful for probing complex protein-nucleic acid interactions.

Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes.

NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3-0.6 ppm and correlation coefficients (r(2)) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.

An information-based approach to change-point analysis with applications to biophysics and cell biology.

This article describes the application of a change-point algorithm to the analysis of stochastic signals in biological systems whose underlying state dynamics consist of transitions between discrete states. Applications of this analysis include molecular-motor stepping, fluorophore bleaching, electrophysiology, particle and cell tracking, detection of copy number variation by sequencing, tethered-particle motion, etc. We present a unified approach to the analysis of processes whose noise can be modeled by Gaussian, Wiener, or Ornstein-Uhlenbeck processes. To fit the model, we exploit explicit, closed-form algebraic expressions for maximum-likelihood estimators of model parameters and estimated information loss of the generalized noise model, which can be computed extremely efficiently. We implement change-point detection using the frequentist information criterion (which, to our knowledge, is a new information criterion). The frequentist information criterion specifies a single, information-based statistical test that is free from ad hoc parameters and requires no prior probability distribution. We demonstrate this information-based approach in the analysis of simulated and experimental tethered-particle-motion data.

Chemical Interactions of Polyethylene Glycols (PEGs) and Glycerol with Protein Functional Groups: Applications to Effects of PEG and Glycerol on Protein Processes.

In this work, we obtain the data needed to predict chemical interactions of polyethylene glycols (PEGs) and glycerol with proteins and related organic compounds and thereby interpret or predict chemical effects of PEGs on protein processes. To accomplish this, we determine interactions of glycerol and tetraEG with >30 model compounds displaying the major C, N, and O functional groups of proteins. Analysis of these data yields coefficients (α values) that quantify interactions of glycerol, tetraEG, and PEG end (-CH2OH) and interior (-CH2OCH2-) groups with these groups, relative to interactions with water. TetraEG (strongly) and glycerol (weakly) interact favorably with aromatic C, amide N, and cationic N, but unfavorably with amide O, carboxylate O, and salt ions. Strongly unfavorable O and salt anion interactions help make both small and large PEGs effective protein precipitants. Interactions of tetraEG and PEG interior groups with aliphatic C are quite favorable, while interactions of glycerol and PEG end groups with aliphatic C are not. Hence, tetraEG and PEG300 favor unfolding of the DNA-binding domain of lac repressor (lacDBD), while glycerol and di- and monoethylene glycol are stabilizers. Favorable interactions with aromatic and aliphatic C explain why PEG400 greatly increases the solubility of aromatic hydrocarbons and steroids. PEG400-steroid interactions are unusually favorable, presumably because of simultaneous interactions of multiple PEG interior groups with the fused ring system of the steroid. Using α values reported here, chemical contributions to PEG m-values can be predicted or interpreted in terms of changes in water-accessible surface area (ΔASA) and separated from excluded volume effects.

Amino acid positions subject to multiple coevolutionary constraints can be robustly identified by their eigenvector network centrality scores.

As proteins evolve, amino acid positions key to protein structure or function are subject to mutational constraints. These positions can be detected by analyzing sequence families for amino acid conservation or for coevolution between pairs of positions. Coevolutionary scores are usually rank-ordered and thresholded to reveal the top pairwise scores, but they also can be treated as weighted networks. Here, we used network analyses to bypass a major complication of coevolution studies: For a given sequence alignment, alternative algorithms usually identify different, top pairwise scores. We reconciled results from five commonly-used, mathematically divergent algorithms (ELSC, McBASC, OMES, SCA, and ZNMI), using the LacI/GalR and 1,6-bisphosphate aldolase protein families as models. Calculations used unthresholded coevolution scores from which column-specific properties such as sequence entropy and random noise were subtracted; "central" positions were identified by calculating various network centrality scores. When compared among algorithms, network centrality methods, particularly eigenvector centrality, showed markedly better agreement than comparisons of the top pairwise scores. Positions with large centrality scores occurred at key structural locations and/or were functionally sensitive to mutations. Further, the top central positions often differed from those with top pairwise coevolution scores: instead of a few strong scores, central positions often had multiple, moderate scores. We conclude that eigenvector centrality calculations reveal a robust evolutionary pattern of constraints-detectable by divergent algorithms--that occur at key protein locations. Finally, we discuss the fact that multiple patterns coexist in evolutionary data that, together, give rise to emergent protein functions.

Transcription closed and open complex dynamics studies reveal balance between genetic determinants and co-factors.

In E. coli, promoter closed and open complexes are key steps in transcription initiation, where magnesium-dependent RNA polymerase catalyzes RNA synthesis. However, the exact mechanism of initiation remains to be fully elucidated. Here, using single mRNA detection and dual reporter studies, we show that increased intracellular magnesium concentration affects Plac initiation complex formation resulting in a highly dynamic process over the cell growth phases. Mg2+ regulates transcription transition, which modulates bimodality of mRNA distribution in the exponential phase. We reveal that Mg2+ regulates the size and frequency of the mRNA burst by changing the open complex duration. Moreover, increasing magnesium concentration leads to higher intrinsic and extrinsic noise in the exponential phase. RNAP-Mg2+ interaction simulation reveals critical movements creating a shorter contact distance between aspartic acid residues and Nucleotide Triphosphate residues and increasing electrostatic charges in the active site. Our findings provide unique biophysical insights into the balanced mechanism of genetic determinants and magnesium ion in transcription initiation regulation during cell growth.

Effect of helix length on the stability of the Lac repressor antiparallel coiled coil.

The helix length dependence of the stability of antiparallel four-chain coiled coils is investigated using eight synthetic peptides (Lac21-Lac28) whose sequences are derived from the tetramerization domain of the Lac repressor protein. Previous studies using analytical ultracentrifugation sedimentation equilibrium experiments to characterize Lac21 and Lac28 justifies the use of a two state model to describe the unfolding behavior of these two peptides. Using circular dichroism spectropolarimetry as a measure of tetramer assembly, both chemical and thermal denaturation experiments were carried out to determine thermodynamic parameters. We found that the hydrophobic core residues provide the greatest impact on stability and, as a consequence, must reorganize the register of the antiparallel helices to accommodate the burial of the nonpolar amino acids. Addition of noncore residues appears to have only a minor effect on stability, and in some cases, show a slight destabilization.

Single molecule study of non-specific binding kinetics of LacI in mammalian cells.

Many key cellular processes are controlled by the association of DNA-binding proteins (DBPs) to specific sites. The kinetics of the search process leading to the binding of DBPs to their target locus are largely determined by transient interactions with non-cognate DNA. Using single-molecule microscopy, we studied the dynamics and non-specific binding to DNA of the Lac repressor (LacI) in the environment of mammalian nuclei. We measured the distribution of the LacI-DNA binding times at non-cognate sites and determined the mean residence time to be τ(1D) = 182 ms. This non-specific interaction time, measured in the context of an exogenous system such as that of human U2OS cells, is remarkably different compared to that reported for the LacI in its native environment in E. coli (<5 ms). Such a striking difference (more than 30 fold) suggests that the genome, its organization, and the nuclear environment of mammalian cells play important roles on the dynamics of DBPs and their non-specific DNA interactions. Furthermore, we found that the distribution of off-target binding times follows a power law, similar to what was reported for TetR in U2OS cells. We argue that a possible molecular origin of such a power law distribution of residence times is the large variability of non-cognate sequences found in the mammalian nucleus by the diffusing DBPs.