Molecular basis for allosteric specificity regulation in class Ia ribonucleotide reductase from Escherichia coli
Abstract
Ribonucleotide reductase (RNR) converts ribonucleotides to deoxyribonucleotides, a reaction that is essential for DNA biosynthesis and repair. This enzyme is responsible for reducing all four ribonucleotide substrates, with specificity regulated by the binding of an effector to a distal allosteric site. In all characterized RNRs, the binding of effector dATP alters the active site to select for pyrimidines over purines, whereas effectors dGTP and TTP select for substrates ADP and GDP, respectively. Here, we have determined structures of Escherichia coli class Ia RNR with all four substrate/specificity effector-pairs bound (CDP/dATP, UDP/dATP, ADP/dGTP, GDP/TTP) that reveal the conformational rearrangements responsible for this remarkable allostery. These structures delineate how RNR ‘reads’ the base of each effector and communicates substrate preference to the active site by forming differential hydrogen bonds, thereby maintaining the proper balance of deoxynucleotides in the cell.
https://doi.org/10.7554/eLife.07141.001eLife digest
DNA contains the instructions required to make proteins and other molecules in cells. DNA is made of four building blocks called deoxyribonucleotides, which are in turn made from molecules called ribonucleotides by enzymes known as ribonucleotide reductases (RNRs for short). RNR enzymes are responsible for maintaining a good balance in the levels of the different deoxyribonucleotides in cells, which is essential for DNA to be made and repaired correctly.
Previous work has shown that each RNR can act on all four ribonucleotides. However, these enzymes become more selective for certain ribonucleotides depending on which deoxyribonucleotide is most common within the cell. For example, when a deoxyribonucleotide called dGTP is plentiful, it binds to a so-called “specificity site” on the enzyme and alters the shape of the enzyme’s active site. This then means that a ribonucleotide called ADP will bind in preference to the other ribonucleotides. However, it was not clear how the binding of deoxyribonucleotides to the enzyme influences the shape of the active site.
Zimanyi et al. used a technique called X-ray crystallography to determine the three-dimensional structures of a bacterial RNR enzyme when it is bound to all four different combinations of deoxyribonucleotides and ribonucleotides. In the absence of nucleotides, the active site adopts a shape that resembles an open barrel. However, when RNR is bound to a deoxyribonucleotide at the specificity site and a ribonucleotide at the active site, the barrel clamps down, bringing the specificity site and the active site closer together. Additionally, a loop of the protein interacts with each of the deoxyribonucleotides in a different way and communicates their identity directly to the active site, which rearranges itself to hold on to the corresponding preferred ribonucleotide.
Zimanyi et al.’s findings provide an explanation for how RNRs can select between ribonucleotides so that they produce a good balance of deoxyribonucleotides in cells. This will inform future efforts to develop molecules that inhibit RNRs, which may have the potential to be used to treat bacterial infections or to kill cancer cells.
https://doi.org/10.7554/eLife.07141.002Introduction
Deoxyribonucleotides, the building blocks for DNA biosynthesis, are produced in the cell from ribonucleotide precursors by ribonucleotide reductase (RNR) (Figure 1). Three classes of RNRs are known, categorized by the cofactor they use to generate a protein radical required for catalysis. The best characterized of all RNRs is the class Ia enzyme from Escherichia coli that employs a di-iron-tyrosyl-radical cofactor to initiate chemistry and requires two dimeric protein subunits for enzymatic activity. The α2 subunit contains two (β/α)10 barrels, which house the active sites at the barrel centers (Eriksson et al., 1997; Uhlin and Eklund, 1994), and the β2 subunit utilizes a largely helical secondary structure to house the radical cofactor (Sjöberg and Reichard, 1977) (Figure 1B–C). As a central controller of nucleotide metabolism, RNR uses multiple allosteric mechanisms to maintain the balanced deoxyribonucleoside triphosphate (dNTP) pools that are required for accurate DNA replication. First, allosteric activity regulation modulates the overall size of dNTP pools. ATP or dATP binding at an allosteric activity site, found at the N-terminus of α2 (Figure 1D), leads to up-regulation or down-regulation of enzyme activity, respectively (Brown and Reichard, 1969). In E. coli class Ia RNR, this regulation is achieved by changes in the oligomeric arrangement of the α2 and β2 subunits (Brown and Reichard, 1969; Rofougaran et al., 2008; Ando et al., 2011). When ATP is bound at the activity site, an α2β2 complex is favored. Although no X-ray structure of the active complex has been determined, low resolution models have been generated using small-angle X-ray scattering (Ando et al., 2011), electron microscopy (Minnihan et al., 2013), and distance measurements made through spectroscopic analyses (Seyedsayamdost et al., 2007) (Figure 1D). This active α2β2 complex is capable of a long-range proton coupled electron transfer from β2 to α2, forming a transient thiyl radical on Cys439 to initiate catalysis (Licht et al., 1996). Alternatively, when concentrations of dATP become too high in the cell, dATP binds at the allosteric activity site and formation of an α4β4 complex is promoted. The structure of this complex was recently solved (Ando et al., 2011), revealing a ring of alternating α2 and β2 units that cannot form a productive electron transfer path, thus inhibiting the enzyme (Figure 1D).
The second form of allosteric regulation is specificity regulation, which maintains the proper relative ratios of dNTPs in the cell. Briefly, the binding of (d)NTP effectors to an allosteric specificity site in α2 influences the preference of RNR for its four nucleoside diphosphate (NDP) substrates. Whereas high levels of dATP inhibit class Ia RNR, at lower levels, dATP promotes CDP or UDP reduction. Likewise, TTP promotes GDP reduction, and dGTP promotes ADP reduction (Figure 1E) (Brown and Reichard, 1969; Rofougaran et al., 2008; von Döbeln and Reichard, 1976). Importantly, the affinity of the α2 and β2 subunits for each other is weak (~0.4 μM) in the absence of effectors, whereas the binding of a complementary substrate/specificity effector pair increases the affinity of the class Ia RNR subunits fivefold (Crona et al., 2010; Hassan et al., 2008). Previous structural work, which includes: X-ray structures of GDP and TTP bound to E. coli α2 (Eriksson et al., 1997), structures of all four substrate/effector pairs bound to class Ia α2 from Saccharomyces cerevisiae (Xu et al., 2006), and class II α2 from Thermotoga maritima (Larsson et al., 2004), revealed the location of the allosteric specificity sites at the ends of a four helix bundle at the dimer interface (Figure 1B). These data, and accompanying in vitro and in vivo studies on S. cerevisiae (Ahmad et al., 2012; Kumar et al., 2010; Kumar et al., 2011), also implicated which residues (Gln294 and Arg298, E. coli numbering) and which regions of the structure are involved in the communication between the specificity site and the active site. A flexible loop, termed loop 2 (residues 292–301 in E. coli) (Figure 1B), bridges the two sites, which are approximately 15 Å apart, and becomes more ordered upon effector and substrate binding (Eriksson et al., 1997). An additional loop, termed loop 1 (residues 259–278), is also near the effector-binding site and is stabilized upon effector binding (Eriksson et al., 1997). Thus, previous work established the players involved in allosteric specificity regulation in RNR. This work supports those assignments, and goes on to provide a novel set of crystallographic snapshots that reveal how these residues in the prototypic RNR from E. coli are able to communicate and thereby regulate substrate preference.
Results
We have utilized an α4β4 crystal form of the E. coli class Ia RNR (Figure 1D, right) (Ando et al., 2011; Zimanyi et al., 2012) to obtain structures of all four substrate/specificity effector pairs (CDP/dATP, UDP/dATP, ADP/dGTP or GDP/TTP) bound to the α2 subunit (Table 2). Each model consists of two α2 dimers, two β2 dimers, and activity effector dATP in all four allosteric activity sites. The residues modeled for each chain of each structure are given in Table 3. Although the α4β4 form of the enzyme cannot turn over substrates due to abrogation of the radical transfer pathway, the substrate and effector binding sites are not affected (Video 1). In all structures, substrates and effectors were found ordered at their respective binding sites (loop 2 electron density shown in Figure 2).
Substrate/effector binding induces a clamping of the (β/α)10 barrel in the α4β4 crystal form
Each substrate/effector bound α4β4 structure reveals a closing of the (β/α)10 barrel around the bound substrate that is not observed in the absence of substrate. An initial superposition of one α chain from the previously reported free α2 soaked with AMP-PNP (PDB ID: 3R1R (Eriksson et al., 1997)) onto the α4β4 crystal form was performed using only residues 432–446, that comprise the conserved active site loop (‘finger loop’) found in the center of the (β/α)10 barrel. In this structural alignment, one half of the active site barrel overlays exactly, whereas the second half of the barrel and the N-terminal portion of the structure has undergone a clamping movement (Figure 3). To better characterize this movement, a difference distance matrix (DDM) was calculated. This plot reveals relative movements in a reference-independent manner. Briefly, the distance between every possible pair of atoms in one chain of α from the α4β4 structure is subtracted from the same distance measured in the free α2 structure, giving a representation of how much the distances between pairs of atoms change between the two structures. The pattern of movements seen in the plot indicates that the N-terminal 220 residues (region 1 of Figure 3B,C) shift substantially (>3Å), bringing them closer to half of the active site barrel (residues ~224–439, region 2 of Figure 3B,C). We also observe that loop 2 (region 3 of Figure 3B,C) and a flexible β hairpin that sits adjacent to the N-terminus (residues 646–651, region 5 of Figure 3B,C) undergo substantial motions.
Closer inspection of the active site in the superposition (again, aligned by the active site finger loop) reveals that the movements described above directly affect substrate binding (Figure 4A). Residues Ser622, Ser625, and Thr209 form hydrogen bonds with the NDP substrate, whereas in the substrate-free state (pink ribbon in Figure 4A), these residues are shifted away from the center of the barrel. Thr209 is part of the N-terminal region that is shown to move in the DDM analysis. The structure superposition and DDM analysis also indicate that loop 2 shifts, bringing Arg298 closer to the active site. The Arg298 side chain now forms a cation-π interaction with the substrate base and a charge-charge interaction with the substrate β-phosphate. In the substrate-free structure, this side chain is pointing in the opposite direction. Protein contacts like Arg298 to the substrate diphosphate are expected to be particularly important for charge neutralization in enzymes, such as E. coli RNR, that do not employ metal cations such as Mg2+ for this purpose. The positions of residues that contact the ribose O3' (Glu441 and Asn437) are virtually unchanged between the NDP-free and NDP-bound structures (Figure 4A). All four of our substrate-bound structures show almost identical contacts between the enzyme and the substrate ribose and phosphates as described here and shown in Figure 4A,C and Figure 5 for GDP and CDP, respectively.
In contrast to the compact α2 barrel observed in our substrate/effector-bound α4β4 structures, no clamping motion of the barrel is observed between substrate-free α2 and a previously reported structure of α2 with GDP substrate and TTP effector bound (Eriksson et al., 1997). In this structure, Ser622, Ser625, and Thr209 are all 3.0 Å or greater from the substrate diphosphate, and Arg298 points away from the active site (Figure 4B). These weaker (or missing) interactions between protein and substrate in this previously determined structure are consistent with the fact that substrate is observed in an altered orientation in the active site (Figure 4C,D) and is only present at half occupancy. We attribute the differences in substrate positioning in this structure to the inability of the active site barrel to clamp around substrate due to crystal lattice contacts. Symmetry-related molecules are closer to the active site in the α2-only crystalline state than they are in the α4β4 crystal form. Importantly, differences in flexibility cannot be attributed to the fact that α4β4 is an inactive state of E. coli RNR because α2 is also an inactive state; the former because β2 is held at arm’s length from α2, and the latter because β2 is absent altogether. Additionally, other than contributing to an amenable crystal lattice, there is nothing to suggest that β2 must be present for α2 to clamp. When lattice contacts do not restrain α2 movement, α2 should be able to clamp in the presence of a cognate substrate/effector pair regardless of β2. Thus, it seems that different crystal forms have given us two distinct snapshots of substrate-bound states for the E. coli class Ia enzyme: a high affinity state which is closed and ready for radical-based chemistry when β2 becomes available (Figure 4C) and a lower affinity state in which substrate is still exposed to solvent and not ready to undergo catalysis, regardless of the availability of β2 (Figure 4D).
Specificity effectors are anchored through conserved ribose and phosphate interactions
Regardless of the identity of the dNTP base, the ribose and phosphates of the three specificity effectors (dATP, dGTP, and TTP) form almost identical contacts with RNR (Figure 6). Despite the presence of the specificity effector site at the dimer interface, these common interactions are made with only one chain of the α2 dimer. Asp232 and His275 are within hydrogen-bonding distance of O3' of the ribose. One Mg2+ ion appears to form an octahedral coordination complex with three phosphate oxygens of the specificity effector and three water molecules (Figure 6). There are also hydrogen bonds between the phosphates and the backbone amides of residues Arg269 and Leu234 and charge-charge interactions with the side chains of Arg269 and/or Arg262.
Three different loop 2 conformations are responsible for substrate specificity
Loop 2 is ordered in all four structures as indicated by 2Fo–Fc composite omit electron density (Figure 2). As noted above, Arg298 of loop 2 forms a conserved interaction with the β-phosphate of all four substrates (Figures 7, 8). On the effector side of loop 2, Cys292 packs against the base of all three effectors (Figures 7, 8), consistent with binding studies that show the affinity for all specificity effectors is decreased when Cys292 is mutated to alanine (Ormö and Sjöberg, 1996). The remaining interactions with loop 2 are distinct for each substrate/effector pair (Figures 7, 8). dATP, which preferentially increases RNR activity for CDP and UDP substrates, makes two hydrogen bonds with loop 2. The backbone amide and carbonyl of Ser293 are in position to hydrogen bond to N1 and N6, respectively, of the adenine base of dATP. These hydrogen bonds hold Ser293 in place and position the adjacent Gln294 in to a conformation that points into the substrate-binding site (Figure 8A,B). In this conformation, Gln294 can form a hydrogen bond with O2 of CDP or UDP. There is no difference in substrate positioning between CDP and UDP, as the unique positions of the pyrimidine ring (N2 and N4/O4) are not involved in any contacts with the protein. Furthermore, the positioning of Gln294 by dATP selects for pyrimidine binding, as binding of a purine base would be disfavored due to steric occlusion by Gln294.
dGTP, which preferentially increases RNR activity for substrate ADP, forms one hydrogen bond with loop 1 residue Thr276 (Figure 6C) in addition to its loop 2 contacts. It is the only effector for which a hydrogen bond is formed between the base and loop 1. Because the position of hydrogen bond donors and acceptors at N1 and O6 in dGTP is opposite from that of N1 and N6 of dATP, the interaction with the backbone amide of Ser293 in loop 2 is unfavorable in this structure. Instead, Ser293 is shifted away from the base, and N1 and N2 – both hydrogen bond donors – form an interaction with the backbone carbonyl of Gln294, which has shifted away from the substrate (Figure 7A–C and 8A–C). Ser293 thus makes no hydrogen bond contacts to dGTP and instead provides hydrogen bonds to Gly300 on the other side of loop 2 through both its side chain and backbone. Stabilization of Gly300 in this position in turn allows for a hydrogen bond between the carbonyl of Gly299 with the N6 position of the preferred substrate, ADP. The movement of Gln294 away from the substrate-binding site allows room for the larger ADP to bind. Thus, Ser293 and Gln294 are responsible both for differentiating between the purine bases adenine and guanine at the effector-binding site and for differentiating pyrimidines from purines at the substrate site, either directly (in the case of Gln294 hydrogen bonding to CDP or UDP) or indirectly through stabilization of Gly299 and Gly300. With dATP bound, Ser293 plays a role in effector recognition and Gln294 in substrate recognition; these roles are reversed for ADP/dGTP.
TTP, which preferentially increases RNR activity for substrate GDP, binds to the specificity site such that only the carbonyl of Cys292 from loop 2 is positioned to form a hydrogen-bonding interaction with the thymine base (Figure 8D). As in the ADP/dGTP structure, Gln294 is swung away from the active site to make room for the binding of a purine substrate (Figure 8D). On the substrate side of loop 2, N1 and N2 of GDP are both within hydrogen-bonding distance of the backbone carbonyl of Gly299. N2 of GDP is also within hydrogen-bonding distance of the backbone carbonyl of Ala252 (not shown), giving GDP one more hydrogen bond than the other substrates
For the purine substrates ADP and GDP, preference appears to rely on a shift of the carbonyl of Gly299 in loop 2 forward or backward in the active site, enabled by the binding of dGTP or TTP, respectively, at the specificity site. With dGTP bound, Ser293 is positioned such that it pushes the carbonyl of Gly299 'forward' (away from the effector-side of loop 2, see Figure 8C), allowing a hydrogen bond to form between the carbonyl of Gly299 and N6 of adenine. With the smaller TTP in the effector site, Ser293 rearranges such that the carbonyl of Gly299 can relax ‘backward’ (toward the effector-side of loop 2, see Figure 8D), allowing for a hydrogen bond between the carbonyl of Gly299 and N1 and N2 of guanine. Specificity is imparted by the fact that the ‘forward’ position of the Gly299 carbonyl interacts with the 6-position of the purine, which is a hydrogen bond donor (NH2) in adenine and a hydrogen bond acceptor (O) in guanine; and the ‘backward’ position of the Gly299 carbonyl interacts with the 1-position and 2-position of the purine, which are both hydrogen bond donors (NH and NH2) in guanine but not in adenine (Figure 8C,D). To briefly summarize purine specificity, purines are favored over pyrimidines when Gln294 is positioned away from the active site creating a larger substrate-binding pocket, and ADP is favored over GDP when dGTP binding stabilizes Ser293 close to the Gly299 backbone pushing its carbonyl forward, and GDP is favored over ADP when TTP allows Ser293 to fall back.
Activity of mutant proteins Gln294Ala-α2 and Arg298Ala-α2 support structural observations
Our sets of structures implicate Arg298 and Gln294 as key residues in the specificity regulation of E. coli RNR. To confirm this prediction, we prepared Gln294Ala and Arg298Ala mutant enzymes and tested the activity for each substrate/effector pair. Because dATP at high concentrations inhibits RNR activity by binding to the allosteric activity site, whereas at lower concentrations it is a specificity effector, we first carried out control experiments to determine activity levels for CDP reduction in the absence of dATP and at an inhibitory dATP concentration (175 μM) (Figure 9). Results show that wild-type RNR is as active (CDP/dATP) or more active (UDP/dATP) at 1 μM dATP than it is in the positive control (CDP/ATP), indicating that at 1 μM, dATP is acting as a specificity effector and not as an allosteric inhibitor of activity, consistent with previous studies (Birgander et al., 2004). The negative control shows the substantial decrease in wild-type RNR activity on CDP when the concentration of dATP is high enough (175 μM) for dATP to bind to the allosteric activity site and inhibit the enzyme. With these controls in place, the activity of mutant RNRs with all four substrate/effector pairs can be established. Our structures predict that mutation of Arg298 to alanine should decrease activity for all substrates, and that is exactly what we observe (Figure 9). The effect is dramatic for this mutant, with negligible activity observed in all cases. In fact, the small amount of activity detected could be due to low levels of contamination of wild-type E. coli RNR instead of the Arg298Ala mutant protein. Thus, Arg298Ala-RNR is either mostly or completely inactive. In contrast, for Gln294Ala-RNR, we would expect purine and pyrimidine substrates to be differentially affected, with minimal or no loss of activity on ADP/GDP and measurable loss on CDP/UDP, and that is again what we observe. For Gln294Ala, activity is the same as wild-type for ADP/dGTP and perhaps even higher than wild-type for GDP/TTP, whereas for both CDP and UDP with 1 μM dATP, activity is decreased (Figure 9). These results are consistent with the removal of the Gln294 side chain from the active site to create room for the larger substrates (ADP/GDP) to bind and also with the positioning of Gln294 into the active site to stabilize CDP/UDP binding through hydrogen-bonding to the respective bases.
Discussion
In this work, we have sought to decipher the molecular rules of substrate/effector recognition in RNR and determine the basis for increased substrate affinity in the presence of a cognate effector. We observe that in the absence of substrates or effectors (for example, in the previously reported structure of α2 (PDB ID: 3R1R (Eriksson et al., 1997)), the barrel is not clamped and loop 2 is not found in any of the conformations seen in our effector-bound structures (Figure 7E). However, in the presence of substrate/effector pairs, we find a clamped down barrel and a stabilized loop 2 that adopts three different conformations depending on the allosteric effector that is bound (Figure 7A–D) (Video 2). Loop 2 residues Cys292, Ser293, and Gln294 appear to be involved in specific interactions that read out the identity of the base and communicate that identity to the active site. Briefly, dATP, which enhances reduction of CDP/UDP (von Döbeln and Reichard, 1976), makes specific contacts with Ser293, orienting Gln294 toward the active site and stabilizing the binding of both CDP and UDP (Figure 8A–B). The lack of discrimination between CDP and UDP substrates by RNR is not problematic for the cell as cytidine deaminase provides another level of control for dCTP/TTP ratios (Wang and Weiss, 1992; O'Donovan et al., 1971). In contrast, dGTP binding stabilizes Gln294 away from the RNR active site, so that CDP/UDP binding is not stabilized (although it is not prohibited) and there is room for the larger purine substrates to bind and to hydrogen bond to the carbonyl of Gly299 (Figure 8C–D). The creation of a more expansive active site is consistent with the ability of RNR to reduce both purine substrates in the presence of dGTPγS, albeit with a much greater fold activity increase for preferred substrate ADP (von Döbeln and Reichard, 1976). Specificity for ADP versus GDP appears to be modulated by whether effector-loop contacts stabilize the carbonyl of Gly299 in a forward or a backward position, respectively (Figure 8C–D). Thus, our structural data suggest that movement of Gln294 in and out of the active site alternatively shrinks and expands the active site for CDP/UDP versus ADP/GDP, whereas a more modest shift of the Gly299 carbonyl is involved in ADP/GDP selectivity. Mutagenesis of Gln294 to Ala is consistent with this proposal, showing decreased activity on CDP/UDP and no change or increased activity on ADP/GDP (Figure 9).
When a cognate substrate/effector pair is bound and the barrel is clamped, Arg298 is able to reach across the active site, stack against the NDP base, and hydrogen bond to the NDP β-phosphate, thereby sequestering the substrate in the active site, and neutralizing the substrate negative charge. These interactions yield what appears to be a high affinity substrate-bound state that is sequestered by solvent and thus amenable to radical-based chemistry. To invoke an analogy, Arg298 is like a latch on a suitcase, locking the active site barrel in cases in which the active site is appropriately packed (Figure 10A). When a suitcase is packed with too many clothes, the latch cannot reach the lock, and the clothes are not secured. Similarly, a mismatched substrate/effector pair such as a purine nucleotide with dATP would not be expected to allow the barrel to clamp and loop 2 to rearrange such that Arg298 can reach the substrate phosphate and thus form a ‘latched’ complex (Figure 10B). In the latter case, one would expect the mismatched substrate to be released from the enzyme, so that a complementary one can bind. Chemical logic tells us that allosteric regulation of specificity would require both low and high affinity substrate/effector-enzyme states; low affinity states to sample substrate/effector pairs, and high affinity states to capture the correct pair, recruit β2 (class I) or adenosylcobalamin (class II), and initiate catalysis (Figure 10). Although the molecular basis for the fivefold increase in binding affinity of β2 for α2 in E. coli class Ia RNR in the presence of bound substrate/effector pairs is not understood, the clamping of the barrel that we observe here may be part of the molecular explanation. The use of Arg298 as a molecular latch to secure the barrel in this closed state is beautiful in its simplicity. It is also consistent with both our structural and our biochemical data that show almost complete loss of activity in the Arg298Ala variant (Figure 9).
Though it is generally accepted that the rules of specificity regulation are conserved among class I, II, and III RNRs, this hypothesis is based on a relatively small number of characterized RNRs. Loop 2 is clearly an important player in specificity regulation, and yet it is not highly conserved, even among characterized RNRs (Figure 11). This being said, our structures show that backbone atoms of loop 2 make the vast majority of contacts to the substrate base and the effector base, limiting the number of residues that need to be strictly conserved to allow for the same molecular mechanism. Also, despite the fairly low sequence identities between the E. coli and the T. maritima α subunits (20.5%) and the E. coli and the S. cerevisiae α subunits (29.3%), a structural comparison shows that all three enzymes interact with the phosphate groups and ribose moieties of their substrates and effectors in similar ways (Table 4).
The only noteworthy difference is that Arg298, critical for reduction of all four NDPs in E. coli, directly contacts the β-phosphate of all four substrates. In contrast, none of the four S. cerevisiae RNR structures show a direct contact to the β-phosphate by the equivalent arginine, although one structure shows a through-water contact, and only one of the substrate/effector bound structures of T. maritima RNR (GDP/TTP) shows direct contact. This GDP/TTP-bound structure is the only T. maritima RNR structure in which loop 2 is fully ordered and it is remarkably similar to our E. coli structures (Figure 10H,I). All contacts made to GDP are the same: Arg side chain to β-phosphate, and loop 2 carbonyl to N1 and N2 of base (Figure 10H, 4C). The TTP contacts are not identical, but in both cases, the resulting loop 2 conformations allow the carbonyl of 299 (E. coli numbering) and 208 (T. maritima numbering) to ‘fall back’ to contact N1 and N2 of GDP.
Excitingly, the T. maritima class II RNR structures also reveal a movement of Gln294 in and out of the active site in response to effector binding. In particular, the GDP/TTP and CDP/dATP bound structures from T. maritima class II RNR, which have an ordered and semi-ordered loop 2, respectively, show that Gln203 (equivalent to Gln294 in E. coli) interacts directly with CDP when dATP is bound, and moves away from the active site when GDP/TTP are bound (Larsson et al., 2004). Unfortunately, many loop 2 residues are disordered and thus not modeled in the T. maritima structures with CDP/dATP, UDP/dATP, and ADP/dGTP, preventing further comparison to our E. coli structures.
In stark contrast to the behavior of Gln294 in E. coli and Gln203 in T. maritima, Gln288 in S. cerevisiae class Ia RNR is observed to point into the active site in all structures, regardless of which substrate/effector pairs are bound. Structural comparisons show that the protein does not pack as tightly around purine substrates when Gln points into the active site as it does when Gln is flipped away (Figure 10E,F compared with 10H,I). It is possible that this difference in Gln positioning in S. cerevisiae, and differences described above for Arg298 (E. coli numbering), indicate that the mechanism of allosteric regulation is not conserved among RNRs. However, it is also possible that these observed structural differences are a result of crystal packing variations or other deviations in how crystals were prepared. Additional structures of high-affinity RNR-substrate complexes will help to clarify the degree to which the roles of Gln294 and Arg298 are conserved across RNR species. Regardless of their exact roles in S. cerevisiae RNR, there are a number of studies that support the importance of these residues (Ahmad et al., 2012; Kumar et al., 2010; Kumar et al., 2011). In particular, when mutant RNR is the only RNR being expressed in S. cerevisiae, mutation of Arg293 (the Arg298 equivalent) to Ala is lethal, and mutation of Gln288 (the Gln294 equivalent) to Ala yields a severe S phase defect (Ahmad et al., 2012). Gln288Ala mutation in S. cerevisiae also leads to substantially elevated dGTP/dATP levels compared with dCTP/TTP levels when one compares wild-type S. cerevisiae expressing two wild-type RNRs with mutant S. cerevisiae expressing one wild-type RNR and one Gln288Ala mutant RNR (Kumar et al., 2010). Under the same experimental conditions, all four dNTPs are elevated by similar amounts (within 1–3%) for the Arg293Ala mutation in S. cerevisiae (Kumar et al., 2010). Although more studies are clearly needed to confirm or refute that the molecular basis of allosteric specificity regulation among RNRs is conserved, it is interesting to note that our in vitro data on E. coli (Figure 9) is consistent with the more severe phenotype in S. cerevisiae for Arg than Gln mutation and is consistent with the observation that mutation of Gln in S. cerevisiae elevates purines over pyrimidines whereas mutation of Arg does not show differential elevation.
In conclusion, the work presented here provides a unifying mechanism for substrate specificity regulation in the most studied RNR, the E. coli class Ia enzyme. Our structures show how each specificity effector is read out at a distal allosteric site and how that information is communicated to the active site where residues rearrange such that specific hydrogen bonds can be formed with the cognate substrate base. When an effector/substrate match is discovered, the barrel is clamped and latched in preparation for catalysis. Just as DNA replication and transcription take advantage of the unique hydrogen-bonding properties of each nucleotide base, enzymatic ribonucleotide reduction also employs these unique hydrogen-bonding properties for specificity regulation. Through an elegant set of protein rearrangements, E. coli RNR screens and selects its substrate from the four potential NDPs, ensuring appropriate pools of deoxynucleotides are available for DNA biosynthesis and repair.
Materials and methods
Effector and substrate preparation
Request a detailed protocolFor the enzyme assays, sodium salts of CDP, ADP, UDP, GDP, dATP, dGTP, and TTP were purchased from Sigma-Aldrich (St. Louis, MO) and dissolved into assay buffer (50 mM HEPES pH 7.6, 15 mM MgCl2, 1 mM ethylenediaminetetraacetic acid (EDTA)). The pH of each solution was slowly adjusted to 7–8 with NaOH, and the nucleotide concentrations were determined spectroscopically, using ε271 of 9.1 mM-1 cm-1 for CDP, ε259 of 15.4 mM-1 cm-1 for ADP, ε262 of 10.0 mM-1 cm-1 for UDP, ε253 of 13.7 mM-1 cm-1 for GDP, ε259 of 15.4 mM-1 cm-1 for ATP, ε259 of 15.2 mM-1 cm-1 for dATP, ε253 of 13.7 mM-1 cm-1 for dGTP, and ε262 of 9.6 mM-1 cm-1 for TTP. Preparation of nucleotides for crystallography was the same as described above except that 100 mM solutions of sodium salts of dATP, dGTP, and TTP were purchased from USB Corporation (Cleveland, OH) or Invitrogen (Carlsbad, CA).
Protein preparation for crystallography
Request a detailed protocolThe α2 and β2 proteins were prepared as described (Salowe et al., 1987; Salowe and Stubbe, 1986). The concentrations of α2 and β2 were determined using ε280 of 189 and 131 mM-1cm-1, respectively; unless noted otherwise, all molar concentrations are dimer concentrations (i.e. α2 or β2). For all structures, hydroxyurea-inactivated β2 (met-β2) was used in place of active β2. Met-β2 was prepared from purified active β2 as described (Ando et al., 2011). α2 had a specific activity of 3800 nmol min-1 mg-1. Prior to hydroxyurea treatment, β2 had a specific activity of 7700 nmol min-1 mg-1 as determined by a coupled spectrophotometric assay (Ge et al., 2003).
Crystallization
Request a detailed protocolCrystals were grown using the hanging drop vapor diffusion technique by mixing 1 μL of protein (25 μM α2 and 50 μM met-β2 in 50 mM HEPES, 15 mM MgCl2, and 1 mM EDTA, pH 7.6, supplemented with 10 mM dATP and 5 mM DTT, 1% (w/v) isopropyl-β-thiogalactopyranoside) with 1 μL of precipitant solution and equilibrating over a reservoir of 500 μL of precipitant at room temperature (~25°C). For the CDP/dATP structure, the precipitant solution was 9.5% (w/v) PEG 3350, 100 mM MOPS pH 7.5, 250 mM Mg(CH3COO)2, 25 mM MgCl2, and 5% (v/v) glycerol. For the UDP/dATP, ADP/dGTP, and GDP/TTP structures, the precipitant solution was 12% (w/v) PEG 3350, 100 mM MOPS pH 7.5, 300 mM Mg(CH3COO)2, 30 mM MgCl2, and 5% (v/v) glycerol. Freshly prepared drops were streak seeded with microcrystals grown under the same conditions with the exception of the MgCl2 concentration in the precipitant being 100 mM MgCl2. Streak seeding was used to obtain larger single crystals.
After 2 days of growth, crystals were transferred to a drop of soaking solution. For the CDP/dATP structure, dATP was co-crystallized, and CDP was soaked into the structure using a solution containing 10.5% (w/v) PEG 3350, 100 mM MOPS pH 7.5, 25 mM MgCl2, 250 mM Mg(CH3COO)2, 5% (v/v) glycerol, 5 mM DTT, and 10 mM CDP. For the GDP/TTP, UDP/dATP and ADP/dGTP structures, the soaking solutions were the same as for CDP except for having 13% (w/v) PEG 3350, and each nucleotide at 10 mM. Crystals were left in the soaking solution for 2 min and then cryoprotected. CDP-soaked crystals were looped through a solution of 12% (w/v) PEG 3350, 100 mM MOPS pH 7.5, 250 mM Mg(CH3COO)2, 60 mM MgCl2 and 10%, 15%, and 20% (v/v) glycerol in succession and then plunged directly into liquid N2. GDP/TTP, UDP/dATP, and ADP/dGTP soaked crystals were looped through a solution of 14% (w/v) PEG 3350, 100 mM MOPS pH 7.5, 300 mM Mg(CH3COO)2, 100 mM MgCl2, and 10%, 15%, and 20% (v/v) glycerol in succession and then plunged directly into liquid N2.
Structure solution and refinement
Request a detailed protocolDiffraction data were collected at the Advanced Photon Source at Argonne National Laboratory. The CDP/dATP data set was collected on beamline 24ID-C at 100 K on an ADSC Q315 detector. The UDP/dATP, ADP/dGTP, and GDP/TTP data sets were collected on beamline 24ID-C at 100 K on a Pilatus 6M detector. All data were processed using HKL2000 (Otwinowski and Minor, 1997) (Table 2).
All four α4β4 complex structures soaked with substrates and effectors were solved to the full extent of the data resolution using the previously published 3.95-Å resolution structure (Zimanyi et al., 2012) with all nucleotides removed. Rfree test sets were chosen to contain the same reflections across all data sets. For all structures, initial refinement was carried out in CNS 1.3 (Brünger et al., 1998) and later in Phenix (Afonine et al., 2012) with model building performed in COOT (Emsley et al., 2010). Refinement consisted of rigid body, simulated annealing, positional and individual B-factor refinement with no sigma cutoff. Loose non-crystallographic symmetry restraints were used throughout refinement. Simulated annealing composite omit maps generated in CNS 1.3 and Phenix were used to verify the models. All figures were made using the PyMOL Molecular Graphics System, version 1.5.0.4 (Schrödinger, LLC). Final refinement statistics are shown in Table 2. In addition to the substrate/specificity effector pairs bound in each structure, all structures have dATP or its hydrolysis product dADP in the allosteric activity site. In some cases, the high concentration of nucleotide used (10 mM) has resulted in more than one nucleotide molecule bound near this allosteric site.
Calculation of difference distance matrix plots
Request a detailed protocolDifference distance matrix plots were produced using the DDMP program from the Center for Structural Biology at Yale University (New Haven, CT). Residues 4 to 737 of chain A from each structure were used for the analysis.
Mutagenesis and protein purification for activity assays
Request a detailed protocolGln294Ala-α2 and Arg298Ala-α2 were constructed from the previously published vector pET-nrdA (Minnihan et al., 2011) through QuikChange site-directed mutagenesis (Agilent, Santa Clara, CA) using the following primers from Integrated DNA Technologies (Coralville, IA): Gln294Ala forward primer (5'’-AAATCCTGCTCTGCGGGCGGTGTGC-3’), Gln294Ala reverse primer (5’-GCACACCGCCCGCAGAGCAGGATTT-3’), Arg298Ala forward primer (5’-CGTTGCCGCACCGCCAGCCACACCGC-3’), and Arg298Ala reverse primer (5’-GCGGTGTGGCTGGCGGTGCGGCAACG-3’). All mutations were confirmed by DNA sequencing performed by Genewiz (South Plainfield, NJ).
Expression and purification of N-terminal hexahistidine tagged E. coli wild-type α2 and N-terminal hexahistidine tagged mutant α2 variant proteins (Gln294Ala-α2 and Arg298Ala-α2) were carried out as previously described (Minnihan et al., 2011). Briefly, cells were resuspended in 40 mL of Buffer A (50 mM Tris pH 7.6, 300 mM NaCl, 1 mM TCEP), lysed by sonication, and clarified by centrifugation at 29,000 × g. Lysate was applied to a 5 mL HisTrap HP column (GE Healthcare Life Sciences, Pittsburg, PA), washed with Buffer A supplemented with 30 mM imidazole, and eluted with Buffer A supplemented with 300 mM imidazole. Protein was further purified on a Superdex 200 16/60 size exclusion column (GE Healthcare Life Sciences) and transferred to a final storage buffer of 20 mM HEPES 7.6, 100 mM NaCl, and 5% glycerol. A final yield of ~25–50 mg/L of culture for wild-type α2 is typical. The purification for the mutants was identical, with similar yields. All proteins were judged as purified to homogeneity by sodium dodecyl sulfate/polyacrylamide gel electrophoresis (SDS/PAGE), and their concentrations were determined using ε280 of 189 mM-1cm-1. Hexahistidine tags were not removed since previous studies showed that these tags on E. coli α2 do not significantly alter activity (Minnihan et al., 2011). Untagged E. coli β2 was purified as previously described (Salowe and Stubbe, 1986) and contained ~1.1 radicals per dimer as estimated by UV-visible spectroscopy of the Y122 radical (ε411 of 1760 mM-1 cm-1) after drop-line subtraction of the diferric cluster absorbance (Bollinger et al., 1995). Untagged E. coli β2 was exchanged into a storage buffer containing 50 mM HEPES 7.6 and 5% glycerol, and its concentration determined using ε280 of 131 mM-1 cm-1.
For the coupled assay described below, E. coli thioredoxin reductase (TrxR) and E. coli thioredoxin (Trx) were prepared. The gene for TrxR was amplified from E. coli genomic DNA using primers with ends suitable for a second round of PCR to generate the plasmid borne gene in the pRham SUMO fusion vector (Lucigen, Middleton, WI). The primers, ordered from Integrated DNA Technologies, were: forward 5’-CGCGAACAGATTGGAGGTGGCACGACCAAACACAGTAAACTG-3’, and reverse 5’-GTGGCGGCCGCTCTATTATTTTGCGTCAGCTAAACCATCGAG-3’. DNA sequencing performed by Genewiz was used to confirm the sequence of the resulting construct. This construct contains a hexahistidine tag and SUMO protein fusion at the N-terminus under a rhamnose promoter. The resulting protein was expressed according to the plasmid manufacturer protocol (Lucigen) and purified as described above for hexahistidine-tagged α2, using the same buffers. The amount of flavin cofactor was quantified by absorption at 440 nm (Thelander, 1967; Gleason et al., 1990), and a specific activity of 38,200 nmol min-1 mg-1 was determined by a coupled assay with 30 µM TR and 150 µM Ellman’s reagent (5,5'-dithiobis[2-nitrobenzoic acid]) (Pierce, Rockford, IL) (Luthman and Holmgren, 1982). The final concentration for assays (0.5 µM) was determined by flavin concentration as the as-expressed protein was only 60% loaded with flavin as determined by the A280/A440 ratio. The resulting protein behaved identically in assays as untagged TrxR purified from an overexpressing E. coli strain (Russel and Model, 1985). E. coli Trx was purified as previously described (Chivers et al., 1997).
Activity assays
Request a detailed protocolActivity assays for wild-type α2 and the Gln294Ala and Arg298Ala mutants were performed using a continuous, coupled, spectrophotometric assay monitoring the consumption of NADPH by the Trx/TrxR system (Ge et al., 2003). All experiments were performed on a Cary Bio300 spectrometer (Agilent) with data analysis performed using the Cary WinUV Kinetics program (Agilent) and Microsoft Excel. The assay buffer consisted of 50 mM HEPES pH 7.6, 15 mM MgCl2, 1 mM EDTA, and the following substrate and effector concentrations were used: 1 mM CDP and 3 mM ATP (control for presence of activity), 1 mM CDP and 175 μM dATP (control for inactivation by dATP), 1 mM ADP and 120 μM dGTP, 1 mM GDP and 250 μM TTP, and 1 mM CDP/UDP and 1 μM dATP. Substrate, effector, E. coli Trx (30 μM), E. coli TrxR (0.5 μM), and NADPH from Sigma-Aldrich (200 μM) were mixed in assay buffer, and the reaction was initiated by the addition of α2 (0.1 μM) and wild-type β2 (1 μM) to a final volume of 120 μL. The decrease in NADPH absorbance at 340 nm was monitored over 1 min. The basal level of NADPH oxidation was monitored over 30 s prior to the addition of enzyme.
Data availability
References
-
Towards automated crystallographic structure refinement with phenix.refineActa Crystallographica Section D Biological Crystallography 68:352–367.https://doi.org/10.1107/S0907444912001308
-
Structural interconversions modulate activity of Escherichia coli ribonucleotide reductaseProceedings of the National Academy of Sciences of the United States of America 108:21046–21051.https://doi.org/10.1073/pnas.1112715108
-
Mutant R1 proteins from Escherichia coli class Ia ribonucleotide reductase with altered responses to dATP inhibitionJournal of Biological Chemistry 279:14496–14501.https://doi.org/10.1074/jbc.M310142200
-
Role of effector binding in allosteric control of ribonucleoside diphosphate reductaseJournal of Molecular Biology 46:39–55.https://doi.org/10.1016/0022-2836(69)90056-4
-
Crystallography & NMR system: a new software suite for macromolecular structure determinationActa Crystallographica. Section D, Biological Crystallography 54:905–921.
-
Microscopic pka values of Escherichia coli thioredoxinBiochemistry 36:14985–14991.https://doi.org/10.1021/bi970071j
-
Subunit and small-molecule interaction of ribonucleotide reductases via surface plasmon resonance biosensor analysesProtein Engineering Design and Selection 23:633–641.https://doi.org/10.1093/protein/gzq035
-
Features and development of CootActa Crystallographica Section D Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Methodology to probe subunit interactions in ribonucleotide reductasesBiochemistry 47:13046–13055.https://doi.org/10.1021/bi8012559
-
Highly mutagenic and severely imbalanced dNTP pools can escape detection by the S-phase checkpointNucleic Acids Research 38:3975–3983.https://doi.org/10.1093/nar/gkq128
-
Mechanisms of mutagenesis in vivo due to imbalanced dNTP poolsNucleic Acids Research 39:1360–1371.https://doi.org/10.1093/nar/gkq829
-
Enzymatic synthesis of deoxyribonucleotides. IX. Allosteric effects in the reduction of pyrimidine ribonucleotides by the ribonucleoside diphosphate reductase system of Escherichia coliThe Journal of Biological Chemistry 241:2533–2539.
-
Structural mechanism of allosteric substrate specificity regulation in a ribonucleotide reductaseNature Structural & Molecular Biology 11:1142–1149.https://doi.org/10.1038/nsmb838
-
Kinetics of radical intermediate formation and deoxynucleotide production in 3-aminotyrosine-substituted Escherichia coli ribonucleotide reductasesJournal of the American Chemical Society 133:9430–9440.https://doi.org/10.1021/ja201640n
-
Generation of a stable, aminotyrosyl radical-induced a2α2bβ2 complex of Escherichia coli class Ia ribonucleotide reductaseProceedings of the National Academy of Sciences of the United States of America 110:3835–3840.https://doi.org/10.1073/pnas.1220691110
-
Deoxycytidine triphosphate deaminase: characterization of an Escherichia coli mutant deficient in the enzymeJournal of Bacteriology 105:666–672.
-
Processing of x-ray diffraction data collected in oscillation modeMacromolecular Crystallography, Pt A 276:307–326.https://doi.org/10.1016/S0076-6879(97)76066-X
-
Oligomerization status directs overall activity regulation of the Escherichia coli class Ia ribonucleotide reductaseJournal of Biological Chemistry 283:35310–35318.https://doi.org/10.1074/jbc.M806738200
-
Direct cloning of the TrxB gene that encodes thioredoxin reductaseJournal of Bacteriology 163:238–242.
-
Cloning, overproduction, and purification of the Bβ2 subunit of ribonucleoside-diphosphate reductaseJournal of Bacteriology 165:363–366.
-
PELDOR spectroscopy with DOPA-β2 and NH2Y-α2s: distance measurements between residues involved in the radical propagation pathway of E. coli ribonucleotide reductaseJournal of the American Chemical Society 129:15748–15749.https://doi.org/10.1021/ja076459b
-
Nature of the free radical in ribonucleotide reductase from Escherichia coliThe Journal of Biological Chemistry 252:536–541.
-
Thioredoxin reductase. characterization of a homogenous preparation from Escherichia coli BThe Journal of Biological Chemistry 242:852–859.
-
Binding of substrates to Escherichia coli ribonucleotide reductaseThe Journal of Biological Chemistry 251:3616–3622.
-
Dcd (dCTP deaminase) gene of Escherichia coli: mapping, cloning, sequencing, and identification as a locus of suppressors of lethal dut (dUTPase) mutationsJournal of Bacteriology 174:5647–5653.
-
Structures of eukaryotic ribonucleotide reductase I provide insights into dNTP regulationProceedings of the National Academy of Sciences of the United States of America 103:4022–4027.https://doi.org/10.1073/pnas.0600443103
Article and author information
Author details
Funding
National Institutes of Health (T32GM08334)
- Christina M Zimanyi
Howard Hughes Medical Institute
- Catherine L Drennan
National Institutes of Health (P30-ES002109)
- Catherine L Drennan
National Science Foundation (0645960)
- Michael A Funk
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. EJ Brignole for a critical reading of this manuscript. This work is based upon research conducted at the Advanced Photon Source on the Northeastern Collaborative Access Team beamlines, which are supported by award GM103403 from the National Center for Research Resources at the National Institutes of Health. Use of the Advanced Photon Source is supported by the U.S. Department of Energy, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. This work was supported by the National Institutes of Health grants T32GM08334 (CMZ), P30-ES002109 (CLD), and a National Science Foundation Graduate Research Fellowship under Grant No. 0645960 (MAF). CLD is a Howard Hughes Medical Institute Investigator. The Protein Data Bank accession numbers for the CDP/dATP, UDP/dATP, ADP/dGTP, and GDP/TTP E. coli class Ia RNR structures are 5CNS, 5CNT, 5CNU, and 5CNV, respectively.
Copyright
© 2016, Zimanyi 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
-
- 5,784
- views
-
- 882
- downloads
-
- 65
- 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
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Pre-mRNA splicing is catalyzed in two steps: 5ʹ splice site (SS) cleavage and exon ligation. A number of proteins transiently associate with spliceosomes to specifically impact these steps (first and second step factors). We recently identified Fyv6 (FAM192A in humans) as a second step factor in Saccharomyces cerevisiae; however, we did not determine how widespread Fyv6’s impact is on the transcriptome. To answer this question, we have used RNA sequencing (RNA-seq) to analyze changes in splicing. These results show that loss of Fyv6 results in activation of non-consensus, branch point (BP) proximal 3ʹ SS transcriptome-wide. To identify the molecular basis of these observations, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of a yeast product complex spliceosome containing Fyv6 at 2.3 Å. The structure reveals that Fyv6 is the only second step factor that contacts the Prp22 ATPase and that Fyv6 binding is mutually exclusive with that of the first step factor Yju2. We then use this structure to dissect Fyv6 functional domains and interpret results of a genetic screen for fyv6Δ suppressor mutations. The combined transcriptomic, structural, and genetic studies allow us to propose a model in which Yju2/Fyv6 exchange facilitates exon ligation and Fyv6 promotes usage of consensus, BP distal 3ʹ SS.
-
- Biochemistry and Chemical Biology
- Neuroscience
TIMM50, an essential TIM23 complex subunit, is suggested to facilitate the import of ~60% of the mitochondrial proteome. In this study, we characterized a TIMM50 disease-causing mutation in human fibroblasts and noted significant decreases in TIM23 core protein levels (TIMM50, TIMM17A/B, and TIMM23). Strikingly, TIMM50 deficiency had no impact on the steady-state levels of most of its putative substrates, suggesting that even low levels of a functional TIM23 complex are sufficient to maintain the majority of TIM23 complex-dependent mitochondrial proteome. As TIMM50 mutations have been linked to severe neurological phenotypes, we aimed to characterize TIMM50 defects in manipulated mammalian neurons. TIMM50 knockdown in mouse neurons had a minor effect on the steady state level of most of the mitochondrial proteome, supporting the results observed in patient fibroblasts. Amongst the few affected TIM23 substrates, a decrease in the steady state level of components of the intricate oxidative phosphorylation and mitochondrial ribosome complexes was evident. This led to declined respiration rates in fibroblasts and neurons, reduced cellular ATP levels, and defective mitochondrial trafficking in neuronal processes, possibly contributing to the developmental defects observed in patients with TIMM50 disease. Finally, increased electrical activity was observed in TIMM50 deficient mice neuronal cells, which correlated with reduced levels of KCNJ10 and KCNA2 plasma membrane potassium channels, likely underlying the patients’ epileptic phenotype.