The sixth transmembrane segment (TM6) of the CFTR chloride channel has been intensively investigated. The effects of amino acid substitutions and chemical modification of engineered cysteines (cysteine scanning) on channel properties strongly suggest that TM6 is a key component of the anion-conducting pore, but previous cysteine-scanning studies of TM6 have produced conflicting results. Our aim was to resolve these conflicts by combining a screening strategy based on multiple, thiol-directed probes with molecular modeling of the pore. CFTR constructs were screened for reactivity toward both channel-permeant and channel-impermeant thiol-directed reagents, and patterns of reactivity in TM6 were mapped onto two new, molecular models of the CFTR pore: one based on homology modeling using Sav1866 as the template and a second derived from the first by molecular dynamics simulation. Comparison of the pattern of cysteine reactivity with model predictions suggests that nonreactive sites are those where the TM6 side chains are occluded by other TMs. Reactive sites, in contrast, are generally situated such that the respective amino acid side chains either project into the predicted pore or lie within a predicted extracellular loop. Sites where engineered cysteines react with both channel-permeant and channel-impermeant probes occupy the outermost extent of TM6 or the predicted TM5−6 loop. Sites where cysteine reactivity is limited to channel-permeant probes occupy more cytoplasmic locations. The results provide an initial validation of two, new molecular models for CFTR and suggest that molecular dynamics simulation will be a useful tool for unraveling the structural basis of anion conduction by CFTR.

Mitochondrial DNA levels are significantly decreased in patients with symptomatic, nucleoside-related hyperlactatemia, an effect that resolves on the discontinuation of therapy.

We developed molecular models for the CFTR chloride channel based on the prokaryotic ABC transporter, Sav1866. Here we analyze predicted pore geometry and side-chain orientations for TMs 3, 6, 9 and 12; with particular attention to the location of the rate-limiting barrier for anion conduction. Side-chain orientations assayed by cysteine scanning were found to be from 77% to 90% in accord with model predictions. The predicted geometry of the anion conduction path was defined by a space-filling model of the pore and confirmed by visualizing the distribution of water molecules from a molecular dynamics (MD) simulation. Pore shape is that of an asymmetric hour glass, comprising a shallow outward-facing vestibule that tapers rapidly toward a narrow “bottleneck” linking the outer vestibule to a large inner cavity extending toward the cytoplasmic extent of the lipid bilayer. The junction between the outer vestibule and the bottleneck features an outward–facing rim marked by T338 in TM6 and I1131 in TM12, consistent with the observation that cysteines at both of these locations reacted with both channel-permeant and channel-impermeant, thiol-directed reagents. Conversely, cysteines substituted for S341 in TM6 or T1134 in TM12, predicted by the model to lie below the rim of the bottleneck, were found to react exclusively with channel-permeant reagents applied from the extracellular side. The predicted dimensions of the bottleneck are consistent with the demonstrated permeation of Cl− pseudohalide anions, water and urea.