High-resolution structures with bound Mn2+ and Cd2+ map the metal import pathway in an Nramp transporter

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Appendix 1
Data availability
References
Article and author information
Metrics

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, Mn²⁺, across conformations, and by conserved networks of polar residues lining the inner and outer gates. In addition, a high-resolution Cd²⁺-bound structure highlights differences in how Cd²⁺ and Mn²⁺ 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 Mn²⁺ is different and more robust to perturbation than for transporting the toxic Cd²⁺ 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.sa0

Introduction

Transition metal ions like Mn²⁺ and Fe²⁺ 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, Fe²⁺ deficiency causes anemia and neurodegenerative diseases, whereas Fe²⁺ overload increases the risk of cancer by generating toxic reactive oxygen species (ROS) and mutations (Ekiz et al., 2005; Jung et al., 2015). Mn²⁺ 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 Cd²⁺ and Hg²⁺, 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 Fe²⁺ and Mn²⁺ 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 Mn²⁺ and Fe²⁺ into the cytosol (Pujol-Giménez et al., 2017; Skamene et al., 1998; Vidal et al., 1993). Plant and fungal Nramps aid in Fe²⁺ and Mn²⁺ uptake and trafficking, and bacterial Nramps are involved in the acquisition of Mn²⁺, an essential nutrient (Bozzi and Gaudet, 2021). In addition to the physiological substrates Fe²⁺ and Mn²⁺, Nramps can also transport toxic metals like Cd²⁺ and Hg²⁺ but exclude the abundant alkaline earth metals like Mg²⁺ and Ca²⁺ (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 Mg²⁺ transporter), can transport Mg²⁺ (Ramanadane et al., 2022). However, little is known about whether the canonical Nramps can mechanistically distinguish between their physiological substrates (Fe²⁺ and Mn²⁺) from non-essential ones like Cd²⁺ within their broad spectrum of transition metal substrates. Functional studies on Deinococcus radiodurans (Dra)Nramp revealed that Mn²⁺ and Cd²⁺ 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 Mn²⁺-bound and metal-free states, providing the first molecular map of the entire Mn²⁺ transport cycle. The structures along with molecular simulations reveal that Nramps achieve alternate access during transport by adopting distinct Mn²⁺-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 Mn²⁺ transport cycle also informs on the transport of Fe²⁺, the other common physiological Nramp substrate, because Fe²⁺ and Mn²⁺ 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, Cd²⁺, 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 Mn²⁺ 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 Mn²⁺ at 2.38 Å by soaking crystals with Mn²⁺ (WT•Mn²⁺; Table 1, Figure 1A–B). We resolved a comparable structure using co-crystallization with Mn²⁺ and the inward-locking mutation A47W (A47W•Mn²⁺; 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 Mn²⁺-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•Mn²⁺ and A47W•Mn²⁺ differ in their TM1a position, with WT•Mn²⁺ nearly identical to G45R whereas the A47W•Mn²⁺ TM1a is displaced within the inner vestibule, likely to accommodate the bulky tryptophan sidechain. Therefore, we generally used the WT•Mn²⁺ structure for analysis of the occluded state. As in the metal-free G45R, the Mn²⁺-bound occluded structures have a completely sealed outer vestibule and a partially closed inner vestibule, with the Mn²⁺ occluded from bulk solvent (Figure 1A).

Figure 1 with 4 supplements see all

Download asset Open asset

The occluded structure of DraNramp reveals a largely dehydrated Mn²⁺-coordination sphere.

(A) Cartoon representation of WT•Mn²⁺ in an occluded state. Anomalous signal confirmed the presence of Mn²⁺ in both the orthosteric metal-binding site and an additional site at the mouth of the external vestibule (Figure 1—figure supplement 1A) which is less conserved across the Nramp family (Figure 1—figure supplement 3). TM1 and TM6 are labeled. (B) Detail of the orthosteric metal-binding site of WT•Mn²⁺ where D56, N59, M230, and the pseudo-symmetrically related carbonyls of A53 and A227 coordinate the Mn²⁺ ion (Figure 1—figure supplement 1B). A water molecule completes the six-ligand coordination sphere. Coordinating residues are shown as sticks, and coordinating distances are indicated in Å. (C) ITC measurement of the affinity of WT DraNramp for Mn²⁺. Top graph shows heat absorbed upon injection of Mn²⁺ solution to the protein solution. Bottom graph shows the fit of the integrated and corrected heat to a binding isotherm. The data show an endothermic mode of binding and fits best with a two-site sequential binding model. The figure shows one of three measurements and the average K_d values ± SEM (K_d1=190±30 µM, K_d2=1970±520 µM; see Appendix 1). Based on ITC experiments comparing Mn²⁺ binding to WT or DraNramp constructs with mutations at the external site (Figure 1—figure supplement 2A), we assigned K_d1 to the orthosteric site. In all figures, unless otherwise noted, TMs 1, 5, 6, and 10 are pale yellow, TMs 2, 7, and 11 gray, TMs 3, 4, 8, and 9 light blue, and Mn²⁺ atoms are magenta spheres.

Figure 1—source data 1 Multiple sequence alignment of 6172 Nramp homologs.: https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data1-v2.zip
Download elife-84006-fig1-data1-v2.zip
Figure 1—source data 2 Maximum likelihood phylogenetic tree of Nramp homologs built with RAxML-NG.: https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data2-v2.zip
Download elife-84006-fig1-data2-v2.zip
Figure 1—source data 3 Raw data of metal ion uptake into proteoliposomes measured at four ΔΨ values for each DraNramp construct.: https://cdn.elifesciences.org/articles/84006/elife-84006-fig1-data3-v2.xlsx
Download elife-84006-fig1-data3-v2.xlsx

Table 1

Data collection and refinement statistics for four new DraNramp structures.

StructureConformationBound metal ion substratePDB ID	WT_soakOccluded none8E5V	WT•Mn²⁺OccludedMn²⁺8E60	M230A•Mn²⁺Inward openMn²⁺8E6I	WT•Cd²⁺Inward openCd²⁺8E6M
Data Collection
Beamline	GMCA 23IDB	GMCA 23IDB	GMCA 23IDB	NECAT 24IDC
Wavelength (Å)	1.033	1.033	1.033	0.984
Resolution range (Å)	41.23–2.36 (2.44–2.36)	41.28–2.38 (2.46–2.38)	45.32–2.52 (2.61–2.52)	45.54–2.48 (2.57–2.48)
Space group	P 2 21 21	P 2 21 21	P 2 21 21	P 2 21 21
Unit cell (a, b, c)	58.95, 71.04, 98.77	59.08, 71.10, 98.75	58.67, 71.35, 98.59	59.14, 71.37, 99.05
Unit cell (α, β, γ)	90, 90, 90	90, 90, 90	90, 90, 90	90, 90, 90
Number of crystals	1	1	3	1
Total reflections	58744 (5928)	57472 (5733)	146077 (14913)	76829 (7275)
Unique reflections	17477 (1718)	16468 (1646)	14548 (1427)	15351 (1507)
Redundancy	3.4 (3.4)	3.5 (3.5)	10.0 (10.4)	5.0 (4.8)
Completeness (%)	98.71 (98.85)	95.11 (96.92)	99.90 (99.79)	99.03 (99.47)
Mean I/σ (I)	8.89 (0.97)	8.92 (0.89)	8.47 (0.75)	9.90 (1.12)
Wilson B-factor	49.97	50.55	54.29	49.52
R_merge	0.106 (1.292)	0.109 (1.241)	0.269 (2.475)	0.158 (1.626)
R_meas	0.127 (1.511)	0.127 (1.447)	0.284 (2.603)	0.178 (1.831)
R_pim	0.067 (0.759)	0.063 (0.722)	0.090 (0.800)	0.078 (0.816)
CC1/2	0.99 (0.37)	0.99 (0.39)	0.98 (0.34)	0.99 (0.34)
Refinement
Resolution range (Å)	41.23–2.36 (2.44–2.36)	41.28–2.38 (2.46–2.38)	45.32–2.52 (2.61–2.52)	45.54–2.48 (2.57–2.48)
No. reflections	17441 (1714)	16438 (1636)	14547 (1425)	15291 (1507)
No. reflections in R_free	1743 (171)	1642 (164)	1454 (143)	1530 (151)
R_work	0.217 (0.340)	0.207 (0.316)	0.225 (0.313)	0.202 (0.319)
R_free	0.245 (0.350)	0.259 (0.358)	0.266 (0.349)	0.250 (0.354)
Number of atoms	3449	3385	3451	3321
Protein	2945	2933	2934	2905
Ligand	443	405	448	362
Water	61	47	69	54
Protein Residues	392	393	392	388
Ramachandran plot
Favored (%)	98.46	98.47	98.21	98.96
Allowed (%)	1.54	1.53	1.79	1.04
Outliers (%)	0	0	0	0
Rotamer outliers (%)	0.33	1.00	1.01	1.01
Clashscore	8.25	8.97	7.15	5.57
RMS (bonds)	0.002	0.002	0.002	0.002
RMS (angles)	0.43	0.46	0.43	0.46
Average B-factor	65.12	64.98	66.61	64.82
Protein	63.36	63.23	65.04	62.68
Ligand	77.99	78.48	77.70	83.11
Water	56.34	58.24	61.33	57.45
No. of TLS groups	9	8	3	3

Values in parentheses are for highest-resolution shell. Data for M230A•Mn²⁺ merge reflections from three crystals. Data for the other structures were obtained from a single crystal. See Supplementary file 1a for details on soaking or co-crystallization procedures for bound metal ion substrates.

Anomalous difference Fourier maps confirmed presence of Mn²⁺ 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•Mn²⁺, the Mn²⁺ 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 Mn²⁺ with Q378, a residue previously proposed to directly coordinate Mn²⁺ in the occluded state (Bozzi et al., 2019b). This yields a coordination number of 6, typical for Mn²⁺, and a largely dehydrated metal-binding site with a distorted octahedral Mn²⁺-coordination geometry (Supplementary file 1f), as often observed in other Mn²⁺-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 Mn²⁺ 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 (K_d1=190±30 µM, K_d2=1970±520 µM; Figure 1C) and A47W (K_d1=125±5 µM, K_d2=2450±650 µM; Figure 1—figure supplement 2A; all K_d 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 Mn²⁺ transport (Figure 1—figure supplement 2B). The D296A and D369A variants each bind one Mn²⁺ with K_d=370±30 µM and 420±30 µM respectively (Figure 1—figure supplement 2A), which is closest to K_d1 of WT. Hence, the affinity for Mn²⁺ 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 Mn²⁺ transport cycle by DraNramp

We also determined high-resolution DraNramp structures in metal-free occluded (WT) and Mn²⁺-bound inward-open (M230A•Mn²⁺) states and re-refined a Mn²⁺-bound outward-open conformation (G223W•Mn²⁺; 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 Mn²⁺ transport cycle to three major conformations, each in Mn²⁺-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.

Figure 2 with 2 supplements see all

Download asset Open asset

Structures in new conformations complete the Mn²⁺ transport cycle by DraNramp.

(A) Schematic of the conformational states that DraNramp traverses to import Mn²⁺. The mobile and stationary parts are pale yellow and light blue, respectively. (B) Corresponding structures of DraNramp, showing TMs 1 and 6 in green, TMs 5 and 10 in pale yellow, and TMs 2, 7, and 11 gray. Stationary TMs 3, 4, 8, and 9 are omitted to highlight the key motions in the mobile parts. Mn²⁺ ions are magenta. Black arrows indicate the key motions in TMs 1a and 6a detailed in panel (C), and TMs 5 and 10 detailed in Figure 2—figure supplement 1. Full structures and the electron density for TMs 1 and 6 are illustrated in Figure 2—figure supplement 2. (C) Pairwise superpositions of whole Mn²⁺-bound structures highlight the motions of TMs 1 and 6. Conformations are indicated at the bottom. The distance between residues 46 and 240 in TMs 1a and 6b, indicated for the green structures, increases from 5.9 Å to 13.8 Å from outward open to inward open. The large angular motions of TM1a and TM6b are also indicated. (D) Plots of B-factor by residue for the TM1 region (residues 40–70) and the TM6 region. The B-factors are highest for the inward-open state in which the interaction between TMs 1a and 6b (both in the inner leaflet) is broken.

TMs 1, 5, 6, and 10 move the most as the Mn²⁺-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 Mn²⁺ 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 Mn²⁺-bound occluded and inward-open states, allowing release of the Mn²⁺ into the inner vestibule. TM1a swings to a similar angle in the new inward-open M230A•Mn²⁺ 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 Mn²⁺-bound structures, suggesting that once Mn²⁺ 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 Mn²⁺ (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 Mn²⁺ to close the outer gate. The central unwound regions of TM1 and TM6 are closest in the occluded state, resulting in an almost dehydrated Mn²⁺-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 Mn²⁺ coordinations

Our Nramp structures provide snapshots of the complete Mn²⁺-coordination sphere geometries in each conformation (Figure 3). Two new structures of DraNramp point-mutants reveal the coordination of Mn²⁺ in the inward-open state: M230A•Mn²⁺ and D296A•Mn²⁺ (Cα RMSD of 0.42 Å). The Mn²⁺ is in the same location of the orthosteric site as in the occluded state, as confirmed by anomalous diffraction for D296A•Mn²⁺ (Supplementary file 1d, Figure 3—figure supplement 1). As in the occluded state, the Mn²⁺ binds D56, N59, and the A227 carbonyl, with a water replacing M230 in M230A•Mn²⁺ (Figure 3A, Supplementary file 1e). However, with TM1a displaced, the A53 carbonyl no longer coordinates Mn²⁺; instead, the Y54 carbonyl approaches Mn²⁺ 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 Mn²⁺-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 Mn²⁺ transport cycle.

Figure 3 with 3 supplements see all

Download asset Open asset

Coordination sphere changes across the Mn²⁺ transport cycle of DraNramp.

(A) Structures of the orthosteric metal-binding site in six conformations reveal the differences in coordination geometry and illustrate that the bound Mn²⁺ is more hydrated in the outward-open and inward-open states than the occluded state. In the occluded structure of metal-free WT DraNramp a density we have assigned as water replaces Mn²⁺. Y54 in TM1a progressively moves to open the inner vestibule in the transition from outward to inward open, shown by black curved arrows. (B) TM1 and TM6 from a superposition of the three Mn²⁺-bound structures in panel a illustrate the swing of the Y54 sidechain as sticks. The view is rotated 180° along the vertical axis from Figure 2C. (C) Initial Mn²⁺ uptake rates for DraNramp variants Y54A and Y54F at membrane potentials ranging from ΔΨ=0 to −120 mV (n=2–3; each data point is on the scatter plot and black bars are the mean values). The Mn²⁺ concentration was 750 μM, and the pH was 7 on both sides of the membrane. Y54A nearly abolishes transport whereas Y54F has near-wildtype initial transport rates. Corresponding time traces are plotted in Figure 1—figure supplement 4. (D) ITC measurements of G223W (left; one-site binding model with fixed n=1) and M230A (right; two-site sequential binding model) binding to Mn²⁺. One isotherm is shown of two measured, and the listed K_d values are the average ± SEM (see Appendix 1 for ITC analysis).

The new occluded and inward-open Mn²⁺-bound structures have a monodentate coordination of D56 with Mn²⁺. For consistency, we reinterpreted the outward-open G223W•Mn²⁺ map (PDB ID: 6BU5) with a monodentate coordination of D56 with Mn²⁺ 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•Mn²⁺ has six Mn²⁺-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 2F_o-F_c maps of the closeup snapshots of their Mn²⁺-bound orthosteric site (Figure 3—figure supplement 1C, Figure 5—figure supplement 2).

Comparing the three Mn²⁺-bound conformations (Figure 3A), the carbonyls of A53 in TM1a and A227 in TM6b alternately coordinate Mn²⁺ in the outward- and inward-open structures respectively, and both residues interact with Mn²⁺ in the occluded structure. The pseudosymmetrically related A53 and A227 may thus act as hinges altering the Mn²⁺-coordination sphere as TM1 and TM6 move in turn to open the gates during Mn²⁺ 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 Mn²⁺ 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 Mn²⁺ transport, we purified and reconstituted into proteoliposomes the Y54A and Y54F variants. While Y54F has near-wildtype Mn²⁺-transport activity, Y54A nearly eliminates Mn²⁺ transport (Figure 3C), indicating that an aromatic ring is essential for the gating motion required for transport.

We also measured the Mn²⁺-binding affinity of the constructs that yielded inward- or outward- open structures, M230A and G223W, respectively. Mn²⁺ is present at the external site in the inward-open M230A•Mn²⁺ (Figure 1—figure supplement 3A), and consistently, the ITC data fit a two-site model (K_d1=215±65 µM assigned to the orthosteric site, K_d2=4600±1300 µM for the external site; Figure 3D and Appendix 1). The ITC data for G223W with Mn²⁺ fits only in a one-site model (K_d=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).

Mn²⁺ 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 WT_soak (2.36 Å; Table 1). The crystal used for WT_soak was mock-soaked (with no metal in the soaking solution). WT and WT_soak are nearly identical, confirming that the soaking process does not influence the conformational state. We analyzed WT_soak, unless otherwise noted. In WT_soak, we observed density but no anomalous signal at the orthosteric site and modeled a water molecule at the position where Mn²⁺ sits in the occluded state (Figure 3A). WT_soak is nearly identical to WT•Mn²⁺, 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 Mn²⁺ 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 Mn²⁺ transport (Bozzi et al., 2020; Bozzi et al., 2019b). To pinpoint protein features that enable these motions, we used our Mn²⁺-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).

Figure 4 with 3 supplements see all

Download asset Open asset

Networks of polar residues lining the outer and inner vestibules rearrange through the conformational transitions needed for Mn²⁺ transport.

(A) The Cα positions of residues in the Q378 (orange) and T228 (pink) networks lining the outer gate, R244 (cyan) and Q89 (blue) networks within the inner gate, and the H232 (orange) network coordinating with the proton pathway, are mapped on the occluded WT•Mn²⁺ structure. Metal-binding and proton pathway residues (Bozzi et al., 2019b) are represented as brown spheres. (B) The Q378 network forms as TM10 moves when DraNramp transitions from outward-open to occluded to close the outer vestibule. Water-mediated interactions form between Q378, D56, and T130. The other D56 carbonyl interacts directly with Mn²⁺. TM1 is transparent. (C) The T228 network forms with N275, N82, and T228 coordinating a water as TM6a moves to close the outer vestibule. TM1 is transparent. (D) In the R244 network, interactions between R244, E176, and D263 break as TM5 moves in the transition from outward-open to occluded state to initiate the opening of the inner vestibule. (E) In the Q89 network, Y54, Q89, and H237 rearrange from occluded to inward-open state as TM1a swings up to allow for metal release. (F) H232, which abuts the orthosteric Mn²⁺-binding site, interacts with E134 and T130 through waters conserved in all conformations. In the M230A•Mn²⁺ structure, H232 flips and is replaced by a water, retaining the interaction with T130 but breaking the connection with E134, suggesting that M230 helps stably position H232. TM8 is transparent. In panels c-e, TMs 3, 4, 8, and 9 are omitted to better visualize the interactions. The illustrated structures are indicated on the figure.

Two networks seal the outer gate. In the occluded and inward-open conformations, Q378 interacts with two waters, one coordinating the orthosteric Mn²⁺ 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 Mn²⁺-coordination sphere (Figure 1B). As the inner gate opens to release Mn²⁺, 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 Mn²⁺ 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 Mn²⁺ 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 Mn²⁺ 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 Mn²⁺ 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•Mn²⁺ 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 Mn²⁺ binding, Cd²⁺ binding to DraNramp is exothermic

Nramps transport divalent transition metals quite promiscuously, including both physiological (Fe²⁺ and Mn²⁺) and non-physiological substrates (Cd²⁺, Zn²⁺, Co²⁺, Ni²⁺, Pb²⁺), but select against alkaline earth metals (Mg²⁺, Ca²⁺) (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 Cd²⁺ 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 Mn²⁺ and Cd²⁺.

In contrast to endothermic binding of Mn²⁺, ITC measurements show exothermic binding of Cd²⁺ to WT DraNramp (Figure 5—figure supplement 1B; see Appendix 1 for details of the ITC analyses). Like for Mn²⁺, the Cd²⁺ isotherm fits best in a two-site model (K_d1=55±15 µM, K_d2=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 (K_d=120±1 µM for D296A and K_d=70±10 µM for D369A), we conclude that the orthosteric site has higher affinity toward Cd²⁺ than the external site (Figure 5A–B). Both the orthosteric and external site shows higher affinity toward Cd²⁺ than Mn²⁺. The affinities are comparable at the orthosteric site, in low micromolar range for both metals (K_d of 190±30 µM for Mn²⁺ vs. 55±15 µM for Cd²⁺), whereas the external site has a much higher affinity toward Cd²⁺ (K_d of 1970±520 µM for Mn²⁺ vs. 220±20 µM for Cd²⁺). Consistent with its inability to transport Mg²⁺, a representative alkaline earth metal (Bozzi and Gaudet, 2021; Bozzi et al., 2019b), DraNramp does not bind Mg²⁺ (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 Cd²⁺ 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.

Figure 5 with 4 supplements see all

Download asset Open asset

DraNramp binds differently to Mn²⁺ and Cd²⁺.

(A) ITC measurements show that DraNramp binds Mn²⁺ in an endothermic mode, Cd²⁺ in an exothermic mode and does not bind Mg²⁺ (6 mM metal). ITC was performed with DraNramp in which the external site is mutated (D296A) to analyze binding at the orthosteric site. One isotherm is shown of 2 or 3 measured. These isotherms were fit using a one-site model with fixed n=1, and the listed K_d values are the average ± SEM (see Appendix 1 for ITC analysis). (B) Schematic showing the ITC-measured K_d values of various DraNramp constructs for Mn²⁺ (top) and Cd²⁺ (bottom). WT DraNramp binds Mn²⁺ and Cd²⁺ at the same two sites, with the external-site affinity ~10-fold higher and the orthosteric-site affinity ~threefold higher for Cd²⁺ than Mn²⁺. The M230A mutation eliminates binding of Cd²⁺ but not Mn²⁺ at the orthosteric site. The D369A mutation eliminates binding of either metal at the external site. A variant with mutations at both the orthosteric and external sites, M230A-D369A, does not bind Cd²⁺ but maintains orthosteric site binding for Mn²⁺. For all ITC data, the K_d values were computed while fixing the number of sites (n) to 1 or 2 and assigned to the external or orthosteric site based on knowledge of the crystal structures and mutational analysis (Appendix 1; Supplementary file 1g). ITC traces are shown in Figure 5—figure supplement 1A, Figure 5—figure supplement 2 and Figure 5—figure supplement 3. (C) Comparison of the inward-open state bound to Mn²⁺ (M230A•Mn²⁺; top) and Cd²⁺ (WT•Cd²⁺; bottom) shows differences in coordination geometry at the orthosteric site. D56 is oriented differently and does not directly coordinate Cd²⁺. Cd²⁺ coordinates N59, M230, carbonyls of A227 and Y54, and two waters, for a total of six ligands compared to seven for Mn²⁺. Mn²⁺ and Cd²⁺ are magenta and brown spheres, respectively. Of note, we use M230A•Mn²⁺ here because soaking of WT crystals with Cd²⁺ yielded an inward-open WT•Cd²⁺ structure (Figure 5—figure supplement 1C) whereas soaking the same kind of crystals with Mn²⁺ yielded an occluded WT•Mn²⁺ structure (Figure 1).

Figure 5—source data 1 Source files (Origin files) of ITC experiments of magnesium binding to each DraNramp construct.: https://cdn.elifesciences.org/articles/84006/elife-84006-fig5-data1-v2.zip
Download elife-84006-fig5-data1-v2.zip

A Cd²⁺-bound structure helps explain functional differences

To understand the differences in binding and transport of Mn²⁺ and Cd²⁺, we determined a 2.5 Å Cd²⁺-bound structure by soaking WT DraNramp crystals with 2 mM Cd²⁺ (Table 1). In agreement with the ITC data, WT•Cd²⁺ shows Cd²⁺ 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 Cd²⁺ than Mn²⁺ at the external site suggests that local geometry favors non-physiological Cd²⁺ over the physiological substrate Mn²⁺ (Figure 1—figure supplement 3A). The low affinity of this site for Mn²⁺, compared to Cd²⁺, 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 Mn²⁺.

Interestingly, WT•Cd²⁺ is inward open (Figure 5—figure supplement 1C), although both mocked-soaked and Mn²⁺-soaked crystals under otherwise equivalent conditions yielded occluded structures (WT_soak and WT•Mn²⁺, respectively; Figures 1 and 3). The larger ionic radius of Cd²⁺ (~0.95 Å) compared to Mn²⁺ (~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 Cd²⁺ bound. Compared with the inward-open M230A•Mn²⁺, D56 adopts a different rotamer in WT•Cd²⁺ and does not coordinate the Cd²⁺ 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. Cd²⁺ has six coordinating ligands, and the distortion from ideal octahedral geometry is more pronounced compared to Mn²⁺ (Supplementary file 1f). The coordination distances are larger for Cd²⁺ than for Mn²⁺ (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. Cd²⁺ 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.

Cd²⁺ binding is prone to perturbations and favors the inward-open state

ITC data analysis of different DraNramp variants highlights additional differences between Mn²⁺ and Cd²⁺ 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, Mn²⁺ and Cd²⁺ binding are entropy- and enthalpy-driven, respectively, and entropy contributions (-TΔS) are smaller for Cd²⁺ than Mn²⁺ (Appendix 1), suggesting that Cd²⁺ 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, Mn²⁺ binding fits with a two-site model and Cd²⁺ fits only with a one-site model (Figure 5B and Figure 5—figure supplement 2, Appendix 1). M230 is thus a crucial ligand for Cd²⁺ but not Mn²⁺, which is further reflected in the transport behavior where M230A affects Cd²⁺ transport drastically but has negligible effect on Mn²⁺ (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 Mn²⁺-binding isotherms, which we assigned as binding to the orthosteric site. However, all four double mutants showed complete loss of Cd²⁺ binding, which is expected for ones with M230A—which eliminates Cd²⁺ 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 Mn²⁺ substrate at the orthosteric site is more robust to perturbations than binding of Cd²⁺, 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 Mn²⁺ and Cd²⁺ (Figure 5—figure supplement 3). ITC with Mn²⁺ 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 Cd²⁺, 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 Cd²⁺ at the orthosteric site of A47W and G223W. This is consistent with our inability to obtain Cd²⁺-bound structures in conformations other than inward-open, despite trying both co-crystallization and soaking with A47W and G223W. These data suggest that Cd²⁺ 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 WT_soak structure) produced an inward-open Cd²⁺-bound structure.

Gating network differences in the Cd²⁺-bound structure are restricted to D56

To better understand the Cd²⁺ transport mechanism given that its binding at the orthosteric site seems less optimal and robust than Mn²⁺, we compared the gating networks described above for the inward-open Mn²⁺- and Cd²⁺-bound structures. The Q89, T228, and R244 networks are essentially identical in the inward-open state, regardless of whether Mn²⁺ or Cd²⁺ 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 Cd²⁺ 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 Cd²⁺ and Q378 to close the outer gate. Most of the H232 network is similar in the Mn²⁺- and Cd²⁺-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. Mn²⁺ and Cd²⁺ follow similar trends although with less voltage dependence for Cd²⁺ than Mn²⁺. These results support the idea that although Mn²⁺ and Cd²⁺ 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 Mn²⁺ import pathway in DraNramp (Video 1). Binding of the Mn²⁺ 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 Mn²⁺ and Cd²⁺ and corresponding analyses of substrate binding demonstrate that the orthosteric site shows similar affinity for physiological (Mn²⁺) and toxic (Cd²⁺) substrates, but Cd²⁺ binding is less robust to various perturbations.

Video 1

Download asset

posterframe for video — Conformational rearrangement and distinct coordination geometry adopted by DraNramp during Mn²⁺ import.

The global conformational changes are highlighted first, with cartoon representations of each of the six structures in the Mn²⁺ transport cycle starting from the outward-open Mn²⁺-bound conformation rotating through the occluded Mn²⁺-bound, inward-open Mn²⁺-bound, and then the metal-free states transitioning back to outward-open Mn²⁺-bound form. The second set of scenes then focusses on the orthosteric site, cycling through the structures in the same order to illustrate the distinct coordination geometries of the bound Mn²⁺, and the corresponding positions of its interacting ligands (protein residues and water) in the metal-free states.

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 Mn²⁺ transport. Starting in the outward-open state, Mn²⁺ 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 Mn²⁺ and restrict solvent access from the extracellular side in the occluded structure. A227 replaces a water of the outward-open Mn²⁺-coordination sphere, thus retaining a six-coordination geometry in both conformations. Mn²⁺ 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).

Figure 6 with 1 supplement see all

Download asset Open asset

A structure-guided model of the conformational cycle and thermodynamic landscape of metal transport by Nramps.

(A) After Mn²⁺ enters through the outer vestibule between TM6a and TM10 in the outward-open state, a bulk conformation change closes the outer gate. The occluded conformation arises though rearrangements of TM6a and TM10 facilitated by formation of the T228 and Q378 networks, respectively. The inner gate partially opens in the occluded state as the R244 network breaks and TM5 moves. To achieve the inward-open conformation, disruption of the Q89 network frees TM1a to swing up to fully open the inner vestibule for Mn²⁺ release into the cytosol. (B) Our data indicate that the most stable Mn²⁺-bound state is the occluded state, and the three main states are readily accessible to facilitate transport. In contrast, Cd²⁺ binding stabilizes the inward-open state.

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 Mn²⁺ in the inward-open structure is less favored and more distorted, which may facilitate Mn²⁺ 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 Mn²⁺ 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 Mn²⁺ but are not metal transporters, like the Mn²⁺ 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 Mn²⁺-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 Mn²⁺-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 Mn²⁺.

While the Mn²⁺ 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 Mn²⁺-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 Mn²⁺ 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•Mn²⁺ and D296A•Mn²⁺). 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•Cd²⁺ retains the M230 sulfur as a coordinating ligand but excludes the D56 carboxylate can be rationalized by the fact that Cd²⁺ is a softer metal than Mn²⁺. 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 Mn²⁺ but not Cd²⁺ (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, Cd²⁺ tends to form weaker complexes than Mn²⁺ 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 Cd²⁺-bound DraNramp conformations will be less stable than the Mn²⁺-bound states. Our ITC data also indicate that binding of the toxic Cd²⁺ to DraNramp is less robust to perturbations compared to its physiological substrate Mn²⁺. Furthermore, Cd²⁺ only binds well to the inward-open state, and its binding is exothermic rather than endothermic for Mn²⁺ (Figure 5A). The crystallographic data similarly suggest that inward-open is the most stable Cd²⁺-bound state (Figure 6B), based on the following observations: (i) soaking crystals of occluded WT yielded inward-occluded WT•Cd²⁺, (ii) co-crystallization efforts yielded no other Cd²⁺-bound structures, and (iii) soaking Cd²⁺ 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 Mn²⁺ (and likely the similar ion Fe²⁺). The distinct interactions of Nramps with Mn²⁺ and Cd²⁺ 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

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Gene (Deinococcus radiodurans)	DraNramp	Genomic DNA	Uniprot: Q9RTP8
Strain, strain background (Escherichia coli)	C41(DE3)	Lucigen	60442–1	Chemically competent cells
Strain, strain background (Escherichia coli)	DH5α	Invitrogen	18265–017	Chemically competent cells
Recombinant DNA reagent	pET21a	Novagen	69740–3	Vector backbone for cloning DraNramp
Chemical compound, drug	n-Dodecyl-β-D-Maltopyranoside	Anatrace	D310S	1% for solubilizing membrane, 0.03% used in wash buffer
Chemical compound, drug	n-Decyl-β-D-Maltopyranoside	Anatrace	D322S	0.1% used in exchange buffer
Chemical compound, drug	Lauryl maltose neopentyl glycol	Anatrace	NG310	0.01% in elution buffer, 0.003% used in SEC buffer
Chemical compound, drug	Monoolein	Anatrace	LCP18	1:1.5 (protein: monoolein) for LCP crystallization
Chemical compound, drug	Fura-2 Pentapotassium Salt, cell impermeant	Life Technologies	F-1200	5 mM mixed with in proteoliposome
Chemical compound, drug	1-Palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG)	Avanti Polar Lipids	850457 C	lipid mixture (3 POPE: 1 POPG) and protein in 400:1 ratio for proteoliposome preparation
Chemical compound, drug	1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE)	Avanti Polar Lipids	840757 C	lipid mixture (3 POPE: 1 POPG) and protein in 400:1 ratio for proteoliposome preparation
Chemical compound, drug	Valinomycin	Sigma-Aldrich	V0627	Creates membrane potential by transporting K⁺ in proteoliposome assay
Commercial assay, kit	Ni Sepharose High Performance Resin	Cytiva	95055–838	IMAC resin
Commercial assay, kit	Superdex 200 10/300 GL	Cytiva	89497–272	SEC column
Commercial assay, kit	PD-10 Desalting Column	Cytiva	95017–001	Buffer exchange
Software, algorithm	XDS	PMID:20124692	RRID:SCR_015652	Data Porcessing
Software, algorithm	PHASER	PMID:19461840	RRID:SCR_014219	Model Builiding
Software, algorithm	PHENIX	PMID:22505256	RRID:SCR_014224	Refinement
Software, algorithm	coot	PMID:20383002	RRID:SCR_014222	Refinement and model building
Software, algorithm	PyMOL	Schrödinger	RRID:SCR_000305	Figure making
Software, algorithm	HMMER	hmmer.org	RRID:SCR_005305	Collected sequences with jackhmmer and aligned with hmmalign
Software, algorithm	MUSCLE	PMID:15318951	RRID:SCR_011812	Aligning sequences
Software, algorithm	RAxML-NG	PMID:31070718	RRID:SCR_022066	Tree building
Software, algorithm	PropKa	PMID:26596171		Identifying protonation states
Software, algorithm	CHARMM-GUI	PMID:25130509		Building MD system
Software, algorithm	NAMD	PMID:20675161	RRID:SCR_014894	Running MD simulations
Software, algorithm	mdtraj	PMID:26488642		Analyzing MD data
Software, algorithm	VMD	PMID:8744570	RRID:SCR_001820	Analyzing MD data

Cloning and protein expression vectors

Request a detailed protocol

The 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 (Mn²⁺, Cd²⁺) 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 protocol

Protein 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 protocol

Cells 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 protocol

To 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 protocol

Crystallization 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 MnCl₂) co-crystals (A47W•Mn²⁺, M230A•Mn²⁺, D296A•Mn²⁺) 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•Cd²⁺, WT•Mn²⁺, WT_soak) 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 protocol

Diffraction 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 (Mn²⁺ or Cd²⁺) in the binding site. For D296A•Mn²⁺ and A47W•Mn²⁺, 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•Cd²⁺, 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 protocol

Models 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•Cd²⁺, ‘anomalous group refinement’ was used to improve the fit to density of the Cd²⁺ ions, with Cd²⁺ 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•Mn²⁺ and WT•Cd²⁺, 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 Mn²⁺ (D296A•Mn²⁺ and A47W•Mn²⁺) and Cd²⁺ (WT•Cd²⁺) 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 Mn²⁺ and Cd²⁺ in all structures. All software were provided by SBGrid (Morin et al., 2013).

Metal binding measurements using ITC

Request a detailed protocol

ITC experiments were performed using MicroCal iTC200 (GE Healthcare) to determine the affinity and thermodynamic parameters of binding of divalent metals (Mn²⁺, Cd²⁺ and Mg²⁺) 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 K_d 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

View detailed protocol

Protein 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 protocol

We 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 protocol

Molecular dynamics (MD) simulations were initialized from three high-resolution structures of DraNramp: the outward-open G223W•Mn²⁺ structure (6BU5) with Mn²⁺ removed and W223 mutated back to the native glycine residue in silico, the inward-open WT•Cd²⁺ structure with Cd²⁺ 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 Å³ 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 protocol

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 (WT_soak), 8E60 (WT•Mn²⁺), 8E6H (A47W•Mn²⁺), 8E6I (M230A•Mn²⁺), 8E6L (D296A•Mn²⁺), 8E6M (WT•Cd²⁺), and 8E6N (re-refined G223W•Mn²⁺). 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 (Mn²⁺ isotherms), Appendix 1—table 2—source data 1 (Cd²⁺ isotherms) and Figure 5—source data 1(Mg²⁺ 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, K_d, (or association constant, K_a) 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 K_d 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 K_d 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

ITC analysis of all Mn²⁺ binding isotherms, with chosen model and values shaded.

Protein construct	Two-site sequential binding model (fixed n=2)								One-site fit with fixed n=1
Protein construct	K_d1(µM)	ΔH₁(kcal/mol)	-TΔS₁(kcal/mol)	ΔG₁(kcal/mol)	K_d2(µM)	ΔH₂(kcal/mol)	-TΔS₂(kcal/mol)	ΔG₂(kcal/mol)	K_d(µM)	ΔH(kcal/mol)	-TΔS(kcal/mol)	ΔG(kcal/mol)
WT	220	3.2	–8.1	-5	961	5.4	–9.2	–3.8	650	8.8	–13.1	–4.3
	125	5.7	–11.1	–5.4	2700	17.1	–20.5	–4.6	380	12.6	–16.9	–4.3
	220	8	–12.8	-5	2250	13.7	–17.3	–3.4	500	15.4	–19.6	–4.2
	190±30	5.6±1.3	–10.6±1.4	–5.1±0.1	1970±520	12.1±3.3	–15.7±3.1	–4.0±0.3	510±80	12.3±1.9	–16.5±1.8	–4.2±0.03
A47W	130	6.3	–11.3	-5	3100	20.6	–23.8	–3.2	380	13.6	–18.1	–4.5
	120	6.2	–11.3	–5.1	1800	11.4	–14.9	–3.5	300	11.9	–16.6	–4.7
	125±5	6.2±0.5	–11.3±0.0	–5.0±0.1	2450±650	16±4.5	–19.3±4.4	–3.3±0.1	340±40	12.7±0.8	–17.3±0.7	–4.6±0.1
D56A	305	2.4	–7.1	–4.7	4800	20.4	–23.5	–3.1	2200	16.3	–19.6	–3.3
	140	1.8	–6.8	-5	2700	22.6	–25.9	–3.3	2500	22.6	–26.2	–3.6
	230±80	2.1±0.3	–6.9±0.1	–4.8±0.1	3800±1,100	21.5±1.1	–24.7±0.1	–3.2±0.1	2350±150	19.4±3	–22.9±3.3	–3.4±0.1
M230A	280	3	–7.7	–4.7	3300	14.9	–18.1	–3.2	1500	13.3	–16.9	–3.6
	150	2.9	-8	–5.1	5800	29.9	–33.1	–3.2	970	12.5	–16.3	–3.8
	215±65	2.9±0.05	–7.8±0.1	–4.9±0.2	4600±1,300	22.4±7.5	–25.6±7.5	–3.2±0	1240±270	12.9±0.4	–16.6±0.3	–3.7±0.1
G223W	Did not fit								460	5.8	–10.1	–4.3
									430	3	–7.4	–4.4
									440±15	4.4±1.4	–8.7±1.3	–4.3±0.1
D296A	Did not fit								430	8.9	–13.4	–4.5
									320	9.5	–14.3	–4.8
									370	10.1	–14.6	–4.5
									370±30	9.5±0.3	–14.1±0.3	–4.6±0.1
D369A	Did not fit								390	9.5	–14	–4.5
									480	8.9	–13.4	–4.5
									380	8.3	–12.8	–4.5
									420±30	8.9±0.3	–13.4±0.3	–4.5±0
A47W-D296A	Did not fit								300	6.1	–10.7	–4.6
									210	13.2	–17.5	–4.3
									255±45	9.6±3.5	–14.1±3.4	–4.4±0.1
A47W-D369A	Did not fit								210	7.5	–12.5	-5
									400	6.5	–11.1	–4.6
									300±95	7±0.5	–11.8±0.7	–4.8±0.1
D56A-D296A	Did not fit								170	5.3	–10.4	–5.1
									330	4.8	–9.5	–4.7
									250±80	5.0±0.2	–9.9±0.4	–4.9±0.2
D56A-D369A	Did not fit								280	4.9	–9.5	–4.6
									145	4	–8.9	–4.9
									210±70	4.4±0.4	–9.2±0.3	–4.7±0.1
M230A-D296A	Did not fit								320	7.1	–11.6	–4.5
									140	3.7	–8.9	–5.2
									230±90	5.4±1.7	–10.2±1.3	–4.8±0.3
M230A-D369A	Did not fit								750	11.5	–15.4	–3.9
									790	7.9	–11.9	-4
									770±20	9.7±1.8	–13.6±1.7	–3.9±0.05

Appendix 1—table 1—source data 1 Source files (Origin files) of ITC experiments of manganese binding to each DraNramp construct.: https://cdn.elifesciences.org/articles/84006/elife-84006-app1-table1-data1-v2.zip
Download elife-84006-app1-table1-data1-v2.zip

Appendix 1—table 2

ITC analysis of all Cd²⁺-binding isotherms, with chosen model and values shaded.

Protein construct	Two-site sequential binding model (fixed n=2)								One-site fit with fixed n=1
Protein construct	K_d1(µM)	ΔH₁(kcal/mol)	-TΔS₁ (kcal mol^–1)	K_d1(µM)	ΔH₁(kcal/mol)	ΔH₂ (kcal mol^–1)	K_d1(µM)	ΔH₁(kcal/mol)	K_d(µM)	K_d1(µM)	ΔH₁(kcal/mol)	ΔG (kcal mol^–1)
WT	85	–4.5	–0.9	–5.5	260	–2.6	-2	–4.6	105	–2.3	–2.9	–5.2
	50	–4.7	–0.9	–5.6	220	–4.1	–0.6	–4.7	105	-9	3.5	–5.5
	30	–2.9	–2.9	–5.8	180	–4.7	–0.3	-5	115	–7.6	2.1	–5.5
	55±15	–4.0±0.5	–1.6±0.6	–5.6±0.1	220±20	–3.8±0.6	–1.0±0.5	–4.8±0.1	110±3	6±2	–2.8±1.9	–5.4±0.1
A47W	220	–4.4	–0.2	–4.6	130	0.6^*	–5.6	-5	165	–4.8	–0.2	-5
	Did not fit								140	–3.3	–1.7	-5
									150±10	–4.0±0.7	–0.9±0.5	–5.0±0.0
D56A	100	–6.1	0.6	–5.5	3000	–11.2	7.7	–3.5	190	–9.2	4.1	–5.1
	55	–4.5	–1.2	–5.7	1500	–6.2	2.3	–3.9	130	–6.9	1.5	–5.4
	80±20	–5.3±0.8	–0.6±0.9	–5.6±0.1	2250±750	–8.7±2.5	5.0±0.5	–3.7±0.1	160±30	–8.0±1	–2.8±1.3	–5.2±0.1
M230A	Did not fit								142	–6.1	0.9	–5.2
									182	–3.1	–1.7	–4.8
									160±20	–4.6±1.5	–0.4±0.3	–5.0±0.2
G223W	No Binding
D296A	Did not fit								122	–3.5	–1.8	–6.8
									120	–3.7	–1.5	–5.2
									120±1	–3.6±0.1	–1.6±0.1	–6.0±0.8
D369A	100	–1.6	–3.5	–5.1	5350	2.8*	–5.9	–3.1	80	–1.4	–4.1	–5.5
	90	–2.5	–2.6	–5.7	4850	5^*	-8	-3	60	–2.1	–3.8	–5.9
	95±5	–5.3±0.8	–3.1±0.4	–5.4±0.1	5100±250	3.9±1.0*	–6.9±1.0	–3.0±0.1	70±10	–1.7±0.3	–3.9±0.1	–5.7±0.1
A47W-D296A	No Binding
A47W-D369A	No Binding
D56A-D296A	No Binding
D56A-D369A	No Binding
M230A-D296A	No Binding
M230A-D369A	No Binding

*

Positive values of ΔH corresponding to an exothermic mode of Cd²⁺ binding to protein indicates inappropriate fit

Appendix 1—table 2—source data 1 Source files (Origin files) of ITC experiments of cadmium binding to each DraNramp construct.: https://cdn.elifesciences.org/articles/84006/elife-84006-app1-table2-data1-v2.zip
Download elife-84006-app1-table2-data1-v2.zip

WT binding to Mn²⁺

The WT•Mn²⁺ crystal structure revealed two Mn²⁺ 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 K_d 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 K_d1 to the orthosteric site (190±30 µM) and the lower affinity K_d2 to the external site (1970±520 µM).

Appendix 1—figure 1

Download asset Open asset

Comparison of one-site and two-site fit to ITC binding isotherms for each case where both types of fits gave results.

We show fits to a two-site sequential binding model (fixed n=2; left) or a one-site model with a fixed n=1 (right) for each case where both types of fits gave a result (Appendix 1—tables 1 and 2). For constructs•ligand pairs WT•Mn²⁺ (A), A47W•Mn²⁺ (B), M230A•Mn²⁺ (C), and D56A•Cd²⁺ (F), the two-site fit is better than the one-site fit, which corroborates findings from the crystal structures. For pairs WT•Cd²⁺ (E) and D56A•Mn²⁺ (D), the fits are similar, and we used additional biochemical and structural information to select the most appropriate fit.

A47W binding to Mn²⁺

The A47W•Mn²⁺ crystal structure revealed two Mn²⁺ ions bound to the same sites as with WT•Mn²⁺, 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 K_d 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 K_d assignments for WT, with the orthosteric site having higher affinity for Mn²⁺ (K_d1=125±5 µM) than the external site (K_d2=2450±650 µM).

M230A binding to Mn²⁺

The M230A•Mn²⁺ crystal structure again revealed two bound Mn²⁺, 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 K_d 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 K_d assignments for WT, with the orthosteric site having higher affinity for Mn²⁺ (K_d1=215±65 µM) than the external site (K_d2=4600±1300 µM).

D56A binding to Mn²⁺

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 Mn²⁺, as it agrees best with the overall analysis, with the orthosteric site having higher affinity for Mn²⁺ (K_d1=230±80 µM) than the external site (K_d2=3800±1100 µM).

WT binding to Cd²⁺

The WT•Cd²⁺ crystal structure reveals two Cd²⁺ ions bound at the same two distinct sites where Mn²⁺ 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 Mn²⁺ binding, the two K_d 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 K_d (K_d1=55±15 µM) to the orthosteric site, and the lower affinity one to the external site (K_d2=220±20 µM).

A47W binding to Cd²⁺

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 Cd²⁺ and we assigned the resulting K_d value (150±10 µM) to the external site. We conclude that Cd²⁺ does not bind to the A47W orthosteric site with measurable affinity.

M230A binding to Cd²⁺

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 K_d value obtained using the one-site model (160±20 µM) to the external site, and we conclude that Cd²⁺ does not bind to the M230A orthosteric site with measurable affinity.

D56A binding to Cd²⁺

The ITC data fit better when using a two-site sequential binding model than a one-site model (Appendix 1—figure 1F). The WT•Cd²⁺ structure shows no direct interaction of D56 with Cd²⁺, suggesting that the D56A mutation may not completely impair Cd²⁺ binding at the orthosteric site. We thus chose to use the two-site sequential binding model and assign higher affinity (K_d1=80±20 µM) to the external site, and lower affinity (K_d2=2250±750 µM) to the orthosteric site.

G223W binding to Mn²⁺

Consistent with the G223W•Mn²⁺ crystal structure in which Mn²⁺ is only bound at the orthosteric site, the ITC data could only be fitted with the one-site model. We assigned the resulting K_d value (440±15 µM) to the orthosteric site.

G223W binding to Cd²⁺

The ITC data show no binding, which is consistent with our inability to obtain a crystal structure of G223W bound to Cd²⁺. This indicates that the outward-locked G223W does not bind Cd²⁺ with measurable affinity at the orthosteric site.

D296A binding to Mn²⁺

Consistent with the D296A•Mn²⁺ crystal structure in which Mn²⁺ 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 K_d value (370±30 µM) to the orthosteric site.

D296A binding to Cd²⁺

The ITC data could only be fitted with the one-site model. We assigned the resulting K_d value (120±10 µM) to the orthosteric site, because the D296A mutation removes an aspartate ligand at the external site.

D369A binding to Mn²⁺

The ITC data could only be fitted with the one-site model. We assigned the resulting K_d value (420±30 µM) to the orthosteric site, because the D369A mutation removes an aspartate ligand at the external site.

D369A binding to Cd²⁺

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 K_d 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•Mn²⁺), 8E6H (A47W•Mn²⁺), 8E6I (M230A•Mn²⁺), 8E6L (D296A•Mn²⁺), 8E6M (WT•Cd²⁺), and 8E6N (re-refined G223W•Mn²⁺). 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 (Mn²⁺ isotherms), Appendix 1-table 2-source data 1 (Cd²⁺ isotherms) and Figure 5-source data 1 (Mg²⁺ isotherms).

The following data sets were generated

1. Wilson EA
2. Berry SP
3. Shekhar M
4. Gaudet R
5. Singharoy A
(2022) Dryad Digital Repository
Molecular 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
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E5S. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an occluded state.

https://www.rcsb.org/structure/8E5S
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E5V. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WTsoak in an occluded state.

https://www.rcsb.org/structure/8E5V
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E60. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an occluded, manganese-bound state.

https://www.rcsb.org/structure/8E60
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E6H. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter A47W mutant in an occluded, manganese-bound state.

https://www.rcsb.org/structure/8E6H
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E6I. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter M230A mutant in an inward-open, manganese-bound state.

https://www.rcsb.org/structure/8E6I
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E6L. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter D296A mutant in an inward-open, manganese-bound state.

https://www.rcsb.org/structure/8E6L
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E6M. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter WT in an inward-open, cadmium-bound state.

https://www.rcsb.org/structure/8E6M
1. Ray S
2. Gaudet R
(2022) RCSB Protein Data Bank
ID 8E6N. X-ray structure of the Deinococcus radiodurans Nramp/MntH divalent transition metal transporter G223W mutant in an outward-open, manganese-bound state.

https://www.rcsb.org/structure/8E6N
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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
1. Ray S
2. Gaudet R
(2022) SBGrid Data Bank
X-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

The following previously published data sets were used

(2019) SBGrid Data Bank
X-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

1. Andrews NC
(2002) Metal transporters and disease
Current Opinion in Chemical Biology 6:181–186.

https://doi.org/10.1016/s1367-5931(02)00307-1
- PubMed
- Google Scholar
(2017) Transition metal binding selectivity in proteins and its correlation with the phylogenomic classification of the cation diffusion facilitator protein family
Scientific Reports 7:16381.

https://doi.org/10.1038/s41598-017-16777-5
- PubMed
- Google Scholar
1. Bateman A
2. Martin MJ
3. Orchard S
4. Magrane M
5. Agivetova R
6. Ahmad S
7. Alpi E
8. Bowler-Barnett EH
9. Britto R
10. Bursteinas B
11. Bye-A-Jee H
12. Coetzee R
13. Cukura A
14. Da Silva A
15. Denny P
16. Dogan T
17. Ebenezer T
18. Fan J
19. Castro LG
20. Garmiri P
21. Georghiou G
22. Gonzales L
23. Hatton-Ellis E
24. Hussein A
25. Ignatchenko A
26. Insana G
27. Ishtiaq R
28. Jokinen P
29. Joshi V
30. Jyothi D
31. Lock A
32. Lopez R
33. Luciani A
34. Luo J
35. Lussi Y
36. MacDougall A
37. Madeira F
38. Mahmoudy M
39. Menchi M
40. Mishra A
41. Moulang K
42. Nightingale A
43. Oliveira CS
44. Pundir S
45. Qi G
46. Raj S
47. Rice D
48. Lopez MR
49. Saidi R
50. Sampson J
51. Sawford T
52. Speretta E
53. Turner E
54. Tyagi N
55. Vasudev P
56. Volynkin V
57. Warner K
58. Watkins X
59. Zaru R
60. Zellner H
61. Bridge A
62. Poux S
63. Redaschi N
64. Aimo L
65. Argoud-Puy G
66. Auchincloss A
67. Axelsen K
68. Bansal P
69. Baratin D
70. Blatter MC
71. Bolleman J
72. Boutet E
73. Breuza L
74. Casals-Casas C
75. de Castro E
76. Echioukh KC
77. Coudert E
78. Cuche B
79. Doche M
80. Dornevil D
81. Estreicher A
82. Famiglietti ML
83. Feuermann M
84. Gasteiger E
85. Gehant S
86. Gerritsen V
87. Gos A
88. Gruaz-Gumowski N
89. Hinz U
90. Hulo C
91. Hyka-Nouspikel N
92. Jungo F
93. Keller G
94. Kerhornou A
95. Lara V
96. Le Mercier P
97. Lieberherr D
98. Lombardot T
99. Martin X
100. Masson P
101. Morgat A
102. Neto TB
103. Paesano S
104. Pedruzzi I
105. Pilbout S
106. Pourcel L
107. Pozzato M
108. Pruess M
109. Rivoire C
110. Sigrist C
111. Sonesson K
112. Stutz A
113. Sundaram S
114. Tognolli M
115. Verbregue L
116. Wu CH
117. Arighi CN
118. Arminski L
119. Chen C
120. Chen Y
121. Garavelli JS
122. Huang H
123. Laiho K
124. McGarvey P
125. Natale DA
126. Ross K
127. Vinayaka CR
128. Wang Q
129. Wang Y
130. Yeh LS
131. Zhang J
132. Ruch P
133. Teodoro D
(2021) UniProt: the universal protein knowledgebase in 2021
Nucleic Acids Research 49:D480–D489.

https://doi.org/10.1093/nar/gkaa1100
- PubMed
- Google Scholar
1. Begg SL
2. Eijkelkamp BA
3. Luo Z
4. Couñago RM
5. Morey JR
6. Maher MJ
7. Ong C-LY
8. McEwan AG
9. Kobe B
10. O’Mara ML
11. Paton JC
12. McDevitt CA
(2015) Dysregulation of transition metal ion homeostasis is the molecular basis for cadmium toxicity in Streptococcus pneumoniae
Nature Communications 6:6418.

https://doi.org/10.1038/ncomms7418
- PubMed
- Google Scholar
1. Bleackley MR
2. Macgillivray RTA
(2011) Transition metal homeostasis: from yeast to human disease
Biometals 24:785–809.

https://doi.org/10.1007/s10534-011-9451-4
- PubMed
- Google Scholar
1. Bozzi AT
2. Bane LB
3. Weihofen WA
4. McCabe AL
5. Singharoy A
6. Chipot CJ
7. Schulten K
8. Gaudet R
(2016a) Conserved methionine dictates substrate preference in nramp-family divalent metal transporters
PNAS 113:10310–10315.

https://doi.org/10.1073/pnas.1607734113
- PubMed
- Google Scholar
1. Bozzi AT
2. Bane LB
3. Weihofen WA
4. Singharoy A
5. Guillen ER
6. Ploegh HL
7. Schulten K
8. Gaudet R
(2016b) Crystal structure and conformational change mechanism of a bacterial nramp-family divalent metal transporter
Structure 24:2102–2114.

https://doi.org/10.1016/j.str.2016.09.017
- PubMed
- Google Scholar
1. Bozzi AT
2. Bane LB
3. Zimanyi CM
4. Gaudet R
(2019a) Unique structural features in an nramp metal transporter impart substrate-specific proton cotransport and a kinetic bias to favor import
The Journal of General Physiology 151:1413–1429.

https://doi.org/10.1085/jgp.201912428
- PubMed
- Google Scholar
1. Bozzi AT
2. Zimanyi CM
3. Nicoludis JM
4. Lee BK
5. Zhang CH
6. Gaudet R
(2019b) Structures in multiple conformations reveal distinct transition metal and proton pathways in an Nramp transporter
eLife 8:e41124.

https://doi.org/10.7554/eLife.41124
- PubMed
- Google Scholar
(2020) Transmembrane helix 6b links proton and metal release pathways and drives conformational change in an nramp-family transition metal transporter
The Journal of Biological Chemistry 295:1212–1224.

https://doi.org/10.1074/jbc.RA119.011336
- PubMed
- Google Scholar
1. Bozzi AT
2. Gaudet R
(2021) Molecular mechanism of nramp-family transition metal transport
Journal of Molecular Biology 433:166991.

https://doi.org/10.1016/j.jmb.2021.166991
- PubMed
- Google Scholar
1. Budinger D
2. Barral S
3. Soo AKS
4. Kurian MA
(2021) The role of manganese dysregulation in neurological disease: emerging evidence
The Lancet. Neurology 20:956–968.

https://doi.org/10.1016/S1474-4422(21)00238-6
- PubMed
- Google Scholar
Conference
1. Cammack R
2. Hughes M
(2008)
Considerations for the Specification of Enzyme Assays Involving Metal Ions

Proceedings of the 3rd Beilstein ESCEC Symposium “Experimental Standard Conditions Of Enzyme Characterization.”.
- Google Scholar
Book
1. Cellier M
2. Gros P
(2004)
The NRAMP Family

New York, NY: Springer.
- Google Scholar
1. Chen PR
2. He C
(2008) Selective recognition of metal ions by metalloregulatory proteins
Current Opinion in Chemical Biology 12:214–221.

https://doi.org/10.1016/j.cbpa.2007.12.010
- PubMed
- Google Scholar
1. Cock PJA
2. Antao T
3. Chang JT
4. Chapman BA
5. Cox CJ
6. Dalke A
7. Friedberg I
8. Hamelryck T
9. Kauff F
10. Wilczynski B
11. de Hoon MJL
(2009) Biopython: freely available python tools for computational molecular biology and bioinformatics
Bioinformatics 25:1422–1423.

https://doi.org/10.1093/bioinformatics/btp163
- PubMed
- Google Scholar
(2016) X-Ray structures and mechanism of the human serotonin transporter
Nature 532:334–339.

https://doi.org/10.1038/nature17629
- PubMed
- Google Scholar
1. Coleman JA
2. Yang D
3. Zhao Z
4. Wen PC
5. Yoshioka C
6. Tajkhorshid E
7. Gouaux E
(2019) Serotonin transporter-ibogaine complexes illuminate mechanisms of inhibition and transport
Nature 569:141–145.

https://doi.org/10.1038/s41586-019-1135-1
- PubMed
- Google Scholar
1. Couñago RM
2. Ween MP
3. Begg SL
4. Bajaj M
5. Zuegg J
6. O’Mara ML
7. Cooper MA
8. McEwan AG
9. Paton JC
10. Kobe B
11. McDevitt CA
(2014) Imperfect coordination chemistry facilitates metal ion release in the PSA permease
Nature Chemical Biology 10:35–41.

https://doi.org/10.1038/nchembio.1382
- PubMed
- Google Scholar
(2004) Weblogo: a sequence logo generator
Genome Research 14:1188–1190.

https://doi.org/10.1101/gr.849004
- PubMed
- Google Scholar
(1989) Identification of transferrin as the major plasma carrier protein for manganese introduced orally or intravenously or after in vitro addition in the rat
The Journal of Nutrition 119:1461–1464.

https://doi.org/10.1093/jn/119.10.1461
- PubMed
- Google Scholar
(2022) Principles of alternating access in leut-fold transporters: commonalities and divergences
Journal of Molecular Biology 434:167746.

https://doi.org/10.1016/j.jmb.2022.167746
- PubMed
- Google Scholar
1. Dudev M
2. Wang J
3. Dudev T
4. Lim C
(2006) Factors governing the metal coordination number in metal complexes from cambridge structural database analyses
The Journal of Physical Chemistry. B 110:1889–1895.

https://doi.org/10.1021/jp054975n
- PubMed
- Google Scholar
1. Dudev T
2. Lim C
(2014) Competition among metal ions for protein binding sites: determinants of metal ion selectivity in proteins
Chemical Reviews 114:538–556.

https://doi.org/10.1021/cr4004665
- PubMed
- Google Scholar
1. Edgar RC
(2004) MUSCLE: a multiple sequence alignment method with reduced time and space complexity
BMC Bioinformatics 5:113.

https://doi.org/10.1186/1471-2105-5-113
- PubMed
- Google Scholar
(2014) Crystal structure of a SLC11 (NRAMP) transporter reveals the basis for transition-metal ion transport
Nature Structural & Molecular Biology 21:990–996.

https://doi.org/10.1038/nsmb.2904
- PubMed
- Google Scholar
(2017) Structural and mechanistic basis of proton-coupled metal ion transport in the SLC11/NRAMP family
Nature Communications 8:14033.

https://doi.org/10.1038/ncomms14033
- PubMed
- Google Scholar
1. Ekiz C
2. Agaoglu L
3. Karakas Z
4. Gurel N
5. Yalcin I
(2005) The effect of iron deficiency anemia on the function of the immune system
The Hematology Journal 5:579–583.

https://doi.org/10.1038/sj.thj.6200574
- PubMed
- Google Scholar
1. El-Gebali S
2. Mistry J
3. Bateman A
4. Eddy SR
5. Luciani A
6. Potter SC
7. Qureshi M
8. Richardson LJ
9. Salazar GA
10. Smart A
11. Sonnhammer ELL
12. Hirsh L
13. Paladin L
14. Piovesan D
15. Tosatto SCE
16. Finn RD
(2019) The pfam protein families database in 2019
Nucleic Acids Research 47:D427–D432.

https://doi.org/10.1093/nar/gky995
- PubMed
- Google Scholar
1. Emsley P
2. Lohkamp B
3. Scott WG
4. Cowtan K
(2010) Features and development of coot
Acta Crystallographica. Section D, Biological Crystallography 66:486–501.

https://doi.org/10.1107/S0907444910007493
- PubMed
- Google Scholar
1. Evans PR
2. Murshudov GN
(2013) 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
- PubMed
- Google Scholar
1. Faham S
2. Watanabe A
3. Besserer GM
4. Cascio D
5. Specht A
6. Hirayama BA
7. Wright EM
8. Abramson J
(2008) The crystal structure of a sodium galactose transporter reveals mechanistic insights into na+/sugar symport
Science 321:810–814.

https://doi.org/10.1126/science.1160406
- PubMed
- Google Scholar
1. Ferrante A
2. Gorski J
(2012) Enthalpy-Entropy compensation and cooperativity as thermodynamic epiphenomena of structural flexibility in ligand-receptor interactions
Journal of Molecular Biology 417:454–467.

https://doi.org/10.1016/j.jmb.2012.01.057
- PubMed
- Google Scholar
(2003) Structure of the manganese-bound manganese transport regulator of Bacillus subtilis
Nature Structural Biology 10:652–657.

https://doi.org/10.1038/nsb951
- PubMed
- Google Scholar
1. Harding MM
(2000) The geometry of metal-ligand interactions relevant to proteins. II. angles at the metal atom, additional weak metal-donor interactions
Acta Crystallographica. Section D, Biological Crystallography 56:857–867.

https://doi.org/10.1107/s0907444900005849
- PubMed
- Google Scholar
1. Harding MM
(2001) Geometry of metal-ligand interactions in proteins
Acta Crystallographica. Section D, Biological Crystallography 57:401–411.

https://doi.org/10.1107/s0907444900019168
- PubMed
- Google Scholar
(2015) Long-time-step molecular dynamics through hydrogen mass repartitioning
Journal of Chemical Theory and Computation 11:1864–1874.

https://doi.org/10.1021/ct5010406
- PubMed
- Google Scholar
1. Huang J
2. Rauscher S
3. Nawrocki G
4. Ran T
5. Feig M
6. de Groot BL
7. Grubmüller H
8. MacKerell AD
(2017) CHARMM36m: an improved force field for folded and intrinsically disordered proteins
Nature Methods 14:71–73.

https://doi.org/10.1038/nmeth.4067
- PubMed
- Google Scholar
1. I Bannon D
2. Portnoy ME
3. Olivi L
4. Lees PSJ
5. Culotta VC
6. Bressler JP
(2002) Uptake of lead and iron by divalent metal transporter 1 in yeast and mammalian cells
Biochemical and Biophysical Research Communications 295:978–984.

https://doi.org/10.1016/s0006-291x(02)00756-8
- PubMed
- Google Scholar
(2012) Substrate profile and metal-ion selectivity of human divalent metal-ion transporter-1
The Journal of Biological Chemistry 287:30485–30496.

https://doi.org/10.1074/jbc.M112.364208
- PubMed
- Google Scholar
1. Irving H
2. Williams RJP
(1953) 637. the stability of transition-metal complexes
Journal of the Chemical Society 3192.

https://doi.org/10.1039/jr9530003192
- Google Scholar
1. Jo S
2. Kim T
3. Iyer VG
4. Im W
(2008) CHARMM-GUI: a web-based graphical user interface for CHARMM
Journal of Computational Chemistry 29:1859–1865.

https://doi.org/10.1002/jcc.20945
- PubMed
- Google Scholar
(2015) Macrophage iron homeostasis and polarization in the context of cancer
Immunobiology 220:295–304.

https://doi.org/10.1016/j.imbio.2014.09.011
- PubMed
- Google Scholar
1. Jurrus E
2. Engel D
3. Star K
4. Monson K
5. Brandi J
6. Felberg LE
7. Brookes DH
8. Wilson L
9. Chen J
10. Liles K
11. Chun M
12. Li P
13. Gohara DW
14. Dolinsky T
15. Konecny R
16. Koes DR
17. Nielsen JE
18. Head-Gordon T
19. Geng W
20. Krasny R
21. Wei G-W
22. Holst MJ
23. McCammon JA
24. Baker NA
(2018) Improvements to the APBS biomolecular solvation software suite
Protein Science 27:112–128.

https://doi.org/10.1002/pro.3280
- PubMed
- Google Scholar
1. Kabsch W
(2010) XDS
Acta Crystallographica. Section D, Biological Crystallography 66:125–132.

https://doi.org/10.1107/S0907444909047337
- PubMed
- Google Scholar
1. Karplus PA
2. Diederichs K
(2012) Linking crystallographic model and data quality
Science 336:1030–1033.

https://doi.org/10.1126/science.1218231
- PubMed
- Google Scholar
1. Karplus PA
2. Diederichs K
(2015) Assessing and maximizing data quality in macromolecular crystallography
Current Opinion in Structural Biology 34:60–68.

https://doi.org/10.1016/j.sbi.2015.07.003
- PubMed
- Google Scholar
1. Kawabata H
(2019) Transferrin and transferrin receptors update
Free Radical Biology & Medicine 133:46–54.

https://doi.org/10.1016/j.freeradbiomed.2018.06.037
- PubMed
- Google Scholar
1. Kozlov AM
2. Darriba D
3. Flouri T
4. Morel B
5. Stamatakis A
6. Wren J
(2019) RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference
Bioinformatics 35:4453–4455.

https://doi.org/10.1093/bioinformatics/btz305
- PubMed
- Google Scholar
(2009) Unlocking the molecular secrets of sodium-coupled transporters
Nature 459:347–355.

https://doi.org/10.1038/nature08143
- PubMed
- Google Scholar
1. Krishnamurthy H
2. Gouaux E
(2012) X-ray structures of leut in substrate-free outward-open and apo inward-open states
Nature 481:469–474.

https://doi.org/10.1038/nature10737
- PubMed
- Google Scholar
(2005) Characterization of the metal ion binding site in the anti-terminator protein, hutp, of Bacillus subtilis
Nucleic Acids Research 33:5494–5502.

https://doi.org/10.1093/nar/gki868
- PubMed
- Google Scholar
1. Le SQ
2. Gascuel O
(2008) An improved general amino acid replacement matrix
Molecular Biology and Evolution 25:1307–1320.

https://doi.org/10.1093/molbev/msn067
- PubMed
- Google Scholar
1. Leavitt S
2. Freire E
(2001) Direct measurement of protein binding energetics by isothermal titration calorimetry
Current Opinion in Structural Biology 11:560–566.

https://doi.org/10.1016/s0959-440x(00)00248-7
- PubMed
- Google Scholar
1. Lee J
2. Cheng X
3. Swails JM
4. Yeom MS
5. Eastman PK
6. Lemkul JA
7. Wei S
8. Buckner J
9. Jeong JC
10. Qi Y
11. Jo S
12. Pande VS
13. Case DA
14. Brooks CL
15. MacKerell AD
16. Klauda JB
17. Im W
(2016) CHARMM-gui input generator for namd, gromacs, amber, openmm, and CHARMM/openmm simulations using the charmm36 additive force field
Journal of Chemical Theory and Computation 12:405–413.

https://doi.org/10.1021/acs.jctc.5b00935
- PubMed
- Google Scholar
1. Liebschner D
2. Afonine PV
3. Baker ML
4. Bunkóczi G
5. Chen VB
6. Croll TI
7. Hintze B
8. Hung LW
9. Jain S
10. McCoy AJ
11. Moriarty NW
12. Oeffner RD
13. Poon BK
14. Prisant MG
15. Read RJ
16. Richardson JS
17. Richardson DC
18. Sammito MD
19. Sobolev OV
20. Stockwell DH
21. Terwilliger TC
22. Urzhumtsev AG
23. Videau LL
24. Williams CJ
25. Adams PD
(2019) Macromolecular structure determination using x-rays, neutrons and electrons: recent developments in phenix
Acta Crystallographica. Section D, Structural Biology 75:861–877.

https://doi.org/10.1107/S2059798319011471
- PubMed
- Google Scholar
(2009)
Metals and membrane metal transporters in biological systems: the role (S) of Nramp in host-parasite interactions

J Argent Chem Soc 97:210–225.
- Google Scholar
(2021) Iron and manganese transport in mammalian systems
Biochimica et Biophysica Acta. Molecular Cell Research 1868:118890.

https://doi.org/10.1016/j.bbamcr.2020.118890
- PubMed
- Google Scholar
1. Loo DD
2. Hirayama BA
3. Gallardo EM
4. Lam JT
5. Turk E
6. Wright EM
(1998) Conformational changes couple Na⁺ and glucose transport
PNAS 95:7789–7794.

https://doi.org/10.1073/pnas.95.13.7789
- PubMed
- Google Scholar
(2009) Metal transporters and metal sensors: how coordination chemistry controls bacterial metal homeostasis
Chemical Reviews 109:4644.

https://doi.org/10.1021/cr900077w
- Google Scholar
(2019) Mechanistic basis of the inhibition of SLC11/NRAMP-mediated metal ion transport by bis-isothiourea substituted compounds
eLife 8:e51913.

https://doi.org/10.7554/eLife.51913
- PubMed
- Google Scholar
1. McGibbon RT
2. Beauchamp KA
3. Harrigan MP
4. Klein C
5. Swails JM
6. Hernández CX
7. Schwantes CR
8. Wang L-P
9. Lane TJ
10. Pande VS
(2015) MDTraj: A modern open library for the analysis of molecular dynamics trajectories
Biophysical Journal 109:1528–1532.

https://doi.org/10.1016/j.bpj.2015.08.015
- PubMed
- Google Scholar
1. Morin A
2. Eisenbraun B
3. Key J
4. Sanschagrin PC
5. Timony MA
6. Ottaviano M
7. Sliz P
(2013) Collaboration gets the most out of software
eLife 2:e01456.

https://doi.org/10.7554/eLife.01456
- PubMed
- Google Scholar
1. Nevo Y
2. Nelson N
(2006) The NRAMP family of metal-ion transporters
Biochimica et Biophysica Acta 1763:609–620.

https://doi.org/10.1016/j.bbamcr.2006.05.007
- PubMed
- Google Scholar
1. Nies DH
2. Grass G
(2009) Transition metal homeostasis
EcoSal Plus 3:.

https://doi.org/10.1128/ecosalplus.5.4.4.3
- PubMed
- Google Scholar
(2008) The thermodynamics of protein-ligand interaction and solvation: insights for ligand design
Journal of Molecular Biology 384:1002–1017.

https://doi.org/10.1016/j.jmb.2008.09.073
- PubMed
- Google Scholar
(2011) PROPKA3: consistent treatment of internal and surface residues in empirical pka predictions
Journal of Chemical Theory and Computation 7:525–537.

https://doi.org/10.1021/ct100578z
- PubMed
- Google Scholar
1. Perez C
2. Koshy C
3. Yildiz O
4. Ziegler C
(2012) Alternating-Access mechanism in conformationally asymmetric trimers of the betaine transporter BetP
Nature 490:126–130.

https://doi.org/10.1038/nature11403
- PubMed
- Google Scholar
Conference
(2008) Adapting a message-driven parallel application to GPU-accelerated clusters
SC - International Conference for High Performance Computing, Networking, Storage and Analysis.

https://doi.org/10.1109/SC.2008.5214716
- Google Scholar
(2009) Basis of substrate binding and conservation of selectivity in the CLC family of channels and transporters
Nature Structural & Molecular Biology 16:1294–1301.

https://doi.org/10.1038/nsmb.1704
- PubMed
- Google Scholar
1. Pittman JK
(2005) Managing the manganese: molecular mechanisms of manganese transport and homeostasis
The New Phytologist 167:733–742.

https://doi.org/10.1111/j.1469-8137.2005.01453.x
- PubMed
- Google Scholar
1. Potter SC
2. Luciani A
3. Eddy SR
4. Park Y
5. Lopez R
6. Finn RD
(2018) HMMER web server: 2018 update
Nucleic Acids Research 46:W200–W204.

https://doi.org/10.1093/nar/gky448
- PubMed
- Google Scholar
(2017) A novel proton transfer mechanism in the SLC11 family of divalent metal ion transporters
Scientific Reports 7:6194.

https://doi.org/10.1038/s41598-017-06446-y
- PubMed
- Google Scholar
(2022) Structural and functional properties of a magnesium transporter of the SLC11/NRAMP family
eLife 11:e74589.

https://doi.org/10.7554/eLife.74589
- PubMed
- Google Scholar
1. Rulísek L
2. Vondrásek J
(1998) Coordination geometries of selected transition metal ions (Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Hg2+) in metalloproteins
Journal of Inorganic Biochemistry 71:115–127.

https://doi.org/10.1016/s0162-0134(98)10042-9
- PubMed
- Google Scholar
(2001) Properties of the mammalian and yeast metal-ion transporters DCT1 and Smf1p expressed in Xenopus laevis oocytes
The Journal of Experimental Biology 204:1053–1061.

https://doi.org/10.1242/jeb.204.6.1053
- PubMed
- Google Scholar
1. Shimamura T
2. Weyand S
3. Beckstein O
4. Rutherford NG
5. Hadden JM
6. Sharples D
7. Sansom MSP
8. Iwata S
9. Henderson PJF
10. Cameron AD
(2010) Molecular basis of alternating access membrane transport by the sodium-hydantoin transporter mhp1
Science 328:470–473.

https://doi.org/10.1126/science.1186303
- PubMed
- Google Scholar
1. Simmons KJ
2. Jackson SM
3. Brueckner F
4. Patching SG
5. Beckstein O
6. Ivanova E
7. Geng T
8. Weyand S
9. Drew D
10. Lanigan J
11. Sharples DJ
12. Sansom MSP
13. Iwata S
14. Fishwick CWG
15. Johnson AP
16. Cameron AD
17. Henderson PJF
(2014) Molecular mechanism of ligand recognition by membrane transport protein, mhp1
The EMBO Journal 33:1831–1844.

https://doi.org/10.15252/embj.201387557
- PubMed
- Google Scholar
(2020)
Stability and Applications of Coordination

Stability and applications of coordination compounds, Stability and Applications of Coordination, InTechOpen, 10.5772/intechopen.83186.
- Google Scholar
(1998) Infection genomics: nramp1 as a major determinant of natural resistance to intracellular infections
Annual Review of Medicine 49:275–287.

https://doi.org/10.1146/annurev.med.49.1.275
- PubMed
- Google Scholar
(2011) Improved treatment of ligands and coupling effects in empirical calculation and rationalization of pKa values
Journal of Chemical Theory and Computation 7:2284–2295.

https://doi.org/10.1021/ct200133y
- PubMed
- Google Scholar
1. Tareen A
2. Kinney JB
(2020) Logomaker: beautiful sequence logos in python
Bioinformatics 36:2272–2274.

https://doi.org/10.1093/bioinformatics/btz921
- PubMed
- Google Scholar
1. Tellinghuisen J
(2008) Isothermal titration calorimetry at very low c
Analytical Biochemistry 373:395–397.

https://doi.org/10.1016/j.ab.2007.08.039
- PubMed
- Google Scholar
1. Turnbull WB
2. Daranas AH
(2003) 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
- PubMed
- Google Scholar
(2016) Different divalent cations alter the kinetics and fidelity of DNA polymerases
The Journal of Biological Chemistry 291:20869–20875.

https://doi.org/10.1074/jbc.R116.742494
- PubMed
- Google Scholar
1. Vidal SM
2. Malo D
3. Vogan K
4. Skamene E
5. Gros P
(1993) Natural resistance to infection with intracellular parasites: isolation of a candidate for bcg
Cell 73:469–485.

https://doi.org/10.1016/0092-8674(93)90135-d
- PubMed
- Google Scholar
1. Watanabe A
2. Choe S
3. Chaptal V
4. Rosenberg JM
5. Wright EM
6. Grabe M
7. Abramson J
(2010) The mechanism of sodium and substrate release from the binding pocket of vsglt
Nature 468:988–991.

https://doi.org/10.1038/nature09580
- PubMed
- Google Scholar
1. Winn MD
2. Ballard CC
3. Cowtan KD
4. Dodson EJ
5. Emsley P
6. Evans PR
7. Keegan RM
8. Krissinel EB
9. Leslie AGW
10. McCoy A
11. McNicholas SJ
12. Murshudov GN
13. Pannu NS
14. Potterton EA
15. Powell HR
16. Read RJ
17. Vagin A
18. Wilson KS
(2011) Overview of the CCP4 suite and current developments
Acta Crystallographica. Section D, Biological Crystallography 67:235–242.

https://doi.org/10.1107/S0907444910045749
- PubMed
- Google Scholar
(1989) Rapid measurement of binding constants and heats of binding using a new titration calorimeter
Analytical Biochemistry 179:131–137.

https://doi.org/10.1016/0003-2697(89)90213-3
- PubMed
- Google Scholar
1. Wu EL
2. Cheng X
3. Jo S
4. Rui H
5. Song KC
6. Dávila-Contreras EM
7. Qi Y
8. Lee J
9. Monje-Galvan V
10. Venable RM
11. Klauda JB
12. Im W
(2014) CHARMM-GUI membrane builder toward realistic biological membrane simulations
Journal of Computational Chemistry 35:1997–2004.

https://doi.org/10.1002/jcc.23702
- Google Scholar
1. Yokoyama H
2. Tsuruta O
3. Akao N
4. Fujii S
(2012) Crystal structure of Helicobacter pylori neutrophil-activating protein with a di-nuclear ferroxidase center in a zinc or cadmium-bound form
Biochemical and Biophysical Research Communications 422:745–750.

https://doi.org/10.1016/j.bbrc.2012.05.073
- PubMed
- Google Scholar

Article and author information

Author details

Shamayeeta Ray

Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States

Contribution
Conceptualization, Data curation, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-7906-0572
Samuel P Berry

Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States

Contribution
Data curation, Formal analysis, Investigation, Visualization, Writing - original draft, Writing – review and editing

Competing interests
No competing interests declared
Eric A Wilson

School of Molecular Sciences, Arizona State University, Tempe, United States

Contribution
Data curation, Formal analysis, Investigation, Methodology, Writing – review and editing

Competing interests
No competing interests declared
Casey H Zhang

Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States

Present address
Grossman School of Medicine, New York University, New York, United States

Contribution
Data curation, Formal analysis, Validation, Methodology, Writing – review and editing

Competing interests
No competing interests declared
Mrinal Shekhar

Broad Institute, Cambridge, United States

Contribution
Data curation, Investigation, Methodology, Writing – review and editing

Competing interests
No competing interests declared
Abhishek Singharoy

School of Molecular Sciences, Arizona State University, Tempe, United States

Contribution
Resources, Supervision, Funding acquisition, Methodology, Project administration, Writing – review and editing

Competing interests
No competing interests declared
Rachelle Gaudet

Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States

Contribution
Conceptualization, Resources, Supervision, Funding acquisition, Validation, Investigation, Visualization, Project administration, Writing – review and editing

For correspondence
gaudet@mcb.harvard.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-9177-054X

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

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,701

views
267

downloads
14

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Shamayeeta Ray
Samuel P Berry
Eric A Wilson
Casey H Zhang
Mrinal Shekhar
Abhishek Singharoy
Rachelle Gaudet

(2023)

High-resolution structures with bound Mn²⁺ and Cd²⁺ map the metal import pathway in an Nramp transporter

eLife 12:e84006.

https://doi.org/10.7554/eLife.84006

Categories and tags

Research organism

E. coli

Share this article

Cite this article

The occluded structure of DraNramp reveals a largely dehydrated Mn2+-coordination sphere.

Figure 1—source data 1

Figure 1—source data 2

Figure 1—source data 3

Data collection and refinement statistics for four new DraNramp structures.

Structures in new conformations complete the Mn2+ transport cycle by DraNramp.

Coordination sphere changes across the Mn2+ transport cycle of DraNramp.

Networks of polar residues lining the outer and inner vestibules rearrange through the conformational transitions needed for Mn2+ transport.

DraNramp binds differently to Mn2+ and Cd2+.

Figure 5—source data 1

Conformational rearrangement and distinct coordination geometry adopted by DraNramp during Mn2+ import.

A structure-guided model of the conformational cycle and thermodynamic landscape of metal transport by Nramps.

ITC analysis of all Mn2+ binding isotherms, with chosen model and values shaded.

Appendix 1—table 1—source data 1

ITC analysis of all Cd2+-binding isotherms, with chosen model and values shaded.

Appendix 1—table 2—source data 1

Comparison of one-site and two-site fit to ITC binding isotherms for each case where both types of fits gave results.

Author details

Shamayeeta Ray

Contribution

Competing interests

Samuel P Berry

Contribution

Competing interests

Eric A Wilson

Contribution

Competing interests

Casey H Zhang

Present address

Contribution

Competing interests

Mrinal Shekhar

Contribution

Competing interests

Abhishek Singharoy

Contribution

Competing interests

Rachelle Gaudet

Contribution

For correspondence

Competing interests

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)

Categories and tags

Research organism

Further reading

The occluded structure of DraNramp reveals a largely dehydrated Mn²⁺-coordination sphere.

Structures in new conformations complete the Mn²⁺ transport cycle by DraNramp.

Coordination sphere changes across the Mn²⁺ transport cycle of DraNramp.

Networks of polar residues lining the outer and inner vestibules rearrange through the conformational transitions needed for Mn²⁺ transport.

DraNramp binds differently to Mn²⁺ and Cd²⁺.

Conformational rearrangement and distinct coordination geometry adopted by DraNramp during Mn²⁺ import.

ITC analysis of all Mn²⁺ binding isotherms, with chosen model and values shaded.

ITC analysis of all Cd²⁺-binding isotherms, with chosen model and values shaded.