Real time dynamics of Gating-Related conformational changes in CorA
Abstract
CorA, a divalent-selective channel in the metal ion transport superfamily, is the major Mg2+-influx pathway in prokaryotes. CorA structures in closed (Mg2+-bound), and open (Mg2+-free) states, together with functional data showed that Mg2+-influx inhibits further Mg2+-uptake completing a regulatory feedback loop. While the closed state structure is a symmetric pentamer, the open state displayed unexpected asymmetric architectures. Using high-speed atomic force microscopy (HS-AFM), we explored the Mg2+-dependent gating transition of single CorA channels: HS-AFM movies during Mg2+-depletion experiments revealed the channel’s transition from a stable Mg2+-bound state over a highly mobile and dynamic state with fluctuating subunits to asymmetric structures with varying degree of protrusion heights from the membrane. Our data shows that at Mg2+-concentration below Kd, CorA adopts a dynamic (putatively open) state of multiple conformations that imply structural rearrangements through hinge-bending in TM1. We discuss how these structural dynamics define the functional behavior of this ligand-dependent channel.
Introduction
Magnesium (Mg2+) is a key divalent cation in biology. It regulates and maintains numerous, physiological functions such as nucleic acid stability, muscle contraction, heart rate and vascular tone, neurotransmitter release, and serves as cofactor in a myriad of enzymatic reactions (Altura, 1991; de Baaij et al., 2015; Jahnen-Dechent and Ketteler, 2012; Romani, 2013; Ryan, 1991). Most importantly, it coordinates with ATP, and is thus crucial for energy production in mitochondria (Altura, 1991; Pilchova et al., 2017; Romani, 2011; Swaminathan, 2003; Yamanaka et al., 2016). In order to store Mg2+ in the mitochondrial lumen it is imported via Mrs2 (Kolisek et al., 2003) and Alr2 (Liu et al., 2002) ion channels that are closely related to CorA, the main Mg2+-importer in bacteria (Guskov and Eshaghi, 2012; Hmiel et al., 1986; Knoop et al., 2005; Schweyen and Froschauer, 2007). Although these Mg2+-transport proteins do not show much sequence conservation, they all share two trans-membrane domains (TMDs) with the signature motif Glycine-Methionine-Asparagine (GMN) at the extracellular loop (Knoop et al., 2005; Schweyen and Froschauer, 2007).
The crystal structure of CorA from Thermotoga maritima in its Mg2+-bound closed state (Eshaghi et al., 2006; Lunin et al., 2006; Payandeh and Pai, 2006) revealed a 5-fold symmetric homo-pentamer forming an ~11 nm funnel-like structure, with large intracellular domains and Mg2+-ions bound between the subunits (Eshaghi et al., 2006; Guskov et al., 2012; Lerche et al., 2017; Lunin et al., 2006; Payandeh and Pai, 2006). Through CorA, Mg2+-homeostasis is achieved by a negative feedback mechanism in which Mg2+ acts as both, charge carrier and ligand (Dalmas et al., 2014b). The binding of Mg2+ at the cytoplasmic subunit interfaces leads to channel closing and thereby limits further Mg2+-influx (Dalmas et al., 2014a; Dalmas et al., 2014b; Palombo et al., 2012; Pfoh et al., 2012; Schindl et al., 2007).
Determining the structure of CorA in its Mg2+-unbound (Apo) form is considered a fundamental step towards understanding its gating mechanism. Yet, the first CorA crystal structure in the absence of divalent cations showed little or no changes compared to the fully Mg2+-bound closed channel structure, with only a slight kink between TMDs and intracellular domains (Pfoh et al., 2012). In contrast, EPR spectroscopy indicated much larger structural rearrangements within the protein upon Mg2+-dissociation, suggestive of a more dramatic closed-to-open channel transition (Dalmas et al., 2010; Dalmas et al., 2014b). These results were confirmed by cryo-electron microscopy (cryo-EM), reporting surprising asymmetric rearrangements of the individual subunits within the homo-pentamer in Mg2+-free condition (Cleverley et al., 2015; Matthies et al., 2016). Based on ~7 Å cryo-EM structures (Matthies et al., 2016), a model was proposed in which Mg2+-unbinding leads to an increase in the inter-subunit conformational flexibility of CorA, thereby resulting in at least two asymmetric structures with, presumably, open gates (Matthies et al., 2016).
Despite progresses in the structural characterization of the Mg2+-free CorA structure(s), the transition from the closed (Mg2+-bound) state to the conductive unliganded conformation is still unclear. At 7.1 Å, the putative open CorA asymmetric cryo-EM structures were solved only from a relatively low number of particles (26,271 and 27,416 of 173,653 particles for open-I and open-II states, respectively) (Matthies et al., 2016), and due to its intrinsic averaging for structure determination, the cryo-EM data cannot inform about time-dependent behavior. We posit that these structures could be sampled from three different dynamic processes: they could represent two different open state conformations with significant state dwell-times (structurally resolved by cryo-EM); they could be intermediates along a conformational trajectory; or finally, they might be mere snapshots of a highly dynamic, fluctuating molecule without a long-lasting and well-defined high-resolution structure.
In order to elucidate the sequence of events underlying gating-related CorA conformational changes, we used high-speed atomic force microscopy (HS-AFM) (Ando et al., 2001; Ando et al., 2014) to capture CorA structural rearrangements during Mg2+-depletion experiments. HS-AFM is unique in its ability to concomitantly characterize molecular structures and dynamics under native-like conditions. In agreement with electrophysiological, spectroscopic and biophysical data (Dalmas et al., 2014b), we found that during Mg2+-depletion experiments individual CorA molecules lost pentameric symmetry and became highly dynamic at Mg2+-concentrations below ~2 mM. Under these conditions, the symmetric state is reversibly adopted by means of dynamic fluctuations among the individual subunits, indicative of spontaneous Mg2+-rebinding. However, in the absence of Mg2+, CorA transitions to an asymmetric state. Thus, the conformational energy landscape likely comprises two deep energy minima represented by the fully symmetric (Mg2+-bound) and the asymmetric (Mg2+-depleted) states interspersed by a wide plateau of conformational fluctuations of a highly flexible molecule, which likely represents the conductive state.
Results
Lipid composition determines reconstitution density and morphology
In recent years, HS-AFM has demonstrated its power to observe molecular mechanisms and structural dynamics of single molecules under physiological conditions with extraordinary detail (Chiaruttini et al., 2015; Kodera et al., 2010; Marchesi et al., 2018; Preiner et al., 2014; Rangl et al., 2016; Ruan et al., 2018; Ruan et al., 2017). Here, we have used HS-AFM to track the conformational changes of CorA as a result of Mg2+-concentration changes in real-time. For this, wild-type (wt) CorA was reconstituted at low lipid-to-protein ratios (LPRs between 0.2 and 0.4) in the presence of Mg2+. The resulting proteo-liposomes were adsorbed onto freshly cleaved mica under saturating Mg2+ (10 mM) condition and imaged at an acquisition rate of 1–2 frames s−1 with a resolution of 0.5 nm pixel−1. CorA-containing vesicles were generated from protein reconstitutions in POPC/POPG (3:1, w:w). When imaged by HS-AFM, vesicles spread on the mica support resulting in large membrane patches with densely packed CorA (Video 1). These patches mostly exposed the periplasmic face. Smaller crowded areas of molecules expose the intracellular side (Figure 1a, Video 1).
Cross-section analysis of the CorA membrane patches showed overall heights of ~12 nm (Figure 1b). High-resolution movies of the periplasmic domains revealed protrusions of 0.9 ± 0.4 nm in height and 3.4 ± 0.4 nm in diameter, and resolved a central indentation where the channel pore is located (Figure 1c). Despite the fast diffusion of CorA molecules within the clusters exposing the intracellular side, some high-resolution snapshots allowed analysis of the surface structure. The intracellular face protruded ~7 nm from the membrane (Figure 1—figure supplement 1a) with top-ring and outer diameters of ~6 nm and ~10 nm, respectively (Figure 1—figure supplement 1b,c; Video 2). These measurements are in good agreement with the molecular dimensions of CorA (Eshaghi et al., 2006; Lunin et al., 2006; Matthies et al., 2016).
However, it became clear that in order to study Mg2+-dependent conformational dynamics, more stably packed molecules were required (Müller et al., 1999; Ramadurai et al., 2010). We found that reconstituting the protein in DOPC/DOPE/DOPS (4:5:1, w:w:w) and adsorbed and imaged at lower pH 6.0 resulted in widespread surface coverage of CorA-crowded membranes with only slowly moving molecules. Moreover, under these conditions densely packed CorA patches were stacked and exposed the intracellular face of the channel (Figure 1d,e; Video 3), thus providing an excellent experimental platform for studying CorA Mg2+-dependent conformational changes. High-resolution HS-AFM topographs displayed a ‘flower-shaped’ surface structure corresponding to the CorA intracellular face. This view allowed to resolve individual subunits in the pentamer with top-ring diameter of 5.0 ± 0.9 nm and a center-to-center distance of the molecules, that is outer diameter, of 10.9 ± 2.1 nm (Figure 1f).
Mg2+-depletion induces large conformational changes of the intracellular face
We then monitored the structural changes of the channel upon Mg2+-depletion in real time at higher magnification, that is at scan sizes < 200 nm and pixel sampling of 0.5 nm pixel−1. First, we studied the periplasmic face of CorA under saturating Mg2+ (10 mM), in which all channels are expected to be in the closed conformation. Subsequently, membranes were imaged at reduced Mg2+-concentration near the reported apparent Kd for Mg2+ (~2 mM), and finally at 0 mM Mg2+ to focus on the dynamically open Apo form.
Mg2+-depletion was achieved by live injection of EDTA, while continuously monitoring CorA. Using the closed state as reference, no structural changes of the periplasmic face of CorA were observed at 10 mM Mg2+, at 2 mM Mg2+ (after ~10 min) and in absence of Mg2+ (after ~12 min) (Figure 2—figure supplement 1). At the present resolution (0.5 nm pixel−1) any putative changes in the conformation of the periplasmic face of the channel are beyond the current resolution limit of HS-AFM (although we cannot rule out that the structures of the small periplasmic loops might not change upon Mg2+-depletion). This is consistent with the ~7 Å cryo-EM structures, which show little or no Mg2+-dependent changes on the periplasmic face (Matthies et al., 2016; Pfoh et al., 2012).
In stark contrast, previous EPR (Dalmas et al., 2010; Dalmas et al., 2014b) and cryo-EM (Matthies et al., 2016) analyses have revealed dramatic conformational rearrangements of the intracellular domains (resulting in a loss of symmetry) in the nominal absence of Mg2+ (conditions favoring the functionally open state) (Figure 2a). First, to measure the CorA Mg2+-affinity (Pilchova et al., 2017) in our experimental set-up, we used a microfluidic system connected to a constant pressure and flow pump (Miyagi et al., 2016), with which we slowly exchanged the complete 10 mM Mg2+ measuring solution to a Mg2+-free buffer (containing additional 2 mM EDTA). Analysis of the distribution of symmetric and dynamically asymmetric CorA particles pointed to a Kd of ~2 mM Mg2+, which is in good agreement with the reported affinity (Figure 2—figure supplement 2). Next, we monitored the structural and dynamical transition of CorA upon Mg2+-depletion by pipetting defined amounts of EDTA into the measurement fluid cell to achieve the following equilibrium concentrations: 10 mM Mg2+, ~2 mM Mg2+ (~Kd), and 0 mM Mg2+ (putative).
Consistent with the Mg2+-liganded cryo-EM structure (Figure 2a, left), at saturating 10 mM Mg2+ condition CorA channels revealed a stable, flower-like 5-fold symmetric structure (Figure 2b, left). However, once Mg2+-concentrations dropped to ~2 mM and below, individual channels start to fluctuate between various structural states, occasionally assuming increased height and, after ~20 min in absence of Mg2+, adopting an ill-defined conformation with significantly increased protrusion height (Figure 2b, right). The Mg2+-dependent loss of symmetry and increased protrusion height are also consistent with expected structural features of the Mg2+-free (open) CorA structures, where individual subunits move towards the former 5-fold axis and thus stand taller (Matthies et al., 2016) (Figure 2a, bottom). The large, ~1.5 nm, protrusion height difference (ΔHeight) of such a conformation change presents a useful signature to detect and follow the conformational dynamics of individual channels in the membrane. However, these two ΔHeight-states must not be mistaken with functional or even structural states. They are merely a way to topographically discriminate between the closed 5-fold symmetric state and any other state where single subunits stand up, that is move towards the 5-fold axis and thus appear higher.
Height section kymographs of individual molecules (Figure 3a, top) over extended imaging periods (Figure 3a, middle) allowed the subsequent computational detection of single molecule dynamics based on the fluctuations in height differences (ΔHeight) as a function of time (Figure 3a, bottom): While a ΔHeight of ~1 nm represents a symmetric CorA structure with a pore in the center, larger ΔHeight values represent molecules with increased protrusion height. ΔHeight values < 1 nm were found for molecules of low height (putative closed state) where the central pore could not be resolved. Hence, the ΔHeight/time traces were idealized to represent two ΔHeight-states using the Step Transition and State Identification (STaSI) algorithm developed for single molecule experiments (Shuang et al., 2014), which we successfully adapted for the analysis of HS-AFM ΔHeight/time traces (Heath and Scheuring, 2018) (Figure 3a, red line).
We monitored molecular transitions from the Mg2+-bound to Mg2+-free states in time-lapse experiments. HS-AFM images of CorA membrane patches during these extended experiments reproducibly revealed the conformational changes of individual molecules over time (Figure 3b, Video 4, Video 5). While CorA maintained the 5-fold symmetric state at saturating Mg2+ condition (Figure 3b, green time stamps), the channels started switching dynamically between (at least) two conformational states at concentrations below ~2 mM Mg2+ (Figure 3b, blue and early red time stamps), and began populating strongly protruding asymmetric states as Mg2+-depletion progressed (Figure 3b, red time stamps). A cumulative height-based state-assignment of CorA as a function of Mg2+-depletion demonstrated that the number of channels in the putatively open state(s) gradually increased (Figure 3b, bottom). In contrast, the number of transitions (i.e. switching back and forth between the high (open) state(s) and the state of lower height) cumulated earlier during the titration experiment (Figure 3c). Thus, we are looking at a transition from the symmetric low (closed) to the asymmetric high (open) structures and a wide range of fluctuations between states at intermediate Mg2+-concentrations, or within the first minutes of complete Mg2+-depletion.
CorA state transition revealed a highly dynamic intermediate
In addition to the gradual state conversion during Mg2+-depletion, the number of transition events varied over time. At saturating Mg2+, all molecules occupied a stable conformation and the state interconversion activity is virtually zero (Figure 3b, green). About 5–8 min into the complete Mg2+-depletion experiment (at 0 mM Mg2+) the dynamics of the channels reach a maximum (Figure 3c), followed by another more stable, strongly protruding and asymmetric state, putatively representing an open conformation (Figure 3b,c, red). The time windows of these observations correlate well with electrophysiology experiments, in which the CorA-driven Mg2+-currents in Xenopus leavis oocytes decay within 15–20 min (Dalmas et al., 2010; Dalmas et al., 2014b). In addition to the number of transitions, the same trend is represented when analyzing the dwell-times of the ΔHeight-states in the nominal absence of Mg2+ (Figure 3c): Within the first 5 min in 0 mM Mg2+ the average time spent in the closed (low height) state is almost twice (~60 s) the averaged dwell-time in one of the putatively open (elevated height) conformations. After ~10 min in Mg2+-free condition, the situation is reversed, and more molecules show long dwell-times in the elevated state (Figure 3c, turquois pentagons and red diamonds). Thus, in this experiment, the CorA gating pathway initially favors the stable Mg2+-bound, low-height state by ~- 0.7kBT while elevated Mg2+-free (open) conformation(s) are favored when Mg2+ is depleted from all binding sites. In between, the channel is in a highly dynamic regime, probably reflecting Mg2+-unbinding and -rebinding events (Figure 3, Videos 4 and 5). We must highlight, however, that a total depletion of Mg2+ probably never occurs under physiological conditions, thus the fully Mg2+-depleted asymmetric stable states observed by cryo-EM and likely adopted here at the end of the depletion experiments might not be visited in the cell. More likely, the 5-fold symmetric closed state interchanges with a highly fluctuating molecule where single subunits dissociate from the quaternary structure of the cytoplasmic ensemble, opening the channel.
CorA fluctuates between several conformations in the open state
Following the general observation that the stable closed and open conformations are interconnected by a regime in which the channels are highly dynamic, we pursued a detailed structural examination of CorA in this intermediate dynamic stage at low Mg2+-concentrations. High-resolution HS-AFM image sequences of individual CorA channels revealed that CorA undergoes conformational rearrangements of the entire cytosolic Mg2+-sensor domain at rates beyond the bandwidth of the current measurements (<550 ms, the imaging rate of our videos). We also acquired movies at 250 ms frame acquisition and found that individual sequential frames displayed different apparently unrelated conformations, an indication that the conformational fluctuations are even faster (Figure 4d).
Movie snapshots were classified into four conformational classes (Figure 4, Video 6): First, the fully-liganded 5-fold symmetric closed conformation (Eshaghi et al., 2006; Lunin et al., 2006; Matthies et al., 2016; Payandeh and Pai, 2006) (Figure 4a, top row, Figure 4b, 3.85 s, 7.15 s, 8,8 s and 14.3 s and Figure 4c, 1.1 s, 5.5 s and 6.05 s, Figures 4d, 0s, 0.75 s and 2.25 s). Second, a structure of reduced diameter with three asymmetrically distributed protrusions of increased height (Figure 4a, second row, Figure 4b, 9.35 s, Figure 4c, 0.55 s, 4.95 s, Figure 4d, 0.5 s), likely corresponding to the cryo-EM open-I state (Matthies et al., 2016). Third, a round- or star-shaped conformation with an elevated center resembling the cryo-EM open-II state where one subunit is displaced towards the channel axis (Matthies et al., 2016) (Figure 4a, third row, Figure 4b, 0s, 6.05 s, 7.7 s, 17.05 s and Figures 4c, 0s, 0.55 s, 2.75 s). However, given that our tentative assignment to open-I and open-II states only covers a fraction of all molecular observations, a fourth category was invoked to include the remaining range of structural variations, which we named open-+ states. This category includes asymmetric molecules of all kinds, even those without increased height (Figure 4a, fourth row, Figure 4b, 16.5 s, Figure 4d, 1.5 s, 1.75 s). It is clear that a key challenge in defining the conformational landscape of CorA will be the unbiased classification of discrete states in this highly flexible molecule, made evident by following a single CorA channel frame by frame (Figure 4c,d). Unfortunately, all our computational classification attempts failed likely due to the presence of protein-protein contacts with neighboring molecules and limited resolution due to the high mobility.
To work out a more detailed picture of the CorA-transitions we visually assigned these four different types of states to high-resolution HS-AFM movie sections recorded at different Mg2+-concentrations. In saturating Mg2+-concentrations (Figure 5a, green bars, Video 7), ~91% of the time molecules adopt the flower-shaped, 5-fold symmetric state, while the other states are barely populated. After Mg2+-reduction to 3 mM (Figure 5a, blue bars) and full Mg2+-depletion (Figure 5a, red bars, Video 8), the probability of finding CorA in any of the asymmetric states is reversed (Figures 5a, 3 classes on the right). Importantly, in 3 mM Mg2+ ~50% of the molecules pool in the open-+ (others) class, representative of a highly mobile and flexible molecule that switches fast between numerous sub-structures that could not be classified to any of the conformations described by cryo-EM. After full depletion, re-addition of Mg2+ could partially recover ~35% the symmetric closed-state structure (Figure 5a, yellow bars). Analysis of the probability of transitions from one structural state to the others (Figure 5b) demonstrated that under saturating Mg2+-concentrations the symmetric closed state is, as expected, stable. Once CorA starts adopting any of the asymmetric conformational states in Mg2+concentrations ~3 mM and lower, the probability to switch back to the symmetric state is only ~14%. This probability is further reduced to ~7% at 0 mM Mg2+.
Interestingly, comparing 3 mM and 0 mM Mg2+-conditions gave a hint of the different Mg2+-loads between the various asymmetric states: while the number of transitions between elongated states (Figure 5b, red outline: open-I and open-+) is virtually identical in these two conditions, the transition into a state with elevated height is strongly favored at 0 mM Mg2+ (Figure 5b, dashed blue outline: open-II and open-I). We suggest that molecules with increased height where one subunit moves up close to the channel axis represent Mg2+-free states.
Discussion
On the basis of real-time HS-AFM imaging, we find that CorA Mg2+-dependent gating can be best described in three phases: Above the apparent Mg2+-affinity of the intracellular Mg2+-sensor sites (~2 mM Mg2+), the channel adopts a stable 5-fold symmetric state reminiscent of the high-resolution structures (phase 1). After several minutes of exposure to a Mg2+-free solution, CorA transitions into a set of asymmetric architectures, which are characterized by elevated cytoplasmic domains (phase 3). These structures are adopted through major conformational changes implicating hinge-bending of TM1 and reorientation of the cytoplasmic domains. This picture is wholly consistent with data derived from functional (electrophysiology), biochemical (crosslinking), biophysical (EPR, Fluorescence) and structural (cryo-EM) approaches. However, these results contrast with recent crystallographic data of CorA Mg2+ binding site mutants that reported limited or no structural changes compared to the closed wt channel, thereby challenging the hypothesis that Mg2+-depletion leads to large conformational changes and channel opening (Kowatz and Maguire, 2019). We argue, that the constrains of the crystal lattice on CorA limits the range of conformational changes that can be observed in X-ray structures. In contrast, HS-AFM imaging under physiological conditions shown here is more consistent with the single particle cryo-EM data reported in Cleverley et al. (2015); Matthies et al. (2016) and the EPR spectroscopic data (Dalmas et al., 2010; Dalmas et al., 2014b), that is when the Mg2+ is investigated under lattice-free conditions. These facts must be considered a critical factor in explaining these differences. Among multiple conformational states within phase 3, we identified CorA structures that resembled the low-resolution cryo-EM open-I and open-II structures. These phases 1 and 3 are characterized by structural stability or at least limitations by the structural freedom of CorA. Here, more work is required to identify and clearly distinguish between different sub-structures in the Mg2+-free elevated state of CorA. In contrast, at intermediate Mg2+-concentrations, and/or after short incubation times at 0 mM Mg2+, the channel reveals a highly mobile and fluctuating state (phase 2). Active subunits change and the entire molecule displays high structural variability, likely adopting many more conformations than the three states so far identified by cryo-EM. Indeed, while we can assign the fully 5-fold symmetric state with certainty, our assignment of molecules to state open-I and open-II must be considered as putative. Perhaps more importantly, at intermediate Mg2+-concentrations ~ 50% of the molecules could not be assigned and were pooled in a class of undefined quaternary structures (open-+). In this context, the details of the cryo-EM study become relevant: open-I and open-II only pooled 26,271 and 27,416 particles, respectively, of the total 173,653 particles, and thus a significant portion of the particles remained outside of these classes (Matthies et al., 2016). We note, however, that this initial assignment is ultimately constrained by the resolution limitations of both cryo-EM structural assignments and the present application of HS-AFM. We expect that given a larger number of particles to process (in cryo-EM) or improvements in the spatial and/or temporal resolution of HS-AFM, additional states are likely to emerge. We propose that the various CorA structural states and the associated subunit movements observed by HS-AFM may reflect variations of number and location of Mg2+-ions bound to the cytosolic domain (Figure 5c).
Among the main structural features of CorA is its long (~11 nm) first transmembrane helix (TM1), ranging all the way from the periplasmic face to the surface of the cytoplasmic Mg2+-sensor. TM1 is actually the only connection between the TM region and the cytoplasmic domain (Figure 5—figure supplement 1). Thus, it seems plausible that the large fluctuations of the cytoplasmic domains might be directly translated by TM1 into fluctuations within the channel pore. We suggest this mechanism as the basis for ion conductance through a particularly long (in comparison to other channels) ~55 Å pore, in which moving TM1s allow the asynchronous progression of ions through the pore. Considering the typical persistence length of an α-helix of ~100 nm (Choe and Sun, 2005), allowed us to estimate how much an 11 nm long thermally activated helix fluctuates (Figure 5 – Information Supplement 1). We find that at 1kBT TM1 could bend at its cytoplasmic end ~2 nm out of axis, in good agreement with observed bending of this helix in the open-I and open-II cryo-EM structures and the HS-AFM movies of subunit displacement. Thus, we propose that in the open state the pore is by far less narrow than assumed based on static structures. Indeed, constantly fluctuating TM1s would provide a considerable channel diameter towards the cytoplasmic side. This model would predict that amino acids close to the periplasmic side of the pore play the key role for gating. In agreement with this hypothesis, we find that the Mg2+-channel signature motif (GMN) is located right at the periplasmic end of the channel, whereas amino acids towards the cytoplasmic face of TM1 are less conserved (Figure 5—figure supplement 1).
The idea that there are multiple conductive conformations of CorA, or, as a paradigm shift, the conductive conformation of CorA is a fluctuating molecule, is favored by both, the subunit movement analysis and state transition analysis. In the long-term absence of Mg2+, CorA molecules adopt a rather stabilized asymmetric state resembling the putatively open state cryo-EM structures. However, a living cell under physiological conditions is unlikely to ever reach Mg2+-concentrations below ~1 mM. Thus we hypothesize that in cellula the physiologically relevant open state is a collection of fluctuating conformational states. Thus, the conformational energy landscape of CorA gating would consist of two deep energy minima representing the stable closed (symmetric) and putatively open (asymmetric) conformations connected by a wide plateau in which a variety of open conformations are adopted through permanent molecular fluctuations. The compelling utility of HS-AFM for the study of macromolecular dynamics at high spatio-temporal resolution and in native-like conditions is demonstrated, elucidating a process so far inaccessible to other structural or biophysical techniques. The present data set and our proposed mechanistic interpretation of the conformational dynamics of Mg2+-dependent CorA gating sets the stage for an unprecedented understanding of CorA as a ‘reverse’ polarity ligand gated ion channel with an unique gating mechanism.
Materials and methods
Protein purification
Request a detailed protocolCorA from T. maritima was expressed and purified as previously described (Dalmas et al., 2010). Briefly, the CorA-Pet15b vector was used to transform, then express CorA in E. coli BL21 DE3. After cell harvesting and disruption, membranes were collected by ultracentrifugation and gently solubilized. The sample was then cleared by ultracentrifugation and purified using cobalt high affinity chromatography column (Clontech Laboratories). The concentrated protein sample (AMICON 100 kDa cutoff membrane filters, EMD Millipore) was homogenized by gel filtration (Superdex 200 10/300 GL column, GE Healthcare Bio Sciences) and equilibrated in 50 mM HEPES, pH 7.3, 200 mM NaCl, 20 mM MgCl2, and 1 mM DDM).
Protein reconstitution
Request a detailed protocolFor CorA reconstitution into liposomes, the protein was gently mixed with freshly prepared solubilized POPC-POPG (3:1) or DOPC-DOPE-DOPS (4:5:1) lipids (Avanti Polar Lipids) at low lipid to protein ratios (LPR) between 0.2–0.4 (w:w) at a total protein concentration of 1 mg/mL. After 4 hr equilibration, detergent was removed by addition of biobeads overnight.
Sample preparation for HS-AFM
Request a detailed protocolA 1.5 mm diameter muscovite mica sheet was glued on a HS-AFM glass rod sample support and mounted on a HS-AFM scanner. Reconstituted CorA membranes were adsorbed on freshly cleaved mica for ~5 min. Subsequently, the sample was rinsed with imaging buffer (50 mM MES, pH6.0, 200 mM NaCl) containing 10 mM Mg2+.
HS-AFM
Request a detailed protocolAll experiments were performed using HS-AFM (Ando et al., 2001) (SS-NEX, Research Institute of Biomolecule Metrology Co.) operated in amplitude modulation mode, using ultra-short cantilevers (8 µm) with a nominal spring constant of ~0.15 N/m and a resonance frequency of ~600 kHz in liquid (USC, NanoWorld). Videos of CorA membranes were recorded with imaging rates of ~1–2 frames s−1 and at a resolution of 0.5 nm pixel−1. The energy input by the AFM tip (estimated to ~1.5 kBT, considering a 90% imaging amplitude of a 1 nm free amplitude)(Miyagi et al., 2016) was minimized by continuously adapting the drive and setpoint amplitude and optimizing the feedback parameters.
Structural titration experiments
Request a detailed protocolMonitoring the transition of CorA from Mg2+-saturated to low/no Mg2+-conditions was achieved by depleting Mg2+ by either adding EDTA into the measuring solution or alternatively by buffer solution exchange to Mg2+-free buffer using an integrated constant-pressure and constant-flow pump system (Miyagi et al., 2016). Experiments were performed with CorA membranes prepared from three different purifications and ~10 different reconstitutions. In total, about 50 transition experiments on different membranes patches were performed on different days (over ~30 days) using two different, but similar HS-AFM systems, and ~10 USC cantilevers. All recordings showed a clear structural transition from stable 5-fold symmetrical proteins to very dynamic molecules with increased height. About 20% of the recorded CorA transition movies were considered for further analysis.
Data analysis
Request a detailed protocolHS-AFM images were first-order flattened and contrast adjusted using laboratory-made routines in Igor Pro software (WaveMetrics). Videos were then drift corrected with respect to the membrane patch, or aligned on individual CorA molecules using an in-house developed analysis software plug-in for ImageJ (Fechner et al., 2009; Husain et al., 2012). Dimensions of the CorA were calculated by height histogram analysis (n = 20,500 height values) and cross section analysis (n = 25) for each condition. Estimation of dwell-times was based on CorA protrusion height: ΔHeight/time traces were generated by subtracting the minimum pixel value from the maximum pixel value in the 5 × 5 nm center area of the molecule. For the molecular height transition detection the ΔHeight traces were analyzed by a Step Transition and State Identification (STaSI) method in a MatLAB (MathWorks) routine (Heath and Scheuring, 2018; Shuang et al., 2014). StaSI indicated a minimum description length (MDL) for fitting the data with two ΔHeight states (where the increased height states comprises all conformational states with elevated topography). For in depth analysis, 25 molecules of two individual experiments were tracked over 900 and 1600 frames at different Mg2+ conditions, resulting in ~30,500 analyzed frames. Molecules in different height states and the corresponding dwell-times were binned either over 20 events along the time axis. Note, the transition analysis only discerned between states where subunits fluctuate between different height levels and did not classify between sub-states that exposed equivalent height levels to the putative closed (symmetric) state or the putatively open, activated (asymmetric and elevated) state. Such, the height/time traces did also not discriminate which of the open states was assumed. Notably however, our experiments showed that structures exposing increased protrusion heights are likely representatives of Mg2+-depleted molecules and thus detection of elevated height is a valuable fingerprint for the state transition. For conformational transition analysis membrane patches of ~20 CorA channels were imaged at Mg2+ concentrations of 10 mM, 3 mM, 0 mM, and after re-addition of 25 mM. Extracted CorA molecules from high-resolution movie sections (1,000–2000 molecular representations extracted from about 20 molecules for each Mg2+ condition) were manually assigned to different structural states of (C) symmetric flower-shaped, (O-1) asymmetric elevated, (O-2) dome-shaped, elevated and (O-+) all others. The significance of the occurrence of each conformation under the four tested Mg2+ concentrations was tested using a two-tailed students test.
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files.
References
-
Basic biochemistry and physiology of magnesium: a brief reviewMagnesium and Trace Elements 10:167–171.
-
Filming biomolecular processes by high-speed atomic force microscopyChemical Reviews 114:3120–3188.https://doi.org/10.1021/cr4003837
-
The elasticity of alpha-helicesThe Journal of Chemical Physics 122:244912.https://doi.org/10.1063/1.1940048
-
The Cryo-EM structure of the CorA channel from Methanocaldococcus jannaschii in low magnesium conditionsBiochimica Et Biophysica Acta (BBA) - Biomembranes 1848:2206–2215.https://doi.org/10.1016/j.bbamem.2015.06.002
-
Molecular mechanism of Mg2+-dependent gating in CorANature Communications 5:3590.https://doi.org/10.1038/ncomms4590
-
Magnesium in man: implications for health and diseasePhysiological Reviews 95:1–46.https://doi.org/10.1152/physrev.00012.2014
-
The mechanisms of Mg2+ and Co2+ transport by the CorA family of divalent cation transportersCurrent Topics in Membranes 69:393–414.https://doi.org/10.1016/B978-0-12-394390-3.00014-8
-
Software for drift compensation, particle tracking and particle analysis of high-speed atomic force microscopy image seriesJournal of Molecular Recognition 25:292–298.https://doi.org/10.1002/jmr.2187
-
Transport of magnesium and other divalent cations: evolution of the 2-TM-GxN proteins in the MIT superfamilyMolecular Genetics and Genomics 274:205–216.https://doi.org/10.1007/s00438-005-0011-x
-
Loss of cytosolic Mg2+ binding sites in the Thermotoga maritima CorA Mg2+ channel is not sufficient for channel openingBiochimica Et Biophysica Acta (BBA) - General Subjects 1863:25–30.https://doi.org/10.1016/j.bbagen.2018.09.001
-
The periplasmic loop provides stability to the open state of the CorA magnesium channelJournal of Biological Chemistry 287:27547–27555.https://doi.org/10.1074/jbc.M112.371484
-
A structural basis for Mg2+ homeostasis and the CorA translocation cycleThe EMBO Journal 25:3762–3773.https://doi.org/10.1038/sj.emboj.7601269
-
The involvement of Mg2+ in regulation of cellular and mitochondrial functionsOxidative Medicine and Cellular Longevity 2017:1–8.https://doi.org/10.1155/2017/6797460
-
IgGs are made for walking on bacterial and viral surfacesNature Communications 5:4394.https://doi.org/10.1038/ncomms5394
-
Cellular magnesium homeostasisArchives of Biochemistry and Biophysics 512:1–23.https://doi.org/10.1016/j.abb.2011.05.010
-
Magnesium in health and diseaseMetal Ions in Life Sciences 13:49–79.https://doi.org/10.1007/978-94-007-7500-8_3
-
The role of magnesium in clinical biochemistry: an overviewAnnals of Clinical Biochemistry: An International Journal of Biochemistry and Laboratory Medicine 28:19–26.https://doi.org/10.1177/000456329102800103
-
Mrs2p forms a high conductance Mg2+ selective channel in mitochondriaBiophysical Journal 93:3872–3883.https://doi.org/10.1529/biophysj.107.112318
-
New Perspectives in Magnesium Research: Nutrition and HealthCorA-Mrs2-Alr1 Superfamily of Mg2+ Channel Proteins, New Perspectives in Magnesium Research: Nutrition and Health, Springer.
-
Fast step transition and state identification (STaSI) for discrete Single-Molecule data analysisThe Journal of Physical Chemistry Letters 5:3157–3161.https://doi.org/10.1021/jz501435p
-
Magnesium metabolism and its disordersThe Clinical Biochemist. Reviews 24:47–66.
Article and author information
Author details
Funding
National Institutes of Health (R01GM120561)
- Eduardo Perozo
National Institutes of Health (DP1AT010874)
- Simon Scheuring
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank George Heath for providing the Matlab scripts for the STaSI trace analysis, and Yi-Chih Lin for valuable comments on data analysis.
Copyright
© 2019, Rangl et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.
-
- Structural Biology and Molecular Biophysics
The functional effects of an RNA can arise from complex three-dimensional folds known as tertiary structures. However, predicting the tertiary structure of an RNA and whether an RNA adopts distinct tertiary conformations remains challenging. To address this, we developed BASH MaP, a single-molecule dimethyl sulfate (DMS) footprinting method and DAGGER, a computational pipeline, to identify alternative tertiary structures adopted by different molecules of RNA. BASH MaP utilizes potassium borohydride to reveal the chemical accessibility of the N7 position of guanosine, a key mediator of tertiary structures. We used BASH MaP to identify diverse conformational states and dynamics of RNA G-quadruplexes, an important RNA tertiary motif, in vitro and in cells. BASH MaP and DAGGER analysis of the fluorogenic aptamer Spinach reveals that it adopts alternative tertiary conformations which determine its fluorescence states. BASH MaP thus provides an approach for structural analysis of RNA by revealing previously undetectable tertiary structures.