Structural Determination of Functional Units of the Nucleotide Binding Domain (NBD94) of the Reticulocyte Binding Protein Py235 of Plasmodium yoelii

Background Invasion of the red blood cells (RBC) by the merozoite of malaria parasites involves a large number of receptor ligand interactions. The reticulocyte binding protein homologue family (RH) plays an important role in erythrocyte recognition as well as virulence. Recently, it has been shown that members of RH in addition to receptor binding may also have a role as ATP/ADP sensor. A 94 kDa region named Nucleotide-Binding Domain 94 (NBD94) of Plasmodium yoelii YM, representative of the putative nucleotide binding region of RH, has been demonstrated to bind ATP and ADP selectively. Binding of ATP or ADP induced nucleotide-dependent structural changes in the C-terminal hinge-region of NBD94, and directly impacted on the RBC binding ability of RH. Methodology/Principal Findings In order to find the smallest structural unit, able to bind nucleotides, and its coupling module, the hinge region, three truncated domains of NBD94 have been generated, termed NBD94444–547, NBD94566–663 and NBD94674–793, respectively. Using fluorescence correlation spectroscopy NBD94444–547 has been identified to form the smallest nucleotide binding segment, sensitive for ATP and ADP, which became inhibited by 4-Chloro-7-nitrobenzofurazan. The shape of NBD94444–547 in solution was calculated from small-angle X-ray scattering data, revealing an elongated molecule, comprised of two globular domains, connected by a spiral segment of about 73.1 Å in length. The high quality of the constructs, forming the hinge-region, NBD94566–663 and NBD94674–793 enabled to determine the first crystallographic and solution structure, respectively. The crystal structure of NBD94566–663 consists of two helices with 97.8 Å and 48.6 Å in length, linked by a loop. By comparison, the low resolution structure of NBD94674–793 in solution represents a chair–like shape with three architectural segments. Conclusions These structures give the first insight into how nucleotide binding impacts on the overall structure of RH and demonstrates the potential use of this region as a novel drug target.


Introduction
Malaria continues to be one of the major public health problems for mankind; it is the transmissible disease with the greatest morbidity around the world [1]. The complex life cycle of the protozoan parasite is characterized by distinct invasive forms of the sporozoite and merozoite that invade hepatocytes and erythrocytes in the vertebrate host, respectively, and the ookinetes inside the insect vector that penetrates the mosquito midgut epithelium [1][2][3][4][5][6]. Invasion of red blood cells (RBC) by the merozoite and the subsequent cyclical replication of the parasite is the cause of malaria associated pathology. Multiple merozoite protein families are implicated in the invasion of RBCs, including the erythrocyte binding proteins (EBPs) and the reticulocyte binding protein homologues (RH), which bind to different RBC membrane receptors [1,[6][7][8][9][10][11][12][13][14][15][16]. Little is known about how the large RH transmembrane proteins mediate their function during erythrocyte invasion, but a crucial step appears to be the proteolytic cleavage during the invasion process [14,17]. Members of RH have been identified in all plasmodium species so far analyzed indicating the conserved function and importance of this protein family to the malaria parasite [1]. The P. yoelii RH protein, termed Py235 (235 kDa in mass), has been shown to be a potential virulence factor that allows the parasite to invade a wider range of host erythrocytes [18][19][20]. Py235 is also involved in clonal phenotypic variation of merozoites [21], enabling the parasite to evade immune responses and adapt to changes in the host environment during the invasion step [22]. Previous studies in P. vivax has indicated that RH may have an initial sensing role preceding and possibly enabling the subsequent interaction of the EBP member with its corresponding receptor [23]. Besides multiple receptor ligand interaction between the merozoite and the RBC, the effects of other intra and/or extra cellular molecules like ATP on invasion have been discussed [24][25][26][27]. ATP, released by erythrocytes under normal physiological conditions, has recently been shown to serve as a signalling molecule regulating vascular tone [24][25][26] and ATP receptors have been implicated in a number of signalling pathways [28,29]. Mechanical deformation of erythrocytes as encountered in capillaries forming the microvasculature lead to increased ATP release [30]. Importantly, intracellular ATP is a requirement for merozoite invasion [31][32][33][34] with erythrocytes that have been depleted of ATP being refractory to invasion. These findings suggest that it would be advantageous for merozoites to sense the intracellular ATP level of the erythrocyte. In this context a 94 kDa domain of Py235 of P. yoelii, which is highly conserved amongst the RBLs and called the nucleotide-binding domain (NBD94), has been found to selectively bind ATP and ADP. The amino acid sequence 488 EKL-KHYNFDDFVK 500 in NBD94 has been identified as nucleotidebinding region as shown by photoaffinity labeling of the nucleotide-analogue 8-N 3 -39-biotinyl-ATP ( Fig. 1 [27]). The preference of MgATP over MgADP recognition is associated with specific structural alterations in the C-terminal domain of NBD94 as depicted by spectroscopic comparison of NBD94 and its C-terminal truncated form, NBD94 1-550 , in which no nucleotide-dependent alteration could be observed [27]. The nucleotide effect in the recombinant protein is confirmed by a strong binding of Py235 to RBCs in the presence of MgATP which is significantly lower in the presence of MgADP or the absence of nucleotides [27]. Based on these traits and the absence of significant ATPase activity of NBD94, this domain was suggested to serve as an ATP/ADP sensor during the invasion process [27]. However, the recombinant NBD94 protein is prone to degradation over a period of time. Whether degradation occurs through autolysis or proteolytic cleavage, as described for a series of proteolytic cleavage events of both the EBP and RH protein families during invasion, resulting in 9 to 13 kDa fragments, is unclear [35,36].
In many cases a protein fragment obtained by cleavage of the full length protein or by expression of part of the protein can retain a functional domain. This is particularly relevant in drug designing. In order to get insight into the functional regions of NBD94 and to define the smallest segment, still able to bind ATP/ADP, three truncated constructs, called NBD94 444-547 , including the 8-N 3 -39-biotinyl-ATP binding sequence, as well as NBD94 566-663 and NBD94 674-793 , described to belong to the Cterminal hinge region and to couple the nucleotide-binding event [27], have been generated (Fig. 1) and purified to homogeneity. As demonstrated by fluorescence correlation spectroscopy (FCS), NBD94 444-547 forms the smallest ATP/ADP-binding segment, which becomes inhibited by the ATPase/-synthase inhibitor 4-Chloro-7-nitrobenzofurazan (NBD-Cl). The high quality of all proteins formed the platform to solve the first low resolution structure of NBD94 444-547 , and NBD94 674-793 in solution and the first crystallographic structure of NBD94 566-663 . The structural features of this ensemble will be discussed in the light of coupling of nucleotide-binding with structural rearrangements in the hinge region, followed by downstream effects, enabling the parasite to continue the invasion process.
The secondary structures of these proteins were determined from circular dichroism spectra, measured between 185-260 nm ( Fig. 2A-C). The CD-spectra of NBD94 444-547 , NBD94 566-663 and NBD94  show that all proteins are mainly a-helical, as reflected by its minima at 208 and 222 nm and as predicted from its amino acid sequence. The a-helical content of NBD94 444-547 , NBD94 566-663 and NBD94 674-793 is determined to be 83%, 86% and 69%, respectively, congruent to the values determined for the entire NBD94 and indicating the correct secondary structure of the recombinant truncated proteins generated, as well as reflecting a proper selection of stable constructs. The molar ellipticity values at 208 nm and at 222 nm of NBD94 444-547 , NBD94 566-663 and NBD94 674-793 are in a ratio of 0.94, 0.99 and 0.97, respectively.
Shape Determination of NBD94 444-547 and NBD94 674-793 in Solution The high purity allowed small-angle X-ray scattering (SAXS) experiments to be performed, with the aim to determine the first low resolution structures of the nucleotide-binding segment NBD94 444-547 and NBD94 674-793 in solution. SAXS patterns from solutions of both proteins were recorded as described in Materials and Methods to yield the final composite scattering curves in figure 3A and 4A, respectively, that both proteins are monodispersed in solution. The radius of gyration R g of NBD94 444-547 is 33.461 Å and the maximum dimension D max of the protein is 13462 Å (Fig. 3B). The gross structure of NBD94 444-547 was restored ab initio from the scattering pattern in figure 3A using the dummy residues modeling program GASBOR [37], which fitted well to the experimental data in the entire scattering range (a typical fit displayed in figure  3A, curve 2, has the discrepancy x = 1.28). Ten independent reconstructions yielded reproducible models and the average model and the most probable model are displayed in figure 3C. NBD94 444-547 appears as an elongated molecule with a length of 134 Å , composed of two more globular domains and a spiral-like segment of about 73.1 Å in length between both domains.
By comparison, the radius of gyration R g of NBD94 674-793 is 22.261 Å nm and the maximum dimension D max of the protein is 11862 Å (Fig. 4B), respectively, suggesting that the protein is rather elongated. The shoulders at larger intraparticle distances of about 70 Å and 100 Å indicate that the molecule consists of three distinct domains. Like for NBD94 444-547 the shape of NBD94 674-793 was determined ab initio from the scattering pattern in figure 4A using the program GASBOR [37]. The experimental data fitted well with the calculated data, reflected by an x value of 1.03 (

Nucleotide-Binding Determined by Fluorescence Correlation Spectroscopy
The proper structural folding and structure formation of NBD94 444-547 enabled us to study the ability of this protein to bind nucleotides by fluorescence correlation spectroscopy using fluorescent ATP and ADP derivatives ATP ATTO-647N and ADP ATTO-647N, respectively. As a reference, the mean count rate per Cyanine 5 (Cy5) fluorophore was determined to be 32.  figure 5A and in supplementary figure S1A, respectively. The increase of the mean diffusion time t D was due to the increase in the mass of the diffusing particle, when fluorescently labelled nucleotide bound to NBD94 444-547 , which is apparent in the displayed autocorrelation curves with increased protein concentrations from left to right. A binding constant (K d ) of 22862.3 mM for MgATP ATTO-647N and 33161.8 mM for MgADP ATTO-647N bound to NBD94 444-547 was determined (Fig. 5B, C). By contrast no nucleotide-binding could be observed in the constructs NBD94 566-663 and NBD94 674-793 .
Most recently, we observed that 4-Chloro-7-nitrobenzofurazan (NBD-Cl) is a potent inhibitor of nucleotide-binding in the entire NBD94 [27]. We tested, whether MgATP ATTO-647N-binding of the minimal nucleotide-core domain, NBD94 444-547 , becomes inhibited by NBD-Cl. The plotted autocorrelation functions in supplementary figure S1B show a change of the diffusion time due to an increase of the concentration of NBD-Cl from right to left for MgATP ATTO-647N. The calculated bound fraction for increasing inhibitor concentrations were plotted to determine the effect of nucleotide binding. The interaction of MgATP ATTO-647N to NBD94 444-547 showed an IC 50 value of 40063 mM (Fig. 5D). A similar effect of NBD-Cl-inhibition has been observed for NBD94 444-547 -ADP-ATTO-647N formation (data not shown).

Crystallographic Structure of NBD94 566-663
The second domain, NBD94 566-663 , covering the amino acids close to the nucleotide binding segment NBD94 444-547 , has been crystallized by vapor diffusion using 2-methyl-2,4 pentanediol as a precipitant in a form suitable for X-ray diffraction analysis. This success was due to controlling the rate of vapor diffusion via the introduction of an oil barrier over the reservoir of a vapordiffusion trial. The two different oils, paraffin and silicone, were applied as barriers, separately and also as mixtures of these two oils. Paraffin oil enabled the increase of crystal size (Supplementary Figure S2A and S2B) as well as improvement in resolution. The native crystal diffracted to 4.1 Å while the SeMet crystals diffracted to 3.9 Å at peak wavelength. A MAD dataset was collected which was indexed, integrated and scaled using the HKL-2000 [38] package. The data collection statistics are given in Table 1. Both the native and the SeMet derivative, crystallizes in the hexagonal space group P6 3 22 with three molecules in the asymmetrical unit and a solvent content and V M value of 35.44% and 1.9 Å 3 per Dalton, respectively [39]. The data collection statistics are given in Table 1.
Both MAD and SAD techniques for structure solution were tried. Various programs like, SHELX [40], SOLVE [41], CNS [42] and CCP4 [43] were used to identify the heavy atom sites (Se). Reasonable statistics and good initial map could be achieved with the SAD technique and using the help of SHELX program [40]. The resolution cut-off used for Se site identification is 4.4 Å and the correlation co-efficient stood at 39.51% with the Pseudo- free CC of 71.47%. Peak wavelength (0.97898 Å ) data was used for the SAD technique in which 12 out of 18 Se sites (6 SeMet in each molecule) could be identified by SHELXD [44] that were further refined by SHELXE [45]. These 12 Se-sites were used to phase the structure factor with SHELXE and the resulting electron densities were improved by solvent flattening with SOLOMON and density modification by DM from the CCP4 package and RESOLVE [46]. The helical region could be readily identified from the initial map and the model of the NBD94 566-663 was built manually using the program COOT [47]. Several cycles of manual building and fitting were carried out by COOT in combination with restrained refinement using REFMAC5 [48] of the CCP4 suite, keeping the temperature factor at overall. Each time a simulated annealing omit map was calculated and was used for model building (Supplementary Figure S3). The electron density map for the main chain atoms is excellent showing the typical sausage like features for the helical regions (Fig. 6A), whereas for the side chain atoms no visible densities could be identified. Se positions from the anomalous map were used to trace the chains. The final refined model has an R-factor of 33.9% and a R-free of 37.32% with good stereochemistry as can be inferred from the Ramachandran plot statistics given by PROCHECK [49]. The detailed summary of the refinement statistics is given in Table 1. The final electron density map for the structure shows good density for most of the backbone residues except for the residues 1-6 and 73-75 in chain A (Fig. 6A), 1-12 and 73-75 in chain B and 1-13, 31-33 and 73-75 in chain C, respectively (Fig. 6B).
The structure of the NBD94 566-663 molecule consists of two helices that are linked by a loop, which is not visible in the electron density map. The length of the long and short helix is 97.8 Å and 48.6 Å , which are formed by 65 and 30 residues, respectively (Fig. 6B). Amino acids Gly61 and Ser85 of the N-and C-terminal helix, respectively, are in a very close proximity of 9.85 Å . Due to the kinks from residues Tyr60-Lys62 and Ser85-Glu87 both helices spread apart by an angle of 40.2u, thereby forming a hingelike feature. The total number of amino acid residues is 107 that includes the N-terminal His-tag residues. The longer helix of chain C in the asymmetric unit is highly distorted when compared to the other molecules (Fig. 6B). Slight kinks could be noted in the long helix from residues Val20-Lys22 and Tyr60-Lys62 in chain A. Whereby in chain B, the residues Asp24-Lys26 show slight helical distortion but for residues Leu37-Lys47 the helix is deviated more. The shorter helices of all the three molecules do show some amount of distortion but is less when compared to the long helices, wherein residues Ser85-Glu87 in chain A, Glu80-Met82 in chain B and Tyr89-Lys91 in chain C show notable deviations.

Discussion
The mechanism by which a merozoite recognizes a suitable host cell is mediated by a cascade of receptor ligand interactions. In addition to the availability of the appropriate receptors, intracellular ATP plays an important role in determining if erythrocytes are suitable for merozoite invasion [24][25][26]. In erythrocyte binding assays (EBA) it has been demonstrated that Py235 of P. yoelii binds strongly to erythrocytes in the presence of ATP, whereby weak interactions have been found in the presence of ADP or the absence of nucleotides [27]. The ATP/ADP modulation of Py235receptor binding suggested a nucleotide-dependent rearrangement, making the binding domain of Py235 more accessible. Such a nucleotide induced change has been observed in the 94 kDa nucleotide-binding domain, NBD94, of Py235, in which ATPbinding causes alterations in the C-terminal hinge region [27]. The recombinant NBD94 444-547 is identified as the smallest segment of NBD94 still able to bind nucleotides with a preference of ATP-over the ADP analogue, important for sensing the signal for receptor binding of Py235. NBD94 444-547 includes the 488 EKLKHYNFDDFVK 500 peptide, observed to bind the nucleotide-analogue 8-N 3 -39-biotinyl-ATP [27]. Tyr493 is the only Tyr residue of this segment and thereby the candidate for covalently binding the potent inhibitor NBD-Cl, which was predicted to act as a nucleotide analogue in F 1 F O ATP synthases [50]. NBD-Cl interacts with the phenolic oxygen of Tyr331 in the catalytic b subunit of bovine F 1 F O ATP synthases and appears to act by preventing the nucleotide empty b subunit from converting to a productive nucleotide-binding state [51][52][53]. The hydrated NBD94 444-547 is a 134 Å long molecule, comprised of two globular segments, connected by a spiral region of about 73.1 Å in length with an a-helical content of 83%. The H 222 /H 208 ratio of 0.94 indicates that many of the residues in NBD94 444-547 are in close neighborhood, since non-interacting helices typically give ratios of around 0.8. As ATP/ADP-binding in NBD94 [27] and NBD94 444-547 is shown to be inhibited by NBD-Cl, proposed to interact with Tyr493, the elongated NBD94 444-547 becomes a new target for modified or novel classes of ATP agents.
Recently, we proposed that the ATP-binding event in NBD94 triggers a conformation rearrangement within the protein, going along with alterations in Py235 and thereby changing the binding strength of this reticulocyte binding-like protein to the red blood cell, as reflected by erythrocyte binding assays performed with Py235 of P. yoelii [27]. From nucleotide-dependent CDspectroscopy data of NBD94 and its C-terminal truncated form NBD94 1-550 it was also concluded that a hinge-like region in the C-terminal region of NBD94 may mediate the nucleotide-binding with further down stream events [27]. The NBD94 566-663 molecule with its long N-terminal and shorter C-terminal a helices, linked by a flexible loop, shows a hinge-like arrangement (Fig. 6B).
In the crystallographic structure, the residues Gly61 and Ser85 of the N-and C-terminal helix, respectively, are in a very close proximity of 9.8 Å and the very N-and C-terminal amino acids of NBD94 566-663 spread apart by an angle of 40.2u. This structural feature of NBD94 566-663 enables this protein to move either like a hinge or like a rattle up and down, thereby transmitting ATP-ADP modulation in NBD94 444-547 with up-and down movements in NBD94 566-663 , which have to be coupled to the a-helical and chair-like NBD94 674-793 fragment of the C-terminal region of NBD94, representing an elongated protein with three subdomains of 96 Å , 35 Å and 20 Å in length. Different orientation of the protein ensemble relative to each other have been speculated, in which the N-terminal segment NBD94 444-547 is linked to the Nterminal helix of NBD94 566-663 and the C-terminal segment NBD94 674-793 is oriented to the C-terminal helix of NBD94 566-663 ( Fig. 7A-C). In all hypothesized models, the structure of the NBD94 566-663 molecule forms the middle element, able to mediate the sensing of ATP-ADP binding in NBD94 444-547 with concerted conformational changes in NBD94 566-663 and NBD94 674-793 , whose subdomains might undergo structural rearrangements, transferred to down-stream events in NBD94, thereby facilitating the linkage of nucleotide signaling and Py235 binding to the blood cell. An alternative prediction is presented in figure 7D revealing the possibility of a direct NBD94 444-547 and NBD9 674-793 contact, in which NDB94 566-663 presents a linking element. The increased nucleotide affinity in presence of the C-terminal domain of NBD94 determined recently [27] may hint at such an assembly. In summary, the data presented demonstrate that NBD94 can be divided into a nucleotide-binding segment, NBD94 444-547 , existing in solution as an elongated molecule. The first crystallographic and solution structure of NBD94 566-663 and the C-terminal domain NBD94 674-793 of the C-terminal hinge region in NBD94 provide the structural basis towards a better understanding of a concerted interaction of the protein ensemble in NBD94, which is triggered by ATP, a requirement for merozoite invasion [31][32][33][34]. Furthermore, this work will provide the foundation for future studies to identify new compounds that directly interfere with the invasion process.

Biochemicals
Pfu DNA polymerase and restriction enzymes were purchased from New England Biolabs (Ipswich, MA, USA). Ni 2+ -Sepharose TM High Performance chromatography resin was obtained from GE Healthcare Bio-Sciences AB (Uppsala, Sweden). Chemicals for gel electrophoresis were received from Biorad (Hercules, USA). The ATP-and ADP-analogues EDA-ATP ATTO-647N and EDA-ADP ATTO-647N were received from ATTO-TEC (Siegen, Germany). All other chemicals were at least of analytical grade and purchased from Sigma-Aldrich (Deisenhofen, Germany), Amersham Bioscience (Buckinghamshire, UK), BIOMOL (Hamburg, Germany) and Merck (Darmstadt, Germany).

Cloning, Expression and Protein Purification
In order to obtain the truncated domains NBD94 444-547 , NBD94 566-663 and NBD94 674-793 of NBD94, the following primers were designed: forward primer 59-AAC ATC GAA CCA TGG TAA TTC CAT-39 and reverse primer 59-CGC TTA TTT GGA GCT CTT ACG TTT-39 for NBD94 444-547 ; forward primer 59-GAT CCA CCA TGG TAA AGG AAA-39 and reverse  TTA TAT TAC TAG-39 and  reverse primer 59-TTG ATT CGA GCT CTT ATA TTT TCG  ATT-39 for NBD94 674-793 . The genomic P. yoelii YM DNA was used as the template. In all the constructs the restriction sites NcoI and SacI were incorporated. Following digestion with NcoI and SacI, the PCR products were ligated into the pET9d1-His 3 vector [54]. The pET9d-His 3 vector, containing the respective gene, was then transformed into E. coli cells (strain BL21 (DE3)) and grown on 30 mg/ml kanamycin-containing Luria-Bertani (LB) agar-plates. To express NBD94 444-547 , NBD94 566-663 and NBD94 674-793 , liquid cultures were shaken in LB medium containing kanamycin (30 mg/ml) for about 20 h at 37uC until an optical density of 0.6-0.7 (OD 600 nm) was reached. To induce production of the recombinant proteins, the cultures were supplemented with isopropyl (thio)-b-D-galactoside (IPTG) to a final concentration of 1 mM. Following incubation for another 4 h at 37uC, the cells were harvested at 6 0006g for 20 min, 4uC. Selenomethionine containing protein was expressed in the same way except that methionine biosynthesis was inhibited by the growth conditions [55]. A 1:1000 preculture in LB medium was used to inoculate 300 ml of M9 minimal medium supplemented with 0.4% glucose, 2 mM MgSO 4 , 30 mg/ml kanamycin, vitamins and trace elements. After overnight growth, the culture was diluted 1:100 into 1 l of minimal medium. After an OD 600 of 0.6 was reached, 100 mg/l DL-selenomethionine (Sigma), 100 mg/l lysine, threonine and phenylalanine and 50 mg/l leucine, isoleucine and valine were added as solids. IPTG (1 mM) was added after 15 min, and cells were grown for 3 h and harvested as described above. Cells were lysed on ice by sonication for 361 min in buffer A (50 mM Tris/HCl, pH 7.5, 200 mM NaCl, and 2 mM PMSF) for NBD94 444-547 , buffer B (50 mM Tris/HCl, pH 9.0, 200 mM NaCl, and 2 mM PMSF) for NBD94 566-663 and buffer C (50 mM Tris/HCl, pH 7.5, 200 mM NaCl, 2 mM PMSF and 0.5 mM DTT) for NBD94 674-793 . The lysate was centrifuged at 10 0006g for 35 min. The supernatant was filtered (0.45 mm; Millipore) and passed over a Ni 2+ -NTA resin column to isolate NBD94 444-547 , NBD94 566-663 and NBD94 674-793 , respectively. The His-tagged protein was allowed to bind to the matrix for 1.5 h at 4uC and eluted with an imidazole-gradient (25-300 mM) in buffer A, B and C, respectively. Fractions containing His-NBD94 444-547 , NBD94 566-663 and NBD94 674-793 , were identified by SDS-PAGE [56], pooled and concentrated using Centricon YM-3 (3 kDa molecular mass cut off) spin concentrators (Millipore). The concentrated samples were applied on a sizeexclusion column (Superdex TM 75 10/300 GL, GE Healthcare) using their respective buffers with additional 10 mM EDTA. The purity of all protein samples were analyzed by SDS-PAGE [56] and the gels were stained with Coomassie Brilliant Blue G250. Protein concentrations were determined by the bicinchoninic acid assay (BCA; Pierce, Rockford, IL., USA).

Circular Dichroism (CD) Spectroscopy
Measurements of steady state CD spectra were carried out in the far UV-light (185-260 nm) using a CHIRASCAN spectropolarimeter (Applied Photophysics, UK). The CD spectroscopy measurements of NBD94 444-547 , NBD94 566-663 and NBD94 674-793 were performed in their corresponding buffers (see above) using a 60 ml quartz cell (Hellma) at 18uC with 1 nm step points and a protein concentration of 2 mg/ml for each recombinant protein. The spectrum for the buffer was subtracted from the spectrum of the protein. The ellipticity values were calculated by the average of triple determinations for each sample with a bandwidth of 1 nm from 185 to 260 nm (1 nm step points). The CD data were converted to mean residue ellipticity (H) in units of deg x dmol 21 6aa 21 using the software Chirascan version 1.2.1 (Applied Photophysics). Baseline corrected spectra were used as input for computer methods to obtain predictions of secondary structure. In order to analyze the CD spectrum, the following algorithms were used: Varselec [57], Selcons [58], Contin [59], K2D [60]; all methods as incorporated into the program Dicroprot [61] and NeuralNet [62].

Fluorescence Correlation Spectroscopy (FCS)
Fluorescence correlation spectroscopy was performed at 25uC on a LSM 510 Meta/ConfoCor 3 (Zeiss, Jena, Germany) using the ATP-analogue EDA-ATP ATTO-647N (ATTO-TEC, Siegen, Germany). The temperature was adjusted to 25uC in an incubation chamber (Zeiss). The 633 nm laser line of a HeNe633 laser was attenuated to 6 mW and focused into the aqueous solution by a water immersion objective (406/1,2 W Korr UL-VIS-IR, Zeiss). The proteins were dissolved in the buffer, containing 50 mM Tris, pH 7.5, and 200 mM NaCl. FCS was measured in 15 ml droplets of the diluted fluorescent derivatives of ATP, which were placed on Nunc 8 well chambered cover glass. Before usage, the cover glasses were treated according to Hunke et al.

Small Angle X-Ray Scattering Experiments
The synchrotron radiation X-ray scattering data for NBD94 444-547 and NBD94 674-793 were collected following standard procedures on the X33 SAXS camera [64,65] of the EMBL Hamburg located on a bending magnet (sector D) on the storage ring DORIS III of the Deutsches Elektronen Synchrotron (DESY). A photon counting Pilatus 1 M pixel detector (676420 mm 2 ) was used at a sample -detector distance of 2.4 m covering the range of momentum transfer 0.1,s,4.5 nm 21 (s = 4p sin(q)/l, where q is the scattering angle and l = 0.15 nm is the X-ray wavelength). The s-axis was calibrated by the scattering pattern of Silverbehenate salt (d-spacing 5.84 nm). The scattering patterns from NBD94 444-547 and NBD94 674-793 were measured at 8.0 mg/ml, respectively. Protein samples were prepared in 50 mM Tris/HCl (pH 7.5), 200 mM NaCl and 1.25 mM DTT as radical quencher and injected automatically using the sample-changing robot for solution scattering experiments at the SAXS station X33 [66]. The data were normalized to the intensity of the incident beam; the scattering of the buffer was subtracted and the difference curves were scaled for concentration. All the data processing steps were performed using the program package PRIMUS [67]. The forward scattering I(0) and the radius of gyration R g were evaluated using the Guinier approximation [68]. The shape of NBD94 444-547 and NBD94 674-793 in solution was built by the program GASBOR [37] as described in [69].

Crystallization of NBD94 566-663
NBD94 566-663 (10 mg/ml) was crystallized by vapor diffusion using 35% (v/v) 2-methyl-2,4-pentanediol as precipitant and acetate buffer of pH 4.5 (Wizard II Screen, Emerald BioSystems, Inc. USA) at 18uC. To avoid oxidation of selenomethionine during purification and crystallization, 1 mM tris-2-carboxyethyl-phosphine (TCEP) was added in the buffer. The crystal quality was improved by controlling the rate of vapor diffusion by the introduction of an oil barrier over the reservoir of a vapordiffusion trial. The two different oils, paraffin and silicone, were applied as barriers, separately and also as mixtures of these two oils.

Structure Determination
Crystals of selenomethionine (SeMet) substituted NBD94 566-663 were prepared for MAD phasing. Initially, crystals were screened at beamline BL12B2 at SPring-8 using a Q4R detector. Native dataset and a complete MAD dataset from single crystals of SeMet incorporated NBD94 566-663 were collected at 100 K on beamline BL26B2 at SPring-8 using a Mar CCD 225 detector. All the data were indexed, integrated and scaled using the program HKL2000 [38]. Selenium sites were identified and refined by the program SHELX [40], and density modification of the experimental map was performed with RESOLVE [46]. Manual model building and refinement of the structure was carried out in iterative cycles using COOT [47] and REFMAC5 [48]. The atomic coordinates and structure factors for the NBD94 566-663 structure has been deposited in the Protein Data Bank under accession code 3HGF.