Real-time imaging of single DNA molecules with fluorescence microscopy.

A fluorescence microscopy technique was used to image the dynamics of individual DNA molecules. Lambda, calf thymus, cosmid (circular), and T4 DNA were studied with the fluorescent dye acridine orange. Experiments with DNAase I were conducted, and the results indicate that these observations correspond to DNA molecules. The results of experiments with circular DNA provide strong evidence that these were single DNA molecules. Molecules were observed free in solution or attached to a glass or copper surface at one or several points. The Brownian motion of these molecules was observed, indicating that DNA in solution exists in a partially supercoiled state. Some molecules appeared stretched and were attached to the surface by their termini; the lengths of these molecules were measured. Such molecules also exhibited elastic behavior upon breaking. The power of this technique is demonstrated in images of cosmid DNA molecules, catenanes, and DNA extending from T4 phage particles. These results suggest immediate applications to molecular biology, such as examining the dynamics of protein-DNA interactions. Areas of ongoing research are discussed.     

Sliding of Proteins Non-specifically Bound to DNA: Brownian Dynamics Studies with Coarse-Grained Protein and DNA Models

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.     

Brownian Dynamics Simulation of DNA Gel Electrophoresis

Abstract

Brownian dynamics computer simulation technique was applied to investigate DNA dynamics in gel electrophoresis. Under a constant electric field of moderate strength, large DNA chains take stretched and contracted conformations alternatively during the migration. The conformation change is quasi-periodic under certain conditions, and its frequency is closely related to the experimentally-found suitable frequency of pulse field gel electrophoresis.     

Stochastic Simulation of Signal Transduction: Impact of the Cellular Architecture on Diffusion

The transduction of signals depends on the translocation of signaling molecules to specific targets. Undirected diffusion processes play a key role in the bridging of spaces between different cellular compartments. The diffusion of the molecules is, in turn, governed by the intracellular architecture. Molecular crowding and the cytoskeleton decrease macroscopic diffusion. This article shows the use of a stochastic simulation method to study the effects of the cytoskeleton structure on the mobility of macromolecules. Brownian dynamics and single particle tracking were used to simulate the diffusion process of individual molecules through a model cytoskeleton. The resulting average effective diffusion is in line with data obtained in the in vitro and in vivo experiments. It shows that the cytoskeleton structure strongly influences the diffusion of macromolecules. The simulation method used also allows the inclusion of reactions in order to model complete signaling pathways in their spatio-temporal dynamics, taking into account the effects of the cellular architecture.     

Simulation of the conformation and dynamics of a double-helical model for DNA

We propose a partially flexible, double-helical model for describing the conformational and dynamic properties of DNA. In this model, each nucleotide is represented by one element (bead), and the known geometrical features of the double helix are incorporated in the equilibrium conformation. Each bead is connected to a few neighbor beads in both strands by means of stiff springs that maintain the connectivity but still allow for some extent of flexibility and internal motion. We have used Brownian dynamics simulation to sample the conformational space and monitor the overall and internal dynamics of short DNA pieces, with up to 20 basepairs. From Brownian trajectories, we calculate the dimensions of the helix and estimate its persistence length. We obtain translational diffusion coefficient and various rotational relaxation times, including both overall rotation and internal motion. Although we have not carried out a detailed parameterization of the model, the calculated properties agree rather well with experimental data available for those oligomers.     

Calculation of self-diffusion coefficients of the [BMIM][TFSA]/water system by molecular dynamics simulation

We performed a molecular dynamics simulation to calculate the self-diffusion coefficients of 1-Butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide and water in a water–ionic liquid mixture. We then compared the simulated self-diffusion coefficients of cation, anion and water molecules with experimental data and with simulated data from the literature. Although the simulation overestimated the self-diffusion coefficients of ions, the simulated results qualitatively reproduced the enhancement of the self-diffusion coefficients of water as the water molar fraction increased. We also calculated the radial distribution functions to investigate the solution structure, i.e. the clustering of water molecules. The clustering of water in ionic liquid was found to play an important role in the enhancement of the diffusion of water molecules in the ionic liquid.     

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.     

Determining Protein-Induced DNA Bending in Force-Extension Experiments: Theoretical Analysis

Computer simulations were used to investigate the possibility of determining protein-induced DNA bend angles by measuring the extension of a single DNA molecule. Analysis of the equilibrium sets of DNA conformations showed that shortening of DNA extension by a single protein-induced DNA bend can be as large as 35 nm. The shortening has a maximum value at the extending force of ∼0.1 pN. At this force, the DNA extension experiences very large fluctuations that dramatically complicate the measurement. Using Brownian dynamics simulation of a DNA molecule extended by force, we were able to estimate the observation time needed to obtain the desired accuracy of the extension measurement. Also, the simulation revealed large fluctuations of the force, acting on the attached magnetic bead from the stretched DNA molecule.     

Brownian dynamics simulation of rigid particles of arbitrary shape in external fields

We have developed a Brownian dynamics simulation algorithm to generate Brownian trajectories of an isolated, rigid particle of arbitrary shape in the presence of electric fields or any other external agents. Starting from the generalized diffusion tensor, which can be calculated with the existing HYDRO software, the new program BROWNRIG (including a case-specific subprogram for the external agent) carries out a simulation that is analyzed later to extract the observable dynamic properties. We provide a variety of examples of utilization of this method, which serve as tests of its performance, and also illustrate its applicability. Examples include free diffusion, transport in an electric field, and diffusion in a restricting environment.     

Hydrodynamics of segmentally flexible macromolecules

Segmentally flexible macromolecules are composed of a few rigid subunits linked by joints which are more or less flexible. The dynamics in solution of this type of macromolecule present special aspects that are reviewed here. Three alternative approaches are described. One is the rigid-body treatment, which is shown to be valid for overall dynamic properties such as translational diffusion and intrinsic viscosity. Another approach is the Harvey-Wegener treatment, which is particularly suited for rotational diffusion. The simplest version of this treatment, which ignores hydrodynamic interaction (HI) effects, is found to be quite accurate when compared to a more rigorous version including HI. A third approach is the Brownian dynamics simulation that, albeit at some computational cost, might describe rigorously cases of arbitrary complexity. This technique has been used to test the approximations in the rigid-body and Harvey-Wegener treatments, thus allowing a better understanding of their validity. Brownian trajectories of simplified models such as the trumbbell and the broken rod have been simulated. The comparison of the decay rates of some correlation functions with the predictions of the two treatments leads to a general conclusion: the Harvey-Wegener treatment determines the initial rate, while the long-time behavior is dominated by the rigid-body relaxation time. As an example of application to a specific biological macromolecule, we present a simulation of an immunoglobulin molecule, showing how Brownian Dynamics can be used to predict rotational and internal dynamics. Another typical example is myosin. Literature data of hydrodynamic properties of whole myosin and the myosin rod are compared with predictions from the Harvey-Wegener and rigid-body treatments. The present situation of the problem on myosin flexibility is analyzed, and some indications are given for future experimental and simulation work.     

Novel display of knotted DNA molecules by two-dimensional gel electrophoresis

We describe a two-dimensional agarose gel electrophoresis procedure that improves the resolution of knotted DNA molecules. The first gel dimension is run at low voltage, and DNA knots migrate according to their compactness. The second gel dimension is run at high voltage, and DNA knots migrate according to other physical parameters such as shape and flexibility. In comparison with one-dimensional gel electrophoresis, this procedure segregates the knotted DNA molecules from other unknotted forms of DNA, and partially resolves populations of knots that have the same number of crossings. The two-dimensional display may allow quantitative and qualitative characterization of different types of DNA knots simply by gel velocity.     

Brownian dynamics simulation of probe diffusion in DNA: effects of probe size, charge and DNA concentration

We have used Brownian dynamics simulation to study probe diffusion in solutions of short chain DNA using our previously developed simulation algorithm. We have examined the effect of probe size, charge, and DNA concentration on the probe diffusion coefficient, with the aim of gaining insight into the diffusion of proteins in a concentrated DNA environment. In these simulations, DNA was modeled as a worm-like chain of hydrodynamically equivalent spherical frictional elements while probe particles were modeled as spheres of given charge and hydrodynamic radius. The simulations allowed for both short range Lennard-Jones interactions and long ranged electrostatic interactions between charged particles. For uncharged systems, we find that the effects of probe size and DNA concentration on the probe diffusion coefficient are consistent with excluded volume models and we interpret our results in terms of both empirical scaling laws and the predictions of scaled particle theory. For charged systems, we observe that the effects of probe size and charge are most pronounced for the smallest probes and interpret the results in terms of the probe charge density. For an ionic strength of 0.1 M we find that, below a critical probe surface charge density, the probe diffusion coefficient is largely independent of probe charge and only weakly dependent on the DNA charge. These effects are discussed in terms of the interactions between the probe and the DNA matrix and are interpreted in terms of both the underlying physics of transport in concentrated solutions and the assumptions of the simulation model.     

Topo IV is the topoisomerase that knots and unknots sister duplexes during DNA replication

DNA topology plays a crucial role in all living cells. In prokaryotes, negative supercoiling is required to initiate replication and either negative or positive supercoiling assists decatenation. The role of DNA knots, however, remains a mystery. Knots are very harmful for cells if not removed efficiently, but DNA molecules become knotted in vivo. If knots are deleterious, why then does DNA become knotted? Here, we used classical genetics, high-resolution 2D agarose gel electrophoresis and atomic force microscopy to show that topoisomerase IV (Topo IV), one of the two type-II DNA topoisomerases in bacteria, is responsible for the knotting and unknotting of sister duplexes during DNA replication. We propose that when progression of the replication forks is impaired, sister duplexes become loosely intertwined. Under these conditions, Topo IV inadvertently makes the strand passages that lead to the formation of knots and removes them later on to allow their correct segregation.     

Brownian dynamics simulations of probe and self-diffusion in concentrated protein and DNA solutions.

We have developed a Brownian dynamics algorithm for simulating probe and self-diffusion in concentrated solutions of DNA and protein. In these simulations, proteins are represented as spheres with radii given by their hydrodynamic radii, while DNA is modeled as a wormlike chain of hydrodynamically equivalent spherical frictional elements. The molecular interaction potentials employed by the program allow for intramolecular stretching and bending motions of the DNA chains, short-range Lennard-Jones interactions, and long-range electrostatic interactions. To test the program, we have carried out simulations of bovine serum albumin (BSA) probe diffusion and DNA self-diffusion in solutions of short-chain DNA as a function of both DNA concentration and solution ionic strength. In addition, we report on simulations of BSA self-diffusion as a function of BSA concentration and ionic strength. Based on a comparison to available experimental data, we find that our simulations accurately predict these transport properties under conditions of physiological salt concentration and predict the stronger concentration dependence observed at lower salt concentrations. These results are discussed in light of the nature of the intermolecular interactions in such systems and the approximations and limitations of the simulation algorithm.     

Multiple Time Step Brownian Dynamics for Long Time Simulation of Biomolecules

We report a multiple time step algorithm applied to an atomistic Brownian dynamics simulation for simulating the long time scale dynamics of biomolecules. The algorithm was based on the original multiple time step method; a short time step was used to keep faster motions in local equilibrium. When applied to a 28-mer # # ! folded peptide, the simulation gave stable trajectories and the computation time was reduced by a factor of 160 compared to a conventional molecular dynamics simulation using explicit water molecules. We applied it for the folding simulation of a 13-mer ! -helical peptide, giving a successful folding simulation. These results indicate that the Brownian dynamics with the multiple time step algorithm is useful for studies of biomolecular motions by long time simulation.     

Brownian dynamics simulation of nucleosome formation and disruption under stretching

Using a Brownian dynamics simulation, we numerically studied the interaction of DNA with histone and proposed an octamer-rotation model to describe the process of nucleosome formation. Nucleosome disruption under stretching was also simulated. The theoretical curves of extension versus time as well as of force versus extension are consistent with previous experimental results.     

Kinetics of nanochain formation in a simplified model of amelogenin biomacromolecules

We show that the kinetics of nanochain formation of amelogenin molecules is well described by a combination of translational and rotational diffusion of a simplified anisotropic bipolar model consisting of hydrophobic spherical colloid particles and a point charge located on each particle surface. The colloid particles interact via a standard depletion attraction whereas the point charges interact through a screened Coulomb repulsion. We study the kinetics via a Brownian dynamics simulation of both translational and rotational motions and show that the anisotropy brought in by the charge dramatically affects the kinetic pathway of cluster formation and our simple model captures the main features of the experimental observations.     

Diffusion-based DNA target colocalization by thermodynamic mechanisms

In eukaryotic cell nuclei, a variety of DNA interactions with nuclear elements occur, which, in combination with intra- and inter-chromosomal cross-talks, shape a functional 3D architecture. In some cases they are organized by active, i.e. actin/myosin, motors. More often, however, they have been related to passive diffusion mechanisms. Yet, the crucial questions on how DNA loci recognize their target and are reliably shuttled to their destination by Brownian diffusion are still open. Here, we complement the current experimental scenario by considering a physics model, in which the interaction between distant loci is mediated by diffusing bridging molecules. We show that, in such a system, the mechanism underlying target recognition and colocalization is a thermodynamic switch-like process (a phase transition) that only occurs if the concentration and affinity of binding molecules is above a threshold, or else stable contacts are not possible. We also briefly discuss the kinetics of this `passive-shuttling' process, as produced by random diffusion of DNA loci and their binders, and derive predictions based on the effects of genomic modifications and deletions.     

Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot type.