Gramicidin S

ABSTRACT The dissociation of gramicidin A (gA) channels into monomers is the simplest example of a channel gating process. The initial steps in this process are studied via a computational model that simulates the reaction coordinate for dimermonomer dissociation. The nonbonded interaction energy between the monomers is determined, allowing for their free relative translational and rotational motion. Lowest energy pathways and reaction coordinates of the gating process are determined. Partial rupture of the six hydrogen bonds (6HB) at the dimer junction takes place by coupling monomer rotation and lateral displacement. Coupling rotation with axial separation is far more expensive energetically. The transition state for channel dissociation occurs when monomers are displaced laterally by ~4-6 [Angstrom], separated by ~1.6-2 [Angstrom], and rotated by ~120°, breaking two hydrogen bonds. In membranes with significant hydrophobic mismatch there is a much greater likelihood of forming 4HB and possibly even 2HB states. In the 4HB state the pore remains fully open and conductive. However, transitions from the 6HB to 4HB and 4HB to 2HB states take place via intermediates in which the gA pore is closed and nonconductive. These lateral monomer displacements give rise to transitory pore occlusion at the dimer junction, which provides a rationale for fast closure events (flickers). Local dynamics of gA monomers also leads to lateral and rotational diffusion of the whole gA dimer, giving rise to diffusional rotation of the dimer about the channel axis.

INTRODUCTION

Ion channel function is characterized by three basic properties: permeability, selectivity, and gating (Hille, 2001). Permeability and selectivity are determined by specific interactions between the channel forming protein (or peptide) and the ion(s) and water(s) occupying its open transmembrane pore. With the availability of atomic level structures for an ever-increasing number of such proteins (Doyle et al., 1998; Murata et al, 2000; Sui et al., 2002; Fu et al., 2000; Dutzler et al., 2002, 2003; Kuo et al., 2003), theoretical study of these processes has become a major research focus (see Kuyucak et al., 2001, for a recent review of various approaches). Even though the rates of ion flow are typically 10^sup 6^-10^sup 7^ s^sup -1^, too slow to be directly simulated, considerable atomic level insight into these processes has been gained by molecular dynamics or Monte Carlo (MC) simulations that establish, with greater or lesser precision, the energetics of permeation. The simplification that permits such analyses to be successful is the existence of (relatively) simple reaction coordinates for permeation, ones that can be (more or less) determined by inspection of the channel structure.

Gating, the process that controls ionic flow across a membrane, is a far more difficult nut to crack. Gating displacements (leading to channel opening and closing) are much slower, typically requiring milliseconds or longer. Unlike protein folding, in which the assembly seeks the minimum of the energy landscape funnel, these conformational changes must involve specific transitions over saddles, for which the corresponding "reaction coordinates" may be far from obvious. Concerted movements of substantial segments of the protein are presumably involved during opening, closing, or inactivation. Cartoon models for such processes pervade channel biophysics (see Hille, 2001), but to date no theoretical study has provided a molecular level analysis of gating in even the simplest system, although recent studies on potassium channels (Jiang et al., 2002, 2003) and the mechanosensitive MscL channel (Martinac and Perozo, 2002) provide data that will surely form the basis of such investigations. Gramicidin A (gA), in its conducting head-to-head conformation, provides a sufficiently tractable system that permits undertaking such a study. But even here there are significant complications. The "reaction coordinate" involves coupled complex relative motion of the gA monomers and membrane influences, mainly due to hydrophobic mismatch, that affect the overall channel dissociation process.

The gA channel is formed by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet (Urry et al., 1971; Andersen, 1984). Upon association the assembly is stabilized by six junctional hydrogen bonds (O'Connell et al., 1990). Cifu et al. (1992) showed that the elementary conducting unit is a dimer; in the absence of intermonomer linkers, supramolecular association (the formation of coupled channels) is unlikely. There is an enormous body of experimental and theoretical work investigating the structure, selectivity, and conductance of gA channels (Neher et al., 1978; Urban et al., 1980; Andersen, 1983; Arseniev et al., 1986; Becker et al., 1992; Partenskii and Jordan, 1992; Roux and Karplus, 1994; Ketchem et al., 1997; Wallace, 1998; Townsley et al., 2001). Less attention has been paid to channel gating: the energetics and the "reaction coordinate" for dissociation and formation of the gA dimer. A number of studies have demonstrated the dependence of the formation and dissociation rates of gA channels on voltage (Bamberg and Benz, 1976; Sandblom et al., 2001), membrane thickness (Kolb and Bamberg, 1977; Elliot et al., 1983), ion concentration (Ring and Sandblom, 1988; Ring, 1992), and elastic properties of lipid bilayers (Neher and Eibl, 1977;Ring, 1996; Goulian et al., 1998; Lundbk and Andersen, 1999). To understand the origin of these various phenomena, a molecular mechanism for dimer formation and dissociation in lipid bilayer membranes is needed.

Gating presumably occurs via dissociation and association of the monomers with the closing transition triggered by breaking the dinner's stabilizing hydrogen bonds (Durkin et al., 1993; Lundbaek and Andersen, 1999). The deletion of a single hydrogen bond (HB) at the junction between the monomers destabilizes the g A dimer (by ~ 10 kJ mol^sup -1^) and reduces its conductance (Durkin et al., 1993). As the monomers can rotate relative to each other, formation of 4HB or 2HB dimers is conceivably possible (Lundbaek and Andersen, 1999). The transition state in dissociation of the normal (6HB) dimer is believed to be reached when two hydrogen bonds are broken and the monomers separate by ~ 1.6 [Angstrom] (Lundbaek and Andersen, 1999). The average lifetime of the conducting state in native gA channels is on the order of milliseconds to seconds (Elliot et al., 1983; Sigworth and Shenkel, 1988; Ring and Sandblom, 1988).

Although many aspects of gA structure and function can be studied using all-atom molecular dynamics simulations aimed at determining the exact microscopic details of channel behavior, determination of transition states and reaction pathways requires special methods. To do this we attempt to capture essential features of gA dynamics by simplifying the model system. Our approach is based on kinetic MC methodologies widely used in statistical and condensed matter physics to study phase transitions, polymer systems, and critical phenomena in alloys and magnets (Binder, 1992; Landau and Binder, 2000). To investigate gA dissociation in detail we present a new computational model, previously described in preliminary fashion (Miloshevsky and Jordan, 2003), to simulate reaction coordinates for the dissociation process. We calculate the nonbonded interaction energy between the monomers, which are permitted to undergo free translational and rotational motion, and determine the lowest energy pathways and the reaction coordinates for the initial steps in the process D(imef) [if and only if] 2M(onomer). The effects that the relative monomer rotation angle, the intermonomer separation distance, the lateral intermonomer displacement, and the monomer-monomer tilt angle have on the intermonomer interaction energy are analyzed. The potential of mean force (PMF) along the minimum energy pathway for dimer dissociation is calculated using the free energy perturbation method. We find that partial rupture of the hydrogen bonds at the dimer junction takes place via coupled monomer rotation and lateral displacement. Coupling rotation with axial separation is energetically far more expensive. Hydrophobic mismatch, due to differences between the length of the nonpolar portion of gA and the membrane bilayer's acyl chains, greatly aids formation of 4HB and possibly even 2HB states. Finally, we discuss rotational diffusion of the gA dimer around its channel axis, observed in MC runs.

COMPUTATIONAL MODEL

Nonbonded interaction energy between gA monomers is calculated using our Monte Carlo ion channel proteins code as follows. We use partial charges and van der Waals parameters from the CHARMM22 all-hydrogen force field (MacKerell et al., 1998), treat bulk water regions as continua with [epsilon] = 80, and immerse the gA monomers in a low dielectric ([epsilon] = 1) membrane slab. The choice of [epsilon] = 1, used for consistency with the CHARMM22 force field, where polarization effects are not explicitly included, somewhat overestimates electrostatic effects; consequently representative computations with a slab [epsilon] of 2, 4, and 10 were also carried out. The reaction field is treated by the method of images (Dorman et al., 1996). Van der Waals and electrostatic interactions are computed with no cutoff. Most simulations were based on two crystallographic structures (Ketchem et al., 1996; Townsley et al., 2001) (pdb entries 1MAG and 1JNO, resolved by solid-state and micellarNMR, respectively). At each MC step the full intermonomer interaction energy is computed directly. Calculating the image charge reaction field energy is time-consuming (with the 552 real gA charges; 10 generations of images yield 11,040 image charges). The total interaction energy arises from the direct interaction of the partially charged atoms of the two monomers (both electrostatic and van der Waals) and the interaction between real and image charges. To accelerate calculation we separate the dimer into monomers. We introduce a crucial simplification, to be justified later, holding atoms within each monomer fixed while the monomers are relatively mobile; thus the interaction energy of atoms in the individual monomers is invariant and need not be computed. Excluding these contributions greatly reduces the computational cost. Although the monomers are held rigid, they undergo free three-dimensional relative rotation and translation. The bulk water-low [epsilon] slab interfaces are chosen to be ~ 13 [Angstrom] from the center of the coordinate frame along the Z axis. Image planes are perpendicular to the Z axis and the outermost protein atoms, including their van der Waals radii, are included in the low [epsilon] membrane slab. The image planes move along the Z axis together with monomers as they separate axially, but their separation is always kept symmetric relative to the membrane slab.

When gramicidin is inserted in a membrane, its hydrophobic length typically differs from that of the lipid bilayer (Elliot et al, 1983). We compute the deformational energy associated with this "hydrophobic mismatch" in an harmonic approximation (Nielsen et al., 1998; Lundbaek and Andersen, 1999; Neustadt and Partenskii, 2002; Partenskii et al., 2003) as E^sub elst^ = cu^sup 2^, where c is an elastic constant, specific to a particular membrane, and u is the vertical displacement of the membrane from its nondeformed (channel free) state. The vertical displacement is related to the monomer-monomer separation distance d, where u = u^sub 0^ -0.5 × d, and u^sub 0^ is the deformation depth (the hydrophobic mismatch) of each monolayer when the channel is in its native (6HB) conformation. The membrane's deformational energy depends not only on d, but also on the contact angle between lipid and channel, conventionally expressed in terms of the contact slope s. The effect that the membrane's elastic deformation has on channel energetics is effectively described by the constrained boundary condition s - 0 (Huang, 1986; Lundbaek and Andersen, 1999; Partenskii et al., 2003). The appropriate effective elastic constants are c = 1.7 kT/[Angstrom]^sup 2^ for a dimyristoyl phosphatidylcholine (DMPC) membrane (u^sub 0^ = 1.65 [Angstrom]) and c = 0.95 kT/[Angstrom]^sup 2^ for a glyceryl monooleate (GMO) membrane (u^sub 0^) = 3.4 [Angstrom]) (Partenskii et al., 2003).

The influence that this elastic energy has on monomer reorientation is implemented using the Metropolis Monte Carlo technique (Metropolis et al., 1953), in which transition to a new state depends upon the energy change, [Delta]E; the acceptance criterion is that exp(-[Delta]E/kT) is greater than a randomly chosen number between 0 and 1. The elastic energy E^sup old^^sub elst^ is added to the total (electrostatic + vdW) interaction energy E^sup old^^sub tot^ between gA monomers. The elastic energy is only a function of their separation distance. Possible effects of monomer tilt or lateral motion on the elastic energy are neglected. As the relative monomer orientation and separation is changed, the intermonomer energy becomes E^sup new^^sub tot^. The elastic energy E^sup new^^sub elst^ is cu^sup 2^; the total energy change is [Delta]E = E^sup new^^sub tot^ + E^sup new^^sub elst^ - E^sup old^^sub tot^ - E^sup old^^sub elst^. We use the Metropolis criterion, with T - 300 K, to determine if the new configuration is accepted. As monomers separate the membrane's elastic energy decreases and the direct monomer-monomer interaction energy increases. Dimer behavior is determined by the competition.

Simulation of complex gating motions in proteins with standard Monte Carlo or molecular dynamics fails to sample conformations separated by high energy barriers. These cannot be overcome by direct methods since they describe very low probability configuration space domains. Special computational techniques such as umbrella sampling (Torrie and Valleau, 1974, 1977), constrained reaction coordinate dynamics for the simulation of rare events (Carter et al., 1989), or high temperature molecular dynamics (Rosso et al., 2002) are required to sample reactive trajectories, the rare but important dynamical pathways that bridge long-lived stable states, for example, native and excited states of the gA channel. In these techniques the reaction coordinate is constrained within narrow windows or fixed at certain values and individual simulations performed for each window or at each fixed coordinate value, a computationally expensive approach. The novelty of our Kinetic Monte Carlo Reaction Path Following technique is its efficiency for following structural evolution unidirectionally along the reaction coordinate, determining the lowest energy pathways and the reaction coordinates for the initial steps of dissociating the dimer into monomers. Monomer B is held fixed and monomer A (see Fig. 1) is allowed both translational (three Cartesian degrees of freedom) and rotational (rotation [phi] and tilt [theta] angles) motion. As long as tilt is small, [phi], and the third Euler angle, [psi], are almost equal; thus, if [psi] were freely variable, motion in that direction would immediately counter the rotation about [phi] and it would be impossible to monitor reaction. Consequently it was forbidden. The Metropolis method (Metropolis et al., 1953) allows moves to states of higher energy. The smaller the energy difference the greater the probability of an acceptable uphill move. Thus the rotation angle [phi] is only allowed to increase. Other degrees of freedom are unconstrained. The reaction coordinate (angle [phi]) evolves slowly relative to the other degrees of freedom. As monomer A approaches a saddle point with rapidly increasing energy, only small changes in [phi] are likely to be accepted. This describes a major aspect of the motion along the reaction pathway. However, it is incomplete since the energy fluctuates around the lowest-energy groove. A large set of closely related paths may be sampled within a minimum energy pathway. At the peak of the energy barriers separating the stable states the reaction pathway may bifurcate. For accepted configurations several MC trials are used to relax monomer A with the new angle [phi] fixed. The remaining degrees of freedom then relax fully in response to movement along the reaction coordinate, an especially important feature as monomer A evolves downhill, since any new configuration with lower energy is always accepted.

We validated this technique against results of a grid search (Leach, 2001). Since grid searching must be restricted to a very few degrees of freedom, comparisons were made with monomer A that only permitted two degrees of freedom, Z and [phi]. Z was unrestricted and only increases in [phi] were permitted; other degrees of freedom (lateral motion and tilt) were frozen. The two approaches yielded identical energy profiles over the whole [phi]-range.

RESULTS

The molecular model of Fig. 1 depicts the 1MAG gA structure derived by Ketchem et al. (1996). The lowest energy pathway was determined by [phi] rotation of monomer A by 360° from the 6HB state to the 4HB state, then to the 2HB state, finally returning to the 6HB state. The 6HB state is one where all the Os and Hs of each monomer's six junctional COs and NHs are all (roughly) equidistant from the complementary Hs and Os of the other monomer (essentially, [phi]

Fig. 2 presents an energy map for the 1MAG structure (the d,[phi]-map) as a function of the separation distance d and rotation angle [phi]. Lateral displacement and tilt are forbidden. The map, calculated by grid search (Leach, 2001), with d and [phi] steps of 0.05 [Angstrom] and 0.5° respectively, describes this two-dimensional energy surface away from the reaction path. The black parts of the energy map are regions sterically inaccessible to monomer A. There are three clear energy wells corresponding to 6HB, 4HB, and 2HB states; these are uncoupled in d,[phi]-space. Large energy barriers separate transition between the wells. For the 1MAG and 1JNO structures used here, the energy, separation distance, and rotation angle differences between the three states determined from d,[phi]-maps are presented in Table 1.

In our approach, the monomers are treated as rigid bodies; effects due to their flexibility are suppressed. To justify this, we rely on solid-state NMR observations (Tian et al., 1996; Tian and Cross, 1999) demonstrating no significant change in polypeptide gA structure even at the highest cation loading of the channel. The hydrogen-bonded COs and NHs in the polypeptide skeleton are not noticeably perturbed (Tian et al., 1996), suggesting that the monomers move as nearly rigid bodies. We considered all four experimental head-to-head dimer structures, corresponding to different conformations and peptide sequences. Fig. 3 shows the lowest energy contours on the d,[phi]-map as a function of [phi]; the monomer axes are coincident with both lateral motion and tilt suppressed. 1MAG (Ketchem et al., 1996) and 1JNO (Townsley et al., 2001) are for gramicidin A; 1JO3 (Townsley et al., 2001) is for gramicidin B, and 1JO4 (Townsley et al., 2001) is for gramicidin C. All exhibit qualitatively similar lowest energy pathways for dissociation; the basic profile is conserved. Differences most likely reflect conformational differences at TRP9 and sequence differences at residue 11, regions quite distant from the dimer junction. These profiles suggest that major aspects of gramicidin dissociation are basically independent of the structural details of the monomer. These observations are unaltered if all backbone carbonyl oxygens (COs) and amino hydrogens (NHs) are mobile. Calculations, not shown here, performed for the 1MAG structure, demonstrated that allowing backbone COs and NHs to librate and stretch has no significant effect on the intermonomer interaction energy. The three wells are still separated by huge energy barriers, and the monomers must undergo large (>4 [Angstrom]) axial separations to migrate from one well to another. As rotation occurs, Os and Hs from the separate strands come into opposition, with strong consequent repulsion, a point discussed at length in what follows.

Fig. 4 illustrates the effect of the membrane dielectric constant, [epsilon], on the 1MAG energy profile for axial displacement of monomer A along the Z axis, starting from the 6HB state. The rotation angle [phi] is fixed at zero. Increasing [epsilon] greatly alters the 6HB energy-well depth, with saturation at [epsilon] ~ 10. This well (~30 kT) is conserved even when electrostatics is suppressed and all partial charges on the protein atoms are turned off. The well remains deep and the profile is steep; at Z ~ 0 [Angstrom] the interaction is repulsive, but becomes attractive at Z > ~0.5 [Angstrom]. Thus, in equilibrium the strong electrostatic attraction between the monomers is counterbalanced by van der Waals repulsion. When the monomers separate >0.5 [Angstrom], van der Waals forces aid intermonomer binding. To estimate the importance of interaction with explicit waters, we augmented the model with two cylindrical regions containing ~300 waters each, with the bulk continuum regions moved appropriately outwards. In the 6HB state the net solvent contribution (both electrostatic and van der Waals) to the intermonomer binding force is very small (~1 kT/[Angstrom]), far less than the ~15 kT/[Angstrom] intermonomer van der Waals force. The elastic energy of a GMO membrane as a function of the intermonomer separation distance, d, is also shown in Fig. 4. In the 6HB state it is ~11 kT, gradually decreasing to zero with increasing d. The magnitude of the slope of the elastic and the intermonomer energies differ greatly in the separation range, d

Fig. 5, A and B, illustrate the lowest energy profiles for the 1MAG and 1JNO structures, corresponding to the reaction coordinate for dimer dissociation. A large number of MC runs were followed to determine the behavior of monomer A and characteristic results are presented. Here monomer A is permitted lateral displacement and tilt. Only reaction coordinates, representing a trajectory in five-dimensional space, were determined. Maps of energy surface topology away from these pathways were not determined; thus energetics away from the lowest energy path remains unknown. The profiles, although different in detail, exhibit substantial qualitative similarities. The importance of the additional degrees of freedom (lateral displacement and tilt) is evident from the difference between the MC curves and the d,[phi]-curves of Fig. 3. Energy barriers decreased significantly. The data shown incorporate hydrophobic mismatch; similar behavior is found with no mismatch. In 1MAG (Fig. 5 A) the barrier separating the 6HB and 4HB wells is ~50 kT; the well depth difference is ~18 kT. Intermonomer hydrogen bonds break much more easily if monomer A undergoes rotary motion with simultaneous lateral displacement rather than with direct axial separation. At the peak of the 6HB-4HB barrier the reaction pathway bifurcates. Path 1 leads to the 4HB state. There is switching between hydrogen bonds at the intermonomer junction (a rotational shift and the breaking of two hydrogen bonds) with an abrupt discontinuity in the separation distance (~1.6 [Angstrom], see Fig. 6 A). Along path 2 the monomers dissociate directly from the native 6HB state by coupled rotary and lateral motion with neither a separation distance discontinuity nor a rotational shift between hydrogen bonds. In 1JNO (Fig. 5 B) the picture is less sharply drawn. The 6HB-4HB barrier is now ~60 kT, and the well-depth difference is in the range of ~20-40 kT. For both structures the probability of realizing the different paths depends on the tilt permitted to monomer A, the membrane thickness, etc. However, in both the 1MAG and 1JNO structures, for [phi] > 150° (escape from the 4HB state) the energy fluctuates between 40 and 60 kT and the barriers between the 4HB and 2HB wells and 2HB and 6HB wells disappear (compare the MC curves with the d,[phi]-curve). Increasing the membrane [epsilon] to as much as 4 naturally lowers the energy differences; however, it has no qualitative effect on the results. In what follows, we only present results for which [epsilon] = 1.

Fig. 6, A and B, illustrate the angular dependence of the axial separation distance between monomers for the 1MAG and 1JNO structures, respectively. Zero separation corresponds to the native state, i.e., the fully open conducting channel. In the 6HB state the hydrogen-bonding pattern is ^sup 1^^sub A^H-^sup 5^^sub B^O (where ^sup 1^^sub A^H is the amino hydrogen of residue 1A, and ^sup 5^^sub B^O is the carbonyl oxygen of residue 5B, etc.), ^sup 1^^sub A^O-^sup 5^^sub B^H, ^sup 3^^sub A^H-^sup 3^^sub B^O, ^sup 3^^sub A^O-3^sup ^^sub B^H, ^sup 5^^sub A^H-^sup 1^^sub B^O, and ^sup 5^^sub A^O-^sup 1^^sub B^H. For the d, [phi]-path (neither lateral displacement nor tilt allowed) 1MAG monomers separate ~4 [Angstrom] in transiting from the 6HB to the 4HB state (Fig. 6 A), due to a large electrostatic barrier. At [phi] ~ 60° five COs and NHs from the individual monomers are in opposition: ^sup 1^^sub A^H-^sup 5^^sub B^H, ^sup 1^^sub A^O-^sup 3^^sub B^O, ^sup 3^^sub A^H-^sup 3^^sub B^H, ^sup 3^^sub A^O-^sup 1^^sub B^O, and ^sup 5^^sub A^H-^sup 1^^sub B^H. This barrier is unaltered even if H-atom and O-atom motions (libration and stretching of the NH and CO bonds) are permitted. The [phi]-dependence of the axial separation distance is very different when lateral displacement and tilt are allowed (MC curves).

For 1MAG in a GMO membrane depicted in MC runs 1 and 2 of Fig. 6 A (mismatch u^sub 0^ = 3.4 [Angstrom]) the axial separation exhibits three distinct levels, ~0 [Angstrom], ~1.6 [Angstrom], and ~3.4 [Angstrom] corresponding to 6HB, 4HB, and 2HB states, respectively. For both the 1MAG and 1JNO structures, the monomers' axial separation is insignificant for [phi]

Fig. 6 B presents data for 1JNO in a GMO membrane (hydrophobic mismatch 3.4 [Angstrom]). Usually in GMO, 1JNO undergoes a transition to the 4HB state with an abrupt separation distance discontinuity at [phi] ~ 90°. However, direct dissociation from the 6HB state is observed in MC run 1. The crucial determinant is the lateral displacement (see Fig. 7 B), which is ~2 [Angstrom] (MC run 1); the permeation pathway, once displaced, never reformed. In the absence of a hydrophobic mismatch, direct dissociation of 1JNO from the 6HB state was also observed (results not shown).

Hydrophobic mismatch aids formation of the 4HB state in both structures and the 2HB state in the 1MAG structure. When hydrophobic mismatch is taken into account, the axial separation abruptly increases by ~1.6-2.0 [Angstrom] (level 2), corresponding to formation of the 4HB state. There is a corresponding rotational shift of junctional hydrogen bonds; residues 5A and 5B no longer contribute to junctional stability. The new bonding pattern is ^sup 1^^sub A^H-^sup 3^^sub B^O, ^sup 1^^sub A^O-^sup 3^^sub B^H, ^sup 3^^sub A^H-^sup 1^^sub B^O, and ^sup 3^^sub A^O-^sup 1^^sub B^H. This pathway (level 2) is also observed at low probability when there is no mismatch. For a large mismatch (GMO), formation of a 2HB state was also observed for 1MAG with an axial intermonomer separation of ~3.4 [Angstrom] (level 3) at [phi] ~ 210°. There is a further rotational shift and breakage of junctional hydrogen bonds. The residual bonding pattern is ^sup 1^^sub A^H-^sup 1^^sub B^O and ^sup 1^^sub A^O-^sup 1^^sub B^H. For 1MAG the 2HB state was not observed in modeling behavior in a DMPC membrane, as the hydrophobic mismatch (u^sub 0^ = 1.65 [Angstrom]) is too small. For GMO the 2HB state did not form in the 1JNO structure (Fig. 6 B). Transient dissociation events from levels 2 and 3 were observed in some MC runs (MC run 2 illustrates dissociation from level 2) where the 1MAG monomers separate by ~6.8 [Angstrom]. In sum, channel behavior is sensitive to bilayer deformation for monomers separated by ~1.6 [Angstrom], in agreement with the assumption of Lundbaek and Andersen (1999).

Fig. 7, A and B, illustrate the lateral displacement between the monomers along the dissociative pathways for the 1MAG and 1JNO structures, respectively. The intermonomer junction is displaced sideways by -3.5 [Angstrom] (1MAG) and ~5 [Angstrom] (1JNO) at the peak of the barrier separating the 6HB and 4HB wells; here [phi] ~ 75° and the dimer pore is fully occluded. In the absence of any hydrophobic mismatch, dissociation takes place laterally for 1MAG with a high probability (MC run 3). Along this pathway, for [phi] > ~75°, separation occurs laterally without the open pore reforming. With increasing hydrophobic mismatch, there is a greatly increased probability that the monomers' lateral displacement again approaches zero (forming the 4HB state) with a correspondingly fully reopened pore. In 1JNO the open pore (the 4HB state) reformed in all but one MC run (MC run 1 of Fig. 7 B, where the lateral separation is ~2 [Angstrom] at the "4HB angle" of ~120°). For GMO the formation of 2HB state with a fully open pore in the 1MAG structure is observed (Fig. 7 A, MC run 1). The transition between the 4HB and 2HB states again involves sideways displacement of the intermonomer junction by ~4 [Angstrom] (Fig. 7 A). Three COs and NHs from the individual monomers are in opposition: ^sup 1^^sub A^H-^sup 3^^sub B^H, ^sup 1^^sub A^O-^sup 1^^sub B^O, and ^sup 3^^sub A^H-^sup 1^^sub B^, producing an electrostatic barrier. Thus, in 1MAG transitions between the 6HB, 4HB, and 2HB states occur with intermediate pore closures. In 1JNO behavior is similar, but without formation of a 2HB state (Fig. 7 B). The 2HB state of 1JNO was observed if the tilt angle was fixed at zero.

Fig. 8, A and B, illustrate axial intermonomer displacement along the dissociative pathways as a function of the lateral separation distance for the two structures. For 1MAG, at lateral displacements 3 [Angstrom].

Fig. 9 illustrates the tilt angle [theta] for monomer A of 1MAG. It fluctuates between 7° and 10°. In the absence of hydrophobic mismatch either path 1 or 2 (Fig. 5) can occur depending on the tilt angle constraints. If tilt is forbidden or near zero, path 1 is followed frequently (9 out of 10 MC runs). When tilt is freely variable direct dissociation from the 6HB state (path 2) is most probable. For the 1JNO structure, tilt angle fluctuations are substantially larger, fluctuating as much as 25° (results not shown).

The potential of mean force (PMF) for dimer dissociation was calculated via free energy perturbation (Leach, 2001). With [phi] as reaction coordinate, varied from [phi]^sub 0^ = 0° to [phi]^sub 1^ = 360°, intermediate angular states were described by a coupling parameter [lambda] ranging from 0 to 1. The intermediate angle [phi], corresponding to [lambda]^sub i^, is [phi]([lambda]^sub i^) = [lambda]^sub i^[phi]^sub 1^ + (1 - [lambda]^sub i^)[phi]^sub 0^. In each MC simulation [phi] was fixed at [phi]([lambda]^sub i^). The free energy differences [Delta]A([lambda]^sub i^ [arrow right] l^sub i+1^) and [Delta]A([lambda]^sub i^ [arrow right] l^sub i-1^) were obtained from a single MC simulation. Three configurations, [phi]([lambda]^sub i-1^), [phi]([lambda]^sub i^), and [phi]([lambda]^sub i+1^), are treated simultaneously by calculating the energy difference between the configuration [phi]([lambda]^sub i^) and its [delta][lambda]-incremented neighbors [phi]([lambda]^sub I^ - [delta][lambda]) and [phi]([lambda]^sub I^ + [delta][lambda]) at each perturbation step. Fig. 10 illustrates the potential of mean force along the reaction coordinate for dissociation of 1MAG. The reaction coordinate [phi] was varied in 0.3° steps. In the main the PMF mimics the total energy profiles along reaction pathways. In the 6HB region there is no difference between the PMF and total energy profiles. The PMF barriers separating the 6HB and 4HB states are lower than the total energy barriers.

DISCUSSION

Our results demonstrate that rupture of the six very strong hydrogen bonds between the monomers is not immediately dissociative. The monomers undergo rotary motion with simultaneous relative lateral displacement rather than direct axial separation. The reaction paths are five-dimensional trajectories that cannot be summarized in a single picture. Movies (viewed from within the membrane slab and along the channel axis) showing the monomers' dissociation behavior in a GMO membrane (mismatch 3.4 [Angstrom]) are available at the URL (http://people.brandeis.edu/gennady/ gA.html). The hydrogen bonds break slowly as the monomers undergo rotational and lateral shifts. Only after substantial rupture is the transition state reached, at which point there is the abrupt shift leading to hydrogen-bond reduction. After the first transition state is reached, the [phi]-dependence of the hydrogen-bonding pattern depends upon the path followed, i.e., upon the channel-bilayer hydrophobic mismatch. When two hydrogen bonds are broken, membrane elasticity greatly aids the transition from the native 6HB state to the 4HB state. This switching between hydrogen bonds at the intermonomer junction (a rotational shift by two hydrogen bonds) occurs abruptly with an axial separation distance discontinuity of ~1.6 [Angstrom], a model prediction in excellent agreement with the conclusion of Lundbaek and Andersen (1999).

Hydrophobic mismatch has no discernable effect on dimer behavior in the 6HB state as the d-dependence of the intermonomer energy is much larger than that of the elastic energy (see Fig. 4). Changes in the intermonomer interaction energy and the elastic energy of the membrane deformation do not counterbalance one another. For a 3.4 [Angstrom] mismatch in 1MAG the intermonomer energy increases by ~31 kT and the elastic energy decreases by ~3 kT as d increases from 0 to 1 [Angstrom]. Consequently the probability of observing transient dissociation events originating from the 6HB state is too low to be sampled. Even for the artificial case with electrostatics ignored ([epsilon] = [infinity]), the intermonomer energy varies much more rapidly than the elastic energy; the 6HB energy well remains deep and steep. Thus, in the 6HB state, relative monomer motion behavior is overwhelmingly likely to be governed by the intermonomer interaction energy. We calculated PDFs in 1MAG for intermonomer separation distances d, lateral displacements R, rotation angles [phi], and tilt angles [theta] for both monomers. Both monomers were totally mobile, permitted free translational and rotational motion. The PDFs determined the most probable values of d, R, [phi], and [theta], and their fluctuation range in the 6HB state. We studied the 1MAG dimer's 6HB state in three cases: 1), no mismatch; 2), mismatch, u^sub 0^ = 1.6 [Angstrom], DMPC-like; and 3), mismatch, u^sub 0^ = 3.4 [Angstrom], GMO-like. The channel-bilayer hydrophobic mismatch had no noticeable effect on the PDFs in the 6HB state. Accounting for membrane deformation (mismatches of 1.6 or 3.4 [Angstrom]) all fluctuations in the PDFs are within the MC error. This suggests that the gA channels prefer to remain dimeric (if initially dimeric) in membranes of reasonable thickness.

In the PDF determinations with both monomers freely mobile, direct fluctuations leading to the transition state (two broken hydrogen bonds and monomers separated by ~1.6 [Angstrom]) were never observed. This contrasts sharply with the transitory dissociative event seen in Fig. 6 A from level 2 (the 4HB state); nothing similar is seen originating from level 1 in either 1MAG or 1JNO. We conclude that direct dissociation (or attaining the transition state) from the 6HB state is a very rare event, too improbable to be observed by standard MC sampling of PDFs; thus the average lifetime of the open state of the gA dimer corresponds to the 6HB state transiting to unbound monomers. Thermally driven relative motion of the monomers coupled to the hydrophobic mismatch is not directly dissociative. Possibly the local concentration of bilayer's undulatory fluctuations that are eventually released leading to the gA transition state play a role in this process. Other factors encouraging fluctuations leading to instability might be the presence of ions in the channel, the transmembrane potential, etc. The channel-bilayer hydrophobic mismatch clearly affects the dissociation pathway when two hydrogen bonds are broken and the monomers separated at ~1.6-2.0 [Angstrom]. Full dissociation events were observed in a GMO-like membrane when the 1MAG monomers are no longer in the 6HB energy well (Fig. 6 A). However, instantaneous dimerization follows these dissociation events, which suggests that transition states with two broken hydrogen bonds are very short-lived compared to the 6HB state.

For a gA dimer inserted into GMO, hydrophobic mismatch greatly aids the formation of 4HB and even 2HB (in 1MAG) states. At the transition between the 6HB and 4HB states the monomers are laterally and rotationally shifted relative to each other by ~4-6 [Angstrom] and ~75°-90°, respectively. Their axial separation gradually increases above zero. In the absence of a mismatch, there would be a high probability for either the 1MAG or 1JNO dimers to dissociate laterally (without forming a 4HB state). When there is mismatch just the opposite happens; 4HB state formation is highly probable in both cases. In the 4HB state the monomers are rotationally shifted by ~140° and the corresponding axial separation distance is ~1.6 [Angstrom]. However, there is now no lateral displacement. The pore is perfectly open, as in the 6HB state. During the 6HB-to-4HB transition the pore is closed fleetingly due to lateral displacement (Fig. 8, A and B). There is experimental evidence (Mobashery et al., 1997) that increased membrane-channel mismatch leads to a new conducting state involving the same pair of monomers. We concluded, after discussion with O. S. Andersen (private communication), that this new conducting state could well be the 4HB state. Our analysis predicts similar behavior leading to formation of a 2HB state in a thick membrane; the pore would again be fully open at an axial separation distance of ~3.4 A. Again in the 4HB-to-2HB transition in 1MAG there is intermediate pore closure due to lateral displacement (Fig. 8 A). Thus, there can be three levels in the monomers' axial separation distance: ~0 [Angstrom], ~1.6 [Angstrom], and ~3.4 [Angstrom]. Direct dissociation due to the membrane contribution was observed only from the two last levels, not from the first level. Any differences between the results for 1MAG and 1JNO, seen in Figs. 5-8, probably reflect the latter conformer's less "stressed" backbone (Allen et al, 2003).

There is experimental evidence (Lee et al., 1993) that gA dimers undergo rotational diffusion around the bilayer normal. Our MC results demonstrate that rotational diffusion of the dimers is a consequence of local monomer dynamics. Movies of MC simulations performed for a 1MAG dimer inserted into a GMO membrane (mismatch 3.4 [Angstrom]) are available at the URL (http://people.brandeis.edu/gennady/ gARotation.html). Each monomer rotates and translates freely in three-dimensional space. These degrees of freedom (local dynamics) represent possible large amplitude motions of the rigid monomers. In the 6HB state the monomers only undergo rotary and sliding motion, with negligible fluctuations in axial separation. Rotational diffusion of the dimer arises naturally from local monomer dynamics. There is no preferred clockwise or anticlockwise rotation; dimer rotation occurs randomly in either direction. However, during this process the channel can make a full 360° turn around its axis, equivalent to rotational diffusion of the dimer as a whole. In addition to this uniaxial rotational diffusion, the dimer diffuses laterally as a unit. A recent solid-state NMR study (Mo et. al., 2004) of the monomeric gA closed state shows that the free monomer does not exhibit global axial rotation about the bilayer normal. This is consistent with our study; each monomer needs a partner for rotational diffusion to occur. For rotation axial angular momentum is required. Unless the membrane (or an external field) applies a torque, an isolated monomer will not rotate. However, in the dimer, the monomer motions are coupled. Their relative translational and rotational motions are highly restricted, manifesting themselves in fluctuations of the rotational angle [phi] and the lateral displacement R. Both dimer rotation around its axis and its translational diffusion in the membrane plane arise from the restricted translational and rotational monomer motions.

Our results provide a natural explanation for flickering during disappearance of double-barreled gA channels (Goforth et al., 2003; Rokitskaya et al., 2003). Lateral motion of either the 1MAG or 1JNO monomers along the dissociation pathway leads naturally to fluctuations in the pore cross-section at the dimer junction; this may rationalize abrupt transitions (flickers) on the sub-ms timescale. In going from the 6HB to the 4HB state, ionic conduction is partially interrupted and restored as the monomers undergo lateral fluctuations at the intermonomer junction. Armstrong and Cukierman (2002) have observed fast closure flickers, promoted by greater hydrophobic mismatch, for dioxolane-linked gramicidins. As these channels cannot separate, flickering was viewed as arising from bilayer undulations blocking the mouths of the channels. However, the dioxolane linkage does not preclude flicker due to rotation, lateral motion, and tilt; the membrane thickness effect is consistent with our observation of increased probability of occlusion with increased mismatch.

In our study both membrane and bulk water were treated as dielectric continua. For bulk water this is a good approximation. We carried out a series of MC runs and investigated the effect that explicit water molecules within the gA pore have on dissociation behavior for the d,[phi]-path (no lateral motion and tilt). The explicit water molecules reduce the energy difference between the wells and barriers; however, they do not alter the profile's basic shape. We also considered the possible influence of a permeant cation. A K+ ion located at mid-dimer increases the intermonomer attraction in agreement with numerical estimates and experimental results demonstrating channel stabilization by ion occupancy (Ring and Sandblom, 1988; Ring, 1992). With two K+ ions located at Z = 9 [Angstrom] and Z = -9 [Angstrom] (the binding sites close to the channel mouths) the intermonomer energy profile is slightly altered but its basic shape is conserved. In addition, the energy profile is sensitive to the value chosen for the dielectric constant of the membrane slab. However, in all the cases the shape of the intermonomer energy profile is conserved.

The effect that explicit water molecules and ions located within the pore had on the dissociation pathways, accounting for the monomer's lateral motion and tilt, could not be studied in full detail due to various technical problems. However, we monitored the initial stages of the dissociation process permitting both lateral motion and tilt while incorporating two K+ ions and seven water molecules in the gA pore, for 6HB, 4HB, and 2HB states, using the 1MAG structure. In the 6HB state, for [phi] up to ~75° the most favorable pathways are identical; with the pore empty or occupied, monomers separate laterally, not axially. By [phi] ~ 75°, the lateral separation, whether or not water is present, is almost 4 [Angstrom] (data not shown), so that here too the conductance pathway is occluded at the junction. For [phi] > ~75°, two or three water molecules leak into the membrane slab through the gap formed by the rotational shift and corresponding lateral displacement of the monomers, but the ions stay at their binding sites with the remaining water molecules. We found similar behavior for the 4HB and 2HB states, constructed with [phi] ~ 120° and [phi] ~ 240°, respectively. In both cases as [phi] increases from these initial values, the monomers again separate laterally, not axially; the conductance pathways become occluded at the junction as the lateral separations are again ~4 [Angstrom]. Although it is not possible to fully track the dissociative process, these studies confirmed the results presented in Figs. 6-8 for the empty channel. Even with an occupied pore, the initial steps in the dissociation mechanism are essentially the same, involving rotation and lateral displacement of the monomers. Neither the presence of water or ions favors dissociation via rotation and direct axial displacement. These results complement the study of Durkin et al. (1993) indicating that the gA monomers cannot be stabilized by the water molecules located in a defect at the monomer's junction and forming a hydrogen-bonded bridge between the free amide and carbonyl groups.

Interaction between the dissociating gA monomers and the surrounding explicit lipid molecules, not included in our study, may affect the dissociation pathways. The gA dimer dissociating by lateral monomer displacement shifts the lipids adjacent to the channel. The monomers may experience a lateral resistance from the surrounding lipid. However, experimental results (Tank et al., 1982; Borisenko et al., 2003) indicate that lateral diffusion of gramicidin monomers in the plane of the membrane is quite fast (D ~ 3 × 10^sup -8^ cm^sup 2^ s^sup -1^). This suggests that the monomers are relatively free to move sideways in the lipid bilayer (during 1 µs a gA monomer would travel ~35 [Angstrom]) so that the influence that the lipids adjacent to the channel have on the reaction pathways is likely insignificant.

Our fundamental observation, from which all results flow, is that for both 1MAG and 1JNO the initial step in the dissociation process involves lateral, not axial, separation of the monomers. The choice of the slab [epsilon] is unimportant; qualitative behavior is the same for all [epsilon]

CONCLUSIONS

We have presented a Monte Carlo study of the simplest channel gating process: the mechanism of the gA dimer to monomer dissociation reaction. By using a simplified model capturing essential features of the monomer's kinetics we have found the reaction pathways for gA dimer dissociation and investigated the effect of hydrophobic mismatch on these lowest energy paths. Our results reveal numerous important mechanistic features that are likely to be generally characteristic of the dimer dissociation process. Our predictions are in good agreement with much experimental data.

A large electrostatic barrier, separating the 6HB and 4HB states, arises as the monomers rotate and the interjunctional H-bonding scheme is dephased. To reduce this barrier, the monomers are displaced laterally. Axial separation is energetically far more expensive than this lateral motion. The transition state for channel dissociation is reached when the monomers are displaced laterally ~4-6 [Angstrom] and separated by ~1.6-2.0 [Angstrom] with a rotational shift by two hydrogen bonds. This prediction is supported by the study of Lundbaek and Andersen (1999). Lateral displacement fluctuations in the 6HB state provide a rationale for closure events (flickers) at sub-ms times.

Hydrophobic mismatch substantially stabilizes formation of 4HB and even possibly 2HB states. In transiting between the 6HB, 4HB, and 2HB states ionic conductance is interrupted as the monomers undergo a relative lateral shift at the intermonomer junction by ~4-6 [Angstrom]. However, when the 4HB state forms the pore is fully open and ionic conductance is restored. This prediction of our model is also supported by experimental evidence (Mobashery et al., 1997).

The dimer undergoes uniaxial rotational and lateral diffusion due to local monomer dynamics. This local dynamics represents possible large amplitude motions of the gA monomers leading to global axial rotation of the gA dimer. Experimental evidence for a uniaxial global gA dimer rotation has been demonstrated by Lee et al. (1993).

Note added in proof: Since this manuscript was submitted, Harms et al., (2003) used a new approach, single-molecule patch-clamp fluorescence microscopy, to probe conformational changes of gA channels. Their results strongly suggest that gA channel dynamics involves multiple open and closed states associated with different conformations of the gA dimer.

We thank Olaf S. Andersen (Weill Medical College of Cornell University) and Yuri N. Antonenko (Belozersky Institute of Moscow State University) for valuable suggestions and comments. We are grateful to Michael B. Partenskii (Brandeis University) for stimulating discussions. We also appreciate the care taken by the anonymous referees in assessing the article, and thank them for their valuable suggestions.

Work was supported by National Institutes of Health grant GM-28643.

REFERENCES

Allen, T. W., O. S. Andersen, and B. Roux. 2003. Structure of gramicidin A in a lipid bilayer environment determined using molecular dynamics simulations and solid-state NMR. J. Am. Chem. Soc. 125:9868-9877.

Andersen, O.S. 1983. Ion movement through gramicidin A channels. Studies on the diffusion-controlled association step. Biophys. J. 41:147-165.

Andersen, O. S. 1984. Gramicidin channels. Annu. Rev. Physiol. 46: 531-548.

Armstrong, K. M., and S. Cukierman. 2002. On the origin of closing flickers in gramicidin channels: a new hypothesis. Biophys. J. 82:1329-1337.

Arseniev, A. S., A. L. Lomize, I. L. Barsukov, and V. F. Bystrov. 1986. Gramicidin A transmembrane ion-channel. Three-dimensional structure reconstruction based on NMR spectroscopy and energy refinement. Biol. Membr. 3:1077-1104.

Bamberg, E., and R. Benz. 1976. Voltage-induced thickness changes of lipid bilayer membranes and the effect of an electric field on gramicidin A channel formation. Biochim. Biophys. Acta. 426:570-580.

Becker, M. D., R. E. Koeppe II, and O. S. Andersen. 1992. Amino acid substitutions and ion channel function. Model-dependent conclusions. Biophys. J. 62:25-27.

Binder, K. 1992. The Monte Carlo Method in Condensed Matter Physics. Springer-Verlag, Berlin, Germany.

Borisenko, V., T. Lougheed, J. Hesse, E. Fureder-Kitzmuller, N. Fertig, J. C. Behrends, G. A. Woolley, and G. Schutz. 2003. Simultaneous optical and electrical recording of single gramicidin channels. Biophys. J. 84: 612-622.

Carter, E. A., G. Ciccolti, J. T. Hynes, and R. Kapral. 1989. Constrained reaction coordinate dynamics for the simulation of rare events. Chem. Phys. Lett. 156:472-477.

Cifu, A. S., R. E. Koeppe II, and O. S. Andersen. 1992. On the supramolecular organization of gramicidin channels. The elementary conducting unit is a dimer. Biophys. J. 61:189-203.

Dorman, V., M. B. Partenskii, and P. C. Jordan. 1996. A semi-microscopic Monte Carlo study of permeation energetics in a gramicidin-like channel: the origin of cation selectivity. Biophys. J. 70:121-134.

Doyle, D. A., J. M. Cabrai, R. A. Pfuetzner, A. Kuo, J. M. Gulbis, S. L. Cohen, B. T. Chait, and R. MacKinnon. 1998. The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science. 280:69-77.

Dutzler, R., E. B. Campbell, M. Cadene, B. T. Chait, and R. MacKinnon. 2002. X-ray structure of a ClC chloride channel at 3.0 [Angstrom] reveals the molecular basis of anion selectivity. Nature. 415:287-294.

Dutzler, R., E. B. Campbell, and R. MacKinnon. 2003. Gating the selectivity filter in ClC chloride channels. Science. 300:108-112.

Durkin, J. T., L. L. Providence, R. E. Koeppe II, and O. S. Andersen. 1993. Energetics of heterodimer formation among gramicidin analogues with an NH^sub 2^-terminal addition or deletion: consequences of a missing residue at the join in channel. J. Mol. Biol. 231:1102-1121.

Elliot, J. R., D. Neddham, J. P. Dilger, and D. A. Haydon. 1983. The effects of bilayer thickness and tension on gramicidin single-channel lifetime. Biochim. Biophys. Acta. 735:95-103.

Fu, D., A. Libson, L. J. W. Miercke, C. Weitzman, P. Nollert, J. Krucinski, and R. M. Stroud. 2000. Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 290:481-486.

Goforth, R. L., K. C. Aung, V. G. Denise, L. P. Lyndon, R. E. Koeppe II, and O. S. Andersen. 2003. Hydrophobic coupling of lipid bilayer energetics to channel function. J. Gen. Physiol. 121:477-493.

Goulian, M., O. N. Mesquita, D. K. Fygenson, C. Nielsen, O. S. Andersen, and A. Libacher. 1998. Gramicidin channel kinetics under tension. Biophys. J. 74:328-337.

Harms, G. S., G. Orr, M. Montal, B. D. Thrall, S. D. Colson, and H. P. Lu. 2003. Probing conformational changes of gramicidin ion channels by single-molecule patch-clamp fluorescence microscopy. Biophys. J. 85:1826-1838.

He, K., S. J. Ludtke, Y. Wu, H. W. Huang, O. S. Andersen, D. Greathouse, and R. E. Koeppe. 1994. Closed state of gramicidin channel detected by x-ray in-plane scattering. Biophys. Chem. 49:83-89.

Hille, B. 2001. Ionic Channels of Excitable Membranes, 3rd ed. Simmer Associates, Sunderland, MA.

Huang, H. W. 1986. Deformation free energy of bilayer membrane and its effect on gramicidin channel lifetime. Biophys. J. 50:1061-1070.

Jiang, Y., A. Lee, J. Chen, M. Cadene, B. T. Chait, and R. MacKinnon. 2002. The open pore conformation of potassium channels. Nature. 417:515-522.

Jiang, Y., A. Lee, J. Chen, V. Ruta, M. Cadene, B. T. Chait, and R. MacKinnon. 2003. X-ray structure of a voltage-dependent K+ channel. Nature. 423:33-41.

Ketchem, R. R., K. C. Lee, S. Huo, and T. A. Cross. 1996. Macromolecular structural elucidation with solid-state NMR-derived orientational constraints. J. Biomol. NMR. 8:1-14.

Ketchem, R. R., B. Roux, and T. A. Cross. 1997. High-resolution polypeptide structure in a lamellar phase lipid environment from solidstate NMR-derived orientational constraints. Structure. 5:1655-1669.

Kolb, H. A., and E. Bamberg. 1977. Influence of membrane thickness and ion concentration on the properties of the gramicidin A channel. Biochim. Biophys. Acta. 464:127-141.

Kuo, A., J. M. Gulbis, J. F. Antcliff, T. Rahman, E. D. Lowe, J. Zimmer, J. Cuthbertson, F. M. Ashcroft, T. Ezaki, and D. A. Doyle. 2003. Crystal structure of the potassium channel KirBac1.1 in the closed state. Science. 300:1922-1926.

Kuyucak, S., O. S. Andersen, and S. H. Chung. 2001. Models of permeation in ion channels. Rep. Prog. Phys. 64:1427-1472.

Landau, D. P., and K. Binder. 2000. A Guide to Monte Carlo Simulations in Statistical Physics. University Press, Cambridge, UK.

Leach, A. R. 2001. Molecular Modelling: Principles and Applications. Prentice Hall, Harlow, England and New York.

Lee, K.-C., W. Hu, and T. A. Cross. 1993. ^sup 2HB^NMR determination of the global correlation time of the gramicidin channel in a lipid bilayer. Biophys. J. 65:1162-1167.

Li, G., and Q. Cui. 2002. A coarse-grained normal mode approach for macromolecules: an efficient implementation and application to Ca^sup 2+^ ATPase. Biophys. J. 83:2457-2474.

Lundbaek, J. A., and O. S. Andersen. 1999. Spring constants for channel-induced lipid bilayer deformations. Estimates using gramicidin channels. Biophys. J. 76:889-895.

MacKerell, A. D. Jr., D. Bashford, M. Bellolt, R. L. Dunbrack Jr., J. D. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph-McCarthy, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, W. E. Reiher, B. Roux, M. Schlenkrich, J. C. Smith, R. Stole, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, D. Yin, and M. Karplus. 1998. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B. 102:3586-3616.

Martinac, B., and E. Perozo. 2002. Physical principles of mechanosensitive channel gating by bilayer deformation forces. J. Physiol. (London). 544:16S-17S.

Metropolis, N., A. Rosenbluth, M. Rosenbluth, A. Teller, and E. Teller. 1953. Equations of state calculations by fast computing machines. J. Chem. Phys. 21:1087-1091.

Miloshevsky, G. V., and P. C. Jordan. 2003. Interaction energy between gramicidin A monomers. Biophys. J. 84:247a.

Mo, Y., T. A. Cross, and W. Nerdal. 2004. Structural restraints and heterogeneous orientation of the gramicidin A channel closed state in lipid bilayers. Biophys. J. In press.

Mobashery, N., C. Nielsen, and O. Andersen. 1997. The conformational preference of gramicidin channels is a function of lipid bilayer thickness. FEBS Lett. 412:15-20.

Murata, K., K. Mitsuoka, T. Hirai, T. Walz, P. Agre, J. B. Heymann, A. Engel, and Y. Fujiyoshi. 2000. Structural determinants of water permeation through aquaporin-1. Nature. 407:599-605.

Neher, E., and H. J. Eibl. 1977. The influence of phospholipid polar groups on gramicidin A channels. Biochim. Biophys. Acta. 464:37-44.

Neher, E., J. Sandblom, and G. Eisenman. 1978. Ionic selectivity saturation and block in gramicidin A channels. II. Saturation behavior of single channel conductance and evidence for the existence of multiple binding sites in the channel. J. Membr. Biol. 40:97-116.

Neustadt, Y. S., and M. B. Partenskii. 2002. Effective spring constants for the elastically coupled insertions in membranes. ArXiv:physics. 0212038 V4. Published online December 19, 2002.

Nielsen, C., M. Goulian, and O. S. Andersen. 1998. Energetics of inclusion-induced bilayer deformations. Biophys. J. 74:1966-1983.

O'Connell, A. M., R. E. Koeppe II, and O. S. Andersen. 1990. Kinetics of gramicidin channel formation in lipid bilayers: transmembrane monomer association. Science. 250:1256-1259.

Partenskii, M. B., and P. C. Jordan. 1992. Theoretical perspectives on ion-channel electrostatics: continuum and microscopic approaches. Q. Rev. Biophys. 25:477-510.

Partenskii, M. B., G. V. Miloshevsky, and P. C. Jordan. 2003. Stabilization of ion channels due to membrane-mediated elastic interaction. J. Chem. Phys. 118:10306-10311.

Ring, A., and J. Sandblom. 1988. Modulation of gramicidin A open channel life time by ion occupancy. Biophys. J. 53:549-559.

Ring, A. 1992. Influence of ion occupancy and membrane deformation on gramicidin A stability in lipid membranes. Biophys. J. 61:1306-1315.

Ring, A. 1996. Gramicidin channel-induced lipid membrane deformation energy: influence of chain length and boundary conditions. Biochim. Biophys. Acta. 1278:147-159.

Rokitskaya, T. I., E. A. Kotova, and Y. N. Antonenko. 2003. Tandem gramicidin channels cross-linked by streptavidin. J. Gen. Physiol. 121:463-476.

Rosso, L., P. Minary, Z. Zhu, and M. E. Tuckerman. 2002. On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles. J. Chem. Phys. 116:4389-4402.

Roux, B., and M. Karplus. 1994. Molecular dynamics simulations of the gramicidin channel. Annu. Rev. Biophys. Biomol Struct. 23:731-761.

Sandblom, J., J. Galvanovskis, and B. Jilderos. 2001. Voltage-dependent formation of gramicidin channels in lipid bilayers. Biophys. J. 81: 827-837.

Sigworth, F. J., and S. Shenkel. 1988. Rapid gating events and current fluctuation in gramicidin A channels. Curr. Topics Membr. Transp. 33:113-130.

Sui, H., B. G. Han, J. K. Lee, P. Walian, and B. K. Jap. 2002. Structural basis of water-specific transport through AQP1 water channel. Nature. 414:872-878.

Tama, F., F. X. Gadea, O. Marques, and Y. -H. Sanejouand. 2000. Building-block approach for determining low-frequency normal modes of macromolecules. Prot. Struct. Funct. Gen. 41:1-7.

Tank, D. W., E. S. Wu, P. R. Meers, and W. W. Webb. 1982. Lateral diffusion of gramicidin C in phospholipid multibilayers. Effects of cholesterol and high gramicidin concentration. Biophys. J. 40:129-135.

Tian, F., K.-C. Lee, W. Hu, and T. A. Cross. 1996. Monovalent cation transport: lack of structural deformation upon cation binding. Biochemistry. 35:11959-11966.

Tian, F., and T. A. Cross. 1999. Cation transport: an example of structural based selectivity. J. Mol. Biol. 285:1993-2003.

Tirion, M. M. 1996. Large amplitude elastic motions in proteins from a single-parameter, atomic analysis. Phys. Rev. Lett. 77:1905-1908.

Torrie, G. M., and J. P. Valleau. 1974. Monte Carlo free energy estimates using non-Boltzmann sampling: application to the sub-critical Lennard-Jones fluid. Chem. Phys. Lett. 28:578-581.

Torrie, G. M., and J. P. Valleau. 1977. Nonphysical sampling distributions in Monte Carlo free-energy estimation: umbrella sampling. J. Comp. Phys. 23:187-199.

Townsley, L. E., W. A. Tucker, S. Sham, and J. F. Hinton. 2001. Structures of gramicidins A, B, and C incorporated into sodium dodecyl sulfate micelles. Biochemistry. 40:11676-11686.

Urban, B. W., S. B. Hladky, and D. A. Haydon. 1980. Ion movements in gramicidin pores. An example of single-file transport. Biochim. Biophys. Acta. 602:331-354.

Urry, D. W., M. C. Goodall, J. D. Glickson, and D. F. Mayers. 1971. The gramicidin A transmembrane channel characteristics of head to head dimerized (L, D) helices. Prot: Natl. Acad. Sci. USA. 68:1907-1911.

Wallace, B. A. 1998. Recent advances in the high resolution structures of bacterial channels: gramicidin A. J. Struct. Biol. 121:123-141.

Gennady V. Miloshevsky and Peter C. Jordan

Department of Chemistry, Brandeis University, Waltham, Massachusetts

Submitted July 23, 2003, and accepted for publication September 17, 2003.

Address reprint requests to Peter C. Jordan, Dept. of Chemistry, MS-015, Brandeis University, P.O. Box 5491 10, Waltham, MA 02454-9110. Tel: 781-736-2540; E-mail: jordan@brandeis.edu.

© 2004 by the Biophysical Society

0006-3495/04/01/92/13 $2.00

Copyright Biophysical Society Jan 2004
Provided by ProQuest Information and Learning Company. All rights Reserved