A modeling algorithm is presented to compute simultaneously polymer conformations and ionic current, as single polymer molecules undergo translocation through protein channels. The method is based on a combination of Langevin dynamics for coarse-grained models of polymers and the Poisson-Nernst-Planck formalism for ionic current. For the illustrative example of ssDNA passing through the ␣-hemolysin pore, vivid details of conformational fluctuations of the polymer inside the vestibule and -barrel compartments of the protein pore, and their consequent effects on the translocation time and extent of blocked ionic current are presented. In addition to yielding insights into several experimentally reported puzzles, our simulations offer experimental strategies to sequence polymers more efficiently.T ranslocation of polymers through biological channels is very complex involving many machineries and is a fundamental step in many life processes. Although several essential features of translocation are richly documented in systems such as mRNP complex through nuclear pores (1-3), a simple system has only recently been identified for following the single-file passage of one isolated polymer through one channel (4-11). In this system, the channel is constituted by self-assembling heptamers of the Staphylococcus aureus ␣-hemolysin (␣HL) protein. The channel is assembled in a phospholipid bilayer, which offers a physical barrier, and the channel has an opening diameter of Ϸ1.4 nm at the narrowest constriction (12). A single-stranded polynucleotide, such as poly(deoxyadenylate) and poly(deoxycytidylate), is pulled through the channel by an externally applied voltage gradient across the channel in a solution of a strong electrolyte. The idea is that the ionic current through the channel caused by the passage of small ions of the electrolyte is blocked to a certain extent during the event of translocation of the polymer. It has been hoped that the extent and duration of the current blockade are unique signatures of the identity of the polymer, both in terms of the polymer's chemical characteristics and physical length.Even this simplest setup, where identical molecules undergo translocation, has generated several puzzling results. The distribution, P( ), of the duration of blockade of ionic current I b is very broad and appears to exhibit at least two peaks. In addition, there are several levels of ionic current blockade I b for the same molecule. It is standard practice in experimental investigations to combine the histograms of and I b (8). The resultant scatter plots always yield two groups of data even for monodisperse homopolymers.To gain insight into these puzzles, we have developed the following simulation. It is complementary to a flurry of theoretical activity (13-21), based on entropic barrier dynamics (22), all of which lead to a generic P( ) unlike in experiments. Although it is indeed desirable to perform the computation ab initio, the size of the system to be simulated is forbiddingly large to enable such a compu...