Dinucleotides as simple models of the base stacking-unstacking component of DNA 'breathing' mechanisms.

Regulatory protein access to the DNA duplex 'interior' depends on local DNA 'breathing' fluctuations, and the most fundamental of these are thermally-driven base stacking-unstacking interactions. The smallest DNA unit that can undergo such transitions is the dinucleotide, whose structural and dynamic properties are dominated by stacking, while the ion condensation, cooperative stacking and inter-base hydrogen-bonding present in duplex DNA are not involved. We use dApdA to study stacking-unstacking at the dinucleotide level because the fluctuations observed are likely to resemble those of larger DNA molecules, but in the absence of constraints introduced by cooperativity are likely to be more pronounced, and thus more accessible to measurement. We study these fluctuations with a combination of Molecular Dynamics simulations on the microsecond timescale and Markov State Model analyses, and validate our results by calculations of circular dichroism (CD) spectra, with results that agree well with the experimental spectra. Our analyses show that the CD spectrum of dApdA is defined by two distinct chiral conformations that correspond, respectively, to a Watson-Crick form and a hybrid form with one base in a Hoogsteen configuration. We find also that ionic structure and water orientation around dApdA play important roles in controlling its breathing fluctuations.

Stability assessment of chitosan-sodium hexametaphosphate capsules.

The assessment of the stability of capsules based on chitosan-sodium hexametaphosphate complex formation has been carried out using two independent methods--compression and osmotic swelling, and the influence of the preparation variables was evaluated. The formulation containing 1.5% core polymer (chitosan) and 1.5% oligophosphate, in the absence of salt or at low ionic strength (0.15% NaCl) was found to provide the best membrane resistance. A higher concentration of cross-linker (2.25%) produced stable capsules only in absence of electrolyte. Mannitol, a porogen added to the preparation solutions, did not affect the stability of the obtained membranes. At elevated polyol (1%) and cross-linker levels (2.25%), and at 0% salt, membranes with decreased elasticity were obtained, having lower compression and osmotic bursting values and lower deformation at the breaking points. A significant influence of salt amount on the capsule stability was also found. This was attributed to changes in the membrane formation process resulting in membranes with different thickness and structure. Membrane compression stability was found to be dependent on the pH of both oligophosphate and chitosan solutions, as well as on the reaction time. The bursting force values decreased for capsule diameters below 1.6 mm. The increased membrane/capsule volume ratio for the small capsules decreased the capsule deformation freedom and caused capsule rupture at low force values. The capsules made at low salt amounts showed very good storage stability over time and at elevated temperatures. The results demonstrated that the capsules could be formulated with controlled properties for various biomedical applications.

The global intrinsic curvature of archaeal and eubacterial genomes is mostly contained in their dinucleotide composition and is probably not an adaptation.

Until now, the genomic DNA of all eubacteria analyzed has been hyper-curved, its global intrinsic curvature being higher than that of a random sequence. In contrast, that rule failed for archaea or eukaryotes, which could be either hypo- or hyper-curved. The existence of the rule suggested that, at least for eubacteria, global intrinsic curvature is adaptive. However, the present results from analyzing 21 eubacterial and six archaeal genomes argue against adaptation. First, there are two eubacterial exceptions to the former rule. More significantly, we found that the dinucleotide composition of the genome alone (which lacks all sequence information) is enough to determine the genome curvature. Additional evidence against adaptation came from showing that the global curvature of bacterial genomes could not have evolved under either of two complementary models of curvature selection: (i) that curvature is selected locally from unbiased variability; (ii) that curvature is established globally through the selection of a curvature-altering mutational bias. We found that the observed relationship between curvature and dinucleotide composition is incompatible with model (i). We also found that, contrary to the predictions of model (ii), the dinucleo-tide compositions of bacterial genomes were not statistically special in their curvature-related properties (when compared to stochastically generated dinucleotide compositions).

A study of the hydration of deoxydinucleoside monophosphates containing thymine, uracil and its 5-halogen derivatives: Monte Carlo simulation.

An extensive Monte Carlo simulation of hydration of various conformations of the dinucleoside monophosphates (DNP), containing thymine, uracil and its 5-halogen derivatives has been performed. An anti-anti conformation is the most energetically stable one for each of the DNPs. In the majority of cases the energy preference is determined by water-water interaction. For other dimers conformational energy is the most important factor, or both the factors are of nearly equal importance. The introduction of the methyl group into the 5-position of uracil ring most noticeably influences the conformational energy and leads to the decrease of its stabilizing contribution to the total interaction energy. The introduction of halogen atoms increases the relative content of anti-syn and syn-anti conformations of DNPs as compared to the parent ones due to the formation of an energetically more favorable water structure around these conformations. A correlation is observed between the Monte Carlo results for the halogenated DNPs and their experimental photoproduct distribution. The data obtained demonstrates a sequence dependence in the photochemistry of the halogenated dinucleoside monophosphates.

Grand canonical ensemble Monte Carlo simulation of the dCpG/proflavine crystal hydrate.

The grand canonical ensemble Monte Carlo molecular simulation method is used to investigate hydration patterns in the crystal hydrate structure of the dCpG/proflavine intercalated complex. The objective of this study is to show by example that the recently advocated grand canonical ensemble simulation is a computationally efficient method for determining the positions of the hydrating water molecules in protein and nucleic acid structures. A detailed molecular simulation convergence analysis and an analogous comparison of the theoretical results with experiments clearly show that the grand ensemble simulations can be far more advantageous than the comparable canonical ensemble simulations.

A Monte Carlo method for generating structures of short single-stranded DNA sequences.

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by +/- 10 degrees in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3(7)) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150-325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CXnG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5' end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10(6) conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar-phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques.