We introduce an efficient computational methodology employing the potentials of mean force approach for estimating the detailed three-dimensional ionic distributions around arbitrarily complex charged biomolecular structures for all monovalent salt concentrations of practical interest (e.g., 0.1-5.0 M NaCl). Such distributions are required for specifying thermodynamic and structure-specific features of ion-mediated interactions of charged proteins, DNA and RNA, membranes, and macromolecular assemblies. As a first application, we present results for distributions around the B and Z, conformers of the DNA oligomer d(C-G)j8sd(C-G)js. The ionic microenvironment depends strongly on the DNA conformation, sequence, and bulk salt concentration.The computation of ionic distributions around highly charged biomolecules has been a subject ofcontinuous interest for the past 20 years due to the importance of the ionic environment for structural stability and transitions as well as for intermolecular.interactions and binding equilibria. For example, in the case of DNA, ionic effects play a central role in the extensively studied B -+ Z transition (for recent reviews, see refs. 1 and 2), at all levels of DNA structural organization (refs. 3 and 4 and references therein), and in ligand-DNA binding equilibria (5).The approaches proposed to date for computing ionic distributions around DNA vary widely with respect to (i) the structural modeling of the DNA-water-ions system, (ii) the statistical approximations employed, (iii) the computational complexity, (iv) the range of bulk ionic concentrations, and (v) the DNA conformations studied. If the water, ions, and DNA components are described microscopically, the only currently feasible way to determine ionic distributions is by large-scale computer simulations (i.e., Monte Carlo and molecular dynamics), as reported by Clementi and Corongiu (6-8) and more recently by van Gunsteren et al. (9). These calculations are extremely costly and can only be performed practically in the absence of added salt. For example, to simulate a 0.1 M NaCl electrolyte around DNA, one has to consider at least 50 anion-cation pairs and 25,000 water molecules in the central simulation box, an effort that is at present impossible.In previous studies water has been represented as a dielectric continuum and the ions have been considered either as hard or soft spheres or simply point particles. The treatments have substantial differences in the structural modeling of DNA and the use of statistical approximations. If one pictures a DNA structure as an infinitely long uniformly charged cylinder, it is straightforward to compute the average cylindrically symmetric ionic distributions through the numerical solution of the Poisson-Boltzmann (10-12) and hypernetted chain equations (13,14) This brief overview clearly shows that a general, fairly accurate, computationally efficient method to approximately solve the problem at hand for arbitrarily complex biomolecular structures over a wide range of bulk ionic ...