Hierarchies in eukaryotic genome organization: Insights from polymer theory and simulations.

Eukaryotic genomes possess an elaborate and dynamic higher-order structure within the limiting confines of the cell nucleus. Knowledge of the physical principles and the molecular machinery that govern the 3D organization of this structure and its regulation are key to understanding the relationship between genome structure and function. Elegant microscopy and chromosome conformation capture techniques supported by analysis based on polymer models are important steps in this direction.

Twist propagation in dinucleosome arrays.

We present a Monte Carlo simulation study of the distribution and propagation of twist from one DNA linker to another for a two-nucleosome array subjected to externally applied twist. A mesoscopic model of the array that incorporates nucleosome geometry along with the bending and twisting mechanics of the linkers is employed and external twist is applied in stepwise increments to mimic quasistatic twisting of chromatin fibers. Simulation results reveal that the magnitude and sign of the imposed and induced twist on contiguous linkers depend strongly on their relative orientation.

Coarse-Grained Brownian Dynamics Simulations of the 10-23 DNAzyme.

Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A.

A boundary element method/Brownian dynamics approach for simulating DNA electrophoresis in electrically insulating microfabricated devices.

We present an approach for merging boundary element method (BEM) solutions of the electric field in electrically insulating complex geometries with Brownian dynamics (BD) simulations of DNA electrophoresis therein. Although a rote application of the standard BEM algorithm proves inaccurate and prohibitively expensive, we show that regularization of the near-wall electric field and an updating scheme commensurate with the characteristic length scale of the BD simulation furnishes a robust, efficient simulation protocol. The accuracy of the BEM-BD method is verified by simulating lambda-DNA collisions with an isolated, insulating cylindrical obstacle and comparing the results with equivalent BD simulations that employ the exact solution for the electric field.

Brownian dynamics simulations of single-stranded DNA hairpins.

We present a Brownian dynamics model which we use to study the kinetics and thermodynamics of single-stranded DNA hairpins, gaining insights into the role of stem mismatches and the kinetics rates underlying the melting transition. The model is a base-backbone type in which the DNA bases and sugar-phosphate backbone are represented as single units (beads) in the context of the Brownian dynamics simulations. We employ a minimal number of bead-bead interactions, leading to a simple computational scheme.

Polymer deformation in Brownian ratchets: theory and molecular dynamics simulations.

We examine polymers in the presence of an applied asymmetric sawtooth (ratchet) potential which is periodically switched on and off, using molecular dynamics (MD) simulations with an explicit Lennard-Jones solvent. We show that the distribution of the center of mass for a polymer in a ratchet is relatively wide for potential well depths U0 on the order of several kBT. The application of the ratchet potential also deforms the polymer chains.

Combinatorial design of passive drug delivery platforms.

We introduce a novel computational approach to designing passive drug delivery systems based on porous materials such as hydrogels. Our approach uses three tools: a method to establish the exact release pattern from all possible loading sites inside a given hydrogel; a method to generate a large number of hydrogel structures to be tested numerically, and finally an optimization algorithm which leads to the selection of optimal hydrogel structures. Using this approach, we show that controlled release curves can be obtained by using a genetic algorithm for the optimization step.

The theory of DNA separation by capillary electrophoresis.

The Human Genome has been sequenced in large part owing to the invention of capillary electrophoresis. Although this technology has matured enough to allow such amazing achievements, the physical mechanisms at play during separation have yet to be completely understood and optimized. Recently, new separation regimes and new physical mechanisms have been investigated.

Theory of DNA electrophoresis (approximately 1999-2002(1/2)).

Over the last two decades, the introduction of new methods such as pulsed-field gel electrophoresis and capillary array electrophoresis has made it possible to map and sequence entire genomes, including our own. The development of these experimental methods has been helped by the progress of theoretical and computational sciences, and the interactions between these three modi operandi of modern science are still pushing the limits of our technologies. We now see a clear trend towards proteomics and microfluidic (even nanofluidic!) devices.