Structural basis for the reaction cycle of DASS dicarboxylate transporters
Abstract
Citrate, α-ketoglutarate and succinate are TCA cycle intermediates that also play essential roles in metabolic signaling and cellular regulation. These di- and tricarboxylates are imported into the cell by the divalent anion sodium symporter (DASS) family of plasma membrane transporters, which contains both cotransporters and exchangers. While DASS proteins transport substrates via an elevator mechanism, to date structures are only available for a single DASS cotransporter protein in a substrate-bound, inward-facing state. We report multiple cryo-EM and X-ray structures in four different states, including three hitherto unseen states, along with molecular dynamics simulations, of both a cotransporter and an exchanger. Comparison of these outward- and inward-facing structures reveal how the transport domain translates and rotates within the framework of the scaffold domain through the transport cycle. Additionally, we propose that DASS transporters ensure substrate coupling by a charge-compensation mechanism, and by structural changes upon substrate release.
Introduction
Citrate and dicarboxylates such as α-ketoglutarate (αKG), succinate and malate are intermediates of the TCA cycle. Furthermore, these molecules play essential roles in metabolic signaling and cellular regulation (Tannahill et al., 2013; Mills et al., 2018; Huergo and Dixon, 2015). In particular, citrate acts as a precursor of fatty acid synthesis and allosterically regulates both fatty acid synthesis and glycolysis. Citrate is also used to synthesize acetyl-CoA for histone acetylation and is therefore essential for the regulation of DNA transcription and replication (Wellen et al., 2009). Similarly, cytoplasmic αKG and succinate are important in controlling cell fate. Naive embryonic stem cells that exhibit an elevated αKG-to-succinate ratio maintain pluripotency (Carey et al., 2015). In contrast, pancreatic ductal adenocarcinoma cells with p53-deficiency have a lowed αKG-to-succinate ratio, and increasing the cellular concentration of αKG leads to a phenotype similar to that of tumor suppression by p53 restoration (Morris et al., 2019).
Mammalian cells import di- and tricarboxylates from the bloodstream via the Na+-dependent citrate transporter (NaCT) and the Na+-dependent dicarboxylate transporters 1 and 3 (NaDC1 and NaDC3). These plasma membrane proteins belong to the solute carrier 13 (SLC13) gene family (Markovich and Murer, 2004; Bergeron et al., 2013; Pajor, 2014). Loss-of-function mutations in the human NaCT transporter cause a type of encephalopathy (SLC13A5 Deficiency) (Thevenon et al., 2014; Hardies et al., 2015; Klotz et al., 2016). In contrast, knocking out NaCT in mice leads to protection from obesity and insulin resistance, while mutations in the homologous fly gene extends their lifespan (Birkenfeld et al., 2011; Rogina et al., 2000). Variants in the dicarboxylate transporter NaDC3 cause acute reversible leukoencephalopathy, with accumulation of αKG in the cerebrospinal fluid and urine (Dewulf et al., 2019). These central roles of SLC13 proteins in cell metabolism and signaling make them particularly attractive targets for treating obesity, diabetes, cancer and epilepsy (Huard et al., 2015; Pajor et al., 2016).
The mammalian SLC13 proteins are members of the larger divalent-anion sodium symporter (DASS) family (Markovich and Murer, 2004; Bergeron et al., 2013; Pajor, 2014; Prakash et al., 2003). DASS proteins typically have a molecular weight of 45–65 kDa and form obligate homodimers. The majority of DASS proteins are Na+-coupled cotransporters with a transport stoichiometry of one substrate to 2–4 Na+ ions (Figure 1—figure supplement 1a). However, other DASS members are exchangers (Pos et al., 1998), typically exchanging succinate for other dicarboxylates (Figure 1—figure supplement 1b). While sequence analysis strongly suggests a shared fold (Lolkema and Slotboom, 1998), these two groups separate into distinct clades on a phylogenetic tree and are thereby named DASS-C and DASS-E for cotransporter/symporter and exchanger/antiporter, respectively (Figure 1—figure supplement 1c).
All available DASS structures are of the Na+-driven dicarboxylate transporter VcINDY from Vibrio cholerae, in an inward-facing (Ci-Na+-S) state (Figure 1—figure supplement 1d–f; Mancusso et al., 2012; Nie et al., 2017). Each protomer consists of a scaffold and a transport domain. In the transport domain, the two carboxylate moieties of the bound substrate are coordinated by two conserved Ser-Asn-Thr (SNT) motifs, which also form part of the Na+-binding sites Na1 and Na2. Structural information, along with cross-linking and computer modeling (Mancusso et al., 2012; Mulligan et al., 2016), indicates that DASS proteins operate in an elevator-type transport mechanism (Figure 1a; Reyes et al., 2009; Drew and Boudker, 2016; Garaeva and Slotboom, 2020). However, without the structure of an outward-facing (Co) conformation structure, a detailed description of the changes between the Co and Ci conformations for DASS proteins remains lacking. How DASS exchangers translocate the substrates across the membrane is completely unknown.
Furthermore, it is not known how the transporter couples substrate binding and release to the Co to Ci interconversion (Stein, 1986). Studies on multiple DASS-C cotransporters using transport kinetics, electrophysiology and chemical cross-linking have shown that these proteins sequentially bind sodium and then substrate, with this order being reversed during the release process (Figure 1—figure supplement 1a; Wright et al., 1983; Yao and Pajor, 2000; Hall and Pajor, 2005; Pajor et al., 2013; Mulligan et al., 2014). The arrangement of the observed substrate and Na+-sites in the known VcINDY Ci-Na+-S state is also consistent with such a mechanism. However, it is unclear how Na+ slippage is avoided, namely, how the transporter undergoes interconversion between the apo Ci and Co states only, and between the fully loaded Ci-Na+-S and Co-Na+-S states, but not between the Na+-only loaded Ci-Na+ and Co-Na+ states. This information is essential to understanding the mechanism of SLC13 family disease mutations and the design of drugs to manipulate di- and tricarboxylate import. For the DASS-E exchangers, which we hypothesize follow a typical antiport kinetic mechanism (Figure 1—figure supplement 1b), it is even less clear how they translocate the substrates across the membrane.
We aimed to characterize the structural basis of the entire transport cycle of the DASS family using a combination of single particle cryo-EM, X-ray crystallography, molecular dynamics (MD) simulations and transport activity assays.
Results
Structure determination of multiple states
The only structurally characterized DASS transporter, VcINDY, has exclusively been observed in a Ci-Na+-S state (Figure 1—figure supplement 1e; Mancusso et al., 2012; Nie et al., 2017). To obtain the structure of a DASS protein in a substrate-free, Na+-bound state, we purified VcINDY in 100 mM Na+, but in the absence of a substrate, for single-particle cryo-electron microscopy (cryo-EM). The total mass of a VcINDY dimer is only 96 kDa, and almost that entire mass is embedded in the membrane. This posed a challenge for single-particle analysis, and therefore we used two separate strategies for cryo-EM sample preparation. The first was to increase the effective particle mass using a synthetic antibody fragment (Fab). Using VcINDY reconstituted in lipid nanodiscs, we screened for Fabs using phage display technology (Figure 1—figure supplement 2a; Fellouse et al., 2007; Miller et al., 2012). In this way we identified Fab84, which was subsequently overexpressed, purified, and mixed with VcINDY in nanodiscs (Figure 1—figure supplement 2b and c). The VcINDY-Fab84 structure was determined to 3.15 Å resolution (Figure 1b, Figure 1—figure supplement 2d,e and h, and Table 1). In a parallel approach, we used amphipol polymer (Huynh et al., 2018) to preserve the transporter protein (Figure 1c, Figure 1—figure supplement 2f,g & k, and Table 1), which allowed data collection on particularly thin ice to maximize signal-to-noise. This VcINDY structure in amphipol, without Fab, was determined to 3.16 Å resolution. Both maps are at a sufficiently high resolution and quality to allow direct model building, and the structures are nearly identical (r.m.s.d. 0.750 Å). Finally, to obtain clear density of substrate in the binding site of VcINDY, we crystallized the protein in complex with sodium and terephthalate, and solved the X-ray structure to 3.92 Å resolution by molecular replacement (Figure 1d and Table 2). With two carboxylate moieties at approximately the same distance as those of native substrates, we hypothesized the terephthalate would bind within the VcINDY binding site, while its large benzene ring would provide stronger electron density than succinate, fumarate, or malate.
To obtain the structure of a DASS protein in its Co conformation, we screened various VcINDY homologs and chose a DASS-E protein from Lactobacillus acidophilus (LaINDY) (Figure 1—figure supplement 1c). To characterize the function of LaINDY we performed in vivo complementation and transport assays. LaINDY was able to support aerobic growth on αKG as the sole carbon source for an E. coli strain in which the only aerobically expressed αKG transporter was knocked out (Baba et al., 2006; Figure 2—figure supplement 1a). E. coli transformed with LaINDY accumulated radioactive succinate from an external solution, likely in exchange for endogenous internal dicarboxylates (Figure 2—figure supplement 1b). Subsequent addition of excess external succinate or αKG reduced the amount of internalized radioactive succinate due to LaINDY-facilitated dicarboxylate exchange (Figure 2a).
Like other DASS proteins, purified LaINDY formed a dimer in detergent solution (Figure 2—figure supplement 1d & e). Observing that the VcINDY structures were nearly identical in amphipol and nanodiscs, we surmised that amphipol was suitable for structure determination of another DASS protein, LaINDY. The LaINDY map obtained in this way was at 3.09 Å resolution, allowing for direct model building (Figure 2b, Figure 2—figure supplement 1g–i, and Table 1). Notably, the LaINDY structure was in an outward-facing Co state.
Next, we attempted to determine the structure of LaINDY in its substrate-bound state. To identify substrates suitable for structure determination of such a complex, we examined the thermostability and monodispersity of LaINDY at elevated temperatures in the presence of various potential substrates and substrate analogs (Figure 2—figure supplement 1f). Following identification of malate and αKG as stabilizers, we determined cryo-EM maps of LaINDY in complex with each to 3.36 Å and 3.71 Å, respectively (Figure 2c & d, Figure 2—figure supplement 1j & k, and Table 1). We also solved a 2.85 Å X-ray structure of LaINDY in the presence of both malate and αKG (Figure 2e and Table 2), using the LaINDY-αKG cryo-EM structure as the search model for molecular replacement. As with the apo LaINDY structure, all three LaINDY complex structures were in an outward-facing Co-S state.
The LaINDY and VcINDY maps were all of sufficient quality to build side chains (Figure 2—figure supplement 2a–g), though side-chain rotamers were not always clear at these resolutions. However, we recognized that the significance of comparing structures would be dependent upon the models’ accuracy. This accuracy is a concern particularly as two structural methods and multiple instruments were used. We therefore used the models’ Cα to C bond lengths as an internal benchmark, as we expect this bond to be insensitive to variations in amino acid sequence, secondary structure, or local environment. The Cα – C bond lengths in our cryo-EM structures, as well as those in the PDB, were typically smaller than the ideal and those determined by X-ray crystallography (Figure 2—figure supplement 2h–j). This discrepancy was systematic and independent of map resolution, suggesting it is not the result of microscope calibration errors. Still, the cryo-EM models’ accuracy are sufficient to allow us to interpret the observed side chain movements of 1–2 Å.
Finally, to characterize the Co-S to Ci-S transition of a DASS protein, we used a custom biased MD simulation protocol (Figure 6—figure supplement 1a and Supplementary Note) to drive LaINDY to a Ci-S target model suggested by LaINDY’s internal inverted repeat topological symmetry (Video 1 and Supplementary Note). We bookended this induced transition with unbiased MD simulations of LaINDY’s Co-S and Ci-S states (Video 2), which enabled us to characterize LaINDY’s equilibrium structural dynamics.
Along with the previous VcINDY structures (Mancusso et al., 2012; Nie et al., 2017), the newly determined structures of the Co and Ci conformations in apo and substrate bound states, along with the MD simulations of the Co to Ci transition, allow us to examine the reaction cycle of DASS transporters. We will begin by describing the outward-facing apo state of LaINDY, and subsequently characterize the structural changes associated with substrate binding to the Co state, the Co to Ci transition that carries substrate across the membrane and, finally, substrate release into the cytosol (Figure 1a, and Figure 1—figure supplement 1a and b).
LaINDY is in a Co state
The 3.09 Å cryo-EM map of LaINDY determined in the absence of substrate shows the transporter in its Co apo state (Figure 3a and b). In agreement with LaINDY’s apparent mass in detergent solution (Figure 2—figure supplement 1e), the map shows a transporter dimer. Each protomer consists of a scaffold domain and a transport domain. The transmembrane topology and domain organization of LaINDY resemble that of VcINDY (Mancusso et al., 2012), with the scaffold domain being formed by transmembrane α-helices TMs 1–4 and 7–9, while the transport domain consists of TMs 5, 6, 10 and 11, as well as the helix hairpins HPin and HPout (Figure 3a-c). However, for both hairpins, a bend is found in the second helix, at Val163 in HPinb and Ala404 in HPoutb. Another new structural feature of LaINDY is an extrusion near the dimer interface into the periplasmic space formed by a sequence insertion between TM3 and TM4.
The scaffold and transport domains are linked via two horizontal helices and two loops near the membrane surface, H4c and L4-HPin on the cytosolic side and H9c and L9-HPout on the extracellular side (Figures 3c and 4a). A third linker is formed by a cytoplasmic helix between TM6 and TM7 (H6b), where the equivalent region in VcINDY exists as a long loop (Figure 4b). The domain interface is largely formed by branched or short hydrophobic residues, with only two hydrogen bonds. This results in a smooth domain interface, similar to that observed in the elevator transporter GltPh (Reyes et al., 2009).
While LaINDY’s topology is similar to that of VcINDY, the relative domain positions are different. The transport domain is oriented toward the extracellular side with its substrate binding site facing the periplasm, yielding a Co conformation (Figure 3a). The individual LaINDY and VcINDY domains exhibit strong structural homology, with backbone r.m.s.d.s of 2.897 Å and 2.044 Å for the scaffold and transport domains, respectively (Figure 3d and e). Compared with VcINDY, the transport domain in LaINDY is repositioned by 13.0 Å towards the periplasm with a 37.4° rotation (Figure 4). The domains’ structural conservation and relative positions agree with the notion that DASS proteins operate via a rigid-body, elevator-type movement of the transport domain (Figure 1a).
LaINDY has Na+ surrogate side-chains near the substrate binding site
The structures of LaINDY determined in malate, αKG and in the malate/αKG mixture are in an outward-facing, substrate-bound (Co-S) state (Figure 2e and Figure 2—figure supplement 1j & k). The architecture of the LaINDY binding site is similar to that of VcINDY. Notably, within the binding site, density corresponding to substrate is seen in the X-ray omit map of LaINDY-malate-αKG, and in the LaINDY-αKG cryo-EM difference map, in a similar mode to substrate binding in VcINDY (Figure 5—figure supplement 1e & f). All three Co-S structures are not only similar to each other but, importantly, also to the LaINDY Co apo structure, with pairwise r.m.s.d.s of 0.41–0.52 Å (Figure 3f), indicating that substrate binding introduces little conformational change.
In contrast to the similarity at the substrate binding site, major differences between LaINDY and VcINDY are seen at the cation binding sites. At the Na1 site of VcINDY (Figure 5—figure supplement 1f), an arginine (Arg159) is found in LaINDY, stabilized by a salt bridge with Glu146 (Figure 5a and Figure 5—figure supplement 1a & c). Similarly, a histidine (His392) in LaINDY is located at the equivalent of the VcINDY Na2 site, with another histidine (His401) located 4 Å on its extracellular side. Indeed, these residues are conserved in DASS exchangers but absent in cotransporters (Figure 1—figure supplement 1c). Furthermore, the cation binding sites are more completely enclosed by the surrounding loops in LaINDY compared to VcINDY, facilitated by DASS exchanger specific insertions in the HPin and HPout and L10ab regions (Figure 1—figure supplement 1d). Based on these well enclosed and sterically occupied cation binding sites, it is unlikely that sodium ions bind to LaINDY as Na1 and Na2 do in VcINDY. Rather, we hypothesize these two positively-charged residues in exchangers act as permanent surrogates of the Na+ ions in cotransporters. Such Na+ ion substitutions have previously been observed in other transporters, which makes their substrate transport independent of sodium (Shaffer et al., 2009; Kalayil et al., 2013).
The transport domain moves as a rigid body within the scaffold domain’s arms
The structure determination of LaINDY and VcINDY in multiple states provides an opportunity to characterize conformational changes of a DASS protein during the Co-S to Ci-S transition. MD simulations of succinate-bound LaINDY revealed how the transition between its Co-S and Ci-S states is realized. The transition consists of a 39° rigid-body rotation and an 8.3 Å translation of the transport domain relative to the scaffold domain (Figure 6—figure supplement 1b–d and Video 2). The two ‘arm’ helices, H4c and H9c, are fixed in space with respect to the scaffold domain and cradle the transport domain during the Co-S to Ci-S transition (Figure 6c–e and Video 3). Such arm rigidity agrees with the conserved salt bridge, between Arg122 and Glu283, and bulky-residue interactions at the elbows connecting H4c and H9c to the scaffold domain (Figure 6a & b). The importance of the Arg122 to Glu283 salt bridge is also consistent with recent observations in NaCT, where transport activity was abolished by mutations of the equivalent arginine (Khamaysi et al., 2020), probably by disrupting the conserved salt bridge between H4c and TM7. In contrast to the rigidity at the elbows, flexibility of the hinge loops L4-HPin and L9-HPout at the other end of the arm helices allows transport domain movement (Figure 4). During the Co-S to Ci-S transition the angle between the arm helices and the hairpin helices changes by approximately 30° at both hinges (Figure 6f), allowing the overall rotation and translation of the transport domain.
Despite the large domain movements of LaINDY during MD simulations of the Co-S to Ci-S transition, succinate stably binds without significant changes within the binding site (Figure 5b), though the depth of substrate binding correlates with the orientation of Asn156 (Figure 6—figure supplement 1g,h & l). Furthermore, the conformations and interactions of the Na+ surrogate side chains, Arg159 and His392, are generally stable through the conformational transition (Figure 5—figure supplement 1b & d).
Substrate release from VcINDY causes significant structural changes
The two cryo-EM structures of VcINDY are in a substrate-free, inward-facing (Ci-Na+) state (Figure 1b & c) and provide an opportunity to describe the conformational changes that occur upon substrate release. The Ci-Na+ VcINDY structures are very similar to each other, with an r.m.s.d of 0.75 Å (Figure 7—figure supplement 1b). However, compared with the Ci-Na+-S structures, the substrate-free structures display local changes throughout the protein.
At the substrate binding site of both cryo-EM VcINDY structures, Pro422 at the N-terminus of TM10b in the substrate free structures moved by up to 1.5 Å (Figure 7a). Concurrently, on the cytoplasmic surface, His432 at the C-terminus of the same helix is rotated by 73° (Figure 7—figure supplement 1a). This movement of His432 causes a steric clash with its neighbor, Tyr178 of TM5a, inducing a rotation in that side chain by 41°. One helix turn away, the side chain of Arg175 moves closer to Glu437 of TM11, forming a salt bridge. These rearrangements lead to movements of HPina, HPinb and TM5a by 1.4 Å toward the scaffold domain. Supporting the importance of these interactions, human NaCT’s transport activity is abolished when the equivalent of Glu437 is mutated to histidine, resulting in SLC13A5 Deficiency (Hardies et al., 2015).
Comparison of the previous Ci-Na+-S (Mancusso et al., 2012; Nie et al., 2017) and amphipol preserved Ci-Na+ VcINDY structures also revealed prominent changes on the periplasmic surface (Figure 7b and Video 4). The loop connecting HPoutb and TM10a, from Ala395 to Pro400, has moved on the periplasmic surface. The C-terminus of HPoutb also unwinds by one turn (Val392 – Glu394) in the amphipol structure. As a result, the conserved salt bridge between Glu394 with Lys337 of H9b from the scaffold domain breaks, and Phe396 moves away from its contact with H9c. Such structural changes agree with observations of human NaCT mutations, at positions equivalent to Pro400 and Val401 in VcINDY, which abolish transport and cause SLC13A5 Deficiency (Thevenon et al., 2014; Klotz et al., 2016). Finally, the side chain of Trp461 is inserted between TM6, TM10a and TM11, packing against another conserved aromatic residue, Phe220. Notably, when bound to Fab84, the loop connecting HPoutb and TM10a in VcINDY is midway between the Ci-Na+-S and Ci-Na+ structures. This agrees with Fab84’s epitope covering the periplasmic surface of VcINDY and binding both apo and substrate-bound states (Figure 7—figure supplement 1c & d). Also, a complete lipid molecule from the periplasmic leaflet was found at the interface between the transport and scaffold domains, interacting with TM1, TM2 and HPouta, and therefore may be involved in regulating conformational changes (Figure 7—figure supplement 1e).
Discussion
In this work, we report cryo-EM and X-ray structures of two DASS proteins in four states as well as MD simulations of the Co-S to Ci-S transition. Of these structures, the Ci-Na+, Co, and Co-S are previously un-observed states of the DASS transport cycle. Together with the previously published structures of VcINDY in its Ci-Na+-S state, these results give us a much clearer understanding of the structural basis of transport. The gallery of structures allows us to more completely characterize the reaction cycle of DASS transporters. Thereby, we advance the elevator mechanism of the family from a conceptual model into an atomic description of the transport domain’s movement within the framework of the scaffold domain.
The structures of LaINDY that we have determined represent the first outward-facing structures of any DASS family protein. These structures generally agree with a previous model of VcINDY in its Co conformation (Mulligan et al., 2016), proposed based on the inverted-topology structural repeat and cross-linking distance constraints, with an r.m.s.d of 3.7 Å for the backbone atoms. As noted, the occupation of the Na1 and Na2 sites by conserved basic residues, and absence of any other apparent sodium densities, suggest LaINDY is a DASS exchanger. This is supported by its phylogeny and ability to catalyze succinate-dicarboxylate exchange. However, further experiments, preferably in reconstituted proteoliposomes, will be needed to examine the sodium and proton dependence of transport, and confirm strict substrate coupling in the exchange reaction.
Comparison of the outward-facing and the inward-facing structures, along with MD simulation results, immediately suggests how a DASS protein operates through an elevator-type movement of the transport domain within each protomer (Figure 6c & d). The transport domain moves within the framework formed by the two horizontal α-helix arms on opposing membrane surfaces during the reaction cycle, alternating the substrate binding site between the two sides of the membrane. This mechanism is similar to that proposed for the glutamate transporter GltPh, although the two arms of GltPh are not on the membrane surface but rather transmembrane helices (Reyes et al., 2009).
Analyzing the various conformations also enables us to suggest how substrate binding leads to transporter conformational changes while preventing slippage, or unproductive Co to Ci transitions. In this regard, DASS cotransporters and exchangers appear to employ both unique and shared mechanisms.
Previous experimental data support that Na+-driven DASS cotransporters operate via an ordered sequence (Wright et al., 1983; Yao and Pajor, 2000; Hall and Pajor, 2005; Pajor et al., 2013; Mulligan et al., 2014), namely, Na+ binding induces substrate binding, while substrate release precedes Na+ release. For VcINDY, we have now observed that substrate release in the cytoplasm induces conformational changes as the transporter transitions between the Ci-Na+-S and Ci-Na+ states. Specifically, substrate release leads to an unwinding of the C-terminus of TM11, loop movements, and side chain rotations, resulting in significant changes in local helix packing and protein compactness (Figure 7). This is distinguished from GltPh and homologs in which one or two hairpin gates directly pack against the scaffold domain and block unproductive conformational changes (Reyes et al., 2009; Garaeva et al., 2019; Arkhipova et al., 2020).
In contrast, the exchanger LaINDY exhibits no major conformational changes between the apo and substrate-bound outward-facing states. As two positively-charged residues are found to occupy the cation binding sites, we propose that the DASS exchangers limit slippage, or unproductive conformational changes, via a charge compensation mechanism. In the apo state, where the binding sites have a net positive charge, the hydrophobic surface of the scaffold domain would be an electrostatic barrier to the movement of the transport domain. Only when divalent substrate binds and the net charge is neutralized can the transporter freely exchange between the Ci and Co conformations, ensuring a one-to-one stoichiometric exchange of substrates.
In addition to the noted substrate-release induced conformational changes, charge compensation is also essential to avoid slippage, ensuring Na+-substrate coupling in DASS cotransporters. In fact, the Na+ ions at the Na1 and Na2 sites in cotransporters can be regarded as equivalent to the cationic side chains of exchangers, though Na+ reversibly binds. In cotransporters, the charge compensation model predicts that the transporter can only transition between Co and Ci conformations when the transport domain is either fully-loaded or fully-unloaded. Such a mechanism ensures tightly coupled import of Na+ and substrate in cotransporters, while the reversible binding of Na+ allows for the concentration of the divalent substrate against its electrochemical gradient. Similar charge compensation mechanisms have been proposed for the citrate transporter CitS and the glutamate transporter EAAC1 (Lolkema and Slotboom, 2017; Grewer et al., 2012).
Materials and methods
LaINDY transport activity assays in whole cells
Request a detailed protocolAs a target seed VcINDY was nominated to the cloning core of the New York Consortium of Membrane Protein Structure for the cloning of its homologs (Love et al., 2010). The homologous protein from Lactobacillus acidophilus (Uniprot: Q5FKK5_LACAC, LaINDY) was found to give the highest expression levels in E. coli BL21 DE3 cells when transformed with pET-LaINDY. The transport activity of LaINDY was characterized in E. coli whole cells following a published protocol with minor modifications (Pos et al., 1998; Kim and Unden, 2007). E. coli BL21 DE3 cells were transformed with a modified pET vector (Love et al., 2010) encoding N-terminal 10x His-tagged LaINDY (pET-LaINDY) and grown at 37°C until cells reached OD595 of 0.7. Cells were induced with IPTG and growth was continued for 3 hr. Cells were harvested by centrifugation and resuspended at OD595 = 10 in 10 mM NaCl, 100 mM choline chloride, 50 mM Tris pH 7.5. Cells were kept on ice until needed, and warmed to 30°C five mins prior to transport assay. The uptake reaction was initiated by addition of reaction buffer 10 mM NaCl, 100 mM choline chloride, 50 mM Tris pH 7.5, 1 μM 3H-succinate to the cell suspension at 1:10 (reaction buffer:cells) volumetric ratio. For the exchange reaction, 10 mM succinate or 10 mM α-ketoglutarate (αKG) was added 90 s after initiation of the uptake reaction. Aliquots were collected at fixed time points, with the reaction terminated by collecting the cells on pre-wetted 0.45 mm nitrocellulose filters mounted on a Hoeffer vacuum manifold and immediately washing with 4 mL ice cold 100 mM potassium phosphate buffer pH 7.5. The filters were incubated for 10 mins in scintillation fluid before measuring radioactivity using a Wallac 1450 Microbeta Plus liquid scintillation counter (Shelton, CT) (Law et al., 2007; Law et al., 2008; Law et al., 2009).
E. coli growth assays with LaINDY complementation
Request a detailed protocolThe E. coli strain JW2571 (Baba et al., 2006), in which the only endogenous αKG transporter under aerobic conditions KgtP (Seol and Shatkin, 1991) was knocked out, was transformed with a pGEM-5Zf(+) vector encoding LaINDY. Transformed cells were grown in LB broth with ampicillin to an OD595 ~0.8, then diluted to an OD595 of 0.1. Dilute cells were inoculated into eM9 media supplemented with 50 mM αKG (Rhie et al., 2014), and cell growth monitored using a Tecan SPECTRAFluor Plus microplate reader (Männedorf, Switzerland) incubated at 37°C.
LaINDY expression and purification
Request a detailed protocolLaINDY was expressed by autoinduction (Studier and Moffatt, 1986) at 25°C overnight in the E. coli strain BL21 DE3 transformed with the pET-LaINDY plasmid. Cells were harvested and lysed in a buffer of 50 mM Tris pH 8.0, 400 mM NaCl, 10 mM Imidazole. Ligand, either 10 mM Na+ αKG or 10 mM Na+ DL-malate, was added to the lysis buffer and all subsequent purification steps. Membranes were resuspended in a buffer of 50 mM Tris pH 8.0, 200 mM NaCl, 10 mM Imidazole, and solubilized in 1.2% dodecyl-maltoside (DDM), and protein was purified on a Ni2+-NTA affinity column. The decahistidine tag was removed by overnight digestion at 25°C with TEV protease, followed by preparative size exclusion chromatography (SEC) in a buffer containing 25 mM Tris pH 8.0, 150 mM NaCl, 20% glycerol, and 0.075% DDM.
Thermostabilization assay
Request a detailed protocolA theremostability assay was used to search for compounds that stabilized LaINDY (Auer et al., 2001; Mancusso et al., 2011). Size exclusion chromatography purified LaINDY was dialyzed into a buffer of 20 mM Tris pH 8.0, 150 mM NaCl, 10% glycerol, 0.15% DM. Aliquots of 100 mg LaINDY were incubated with 100 mM test compounds at 42°C for 2 hr, and subsequently injected onto a Shodex KW804 analytical SEC column (Thomson, Clear Brook, VA) on HPLC (Shimadzu, Columbia, MA) in a buffer containing 200 mM Na2SO4, 50 mM Tris 7.5, 3 mM NaN3, and 0.05% DDM. The height of the SEC peak for the 4°C control was used as a reference for normalization.
Multi-angle dynamic light scattering
Request a detailed protocolPurified LaINDY sample (50 μL) was injected onto a Shodex KW803 analytical SEC column on a Waters HPLC (Milford, MA) and eluted with the buffer containing 0.05% DDM at a rate of 0.5 mL/min. The mass of the LaINDY protein was determined using a Wyatt miniDAWN TREOS three angle-static light scattering detector (Santa Barbara, CA), a Wyatt Optilab rEX refractive index detector and a Waters 2489 UV absorbance detector (Slotboom et al., 2008; Waight et al., 2010). The differential refractive index (dn/dc) for DDM, 0.128 mL/g, was calculated using the refractive index detector. The size of the protein–detergent conjugate was deconvoluted following the published method (Kendrick et al., 2001), in which contributions from co-purifying lipids were not distinguished from those of the detergent.
VcINDY expression and purification
Request a detailed protocolExpression and purification of VcINDY was carried out according to our previous protocol (Mancusso et al., 2012). Briefly, E. coli BL21-AI cells (Invitrogen) were transformed with a modified pET vector (Love et al., 2010) encoding N-terminal 10x His tagged VcINDY. Cells were grown at 30°C until OD595 reached 0.8, protein expression occurred at 19°C following IPTG induction, and cells were harvested 16 hr post-induction. Cell membranes were solubilized in 1.2% DDM and the protein was purified on a Ni2+-NTA column. VcINDY was further purified by SEC in buffer containing 25 mM Tris pH 8, 100 mM NaCl, 10 mM Na+-succinate, 5% glycerol and 0.075% DDM unless otherwise indicated.
Nanodisc reconstitution
Request a detailed protocolSEC purified His-tagged VcINDY protein was reconstituted into MSP2N2 nanodiscs (Grinkova et al., 2010) at a molar ratio of VcINDY: MSP: lipids of 1 (9 mM): 8 (72 mM): 277 (2.5 mM) (Bayburt and Sligar, 2010; Schuler et al., 2013). E coli polar lipids in chloroform (Avanti) were vacuum-dried and rehydrated in nanodisc buffer containing 20 mM Tris pH 7.4, 100 mM NaCl, 0.5 mM EDTA and 0.1 mM TCEP to a concentration of 20 mg/mL. VcINDY, lipids and MSP2N2 protein were mixed in the nanodisc buffer to a final volume of 500 mL. Subsequently, 600 mg of Bio-Beads were added and incubated at 4° C overnight. The Bio-Beads were then removed and the solution was filtered through a 0.2 mm centrifugal filter. The nanodiscs containing VcINDY were purified by Ni2+-NTA chromatography. The sample was incubated with a 1:4 molar ratio of Fab84 for 1 hr at 4°C before SEC.
Fab development
Request a detailed protocolThe plasmid pFab101, a modified pFab007 (Miller et al., 2012) which includes a N31T mutation in the VH region, was used as the template to construct a Fab phage-display library termed NYC1 containing 1.1 × 1011 sequences, following a previously described design (Miller et al., 2012). Fab library sorting was performed as previously described (Fellouse et al., 2007; Dominik and Kossiakoff, 2015) with minor modifications. In each round, phage solution was prepared in 50 mM Tris HCl buffer pH 7.5 containing 100 mM NaCl, 1% bovine serum albumin and 0.1 mM TCEP with or without 2 mM succinate, and was first incubated with streptavidin-coated magnetic beads harboring biotinylated nanodisc that did not contain an embedded protein (‘empty’ nanodisc). The supernatant of this reaction was incubated with 100 nM VcINDY embedded in biotinylated nanodisc, and phages bound to the VcINDY-nanodisc complex were captured using streptavidin-coated magnetic beads. A total of four rounds of library sorting were performed. Enriched clones were individually tested using phage ELISA (Fellouse et al., 2007; Sidhu et al., 2000).
Fab expression and purification
Request a detailed protocolTarget Fab84 was subcloned into the Fab expression vector Ptac_Fab_accept_tagless (Burioni et al., 1998). After rigidification (Bailey et al., 2018), the subsequent Fab-containing plasmid was transformed into E. coli 55244 cells. Cells were grown in TGB media at 30°C for 22 hr and harvested. Fab84 protein was purified by injection onto a Protein G column (GE healthcare) and eluted with 100 mM glycine-pH 2.7. Protein was collected in 1 mL fractions containing 2 M Tris pH 8. Fractions containing protein were dialyzed against 50 mM sodium acetate pH 5.3 before purification using a Resource S column (GE healthcare).
Amphipol exchange
Request a detailed protocolAfter nickel column purification PMAL-C8 (Anatrace, Maumee, OH) was added to the detergent-purified transporter protein at a 1:5 protein:amphipol weight ratio (Huynh et al., 2018). The mixture was incubated at 4° C overnight with nutating. To remove detergent, Bio-Beads were incubated with sample at a 1:100 detergent:Bio-Beads weight ratio for 2 hr at 4°C with gentle agitation. The Bio-Beads were then removed by centrifugation at 4,500 rpm and the sample was further purified by SEC in buffer containing 25 mM Tris pH 7.5, 100 mM NaCl and 0.1 mM TCEP.
Cryo-EM sample preparation and data collection
Request a detailed protocolAll cryo-EM grids were prepared by applying 3 mL of protein at ~3 mg/mL to a glow-discharged QuantiAuFoil R1.2/1.3 300-mesh grid (Quantifoil) and blotted for 2.5 to 4 s under 100% humidity at 4°C before plunging into liquid ethane using a Mark IV Vitrobot (FEI).
Cryo-EM data of apo LaINDY were acquired on a Titan Krios microscope with a K2 direct electron detector (Gatan), using a GIF-Quantum energy filter with a 15 eV slit width. Leginon (Suloway et al., 2005) was used for ice thickness targeting and automated data collection. Each micrograph was dose-fractioned over 50 frames, with an accumulated dose of 75 e-/Å2.
Cryo-EM data of LaINDY-malate were acquired on a Titan Krios microscope with a K2 direct electron detector, using a GIF-Quantum energy filter with a 20 eV slit width. SerialEM was used for automated data collection (Schorb et al., 2019). Each micrograph was dose-fractioned over 60 frames, with an accumulated dose of 50 e-/Å2.
Cryo-EM data of LaINDY-αKG were acquired on a Talos Arctica microscope with a K3 direct electron detector. Leginon was used for ice thickness targeting and automated data collection (Rice et al., 2018). Each micrograph was dose-fractioned over 56 frames, with an accumulated dose of 46 e-/Å2.
Cryo-EM data of VcINDY-Na+-Fab84 in nanodiscs were acquired on a Talos Arctica microscope with a K2 direct electron detector. SerialEM was used for automated data collection. Each micrograph was dose-fractioned over 50 frames, with an accumulated dose of 44 e-/Å2.
Cryo-EM data of VcINDY-Na+ in amphipol were acquired on a Talos Arctica microscope with a K2 direct electron detector. SerialEM was used for automated data collection. Each micrograph was dose-fractioned over 48 frames, with an accumulated dose of 40 e-/Å2.
Cryo-EM image processing
Request a detailed protocolMotion correction, CTF estimation, particle picking, 2D classification, ab initio model generation, heterogenous and non-uniform refinement, and per particle CTF refinement were all performed with cryoSPARC (Punjani et al., 2017). Each dataset was processed using the same protocol, except as noted. All maps were sharpened using Auto-sharpen Map in Phenix (Adams et al., 2010).
Micrographs underwent motion correction and the CTF estimated, and those with an overall resolution worse than 8 Å excluded from subsequent steps. An ellipse-based particle picker identified particles used to generate initial 2D classes. These classes were used for template-based particle picking. Template identified particles were extracted and subjected to 2D classification. A subset of well resolved 2D classes were used for the initial ab initio model building, while all picked particles were subsequent used for heterogeneous 3D refinement. After multiple rounds of 3D classification (ab initio model generation and heterogeneous 3D refinement with two or more classes), a single class was selected for non-uniform 3D refinement, resulting in the final map.
For the VcINDY-Na+-Fab84 sample in nanodiscs, refinement of an initial 3.36 Å map indicated oblique views were rare in the particle set. Projections of this map were therefore used as templates for an additional round of particle picking, with subsequent 3D classification, per particle CTF refinement and non-uniform 3D refinement.
LaINDY crystallization and X-ray diffraction data collection
Request a detailed protocolLaINDY protein purified in NG and αKG was concentrated to 2 mg/mL using a centrifugal filtration device and crystallized at 18°C using the hanging drop vapor diffusion method by mixing equal volumes of concentrated protein and well solution of 30% Jeffamine ED-2001, 100 mM HEPES pH 6.8. Crystal quality was improved by the addition of 10 mM malic acid prior to size exclusion chromatography. Crystals were flash frozen in liquid nitrogen with crystallization solution serving as cryoprotectant. Crystals were of the P21 space group with unit cell dimensions of a = 91.3 Å, b = 76.6 Å, c = 96.9 Å, β = 90.5° and contained two molecules in the asymmetric unit, which form a single transporter dimer. Data were collected at the Advanced Light Source beamline 5.0.2. Diffraction data were processed and scaled in HKL2000 (Otwinowski and Minor, 1997).
VcINDY crystallization and X-ray diffraction data collection
Request a detailed protocolVcINDY purified in DM was incubated with 10 mM sodium terephthalate (TTP) on ice for 30 min prior to crystallization with a protein concentration of 4 mg/mL at 4°C by hanging-drop vapor diffusion. Crystals were grown in 30% PEG 550, 50 mM NaCl and 100 mM Tris pH 8.0, and were subsequently frozen in liquid nitrogen with the crystallization solution serving as the cryoprotectant. X-ray data were collected at the Advanced Light Source Beamline 5.0.2. Data processing and scaling were performed using HKL2000. Crystals were of space group P21 with unit cell dimensions around a = 108 Å, b = 103 Å, c = 174 Å, β = 96°, and contained four molecules per asymmetric unit.
Model building and refinement
Request a detailed protocolCryo-EM models were built in Coot (Emsley and Cowtan, 2004) and refined in Phenix real space refine (Adams et al., 2010). The LaINDY-apo model was manually built, and subsequently used as an initial model for the LaINDY-malate and LaINDY-αKG structures. VcINDY models were built using the VcINDY X-ray structure (PDB: 5UL9), with sodium and citrate removed, as an initial model. The Fab was rebuilt from an initial model of a homologous X-ray structure (PDB: 5EII) with the variable loops removed. In both VcINDY and LaINDY maps some discontinuous densities could be attributed to the aliphatic chains of lipids, detergents, or amphipols and were modeled as alkanes for simplicity.
The LaINDY X-ray structure was determined by molecular replacement using the LaINDY-apo cryo-EM model as a starting model, followed by model building in Coot and refinement in Phenix.
The VcINDY-TTP X-ray structure was determined by molecular replacement using the structure of VcINDY (PDB ID: 5UL9), with sodium and citrate removed, as the initial search model. The model was refined using group B-factors, and NCS, secondary structure and reference model restraints.
Phylogenetic tree and sequence conservation
Request a detailed protocolDASS family orthologs were aligned in Promals3D (Pei et al., 2008), including VcINDY (PDB ID: 5UL9) and CitS (PDB ID: 5A1S) as structural models. As a non-member of the DASS family, CitS was used as the outgroup when calculating distances using FastTree (Price et al., 2009).
Bond length analysis
Request a detailed protocolThe distance between Cα and C carbons was calculated using a custom Biopython script (Cock et al., 2009). Median Cα to C distances were calculated for 2045 single particle cryo-EM and 64,676 X-ray structures deposited in the Protein Data Bank between 2013 and 2020, containing only L-type polypeptides, and resolved to at least 4 Å.
System preparation
Request a detailed protocolMolecular dynamics (MD) simulation systems of the Co-S and Co states of LaINDY were constructed from its X-ray crystal structure. Succinate was added to the X-ray crystal structure of LaINDY by aligning the cryo-EM structure of LaINDY with α-ketoglutarate bound and mutating the α-ketoglutarate into succinate. Force field parameters for succinate were determined using the online CGenFF interface (Vanommeslaeghe and MacKerell, 2012; Vanommeslaeghe et al., 2012). The protonation states of the titratable residues of LaINDY were determined using PROPKA (Søndergaard et al., 2011; Olsson et al., 2011). Two independent POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) membranes were prepared using the Membrane Builder (Wu et al., 2014; Jo et al., 2009; Jo et al., 2007; Lee et al., 2019) tool of CHARMM-GUI (Jo et al., 2008). The orientation of LaINDY in the membrane was determined by aligning the protein’s symmetry axis with the membrane-normal axis, and the insertion depth in the membrane was determined by applying the PPM (Positioning of Proteins in Membrane) Server (Lomize et al., 2012) to a model (Mulligan et al., 2016) of the Co-Na+-S state of VcINDY. Once the membrane and protein were combined, all lipids within 1 Å of the protein were removed from the systems. Water was added to the systems using the VMD (Humphrey et al., 1996) Solvate plugin. The systems were neutralized with Cl−, and 0.15 M NaCl was added to the systems using the VMD Autoionize plugin.
Simulation parameters
Request a detailed protocolAll simulations were performed with NAMD (Phillips et al., 2005) using the CHARMM36m (Huang et al., 2017; Klauda et al., 2010; Vanommeslaeghe et al., 2010; Yu et al., 2012; Jorgensen et al., 1983; Venable et al., 2013) force field at a constant temperature of 310 K and a constant pressure of 1 atm. Constant temperature was maintained using Langevin dynamics for all non-hydrogen atoms with a damping coefficient of 1 ps−1. Constant pressure was maintained using a Nosé-Hoover Langevin piston (Martyna et al., 1994; Feller et al., 1995) with a period of 100 fs and a damping time scale of 50 fs. The cut-off distance for both electrostatic and van der Waals interactions was set to 12 Å, and a switching function was applied at 10 Å. Periodic boundary conditions were applied in all simulations, and long-range electrostatic interactions were calculated using the particle mesh Ewald method (Darden et al., 1993; Essmann et al., 1995) with a grid point density of 1 Å−1.
Simulations performed
Request a detailed protocolAfter building the systems, harmonic restraints were applied to the constituent atoms of the protein, membrane, and bound succinate (i.e., only the bulk water and ions were left unrestrained). A force constant of 1.0 kcal/mol/Å2 was used, and an integration time step of 1 fs was used. Three thousand steps of conjugate gradient energy minimization were performed, after which the systems were simulated with MD for 1 ns. The restraints on the membrane were then released, and the systems were minimized for three thousand steps and then simulated for 10 ns. Next, the restraints on the protein and bound succinate were released, and the systems were minimized for three thousand steps and then simulated for 10 ns. With initial equilibration complete, the integration time step was changed to 2 fs, and an additional 100 ns of unbiased equilibrium simulation was performed for each system to sample the dynamics of the Co-S and Co states of LaINDY.
Biases were then applied with the goal of inducing a transition from the Co-S state to an approximate target model of the Ci-S state. These biases and the approximate Ci-S target are described in the ‘Biasing Protocol’ section below. To simplify comparison between the transitions sampled by the two unique protomers of LaINDY and between independent trials of the simulation, a 10 ns transition was first induced to a convenient reference Co-S state: the X-ray crystal structure. The system was then restrained at this reference Co-S state for 10 ns to allow the membrane and aqueous environment to equilibrate around this new protein state. Next, a 100 ns transition to the approximate Ci-S target was induced, and the system was subsequently restrained at the approximate Ci-S target for 50 ns again to allow the environment to equilibrate around the new protein state. Finally, all restraints were released, and a 100 ns unbiased equilibrium simulation was performed to examine the stability of and to sample the dynamics of the newly generated model of the Ci-S state of LaINDY. See Figure 6—figure supplement 1a for a graphical overview of the simulations performed.
The transition to the Ci-S state was first induced without any biases applied to the bound succinate. A successful transition was sampled by protomer B, but succinate spontaneously unbound during the early stages of the transition sampled by protomer A. To increase sampling, an additional trial of this transition was performed with biases applied to the bound succinate (see ‘Biasing Protocol’). Throughout this article, data from the simulation without biases applied to succinate are labeled ‘protomer B1.’ Data from the simulation with biases applied to succinate are labeled ‘protomer A2’ and ‘protomer B2.’.
Biasing protocol
Request a detailed protocolIn transporters with an inverted repeat topology, it has been shown that the Co-S and Ci-S states can be generated from each other by swapping the conformations of the transporter’s repeats (Forrest et al., 2011). The internal symmetry in these transporters means that the fundamental global structural differences between their Co-S and Ci-S states can be observed simply by rotating one state 180° about any axis in the plane of the membrane. When this approach is applied to LaINDY, it is immediately apparent that the fundamental global differences between the Co-S X-ray crystal structure and the newly generated approximate Ci-S target model are: (1) the orientation of the transport domains relative to the scaffold domain and (2) the position of the transport domains relative to the scaffold domain along the membrane-normal axis. To describe these differences quantitatively and to induce a transition to the approximate Ci-S target, the Collective Variables (Colvars) module (Fiorin et al., 2013) in NAMD was used.
In the Colvars module, collective variables are mathematical functions used to describe collective properties of interest of selections of atoms within a system, and transitions are induced by applying moving harmonic potentials to the collective variables. In this case, one collective variable was defined for each of the fundamental global differences observed between the Co-S state and approximate Ci-S target. Specifically, for the orientations of the transport domains, spinAngle collective variables were used, and distanceZ collective variables were used for the positions of the transport domains relative to the scaffold domain along the membrane-normal axis. Helix Cα atoms of the transport domains and transmembrane helix Cα atoms of the scaffold domain were used in the definitions of these collective variables. Schematic depictions of these collective variables are shown in Figure 6—figure supplement 1b.
To determine the rotation axis used in the definition of the spinAngle collective variables, the rotation matrix needed to structurally align the transport domains of the Co-S state with the transport domains of the approximate Ci-S target was calculated using the measure fit command of VMD. Importantly, the order option was used to allow the first repeat of the Co-S state to be compared to the second repeat of the approximate Ci-S target at the same time that the second repeat of the Co-S state was compared to the first repeat of the approximate Ci-S target. This calculation also determined the change in angle (i.e., amount of spin about the spin axis) that was needed to induce the transition to the approximate Ci-S target. To induce the transition to the approximate Ci-S target using biased MD simulations, moving harmonic restraints with a force constant of 8 kcal/mol/(°)2 were applied to the spinAngle collective variables.
Calculating the change in the distanceZ collective variables required to induce the transition to the approximate Ci-S target was straightforward, as the approach used to generate the approximate Ci-S target simply changes the direction of the displacement along the membrane-normal axis of the transport domains relative to the scaffold domain without affecting its magnitude. Moving harmonic restraints with a force constant of 80 kcal/mol/Å2 were applied to the distanceZ collective variables. Video 1 demonstrates how the approximate Ci-S target was generated, and it shows the VMD-based quantification (described above) of the global structural differences between the Co-S state and the approximate Ci-S target used to define the spinAngle and distanceZ collective variables.
In addition to the spinAngle and distanceZ collective variables, harmonic restraints were applied to up to two additional collective variables during the induced transition. First, the orientation of the scaffold domain was restrained to prevent the overall orientation of LaINDY from changing. Specifically, orientation collective variables were used on the scaffold domains (defined as above), and harmonic restraints with a force constant of 100,000 kcal/mol were applied to them. Second, in one trial of the transition (see ‘Simulations Performed’), succinate was restrained to the binding sites of the transport domains using distance collective variables. One distance collective variable per succinate was defined using the non-hydrogen atoms of the succinate and the Cα atoms of the residues that were consistently closest to the bound succinate during the unbiased simulations of the Co-S state of LaINDY (i.e., Pro154-Arg159, Thr208-Pro212, Ala396-Thr399, and Asn442-Pro444). A half-harmonic restraint centered at 2.0 Å and with a force constant of 1.0 kcal/mol/Å2 was applied to these collective variables.
Data and code availability
Request a detailed protocolCryo-EM maps and models have been deposited in the Protein Data Bank and EMDB database, respectively, for VcINDY-Na+ in amphipol (6WU3, EMD-21904), VcINDY-Na+-Fab84 in nanodisc (6WW5, EMD-21928), LaINDY-apo (6WU1, EMD-21902), LaINDY-αKG (6WU4, EMD-21905), and LaINDY-malate (6WU2, EMD-21903). X-ray derived models and diffraction data have been deposited in the Protein Data Bank for LaINDY-malate-αKG (6WTW) and VcINDY-TTP (6WTX). Coordinates of representative LaINDY structures from MD simulations of the Co-S and Ci-S states have been made publicly available on Zenodo (DOI 10.5281/zenodo.3965996). Bond length analysis code is available at Sauer, 2020; https://github.com/DavidBSauer/bond_length_analysis (copy archived at https://github.com/elifesciences-publications/bond_length_analysis).
Appendix 1
Supplementary Note
Molecular Dynamics Simulations
Using MD, we simulated the unbiased equilibrium dynamics of the Co-S state of LaINDY, we induced a transition to an approximate target model of the Ci-S state using a custom biased MD simulation protocol, and we simulated the equilibrium dynamics of the new putative Ci-S state (Video 2). In the custom biased MD simulation protocol, biases were applied to collective variables representing the displacement z and orientation θ of the transport domains of LaINDY relative to its scaffold domain (Figure 6—figure supplement 1b). From the ensembles generated by the unbiased equilibrium MD simulations of the Co-S and putative Ci-S states, we selected as representative the structures with the most representative values of z and θ concurrently in both protomers (Figure 6—figure supplement 1i). Water accessibility profiles of these structures demonstrate that the new state is definitively a Ci-S state (Figure 6—figure supplement 1j), and z and θ time series indicate that the new Ci-S state is stable (Figure 6—figure supplement 1c&d).
To generate a model of a transporter’s missing state based on its internal inverted repeat topological symmetry, standard procedure is to apply homology modeling using the two repeats as structural templates for each other (Vergara-Jaque et al., 2015). To generate our model of the Ci-S state of LaINDY, we instead used a novel MD-based approach that only indirectly drives a swap of the conformations of LaINDY’s inverted repeats (Figure 6—figure supplement 1e,f &k), and we favor this approach for several reasons. Most importantly, the model is generated in the presence of its natural membrane environment, and confirmation of its stability in this environment via MD is built in to the approach. This approach also enables the transporter dynamics to be assessed before, after, and during the transition between states, and it lays the groundwork needed to apply an advanced array of MD-based techniques to characterize detailed structural and thermodynamic properties of the transition (Moradi and Tajkhorshid, 2013; Moradi and Tajkhorshid, 2014; Moradi et al., 2015).
Because the structure of VcINDY’s Ci-Na+-S state has been determined, it also would have been possible to generate a model of the Ci-S state of LaINDY using homology modeling with VcINDY as a template. We favor our MD-based approach over this second homology modeling-based approach for the same reasons listed in the previous paragraph. That said, it would have been possible to use information from the Ci-Na+-S state of VcINDY, rather than information from LaINDY’s inverted repeats, to generate Ci-S targets for the biased MD simulations. Within the collective variable space, the targets suggested by LaINDY’s inverted repeats are −4.16 Å for z and −19.5° for θ, and the targets suggested by VcINDY’s Ci-Na+-S state are −8.34 Å for z and −22.1° for θ. Both sets of targets are reasonable, but we ultimately chose to use the targets provided by the inverted repeats because LaINDY and VcINDY, while homologous, are different proteins. The stability of the LaINDY Ci-S model at the z and θ targets (Figure 6—figure supplement 1c&d) and the accessibility of bound succinate from the cytosol (Figure 6—figure supplement 1j) demonstrate the quality of the inverted repeat-based model, but we cannot comment on what the results would have been if the VcINDY-based targets had been used instead.
The biased MD simulations we have performed are a powerful tool for assessing LaINDY’s dynamics during the transition between its Co-S and Ci-S states, but they also have limitations. Importantly, the pathway taken through the collective variable space during the transition has not been optimized in these simulations. As such, we restricted analyses of our biased MD simulations to the global structural differences between the end states and the internal dynamics of the transport domains, which this work suggests move as rigid bodies relative to the scaffold domain. Because restraints were applied to succinate in simulation trials A2 and B2, the internal structural dynamics of the transport domains could have been perturbed in these simulations, but we suggest that this did not happen because of the consistency of the data from these trials with the data from simulation trial B1, in which no restraints were applied to succinate. Additionally, the biased MD simulations we have performed cannot be used to generate free energy profiles associated with the transitions, also limiting the kinds of analyses we were able to perform. To overcome these limitations and expand on this work, an advanced array of MD-based techniques would need to be applied to the transition (Moradi and Tajkhorshid, 2013; Moradi and Tajkhorshid, 2014; Moradi et al., 2015). We again note that the biased MD simulations we have performed lay the groundwork needed to apply these advanced MD-based techniques.
Data availability
Cryo-EM maps and models have been deposited in the Protein Data Bank and EMDB database, respectively, for VcINDY-Na+ in amphipol (6WU3, EMD-21904), VcINDY-Na+-Fab84 in nanodisc (6WW5, EMD-21928), LaINDY-apo (6WU1, EMD-21902), LaINDY-aKG (6WU4, EMD-21905), and LaINDY-malate (6WU2, EMD-21903). X-ray derived models and diffraction data have been deposited in the Protein Data Bank for LaINDY-malate-aKG (6WTW) and VcINDY-TTP (6WTX). Coordinates of representative LaINDY structures from MD simulations of the Co-S and Ci-S states have been made publicly available on Zenodo (DOI: 10.5281/zenodo.3965996). Bond length analysis code is available at https://github.com/DavidBSauer/bond_length_analysis (copy archived at https://github.com/elifesciences-publications/bond_length_analysis).
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
Locking the elbow: improved antibody Fab fragments as chaperones for structure determinationJournal of Molecular Biology 430:337–347.https://doi.org/10.1016/j.jmb.2017.12.012
-
Membrane protein assembly into nanodiscsFEBS Letters 584:1721–1727.https://doi.org/10.1016/j.febslet.2009.10.024
-
SLC13 family of Na⁺-coupled di- and tri-carboxylate/sulfate transportersMolecular Aspects of Medicine 34:299–312.https://doi.org/10.1016/j.mam.2012.12.001
-
A vector for the expression of recombinant monoclonal Fab fragments in bacteriaJournal of Immunological Methods 217:195–199.https://doi.org/10.1016/S0022-1759(98)00102-1
-
Particle mesh Ewald - an N.log(N) method for Ewald sums in large systemsThe Journal of Chemical Physics 98:10089–10092.https://doi.org/10.1063/1.464397
-
Phage display selections for affinity reagents to membrane proteins in nanodiscsMethods in Enzymology 557:219–245.https://doi.org/10.1016/bs.mie.2014.12.032
-
Shared molecular mechanisms of membrane transportersAnnual Review of Biochemistry 85:543–572.https://doi.org/10.1146/annurev-biochem-060815-014520
-
Coot: model-building tools for molecular graphicsActa Crystallographica. Section D, Biological Crystallography 60:2126–2132.https://doi.org/10.1107/S0907444904019158
-
A smooth particle mesh Ewald methodThe Journal of Chemical Physics 103:8577–8593.https://doi.org/10.1063/1.470117
-
Constant pressure molecular dynamics simulation: the Langevin piston methodThe Journal of Chemical Physics 103:4613–4621.https://doi.org/10.1063/1.470648
-
High-throughput generation of synthetic antibodies from highly functional minimalist phage-displayed librariesJournal of Molecular Biology 373:924–940.https://doi.org/10.1016/j.jmb.2007.08.005
-
Using collective variables to drive molecular dynamics simulationsMolecular Physics 111:3345–3362.https://doi.org/10.1080/00268976.2013.813594
-
The structural basis of secondary active transport mechanismsBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1807:167–188.https://doi.org/10.1016/j.bbabio.2010.10.014
-
Elevator-type mechanisms of membrane transportBiochemical Society Transactions 48:1227–1241.https://doi.org/10.1042/BST20200290
-
Charge compensation mechanism of a Na+-coupled, secondary active glutamate transporterThe Journal of Biological Chemistry 287:26921–26931.https://doi.org/10.1074/jbc.M112.364059
-
Engineering extended membrane scaffold proteins for self-assembly of soluble nanoscale lipid bilayersProtein Engineering Design and Selection 23:843–848.https://doi.org/10.1093/protein/gzq060
-
The emergence of 2-oxoglutarate as a master regulator metaboliteMicrobiology and Molecular Biology Reviews 79:419–435.https://doi.org/10.1128/MMBR.00038-15
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
CHARMM-GUI: a web-based graphical user interface for CHARMMJournal of Computational Chemistry 29:1859–1865.https://doi.org/10.1002/jcc.20945
-
Comparison of simple potential functions for simulating liquid waterThe Journal of Chemical Physics 79:926–935.https://doi.org/10.1063/1.445869
-
Online size-exclusion high-performance liquid chromatography light scattering and differential refractometry methods to determine degree of polymer conjugation to proteins and protein-protein or protein-ligand association statesAnalytical Biochemistry 299:136–146.https://doi.org/10.1006/abio.2001.5411
-
A dynamic anchor domain in slc13 transporters controls metabolite transportJournal of Biological Chemistry 295:8155–8163.https://doi.org/10.1074/jbc.RA119.010911
-
The L-tartrate/succinate antiporter TtdT (YgjE) of L-tartrate fermentation in Escherichia coliJournal of Bacteriology 189:1597–1603.https://doi.org/10.1128/JB.01402-06
-
Update of the CHARMM all-atom additive force field for lipids: validation on six lipid typesThe Journal of Physical Chemistry B 114:7830–7843.https://doi.org/10.1021/jp101759q
-
Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpTJournal of Molecular Biology 378:828–839.https://doi.org/10.1016/j.jmb.2008.03.029
-
CHARMM-GUI Membrane Builder for complex biological membrane simulations with glycolipids and lipoglycansJournal of Chemical Theory and Computation 15:775–786.https://doi.org/10.1021/acs.jctc.8b01066
-
Hydropathy profile alignment: a tool to search for structural homologues of membrane proteinsFEMS Microbiology Reviews 22:305–322.https://doi.org/10.1111/j.1574-6976.1998.tb00372.x
-
Structure and elevator mechanism of the Na+-citrate transporter CitSCurrent Opinion in Structural Biology 45:1–9.https://doi.org/10.1016/j.sbi.2016.10.004
-
OPM database and PPM web server: resources for positioning of proteins in membranesNucleic Acids Research 40:D370–D376.https://doi.org/10.1093/nar/gkr703
-
The New York Consortium on Membrane Protein Structure (NYCOMPS): a high-throughput platform for structural genomics of integral membrane proteinsJournal of Structural and Functional Genomics 11:191–199.https://doi.org/10.1007/s10969-010-9094-7
-
The SLC13 gene family of sodium sulphate/carboxylate cotransportersPflugers Archiv European Journal of Physiology 447:594–602.https://doi.org/10.1007/s00424-003-1128-6
-
Constant pressure molecular dynamics algorithmsThe Journal of Chemical Physics 101:4177–4189.https://doi.org/10.1063/1.467468
-
Computational recipe for efficient description of large-scale conformational changes in biomolecular systemsJournal of Chemical Theory and Computation 10:2866–2880.https://doi.org/10.1021/ct5002285
-
Functional characterization of a Na+-dependent dicarboxylate transporter from Vibrio choleraeThe Journal of General Physiology 143:745–759.https://doi.org/10.1085/jgp.201311141
-
The bacterial dicarboxylate transporter VcINDY uses a two-domain elevator-type mechanismNature Structural & Molecular Biology 23:256–263.https://doi.org/10.1038/nsmb.3166
-
PROPKA3: consistent treatment of internal and surface residues in empirical pKa predictionsJournal of Chemical Theory and Computation 7:525–537.https://doi.org/10.1021/ct100578z
-
Processing of X-ray diffraction data collected in oscillation modeMethods in Enzymology 276:307–326.https://doi.org/10.1016/S0076-6879(97)76066-X
-
Functional characterization of SdcF from Bacillus licheniformis, a homolog of the SLC13 Na+/dicarboxylate transportersThe Journal of Membrane Biology 246:705–715.https://doi.org/10.1007/s00232-013-9590-3
-
Sodium-coupled dicarboxylate and citrate transporters from the SLC13 familyPflügers Archiv - European Journal of Physiology 466:119–130.https://doi.org/10.1007/s00424-013-1369-y
-
PROMALS3D: a tool for multiple protein sequence and structure alignmentsNucleic Acids Research 36:2295–2300.https://doi.org/10.1093/nar/gkn072
-
UCSF Chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
The ion transporter superfamilyBiochimica Et Biophysica Acta (BBA) - Biomembranes 1618:79–92.https://doi.org/10.1016/j.bbamem.2003.10.010
-
FastTree: computing large minimum evolution trees with profiles instead of a distance matrixMolecular Biology and Evolution 26:1641–1650.https://doi.org/10.1093/molbev/msp077
-
Routine determination of ice thickness for cryo-EM gridsJournal of Structural Biology 204:38–44.https://doi.org/10.1016/j.jsb.2018.06.007
-
Nanodiscs as a new tool to examine lipid-protein interactionsMethods in Molecular Biology 974:415–433.https://doi.org/10.1007/978-1-62703-275-9_18
-
Phage display for selection of novel binding peptidesMethods in Enzymology 328:333–363.https://doi.org/10.1016/s0076-6879(00)28406-1
-
Improved treatment of ligands and coupling effects in empirical calculation and rationalization of pKa valuesJournal of Chemical Theory and Computation 7:2284–2295.https://doi.org/10.1021/ct200133y
-
Use of bacteriophage T7 RNA polymerase to direct selective high-level expression of cloned genesJournal of Molecular Biology 189:113–130.https://doi.org/10.1016/0022-2836(86)90385-2
-
Automated molecular microscopy: the new leginon systemJournal of Structural Biology 151:41–60.https://doi.org/10.1016/j.jsb.2005.03.010
-
Mutations in SLC13A5 cause autosomal-recessive epileptic encephalopathy with seizure onset in the first days of lifeThe American Journal of Human Genetics 95:113–120.https://doi.org/10.1016/j.ajhg.2014.06.006
-
CHARMM general force field: a force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fieldsJournal of Computational Chemistry 31:671–690.https://doi.org/10.1002/jcc.21367
-
Automation of the CHARMM general force field (CGenFF) II: assignment of bonded parameters and partial atomic chargesJournal of Chemical Information and Modeling 52:3155–3168.https://doi.org/10.1021/ci3003649
-
Automation of the CHARMM general force field (CGenFF) I: bond perception and atom typingJournal of Chemical Information and Modeling 52:3144–3154.https://doi.org/10.1021/ci300363c
-
Simulations of anionic lipid membranes: development of interaction-specific ion parameters and validation using NMR dataThe Journal of Physical Chemistry B 117:10183–10192.https://doi.org/10.1021/jp401512z
-
Structure and mechanism of a pentameric formate channelNature Structural & Molecular Biology 17:31–37.https://doi.org/10.1038/nsmb.1740
-
Kinetics of sodium succinate cotransport across renal brush-border membranesThe Journal of Biological Chemistry 258:5456–5462.
-
CHARMM-GUI Membrane Builder toward realistic biological membrane simulationsJournal of Computational Chemistry 35:1997–2004.https://doi.org/10.1002/jcc.23702
-
The transport properties of the human renal Na+-dicarboxylate cotransporter under voltage-clamp conditionsAmerican Journal of Physiology-Renal Physiology 279:F54–F64.https://doi.org/10.1152/ajprenal.2000.279.1.F54
-
Extension of the CHARMM general force field to sulfonyl-containing compounds and its utility in biomolecular simulationsJournal of Computational Chemistry 33:2451–2468.https://doi.org/10.1002/jcc.23067
Article and author information
Author details
Funding
National Institutes of Health (R01NS108151)
- Da-Neng Wang
National Institutes of Health (R01GM121994)
- Da-Neng Wang
National Institutes of Health (R01DK099023)
- Da-Neng Wang
National Institutes of Health (U54GM095315)
- Da-Neng Wang
National Institutes of Health (P41GM104601)
- Emad Tajkhorshid
National Institutes of Health (R01GM067887)
- Emad Tajkhorshid
TESS Research Foundation
- Da-Neng Wang
American Epilepsy Society (AES2017SD3)
- Da-Neng Wang
American Cancer Society (129844-PF-17-135-01-TBE)
- David B Sauer
Department of Defense (W81XWH-16-1-0153)
- David B Sauer
National Science Foundation (1746047)
- Noah Trebesch
National Institutes of Health (T32GM088118)
- Nicolette Cocco
XSEDE (TG-MCA06N060)
- Emad Tajkhorshid
National Institutes of Health (R01GM123455)
- Emad Tajkhorshid
Blue Waters
- Emad Tajkhorshid
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was financially supported by the NIH (R01NS108151, R01GM121994 and R01DK099023 to DNW; U54GM095315 to WA Hendrickson; P41GM104601, R01GM067887 and R01GM123455 to ET), the TESS Research Foundation and the American Epilepsy Society (to DNW). DBS was supported by the American Cancer Society Postdoctoral Fellowship (129844-PF-17-135-01-TBE) and Department of Defense Horizon Award (W81XWH-16-1-0153). NT was supported by a NSF Graduate Research Fellowship (1746047). NC was supported by an NIH Predoctoral Training Grant (T32-GM088118). VcINDY homologs were cloned by B Kloss at the New York Consortium of Membrane Protein Structures. We thank the following colleagues for reagent, technical assistance, and helpful discussions: JP Armache, JG Belasco, N Coudray, T Hattori, J Jiang, NK Karpowich, M Lopez Redondo, Z Liu, R Mancusso, AB Rejto and SG Sligar. We are also grateful to the staff at the following facilities for assistance in screening and data collection in cryo-EM and X-ray diffraction: the NYU Cryo-EM Facility, the Pacific Northwest Center for Cryo-EM, Advanced Light Source Beamline 5.0.2 at the Berkeley National Laboratory, beamlines AMX and FMX at NSLS-II and 19-BM and 19-ID at the Advanced Photon Source. EM data processing used computing resources at the HPC Facility of NYULMC, and we were assisted by A Siavosh-Haghighi and M Costantino. MD simulations were performed using supercomputing allocations from Blue Waters and XSEDE (TG-MCA06N060) to ET. Blue Waters is supported by the NSF (OCI-0725070 and ACI-1238993), the State of Illinois and the National Geospatial-Intelligence Agency. XSEDE is supported by the NSF (ACI-1548562).
Copyright
© 2020, Sauer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,475
- views
-
- 549
- downloads
-
- 48
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
The low-complexity domain of hnRNPA1 (A1-LCD) phase separates in a salt-dependent manner. Unlike many intrinsically disordered proteins (IDPs) whose phase separation is suppressed by increasing salt concentrations, the phase separation of A1-LCD is promoted by >100 mM NaCl. To investigate the atypical salt effect on A1-LCD phase separation, we carried out all-atom molecular dynamics simulations of systems comprising multiple A1-LCD chains at NaCl concentrations from 50 to 1000 mM NaCl. The ions occupy first shell as well as more distant sites around the IDP chains, with Arg sidechains and backbone carbonyls the favored partners of Cl– and Na+, respectively. They play two direct roles in driving A1-LCD condensation. The first is to neutralize the high net charge of the protein (+9) by an excess of bound Cl– over Na+; the second is to bridge between A1-LCD chains, thereby fortifying the intermolecular interaction networks in the dense phase. At high concentrations, NaCl also indirectly strengthens π–π, cation–π, and amino–π interactions, by drawing water away from the interaction partners. Therefore, at low salt, A1-LCD is prevented from phase separation by net charge repulsion; at intermediate concentrations, NaCl neutralizes enough of the net charge while also bridging IDP chains to drive phase separation. This drive becomes even stronger at high salt due to strengthened π-type interactions. Based on this understanding, four classes of salt dependence of IDP phase separation can be predicted from amino-acid composition.
-
- Structural Biology and Molecular Biophysics
Many voltage-gated potassium (Kv) channels display a time-dependent phenomenon called C-type inactivation, whereby prolonged activation by voltage leads to the inhibition of ionic conduction, a process that involves a conformational change at the selectivity filter toward a non-conductive state. Recently, a high-resolution structure of a strongly inactivated triple-mutant channel kv1.2-kv2.1-3m revealed a novel conformation of the selectivity filter that is dilated at its outer end, distinct from the well-characterized conductive state. While the experimental structure was interpreted as the elusive non-conductive state, our molecular dynamics simulations and electrophysiological measurements show that the dilated filter of kv1.2-kv2.1-3m is conductive and, as such, cannot completely account for the inactivation of the channel observed in the structural experiments. The simulation shows that an additional conformational change, implicating isoleucine residues at position 398 along the pore lining segment S6, is required to effectively block ion conduction. The I398 residues from the four subunits act as a state-dependent hydrophobic gate located immediately beneath the selectivity filter. By mutating I398 to Asparagine, ion permeation can be resumed in the kv1.2-kv2.1-3m channel, which was not a reversion C-type inactivation, since AgTxII fails to block the ionic permeation of kv1.2-kv2.1-3m_I398N. As a critical piece of the C-type inactivation machinery, this structural feature is the potential target of a broad class of QA blockers and negatively charged activators thus opening new research directions towards the development of drugs that specifically modulate gating-states of Kv channels.