DNA renaturation at the water-phenol interface

Effect of Topology on the Statics and Dynamics of a Polymer Chain at the Fluid-Fluid Interface: A Molecular Dynamics Simulation Study

The effect of chain topology on the statics and dynamics of a chain at the interface of two immiscible fluids is studied by means of molecular dynamic simulations. For three topologically different chains, namely, linear, ring, and trefoil-knot, of the same molecular weight, the effect of varying both polymer-fluid and fluid-fluid interaction nature on the width of the fluid interface, chain conformation, shape, and chain dynamics is investigated. For a sharp-interface binary-fluid system, the interface width is insensitive to both topology and polymer-fluid interaction nature, while a weak nonmonotonic variation is seen for a broad-interface system. Chain extension normal to the interface plane is significantly affected by the topology with a trefoil-knot chain, due to the additional constraint, which has the largest value compared to both linear and ring polymers. Instantaneous shapes are also quantified through shape parameters. Furthermore, it is observed that the qualitative behavior of the center-of-mass mean-square displacement (MSD) is independent of topology, i.e., all the chain types show the same diffusion exponent α ( ∼ 1). However, the self-diffusion constant depends on the topology and it is the largest for the trefoil-knot chain. An interesting observation pertaining to the early time behavior of monomeric MSD is that, within the subdiffusive regime, the values of α for different parameters (independent of topology) are grouped into two distinct ranges (0.52 - 0.59 and 0.62 - 0.67), which are related to the different chain conformation for the polymer-fluid interaction range below and above a threshold value equal to that of the self-interaction of the pure fluid phase.

Multichain adsorption at fluid interfaces: Amphiphilic homopolymers vs copolymers

HYPOTHESIS

At selective liquid-liquid interface, amphiphilic homopolymers, having groups with different affinity for the liquids in each monomer unit, would demonstrate higher occupation of the interfacial layer than copolymers with various distributions of groups and be advantageous as interface stabilizers.

EXPERIMENTS

By means of Langevin dynamics computer simulation, conformations of multiple chains of amphiphilic macromolecules adsorbed at the liquid-liquid interface were studied. Monomer units having different affinity for the liquids were distributed variously along the polymer chains. Homopolymers, amphiphilic at the level of an individual monomer unit, and copolymers with random, altermating and multiblock distribution of groups were considered. The surface coverage, structure of the layer, and spatial distribution of monomer units were investigated depending on the polymer concentration.

FINDINGS

Compared to copolymers with random, alternating and multiblock distributions of the groups, the interfacial layer concentration of amphiphilic homopolymer is about 1.5 times higher, the adsorbed layer is remarkably thinner, has membrane-like structure and is asymmetric with respect to interface boundary. Also, the adsorbed amphiphilic homopolymers form fewer loops and tails, most located on one side of the interface. This combination of properties is promising for practical application in modern self-assembling molecular devices.

Lamellae and parking garage structures in amphiphilic homopolymer brushes with different grafting densities

This article is devoted to the study of polymer layers of amphiphilic homopolymers tightly grafted to a flat surface at the nodes of a square lattice. It was shown that, due to the amphiphilicity of monomer units containing groups with different affinities, in a selective solvent, such layers form lamellae perpendicular to the grafting surface. The period of the lamellae depends on the grafting density and the quality of the solvent. The results are presented in the form of a state diagram in variables "the energy of attraction of the side groups" (effective solvent quality) and "the distance between the grafting points" (inversely proportional to the square root of the grafting density). The diagram contains the regions of stability of lamellae with significantly different periods, and a transitional area with a parking garage structure. The diagram is constructed by calculating the layer-by-layer structure factor and the angle of inclination of the lamellae in the slice. The calculations were performed for different sizes of the simulation box, and the most commensurate size was determined by a special procedure for each grafting density. The results may be interesting not only to specialists in polymer science but also to all those who investigate the processes of self-organization and rearrangement in dense systems.

Ultrasensitive detection of serum hepatitis B virus by coupling ultrafiltration DNA extraction with real-time PCR

Background

A simple and reliable DNA extraction of hepatitis B virus (HBV) is critical in developing an ultrasensitive detection method for HBV infection. Current commercially available serum Hepatitis B Virus (HBV) DNA extraction methods are time-consuming, expensive and/or require specialized equipment, which hinders wide adoption of clinical laboratories. This study offers a report on an ultrasensitive HBV DNA detection method by coupling serum HBV DNA extraction by ultrafiltration (UF) with real-time PCR (qPCR) detection.

Methods

Serum proteins were precipitated by phenol to release HBV DNA in the supernatant which was then transferred to the UF devices. The resultant DNA concentrate was eluted and released into qPCR pre-mixture. The UF-qPCR assay performance, including recovery rate, linearity, detection sensitivity, precision and diagnostic accuracy that compared to the CAP-CTM V2.0 assay by analyzing batched low viral load clinical samples was evaluated.

Results

The recovery rate of the UF-based HBV DNA extraction method was above 80%. The assay linearity was demonstrated with a slope of 0.95 and R² values of 0.99. Limit-of-detection (LOD) of the UF-qPCR assay was determined to be 12.1IU/ml. The coefficient of variation (CV) of HBV quantitation for high, low and limit titer samples was 2.28%, 5.77% and 25.59%, respectively. Accuracy of the UF-qPCR assay was confirmed with the reference panel, and there was a strong correlation between these two methods (R² = 0.55, p < 0.01).

Conclusions

The UF-qPCR assay is reliable, highly sensitive, affordable and time-saving, and the method can be used for ultrasensitive detection of serum HBV.

Theory on the Mechanism of DNA Renaturation: Stochastic Nucleation and Zipping

Renaturation of the complementary single strands of DNA is one of the important processes that requires better understanding in the view of molecular biology and biological physics. Here we develop a stochastic dynamical model on the DNA renaturation. According to our model there are at least three steps in the renaturation process viz. nonspecific-contact formation, correct-contact formation and nucleation, and zipping. Most of the earlier two-state models combined nucleation with nonspecific-contact formation step. In our model we suggest that it is considerably meaningful when we combine the nucleation with the zipping since nucleation is the initial step of zipping and nucleated and zipping molecules are indistinguishable. Nonspecific contact formation step is a pure three-dimensional diffusion controlled collision process. Whereas nucleation involves several rounds of one-dimensional slithering and internal displacement dynamics of one single strand of DNA on the other complementary strand in the process of searching for the correct-contact and then initiate nucleation. Upon nucleation, the stochastic zipping follows to generate a fully renatured double stranded DNA. It seems that the square-root dependency of the overall renaturation rate constant on the length of reacting single strands originates mainly from the geometric constraints in the diffusion controlled nonspecific-contact formation step. Further the inverse scaling of the renaturation rate on the viscosity of reaction medium also originates from nonspecific contact formation step. On the other hand the inverse scaling of the renaturation rate with the sequence complexity originates from the stochastic zipping which involves several rounds of crossing over the free-energy barrier at microscopic levels. When the sequence of renaturing single strands of DNA is repetitive with less complexity then the cooperative effects will not be noticeable since the parallel zipping will be a dominant enhancing factor. However for DNA strands with high sequence complexity and length one needs to consider the underlying cooperative effects both at microscopic and macroscopic levels to explain various scaling behaviours of the overall renaturation rate.

New strategy to create ultra-thin surface layer of grafted amphiphilic macromolecules

It was found first that macromolecules made of amphiphilic monomer units could form spontaneously an ultra-thin layer on the surface which the macromolecules are grafted to. The width of such layer is about double size of monomer unit consisting of hydrophilic A (repulsive) and hydrophobic (attractive) B beads. The hydrophilic A beads are connected in a polymer chain while hydrophobic B beads are attached to A beads of the backbone as side groups. Three characteristic regimes are distinguished. At low grafting density, the macromolecules form ultra-thin micelles of the shape changing with decrease of distance d between grafting points as following: circular micelles-prolonged micelles-inverse micelles-homogeneous bilayer. Those micelles have approximately constant height and specific top-down A-BB-A structure. At higher grafting density, the micelles start to appear above the single bilayer of amphiphilic macromolecules. The thickness of grafted layer in these cases is different in different regions of grafting surface. Only at rather high density of grafting, the height of macromolecular layer becomes uniform over the whole grafting surface. The study was performed by computer modeling experiments and confirmed in framework of analytical theory.

DNA base pair hybridization and water-mediated metastable structures studied by molecular dynamics simulations

The base pair hybridization of a DNA segment was studied using molecular dynamics simulation. The results show the obvious correlation between the probability of successful hybridization and the accessible surface area to water of two successive base pairs, including the unpaired base pair adjacent to paired base pair and this adjacent paired base pair. Importantly, two metastable structures in an A-T base pair were discovered by the analysis of the free energy landscape. Both structures involved addition of a water molecule to the linkage between the two nucleobases in one base pair. The existence of the metastable structures provide potential barriers to the Watson-Crick base pair, and numerical simulations show that those potential barriers can be surmounted by thermal fluctuations at higher temperatures. These studies contribute an important step toward the understanding of the mechanism in DNA hybridization, particularly the effect of temperature on DNA hybridization and polymerase chain reaction. These observations are expected to be helpful for facilitating experimental bio/nanotechnology designs involving fast hybridization.

A simple DNA gate motif for synthesizing large-scale circuits

The prospects of programming molecular systems to perform complex autonomous tasks have motivated research into the design of synthetic biochemical circuits. Of particular interest to us are cell-free nucleic acid systems that exploit non-covalent hybridization and strand displacement reactions to create cascades that implement digital and analogue circuits. To date, circuits involving at most tens of gates have been demonstrated experimentally. Here, we propose a simple DNA gate architecture that appears suitable for practical synthesis of large-scale circuits involving possibly thousands of gates.

Fast molecular beacon hybridization in organic solvents with improved target specificity

DNA hybridization is of tremendous importance in biology, bionanotechnology, and biophysics. Molecular beacons are engineered DNA hairpins with a fluorophore and a quencher labeled on each of the two ends. A target DNA can open the hairpin to give an increased fluorescence signal. To date, the majority of molecular beacon detections have been performed only in aqueous buffers. We describe herein DNA detection in nine different organic solvents, methanol, ethanol, isopropanol, acetonitrile, formamide, dimethylformamide (DMF), dimethyl sulfoxide (DMSO), ethylene glycol, and glycerol, varying each up to 75% (v/v). In comparison with detection in water, the detection in organic solvents showed several important features. First, the molecular beacon hybridizes to its target DNA in the presence of all nine solvents up to a certain percentage. Second, the rate of this hybridization was significantly faster in most organic solvents compared with water. For example, in 56% ethanol, the beacon showed a 70-fold rate enhancement. Third, the ability of the molecular beacon to discriminate single-base mismatch is still maintained. Lastly, the DNA melting temperature in the organic solvents showed a solvent concentration-dependent decrease. This study suggests that molecular beacons can be used for applications where organic solvents must be involved or organic solvents can be intentionally added to improve the molecular beacon performance.

Condensed DNA: condensing the concepts

DNA is stored in vivo in a highly compact, so-called condensed phase, where gene regulatory processes are governed by the intricate interplay between different states of DNA compaction. These systems often have surprising properties, which one would not predict from classical concepts of dilute solutions. The mechanistic details of DNA packing are essential for its functioning, as revealed by the recent developments coming from biochemistry, electrostatics, statistical mechanics, and molecular and cell biology. Different aspects of condensed DNA behavior are linked to each other, but the links are often hidden in the bulk of experimental and theoretical details. Here we try to condense some of these concepts and provide interconnections between the different fields. After a brief description of main experimental features of DNA condensation inside viruses, bacteria, eukaryotes and the test tube, main theoretical approaches for the description of these systems are presented. We end up with an extended discussion of the role of DNA condensation in the context of gene regulation and mention potential applications of DNA condensation in gene therapy and biotechnology.

Mechanism of thermal renaturation and hybridization of nucleic acids: Kramers' process and universality in Watson-Crick base pairing

Renaturation and hybridization reactions lead to the pairing of complementary single-stranded nucleic acids. We present here a theoretical investigation of the mechanism of these reactions in vitro under thermal conditions (dilute solutions of single-stranded chains, in the presence of molar concentrations of monovalent salts and at elevated temperatures). The mechanism follows a Kramers' process, whereby the complementary chains overcome a potential barrier through Brownian motion. The barrier originates from a single rate-limiting nucleation event in which the first complementary base pairs are formed. The reaction then proceeds through a fast growth of the double helix. For the DNA of bacteriophages T7, T4, and phiX174, as well as for Escherichia coli DNA, the bimolecular rate k2 of the reaction increases as a power law of the average degree of polymerization <N> of the reacting single-strands: k2 is proportional to <N> alpha. This relationship holds for 100 < or = <N> < or = 50,000 with an experimentally determined exponent alpha = 0.51 +/- 0.01. The length dependence results from a thermodynamic excluded-volume effect. The reacting single-stranded chains are predicted to be in universal good solvent conditions, and the scaling law is determined by the relevant equilibrium monomer contact probability. The value theoretically predicted for the exponent is alpha = 1 - nutheta2, where nu is Flory's swelling exponent (nu approximately equal 0.588), and theta2 is a critical exponent introduced by des Cloizeaux (theta2 approximately equal 0.82), yielding alpha = 0.52 +/- 0.01, in agreement with the experimental results.

Sequence effects in the melting and renaturation of short DNA oligonucleotides: structure and mechanistic pathways

The renaturation/denaturation of DNA oligonucleotides is characterized in the context of expanded ensemble (EXE) and transition path sampling (TPS) simulations. Free energy profiles have been determined from EXE for DNA sequences of varying composition, chain length, and ionic strength. TPS simulations within a Langevin dynamics formalism have been carried out to obtain further information of the transition state for renaturation. Simulation results reveal that free energy profiles are strikingly similar for the various DNA sequences considered in this work. Taking intact double-stranded DNA to have an extent of reaction ξ = 1.0, the maximum of the free energy profile appears at ξ≈0.15, corresponding to ∼2 base pairs. In terms of chain length, the free energy barrier of longer oligonucleotides (30 versus 15 base pairs) is higher and slightly narrower, due to increased sharpness associated with the transition. Low ionic strength tends to decrease free energy barriers, whereby increasing strand rigidity facilitates reassociation. Two mechanisms for DNA reassociation emerge from our analysis of the transition state ensemble. Repetitive sequences tend to reassociate through a non-specific pathway involving molecular slithering. In contrast, random sequences associate through a more restrictive pathway involving the formation of specific contacts, which then leads to overall molecular zippering. In both random and repetitive sequences, the distribution of contacts suggests that nucleation is favored for sites located within the middle region of the chain. The prevalent extent of reaction for the transition state is ξ≈0.25, and the critical size of the nucleus as obtained from our analysis involves ∼4 base pairs.

Faceting ionic shells into icosahedra via electrostatics

Shells of various viruses and other closed packed structures with spherical topology exhibit icosahedral symmetry because the surface of a sphere cannot be tiled without defects, and icosahedral symmetry yields the most symmetric configuration with the minimum number of defects. Icosahedral symmetry is different from icosahedral-shaped structures, which include some large viruses, cationic-anionic vesicles, and fullerenes. We present a faceting mechanism of ionic shells into icosahedral shapes that breaks icosahedral symmetry resulting from different arrangements of the charged components among the facets. These self-organized ionic structures may favor the formation of flat domains on curved surfaces. We show that icosahedral shapes without rotational symmetry can have lower energy than spheres with icosahedral symmetry caused by preferred bending directions in the planar ionic lattice. The ability to create icosahedral shapes without icosahedral symmetry may lead to the design of new functional materials. The electrostatically driven faceting mechanism we present here suggests that we can design faceted polyhedra with diverse symmetries by coassembling oppositely charged molecules of different stoichiometric ratios.

Compressing a rigid filament: buckling and cyclization

We study elastic properties of rigid filaments modeled as stiff chains shorter than their persistence length. By rigid filaments we mean that fluctuations around the optimal filament shape are weak and that low-order expansions (quadratic or quartic) in the deviation from the optimal shape are sufficient to describe them. Our main interest lies in the profiles of force vs. projected filament length, closure probability and weakly buckled states. Results may be relevant to experiments on self-assembled biological (microtubules, actin filaments) and synthetic (organo-gelators) filaments, carbon nanotubes and polymers grafted with strongly repelling side chains, some of which are discussed here.

Aggregation and adsorption at the air-water interface of bacteriophage phiX174 single-stranded DNA

We study the phase behavior of phage phiX174 single-stranded DNA in very dilute solutions in the presence of monovalent and multivalent salts, in both water (H(2)O) and heavy water (D(2)O). DNA solubility depends on the nature of the salts, their concentrations, and the nature of the solvent. The appearance of attractive interactions between the monomers of the DNA chains in the bulk of the solution is correlated with an adsorption of the chains at the air-water interface. We characterize this correlation in two types of aggregation processes: the condensation of DNA induced by the trivalent cation spermidine and its salting out in the presence of high concentrations (molar and above) of monovalent (sodium) cations, both in water and in heavy water. The overall solubility of single-stranded DNA is decreased in D(2)O compared to H(2)O, pointing to a role of DNA hydration in addition to electrostatic factors in the observed phase separations. DNA adsorption involves attractive van der Waals forces, and these forces are also operating in the bulk aggregation process.

Charged particles on surfaces: coexistence of dilute phases and periodic structures at interfaces

We consider a mixture of two immiscible oppositely charged molecules strongly adsorbed to an interface, with a neutral nonselective molecular background. We determine the coexistence between a high density ionic periodic phase and a dilute isotropic ionic phase. We use a strong segregation approach for the periodic phase and determine the one-loop free energy for the dilute phase. Lamellar and hexagonal patterns are calculated for different charge stoichiometries of the mixture. Molecular dynamics simulations exhibit the predicted phase behavior. The periodic length scale of the solid phase is found to scale as epsilon/(lB psi3/2), where psi is the effective charge density, lB is the Bjerrum length, and epsilon is the cohesive energy.

Effect of external fluctuations on the affinity-specificity negative correlation in DNA-probe interactions

We show that the site-specific interaction of a probe DNA with the template DNA can be well modeled as an unbiased random jump process, where the probe molecule first nonspecifically binds to the template DNA and then searches for the specific site via unbiased random jump motion on the template DNA. By investigating the effects of increasing the jump size, and the fluctuations in the position of the specific site and the fluctuations in the specific site interval on the affinity-specificity negative correlation, we show that (1) increasing the jump size will in turn increase the affinity of the probe toward its target site on the template DNA, however, with a limiting value--the maximum affinity condition; (2) the degree of supercoiling or condensation of the template DNA as well as the electrostatic interactions between the probe and the template in turn control the jump size associated with the dynamics of the probe on the template DNA; (3) under a maximum specificity condition (therefore with minimum affinity), by introducing an external fluctuation in the relative position of the target site on the template DNA with respect to the probe, one can still improve the affinity rate; (4) on the other hand, one can improve the specificity of the probe toward the target site on the template DNA by introducing external fluctuations in the target-site interval. Finally, we propose the design strategies and optimum experimental conditions to simultaneously enhance the affinity as well as the specificity of probe toward its target site on the template DNA.

Semiflexible amphiphilic polymers: Cylindrical-shaped, collagenlike, and toroidal structures

A coarse-grained model is used to study the conformational properties of semiflexible polymers with amphiphilic monomer units containing both hydrophilic and hydrophobic interaction sites. The hydrophobically driven conformational transitions are studied using molecular dynamics simulations for the chains of varying stiffness, as characterized by intrinsic Kuhn segment lengths that vary over a decade. It is shown that the energy of hydrophobic attraction required for the realization of the coil-to-globule transition increases with increasing chain stiffness. For rather stiff backbone, the coil-to-globule transition corresponds to a first order phase transition. We find that depending on the chain stiffness, a variety of thermodynamically stable anisometric chain morphologies are possible in a solvent selectively poor for hydrophobic sites of amphiphilic monomer units. For flexible chains, the amphiphilic polymer forms a cylindrical globule having blob structure with nearly spherical blobs. With increasing stiffness, the number of blobs composing the globule decreases and the shape of blobs transforms into elongated cylinder. Further increase in stiffness leads to compaction of macromolecules into a collagenlike structure when the chain folds itself several times and different strands wind round each other. In this state, the collagenlike structures coexist with toroidal globules, both conformations having approximately equal energies.

How proteins search for their specific sites on DNA: the role of DNA conformation

It is known since the early days of molecular biology that proteins locate their specific targets on DNA up to two orders-of-magnitude faster than the Smoluchowski three-dimensional diffusion rate. An accepted explanation of this fact is that proteins are nonspecifically adsorbed on DNA, and sliding along DNA provides for the faster one-dimensional search. Surprisingly, the role of DNA conformation was never considered in this context. In this article, we explicitly address the relative role of three-dimensional diffusion and one-dimensional sliding along coiled or globular DNA and the possibility of correlated readsorption of desorbed proteins. We have identified a wealth of new different scaling regimes. We also found the maximal possible acceleration of the reaction due to sliding. We found that the maximum on the rate-versus-ionic strength curve is asymmetric, and that sliding can lead not only to acceleration, but also in some regimes to dramatic deceleration of the reaction.

Hydrogel-Encapsulated Lipid Membranes

We have encapsulated free-standing planar lipid membranes in a polymer hydrogel via in situ photopolymerization. These encapsulated membranes were measured to last at least 5 days and provided a stable environment for single-molecule measurements of incorporated alpha-hemolysin channel proteins. The translocation of single-stranded DNA through alpha-hemolysin incorporated into these encapsulated membranes is slowed over 100 times relative to that of the gel-free state. These stabilized membranes, therefore, have application to the single-molecule analysis of DNA as well as broader application to biomolecular sensors, membrane protein biophysical studies, and drug discovery.

DNA reassociation using oscillating phenol emulsions

Reassociating double-stranded DNA from single-stranded components is necessary for many molecular genetics experiments. The choice of a DNA reassociation method is dictated by the complexity of the starting material. Reassociation of simple oligomers needs only slow cooling in an aqueous environment, whereas reannealing the many single-stranded DNAs of complex genomic mixtures requires both a phenol emulsion to accelerate DNA reassociation and dedicated equipment to maintain the emulsion. We present a method that is equally suitable for reassociating either simple or complex DNA mixtures. The Oscillating Phenol Emulsion Reassociation Technique (OsPERT) was primarily developed to prepare heteroduplex DNA from alkali-denatured high molecular weight human genomic DNA samples in which hundreds of thousands of fragments need to be reannealed, but the simplicity of the technique makes it practical for less demanding DNA reassociation applications.

Ligand-induced DNA condensation: choosing the model

We test and compare different models for ligand-induced DNA condensation. Using 14C-labeled spermidine3+, we measure the binding to condensed DNA at micromolar to molar polyamine concentrations. DNA aggregates at a critical polyamine concentration. Spermidine3+ binding becomes highly cooperative at the onset of aggregation. At higher concentrations, spermidine3+ binding to condensed DNA reaches a plateau with the degree of binding equal to 0.7 (NH(4+)/PO3-). Condensed DNA exists in a wide range of spermidine concentrations with the roughly constant degree of ligand binding. At greater concentrations, the degree of binding increases again. Further spermidine penetration between the double helices causes DNA resolubilization. We show that a simple two-state model without ligand-ligand interactions qualitatively predicts the reentrant aggregation-resolubilization behavior and the dependence on the ligand, Na+, and DNA concentrations. However, such models are inconsistent with the cooperative ligand binding to condensed DNA. Including the contact or long-range ligand-ligand interactions improves the coincidence with the experiments, if binding to condensed DNA is slightly more cooperative than to the starting DNA. For example, in the contact interaction model it is equivalent to an additional McGhee-von Hippel cooperativity parameter of approximately 2. Possible physical mechanisms for the observed cooperativity of ligand binding are discussed.