Thermodynamics of the interactions of lac repressor with variants of the symmetric lac operator: effects of converting a consensus site to a non-specific site

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

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

Structural and mutational analysis of MazE6-operator DNA complex provide insights into autoregulation of toxin-antitoxin systems

Of the 10 paralogs of MazEF Toxin-Antitoxin system in Mycobacterium tuberculosis, MazEF6 plays an important role in multidrug tolerance, virulence, stress adaptation and Non Replicative Persistant (NRP) state establishment. The solution structures of the DNA binding domain of MazE6 and of its complex with the cognate operator DNA show that transcriptional regulation occurs by binding of MazE6 to an 18 bp operator sequence bearing the TANNNT motif (-10 region). Kinetics and thermodynamics of association, as determined by NMR and ITC, indicate that the nMazE6-DNA complex is of high affinity. Residues in N-terminal region of MazE6 that are key for its homodimerization, DNA binding specificity, and the base pairs in the operator DNA essential for the protein-DNA interaction, have been identified. It provides a basis for design of chemotherapeutic agents that will act via disruption of TA autoregulation, leading to cell death.

The dimeric MazE6 antitoxin binds to a specific sequence in its cognate operator DNA for autoregulation, and the key residues for dimerization and DNA binding are identified.

Flexible Target Recognition of the Intrinsically Disordered DNA-Binding Domain of CytR Monitored by Single-Molecule Fluorescence Spectroscopy

The intrinsically disordered DNA-binding domain of cytidine repressor (CytR-DBD) folds in the presence of target DNA and regulates the expression of multiple genes in E. coli. To explore the conformational rearrangements in the unbound state and the target recognition mechanisms of CytR-DBD, we carried out single-molecule Förster resonance energy transfer (smFRET) measurements. The smFRET data of CytR-DBD in the absence of DNA show one major and one minor population assignable to an expanded unfolded state and a compact folded state, respectively. The population of the folded state increases and decreases upon titration with salt and denaturant, respectively, in an apparent two-state manner. The peak FRET efficiencies of both the unfolded and folded states change continuously with denaturant concentration, demonstrating the intrinsic flexibility of the DNA-binding domain and the deviation from a strict two-state transition. Remarkably, the CytR-DBD exhibits a compact structure when bound to both the specific and nonspecific DNA; however, the peak FRET efficiencies of the two structures are slightly but consistently different. The observed conformational heterogeneity highlights the potential structural changes required for CytR to bind variably spaced operator sequences.

Analysis of Protein-DNA Interactions Using Isothermal Titration Calorimetry: Successes and Failures

Isothermal titration calorimetry (ITC) is a golden standard for the characterization of protein-DNA binding affinities and allows direct assessment of the accompanying thermodynamic driving forces. Their interpretation can give insight into role of electrostatics, specificity of the DNA recognition, contribution of protein folding upon DNA binding and help to distinguish between minor and major groove binders. The main advantages of ITC are that the binding is measured in solution, and it requires no labeling of the samples, however, the method is not well suited for high-performance studies. Here we describe the sample preparation, a procedure to perform a typical ITC experiment, data analysis, and lastly discuss how to interpret the obtained thermodynamic parameters. In conclusion, we show examples of several unsuccessful ITC experiments and identify the underlying reasons for failed experiments. In most cases with a proper adjustment of the experimental setup, it was possible to obtain data appropriate for further analysis.

Biochemical and structural studies of target lectin SapL1 from the emerging opportunistic microfungus Scedosporium apiospermum

Scedosporium apiospermum is an emerging opportunistic fungal pathogen responsible for life-threatening infections in humans. Host-pathogen interactions often implicate lectins that have become therapeutic targets for the development of carbohydrate mimics for antiadhesive therapy. Here, we present the first report on the identification and characterization of a lectin from S. apiospermum named SapL1. SapL1 was found using bioinformatics as a homolog to the conidial surface lectin FleA from Aspergillus fumigatus known to play a role in the adhesion to host glycoconjugates present in human lung epithelium. In our strategy to obtain recombinant SapL1, we discovered the importance of osmolytes to achieve its expression in soluble form in bacteria. Analysis of glycan arrays indicates specificity for fucosylated oligosaccharides as expected. Submicromolar affinity was measured for fucose using isothermal titration calorimetry. We solved SapL1 crystal structure in complex with α-methyl-L-fucoside and analyzed its structural basis for fucose binding. We finally demonstrated that SapL1 binds to bronchial epithelial cells in a fucose-dependent manner. The information gathered here will contribute to the design and development of glycodrugs targeting SapL1.

Biophysical Characterization of Interaction between E. coli Alanyl-tRNA Synethase with its Promoter DNA

Background:

Aminoacyl-tRNA Synthetases (aaRSs) are well known for their role in the translation process. Lately investigators have discovered that this family of enzymes are also capable of executing a broad repertoire of functions that not only impact protein synthesis, but extend to a number of other activities. Till date, transcriptional regulation has so far only been described in E. coli Alanyl-tRNA synthetase and it was demonstrated that alaRS binds specifically to the palindromic DNA sequence flanking the gene’s transcriptional start site and thereby regulating its own transcription.

Objective:

In the present study, we have characterized some of the features of the alaRS-DNA binding using various biophysical techniques.

Methods:

To understand the role of full length protein and oligomerization of alaRS in promoter DNA binding, two mutants were constructed, namely, N700 (a monomer, containing the N-terminal aminoacylation domain but without the C-terminal part) and G674D (previously demonstrated to form full-length monomer). Protein-DNA binding study using fluorescence spectroscopy, analytical ultracentrifugation, Isothermal Titration Calorimetry was conducted.

Results:

Sedimentation equilibrium studies clearly demonstrated that monomeric variants were unable to bind promoter DNA. Isothermal Calorimetry (ITC) experiment was employed for further characterization of wild type alaRS-DNA interaction. It was observed that full length E. coli Alanyl-tRNA synthetase binds specifically with its promoter DNA and forms a dimer of dimers. On the other hand the two mutant variants were unable to bind with the DNA.

Conclusion:

In this study it was concluded that full length E. coli Alanyl-tRNA synthetase undergoes a conformational change in presence of its promoter DNA leading to formation of higher order structures. However, the exact mechanism behind this binding is currently unknown and beyond the scope of this study.

The lac repressor hinge helix in context: The effect of the DNA binding domain and symmetry

The Lac system of genes has been an important model system in understanding gene regulation. When the dimer lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The hinge region is disordered when binding to nonoperator sequences. This region of the protein must pay a conformational entropic penalty to order when it is bound to operator DNA. Structural studies show that this region is flexible. Previous simulations showed that this region is disordered when free in solution without the DNA binding domain present. Our simulations corroborate that this region is extremely flexible in solution, but we find that the presence of the DNA binding domain proximal to the hinge helix and salt make the ordered conformation more favorable even without DNA present.

Signatures of Specific DNA Binding by the AT-Rich Interaction Domain of BAF250a

The AT-rich interaction domain (ARID) containing BAF250a is a subunit of the BAF-A class of SWI/SNF chromatin remodeling complexes. The ARID belongs to a family of conserved DNA binding domains found in several eukaryotic proteins; however, its exact contribution to BAF250a function and the mechanism of its DNA binding are not well understood. Here we have probed the interaction of the BAF250a ARID with three different double-stranded DNA (dsDNA) sequences to understand its DNA binding properties. A comprehensive biophysical and thermodynamic study using nuclear magnetic resonance (NMR) spectroscopy and isothermal titration calorimetry revealed the complex nature of BAF250a ARID-DNA interactions. The thermodynamic signatures of the BAF250a ARID with 12 A-T bp dsDNA (AT-12) are distinct from those of 12 G-C bp dsDNA (GC-12) or 12 bp Dickerson dodecamer DNA (DD-12) sequences. We observed that the binding of the BAF250a ARID with AT-12 DNA is enthalpically driven in a tested temperature range of 5-25 °C. BAF250a ARID/AT-12 DNA interaction exhibited a larger negative calorimetric specific heat change (ΔC_p) compared to that of BAF250a ARID/GC-12 DNA or BAF250a ARID/DD-12 DNA interactions. In the presence of salt (NaCl), ARID/AT-12 DNA binding was less perturbed than ARID/GC-12 DNA or ARID/DD-12 DNA binding. Overall, these results show that BAF250a ARID/AT-12 DNA interaction has signatures of "specific" binding. Furthermore, using NMR chemical shift perturbation experiments, we have identified DNA binding residues on the BAF250a ARID and generated a data-driven HADDOCK model of the ARID/DNA complex that was further supported by mutating key lysine residues that were found to be important for DNA binding.

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.

Tethered Particle Motion: An Easy Technique for Probing DNA Topology and Interactions with Transcription Factors

Tethered Particle Motion (TPM) is a versatile in vitro technique for monitoring the conformations a linear macromolecule, such as DNA, can exhibit. The technique involves monitoring the diffusive motion of a particle anchored to a fixed point via the macromolecule of interest, which acts as a tether. In this chapter, we provide an overview of TPM, review the fundamental principles that determine the accuracy with which effective tether lengths can be used to distinguish different tether conformations, present software tools that assist in capturing and analyzing TPM data, and provide a protocol which uses TPM to characterize lac repressor-induced DNA looping. Critical to any TPM assay is the understanding of the timescale over which the diffusive motion of the particle must be observed to accurately distinguish tether conformations. Approximating the tether as a Hookean spring, we show how to estimate the diffusion timescale and discuss how it relates to the confidence with which tether conformations can be distinguished. Applying those estimates to a lac repressor titration assay, we describe how to perform a TPM experiment. We also provide graphically driven software which can be used to speed up data collection and analysis. Lastly, we detail how TPM data from the titration assay can be used to calculate relevant molecular descriptors such as the J factor for DNA looping and lac repressor-operator dissociation constants. While the included protocol is geared toward studying DNA looping, the technique, fundamental principles, and analytical methods are more general and can be adapted to a wide variety of molecular systems.

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

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

Thermodynamic Additivity for Impacts of Base-Pair Substitutions on Association of the Egr-1 Zinc-Finger Protein with DNA

The transcription factor Egr-1 specifically binds as a monomer to its 9 bp target DNA sequence, GCGTGGGCG, via three zinc fingers and plays important roles in the brain and cardiovascular systems. Using fluorescence-based competitive binding assays, we systematically analyzed the impacts of all possible single-nucleotide substitutions in the target DNA sequence and determined the change in binding free energy for each. Then, we measured the changes in binding free energy for sequences with multiple substitutions and compared them with the sum of the changes in binding free energy for each constituent single substitution. For the DNA variants with two or three nucleotide substitutions in the target sequence, we found excellent agreement between the measured and predicted changes in binding free energy. Interestingly, however, we found that this thermodynamic additivity broke down with a larger number of substitutions. For DNA sequences with four or more substitutions, the measured changes in binding free energy were significantly larger than predicted. On the basis of these results, we analyzed the occurrences of high-affinity sequences in the genome and found that the genome contains millions of such sequences that might functionally sequester Egr-1.

Aspects of protein-DNA interactions: a review of quantitative thermodynamic theory for modelling synthetic circuits utilising LacI and CI repressors, IPTG and the reporter gene lacZ

Protein-DNA interactions are central to the control of gene expression across all forms of life. The development of approaches to rigorously model such interactions has often been hindered both by a lack of quantitative binding data and by the difficulty in accounting for parameters relevant to the intracellular situation, such as DNA looping and thermodynamic non-ideality. Here, we review these considerations by developing a thermodynamically based mathematical model that attempts to simulate the functioning of an Escherichia coli expression system incorporating two of the best characterised prokaryotic DNA binding proteins, Lac repressor and lambda CI repressor. The key aim was to reproduce experimentally observed reporter gene activities arising from the expression of either wild-type CI repressor or one of three positive-control CI mutants. The model considers the role of several potentially important, but sometimes neglected, biochemical features, including DNA looping, macromolecular crowding and non-specific binding, and allowed us to obtain association constants for the binding of CI and its variants to a specific operator sequence.

A quantitative understanding of lac repressor's binding specificity and flexibility

Lac repressor, the first discovered transcriptional regulator, has been shown to confer multiple-modes of binding to its operator sites depending on the central spacer length. Other homolog members in the LacI/GalR family (PurR and YcjW) cannot bind their operator sites with similar structural flexibility. To decipher the underlying mechanism for this unique property, we used Spec-seq approach combined with site-directed mutagenesis to quantify the DNA binding specificity of multiple hybrids of lacI and PurR. We find that lac repressor's recognition di-residues YQ and its hinge helix loop regions are both critical for its structural flexibility. Also, specificity profiling of the whole lac operator suggests that a simple additive model from single variants suffice to predict other multivariant sites' energy reasonably well, and the genome occupancy model based on this specificity data correlates well with in vivo lac repressor binding profile.

Engineering a reagentless biosensor for single-stranded DNA to measure real-time helicase activity in Bacillus

Single-stranded DNA-binding protein (SSB) is a well characterized ubiquitous and essential bacterial protein involved in almost all aspects of DNA metabolism. Using the Bacillus subtilis SSB we have generated a reagentless SSB biosensor that can be used as a helicase probe in B. subtilis and closely related gram positive bacteria. We have demonstrated the utility of the probe in a DNA unwinding reaction using a helicase from Bacillus and for the first time, characterized the B. subtilis SSB's DNA binding mode switching and stoichiometry. The importance of SSB in DNA metabolism is not limited to simply binding and protecting ssDNA during DNA replication, as previously thought. It interacts with an array of partner proteins to coordinate many different aspects of DNA metabolism. In most cases its interactions with partner proteins is species-specific and for this reason, knowing how to produce and use cognate reagentless SSB biosensors in different bacteria is critical. Here we explain how to produce a B. subtilis SSB probe that exhibits 9-fold fluorescence increase upon binding to single stranded DNA and can be used in all related gram positive firmicutes which employ drastically different DNA replication and repair systems than the widely studied Escherichia coli. The materials to produce the B. subtilis SSB probe are commercially available, so the methodology described here is widely available unlike previously published methods for the E. coli SSB.

Protein-DNA binding in the absence of specific base-pair recognition

Until now, it has been reasonably assumed that specific base-pair recognition is the only mechanism controlling the specificity of transcription factor (TF)-DNA binding. Contrary to this assumption, here we show that nonspecific DNA sequences possessing certain repeat symmetries, when present outside of specific TF binding sites (TFBSs), statistically control TF-DNA binding preferences. We used high-throughput protein-DNA binding assays to measure the binding levels and free energies of binding for several human TFs to tens of thousands of short DNA sequences with varying repeat symmetries. Based on statistical mechanics modeling, we identify a new protein-DNA binding mechanism induced by DNA sequence symmetry in the absence of specific base-pair recognition, and experimentally demonstrate that this mechanism indeed governs protein-DNA binding preferences.

Tethered particle motion reveals that LacI·DNA loops coexist with a competitor-resistant but apparently unlooped conformation

The lac repressor protein (LacI) efficiently represses transcription of the lac operon in Escherichia coli by binding to two distant operator sites on the bacterial DNA and causing the intervening DNA to form a loop. We employed single-molecule tethered particle motion to observe LacI-mediated loop formation and breakdown in DNA constructs that incorporate optimized operator binding sites and intrinsic curvature favorable to loop formation. Previous bulk competition assays indirectly measured the loop lifetimes in these optimized DNA constructs as being on the order of days; however, we measured these same lifetimes to be on the order of minutes for both looped and unlooped states. In a range of single-molecule DNA competition experiments, we found that the resistance of the LacI-DNA complex to competitive binding is a function of both the operator strength and the interoperator sequence. To explain these findings, we present what we believe to be a new kinetic model of loop formation and DNA competition. In this proposed new model, we hypothesize a new unlooped state in which the unbound DNA-binding domain of the LacI protein interacts nonspecifically with nonoperator DNA adjacent to the operator site at which the second LacI DNA-binding domain is bound.

High-resolution specificity from DNA sequencing highlights alternative modes of Lac repressor binding

Knowing the specificity of transcription factors is critical to understanding regulatory networks in cells. The lac repressor-operator system has been studied for many years, but not with high-throughput methods capable of determining specificity comprehensively. Details of its binding interaction and its selection of an asymmetric binding site have been controversial. We employed a new method to accurately determine relative binding affinities to thousands of sequences simultaneously, requiring only sequencing of bound and unbound fractions. An analysis of 2560 different DNA sequence variants, including both base changes and variations in operator length, provides a detailed view of lac repressor sequence specificity. We find that the protein can bind with nearly equal affinities to operators of three different lengths, but the sequence preference changes depending on the length, demonstrating alternative modes of interaction between the protein and DNA. The wild-type operator has an odd length, causing the two monomers to bind in alternative modes, making the asymmetric operator the preferred binding site. We tested two other members of the LacI/GalR protein family and find that neither can bind with high affinity to sites with alternative lengths or shows evidence of alternative binding modes. A further comparison with known and predicted motifs suggests that the lac repressor may be unique in this ability and that this may contribute to its selection.

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.

Unusual characteristics of the DNA binding domain of epigenetic regulatory protein MeCP2 determine its binding specificity

The protein MeCP2 mediates epigenetic regulation by binding methyl-CpG (mCpG) sites on chromatin. MeCP2 consists of six domains of which one, the methyl binding domain (MBD), binds mCpG sites in duplex DNA. We show that solution conditions with physiological or greater salt concentrations or the presence of nonspecific competitor DNA is necessary for the MBD to discriminate mCpG from CpG with high specificity. The specificity for mCpG over CpG is >100-fold under these solution conditions. In contrast, the MBD does not discriminate hydroxymethyl-CpG from CpG. The MBD is unusual among site-specific DNA binding proteins in that (i) specificity is not conferred by the enhanced affinity for the specific site but rather by suppression of its affinity for generic DNA, (ii) its specific binding to mCpG is highly electrostatic, and (iii) it takes up as well as displaces monovalent cations upon DNA binding. The MBD displays an unusually high affinity for single-stranded DNA independent of modification or sequence. In addition, the MBD forms a discrete dimer on DNA via a noncooperative binding pathway. Because the affinity of the second monomer is 1 order of magnitude greater than that of nonspecific binding, the MBD dimer is a unique molecular complex. The significance of these results in the context of neuronal function and development and MeCP2-related developmental disorders such as Rett syndrome is discussed.

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.

Mutual enrichment in ranked lists and the statistical assessment of position weight matrix motifs

Background

Statistics in ranked lists is useful in analysing molecular biology measurement data, such as differential expression, resulting in ranked lists of genes, or ChIP-Seq, which yields ranked lists of genomic sequences. State of the art methods study fixed motifs in ranked lists of sequences. More flexible models such as position weight matrix (PWM) motifs are more challenging in this context, partially because it is not clear how to avoid the use of arbitrary thresholds.

Results

To assess the enrichment of a PWM motif in a ranked list we use a second ranking on the same set of elements induced by the PWM. Possible orders of one ranked list relative to another can be modelled as permutations. Due to sample space complexity, it is difficult to accurately characterize tail distributions in the group of permutations. In this paper we develop tight upper bounds on tail distributions of the size of the intersection of the top parts of two uniformly and independently drawn permutations. We further demonstrate advantages of this approach using our software implementation, mmHG-Finder, which is publicly available, to study PWM motifs in several datasets. In addition to validating known motifs, we found GC-rich strings to be enriched amongst the promoter sequences of long non-coding RNAs that are specifically expressed in thyroid and prostate tissue samples and observed a statistical association with tissue specific CpG hypo-methylation.

Conclusions

We develop tight bounds that can be calculated in polynomial time. We demonstrate utility of mutual enrichment in motif search and assess performance for synthetic and biological datasets. We suggest that thyroid and prostate-specific long non-coding RNAs are regulated by transcription factors that bind GC-rich sequences, such as EGR1, SP1 and E2F3. We further suggest that this regulation is associated with DNA hypo-methylation.

Thermodynamic and structural basis for relaxation of specificity in protein-DNA recognition

As a novel approach to the structural and functional properties that give rise to extremely stringent sequence specificity in protein-DNA interactions, we have exploited "promiscuous" mutants of EcoRI endonuclease to study the detailed mechanism by which changes in a protein can relax specificity. The A138T promiscuous mutant protein binds more tightly to the cognate GAATTC site than does wild-type EcoRI yet displays relaxed specificity deriving from tighter binding and faster cleavage at EcoRI* sites (one incorrect base pair). AAATTC EcoRI* sites are cleaved by A138T up to 170-fold faster than by wild-type enzyme if the site is abutted by a 5'-purine-pyrimidine (5'-RY) motif. When wild-type protein binds to an EcoRI* site, it forms structurally adapted complexes with thermodynamic parameters of binding that differ markedly from those of specific complexes. By contrast, we show that A138T complexes with 5'-RY-flanked AAATTC sites are virtually indistinguishable from wild-type-specific complexes with respect to the heat capacity change upon binding (∆C°P), the change in excluded macromolecular volume upon association, and contacts to the phosphate backbone. While the preference for the 5'-RY motif implicates contacts to flanking bases as important for relaxed specificity, local effects are not sufficient to explain the large differences in ∆C°P and excluded volume, as these parameters report on global features of the complex. Our findings therefore support the view that specificity does not derive from the additive effects of individual interactions but rather from a set of cooperative events that are uniquely associated with specific recognition.

Allostery without conformation change: modelling protein dynamics at multiple scales

The original ideas of Cooper and Dryden, that allosteric signalling can be induced between distant binding sites on proteins without any change in mean structural conformation, has proved to be a remarkably prescient insight into the rich structure of protein dynamics. It represents an alternative to the celebrated Monod-Wyman-Changeux mechanism and proposes that modulation of the amplitude of thermal fluctuations around a mean structure, rather than shifts in the structure itself, give rise to allostery in ligand binding. In a complementary approach to experiments on real proteins, here we take a theoretical route to identify the necessary structural components of this mechanism. By reviewing and extending an approach that moves from very coarse-grained to more detailed models, we show that, a fundamental requirement for a body supporting fluctuation-induced allostery is a strongly inhomogeneous elastic modulus. This requirement is reflected in many real proteins, where a good approximation of the elastic structure maps strongly coherent domains onto rigid blocks connected by more flexible interface regions.

Cooperativity in RNA-protein interactions: global analysis of RNA binding specificity

The control and function of RNA are governed by the specificity of RNA binding proteins. Here, we describe a method for global unbiased analysis of RNA-protein interactions that uses in vitro selection, high-throughput sequencing, and sequence-specificity landscapes. The method yields affinities for a vast array of RNAs in a single experiment, including both low- and high-affinity sites. It is reproducible and accurate. Using this approach,we analyzed members of the PUF (Pumilio and FBF) family of eukaryotic mRNA regulators. Our data identify effects of a specific protein partner on PUF-RNA interactions, reveal subsets of target sites not previously detected, and demonstrate that designer PUF proteins can precisely alter specificity. The approach described here is, in principle, broadly applicable for analysis of any molecule that binds RNA, including proteins, nucleic acids, and small molecules.

DNA looping in prokaryotes: experimental and theoretical approaches

Transcriptional regulation is at the heart of biological functions such as adaptation to a changing environment or to new carbon sources. One of the mechanisms which has been found to modulate transcription, either positively (activation) or negatively (repression), involves the formation of DNA loops. A DNA loop occurs when a protein or a complex of proteins simultaneously binds to two different sites on DNA with looping out of the intervening DNA. This simple mechanism is central to the regulation of several operons in the genome of the bacterium Escherichia coli, like the lac operon, one of the paradigms of genetic regulation. The aim of this review is to gather and discuss concepts and ideas from experimental biology and theoretical physics concerning DNA looping in genetic regulation. We first describe experimental techniques designed to show the formation of a DNA loop. We then present the benefits that can or could be derived from a mechanism involving DNA looping. Some of these are already experimentally proven, but others are theoretical predictions and merit experimental investigation. Then, we try to identify other genetic systems that could be regulated by a DNA looping mechanism in the genome of Escherichia coli. We found many operons that, according to our set of criteria, have a good chance to be regulated with a DNA loop. Finally, we discuss the proposition recently made by both biologists and physicists that this mechanism could also act at the genomic scale and play a crucial role in the spatial organization of genomes.

Sequence dependence of transcription factor-mediated DNA looping

DNA is subject to large deformations in a wide range of biological processes. Two key examples illustrate how such deformations influence the readout of the genetic information: the sequestering of eukaryotic genes by nucleosomes and DNA looping in transcriptional regulation in both prokaryotes and eukaryotes. These kinds of regulatory problems are now becoming amenable to systematic quantitative dissection with a powerful dialogue between theory and experiment. Here, we use a single-molecule experiment in conjunction with a statistical mechanical model to test quantitative predictions for the behavior of DNA looping at short length scales and to determine how DNA sequence affects looping at these lengths. We calculate and measure how such looping depends upon four key biological parameters: the strength of the transcription factor binding sites, the concentration of the transcription factor, and the length and sequence of the DNA loop. Our studies lead to the surprising insight that sequences that are thought to be especially favorable for nucleosome formation because of high flexibility lead to no systematically detectable effect of sequence on looping, and begin to provide a picture of the distinctions between the short length scale mechanics of nucleosome formation and looping.

Determinants of bacteriophage 933W repressor DNA binding specificity

We reported previously that 933W repressor apparently does not cooperatively bind to adjacent sites on DNA and that the relative affinities of 933W repressor for its operators differ significantly from that of any other lambdoid bacteriophage. These findings indicate that the operational details of the lysis-lysogeny switch of bacteriophage 933W are unique among lambdoid bacteriophages. Since the functioning of the lysis-lysogeny switch in 933W bacteriophage uniquely and solely depends on the order of preference of 933W repressor for its operators, we examined the details of how 933W repressor recognizes its DNA sites. To identify the specificity determinants, we first created a molecular model of the 933W repressor-DNA complex and tested the predicted protein-DNA interactions. These results of these studies provide a picture of how 933W repressor recognizes its DNA sites. We also show that, opposite of what is normally observed for lambdoid phages, 933W operator sequences have evolved in such a way that the presence of the most commonly found base sequences at particular operator positions serves to decrease, rather than increase, the affinity of the protein for the site. This finding cautions against assuming that a consensus sequence derived from sequence analysis defines the optimal, highest affinity DNA binding site for a protein.

Nonspecific DNA binding and bending by HUαβ: interfaces of the three binding modes characterized by salt-dependent thermodynamics

Previous isothermal titration calorimetry (ITC) and Förster resonance energy transfer studies demonstrated that Escherichia coli HU(αβ) binds nonspecifically to duplex DNA in three different binding modes: a tighter-binding 34-bp mode that interacts with DNA in large (>34 bp) gaps between bound proteins, reversibly bending it by 140(o) and thereby increasing its flexibility, and two weaker, modestly cooperative small site-size modes (10 bp and 6 bp) that are useful for filling gaps between bound proteins shorter than 34 bp. Here we use ITC to determine the thermodynamics of these binding modes as a function of salt concentration, and we deduce that DNA in the 34-bp mode is bent around-but not wrapped on-the body of HU, in contrast to specific binding of integration host factor. Analyses of binding isotherms (8-bp, 15-bp, and 34-bp DNA) and initial binding heats (34-bp, 38-bp, and 160-bp DNA) reveal that all three modes have similar log-log salt concentration derivatives of the binding constants (Sk(i)) even though their binding site sizes differ greatly; the most probable values of Sk(i) on 34-bp DNA or larger DNA are -7.5±0.5. From the similarity of Sk(i) values, we conclude that the binding interfaces of all three modes involve the same region of the arms and saddle of HU. All modes are entropy-driven, as expected for nonspecific binding driven by the polyelectrolyte effect. The bent DNA 34-bp mode is most endothermic, presumably because of the cost of HU-induced DNA bending, while the 6-bp mode is modestly exothermic at all salt concentrations examined. Structural models consistent with the observed Sk(i) values are proposed.

Fluorescence polarization anisotropy to measure RNA dynamics

RNA requires helical motion to fold and carry out its function. As RNA helical motion occurs on the nanosecond timescale, the timescale probed by fluorescence dyes, fluorescence polarization anisotropy (FPA) is a simple, yet powerful, technique to study helical dynamics in RNA. With the recent development of several fluorescent base analogs that have a nanosecond timescale lifetime in a duplex, FPA has begun to be used for characterizing RNA dynamics. Using the probe 6-methylisoxanthopterin (6-MI) as an example, we describe the procedure for carrying out FPA experiments on model oligonucleotide systems and in a complex RNA, the Tetrahymena group I intron. For smaller RNA systems, isolating the motion of the target helix from the overall tumbling of the whole RNA system is necessary, and nucleic acids binding proteins can be incorporated into the RNA system to increase the overall size of the system, slow the overall tumbling, and thereby reduce the anisotropy contribution from the overall tumbling to negligible. The procedure for incorporating one such protein, the Lac Repressor, is given as an example.

Coulombic free energy of polymeric nucleic acid: low- and high-salt analytical approximations for the cylindrical Poisson-Boltzmann model

An accurate analytical expression for the Coulombic free energy of DNA as a function of salt concentration ([salt]) is essential in applications to nucleic acid (NA) processes. The cylindrical model of DNA and the nonlinear Poisson-Boltzmann (NLPB) equation for ions in solution are among the simplest approaches capable of describing Coulombic interactions of NA and salt ions and of providing analytical expressions for thermodynamic quantities. Three approximations for Coulombic free energy G(u,infinity)(coul) of a polymeric nucleic acid are derived and compared with the numerical solution in a wide experimental range of 1:1 [salt] from 0.01 to 2 M. Two are obtained from the two asymptotic solutions of the cylindrical NLPB equation in the high-[salt] and low-[salt] limits: these are sufficient to determine G(u,infinity)(coul) of double-stranded (ds) DNA with 1% and of single-stranded (ss) DNA with 3% accuracy at any [salt]. The third approximation is experimentally motivated Taylor series up to the quadratic term in ln[salt] in the vicinity of the reference [salt] 0.15 M. This expression with three numerical coefficients (Coulombic free energy and its first and second derivatives at 0.15 M) predicts dependence of G(u,infinity)(coul) on [salt] within 2% of the numerical solution from 0.01 to 1 M for ss (a = 7 A, b = 3.4 A) and ds (a = 10 A, b = 1.7 A) DNA. Comparison of cylindrical free energy with that calculated for the all-atom structural model of linear B-DNA shows that the cylindrical model is completely sufficient above 0.01 M of 1:1 [salt]. The choice of two cylindrical parameters, the distance of closest approach of ion to cylinder axis (radius) a and the average axial charge separation b, is discussed in application to all-atom numerical calculations and analysis of experiment. Further development of analytical expression for Coulombic free energy with thermodynamic approaches accounting for ionic correlations and specific effects is suggested.

Probing the DNA-binding affinity and specificity of designed zinc finger proteins

Engineered transcription factors and endonucleases based on designed Cys(2)His(2) zinc finger domains have proven to be effective tools for the directed regulation and modification of genes. The introduction of this technology into both research and clinical settings necessitates the development of rapid and accurate means of evaluating both the binding affinity and binding specificity of designed zinc finger domains. Using a fluorescence anisotropy-based DNA-binding assay, we examined the DNA-binding properties of two engineered zinc finger proteins that differ by a single amino acid. We demonstrate that the protein with the highest affinity for a particular DNA site need not be the protein that binds that site with the highest degree of specificity. Moreover, by comparing the binding characteristics of the two proteins at varying salt concentrations, we show that the ionic strength makes significant and variable contributions to both affinity and specificity. These results have significant implications for zinc finger design as they highlight the importance of considering affinity, specificity, and environmental requirements in designing a DNA-binding domain for a particular application.

Specificity landscapes of DNA binding molecules elucidate biological function

Evaluating the specificity spectra of DNA binding molecules is a nontrivial challenge that hinders the ability to decipher gene regulatory networks or engineer molecules that act on genomes. Here we compare the DNA sequence specificities for different classes of proteins and engineered DNA binding molecules across the entire sequence space. These high-content data are visualized and interpreted using an interactive "specificity landscape" which simultaneously displays the affinity and specificity of a million-plus DNA sequences. Contrary to expectation, specificity landscapes reveal that synthetic DNA ligands match, and often surpass, the specificities of eukaryotic DNA binding proteins. The landscapes also identify differential specificity constraints imposed by diverse structural folds of natural and synthetic DNA binders. Importantly, the sequence context of a binding site significantly influences binding energetics, and utilizing the full contextual information permits greater accuracy in annotating regulatory elements within a given genome. Assigning such context-dependent binding values to every DNA sequence across the genome yields predictive genome-wide binding landscapes (genomescapes). A genomescape of a synthetic DNA binding molecule provided insight into its differential regulatory activity in cultured cells. The approach we describe will accelerate the creation of precision-tailored DNA therapeutics and uncover principles that govern sequence-specificity of DNA binding molecules.

Cohen's Conjecture, Howard's Hypothesis, and Ptashne's Ptruth: an exploration of the relationship between affinity and specificity

Both affinity and specificity for ligands directly influence the functions of biological macromolecules. Some investigators assume that there is a consistent relationship between the affinity of a receptor molecule for its cognate ligand(s) and the specificity of that same receptor (affinity for cognate versus non-cognate ligands). However, analysis of the range of physical factors that account for changes in affinity, in any particular direction and to any particular degree, of a receptor for a cognate ligand suggests strongly that such factors can have disparate effects on the affinities of the receptor for different non-cognate ligands. Therefore, there can be no simple relationship between affinity and specificity as defined by relative binding of the receptor to cognate and non-cognate ligands.

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.

Interactions of the osmolyte glycine betaine with molecular surfaces in water: thermodynamics, structural interpretation, and prediction of m-values

Noncovalent self-assembly of biopolymers is driven by molecular interactions between functional groups on complementary biopolymer surfaces, replacing interactions with water. Since individually these interactions are comparable in strength to interactions with water, they have been difficult to quantify. Solutes (osmolytes, denaturants) exert often large effects on these self-assembly interactions, determined in sign and magnitude by how well the solute competes with water to interact with the relevant biopolymer surfaces. Here, an osmometric method and a water-accessible surface area (ASA) analysis are developed to quantify and interpret the interactions of the remarkable osmolyte glycine betaine (GB) with molecular surfaces in water. We find that GB, lacking hydrogen bond donors, is unable to compete with water to interact with anionic and amide oxygens; this explains its effectiveness as an osmolyte in the Escherichia coli cytoplasm. GB competes effectively with water to interact with amide and cationic nitrogens (hydrogen bonding) and especially with aromatic hydrocarbon (cation-pi). The large stabilizing effect of GB on lac repressor-lac operator binding is predicted quantitatively from ASA information and shown to result largely from dehydration of anionic DNA phosphate oxygens in the protein-DNA interface. The incorporation of these results into theoretical and computational analyses will likely improve the ability to accurately model intra- and interprotein interactions. Additionally, these results pave the way for development of solutes as kinetic/mechanistic and thermodynamic probes of conformational changes and formation/disruption of molecular interfaces that occur in the steps of biomolecular self-assembly processes.

Concentration and length dependence of DNA looping in transcriptional regulation

In many cases, transcriptional regulation involves the binding of transcription factors at sites on the DNA that are not immediately adjacent to the promoter of interest. This action at a distance is often mediated by the formation of DNA loops: Binding at two or more sites on the DNA results in the formation of a loop, which can bring the transcription factor into the immediate neighborhood of the relevant promoter. These processes are important in settings ranging from the historic bacterial examples (bacterial metabolism and the lytic-lysogeny decision in bacteriophage), to the modern concept of gene regulation to regulatory processes central to pattern formation during development of multicellular organisms. Though there have been a variety of insights into the combinatorial aspects of transcriptional control, the mechanism of DNA looping as an agent of combinatorial control in both prokaryotes and eukaryotes remains unclear. We use single-molecule techniques to dissect DNA looping in the lac operon. In particular, we measure the propensity for DNA looping by the Lac repressor as a function of the concentration of repressor protein and as a function of the distance between repressor binding sites. As with earlier single-molecule studies, we find (at least) two distinct looped states and demonstrate that the presence of these two states depends both upon the concentration of repressor protein and the distance between the two repressor binding sites. We find that loops form even at interoperator spacings considerably shorter than the DNA persistence length, without the intervention of any other proteins to prebend the DNA. The concentration measurements also permit us to use a simple statistical mechanical model of DNA loop formation to determine the free energy of DNA looping, or equivalently, the for looping.

Specificity and affinity of Lac repressor for the auxiliary operators O2 and O3 are explained by the structures of their protein-DNA complexes

The structures of a dimeric mutant of the Lac repressor DNA-binding domain complexed with the auxiliary operators O2 and O3 have been determined using NMR spectroscopy and compared to the structures of the previously determined Lac-O1 and Lac-nonoperator complexes. Structural analysis of the Lac-O1 and Lac-O2 complexes shows highly similar structures with very similar numbers of specific and nonspecific contacts, in agreement with similar affinities for these two operators. The left monomer of the Lac repressor in the Lac-O3 complex retains most of these specific contacts. However, in the right half-site of the O3 operator, there is a significant loss of protein-DNA contacts, explaining the low affinity of the Lac repressor for the O3 operator. The binding mode in the right half-site resembles that of the nonspecific complex. In contrast to the Lac-nonoperator DNA complex where no hinge helices are formed, the stability of the hinge helices in the weak Lac-O3 complex is the same as in the Lac-O1 and Lac-O2 complexes, as judged from the results of hydrogen/deuterium experiments.

Prevalence of temperature-dependent heat capacity changes in protein-DNA interactions

A large, negative DeltaCp of DNA binding is a thermodynamic property of the majority of sequence-specific DNA-protein interactions, and a common, but not universal property of non-sequence-specific DNA binding. In a recent study of the binding of Taq polymerase to DNA, we showed that both the full-length polymerase and its "Klentaq" large fragment bind to primed-template DNA with significant negative heat capacities. Herein, we have extended this analysis by analyzing this data for temperature-variable heat capacity effects (DeltaDeltaCp), and have similarly analyzed an additional 47 protein-DNA binding pairs from the scientific literature. Over half of the systems examined can be easily fit to a function that includes a DeltaDeltaCp parameter. Of these, 90% display negative DeltaDeltaCp values, with the result that the DeltaCp of DNA binding will become more negative with rising temperature. The results of this collective analysis have potentially significant consequences for current quantitative theories relating DeltaCp values to changes in accessible surface area, which rely on the assumption of temperature invariance of the DeltaCp of binding. Solution structural data for Klentaq polymerase demonstrate that the observed heat capacity effects are not the result of a coupled folding event.

Thermodynamic analysis of acetylation-dependent Pb1 bromodomain-histone H3 interactions

An acetyl-histone peptide library was used to determine the thermodynamic parameters that define acetylation-dependent bromodomain-histone interactions. Bromodomains interact with histones by binding acetylated lysines. The bromodomain used in this study, BrD3, is derived from the polybromo-1 protein, which is a subunit of the PBAF chromatin remodeling complex. Steady-state fluorescence anisotropy was used to examine the variations in specificity and affinity that drive molecular recognition. Temperature and salt concentration dependence studies demonstrate that the hydrophobic effect is the primary driving force, consistent with lysine acetylation being required for binding. An electrostatic effect was observed in only two complexes where the acetyl-lysine was adjacent to an arginine. The large change in heat capacity determined for the specific complex suggests that the dehydrated BrD3-histone interface forms a tightly bound, high-affinity complex with the target site. These explorations into the thermodynamic driving forces that confer acetylation site-dependent BrD3-histone interactions improve our understanding of how individual bromodomains work in isolation. Furthermore, this work will permit the development of hypotheses regarding how the native Pb1, and the broader class of bromodomain proteins, directs multisubunit chromatin remodeling complexes to specific acetyl-nucleosome sites in vivo.

Effect of somatic mutation on DNA binding properties of anti-DNA autoantibodies

Autoantibodies that bind DNA are a hallmark of systemic lupus erythematosus. A subset of autoantibody*DNA complexes localize to kidney tissue and lead to damage and even death. 11F8, 9F11, and 15B10 are clonally related anti-DNA autoantibodies isolated from an autoimmune mouse. 11F8 binds ssDNA in a sequence-specific manner and causes tissue damage, while 9F11 and 15B10 bind ssDNA non-specifically and are benign. Among these antibodies, DNA binding properties are mediated by five amino acid differences in primary sequence. Thermodynamic and kinetic parameters associated with recognition of structurally different DNA sequences were determined for each antibody to provide insight toward recognition strategies, and to explore a link between binding properties and disease pathogenesis. A model of 11F8 bound to its high affinity consensus sequence provides a foundation for understanding the differences in thermodynamic and kinetic parameters between the three mAbs. Our data suggest that 11F8 utilizes the proposed ssDNA recognition motif including (Y32)V(L), a hydrogen bonding residue at (91)V(L), and an aromatic residue at the tip of the third heavy chain complementarity determining region. Interestingly, a somatic mutation to arginine at (31)V(H) in 11F8 may afford additional binding site contacts including (R31)V(H), (R96)V(H), and (R98)V(H) that could determine specificity.

Conformational proofreading: the impact of conformational changes on the specificity of molecular recognition

To perform recognition, molecules must locate and specifically bind their targets within a noisy biochemical environment with many look-alikes. Molecular recognition processes, especially the induced-fit mechanism, are known to involve conformational changes. This raises a basic question: Does molecular recognition gain any advantage by such conformational changes? By introducing a simple statistical-mechanics approach, we study the effect of conformation and flexibility on the quality of recognition processes. Our model relates specificity to the conformation of the participant molecules and thus suggests a possible answer: Optimal specificity is achieved when the ligand is slightly off target; that is, a conformational mismatch between the ligand and its main target improves the selectivity of the process. This indicates that deformations upon binding serve as a conformational proofreading mechanism, which may be selected for via evolution.

A single, specific thymine mutation in the ComK-binding site severely decreases binding and transcription activation by the competence transcription factor ComK of Bacillus subtilis

The competence transcription factor ComK plays a central role in competence development in Bacillus subtilis by activating the transcription of the K regulon. ComK-activated genes are characterized by the presence of a specific sequence to which ComK binds, a K-box, in their upstream DNA region. Each K-box consists of two AT-boxes with the consensus sequence AAAA-(N)(5)-TTTT, which are separated by a flexible spacer resulting in either two, three, or four helical turns between the starting nucleotides of the repeating AT-box units. In this study, the effects of potential determinants of ComK regulation in K-boxes were investigated by testing ComK's transcription activation and DNA-binding affinity on altered K-boxes with mutations either in the spacer between the AT-boxes or in the consensus sequence of the AT-boxes. The most striking result demonstrates the importance of the second thymine base in the AT-boxes. Mutation of this T into a guanine resulted in a threefold reduction in transcription activation and DNA binding by ComK. Transcription activation, as well as DNA binding, was almost completely abolished when both AT-boxes contained a T(2)-to-G mutation. This result was corroborated by in silico analyses demonstrating that a combination of mutations at the T(2) positions of both AT-boxes is not found among any ComK-activated K-boxes, indicating that at least one consensus T(2) position is required to maintain a functional K-box. The results suggest an important structural role for T(2) in ComK binding, probably by its specific position in the minor groove of the DNA.

Optimization of a stochastically simulated gene network model via simulated annealing

By rearranging naturally occurring genetic components, gene networks can be created that display novel functions. When designing these networks, the kinetic parameters describing DNA/protein binding are of great importance, as these parameters strongly influence the behavior of the resulting gene network. This article presents an optimization method based on simulated annealing to locate combinations of kinetic parameters that produce a desired behavior in a genetic network. Since gene expression is an inherently stochastic process, the simulation component of simulated annealing optimization is conducted using an accurate multiscale simulation algorithm to calculate an ensemble of network trajectories at each iteration of the simulated annealing algorithm. Using the three-gene repressilator of Elowitz and Leibler as an example, we show that gene network optimizations can be conducted using a mechanistically realistic model integrated stochastically. The repressilator is optimized to give oscillations of an arbitrary specified period. These optimized designs may then provide a starting-point for the selection of genetic components needed to realize an in vivo system.