Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg2+ into cells, we examine early events in the relaxation of Mg2+ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg2+ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.
Magnesium homeostasis is essential for life. In humans, the misregulation of magnesium is implicated in stroke [1], heart disease [2], and diabetes [3]. Magnesium transport is also crucial for bacteria [4]. The movement of magnesium through cell membranes, like that of other ions, is accomplished by integral membrane proteins that provide selective permeability across the dielectric barrier of the lipid bilayer [5]. In bacteria, magnesium uptake is mediated by the CorA protein [6–9], which can substitute for its functional homologue in yeast mitochondria [10]. Electrophysiological data suggests that TmCorA is a channel, not a transporter [8]. Seven crystallographic structures exist for CorA, six of which are from Thermotoga maritima (TmCorA) [11–16]. These structures reveal a homopentamer in which 10-nm-long protomeric α-helices (the “stalk” helices) form a transmembrane (TM) pore through which magnesium is presumed to flow. This pore contains two hydrophobic constrictions: the “MM stretch” (MM), a 1. 9-nm-long constriction formed by pore-lining residues M291, L294, A298, and M302; and the “lower leucine constriction” (LC), a shorter steric bottleneck formed by the sidechain of L280 (Fig 1). Mutagenesis studies suggest that the MM, but not the LC, is involved in channel gating [17,18]. Hydrophobic gates are important for the function of many ion channels, including ligand-gated [19–22], voltage-gated [23,24], phosphorylation-gated [25], and mechanosensitive channels [26,27]. In all crystal structures of CorA, both hydrophobic constrictions are too narrow to be hydrated, suggesting that the channel is in its closed state. Remarkably, the pore extends beyond the relatively small TM domain into a much larger funnel-shaped domain that protrudes by 6 nm into the cytosol (Fig 1). At the far rim of this funnel, divalent cations are bound between cytosolic protomer interfaces (the divalent cation sensor or DCS; Fig 1). Based on crystallographic structures, it was hypothesized that divalent cation occupancy of the DCS regulates magnesium transport by controlling the pore' s diameter or electrostatic profile [11–13]. Recent studies suggest that divalent cation binding to the DCS locks TmCorA in a transport incompetent conformation and that loss of these cations leads to an open conformation of the channel [14,28], which may be asymmetric [14]. To investigate the allosteric regulation of pore opening, we previously conducted a molecular dynamics (MD) study of TmCorA in a hydrated lipid bilayer, either with or without Mg2+ ions in the DCS [29]. The MM remained dehydrated throughout a 110-ns MD simulation in the presence of regulatory ions, but became dilated and hydrated in one of two trajectories generated after these ions were removed [29]. Wetting of the MM involved an iris-like mechanism initiated by the rearrangement of the cytosolic domain interfaces and transmitted to the MM by the long pore-lining stalk helices. These findings suggest a model of allosterically-regulated hydrophobic gating whereby decreasing cytosolic magnesium concentration reduces magnesium occupancy in the DCS, leading to sudden wetting of the pore. Here, we use massive sampling to examine the statistical significance and the kinetics of this apparent hydrophobic gating process, and to assess whether it results in the open state of the channel. Scaling up computational sampling by two orders of magnitude, we compare hundreds of 35-ns MD simulations in which the DCS are either fully occupied or empty. From a total of 54 microseconds of sampling, we observe many events in which completely connected columns of water condense and evaporate in the MM. We show that this hydrophobic gate is more likely to become hydrated in the absence of regulatory ions and we quantify the kinetics of wetting and dewetting transitions. Finally, we show that the extent of hydration of the MM is correlated to a reduction in the free energy barrier for the permeation of Mg2+ ions, supporting the hypothesis that the MM is an allosterically-regulated hydrophobic gate.
To determine the kinetics of allosterically-regulated wetting in the MM (the presumed hydrophobic gate), we massively repeated MD simulations of TmCorA in the two limiting states of ionic regulation, namely, with all ten regulatory (DCS) binding sites either occupied or empty. The PDB: 2HN2 crystallographic structure of TmCorA [13] was embedded in a hydrated lipid bilayer and seven-hundred 35-ns MD simulations were conducted for each regulatory state. Initially, water filled most of the pore but the MM and the LC were completely dehydrated. Although the MM remained dehydrated in most of the simulations, it was hydrated at least part of the time in 12 and 25% (respectively 83 and 175) of the 700 simulations conducted respectively with and without regulatory ions (Table 1), allowing us to quantify the effect of regulatory ions on pore hydration. MM wetting is depicted as S1 Movie. In wetting transitions, water molecules entered the MM along the pore axis and, occasionally, through transient packing defects between pore-lining helices. The latter defects connected the middle of the pore directly to bulk water at the cytosolic membrane-water interface (S1 Fig). To characterize wetting, we computed the maximum distance between consecutive water oxygen atoms along the pore axis, zgap. Wetting in the MM or LC is defined as end-to-end hydration by a completely connected water column (zgap ≤ 0. 38 nm). The fraction of simulations in which the MM was wetted, P (zgap ≤ 0. 38 nm), is shown as a function of simulation time in Fig 2A. In the presence of regulatory ions, this fraction appears to stabilize around 0. 3% after 15 ns (Fig 2A). In contrast, the MM was wetted five times more often in the absence of regulatory ions, and the wetting probability increased from 0 to 2% through the entire length of the simulations (Fig 2A). Conversely, hydration of the LC increased from 0 to 50% independently of ionic occupancy of the DCS, although equilibration occurred more rapidly in the absence of regulatory ions (Fig 2B). To quantify the extent of hydration of the hydrophobic constrictions, we used two metrics, the length of the longest dehydrated stretch, zgap, and the number of water molecules, Nwat. As expected, these two metrics are strongly anti-correlated (S2A and S2E Fig). The normalized distributions of zgap (S2B and S2F Fig) and Nwat (S2C and S2G Fig) in the MM were similar for the two sets of simulations. Hydration defects were most common near M291, at the cytosolic end of the MM (S2D and S2H Fig). These findings are consistent with our previous simulations [29] and with the hypothesis of Lunin et al. that M291 and L294 form the main hydrophobic gate [11]. Finally, S3 Fig shows that hydration of the MM and LC pore constrictions are not correlated with each other. In both regulatory states, many of the wetting events were transient, half of them lasting less than 200 ps. For each simulation, Fig 3 shows the fraction of time that the MM was wetted as a function of the fraction of time that the MM contained >20 water molecules. This analysis leads us to identify eight simulations in which the MM was both stably and highly hydrated, all of which occurred in the absence of regulatory ions. These simulations are henceforth referred to as “stably superhydrated” (SSH). In the most hydrated SSH simulation, the initially dehydrated MM abruptly and durably filled with ≥20 water molecules, at one point surging to twice this amount (Fig 4A). This hydration was concurrent with an increase in mean MM pore diameter from 0. 4 to 0. 74 nm (Fig 4A–4C) and involved asymmetric reorientation of pore-lining sidechains away from the lumen, toward neighboring protomers (Fig 4D and 4E). To assess whether hydration was coupled to reorganization of the pore, we computed the pore diameter along its long axis. On average, the volume available in the MM did not substantially depend on regulatory ion occupancy in the DCS (Fig 5A). However, SSH simulations displayed an increase in MM diameter concurrent with wetting. The mean pore diameter in the MM occasionally fluctuated above 0. 68 nm (the diameter of a hexahydrated Mg2+ ion) in the absence of regulatory ions, but never in their presence (Fig 5A). Upon removal of these ions, there was, on average, a tightening of the pore at its extracellular end, near the proposed selectivity filter of the GMN motif [30] (Fig 5A). Regardless of regulatory ion occupancy, the mean pore diameter in the MM was linearly correlated to the extent of hydration (Fig 5B). It has been suggested that MM wetting could be triggered by the exposure of previously-hidden hydrophilic moieties to the pore lumen following axial rotation of the stalk helices [15,16]. To test this hypothesis, we analyzed protein-water hydrogen bond interactions in the MM (S1 Text and S7 and S8 Figs). Taken together, our analysis indicates that pore wetting is not due to a decrease in the hydrophobic character of the MM but rather that its likelihood increases with the volume of the hydrophobic stretch. Accordingly, the increase in both wetting probability and hydration number with CorA pore size is consistent with hydrophobic wetting in simple nanoscopic systems [31] and in hydrophobic pore analogs such as carbon nanotubes [32]. The sharp dependence of hydrophobic wetting on dynamic fluctuations in pore diameter makes this process well suited to gating in a biological channel. Distributions of Nwat are shown in Fig 6A. Despite the presence of wetting transitions, these distributions are unimodal and values of Nwat greater than 10 are rare. The most likely value of Nwat in the MM is 2 for simulations both with and without regulatory ions (Fig 6A). The distribution of Nwat is shifted to slightly larger values in the absence of regulatory ions, a shift that persists after the removal of the SSH simulations from this analysis (Fig 6A). The distribution of Nwat computed exclusively from the SSH simulations is bimodal, with maxima at Nwat = 3 and 20 (Fig 6A). Fig 6B shows the distribution of hydration number computed separately for the wetted state. In the SSH simulations, Nwat is distributed asymmetrically around 22, while the remaining simulations, with and without regulatory ions, sampled values distributed normally around 15 (σ = 3). Thus, in the wet state, the dependence of the extent of hydration on regulatory ions was entirely due to the eight SSH simulations. To characterize the kinetics of wetting and dewetting of the hydrophobic MM, survival probabilities, S (t), were computed separately for wet and dry states, either with or without regulatory ions (Fig 7). Although wet states are more stable in the absence of regulatory ions, the difference between survival probabilities obtained with and without regulatory ions is abrogated by removing the SSH simulations from this analysis (Fig 7A). Thus, an essential difference between the two sets of simulations is that removal of the regulatory ions led to a hydration state not observed on the 35-ns time-scale in the presence of regulatory ions. The two wetted states differ both in the extent of hydration and in their kinetic properties. Conversely, the dewetted state of the MM was destabilized by the removal of regulatory ions, an effect that persists after removing the SSH simulations from this analysis (Fig 7B). The survival probabilities of the wetted and dewetted MM states were fit to a double exponential decay function from which we computed the rate of channel wetting to be 8. 2×106 ± 7×105 s-1 and 4. 1×106 ± 3×105 s-1 in the absence and presence of regulatory ions, respectively (S1 Table). The estimates of wetting and dewetting rates were used to compute free energy differences between hydrated and dehydrated states (S2 Table), which we used to construct a model of hydrophobic gating in CorA (see Discussion). To assess whether the dilation and hydration of the hydrophobic stretch led to the open state of the channel, we used umbrella sampling (US) [33,34] to compute the free energy associated with translocation of a divalent cation throughout the pore of TmCorA based on multiple starting conformations from each of four distinct states. These states correspond to hydrated and dehydrated MM with and without regulatory ions. Strictly speaking, our US simulations were conducted out of equilibrium since the channel conformation is slowly relaxing (S9 Fig). As such, these profiles provide estimates of the work required for ion permeation given initial states of protein conformation, pore hydration, and regulatory ion occupancy. Multiple free energy profiles from each state are shown together as ensemble averages in Fig 8. The general features of these profiles are similar: a small barrier at z = 2. 6 nm, where the sidechain of N314 partially occludes the pore; a barrier of variable magnitude between −1<z<1 nm, which corresponds to the MM; a smaller barrier at z = −2. 3 nm, which corresponds to the LC; and a minimum at z = −3 nm, which corresponds to the lower of two pore Mg2+ binding sites. The main difference between the four ensemble-averaged free energy profiles shown in Fig 8 is the height of the barrier in the MM. When the MM is fully dehydrated, it presents a very large barrier of 45±5 kcal/mol irrespective of regulatory binding site occupancy (Fig 8). The size of this barrier reduced to 25±5 kcal/mol when the MM was wetted, again irrespective of bound regulatory ions (Fig 8). Thus, the size of the free energy barrier in the MM correlates with its hydration and not with the presence of regulatory ions. This conclusion is supported and clarified by Fig 9, in which the magnitude of the free energy barrier in the MM, ΔG‡, is plotted against values of Nwat from initial US conformations, Nwat (t0), or Nwat (US), the average value of Nwat from US simulations while the lumenal Mg2+ ion was in the MM. The magnitude of ΔG‡ is anticorrelated with both Nwat (t0) and Nwat (US) (Fig 9). Linear fits to all data yield ΔG‡ = 46–0. 82×Nwat (t0) kcal/mol and ΔG‡ = 92–1. 82×Nwat (US) kcal/mol, with r2 coefficients of 0. 80 and 0. 83 respectively. By extrapolating these linear relationships, we predict that reducing the magnitude of the barrier to ΔG‡ = 0 kcal/mol requires 56 or 50 water molecules based on the two above equations, respectively. A linear fit of the relation between average MM pore diameter, dpore¯, and Nwat in the absence of regulatory ions (Fig 5B), dpore¯ = 0. 37+0. 0067×Nwat nm, suggests that the aforementioned values of Nwat correspond to average pore diameters of 0. 71 to 0. 75 nm, slightly larger than the 0. 68 nm diameter of hexahydrated Mg2+. In turn, this result suggests that the barrier at the MM gate is essentially due to Mg2+ desolvation, including effects beyond its first hydration shell. Note that the predicted diameter at which the gating barrier vanishes is significantly smaller than the ~1. 2 nm diameter of the putative open state proposed by Dalmas et al. [28]. Furthermore, note that the requirement ΔG‡ = 0 is excessive as the channel may conduct Mg2+ ions with moderate positive values of ΔG‡. According to the model we present above, a pore diameter equal to that of hexahydrated Mg2+ corresponds to Nwat = 46 and ΔG‡ = 8 kcal/mol. Because pore hydration reduces the free energy barrier to Mg2+ flux, we assessed whether removing a (pentameric) bulky hydrophobic sidechain from the MM enhances its hydration. To this end, we conducted two-hundred 4-ns simulations for wild type and L294A TmCorA in three different wetted conformations extracted from our initial set of fourteen-hundred simulations. Consistent with the role of the MM as a hydrophobic gate, the L294A mutation enhances pore hydration both in the presence and in the absence of regulatory ions (Fig 10). In the presence of regulatory ions, L294A simulations initiated with Nwat = 25 stabilized at Nwat = 32 (σ = 4) over the last 3 ns (Fig 10C). Using the linear fit of Nwat (t0) to ΔG‡ derived above, we predict that the pore remains closed to Mg2+, with ΔG‡ = 20 kcal/mol. Therefore, even with the L294A mutation, further conformational rearrangement is required to attain an open state in the presence of regulatory ions. Similarly, wild type simulations initiated with Nwat = 25 or 42 remained stably hydrated with Nwat = 24 (σ = 5) or spontaneously became less hydrated until Nwat = 32 (σ = 7) at 4 ns, respectively (Fig 10A and 10B). Conversely, L294A simulations initiated with Nwat = 25 or 42 approached Nwat values of 37 and 45, respectively (σ = 5; Fig 10A and 10B). Although these extents of hydration predict ΔG‡ values between 9 and 16 kcal/mol, the standard deviations of five water molecules among our two-hundred repeat simulations indicates that ~33 simulations initiated with Nwat = 42 continued to sample values of Nwat greater than 50 after 4 ns of simulation, and thus are predicted to present a barrier to Mg2+ flux of only 5 kcal/mol. Therefore, in the absence of regulatory ions, the conformational changes that we observe in TmCorA may be sufficient to permit pore opening in the L294A mutant. We previously proposed that the unbinding of regulatory ions induces changes in the arrangement of protomers in the cytosolic domain of TmCorA which increases the probability of hydration of the MM due to allosteric coupling via the long kinked and tilted helices that wrap around the pore in an iris-like fashion [29]. The present work confirms these conformational changes. Specifically, removing all ten regulatory Mg2+ ions results in increased separation between the centers of mass of adjacent Mg2+ binding domains (S11A Fig) and decreases the radial tilt angle of the intracellular region of helix 7 (S11B Fig), which together define the width of the funnel-like cytosolic domain. In addition, both the pore radius and the lateral tilt of helix 7 in the MM region increase upon removing regulatory Mg2+ ions and further increase in the SSH runs (S11C and S11D Fig). The first link in this allosteric mechanism is demonstrated by a clear correlation between fluctuations in Mg2+ binding domain arrangement and changes in the radial tilt of the pore helices; although this correlation exists regardless of regulatory Mg2+ occupancy, these structural fluctuations reach larger amplitudes when the regulatory sites are empty (S12A and S12E and S12I Fig). At the other end of the allosteric transmission, pore diameter and hydration are strongly correlated in the SSH simulations (S12J and S13B Figs), consistent with the above analysis (Fig 5). However, the allosteric mechanism by which we previously proposed radial tilt to cause pore dilation [29] was not strictly maintained throughout our new set of simulations. Specifically, the statistical linkage between average helix tilts and pore diameter is weak when all the SSH runs are considered together (S12L and S12K Fig) and only three out of eight SSH simulations show a significant correlation between average radial and lateral tilts of the pore helices (as defined by a Pearson coefficient magnitude above 0. 4; see labels d, e, and h in S13C Fig). Furthermore, our previous description of the iris mechanism [29] involved a positive correlation between average radial and lateral tilts, whereas this correlation is negative in two of the SSH simulations (see labels d and e in S13C Fig), indicating that both increases and decreases in radial tilt of the cytosolic region of helix 7 can lead to pore dilation. A crystal structure of TmCorA recently obtained in the absence of divalent cations suggests that protomeric asymmetry involving bending of the intracellular domain with respect to the transmembrane region is a precondition for channel opening [14]. Accordingly, the analysis of our simulations shows that asymmetric stalk-helix bending is enhanced by the removal of regulatory ions (S14 Fig). However, there is no correlation between stalk-helix bending and the hydration of the hydrophobic gate (S15 and S16 Figs), suggesting that our 35-ns simulations sample early stages of pore wetting and dilation, and that longer timescales are required for full channel opening.
Like a previous simulation study of protein folding kinetics [35], this study exploited the fact that rare events can be observed from massively-repeated short time trajectories provided that the simulation time exceeds the minimum passage time through the transition state. Our simulations were long enough for wetting and dewetting of the channel to occur, and the number of repeats was large enough to observe hundreds of such transitions, allowing us to derive quantitative estimates of wetting and dewetting kinetics. Significantly, the MM became hydrated more often when regulatory ions were absent from the simulations (Fig 2). Experimentally, Payandeh et al. have shown that a D253F mutation in one of the two pentameric regulatory binding sites enhances bacterial growth in limiting Mg2+ conditions [17], presumably because it reduces the ability of these sites to bind Mg2+, preventing closure of TmCorA. More recently, Dalmas et al. have shown that a D253K mutation abolishes current regulation by the DCS [28]. In addition, our finding that regulatory ions affect hydration of the MM, but not the LC (Figs 2 and 7), suggests that it is the former that is allosterically controlled by regulatory ions in early wetting. However, structural models of the open state generated by Dalmas et al. based on spectroscopic data predict an increased pore diameter throughout the entire channel lumen [28], suggesting that longer-timescale pore opening events involve dilation of the LC. In this context, it remains unclear why a pore-widening mutation at the LC (L280A) compromises bacterial growth in limiting Mg2+ conditions [17]. The magnitude of the free energy barrier in the MM decreased as it became more hydrated (Fig 9). Consistent with this finding, mutations that shorten the sidechain of MM residue L294 enhance bacterial growth in limiting Mg2+ conditions [17,18] and L294A increases the rate of liposomal uptake of Mg2+ by 2. 5-fold in a fluorescence-based flux assay [17]. Simulations of the L294A mutant confirm its increased propensity for MM wetting (Fig 10) and suggest that the largest widening of the TmCorA pore observed in our simulations may be sufficient to permit Mg2+ flux in the L294A mutant. Our results suggest that these mutations enhance channel activity by increasing the water-accessible volume of the MM region, leading to increased hydration and thereby decreasing the barrier to Mg2+ flux. Hydrophobic gating has been discussed extensively in the literature [36–44]. The transient, reversible wetting of the MM revealed in our simulations (Fig 7) is consistent with the concept that water molecules near hydrophobic surfaces are close to a vapor-to-liquid phase transition [45]. Hydrophobic surfaces alter the phase behavior of water [46], especially when it is tightly confined [31], and wetting equilibria can be dramatically altered by subtle changes in hydrophobicity [32]. Specifically, hydrophobic surfaces enhance water density fluctuations in their vicinity [45,47], thereby increasing the probability of cavity formation [47]. When such a cavity spans the gap between hydrophobic surfaces, it can nucleate the formation of a vapor phase [31] via capillary drying [48]. Taken together, the above findings suggest that the hydrophobic gate of TmCorA is poised near the edge of a wetting/dewetting transition that can be pushed toward wetting by a single (pentameric) L→A or L→G mutation. Analogous roles have been ascribed to a single (hexameric) V→A mutation in the hydrophobic bottleneck of a calcium release-activated calcium channel [49] and a single I→G mutation that affects the wetting of hydrophobic surfaces during melittin dimerization [45]. In the present study, condensation of water in the MM was correlated to small-amplitude dilations of the pore-lining helices (Figs 4 and 5, S12B, S12F, and S12J Fig). Apolar surfaces only dry when they are sufficiently large [46] and capillary dehydration occurs more rapidly when the enclosing apolar surfaces are closer and longer [31]. At a modal length of 1. 9 nm (S2B and S2F Fig) and lumenal diameter on the order of 0. 4–0. 6 nm (Fig 5), with an average helical spread of 1. 54 nm (S12C and S12G Fig), the hydrophobic MM appears to be well constructed to employ capillary dehydration as a gating mechanism. We observed the following three distinct hydration states in the MM: First, the majority of our simulations remained in a dry state. Second, a wetted state of low hydration, characterized by an average NWat¯ = 15 (σ = 3), was transiently populated and occurred more often when regulatory ions were absent (Table 1 and Fig 7B and S2 Fig). Although the opening rate of TmCorA is unknown experimentally, the observed rate of wetting in our simulations, 8. 2×106 ± 7×105 s-1 (S1 Table), is 100–1,000 times faster than the opening rate of any known ion channel, including those involved in fast excitatory neurotransmission such as nicotinic acetylcholine receptors (1×104 to 8×104 s-1) [50,51] and AMPA glutamate receptors (7×104 to 8×104 s-1) [52,53]. Thus, together with the fact that this transient hydrated state is impermeable to Mg2+ (Figs 8 and 9), the high rate of wetting suggests that this state is an intermediate. Finally, the third state that we identified, the SSH state (Fig 4 and S4 Fig), differed from transiently-hydrated states in terms of the extent and duration of hydration and in that it was not observed in the presence of regulatory ions (Figs 3,6 and 7 and S6 Fig). The overall effect of removing regulatory ions was to increase the extent of hydration in the MM and to increase the stability of wetted states. We propose the following model of wetting and dewetting transitions in the MM, as depicted in Fig 11. In the presence of regulatory ions, the dry state is favored over the transiently hydrated state by 4. 2 ± 0. 1 kcal/mol (S2 Table). When regulatory ions are removed, the free energy of the transiently-hydrated state, and the barrier to hydration, are reduced by 0. 4 ± 0. 2 kcal/mol. Regulatory ions thus modulate the rate of wetting but not dewetting transitions respectively to and from the transiently-hydrated state. Furthermore, transitions to the SSH state only occurred in the absence of regulatory ions (Figs 3,6 and 7, and S6 Fig). This result, together with the linear increase in MM wetting with time (Fig 2A), indicates that equilibrium has not been reached in the absence of regulatory ions and that longer simulations would lead to increasingly divergent fractions of MM wetting based on regulatory ion occupancy. The steepness of the drop in ΔG‡ vs. the extent of hydration, which was observed both with and without regulatory ions (Figs 8 and 9), indicates switch-like behavior. We hypothesize that the removal of regulatory ions stabilizes the SSH state much more dramatically than it stabilizes transiently hydrated states and that the SSH state plays an important role in TmCorA' s gating mechanism. Extrapolation of the relationship between ΔG‡ and Nwat suggests that the channel may open to divalent cation flux when more than ~50 water molecules condense in the MM. Using large-scale sampling, we have characterized the mechanism and kinetics of wetting and dewetting transitions in a 1. 9-nm-long hydrophobic constriction of the magnesium channel CorA. Massively-repeated simulations reveal essential differences in the dynamic fluctuations and relaxation of this hydrophobic gate in response to the removal of regulatory ions more than 6 nm away. The analysis of rate constants for wetting and dewetting transitions and of free energy profiles of ionic permittivity indicates that the hydrated states observed in the present study are intermediate states rather than the open state of the channel. These results lead to a model of functional gating of CorA in which regulatory ions control pore hydration, which in turn controls the onset of ionic conduction. A key event underpinning this allosteric mechanism, wetting constitutes a requisite but not sufficient step towards channel opening, as it is the extent of hydration of the pore, not the presence of water per se, that determines the barrier to ion permeation. Taken together, these findings demonstrate how capillary wetting transitions help mediate robust, non-linear switching of ionic conduction and underline the relevance of hydrophobic gating to ion channel structure and function.
The simulation systems consisted of TmCorA [13], with and without 10 regulatory Mg2+ ions in the DCS, in a hydrated 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) bilayer. This system comprised 250,000 atoms. Simulations were performed with GROMACS [54]. The water model was TIP3P [55]. TmCorA was modeled by the OPLS-AA/L parameters [56,57]. DMPC was modeled by the Berger parameters [58] using the half-ε double-pairlist method [29]. The Mg2+ parameters were those of Åqvist [59]. For additional details of system setup and massively repeated simulation, and for a description of the algorithms that we used to quantify hydration, see S2 Text. To evaluate the free energy profile for the permeation of a divalent cation throughout the pore of TmCorA, we used umbrella sampling (US) MD simulations [33,34]. In these simulations, we inserted hexahydrated magnesium into the pore, instead of a naked ion, because interactions between magnesium and water molecules in its first hydration shell are very strong [60], suggesting that magnesium is hexahydrated throughout a large part of the conduction pathway. We conducted 261 2-ns simulations in which the axial position of the magnesium ion relative to the center of mass of the hydrophobic gate, z, was harmonically restrained to a specified value, zi0, for restraining potentials (umbrellas), i, distributed every 0. 05 nm in the range -8≤zi0≤5 nm. These 261 simulations constituted one set of US simulations and were used to generate one evaluation of the free energy profile. Altogether, we conducted 18 sets of US simulations, providing 18 independent evaluations of the free energy profile. In so doing, we computed free energy profiles for four different conformational basins of the TmCorA system, all drawn from our massively repeated sampling. The four basins were: those with the largest and smallest values of Nwat in the MM, obtained either in the presence or absence of regulatory ions. For each basin, we conducted 3 sets of US simulations, each using a different conformation extracted from a different simulation, except for the large-Nwat state in the absence of regulatory ions, for which 9 distinct conformations of TmCorA were used to compute 9 free energy profiles. For additional details, see S2 Text.
This study shows how rapid wetting/dewetting transitions in the pores of ion channels participate in the control of biological ion permeation. Ion channels catalyze ionic permeation across non-polar membranes via water-filled pores. However, non-polar stretches or hydrophobic bottlenecks are present in the pores of many ion channels. To clarify the relationship between channel regulation, pore hydration, and ion permeation, we examine how the slow relaxation of magnesium channel CorA from its closed state towards its open state modulates wetting of its hydrophobic bottleneck. Results provide a quantitative description of wetting and dewetting probabilities and kinetics and a quantitative relationship between the extent of pore hydration and the energetics of ion permeation, consistent with a mechanism of hydrophobic gating.
Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA