  We present a semi-quantitative model of condensation of short nucleic acid
(NA) duplexes induced by tri-valent cobalt(III) hexammine (CoHex) ions. The
model is based on partitioning of bound counterion distribution around singleNA
duplex into "external" and "internal" ion binding shells distinguished by the
proximity to duplex helical axis. In the aggregated phase the shells overlap,
which leads to significantly increased attraction of CoHex ions in these
overlaps with the neighboring duplexes. The duplex aggregation free energy is
decomposed into attractive and repulsive components in such a way that they can
be represented by simple analytical expressions with parameters derived from
molecular dynamic (MD) simulations and numerical solutions of Poisson equation.
The short-range interactions described by the attractive term depend on the
fractions of bound ions in the overlapping shells and affinity of CoHex to the
"external" shell of nearly neutralized duplex. The repulsive components of the
free energy are duplex configurational entropy loss upon the aggregation and
the electrostatic repulsion of the duplexes that remains after neutralization
by bound CoHex ions. The estimates of the aggregation free energy are
consistent with the experimental range of NA duplex condensation propensities,
including the unusually poor condensation of RNA structures and subtle sequence
effects upon DNA condensation. The model predicts that, in contrast to DNA, RNA
duplexes may condense into tighter packed aggregates with a higher degree of
duplex neutralization. The model also predicts that longer NA fragments will
condense more readily than shorter ones. The ability of this model to explain
experimentally observed trends in NA condensation, lends support to proposed NA
condensation picture based on the multivalent "ion binding shells".
