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.     

Probing DNA conformational changes with high temporal resolution by tethered particle motion

The tethered particle motion (TPM) technique informs about conformational changes of DNA molecules, e.g. upon looping or interaction with proteins, by tracking the Brownian motion of a particle probe tethered to a surface by a single DNA molecule and detecting changes of its amplitude of movement. We discuss in this context the time resolution of TPM, which strongly depends on the particle-DNA complex relaxation time, i.e. the characteristic time it takes to explore its configuration space by diffusion. By comparing theory, simulations and experiments, we propose a calibration of TPM at the dynamical level: we analyze how the relaxation time grows with both DNA contour length (from 401 to 2080 base pairs) and particle radius (from 20 to 150 nm). Notably we demonstrate that, for a particle of radius 20 nm or less, the hydrodynamic friction induced by the particle and the surface does not significantly slow down the DNA. This enables us to determine the optimal time resolution of TPM in distinct experimental contexts which can be as short as 20 ms.     

Bound Lac repressor protein differentially inhibits the unwinding reactions catalyzed by DNA helicases.

A partial duplex DNA substrate containing the Lac repressor binding site, within the duplex region, was constructed to examine the effect of bound Lac repressor on the unwinding reaction catalyzed by several DNA helicases. The substrate contained 90 base pairs of double-stranded DNA and, in the absence of Lac repressor, was effectively unwound by each of the seven helicases tested. The unwinding reactions catalyzed by Escherichia coli Rep protein, bacteriophage T4 Dda protein and E. coli DNA helicase I were not inhibited by the presence of bound Lac repressor. Both SV40 T antigen and E. coli helicase II were partially inhibited by bound repressor at the highest repressor concentrations tested. The helicase reactions catalyzed by E. coli DnaB protein and helicase IV were substantially inhibited by the presence of bound protein. When the length of the duplex region was increased to 323 base pairs the inhibition spectrum caused by bound Lac repressor on the unwinding reactions catalyzed by DnaB protein, helicase I and helicase II was essentially the same as that observed using the shorter partial duplex molecule. Inhibition of the unwinding reaction was due to the presence of bound Lac repressor as evidenced by the substantially weaker inhibition of helicase IV by Lac repressor in the presence of IPTG. In addition, we have shown that Rep protein displaces the bound repressor protein during the course of an unwinding reaction.     

Three-dimensional characterization of tethered microspheres by total internal reflection fluorescence microscopy

Tethered particle microscopy is a powerful tool to study the dynamics of DNA molecules and DNA-protein complexes in single-molecule experiments. We demonstrate that stroboscopic total internal reflection microscopy can be used to characterize the three-dimensional spatiotemporal motion of DNA-tethered particles. By calculating characteristic measures such as symmetry and time constants of the motion, well-formed tethers can be distinguished from defective ones for which the motion is dominated by aberrant surface effects. This improves the reliability of measurements on tether dynamics. For instance, in observations of protein-mediated DNA looping, loop formation is distinguished from adsorption and other nonspecific events.     

Mechanism for maturation-related reorganization of flavivirus glycoproteins

Flaviviruses, such as dengue, West Nile, and yellow fever viruses, assemble as fusion-incompetent particles and subsequently undergo a large reorganization of their glycoprotein envelope resulting in formation of mature infectious virions. Here we used a combination of three-dimensional cryo-electron tomography and two-dimensional image analysis to study pleomorphic maturation intermediates of dengue virus 2. Icosahedral symmetries of immature and mature regions within one particle were mismatched relative to each other. Furthermore, the orientation of the two regions relative to each other differed among particles. Therefore, there cannot be a specific pathway determining the maturation of all particles. Instead, the region with mature structure expands when glycoproteins on its boundary acquire suitable orientation and conformation to allow them to become a stable part of the mature region. This type of maturation is possible because the envelope glycoproteins are anchored to the phospholipid bilayer that is a part of flavivirus virions and are thus restricted to movement on the two-dimensional surface of the particle. Therefore, compounds that limit movement of the glycoproteins within the virus membrane might be used as inhibitors of flavivirus maturation.     

The solution structure of Lac repressor headpiece 62 complexed to a symmetrical lac operator

BACKGROUND: Lactose repressor protein (Lac) controls the expression of the lactose metabolic genes in Escherichia coli by binding to an operator sequence in the promoter of the lac operon. Binding of inducer molecules to the Lac core domain induces changes in tertiary structure that are propagated to the DNA-binding domain through the connecting hinge region, thereby reducing the affinity for the operator. Protein-protein and protein-DNA interactions involving the hinge region play a crucial role in the allosteric changes occurring upon induction, but have not, as yet, been analyzed in atomic detail. RESULTS: We have used nuclear magnetic resonance (NMR) spectroscopy and restrained molecular dynamics (rMD) to determine the structure of the Lac repressor DNA-binding domain (headpeice 62; HP62) in complex with a symmetrized lac operator. Analysis of the structures reveals specific interactions between Lac repressor and DNA that were not found in previously investigated Lac repressor-DNA complexes. Important differences with the previously reported structures of the HP56-DNA complex were found in the loop following the helix-turn-helix (HTH) motif. The protein-protein and protein-DNA interactions involving the hinge region and the deformations in the DNA structure could be delineated in atomic detail. The structures were also used for comparison with the available crystallographic data on the Lac and Pur repressor-DNA complexes. CONCLUSIONS: The structures of the HP62-DNA complex provide the basis for a better understanding of the specific recognition in the Lac repressor-operator complex. In addition, the structural features of the hinge region provide detailed insight into the protein-protein and protein-DNA interactions responsible for the high affinity of the repressor for operator DNA.     

One-Dimensional Brownian Motion of Charged Nanoparticles along Microtubules: A Model System for Weak Binding Interactions

Various proteins are known to exhibit one-dimensional Brownian motion along charged rodlike polymers, such as microtubules (MTs), actin, and DNA. The electrostatic interaction between the proteins and the rodlike polymers appears to be crucial for one-dimensional Brownian motion, although the underlying mechanism has not been fully clarified. We examined the interactions of positively-charged nanoparticles composed of polyacrylamide gels with MTs. These hydrophilic nanoparticles bound to MTs and displayed one-dimensional Brownian motion in a charge-dependent manner, which indicates that nonspecific electrostatic interaction is sufficient for one-dimensional Brownian motion. The diffusion coefficient decreased exponentially with an increasing particle charge (with the exponent being 0.10 k_BT per charge), whereas the duration of the interaction increased exponentially (exponent of 0.22 k_BT per charge). These results can be explained semiquantitatively if one assumes that a particle repeats a cycle of binding to and movement along an MT until it finally dissociates from the MT. During the movement, a particle is still electrostatically constrained in the potential valley surrounding the MT. This entire process can be described by a three-state model analogous to the Michaelis-Menten scheme, in which the two parameters of the equilibrium constant between binding and movement, and the rate of dissociation from the MT, are derived as a function of the particle charge density. This study highlights the possibility that the weak binding interactions between proteins and rodlike polymers, e.g., MTs, are mediated by a similar, nonspecific charge-dependent mechanism.     

Molecular and genetic characterization of lactose-metabolic genes of Streptococcus cremoris.

Lac+ plasmid DNA from Streptococcus cremoris H2 was subcloned with an Escherichia coli vector on a 3.5-kilobase-pair PstI-AvaI fragment. Genetic analysis of the cloned DNA was possible because linear Lac+ DNA fragments were productive in the S. sanguis transformation system. Complementation of S. sanguis Lac-mutants showed that the 3.5-kilobase-pair fragment included the structural gene for 6-phospho-beta-D-galactosidase and either enzyme II-lac or factor III-lac of the lactose-specific phosphoenolpyruvate-dependent phosphotransferase system. Expression of the S. cremoris-like 40,000-dalton 6-phospho-beta-D-galactosidase in S. sanguis Lac+ transformants, rather than the 52,000-dalton wild-type S. sanguis enzyme, demonstrated the occurrence of gene replacement and not gene repair. The evidence supports chromosomal integration as the mechanism by which S. sanguis Lac- recipients are converted to a Lac+ phenotype after transformation with Lac+ DNA. Southern blot data suggest that the Lac+ DNA does not reside on a transposon, but that integration always occurs within a specific HincII fragment of the recipient chromosome. Hybridization experiments demonstrate homology between the S. cremoris Lac+ DNA and cellular DNA from Lac+ strains of Streptococcus lactis, S. mutans, S. faecalis, and S. sanguis.     

Different characteristics distinguish early versus late arising adaptive mutations in Escherichia coli FC40

The Escherichia coli strain FC40 has frequently been employed to investigate the mechanism of adaptive mutations. The strain cannot utilize lactose due to a +1 frameshift mutation that reduces beta-galactosidase to about 1% of normal levels. Cells undergo a high rate of mutation from Lac- to Lac+ when cells are grown with lactose as the sole energy source. Almost all Lac+ colonies arising 3-6 days after plating result from a base pair deletion in runs of iterated base pairs within a 130-bp target region. In this study we characterized Lac+ colonies arising 3-10 days after plating. Temperature gradient gel electrophoresis (TGGE) was used to detect mutations in the target region as a function of the day a colony appears. TGGE results confirmed the occurrence of mutations within the target region in 36 of 37 FC40 Lac+ colonies arising on days 3-7. However, mutations in this region were not detected in 23 of 37 Lac+ colonies arising from days 8-10. Sequencing data verified the TGGE results. Half of the Lac+ mutants arising on days 8-10 with no base pair change in the target region were unstable and exhibited a Lac- phenotype after successive growth cycles in rich medium. The results suggest that amplification of the lac operon region is a common factor in late arising colonies, and that different characteristics distinguish early and late arising Lac+ colonies.     

Myosin II heavy chain null mutant of Dictyostelium exhibits defective intracellular particle movement

Both cellular motility and intracellular particle movement are compared between normal Dictyostelium amebae of strain AX4 and amebae of a myosin II heavy chain null mutant, HS2215, using the computer assisted "Dynamic Morphology System." In AX4 cells rapidly translocating in buffer, cytoplasmic expansion is apical and the majority of intracellular particles move anteriorly, towards the site of expansion. When these cells are pulsed with 10(-6) M cAMP, the peak concentration of the natural cAMP wave, cells stop translocating and average particle velocity decreases threefold within 2-4 s after cAMP addition. After 8 s, there is a partial rebound both in cytoplasmic expansion and particle velocity, but in both cases, original apical polarity is lost. In HS2215 cells in buffer, both cellular translocation and average particle velocity are already at the depressed levels observed in normal cells immediately after cAMP addition, and no anterior bias is observed in either the direction of cytoplasmic expansion or the direction of particle movement. The addition of cAMP to myosin-minus cells results in no additional effect. The results demonstrate that myosin II is necessary for (a) the rapid rate of intracellular particle movement, (b) the biased anterior directionality of particle movement, and (c) the rapid inhibition of particle movement by cAMP.     

Functional rules for lac repressor�Coperator associations and implications for protein�CDNA interactions

The Lac repressor has been used as a tool to understand protein–DNA recognition for many years. Recent experiments have demonstrated the ability of the Lac repressor to control gene expression in various eukaryotic systems, making the quest for an arsenal of protein–DNA binding partners desirable for potential therapeutic applications. Here, we present the results of the most exhaustive screen of Lac repressor‐DNA binding partners to date, resulting in the elucidation of functional rules for Lac–DNA binding. Even within the confines of a single protein–DNA scaffold, modes of binding of different protein–DNA partners are sufficiently diverse so as to prevent elucidation of generalized rules for recognition for a single protein, much less an entire protein family.     

Evoked potentials of the cat inferior colliculus to acoustic stimuli simulating sound source movement with different velocities in opposite directions

Amplitude changes of inferior colliculus evoked potentials (EPs) in anaesthetized adult cats were studied under presentation of acoustic stimuli simulating both azimuth-moving and stationary sound source. The movement was simulated with gradual changes of interaural time delay between binaurally presented click trains. It was shown that the amplitude of EPs elicited by "moving" signals depended on the velocity of movement. Amplitude differences between EPs to "moving" and stationary stimuli were observed under motion velocities up to 320 deg./s. The greatest response amplitudes in different experiments took place under velocities within the range of 67-320 deg./s with most of them recorded under velocities of 170 and 125 deg./s. Amplitude of the responses to lateral-medial movement with any velocity were always greater than those to opposite direction of movement with the same velocity.     

Lac repressor hinge flexibility and DNA looping: single molecule kinetics by tethered particle motion

The tethered particle motion (TPM) allows the direct detection of activity of a variety of biomolecules at the single molecule level. First pioneered for RNA polymerase, it has recently been applied also to other enzymes. In this work we employ TPM for a systematic investigation of the kinetics of DNA looping by wild-type Lac repressor (wt-LacI) and by hinge mutants Q60G and Q60 + 1. We implement a novel method for TPM data analysis to reliably measure the kinetics of loop formation and disruption and to quantify the effects of the protein hinge flexibility and of DNA loop strain on such kinetics. We demonstrate that the flexibility of the protein hinge has a profound effect on the lifetime of the looped state. Our measurements also show that the DNA bending energy plays a minor role on loop disruption kinetics, while a strong effect is seen on the kinetics of loop formation. These observations substantiate the growing number of theoretical studies aimed at characterizing the effects of DNA flexibility, tension and torsion on the kinetics of protein binding and dissociation, strengthening the idea that these mechanical factors in vivo may play an important role in the modulation of gene expression regulation.     

Inferring resource distributions from Atlantic bluefin tuna movements: an analysis based on net displacement and length of track

We use observed movement tracks of Atlantic bluefin tuna in the Gulf of Maine and mathematical modeling of this movement to identify possible resource patches. We infer bounds on the overall sizes and distribution of such patches, even though they are difficult to quantify by direct observation in situ. To do so, we segment individual fish tracks into intervals of distinct motion types based on the ratio of net displacement to length of track (ΔD/ΔL) over a time window Δt. To find the best segmentation, we optimize the fit of a random-walk movement model to each motion type. We compare results from two distinct movement models: biased turning and biased speed, to check the model-dependence of our inferences, and find that uncertainty in choice of movement model dominates the uncertainties of our conclusions. We find that our data are best described using two motion types: “localized” (ΔD/ΔL small) and “long-ranged” (ΔD/ΔL large). The biased turning model leads to significantly better resolution of localized movement intervals than the biased speed model. We hypothesize that localized movement corresponds to exploitation of resource patches. Comparison with visual behavior observations made during tracking suggests that many inferred intervals of localized motion do indeed correspond to feeding activity. From our analysis, we estimate that, on average, bluefin tuna in the Gulf of Maine encounter a resource patch every 2 h, that those patches have an average radius of 0.7-1.2 km, and that, overall, there are at most 5-9 such patches per in the region studied.     

DNA-binding properties of a lac repressor mutant incapable of forming tetramers

The interaction of proteins bound to sites widely separated on the genome is a recurrent motif in both prokaryotic and eukaryotic regulatory systems. Lac repressor mediates the formation of "DNA loops" by the simultaneous interaction of a single protein tetramer with two DNA-binding sites. The DNA-binding properties of a Lac repressor mutant (LacIadi) deficient in the association of protein dimers to tetramers was investigated. The results of quantitative footprint and gel mobility-shift titrations suggest that the wild-type Lac repressor (LacI+) binds cooperatively to two operator sites separated by 11 helical turns on a linear DNA restriction fragment by the formation of a "looped complex." LacIadi binds to this two-site operator non-cooperatively and without formation of a looped complex. These results demonstrate that the dimer-tetramer association of LacI+ is directly responsible for its cooperative binding and its ability to mediate formation of a looped complex. The Iadi mutation disrupts the monomer-dimer as well as eliminating the dimer-tetramer association equilibria while the DNA binding affinity of LacIadi to a single site is unchanged relative to the wild-type protein. These results suggest that DNA binding and dimer-tetramer association are functionally unlinked. The similarity of the DNA-binding properties of LacIadi and Gal repressor, a protein believed to function by mediating the formation of a looped complex, are discussed.     

A moving–resting process with an embedded Brownian motion for animal movements

Animal movements are of great importance in studying home ranges, migration routes, resource selection, and social interactions. The Global Positioning System provides relatively continuous animal tracking over time and long distances. Nevertheless, the continuous trajectory of an animal’s movement is usually only observed at discrete time points. Brownian bridge models have been used to model movement of an animal between two observed locations within a reasonably short time interval. Assuming that animals are in perpetual motion, these models ignore inactivity such as resting or sleeping. Using the latest developments in applied probability, we propose a moving–resting process model where an animal is assumed to alternate between a moving state, during which it moves in a Brownian motion, and a resting state, during which it does not move. Theoretical properties of the process are studied as a first step towards more realistic models for animal movements. Analytic expressions are derived for the distribution of one increment and two consecutive increments, and are validated with simulations. The induced bridge model conditioning on the starting and end points is used to compute an animal’s probability of occurrence in an observation area during the time of observation, which has wide applications in wildlife behavior research.     

A nonuniform distribution of excision repair synthesis in nucleosome core DNA

We have investigated the distribution in nucleosome core DNA of nucleotides incorporated by excision repair synthesis occurring immediately after UV irradiation in human cells. We show that the differences previously observed for whole nuclei between the DNase I digestion profiles of repaired DNA (following its refolding into a nucleosome structure) and bulk DNA are obtained for isolated nucleosome core particles. Analysis of the differences obtained indicates that they could reflect a significant difference in the level of repair-incorporated nucleotides at different sites within the core DNA region. To test this possibility directly, we have used exonuclease III digestion of very homogeneous sized core particle DNA to "map" the distribution of repair synthesis in these regions. Our results indicate that in a significant fraction of the nucleosomes the 5' and 3' ends of the core DNA are markedly enhanced in repair-incorporated nucleotides relative to the central region of the core particle. A best fit analysis indicates that a good approximation of the data is obtained for a distribution where the core DNA is uniformly labeled from the 5' end to position 62 and from position 114 to the 3' end, with the 52-base central region being devoid of repair-incorporated nucleotides. This distribution accounts for all of the quantitative differences observed previously between repaired DNA and bulk DNA following the rapid phase of nucleosome rearrangement when it is assumed that linker DNA and the core DNA ends are repaired with equal efficiency and the nucleosome structure of newly repaired DNA is identical with that of bulk chromatin. Furthermore, the 52-base central region that is devoid of repair synthesis contains the lowest frequency cutting sites for DNase I in vitro, as well as the only "internal" locations where two (rather than one) histones interact with a 10-base segment of each DNA strand.