High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter
Abstract
Transporters of the Nramp (Natural resistance-associated macrophage protein) family import divalent transition metal ions into cells of most organisms. By supporting metal homeostasis, Nramps prevent diseases and disorders related to metal insufficiency or overload. Previous studies revealed that Nramps take on a LeuT fold and identified the metal-binding site. We present high-resolution structures of Deinococcus radiodurans (Dra)Nramp in three stable conformations of the transport cycle revealing that global conformational changes are supported by distinct coordination geometries of its physiological substrate, Mn2+, across conformations, and by conserved networks of polar residues lining the inner and outer gates. In addition, a high-resolution Cd2+-bound structure highlights differences in how Cd2+ and Mn2+ are coordinated by DraNramp. Complementary metal binding studies using isothermal titration calorimetry with a series of mutated DraNramp proteins indicate that the thermodynamic landscape for binding and transporting physiological metals like Mn2+ is different and more robust to perturbation than for transporting the toxic Cd2+ metal. Overall, the affinity measurements and high-resolution structural information on metal substrate binding provide a foundation for understanding the substrate selectivity of essential metal ion transporters like Nramps.
Editor's evaluation
This manuscript provides fundamental new insight into protein conformational transitions underlying the transport mechanism of Nramps, an important and widespread transporter family that facilitates the uptake and movement of essential transition metals. Eight new crystallographic structures of the prokaryotic homolog DraNramp in a variety of ligand-bound and conformational states, along with companion molecular dynamics simulations and metal binding and transport assays, provide compelling evidence supporting most of the conclusions. These findings will be of broad interest to scientists studying transport mechanisms and ligand recognition.
https://doi.org/10.7554/eLife.84006.sa0Introduction
Transition metal ions like Mn2+ and Fe2+ are essential for various metabolic processes in all living cells and are usually required in low intracellular concentrations for optimal activity (Andrews, 2002; Bozzi and Gaudet, 2021). Excess or deficiency of transition metal ions leads to diseases (Bleackley and Macgillivray, 2011; Nies and Grass, 2009). For example, Fe2+ deficiency causes anemia and neurodegenerative diseases, whereas Fe2+ overload increases the risk of cancer by generating toxic reactive oxygen species (ROS) and mutations (Ekiz et al., 2005; Jung et al., 2015). Mn2+ overload in the brain is linked to neurological disorders and deficiency causes metabolic defects and impairs growth (Budinger et al., 2021; Pittman, 2005). Other transitions metals, like Cd2+ and Hg2+, are toxic and their accumulation affects health by disrupting the physiological levels of essential metals or displacing them in enzyme active sites, thus inhibiting the proteins, or changing their activity (Andrews, 2002; Lin et al., 2009). Cells and organisms have evolved strategies to maintain metal ion homeostasis via highly regulated transport and storage processes (Bleackley and Macgillivray, 2011; Cellier and Gros, 2004; Nies and Grass, 2009).
Natural resistance-associated macrophage proteins (Nramps) are ubiquitous importers of Fe2+ and Mn2+ across cellular membranes into the cytosol (Bozzi and Gaudet, 2021; Cellier and Gros, 2004; Nevo and Nelson, 2006). In humans, Nramp1 extrudes essential metals from phagosomes of macrophages to aid in killing engulfed pathogens, and Nramp2 (DMT1) is expressed at low levels in the endosomes of all nucleated cells and imports Mn2+ and Fe2+ into the cytosol (Pujol-Giménez et al., 2017; Skamene et al., 1998; Vidal et al., 1993). Plant and fungal Nramps aid in Fe2+ and Mn2+ uptake and trafficking, and bacterial Nramps are involved in the acquisition of Mn2+, an essential nutrient (Bozzi and Gaudet, 2021). In addition to the physiological substrates Fe2+ and Mn2+, Nramps can also transport toxic metals like Cd2+ and Hg2+ but exclude the abundant alkaline earth metals like Mg2+ and Ca2+ (Bozzi and Gaudet, 2021).
Recent bacterial Nramp structures reveal a LeuT fold, three stable conformations (outward-open, occluded, and inward-open), and identify the metal-binding site residues, including conserved aspartate, asparagine, and methionine residues (Bozzi et al., 2016b; Bozzi et al., 2019b; Ehrnstorfer et al., 2014; Ehrnstorfer et al., 2017). The metal-binding methionine is essential to select against alkaline earth metals (Bozzi et al., 2016a). This finding is corroborated by the fact that a bacterial Nramp homolog which lacks a metal-binding methionine, NRMT (Nramp-related Mg2+ transporter), can transport Mg2+ (Ramanadane et al., 2022). However, little is known about whether the canonical Nramps can mechanistically distinguish between their physiological substrates (Fe2+ and Mn2+) from non-essential ones like Cd2+ within their broad spectrum of transition metal substrates. Functional studies on Deinococcus radiodurans (Dra)Nramp revealed that Mn2+ and Cd2+ transport differ in their dependence on pH, proton flux, and membrane potential (Bozzi et al., 2019a; Bozzi et al., 2019b). However, we lack high-resolution structural information on binding of different metals to explain these differences.
We present high-resolution structures of DraNramp in three conformations in both Mn2+-bound and metal-free states, providing the first molecular map of the entire Mn2+ transport cycle. The structures along with molecular simulations reveal that Nramps achieve alternate access during transport by adopting distinct Mn2+-coordination spheres in different conformations. These different conformations are also supported by dynamic rearrangements of key polar-residue networks that gate the inner and outer vestibules. This Mn2+ transport cycle also informs on the transport of Fe2+, the other common physiological Nramp substrate, because Fe2+ and Mn2+ have similar coordination preferences and chemical properties (Bozzi and Gaudet, 2021; Davidsson et al., 1989; Kawabata, 2019; Liu et al., 2021). Comparisons with an additional high-resolution structure of DraNramp bound to a non-physiological substrate, Cd2+, and complementary binding and transport measurements and mutational analyses, suggest that Nramps can distinguish physiological from toxic substrates through thermodynamic differences in the conformational landscape of the transport cycle.
Results
DraNramp transports a mostly dehydrated Mn2+ ion
To visualize how the metal substrate is coordinated in Nramps, we determined crystal structures of DraNramp using lipid-mesophase based techniques (Supplementary file 1a). We obtained a structure of wildtype (WT) DraNramp in an occluded state bound to Mn2+ at 2.38 Å by soaking crystals with Mn2+ (WT•Mn2+; Table 1, Figure 1A–B). We resolved a comparable structure using co-crystallization with Mn2+ and the inward-locking mutation A47W (A47W•Mn2+; Supplementary file 1b; Cα RMSD=0.47 Å; all pairwise RMSD values listed in Supplementary file 1c; Bozzi et al., 2016b). The similarity of both structures, including a nearly identical Mn2+-coordination sphere (Figure 1B, Figure 1—figure supplement 1A), suggests that the observations we make based on these two structures are robust. Both structures superimpose best with the published occluded metal-free G45R structure (Bozzi et al., 2019b). Although the inner vestibule is occluded in both structures, WT•Mn2+ and A47W•Mn2+ differ in their TM1a position, with WT•Mn2+ nearly identical to G45R whereas the A47W•Mn2+ TM1a is displaced within the inner vestibule, likely to accommodate the bulky tryptophan sidechain. Therefore, we generally used the WT•Mn2+ structure for analysis of the occluded state. As in the metal-free G45R, the Mn2+-bound occluded structures have a completely sealed outer vestibule and a partially closed inner vestibule, with the Mn2+ occluded from bulk solvent (Figure 1A).
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data1-v2.zip
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data2-v2.zip
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data3-v2.xlsx
Anomalous difference Fourier maps confirmed presence of Mn2+ at the canonical, orthosteric Nramp metal-binding site between the unwound regions of TM1 and TM6 (Supplementary file 1d, Figure 1—figure supplement 1A). In WT•Mn2+, the Mn2+ is coordinated by conserved residues D56, N59, and M230, and backbone carbonyls of A53 and A227 (Figure 1A, Supplementary file 1e). A53 and A227 are pseudosymmetrically related in the inverted repeats of the LeuT fold of DraNramp (Figure 1—figure supplement 1B). The coordination sphere is completed by a water bridging Mn2+ with Q378, a residue previously proposed to directly coordinate Mn2+ in the occluded state (Bozzi et al., 2019b). This yields a coordination number of 6, typical for Mn2+, and a largely dehydrated metal-binding site with a distorted octahedral Mn2+-coordination geometry (Supplementary file 1f), as often observed in other Mn2+-protein complexes (Barber-Zucker et al., 2017; Couñago et al., 2014; Dudev and Lim, 2014).
At the mouth of the outer vestibule, an additional Mn2+ bridges D296 and D369 at the N termini of extracellular helix 2 (EH2) and TM10, respectively (Figure 1—figure supplement 1A). We denote this site as the ‘external site’ and the canonical substrate-binding site as the ‘orthosteric site’. Corroborating the structures, isothermal titration calorimetry (ITC) measurements reveal an endothermic mode of binding and are best fitted with a two-site model for WT (Kd1=190±30 µM, Kd2=1970±520 µM; Figure 1C) and A47W (Kd1=125±5 µM, Kd2=2450±650 µM; Figure 1—figure supplement 2A; all Kd values are in Supplementary file 1g; see Appendix 1 for a description of our ITC data analyses). To determine the affinity of the orthosteric site, we mutated the external-site aspartates. The D296A and D369A substitutions have little impact on Mn2+ transport (Figure 1—figure supplement 2B). The D296A and D369A variants each bind one Mn2+ with Kd=370±30 µM and 420±30 µM respectively (Figure 1—figure supplement 2A), which is closest to Kd1 of WT. Hence, the affinity for Mn2+ at the orthosteric site is higher than at the external site. A metal ion is present at the external site in all inward-open and occluded metal-bound DraNramp structures, but not outward-open structures, as opening the outer vestibule separates D296 and D369 and disrupts the site (Figure 1—figure supplement 3A). D296 and D369 are not conserved across Nramps, but they are more conserved within bacterial clade A, and there is a general abundance of acidic residues in the corresponding loop regions across all clades (Figure 1—figure supplement 3B–C). At present, our results provide little evidence of a biological role for this previously unidentified external site; perhaps the concentration of acidic residues at the mouth of the outer vestibule (Figure 1—figure supplement 3D) provides electrostatic attraction for metal cations.
Snapshots of the complete Mn2+ transport cycle by DraNramp
We also determined high-resolution DraNramp structures in metal-free occluded (WT) and Mn2+-bound inward-open (M230A•Mn2+) states and re-refined a Mn2+-bound outward-open conformation (G223W•Mn2+; Table 1, Supplementary file 1a and b). Along with the published structures of outward-open metal-free G223W and inward-open metal-free ‘Patch’ (which has a patch of mutations in intracellular loops) (Bozzi et al., 2016b; Bozzi et al., 2019b), these new structures allow us to map the entire Mn2+ transport cycle to three major conformations, each in Mn2+-bound and metal-free states (Figure 2A–B). By ordering and comparing these six structures, we outline a molecular mechanism by which metal substrate binds, induces conformational change, and is released.
TMs 1, 5, 6, and 10 move the most as the Mn2+-binding site accessibility switches from outward to inward across the conformations (Figure 2B–C, Figure 2—figure supplement 1), as also highlighted in previous studies (Bozzi and Gaudet, 2021; Bozzi et al., 2019b; Ehrnstorfer et al., 2017). As Mn2+ binds to the outward-open state, TM10 tilts toward TM1b, the upper half of TM5 toward TM7, and TM6a toward TM11 to seal the outer vestibule and yield an occluded state. The lower half of TM5 also moves away from TM1a, allowing it to swing upward to open the inner vestibule in the following transition. This swing of TM1a is the only noteworthy difference between the Mn2+-bound occluded and inward-open states, allowing release of the Mn2+ into the inner vestibule. TM1a swings to a similar angle in the new inward-open M230A•Mn2+ as in the previous low-resolution inward-open metal-free structure (Bozzi et al., 2016b), and structures of the homologous Eggerthella lenta Nramp-related magnesium transporter (EleNRMT), LeuT, and serotonin transporter (Coleman et al., 2019; Krishnamurthy and Gouaux, 2012; Ramanadane et al., 2022). Thus, most of the structural reorganization in Nramps occurs in the shift from outward open to occluded. The three metal-free DraNramp conformations are similar to their corresponding Mn2+-bound structures, suggesting that once Mn2+ is released, the conformational transitions are reversed, including passing through an occluded metal-free intermediate, to reach the outward-open metal-free conformation ready to accept Mn2+ (Figure 2A–B, Figure 2—figure supplements 1 and 2A).
The substrate-binding TM1 and TM6 are well-resolved in our structures (Figure 2—figure supplement 2B) and pairwise superpositions reveal how their motions contribute to the conformational changes across the transport cycle (Figure 2C). TM6a tilts 22° and the unwound region of TM6 becomes more helical as it moves toward the Mn2+ to close the outer gate. The central unwound regions of TM1 and TM6 are closest in the occluded state, resulting in an almost dehydrated Mn2+-coordination sphere. Finally, the inner vestibule opens when TM1a tilts upward by 32°, increasing the distance between TM1a and TM6b by ~8 Å (Figure 2C). TM6b is largely static relative to the protein core, although in the inward-open structure it has high B-factors (Figure 2D), indicating that the interaction with TM1a stabilizes TM6b to close the inner vestibule.
Different conformations have distinct Mn2+ coordinations
Our Nramp structures provide snapshots of the complete Mn2+-coordination sphere geometries in each conformation (Figure 3). Two new structures of DraNramp point-mutants reveal the coordination of Mn2+ in the inward-open state: M230A•Mn2+ and D296A•Mn2+ (Cα RMSD of 0.42 Å). The Mn2+ is in the same location of the orthosteric site as in the occluded state, as confirmed by anomalous diffraction for D296A•Mn2+ (Supplementary file 1d, Figure 3—figure supplement 1). As in the occluded state, the Mn2+ binds D56, N59, and the A227 carbonyl, with a water replacing M230 in M230A•Mn2+ (Figure 3A, Supplementary file 1e). However, with TM1a displaced, the A53 carbonyl no longer coordinates Mn2+; instead, the Y54 carbonyl approaches Mn2+ at a longer distance of 3.1 Å. Two more waters, one bound to Q378 and another from the inner vestibule, complete a seven-coordination sphere resembling a pentagonal bipyramidal geometry with substantial distortion (Supplementary file 1f). Seven coordination is infrequent but found in Mn2+-coordinating proteins like MntR (Chen and He, 2008; Glasfeld et al., 2003). Our inward-open structures provide the first evidence that Y54 participates in the Mn2+ transport cycle.
The new occluded and inward-open Mn2+-bound structures have a monodentate coordination of D56 with Mn2+. For consistency, we reinterpreted the outward-open G223W•Mn2+ map (PDB ID: 6BU5) with a monodentate coordination of D56 with Mn2+ instead of previously modeled bidentate interaction (Bozzi et al., 2019b); the local and global model statistics are very similar to the original structure (Supplementary file 1e). Re-refined G223W•Mn2+ has six Mn2+-coordinating ligands: D56, N59, M230, carbonyl of A53 and two waters (Figure 3A, Supplementary file 1e); the overall geometry resembles a distorted octahedron as in the occluded structure (Supplementary file 1f). The coordination spheres in all conformations of the transport cycle are well defined as confirmed by the 2Fo-Fc maps of the closeup snapshots of their Mn2+-bound orthosteric site (Figure 3—figure supplement 1C, Figure 5—figure supplement 2).
Comparing the three Mn2+-bound conformations (Figure 3A), the carbonyls of A53 in TM1a and A227 in TM6b alternately coordinate Mn2+ in the outward- and inward-open structures respectively, and both residues interact with Mn2+ in the occluded structure. The pseudosymmetrically related A53 and A227 may thus act as hinges altering the Mn2+-coordination sphere as TM1 and TM6 move in turn to open the gates during Mn2+ transport. Furthermore, as DraNramp switches from outward- to inward-open, Y54 progressively swings downward, acting as a gate in concert with TM1a’s upward swing to open the inner vestibule and allow metal release (Figure 3A–B). Our inward-open structures also suggest that the Y54 carbonyl may participate in Mn2+ release through interaction with the metal ion (Figure 3A). All 3796 Nramp sequences in our alignment have either a tyrosine or a phenylalanine at this position and Y54 is completely conserved among bacterial clades A (including DraNramp) and C, while the position is 40% and 100% phenylalanine among eukaryotes and bacterial clade B, respectively (Figure 3—figure supplement 3). To evaluate the significance of Y54 in Mn2+ transport, we purified and reconstituted into proteoliposomes the Y54A and Y54F variants. While Y54F has near-wildtype Mn2+-transport activity, Y54A nearly eliminates Mn2+ transport (Figure 3C), indicating that an aromatic ring is essential for the gating motion required for transport.
We also measured the Mn2+-binding affinity of the constructs that yielded inward- or outward- open structures, M230A and G223W, respectively. Mn2+ is present at the external site in the inward-open M230A•Mn2+ (Figure 1—figure supplement 3A), and consistently, the ITC data fit a two-site model (Kd1=215±65 µM assigned to the orthosteric site, Kd2=4600±1300 µM for the external site; Figure 3D and Appendix 1). The ITC data for G223W with Mn2+ fits only in a one-site model (Kd=440±15 µM; Figure 3D), which we assign to the orthosteric site because the opening of the outer vestibule displaces TM10, disrupting the external site (Figure 1—figure supplement 3A).
Mn2+ binding does not significantly alter the three main DraNramp conformations
To compare the metal-bound states to analogous metal-free states of the transport cycle, we determined two metal-free occluded structures of wildtype DraNramp at a higher resolution than the previously reported G45R structure (Bozzi et al., 2019b), which we refer to as WT (2.38 Å; Supplementary file 1b) and WTsoak (2.36 Å; Table 1). The crystal used for WTsoak was mock-soaked (with no metal in the soaking solution). WT and WTsoak are nearly identical, confirming that the soaking process does not influence the conformational state. We analyzed WTsoak, unless otherwise noted. In WTsoak, we observed density but no anomalous signal at the orthosteric site and modeled a water molecule at the position where Mn2+ sits in the occluded state (Figure 3A). WTsoak is nearly identical to WT•Mn2+, indicating that the occluded conformation is unchanged by the presence of metal ion substrate, and the metal-binding site is instead filled by ordered water molecules.
We used previously reported inward-open ‘Patch’ and outward-open G223W metal-free structures for analysis of the Mn2+ transport cycle (Figures 2B and 3A; Bozzi et al., 2016b; Bozzi et al., 2019b). These structures, resolved at lower resolution than the ones described here, have no density at the orthosteric site. This is consistent with a more flexible organization of a metal-binding site open to bulk aqueous solvent. In contrast, metal-free WT has a more ordered orthosteric site, suggesting a stable occluded intermediate in the switch from inward- to outward-open.
Polar networks latch the gates to achieve alternating access
Vestibules providing access to the orthosteric site from the extracellular or intracellular side alternately open from the motions of DraNramp’s TMs 1, 5, 6, and 10, which form the outer and inner gates during Mn2+ transport (Bozzi et al., 2020; Bozzi et al., 2019b). To pinpoint protein features that enable these motions, we used our Mn2+-bound structures to identify interaction networks with the following attributes: (i) they contain conserved polar residues from at least one of the four mobile helices; (ii) they line the gates; and (iii) they rearrange between the three resolved protein conformations (Figure 4A).
Two networks seal the outer gate. In the occluded and inward-open conformations, Q378 interacts with two waters, one coordinating the orthosteric Mn2+ and the other interacting with D56 and the carbonyl of T130. This network is disrupted in the outward-open conformation as Q378 and the rest of TM10 swing outward to open the outer vestibule (Figure 4B). In the second network, T228, N275, and N82 interact via a water in the occluded and inward-open states, but not in the outward-open state, where the extended unwound region of TM6a positions T228 farther from the orthosteric site and N275 and N82 (Figure 4C). This T228 network helps rearrange TM6a, closing the outer vestibule in the occluded state and generating a nearly dehydrated Mn2+-coordination sphere (Figure 1B). As the inner gate opens to release Mn2+, both networks persist, ensuring that the outer gate remains closed in the inward-open conformation.
All six residues in the Q378 and T228 networks are completely conserved across bacterial clade A and highly conserved across all Nramps (Figure 4—figure supplement 1). To investigate the robustness of these networks, we performed duplicates of molecular dynamics (MD) simulations starting in each of the three conformations and confirmed that within the first 250 ns of these simulations, the T228 and Q378 networks persist in simulations of the occluded and inward-open states and remain broken in simulations starting in the outward-open state (Figure 4—figure supplement 2A–C). This includes the coordinated waters, for example the water at the center of the Q378 network is present in more than 50% of the frames in occluded-state simulations (Figure 4—figure supplement 2E). In line with the X-ray snapshots and simulation outcomes and highlighting their key function in the conformational cycle of Nramps, mutations of any of the six residues across these networks reduces Mn2+ transport by DraNramp (Bozzi et al., 2016a; Bozzi et al., 2019a; Bozzi et al., 2020; Bozzi et al., 2019b).
The inner vestibule is gated by rearrangements of residues in two other polar networks, namely those of R244 and Q89 (Figure 4D–E). In the outward-open state, R244 forms an ion pair with E176 and interacts with the TM1a backbone, as does D263, keeping TM1a and TM5 close and the inner gate closed (Figure 4D; Bozzi et al., 2020). In the occluded state, the E176-R244 interaction breaks and TM5 moves away from TM6b, creating space for TM1a to swing up and open the inner vestibule in the inward-open state. Accordingly, the E176–R244 ion pair is stable in MD simulations of the outward-open state, while these residues are >10 Å apart in simulations of the occluded and inward-open states (Figure 4—figure supplement 2D). Supporting the importance of the R244 network, E176 is 100% and R244 is 85% conserved across all Nramps (Figure 4—figure supplement 1) and mutation of either residue reduces Mn2+ transport by DraNramp (Bozzi et al., 2020).
In the Q89 network, Q89 hydrogen-bonds with Y54 and H237 to seal the inner gate in the outward-open (Bozzi et al., 2020) and occluded states (Figure 4E). In the inward-open structure, the Q89–H237 hydrogen bond is broken and a rearranged Y54–Q89 hydrogen bond buttresses the opening of the inner vestibule (Figure 4E). As discussed above, Y54 is conserved and important for Mn2+ transport by DraNramp (Figure 3C and Figure 3—figure supplement 1C). Similarly, Q89 and H237 are conserved (Figure 4—figure supplement 1), and mutations to either residue impair Mn2+ transport in a cell-based assay (Bozzi et al., 2020). Cysteine accessibility measurements showed that mutations to Q89 or H237 render the outer vestibule solvent-inaccessible (Bozzi et al., 2020), indicating that disrupting the Q89 network likely prevents closing of the inner gate.
H232 (TM6b) is conserved across all Nramps (Figure 4—figure supplement 1), highlighting its importance. H232 sits below the orthosteric site and forms a network conserved across all conformations, with water-mediated hydrogen bonds to E134 (involved in proton transfer to the salt-bridge residues in TMs 3 and 9) (Bozzi et al., 2019a; Bozzi et al., 2019b; Ehrnstorfer et al., 2017) and T130 in TM3 (Figure 4F). Waters occupy these two sites in MD simulations in all states, especially the water coordinated between H232 and T130 (Figure 4—figure supplement 2F–G). However, in M230A•Mn2+ H232 flips to interact with the G330 carbonyl (TM8; Figure 4F), suggesting that it may transiently move during the conformational cycle. Indeed, while the H232 sidechain rotamer is stable in MD simulations of the outward-open state, it explores other rotamers in simulations of the occluded and inward-open states (Figure 4—figure supplement 3).
Unlike Mn2+ binding, Cd2+ binding to DraNramp is exothermic
Nramps transport divalent transition metals quite promiscuously, including both physiological (Fe2+ and Mn2+) and non-physiological substrates (Cd2+, Zn2+, Co2+, Ni2+, Pb2+), but select against alkaline earth metals (Mg2+, Ca2+) (I Bannon et al., 2002; Bozzi et al., 2016a; Bozzi and Gaudet, 2021; Bozzi et al., 2019b; Ehrnstorfer et al., 2017; Illing et al., 2012; Sacher et al., 2001). DraNramp transports Cd2+ well, but without concomitant proton flux and with weaker voltage dependence (Figure 5—figure supplement 1A; Bozzi et al., 2019a; Bozzi and Gaudet, 2021; Bozzi et al., 2019b). To better understand the underlying mechanistic differences, we compared the binding, transport, and structures of DraNramp with Mn2+ and Cd2+.
In contrast to endothermic binding of Mn2+, ITC measurements show exothermic binding of Cd2+ to WT DraNramp (Figure 5—figure supplement 1B; see Appendix 1 for details of the ITC analyses). Like for Mn2+, the Cd2+ isotherm fits best in a two-site model (Kd1=55±15 µM, Kd2=220±20 µM). Both D296A and D369A—containing mutations at the external metal-binding site—showed exothermic binding but fit best in a one-site model (Figure 5A). Based on their affinity (Kd=120±1 µM for D296A and Kd=70±10 µM for D369A), we conclude that the orthosteric site has higher affinity toward Cd2+ than the external site (Figure 5A–B). Both the orthosteric and external site shows higher affinity toward Cd2+ than Mn2+. The affinities are comparable at the orthosteric site, in low micromolar range for both metals (Kd of 190±30 µM for Mn2+ vs. 55±15 µM for Cd2+), whereas the external site has a much higher affinity toward Cd2+ (Kd of 1970±520 µM for Mn2+ vs. 220±20 µM for Cd2+). Consistent with its inability to transport Mg2+, a representative alkaline earth metal (Bozzi and Gaudet, 2021; Bozzi et al., 2019b), DraNramp does not bind Mg2+ (Figure 5A). These are the first ITC measurements comparing the binding of different metals to an Nramp transporter and they show clear differences in the binding mode and affinity of different substrates toward DraNramp (Figure 5A, Figure 5—figure supplements 2 and 3). In contrast to DraNramp, previous ITC studies showed endothermic binding of Cd2+ to the Staphylococcus capitis Nramp homolog (ScaDMT) with 29 µM affinity (Ehrnstorfer et al., 2014). However, in the absence of ITC data with other metals, it is not known whether ScaDMT also shows differences in the mode and affinity of binding to different metals like DraNramp.
A Cd2+-bound structure helps explain functional differences
To understand the differences in binding and transport of Mn2+ and Cd2+, we determined a 2.5 Å Cd2+-bound structure by soaking WT DraNramp crystals with 2 mM Cd2+ (Table 1). In agreement with the ITC data, WT•Cd2+ shows Cd2+ ions at both the external and orthosteric sites as confirmed by anomalous signal (Figure 5—figure supplement 1C–D, Supplementary file 1d). The higher affinity for Cd2+ than Mn2+ at the external site suggests that local geometry favors non-physiological Cd2+ over the physiological substrate Mn2+ (Figure 1—figure supplement 3A). The low affinity of this site for Mn2+, compared to Cd2+, is consistent with the idea that this external site plays a role in general electrostatic attraction of substrate to the orthosteric site, rather than as a finely tuned binding site specific for Mn2+.
Interestingly, WT•Cd2+ is inward open (Figure 5—figure supplement 1C), although both mocked-soaked and Mn2+-soaked crystals under otherwise equivalent conditions yielded occluded structures (WTsoak and WT•Mn2+, respectively; Figures 1 and 3). The larger ionic radius of Cd2+ (~0.95 Å) compared to Mn2+ (~0.82 Å) (Kumarevel et al., 2005; Vashishtha et al., 2016) and different preferred coordination geometry may increase the stability of the inward-open state with Cd2+ bound. Compared with the inward-open M230A•Mn2+, D56 adopts a different rotamer in WT•Cd2+ and does not coordinate the Cd2+ bound at the orthosteric site (Figure 5C). The rest of the coordination sphere is similar and includes N59, M230, the A227 carbonyl, the Y54 carbonyl, a water that coordinates Q378 and another water from inner vestibule. Cd2+ has six coordinating ligands, and the distortion from ideal octahedral geometry is more pronounced compared to Mn2+ (Supplementary file 1f). The coordination distances are larger for Cd2+ than for Mn2+ (Supplementary file 1e), as seen in other proteins (Begg et al., 2015; Yokoyama et al., 2012) and consistent with its larger ionic radius and distinct charge distribution. Cd2+ is a soft metal, likely explaining why it retains coordination by the softer sulfur ligand of M230 (Cammack and Hughes, 2008) but not the hard oxygen of D56, although D56 can still provide favorable electrostatics.
Cd2+ binding is prone to perturbations and favors the inward-open state
ITC data analysis of different DraNramp variants highlights additional differences between Mn2+ and Cd2+ binding and how their binding affects the conformational preferences of DraNramp (Figure 5B, Figure 5—figure supplements 2–3). A general observation is that, at the orthosteric site, Mn2+ and Cd2+ binding are entropy- and enthalpy-driven, respectively, and entropy contributions (-TΔS) are smaller for Cd2+ than Mn2+ (Appendix 1), suggesting that Cd2+ binding conformationally constrains the protein more, leads to less solvent release, or both (Ferrante and Gorski, 2012; Olsson et al., 2008).
ITC with the orthosteric-site mutants, M230A and D56A, shows metal-specific behavior. For M230A, Mn2+ binding fits with a two-site model and Cd2+ fits only with a one-site model (Figure 5B and Figure 5—figure supplement 2, Appendix 1). M230 is thus a crucial ligand for Cd2+ but not Mn2+, which is further reflected in the transport behavior where M230A affects Cd2+ transport drastically but has negligible effect on Mn2+ (Bozzi et al., 2016a; Bozzi et al., 2019a). Data for D56A fit better with a two-site model with both metals, although D56A does not transport either metal (Bozzi et al., 2019b). This suggests that D56 is more important for catalyzing transport than for substrate binding. Data for variants with an additional mutation at the external site (D56A-D296A, D56A-D369A, M230A-D296A, and M230A-D369A) fit only with a one-site model for Mn2+-binding isotherms, which we assigned as binding to the orthosteric site. However, all four double mutants showed complete loss of Cd2+ binding, which is expected for ones with M230A—which eliminates Cd2+ binding at the orthosteric site on its own, as we have also shown previously (Bozzi et al., 2016a)—but more surprising for double mutants with D56A. Thus, binding of the preferred physiological Mn2+ substrate at the orthosteric site is more robust to perturbations than binding of Cd2+, a toxic metal.
The conformation-locking mutants—A47W, which is outward-closed (Bozzi et al., 2016b) and crystallized in an occluded state, and outward-open G223W (Bozzi et al., 2019b)—also show different behavior with Mn2+ and Cd2+ (Figure 5—figure supplement 3). ITC with Mn2+ agrees with the structural data, with two-site fit for A47W (Figure 1—figure supplement 2A), and one-site fit for A47W-D296A, A47W-D369A, and G223W (assigned to the orthosteric site). However, with Cd2+, the ITC data for A47W fit best with a one-site model and we observed no binding with A47W-D296A, A47W-D369A, and G223W, indicating that there is no significant affinity for Cd2+ at the orthosteric site of A47W and G223W. This is consistent with our inability to obtain Cd2+-bound structures in conformations other than inward-open, despite trying both co-crystallization and soaking with A47W and G223W. These data suggest that Cd2+ has highest affinity for the inward-open state and low affinity for other states, in agreement with the fact that soaking of WT crystals (which yielded the metal-free occluded WTsoak structure) produced an inward-open Cd2+-bound structure.
Gating network differences in the Cd2+-bound structure are restricted to D56
To better understand the Cd2+ transport mechanism given that its binding at the orthosteric site seems less optimal and robust than Mn2+, we compared the gating networks described above for the inward-open Mn2+- and Cd2+-bound structures. The Q89, T228, and R244 networks are essentially identical in the inward-open state, regardless of whether Mn2+ or Cd2+ is bound (Figure 5—figure supplement 4A–C).
The other two networks, both of which involve D56, have some differences. In the Q378 network gating the outer vestibule, D56 does not coordinate Cd2+ but does interact with the conserved water that connects the metal ion with Q378 in all outward-closed structures (Figure 5—figure supplement 4D). This preserves a connection to both Cd2+ and Q378 to close the outer gate. Most of the H232 network is similar in the Mn2+- and Cd2+-bound structures, except for the orientation of D56 relative to E134 (Figure 5—figure supplement 4E).
Mutations of several residues in these polar networks reduced transport in proteoliposome-based assays (Figure 5—figure supplement 4F), corroborating results from cell-based assays (Bozzi et al., 2020). Y54A, H237A, and H232A mutations cause the largest decreases, Q89A causes a moderate decrease, and Y54F is similar to WT. Mn2+ and Cd2+ follow similar trends although with less voltage dependence for Cd2+ than Mn2+. These results support the idea that although Mn2+ and Cd2+ bind differently at the orthosteric site, the flexibility of the D56 sidechain enables the networks of polar residues that gate the outer and inner vestibules to engage and enable transport of both metals at similar rates.
Discussion
Our high-resolution structures in different conformations and analysis of substrate-binding affinities reveal a molecular map of the Mn2+ import pathway in DraNramp (Video 1). Binding of the Mn2+ substrate at the orthosteric site takes on different coordination geometries through the three main conformational states in the transport cycle, all of which deviate substantially from ideal, consistent with the moderate binding affinities we measured. We identified several networks of polar interactions that gate both the outer and inner vestibules. Structures of DraNramp bound to Mn2+ and Cd2+ and corresponding analyses of substrate binding demonstrate that the orthosteric site shows similar affinity for physiological (Mn2+) and toxic (Cd2+) substrates, but Cd2+ binding is less robust to various perturbations.
Our structures of one Nramp homolog, DraNramp, in all conformations of the transport cycle provide an opportunity to update the overview of the conformational cycle, focusing on the polar interaction networks that gate the outer and inner vestibules (Figure 6A). The structures also provide the first molecular view of how local changes in the coordination spheres and the global conformational transitions coordinate to facilitate Mn2+ transport. Starting in the outward-open state, Mn2+ entry into the outer vestibule and binding at the orthosteric site triggers the transition to the occluded state by two major rearrangements: (i) closing of the outer gate as the T228 and Q378 networks form and support the reorientation of TM6a and TM10, respectively, and (ii) partial opening of the inner gate through motion of TM5 and breaking of the R244 network. TM6a and TM10 approach the orthosteric site and directly (A227 in TM6a) or indirectly (Q378 in TM10 through a conserved water) coordinate Mn2+ and restrict solvent access from the extracellular side in the occluded structure. A227 replaces a water of the outward-open Mn2+-coordination sphere, thus retaining a six-coordination geometry in both conformations. Mn2+ prefers octahedral (six) coordination (Chen and He, 2008; Dudev et al., 2006; Dudev and Lim, 2014) as in the outward-open and occluded structures, although we observed distortions (Supplementary file 1f).
The inward-open state is achieved when TM1a swings up, rupturing the Q89 network and opening the inner gate. The coordination sphere changes again as a water from the open inner vestibule replaces A53 from TM1a. The TM1a swing introduces Y54 as a long-range seventh ligand. The seven-coordination of Mn2+ in the inward-open structure is less favored and more distorted, which may facilitate Mn2+ release as a solvated ion in the cytosol. In contrast, the more favorable six-coordination in the outward-open and occluded structures could help energize the global conformational changes. Our metal-free structures indicate that after Mn2+ release to the cytosol, the protein resets to the outward-open state through the same occluded state. The residues in the gating networks are highly conserved, suggesting that their role is conserved across the Nramp family.
Compared to other proteins that bind Mn2+ but are not metal transporters, like the Mn2+ regulator MntR (Glasfeld et al., 2003) and PsaA, the solute-binding protein (SBP) domain of an ATP-binding cassette transporter (Couñago et al., 2014), we observe longer Mn2+-coordinating bond lengths for DraNramp (Figure 6—figure supplement 1). Typical manganese-oxygen bonding distances are 2.1–2.5 Å, although ‘weak interactions’ (2.6–3.2 Å) are occasionally part of a Mn2+-coordination sphere (Harding, 2000; Harding, 2001). Metal-sulfur bonding distances are longer owing to the greater van der Waals radius of sulfur (Rulísek and Vondrásek, 1998). The non-ideal metal-ligand bonding distances and angles we observe in DraNramp may allow it to avoid getting trapped in an energy minimum and thus keep moving through the conformational transitions required to transport Mn2+.
While the Mn2+ transport mechanism and associated conformational changes in DraNramp differ appreciably from other LeuT-fold transporters (Bozzi and Gaudet, 2021; Bozzi et al., 2019b), the presence of key residues defining the extracellular and intracellular gates is a common theme (Coleman et al., 2016; Coleman et al., 2019; Faham et al., 2008; Krishnamurthy and Gouaux, 2012; Krishnamurthy et al., 2009; Perez et al., 2012; Shimamura et al., 2010; Watanabe et al., 2010). Comparing the gating networks described for these transporters with DraNramp, the positions and nature of the gating networks are generally not conserved, but two common themes emerge. First, opening (or closing) a particular vestibule often involves a pair of changes, as in the DraNramp intracellular vestibule (Figure 6A). Such paired changes have been associated with ‘thin’ and ‘thick’ gates in Mhp1 (Simmons et al., 2014) and LeuT (Krishnamurthy and Gouaux, 2012), that is, gates based on sidechain and helix motions, respectively. Second, some gating residue positions are shared, but the networks are not, indicating that different families have evolved analogous networks to stabilize equivalent conformations (Coleman et al., 2016; Krishnamurthy and Gouaux, 2012).
The accumulated DraNramp structures also provide clues as to the thermodynamic landscape of the transport process (Figure 6B). Interestingly, both metal-free and Mn2+-bound WT DraNramp crystallized in the occluded conformation, whereas the open conformations were achieved through conformation-locking (Bozzi et al., 2019b) or mutations of functionally important residues (Bozzi et al., 2016b). This occluded state more closely resembles an inward-open conformation (Bozzi et al., 2019b) and is an important intermediate in the Mn2+ transport cycle, with both local changes in coordination geometry and global changes in protein conformation in comparison to the outward-open state. Physiologically, we naïvely expect outward-open, rather than occluded, to be the preferred substrate-free conformation. However, the occluded or inward-open states of other LeuT-fold importers are also more stable, with changes in environmental conditions or the presence of substrate stabilizing specific states and lowering barriers to conformational transitions (Del Alamo et al., 2022). For example, SGLT1 and DraNramp require a negative membrane potential to transport substrates and this negative membrane potential stabilizes the outward-open state of SGLT1 (Bozzi et al., 2019a; Loo et al., 1998). Overall, our DraNramp structures suggest that the occluded state is most stable (at least in the absence of a membrane potential). The occluded and inward-open states may be energetically similar as relatively small perturbations yielded inward-open structures (M230A•Mn2+ and D296A•Mn2+). Furthermore, our previous cysteine accessibility data indicate that the energy barriers between states are low enough for the protein to readily sample the outward- and inward-open states in cell membranes in the absence or presence of metal substrate (Bozzi et al., 2016b; Bozzi et al., 2019b).
In contrast to WT DraNramp, ScaDMT was crystallized in an inward-open state, although its TM1a was deleted from the protein construct (Ehrnstorfer et al., 2014). This deletion would prevent inner vestibule closure, and thus likely affects its energetically preferred conformation. In the case of the Eremococcus coleocola Nramp homolog (EcoDMT), both substrate-free and inhibitor-bound conformations are outward open (Ehrnstorfer et al., 2017; Manatschal et al., 2019), suggesting that the outward-open state is its most stable state. Of note, EcoDMT was crystallized in detergent, whereas the DraNramp structures were obtained in a monoolein lipid bilayer environment. Further studies will be needed to determine to what extent the thermodynamic landscape we begin to outline here for DraNramp is conserved in other Nramp homologs.
Previous metal selectivity studies indicate that Nramps import different substrates with distinct mechanisms (Bozzi et al., 2016a; Bozzi et al., 2019a; Bozzi and Gaudet, 2021; Ehrnstorfer et al., 2014; Nevo and Nelson, 2006). These differences correlate with chemical properties (Irving and Williams, 1953) and preferred coordination chemistry (Ma et al., 2009; Singh et al., 2020) of these transition metals as observed in other transition metal binding proteins like the Psa permease (Begg et al., 2015; Couñago et al., 2014), and cation diffusion facilitators (CDFs) (Barber-Zucker et al., 2017). That WT•Cd2+ retains the M230 sulfur as a coordinating ligand but excludes the D56 carboxylate can be rationalized by the fact that Cd2+ is a softer metal than Mn2+. This difference in coordination likely alters the pKa of nearby residues and water molecules to perturb the proton pathway such that DraNramp co-transports protons with Mn2+ but not Cd2+ (Bozzi et al., 2019a; Bozzi et al., 2019b), although our structures do not yet fully elucidate this mechanistic difference. Moreover, owing to its larger radius, Cd2+ tends to form weaker complexes than Mn2+ with comparable coordination numbers in complexes dominated by harder ligands (Rulísek and Vondrásek, 1998; Singh et al., 2020). The architecture of the orthosteric site in DraNramp appears to facilitate six- or seven-coordinated metal complexes, and thus most Cd2+-bound DraNramp conformations will be less stable than the Mn2+-bound states. Our ITC data also indicate that binding of the toxic Cd2+ to DraNramp is less robust to perturbations compared to its physiological substrate Mn2+. Furthermore, Cd2+ only binds well to the inward-open state, and its binding is exothermic rather than endothermic for Mn2+ (Figure 5A). The crystallographic data similarly suggest that inward-open is the most stable Cd2+-bound state (Figure 6B), based on the following observations: (i) soaking crystals of occluded WT yielded inward-occluded WT•Cd2+, (ii) co-crystallization efforts yielded no other Cd2+-bound structures, and (iii) soaking Cd2+ into crystals of outward-locked G223W yielded very poor diffraction and no structures.
Overall, our data show that the orthosteric metal-binding site of DraNramp, conserved across all Nramps, is best suited to the physiological substrate Mn2+ (and likely the similar ion Fe2+). The distinct interactions of Nramps with Mn2+ and Cd2+ could be leveraged for the design of therapies for metal toxicity and prevention strategies for toxic metal accumulation in crops. These results also lay a foundation for future studies of how metal ion transporters like Nramps evolve their substrate selectivity, for example in response to different environmental conditions. Finally, the first complete set of structures with the same homolog, both in substrate-free and substrate-bound states, suggest a substrate-specific thermodynamic landscape of the transport cycle and provide a framework for future experiments and simulations to fully define this landscape, and for comparisons to other LeuT-fold transporters.
Materials and methods
Cloning and protein expression vectors
Request a detailed protocolThe DraNramp WT and mutant constructs were cloned into pET21a-N8H (Bozzi et al., 2016a). Primer sequences for the mutations are listed in Supplementary file 1h. All constructs used for crystallization had a truncation of 31 residues at the N-terminus (ΔN31, which does not impair metal transport), except for A47W which was full-length and transport deficient (Bozzi et al., 2016b). For proteoliposome-based transport assays, the full-length versions of each construct were used. For ITC, full-length DraNramp constructs were cloned into pET21-NStrep (Bozzi et al., 2016a) to avoid background signal from metals (Mn2+, Cd2+) binding to the His-tag. All mutations were introduced by site-directed mutagenesis using the Quikchange mutagenesis protocol (Stratagene) and confirmed by Sanger DNA sequencing.
Protein expression
Request a detailed protocolProtein expression was performed as previously described (Bozzi et al., 2016a). Briefly, transformed Escherichia coli C41(DE3) (Lucigen) were induced with 0.1 mM isopropyl-β-D-thiogalactopyranoside and cultured at 18 °C for 16 hr. Cell pellets from 10 L of culture were harvested and flash frozen in liquid nitrogen.
Protein purification for crystallography
Request a detailed protocolCells were thawed and resuspended in 50 mL load buffer (20 mM sodium phosphate, pH 7.5, 55 mM imidazole pH 7.5, 500 mM NaCl, 10% (v/v) glycerol) supplemented with 1 mM PMSF, 1 mM benzamidine, 0.3 mg/mL DNAse I and 0.3 mg/mL lysozyme and lysed by sonication on ice (six cycles of 45 s with a Branson Sonifier 450 under duty cycle of 65% and output 10). Lysates were cleared by centrifuging for 20 min at 20,000 rpm (Beckman JA-20) and membranes pelleted from the supernatant by ultracentrifugation at 45,000 rpm (Beckman type 45Ti) for 70 min. Membranes were homogenized in 70 mL load buffer using a glass Potter-Elvehjem grinder, solubilized for 1 hr in 1% (w/v) n-dodecyl-β-D-maltopyranoside (DDM), then ultracentrifuged at 35,000 (Beckman type 45Ti) for 35 min to remove insoluble debris. Pre-equilibrated Ni-Sepharose beads (2 mL; GE Healthcare) were incubated with the supernatant for 90 min at 4 °C, then washed with 20 column volumes (CV) of each of the following buffers sequentially (i) load buffer containing 0.03% DDM, (ii) load buffer containing 0.5% lauryl maltose neopentyl glycol (LMNG), and (iii) load buffer containing 0.1% LMNG. Protein was eluted in 20 mM sodium phosphate, pH 7.5, 450 mM imidazole pH 7.5, 500 mM NaCl, 10% (v/v) glycerol, 0.01% LMNG, concentrated to <0.5 mL in a 50 kDa molecular weight cutoff (MWCO) centrifugal concentrator (EMD Millipore), and purified by size exclusion chromatography (SEC) using a Superdex S200 10/300 (GE Healthcare) pre-equilibrated with SEC buffer (10 mM HEPES pH 7.5, 150 mM NaCl, 0.003% LMNG). Peak protein fractions enriched in DraNramp were combined, concentrated to ~25–40 mg/mL using a 50 kDa MWCO centrifugal concentrator, aliquoted and flash frozen in liquid nitrogen and stored at –80 °C. Purifications of each protein construct were performed at least twice and resulted in similar data.
Purification of DraNramp for ITC
Request a detailed protocolTo purify protein for ITC, harvested cells expressing strep-tagged DraNramp from 10 L of culture were resuspended in 50 mL of buffer W (100 mM Tris, pH 8.0, 150 mM NaCl), and membranes were isolated, homogenized and solubilized in 1% DDM as above. The supernatant was incubated with Strep-Tactin Superflow resin (3 mL; IBA) pre-equilibrated with buffer W+0.03% DDM and washed with the following buffers sequentially: (i) 1 CV buffer W+0.03% DDM, (ii) 2 CV buffer W+0.5% LMNG and (iii) 2 CV buffer W+0.1% LMNG. Protein was eluted with 3 CV buffer W+0.01% LMNG+2.5 mM desthiobiotin. The eluted protein was concentrated up to 2.5 mL using a 50 kDa MWCO centrifugal concentrator and buffer-exchanged into 150 mM NaCl, 10 mM HEPES, pH 7.5, and 0.003% LMNG using disposable PD-10 desalting columns (GE healthcare). Protein was concentrated to ∼2.5 mg/mL a 50 kDa MWCO centrifugal concentrator and flash frozen in liquid nitrogen and stored at –80 °C. Purifications of each protein construct were performed at least twice and resulted in similar data.
DraNramp crystallization
Request a detailed protocolCrystallization of all constructs was performed using lipidic cubic phase (LCP). Protein was mixed with monoolein in 1:1.5 volume ratio using the syringe reconstitution method (Bozzi et al., 2019b). The protein bolus (60 nL) and 720 nL precipitant were dispensed onto custom-made 96 well glass sandwich plates using an NT8 drop-setting robot (Formulatrix). Metal-free crystals (WT) and metal supplemented (5 mM MnCl2) co-crystals (A47W•Mn2+, M230A•Mn2+, D296A•Mn2+) were grown in different precipitant conditions (Supplementary file 1a), harvested within 7–10 days (after reaching their optimal size of 30–40 μm rods) using mesh loops (MiTeGen) and flash-frozen in liquid nitrogen prior to data collection. Some structures (WT•Cd2+, WT•Mn2+, WTsoak) were obtained by soaking WT DraNramp crystals grown in metal-free precipitant (Supplementary file 1a) for 7–10 days: The glass covering the wells was broken without disturbing the bolus and 2 μl soak solution (Supplementary file 1a) were added before resealing with a fresh siliconized glass coverslip, incubating overnight (16–18 hr), then harvesting and flash freezing for data collection.
X-ray diffraction data collection and processing
Request a detailed protocolDiffraction data for structure determination and refinement were collected at beamlines 24-ID-C or 23-ID-B of the Advanced Photon Source at wavelengths of 0.984 Å or 1.033 Å, respectively. We used anomalous signals to confirm the presence of metals (Mn2+ or Cd2+) in the binding site. For D296A•Mn2+ and A47W•Mn2+, the same data we used for structure refinement and collected at 1.033 Å, 0.984 Å, respectively, provided strong anomalous signal in the metal-binding sites. For WT•Cd2+, we were able to collect data at 1.904 Å (near the low-energy boundary for the beamline), to maximize the anomalous signal. Locations of the crystals in the mesh loops were identified by grid scanning with a 20 μm beam at 10% transmission followed by data collection with a 10 μm beam at 15% transmission. Data were indexed in XDS (Kabsch, 2010) and scaled in CCP4 AIMLESS (Version 7.0) (Evans and Murshudov, 2013; Winn et al., 2011). For datasets collected from several crystals, data from each crystal were independently indexed and integrated, then combined during scaling using CCP4 AIMLESS (Evans and Murshudov, 2013; Winn et al., 2011) to obtain complete datasets. The resolution cut-off of each structure was defined based on CC1/2 values of 0.3 and above (Karplus and Diederichs, 2012; Karplus and Diederichs, 2015). Initial phases for all structures were determined by molecular replacement in PHENIX (Version 1.17.1–3660) (Liebschner et al., 2019) using an occluded structure of DraNramp (PDB ID 6C3I chain A) (Bozzi et al., 2019b) as search model. Data statistics are listed in Table 1 and Supplementary file 1b and d.
Model building, refinement, and analysis
Request a detailed protocolModels were built in COOT (Version 7.0) (Emsley et al., 2010) and refined in PHENIX (Liebschner et al., 2019), with macrocycles including reciprocal space, TLS groups, and individual B-factor refinement, and optimization of the X-ray/stereochemistry and X-ray/ADP weights. For WT•Cd2+, ‘anomalous group refinement’ was used to improve the fit to density of the Cd2+ ions, with Cd2+ as an ‘anomalous group’ with the reference f’ and f” values suggested by phenix.form.factor (–0.462 and 2.132, respectively). Ligand restraints for monoolein and spermidine were generated in Phenix.elbow with automatic geometry optimization. All structures contain one protein molecule in the asymmetric unit. The final structures span from residues 45–48 to residues 433–436, except that residues 240–249 and 240–247 were not modeled in D296•Mn2+ and WT•Cd2+, respectively, because of lack of interpretable electron density map. Model refinement statistics are listed in Table 1 and Supplementary file 1b and d. Pairwise RMSD for all structures are listed in Supplementary file 1c. Anomalous difference Fourier maps for Mn2+ (D296A•Mn2+ and A47W•Mn2+) and Cd2+ (WT•Cd2+) were generated in phenix.maps using a high-resolution cutoff of 3.5–4.5 Å. Polder maps omitting the metal ions were generated in phenix.polder to appropriately define the coordination sphere for Mn2+ and Cd2+ in all structures. All software were provided by SBGrid (Morin et al., 2013).
Metal binding measurements using ITC
Request a detailed protocolITC experiments were performed using MicroCal iTC200 (GE Healthcare) to determine the affinity and thermodynamic parameters of binding of divalent metals (Mn2+, Cd2+ and Mg2+) to DraNramp (WT and its mutants) (Leavitt and Freire, 2001; Wiseman et al., 1989). All protein and metal solutions were prepared in ITC buffer (150 mM NaCl, 10 mM HEPES, pH 7.5, and 0.003% LMNG). The sample cell containing 25 μM protein was titrated with 20 2 μL injections of 6 mM metal, with an interval of 120 s between each successive 5 s injection, with a 750 rpm stirring rate at 25 °C. To nullify the heat of dilution, the data from titration of a metal solution into ITC buffer (‘buffer blank’ runs) were subtracted from the metal-protein titration curves prior to model fitting. Data were fitted and analyzed as detailed in Appendix 1. Briefly, as per best practice when c values (association constant × molar protein concentration) are below 1 (Picollo et al., 2009; Tellinghuisen, 2008; Turnbull and Daranas, 2003), all the data reported for each construct are fitted fixing the number of sites (one-site binding model with fixed n=1 or sequential binding model with fixed n=2). The binding stoichiometry was selected based the model fits and knowledge from the crystal structures and mutational analysis. Data were fitted with Origin 7 software. The mean Kd values ± SEM from 2 to 3 repeats (from independent protein purifications) for each sample are reported in Supplementary file 1g.
Proteoliposome-based in vitro transport assays
Request a detailed protocolProtein purification, liposome preparation, and metal transport assays were performed as described (Bozzi et al., 2016a; Bozzi et al., 2019a). Purifications of each protein construct were performed at least twice and resulted in similar data. For each construct, the presented data originates from two to three batches of reconstituted liposomes, each from an independent protein purification, represented as scatter plots of each data point and the corresponding mean.
Sequence alignments
Request a detailed protocolWe used 92 Nramp sequences from Pfam (El-Gebali et al., 2019) to build a seed alignment using MUSCLE (Edgar, 2004). We collected 15,451 sequences from Uniprot (Bateman et al., 2021) using HMMER (Potter et al., 2018). We used HMMER’s hmmalign, with a hidden Markov model profile from the seed alignment as an input, to align all 15,451 sequences. We applied filters to retain sequences 400–600 residues in length and sequences with under 90% pairwise sequence identity, respectively. The final alignment contains 6712 sequences and is well aligned at biologically relevant residues (Figure 1—source data 1). A maximum-likelihood phylogenetic tree was generated via RAxML-NG Kozlov et al., 2019 using the LG substitution model (Le and Gascuel, 2008), with the likeliest final tree selected from 10 parallel optimization trials (Figure 1—source data 2). The canonical Nramp clade in this tree was identified based on conservation of the ‘DPGN’ and ‘MPH’ motifs in transmembrane helices 1 and 6, respectively, and contained 3796 sequences. The Nramp-related magnesium transporters were used to root the canonical Nramp phylogeny. Sequence analysis was done with Biopython (Cock et al., 2009) and sequence logos were generated with logomaker (Tareen and Kinney, 2020) using a ‘chemistry’ color scheme inspired by WebLogo (Crooks et al., 2004).
Molecular dynamics simulation
Request a detailed protocolMolecular dynamics (MD) simulations were initialized from three high-resolution structures of DraNramp: the outward-open G223W•Mn2+ structure (6BU5) with Mn2+ removed and W223 mutated back to the native glycine residue in silico, the inward-open WT•Cd2+ structure with Cd2+ removed, and the inward-occluded WT structure. Crystallographic waters were retained, and protonation states of key titratable residues were selected with PROPKA (Olsson et al., 2011; Søndergaard et al., 2011) assuming a pH of 5.0 for residues exposed to external solvent and a pH of 7.0 for residues exposed to cytosol, a condition under which DraNramp exhibits high activity. All structures were oriented in the membrane with the PPM web server and membrane systems were prepared with CHARMM-GUI (Jo et al., 2008; Lee et al., 2016). A POPC membrane of surface area 99×99 Å was constructed in the XY plane around the protein (Wu et al., 2014), the system was solvated in a 100×100×100 Å3 rectangular box using TIP3 waters and electronically neutralized using potassium and chlorine ions at an overall concentration of 150 mM. The overall system size was approximately 103,000 atoms.
All-atom simulations were run using GPU-accelerated NAMD (Phillips et al., 2008) and the CHARMM36m forcefield (Huang et al., 2017). Prior to simulation, the energy of each system was minimized for 10,000 steps using a conjugate gradient and line search algorithm native to NAMD. To improve simulation stability, the system was initially equilibrated using an NVT-ensemble with harmonic restraints placed on protein and lipid heavy atoms. The harmonic restraints were then incrementally relaxed over a period of 675 ps according to established CHARMM-GUI protocols (Lee et al., 2016). The system was then simulated at a constant pressure, utilizing the Langevin piston method to maintain 1 atm at 303.15 K, from anywhere between 617 and 1176 ns depending on the starting conformation. Simulations were performed using periodic boundary conditions and a time step of 3.0 fs with all bonds to hydrogens being constrained. Large integration timesteps were enabled by employing hydrogen mass repartitioning (Hopkins et al., 2015). Long-range electrostatic interactions were calculated using the particle mesh Ewald (PME) method with nonbonded interactions being cut off at 12 Å. Each simulation was performed in duplicate resulting in approximately 2 µs of total sampling for each system. Simulations are summarized in Supplementary file 1i.
RMSD, residue distance, and dihedral analyses were performed using mdtraj version 1.9.8 (McGibbon et al., 2015). For distance analysis, the minimum interatomic distances were identified between the specified residues across each frame. For water analysis, simulations were centered and wrapped in VMD and a Tcl script was used to produce water density maps. Contour maps of these densities were then visualized in PyMOL. Water occupancies of sites coordinated by specific residues were also calculated in an alignment-agnostic manner by determining for each frame in each simulation whether a water was present within 2.5 Å of both specified residues. Distinct rotamers were identified from dihedrals using spectral clustering as implemented in scikit-learn version 1.0.
Data availability
Request a detailed protocolAtomic coordinates and structure factors for the crystal structures reported in this work have been deposited to the Protein Data Bank under accession numbers 8E5S (WT), 8E5V (WTsoak), 8E60 (WT•Mn2+), 8E6H (A47W•Mn2+), 8E6I (M230A•Mn2+), 8E6L (D296A•Mn2+), 8E6M (WT•Cd2+), and 8E6N (re-refined G223W•Mn2+). Corresponding X-ray diffraction images have been deposited to the SBGrid Data Bank under the respective accession numbers 962 (doi:10.15785/SBGRID/962), 963 (doi:10.15785/SBGRID/963), 964 (doi:10.15785/SBGRID/ 964), 966 (doi:10.15785/SBGRID/966), 967 (doi:10.15785/SBGRID/967), 968 (doi:10.15785/SBGRID/968), 969 (doi:10.15785/SBGRID/969), and previously deposited 564 (doi:10.15785/SBGRID/564). The multiple sequence alignment and phylogenetic tree have been provided as Figure 1—source data 1 and Figure 1—source data 2, respectively. All liposome-based transport data are provided in Figure 1—source data 3. Code for analysis of molecular dynamics data, as well as the raw data plotted in Figure 4—figure supplement 2 and Figure 4—figure supplement 3, can be found at https://github.com/samberry19/nramp-md (MIT license). Raw molecular dynamics trajectory files are available on Dryad (https://doi.org/10.5061/dryad.tx95x6b2b). Source files (origin files) of all ITC experiments are provided in Appendix 1—table 1—source data 1 (Mn2+ isotherms), Appendix 1—table 2—source data 1 (Cd2+ isotherms) and Figure 5—source data 1(Mg2+ isotherms).
Appendix 1
ITC analysis
In this appendix, we describe our choices of models to fit our ITC isotherms. Based on the isotherm shapes and preliminary fits, the c values (association constant × molar protein concentration) for all isotherms are below 1. Curve-fitting is most robust for 5<c < 500, a range which generally allows to fit binding site stoichiometry (n) in addition to the dissociation constant, Kd, (or association constant, Ka) and the enthalpy of association (ΔH) (Turnbull and Daranas, 2003). The combination of low metal ion binding affinity and the achievable protein amounts and concentration for DraNramp makes it impractical to achieve a higher c value. However, previous studies have demonstrated that data with c<1 can still provide useful information on Kd values (Picollo et al., 2009; Tellinghuisen, 2008; Turnbull and Daranas, 2003), Hence, as recommended in those studies, we fit the isotherms with either a one-site model with a fixed n=1, or a two-site sequential binding model (fixed n=2; Appendix 1—tables 1 and 2 and Appendix 1—figure 1). The assumptions inherent to fixing n are that we have accurately determined the concentrations of metal ion ligand and protein, but fixing n has little impact on the resulting Kd value (Tellinghuisen, 2008). We then used the quality of fit and information from our crystal structures, mutational data to select the most appropriate fit to interpret each ITC isotherm dataset. Below we justify our choice of model for pairs of protein construct and metal ion ligand:
-
Appendix 1—table 1—source data 1
- https://cdn.elifesciences.org/articles/84006/elife-84006-app1-table1-data1-v2.zip
-
Appendix 1—table 2—source data 1
- https://cdn.elifesciences.org/articles/84006/elife-84006-app1-table2-data1-v2.zip
WT binding to Mn2+
The WT•Mn2+ crystal structure revealed two Mn2+ ions bound at distinct sites, and the overall fit of the ITC data is better for the two-site sequential binding model than the one-site model (Appendix 1—figure 1A), with the two Kd values differing by ~20 fold. Based on ITC of the constructs with mutations of the external site (D296A and D236A; see below), we assigned the higher affinity Kd1 to the orthosteric site (190±30 µM) and the lower affinity Kd2 to the external site (1970±520 µM).
A47W binding to Mn2+
The A47W•Mn2+ crystal structure revealed two Mn2+ ions bound to the same sites as with WT•Mn2+, and the overall fit of the ITC data is better for the two-site sequential binding model than for the one-site model (Appendix 1—figure 1B). The two Kd values and the ITC data for the constructs with external-site mutations (A47W-D296A and A47W-D236A, both of which can only be fitted with a one-site model) are consistent with the Kd assignments for WT, with the orthosteric site having higher affinity for Mn2+ (Kd1=125±5 µM) than the external site (Kd2=2450±650 µM).
M230A binding to Mn2+
The M230A•Mn2+ crystal structure again revealed two bound Mn2+, and the overall fit of the ITC data is again better for the two-site sequential binding model than the one-site model (Appendix 1—figure 1C). The two Kd values and the ITC data for the constructs with external-site mutations (M230A-D296A and M230A-D236A, both of which only fit with a one-site model) are consistent with the Kd assignments for WT, with the orthosteric site having higher affinity for Mn2+ (Kd1=215±65 µM) than the external site (Kd2=4600±1300 µM).
D56A binding to Mn2+
The ITC data can be fitted using either the two-site sequential model or the one-site model (Appendix 1—figure 1D). In contrast, the data for constructs with external-site mutations (D56A-D296A and D56A-D369A) only fit well using a one-site model. Hence, we selected the two-site sequential binding model to fit D56A binding to Mn2+, as it agrees best with the overall analysis, with the orthosteric site having higher affinity for Mn2+ (Kd1=230±80 µM) than the external site (Kd2=3800±1100 µM).
WT binding to Cd2+
The WT•Cd2+ crystal structure reveals two Cd2+ ions bound at the same two distinct sites where Mn2+ bound. Based on this observation we chose a two-site model rather than a one-site model to fit the ITC data, although the fits look similar (Appendix 1—figure 1E). In contrast to Mn2+ binding, the two Kd values obtained using a two-site sequential binding model are more similar to each other (~4 fold difference). Based on ITC of the constructs with external-site mutations (D296A and D236A; see below), we assigned the higher affinity Kd (Kd1=55±15 µM) to the orthosteric site, and the lower affinity one to the external site (Kd2=220±20 µM).
A47W binding to Cd2+
One-site model-based fits of the ITC data consistently yield better thermodynamic binding parameters; the two-site sequential binding model either does not fit or only fits with a positive ΔH which is inappropriate considering the observed exothermic binding. Furthermore, ITC measurements on the corresponding constructs with external-site mutations (A47W-D296A and A47W-D369A) showed no binding. Therefore, we chose the one-site model to fit the ITC data of A47W binding to Cd2+ and we assigned the resulting Kd value (150±10 µM) to the external site. We conclude that Cd2+ does not bind to the A47W orthosteric site with measurable affinity.
M230A binding to Cd2+
The ITC data could only be fitted with the one-site model. The ITC measurements on the corresponding constructs with external-site mutations (M230A-D296A and M230A-D369A) showed no binding. We thus assigned the Kd value obtained using the one-site model (160±20 µM) to the external site, and we conclude that Cd2+ does not bind to the M230A orthosteric site with measurable affinity.
D56A binding to Cd2+
The ITC data fit better when using a two-site sequential binding model than a one-site model (Appendix 1—figure 1F). The WT•Cd2+ structure shows no direct interaction of D56 with Cd2+, suggesting that the D56A mutation may not completely impair Cd2+ binding at the orthosteric site. We thus chose to use the two-site sequential binding model and assign higher affinity (Kd1=80±20 µM) to the external site, and lower affinity (Kd2=2250±750 µM) to the orthosteric site.
G223W binding to Mn2+
Consistent with the G223W•Mn2+ crystal structure in which Mn2+ is only bound at the orthosteric site, the ITC data could only be fitted with the one-site model. We assigned the resulting Kd value (440±15 µM) to the orthosteric site.
G223W binding to Cd2+
The ITC data show no binding, which is consistent with our inability to obtain a crystal structure of G223W bound to Cd2+. This indicates that the outward-locked G223W does not bind Cd2+ with measurable affinity at the orthosteric site.
D296A binding to Mn2+
Consistent with the D296A•Mn2+ crystal structure in which Mn2+ is bound at the orthosteric site but not at the external site, the ITC data could only be fitted with the one-site model. We assigned the resulting Kd value (370±30 µM) to the orthosteric site.
D296A binding to Cd2+
The ITC data could only be fitted with the one-site model. We assigned the resulting Kd value (120±10 µM) to the orthosteric site, because the D296A mutation removes an aspartate ligand at the external site.
D369A binding to Mn2+
The ITC data could only be fitted with the one-site model. We assigned the resulting Kd value (420±30 µM) to the orthosteric site, because the D369A mutation removes an aspartate ligand at the external site.
D369A binding to Cd2+
Fits of the ITC data using the one-site model consistently yielded better thermodynamic binding parameters; the two-site sequential binding model either resulted in no fit or in a fit with positive ΔH which is inappropriate considering the observed exothermic binding. We thus chose to fit the data with the one-site model, and we assigned the resulting Kd value (70±10 µM) to the orthosteric site.
Data availability
Atomic coordinates and structure factors for the crystal structures reported in this work have been deposited to the Protein Data Bank under accession numbers 8E5S (WT), 8E5V (WTsoak), 8E60 (WT•Mn2+), 8E6H (A47W•Mn2+), 8E6I (M230A•Mn2+), 8E6L (D296A•Mn2+), 8E6M (WT•Cd2+), and 8E6N (re-refined G223W•Mn2+). Corresponding X-ray diffraction images have been deposited to the SBGrid Data Bank under the respective accession numbers 962 (doi:10.15785/SBGRID/962), 963 (doi:10.15785/SBGRID/963), 964 (doi:10.15785/SBGRID/ 964), 966 (doi:10.15785/SBGRID/966), 967 (doi:10.15785/SBGRID/967), 968 (doi:10.15785/ SBGRID/968), 969 (doi:10.15785/SBGRID/969), and previously deposited 564 (doi:10.15785/ SBGRID/564). The multiple sequence alignment and phylogenetic tree have been provided as Figure 1-source data 1 and Figure 1-source data 2, respectively. All liposome-based transport data are provided in Figure 1-source data 3. Code for analysis of molecular dynamics data, as well as the raw data plotted in Figure 4-figure supplement 2 and Figure 4-figure supplement 3, can be found at https://github.com/samberry19/nramp-md (MIT license). Raw molecular dynamics trajectory files are available on Dryad (https://doi.org/10.5061/dryad.tx95x6b2b). Source files (origin files) of all ITC experiments are provided in Appendix 1-table 1-source data 1 (Mn2+ isotherms), Appendix 1-table 2-source data 1 (Cd2+ isotherms) and Figure 5-source data 1 (Mg2+ isotherms).
-
Dryad Digital RepositoryMolecular dynamics simulations in: High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter.https://doi.org/10.5061/dryad.tx95x6b2b
-
RCSB Protein Data BankID 8E5S. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an occluded state.
-
RCSB Protein Data BankID 8E5V. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WTsoak in an occluded state.
-
RCSB Protein Data BankID 8E60. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an occluded, manganese-bound state.
-
RCSB Protein Data BankID 8E6H. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter A47W mutant in an occluded, manganese-bound state.
-
RCSB Protein Data BankID 8E6I. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter M230A mutant in an inward-open, manganese-bound state.
-
RCSB Protein Data BankID 8E6L. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter D296A mutant in an inward-open, manganese-bound state.
-
RCSB Protein Data BankID 8E6M. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an inward-open, cadmium-bound state.
-
RCSB Protein Data BankID 8E6N. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter G223W mutant in an outward-open, manganese-bound state.
-
SBGrid Data BankX-Ray Diffraction data from WT Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E5S structure.https://doi.org/10.15785/SBGRID/962
-
SBGrid Data BankX-Ray Diffraction data from WTsoak Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E5V structure.https://doi.org/10.15785/SBGRID/963
-
SBGrid Data BankX-Ray Diffraction data from Manganese-bound WT Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E60 structure.https://doi.org/10.15785/SBGRID/964
-
SBGrid Data BankX-Ray Diffraction data from Mn-bound A47W mutant Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E6H structure.https://doi.org/10.15785/SBGRID/966
-
SBGrid Data BankX-Ray Diffraction data from Mn-bound M230A mutant Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E6I structure.https://doi.org/10.15785/SBGRID/967
-
SBGrid Data BankX-Ray Diffraction data from Mn-bound D296A mutant Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E6L structure.https://doi.org/10.15785/SBGRID/968
-
SBGrid Data BankX-Ray Diffraction data from Cadmium-bound WT Nramp/MntH divalent transition metal transporter from Deinococcus radiodurans, source of 8E6M structure.https://doi.org/10.15785/SBGRID/969
-
SBGrid Data BankX-Ray Diffraction data from Deinococcus radiodurans Nramp/MntH divalent transition metal transporter in the outward-open, mangan, source of 6BU5 structure.https://doi.org/10.15785/SBGRID/564
References
-
Metal transporters and diseaseCurrent Opinion in Chemical Biology 6:181–186.https://doi.org/10.1016/s1367-5931(02)00307-1
-
UniProt: the universal protein knowledgebase in 2021Nucleic Acids Research 49:D480–D489.https://doi.org/10.1093/nar/gkaa1100
-
Unique structural features in an nramp metal transporter impart substrate-specific proton cotransport and a kinetic bias to favor importThe Journal of General Physiology 151:1413–1429.https://doi.org/10.1085/jgp.201912428
-
Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an nramp-family transition metal transporterThe Journal of Biological Chemistry 295:1212–1224.https://doi.org/10.1074/jbc.RA119.011336
-
Molecular mechanism of nramp-family transition metal transportJournal of Molecular Biology 433:166991.https://doi.org/10.1016/j.jmb.2021.166991
-
The role of manganese dysregulation in neurological disease: emerging evidenceThe Lancet. Neurology 20:956–968.https://doi.org/10.1016/S1474-4422(21)00238-6
-
ConferenceConsiderations for the Specification of Enzyme Assays Involving Metal IonsProceedings of the 3rd Beilstein ESCEC Symposium “Experimental Standard Conditions Of Enzyme Characterization.”.
-
Selective recognition of metal ions by metalloregulatory proteinsCurrent Opinion in Chemical Biology 12:214–221.https://doi.org/10.1016/j.cbpa.2007.12.010
-
Imperfect coordination chemistry facilitates metal ion release in the PSA permeaseNature Chemical Biology 10:35–41.https://doi.org/10.1038/nchembio.1382
-
Principles of alternating access in leut-fold transporters: commonalities and divergencesJournal of Molecular Biology 434:167746.https://doi.org/10.1016/j.jmb.2022.167746
-
Factors governing the metal coordination number in metal complexes from cambridge structural database analysesThe Journal of Physical Chemistry. B 110:1889–1895.https://doi.org/10.1021/jp054975n
-
Crystal structure of a SLC11 (NRAMP) transporter reveals the basis for transition-metal ion transportNature Structural & Molecular Biology 21:990–996.https://doi.org/10.1038/nsmb.2904
-
The effect of iron deficiency anemia on the function of the immune systemThe Hematology Journal 5:579–583.https://doi.org/10.1038/sj.thj.6200574
-
The pfam protein families database in 2019Nucleic Acids Research 47:D427–D432.https://doi.org/10.1093/nar/gky995
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
How good are my data and what is the resolution?Acta Crystallographica. Section D, Biological Crystallography 69:1204–1214.https://doi.org/10.1107/S0907444913000061
-
Structure of the manganese-bound manganese transport regulator of Bacillus subtilisNature Structural Biology 10:652–657.https://doi.org/10.1038/nsb951
-
The geometry of metal-ligand interactions relevant to proteins. II. angles at the metal atom, additional weak metal-donor interactionsActa Crystallographica. Section D, Biological Crystallography 56:857–867.https://doi.org/10.1107/s0907444900005849
-
Geometry of metal-ligand interactions in proteinsActa Crystallographica. Section D, Biological Crystallography 57:401–411.https://doi.org/10.1107/s0907444900019168
-
Long-time-step molecular dynamics through hydrogen mass repartitioningJournal of Chemical Theory and Computation 11:1864–1874.https://doi.org/10.1021/ct5010406
-
Uptake of lead and iron by divalent metal transporter 1 in yeast and mammalian cellsBiochemical and Biophysical Research Communications 295:978–984.https://doi.org/10.1016/s0006-291x(02)00756-8
-
Substrate profile and metal-ion selectivity of human divalent metal-ion transporter-1The Journal of Biological Chemistry 287:30485–30496.https://doi.org/10.1074/jbc.M112.364208
-
637. the stability of transition-metal complexesJournal of the Chemical Society 3192.https://doi.org/10.1039/jr9530003192
-
CHARMM-GUI: a web-based graphical user interface for CHARMMJournal of Computational Chemistry 29:1859–1865.https://doi.org/10.1002/jcc.20945
-
Improvements to the APBS biomolecular solvation software suiteProtein Science 27:112–128.https://doi.org/10.1002/pro.3280
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
Assessing and maximizing data quality in macromolecular crystallographyCurrent Opinion in Structural Biology 34:60–68.https://doi.org/10.1016/j.sbi.2015.07.003
-
Transferrin and transferrin receptors updateFree Radical Biology & Medicine 133:46–54.https://doi.org/10.1016/j.freeradbiomed.2018.06.037
-
Characterization of the metal ion binding site in the anti-terminator protein, hutp, of Bacillus subtilisNucleic Acids Research 33:5494–5502.https://doi.org/10.1093/nar/gki868
-
An improved general amino acid replacement matrixMolecular Biology and Evolution 25:1307–1320.https://doi.org/10.1093/molbev/msn067
-
Direct measurement of protein binding energetics by isothermal titration calorimetryCurrent Opinion in Structural Biology 11:560–566.https://doi.org/10.1016/s0959-440x(00)00248-7
-
CHARMM-gui input generator for namd, gromacs, amber, openmm, and CHARMM/openmm simulations using the charmm36 additive force fieldJournal of Chemical Theory and Computation 12:405–413.https://doi.org/10.1021/acs.jctc.5b00935
-
Macromolecular structure determination using x-rays, neutrons and electrons: recent developments in phenixActa Crystallographica. Section D, Structural Biology 75:861–877.https://doi.org/10.1107/S2059798319011471
-
Metals and membrane metal transporters in biological systems: the role (S) of Nramp in host-parasite interactionsJ Argent Chem Soc 97:210–225.
-
Iron and manganese transport in mammalian systemsBiochimica et Biophysica Acta. Molecular Cell Research 1868:118890.https://doi.org/10.1016/j.bbamcr.2020.118890
-
MDTraj: A modern open library for the analysis of molecular dynamics trajectoriesBiophysical Journal 109:1528–1532.https://doi.org/10.1016/j.bpj.2015.08.015
-
The NRAMP family of metal-ion transportersBiochimica et Biophysica Acta 1763:609–620.https://doi.org/10.1016/j.bbamcr.2006.05.007
-
The thermodynamics of protein-ligand interaction and solvation: insights for ligand designJournal of Molecular Biology 384:1002–1017.https://doi.org/10.1016/j.jmb.2008.09.073
-
PROPKA3: consistent treatment of internal and surface residues in empirical pka predictionsJournal of Chemical Theory and Computation 7:525–537.https://doi.org/10.1021/ct100578z
-
ConferenceAdapting a message-driven parallel application to GPU-accelerated clustersSC - International Conference for High Performance Computing, Networking, Storage and Analysis.https://doi.org/10.1109/SC.2008.5214716
-
Basis of substrate binding and conservation of selectivity in the CLC family of channels and transportersNature Structural & Molecular Biology 16:1294–1301.https://doi.org/10.1038/nsmb.1704
-
HMMER web server: 2018 updateNucleic Acids Research 46:W200–W204.https://doi.org/10.1093/nar/gky448
-
Coordination geometries of selected transition metal ions (Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+) in metalloproteinsJournal of Inorganic Biochemistry 71:115–127.https://doi.org/10.1016/s0162-0134(98)10042-9
-
Properties of the mammalian and yeast metal-ion transporters DCT1 and Smf1p expressed in Xenopus laevis oocytesThe Journal of Experimental Biology 204:1053–1061.https://doi.org/10.1242/jeb.204.6.1053
-
Molecular mechanism of ligand recognition by membrane transport protein, mhp1The EMBO Journal 33:1831–1844.https://doi.org/10.15252/embj.201387557
-
Stability and Applications of CoordinationStability and applications of coordination compounds, Stability and Applications of Coordination, InTechOpen, 10.5772/intechopen.83186.
-
Infection genomics: nramp1 as a major determinant of natural resistance to intracellular infectionsAnnual Review of Medicine 49:275–287.https://doi.org/10.1146/annurev.med.49.1.275
-
Improved treatment of ligands and coupling effects in empirical calculation and rationalization of pKa valuesJournal of Chemical Theory and Computation 7:2284–2295.https://doi.org/10.1021/ct200133y
-
Logomaker: beautiful sequence logos in pythonBioinformatics 36:2272–2274.https://doi.org/10.1093/bioinformatics/btz921
-
Isothermal titration calorimetry at very low cAnalytical Biochemistry 373:395–397.https://doi.org/10.1016/j.ab.2007.08.039
-
On the value of C: can low affinity systems be studied by isothermal titration calorimetry?Journal of the American Chemical Society 125:14859–14866.https://doi.org/10.1021/ja036166s
-
Different divalent cations alter the kinetics and fidelity of DNA polymerasesThe Journal of Biological Chemistry 291:20869–20875.https://doi.org/10.1074/jbc.R116.742494
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
Rapid measurement of binding constants and heats of binding using a new titration calorimeterAnalytical Biochemistry 179:131–137.https://doi.org/10.1016/0003-2697(89)90213-3
-
CHARMM-GUI membrane builder toward realistic biological membrane simulationsJournal of Computational Chemistry 35:1997–2004.https://doi.org/10.1002/jcc.23702
-
Crystal structure of Helicobacter pylori neutrophil-activating protein with a di-nuclear ferroxidase center in a zinc or cadmium-bound formBiochemical and Biophysical Research Communications 422:745–750.https://doi.org/10.1016/j.bbrc.2012.05.073
Article and author information
Author details
Funding
National Institute of General Medical Sciences (R01GM120996)
- Rachelle Gaudet
National Science Foundation (MCB-1942763)
- Abhishek Singharoy
National Science Foundation (1764269)
- Samuel P Berry
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Arghya Deb for discussions of metal coordination chemistry, and previous and current members of the Gaudet lab for discussions and assistance, particularly Gerardo Zavala and Edward Lee for contributions to preliminary molecular dynamics analyses, and José Velilla for help with crystal soaking and fishing prior to data collection. This work was funded by NIGMS grant R01GM120996 (RG), a National Science Foundation (NSF) CAREER award MCB-1942763 (AS), AstraZenaca and ASU-Mayo Foundation (EW), and the NSF-Simons Center for Mathematical and Statistical Analysis of Biology at Harvard (award number 1764269) and the Harvard Quantitative Biology Initiative (SB). Diffraction data reported in this study were collected at NE-CAT beamline 24IDC and GM/CA beamline 23IDB in the Advanced Photon Source. NE-CAT is funded by NIGMS grant P30 GM124165 and GM/CA is funded by the National Cancer Institute (ACB-12002) and the National Institute of General Medical Sciences (AGM-12006, P30GM138396). The Eiger 16 M detector at GM/CA-XSD is funded by NIH grant S10 OD012289. The Advanced Photon Source is a U.S. Department of Energy Facility operated by Argonne National Laboratory under Contract No. DE-AC02-06CH11357. The molecular simulations used the Extreme Science and Engineering Discovery Environment (XSEDE) supported by NSF (ACI-1548562), and Oak Ridge Leadership Computing Facility, supported by the Office of Science, Department of Energy (DE-AC05-00OR22725).
Copyright
© 2023, Ray 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.
Metrics
-
- 1,621
- views
-
- 262
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
Polycystin-1 (PC1) is the protein product of the PKD1 gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Lipopolysaccharides (LPS) confer resistance against harsh conditions, including antibiotics, in Gram-negative bacteria. The lipopolysaccharide transport (Lpt) complex, consisting of seven proteins (A-G), exports LPS across the cellular envelope. LptB2FG forms an ATP-binding cassette transporter that transfers LPS to LptC. How LptB2FG couples ATP binding and hydrolysis with LPS transport to LptC remains unclear. We observed the conformational heterogeneity of LptB2FG and LptB2FGC in micelles and/or proteoliposomes using pulsed dipolar electron spin resonance spectroscopy. Additionally, we monitored LPS binding and release using laser-induced liquid bead ion desorption mass spectrometry. The β-jellyroll domain of LptF stably interacts with the LptG and LptC β-jellyrolls in both the apo and vanadate-trapped states. ATP binding at the cytoplasmic side is allosterically coupled to the selective opening of the periplasmic LptF β-jellyroll domain. In LptB2FG, ATP binding closes the nucleotide binding domains, causing a collapse of the first lateral gate as observed in structures. However, the second lateral gate, which forms the putative entry site for LPS, exhibits a heterogeneous conformation. LptC binding limits the flexibility of this gate to two conformations, likely representing the helix of LptC as either released from or inserted into the transmembrane domains. Our results reveal the regulation of the LPS entry gate through the dynamic behavior of the LptC transmembrane helix, while its β-jellyroll domain is anchored in the periplasm. This, combined with long-range ATP-dependent allosteric gating of the LptF β-jellyroll domain, may ensure efficient and unidirectional transport of LPS across the periplasm.