# Structural Insights into the Nucleotide-Binding Domains of the P1B-type ATPases HMA6 and HMA8 from Arabidopsis thaliana

Hubert Mayerhofer, Emeline Sautron, Norbert Rolland, Patrice Catty, Daphné Seigneurin-Berny, Eva Pebay-Peyroula, Stéphanie Ravaud

► To cite this version:

Hubert Mayerhofer, Emeline Sautron, Norbert Rolland, Patrice Catty, Daphné Seigneurin-Berny, et al. Structural Insights into the Nucleotide-Binding Domains of the P1B-type ATPases HMA6 and HMA8 from Arabidopsis thaliana. PLoS ONE, 2016, 11 (11), pp.e0165666. 10.1371/journal.pone.0165666. hal-01426929

HAL Id: hal-01426929

https://hal.science/hal-01426929v1

Submitted on 27 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH ARTICLE

# Structural Insights into the Nucleotide-Binding Domains of the P₁B-type ATPases HMA6 and HMA8 from *Arabidopsis thaliana*

Hubert Mayerhofer¹, Emeline Sautron², Norbert Rolland², Patrice Catty³, Daphné Seigneurin-Berny², Eva Pebay-Peyroula¹, Stéphanie Ravaud¹*

1. Institut de Biologie Structurale (IBS), Université Grenoble Alpes, CEA, CNRS, Grenoble, France,

2. Laboratoire de Physiologie Cellulaire et Végétale (LPCV), Institut de Biosciences et Biotechnologies de Grenoble (BIG), Université Grenoble Alpes, CEA, CNRS, INRA, Grenoble, France, 3. Laboratoire de Chimie et Biologie des Métaux (LCBM), BIG, Université Grenoble Alpes, CEA, CNRS, INRA, Grenoble, France

* Stephanie.Ravaud@ibs.fr

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Mayerhofer H, Sautron E, Rolland N, Catty P, Seigneurin-Berny D, Pebay-Peyroula E, et al. (2016) Structural Insights into the Nucleotide-Binding Domains of the P₁B-type ATPases HMA6 and HMA8 from *Arabidopsis thaliana*. PLoS ONE 11(11): e0165666. doi:10.1371/journal.pone.0165666

**Editor:** Oleg Y. Dmitriev, University of Saskatchewan, CANADA

**Received:** June 26, 2016

**Accepted:** October 14, 2016

**Published:** November 1, 2016

**Copyright:** © 2016 Mayerhofer et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All structures coordinates and reflections files are available from the PDB database (accession numbers 5LBD and 5LBK).

**Funding:** HM was supported by two post-doc fellowships, first from the European Community's Seventh Framework Programme FP7/2007–2013 under grant agreement no. HEALTH- F4-2007- 201924, EDICT Consortium and then from the Grenoble Alliance for Integrated Structural Cell Biology GRAL (ANR-10-LABX-49-01). EPP was

## Abstract

Copper is a crucial ion in cells, but needs to be closely controlled due to its toxic potential and ability to catalyse the formation of radicals. In chloroplasts, an important step for the proper functioning of the photosynthetic electron transfer chain is the delivery of copper to plastocyanin in the thylakoid lumen. The main route for copper transport to the thylakoid lumen is driven by two P₁B-type ATPases, **Heavy Metal ATPase 6 (HMA6)** and **HMA8**, located in the inner membrane of the chloroplast envelope and in the thylakoid membrane, respectively. Here, the crystal structures of the nucleotide binding domain of HMA6 and HMA8 from *Arabidopsis thaliana* are reported at 1.5Å and 1.75Å resolution, respectively, providing the first structural information on plants Cu⁺-ATPases. The structures reveal a compact domain, with two short helices on both sides of a twisted beta-sheet. A double mutant, aiding in the crystallization, provides a new crystal contact, but also avoids an internal clash highlighting the benefits of construct modifications. Finally, the histidine in the HP motif of the isolated domains, unable to bind ATP, shows a side chain conformation distinct from nucleotide bound structures.

## Introduction

Heavy metals, with their unique chemical properties, play an essential role in cellular processes as cofactors, structural stabilizers or redox partners and plants have evolved to maintain physiological concentrations of crucial metal ions. Their concentration and distribution needs to be tightly regulated as a large fraction of the proteome depends on them, but also to avoid their toxic effects when present in excess. Therefore, a number of cellular mechanisms exist ensuring their proper handling ranging from uptake and trafficking to storage and sequestration. As metal ions cannot diffuse freely across membranes, transporters play an important role in their targeting to different subcellular compartments.
supported by the Institut Universitaire de France.
SR was supported by a chair from the French
Centre National de la Recherche Scientifique and
from the University Grenoble Alpes. These funders
had no role in study design, data collection and
analysis, decision to publish, or preparation of the
manuscript.

**Competing Interests:** The authors have declared
that no competing interests exist.

P-type ATPases are transmembrane proteins that play a key role in the uphill transport of a wide range of cations across membranes using the energy provided by the hydrolysis of ATP, being crucial for the ion homeostasis in nearly all organisms and for heavy metal detoxification [1]. The catalytic cycle follows the classical Post-Albers model [2] and consists of four main steps with an alternating opening of the channel on the different sites of the membrane and the transport of the cation through extensive conformational changes, coupling ATP hydrolysis and cation transport. It includes the binding of ATP by the nucleotide (N) binding domain and its proper positioning for the transient phosphorylation of an aspartate within the invariant cytoplasmic DKTGT motif, part of the Phosphorylation (P) domain, and the subsequent dephosphorylation assisted by the Actuator (A) domain [3]. Conformational changes not only occur at the catalytic part of the transporter but also at its membrane domain affecting the access and the affinity of the transport site and leading to the motion of an ion. The P-type ATPase family can be divided into 5 groups [4]. The P<sub>1B</sub> subfamily, one of the largest, was shown to transfer transition metals (Cu<sup>+</sup>, Ag<sup>+</sup>, Cu<sup>2+</sup>, Zn<sup>2+</sup>...) and can be further divided into five subgroups with varying specificities. They consist of eight transmembrane (TM) helices, three cytoplasmic domains (A, N and P domain) and of a varying number of metal-binding domains at one of the termini. Characteristic features of the subfamily comprise a C/SPC motif found in TM6, cytosolic domains that are considered to be smaller than for other subfamilies [5], a particular membrane topology and conserved residues important for ATP binding, including a HP motif [6], different from the other families, while residues conserved elsewhere are not preserved.

Copper (Cu) is redox-active and plays an important role in electron-transfer reactions including the photosynthetic electron transport. It is also involved in the formation of reactive oxygen radicals and its free form is therefore tightly regulated, with virtually no free copper available in a cell. In plants, about half copper is present in the chloroplast either bound to the Cu/Zn superoxide dismutase (in the stroma) or to the plastocyanin located in the thylakoid lumen, with both proteins acting also as a copper sink in cases of abundant supply [7]. Plastocyanin is a small soluble protein, transferring electrons between cytochrome *f* of the cytochrome *b*<sub>6</sub>/*f* complex and the reaction centre of photosystem I. Two isoforms of plastocyanin are found in *Arabidopsis thaliana* [8]. Eight members of the heavy metal P<sub>1B</sub>-type ATPases subclass are present in *A. thaliana*, three of them being located into the chloroplast [9]. HMA1 and HMA6 are found in the inner membrane of the chloroplast envelope and HMA8 in the thylakoid membrane. HMA1 belongs to subgroup P<sub>1B-4</sub>, while HMA6 and HMA8 belong to subgroup P<sub>1B-1</sub>. HMA6 represents the main route for copper into the chloroplast, also supplying the Cu/Zn superoxide dismutase [10,11], while HMA1 represents a second route for copper, crucial under high light [12]. HMA8 supplies plastocyanin in the thylakoid lumen with copper, which also controls the stability of HMA8, but not of HMA6 [13,14]. HMA6 and HMA8 are both needed for an efficient photosynthetic electron transport [14]. However, they exhibit differing catalytic properties, HMA8 having a higher affinity for Cu and lower dephosphorylation kinetics than HMA6 [15].

Structures of different members of the P-type ATPase family exhibit the same overall fold, but differences exist between the subgroups [3] and, up to date, only few structures of the P<sub>1B</sub> subfamily are available [16–18]. Structural information for plant proteins from the subfamily is basically missing. Here, we present two structures of the nucleotide binding domain of HMA6 and HMA8, the former at the highest resolution, obtained for any P-type ATPase structure. They exhibit a very compact fold forming the minimum ATP binding unit. The conserved residues specific to the subfamily are spatially well aligned to other structures, yet a shift in the position of the conserved His in the HP motif could explain the inability of the isolated domain to bind ATP and can be a way of regulating the nucleotide binding. The surfaces surrounding
the entrance of the nucleotide-binding cavity exhibit notable differences in their charge distribution, which might affect the binding of ATP and subsequently the copper pumping rate.

## Materials and Methods

### Cloning and purification of HMA6<sup>N</sup> and HMA8<sup>N</sup>

The cDNA sequences coding for the nucleotide-binding domains of HMA6 (HMA6<sup>N</sup>, residues 607–734) and HMA8 (HMA8<sup>N</sup>, residues 557–689) were subcloned from the full-length constructs. The sequences of the primers used were: HMA6-forward 5' - CAGGGCGCCAGTCA CCCTGTTGTGACTGAAGTTATTATTCCTG-3', HMA6-reverse 5' - GACCCGACGCGGTTA TTCAAAACGGATGACTGCAGCAAG-3', HMA8-forward 5' - CAGGGCGCCAGTAGACCT GTCGTCTCTGGTGTTG-3' and HMA8-reverse 5' - GACCCGACGCGGTTAATCGGATATT GCAATAGCACCAATGATCC-3'.

The primers contain appropriate extensions for ligation independent cloning (LIC). The gene was inserted into a pETM-11/LIC vector (A. Geerlof, Helmholtz Zentrum München) via LIC cloning. The vector contains an *N*-terminal His6-tag followed by a Tobacco Etch Virus protease (TEV) cleavage site in case of HMA8<sup>N</sup> or an GB1 expression and solubility enhancer [19] flanked by a *N*-terminal His6-tag and a TEV cleavage site in case of HMA6<sup>N</sup>, with a three amino acid linker (GAS) that remains at the *N*-terminus of both proteins after TEV cleavage. Two glutamate residues were mutated to alanine in the HMA6<sup>N</sup> construct (EE709AA) using the QuikChange site-directed mutagenesis kit (Stratagene), yielding the HMA6<sup>N</sup>-AA construct. The final constructs were verified by DNA sequencing.

The plasmids were transformed into the *E. coli* Rosetta strain. A colony from freshly transformed cells was used to inoculate 5 mL lysogeny broth medium containing kanamycin (50 µg. mL⁻¹) and grown overnight at 310K. The overnight HMA8<sup>N</sup> culture was diluted to 1/1000th into 2 L of terrific broth medium containing kanamycin. They were grown to optical densities between 0.6–0.7 at 600 nm after which the temperature was lowered to 293K. The bacteria were induced using 0.05 mM isopropyl β-D-1-thiogalactopyranoside. The overnight HMA6<sup>N</sup>-AA culture was harvested at 2000 g for 5 min followed by resuspension in M9 minimal media. This was used to inoculate 2 L of M9 media supplemented with 5 g.L⁻¹ glucose, 50 µg.L⁻¹ kanamycin, 1 mg.L⁻¹ thiamine and 1 mg.L⁻¹ biotin. At an OD<sub>600</sub> of 0.5, 100 mg.L⁻¹ each of phenylalanine, lysine and threonine, 50 mg.L⁻¹ each of leucine, isoleucine, valine and l-SeMet were added [20]. At an OD<sub>600</sub> of 0.6 the culture was induced with 0.15 mM IPTG after which the temperature was lowered to 301K. Cultivation was continued in both cases for 16 h and cells were harvested by centrifugation at 5000 rpm in a JLA-8.1000 rotor for 20 min at 277K. The cell pellets were stored at 253K. After thawing on ice, pellets were re-suspended in lysis buffer (20 mM Mops pH 7.0, 250 mM NaCl, 10% (w/v) glycerol, 20 mM imidazole, 1 mM complete protease inhibitor cocktail EDTA-free (Roche), 1 mg.mL⁻¹ DNAse (NEB) and 0.1% (w/v) CHAPS (Euromedex)) in a final volume of 20 mL per 5 g of cells and lysed by pulsed sonication on ice for two times 90 sec. The lysates were centrifuged at 37,000 rpm in a Ti-45 rotor at 277K for 60 min. The supernatants were loaded onto a 5 mL HisTrap HP column (GE Healthcare), which had been equilibrated against buffer A (20 mM Mops pH 7.0, 250 mM NaCl, 5% (w/v) Glycerol). The column was washed with 5 column volumes (CV) of buffer A, followed by 5 CV of buffer A with 20% of buffer B (buffer A with 250 mM imidazole). The proteins were eluted with a gradient of 20% to 100% of buffer B in buffer A (50 mM to 250 mM imidazole) within 13 CV. Fractions containing HMA8<sup>N</sup> or HMA6<sup>N</sup>-AA were pooled and dialysed overnight at 277K against buffer C (20 mM Mops pH 7.0, 150 mM NaCl and 6 mM β-ME) or buffer D (20 mM potassium phosphate buffer pH 7.0, 150 mM NaCl and 6 mM β-ME), respectively, in a dialysis bag (molecular-weight cut-off 5 kDa) with 1 mg TEV protease (containing a non-
cleavable *N*-terminal His-tag) added per 30 mg of protein. The solution was collected the next day and passed again over the Ni-NTA column equilibrated against buffer C or D. The flow-through was collected and concentrated using a Vivaspin column (molecular-weight cutoff 5kDa). Size-exclusion chromatography (HiLoad 26/60 Superdex 75, Amersham Biosciences) was used as a final step of purification. The column was pre-equilibrated in buffer E (20 mM Mops pH 7.0, 150 mM NaCl) in case of HMA8N or buffer E (20 mM potassium phosphate buffer pH 7.0, 150 mM KCl) in case of HMA6N-AA. The samples eluted as a single peak and peak fractions were analysed by SDS-PAGE and pooled.

## Crystallization of HMA8N

HMA8N was concentrated to 35 mg.mL⁻¹ using a Vivaspin column (molecular-weight cut-off 5kDa). Initial crystallization trials were carried out with 5 different 96 well screens (Wizard I +II (Rigaku Reagents), Crystal Screen Lite & PEG/Ion and Grid screen (both Hampton Research), JCSG and PACT (both Qiagen) screens) at the High-throughput crystallisation laboratory (HTX Lab, EMBL, PBS, Grenoble). All initial screens were performed at 293K in 96-well Greiner plates using the sitting-drop vapour-diffusion method. 100 nL protein solution was mixed with 100 nL reservoir solution and equilibrated against 50 µL reservoir solution.

Crystallization trials with His tag-free HMA8N resulted in a number of hits. Initial lead conditions were optimized through the use of custom made screens. The best crystals of HMA8N were observed at 293K by mixing 1 µL of protein with 1 µL of 0.1 M sodium acetate trihydrate pH 4.4 and 0.6 M sodium formate. Co-crystallization and soaking trials using 5–10 mM ATP were also performed in the same crystallization conditions.

## Crystallization of HMA6N-AA

HMA6N-AA was concentrated to 10 mg.mL⁻¹ using a Vivaspin column (molecular-weight cut-off 5 kDa). Initial crystallization trials were carried out with 4 different 96 well screens (Wizard I+II (Rigaku Reagents), JCSG and PEGs (both Qiagen)) screens at the High-throughput crystallisation laboratory (HTX Lab, EMBL, PBS, Grenoble) [21]. All initial screens were performed at 293K in 96-well Greiner plates using the sitting-drop vapour-diffusion method. 100 nL protein solution was mixed with 100 nL reservoir solution and equilibrated against 50 µL reservoir solution.

Crystallization trials with His tag-free HMA6N-AA resulted in one hit. Initial lead conditions were optimized through the use of custom made screens, which involved microseeding, using seeds obtained from the initial hit, and varying the protein concentration. The best crystals of HMA6N-AA were observed at 293K by mixing 1 µL of protein at 4 mg.mL⁻¹ with 1 µl of 0.1 M sodium acetate trihydrate pH 4.4 and 0.6 M sodium formate. The plate was incubated at 293K for 4 h after which the wells were reopened and, using a cat whisker, microseeds were added. Crystals appeared after 2 days and reached their final size after about 1 week.

## Data collection and processing

Prior to data collection, crystals of HMA8N were cryoprotected with Paratone-N, while crystals of HMA6N-AA were cryoprotected by adding ethylene glycol to the reservoir solution to a final concentration of 27% (w/v). The crystals were then flash-cooled to 77K in liquid nitrogen. A complete X-ray diffraction data set was collected on beamline BM14 at the European Synchrotron Radiation Facility (ESRF, Grenoble, France) using a MARCCD 225 detector in case of HMA8N. A total of 150 frames were collected with a rotation range of 1° and a detector distance of 187.25 mm. For HMA6N-AA complete X-ray diffraction data sets corresponding to the peak and inflection point were collected on beamline ID23-1 (ESRF, Grenoble, France)
using a PILATUS 6M-F detector. A total of 1210 frames were collected with a rotation range of 0.1° and a detector distance of 271.22 mm. The data were indexed and integrated using *XDS* [22] and scaled with *Aimless* [23].

## Structure solution and refinement

The structure of HMA8<sup>N</sup> was solved by molecular replacement using the N-domain of the CopA structure from *A. fulgidus* as a search model (pdb file 3A1C) in PHASER [24]. The model was manually modified in COOT [25] and refined using Phenix 1.10_2155 [26] by TLS refinement with 22 TLS groups that had been determined by the program [27]. Multiwave-length anomalous dispersion phasing using the SeMet-substituted crystal solved the structure of HMA6<sup>N</sup>-AA. SeMet sites were found with SHELXD and density modifications were performed with the program SHELX E [28]. Using the traceable map automated model building was carried out using buccaneer [29], which traced 80% of the structure. Further building was performed cycling between manual modelling in COOT [25] and refinement in Phenix 1.10_2155 [26], using only the dataset collected at the peak wavelength. MolProbity [30] analysis indicated that the overall geometry of the final model of HMA8<sup>N</sup> and HMA6<sup>N</sup>-AA ranked in the 97<sup>th</sup> and 84<sup>th</sup> percentile, respectively, where the 100<sup>th</sup> percentile is best among structures of comparable resolution. All structure-related figures were generated with PYMOL. Data collection and refinement statistics for HMA6<sup>N</sup>-AA and HMA8<sup>N</sup> are shown in Tables 1 and 2, respectively. The coordinates are deposited in the PDB under the accession codes 5LBD and 5LBK corresponding to HMA6<sup>N</sup> and HMA8<sup>N</sup>, respectively.

## Isothermal titration calorimetry

Calorimetric measurements were tested at 20°C or 25°C under continuous stirring at 286 rpm to ensure rapid mixing, using a VP-ITC instrument with a cell volume of 1.4569 ml (MicroCal, LLC). After concentration HMA6<sup>N</sup> and HMA8<sup>N</sup> were dialysed against the buffer used for size-exclusion chromatography with or without the addition of 1 mM MgCl<sub>2</sub>. The same batch of buffer was used for further dilutions of the proteins and for the dilution of the nucleotide. The protein concentrations used in the cell were between 130 µM–250 µM for HMA6<sup>N</sup> and between 18 µM–250 µM for HMA8<sup>N</sup>. As binding partner ATP was tested at concentrations between 500 µM and 8.6 mM.

## Results and Discussion

### Construct optimization and structure determination

The N-domain of HMA8 (HMA8<sup>N</sup>) was readily expressed and purified, representing monomeric protein, which eluted as a symmetric peak with a calculated mass of 17.2 kDa. The structure of HMA8<sup>N</sup> was solved at 1.75 Å by molecular replacement, using the N-domain of pdb 3A1C as a search model, with two molecules in the asymmetric unit (AU).

The initial construct (HMA6<sup>N</sup>) of the HMA6 N-domain was poorly expressed and required the addition of an N-terminal solubility enhancing tag, which lead to a 15-fold increase in expression. After cleavage of the tag, highly soluble and homogenous protein was obtained, which eluted from the size-exclusion column as a monomer, with a calculated mass of 15.4 kDa, yet did not yield crystals. Modifications of the construct and the buffer remained unsuccessful and therefore the construct sequence was analysed using the Surface-Entropy-Reduction prediction (SERp) server [31] for patches of surface residues, which, once mutated, could potentially promote local stability, and the three highest scoring mutants were tested for expression and crystallisation. One of them, the double mutant EE709AA, dubbed construct
<table><caption>Table 1. Data collection and refinement statistics of HMA6N-AA.</caption><thead><tr><th></th><th>Peak</th><th>Inflection</th><th>Remote</th></tr></thead><tbody><tr><td>Data collection</td><td></td><td></td><td></td></tr><tr><td>Diffraction source</td><td>ID23-1</td><td>ID23-1</td><td>ID23-1</td></tr><tr><td>Wavelength (Å)</td><td>0.97924</td><td>0.97947</td><td>0.97692</td></tr><tr><td>Temperature (K)</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Crystal-detector distance (mm)</td><td>271.40</td><td>271.22</td><td>272.21</td></tr><tr><td>Rotation range per image (°)</td><td>0.1</td><td>0.1</td><td>0.1</td></tr><tr><td>Total rotation range (°)</td><td>121</td><td>121</td><td>121</td></tr><tr><td>Space group</td><td>P 2<sub>1</sub> 2<sub>1</sub></td><td>P 2<sub>1</sub> 2<sub>1</sub></td><td>P 2<sub>1</sub> 2<sub>1</sub></td></tr><tr><td>a, b, c</td><td>38.77 65.15 85.6</td><td>38.79 65.23 85.64</td><td>38.84 65.30 85.68</td></tr><tr><td>α, β, γ</td><td>90 90 90</td><td>90 90 90</td><td>90 90 90</td></tr><tr><td>Mosaicity (°)</td><td>0.062</td><td>0.068</td><td>0.083</td></tr><tr><td>Resolution range</td><td>50–1.50 (1.53–1.50)</td><td>50–1.50 (1.53–1.50)</td><td>50–1.50 (1.53–1.50)</td></tr><tr><td>Total No. of reflections</td><td>155329 (7508)</td><td>155325 (7597)</td><td>155309 (7467)</td></tr><tr><td>No. of unique reflections</td><td>35341 (1689)</td><td>35433 (1707)</td><td>35474 (1689)</td></tr><tr><td>Completeness (%)</td><td>99.5 (99.0)</td><td>99.6 (99.4)</td><td>99.5 (97.9)</td></tr><tr><td>Redundancy</td><td>4.4 (4.4)</td><td>4.4 (4.5)</td><td>4.4 (4.4)</td></tr><tr><td>&lt;I/σ(I)&gt;</td><td>13.4 (1.4)</td><td>13.6 (1.5)</td><td>15.7 (1.1)</td></tr><tr><td>Refinement statistics</td><td></td><td></td><td></td></tr><tr><td>R factor</td><td>18.54 (32.23)</td><td></td><td></td></tr><tr><td>Rfree</td><td>20.72 (34.6)</td><td></td><td></td></tr><tr><td>No. of non-H atoms</td><td></td><td></td><td></td></tr><tr><td>Protein</td><td>1761</td><td></td><td></td></tr><tr><td>Water</td><td>149</td><td></td><td></td></tr><tr><td>Mean B factors (Å<sup>2</sup>)</td><td></td><td></td><td></td></tr><tr><td>Protein</td><td>30.86</td><td></td><td></td></tr><tr><td>Ligands</td><td>44.8</td><td></td><td></td></tr><tr><td>Solvent</td><td>38.4</td><td></td><td></td></tr><tr><td>R.m.s. deviation from ideal</td><td></td><td></td><td></td></tr><tr><td>Bond length (Å)</td><td>0.01</td><td></td><td></td></tr><tr><td>Bond angles (°)</td><td>0.922</td><td></td><td></td></tr><tr><td>Ramachandran plot (%)</td><td></td><td></td><td></td></tr><tr><td>Residues in favoured region</td><td>97.76</td><td></td><td></td></tr><tr><td>Residues in allowed region</td><td>2.24</td><td></td><td></td></tr><tr><td>Outliers</td><td>0</td><td></td><td></td></tr><tr><td>Molprobity Score</td><td>1.52</td><td></td><td></td></tr><tr><td>PDB code</td><td>5LBD</td><td></td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0165666.t001

HMA6N-AA, led, together with a change of the buffer, to an initial crystal hit. Yet the crystals proved difficult to optimize, grew in clusters, required microseeding in order to reproduce the hit and only diffracted to 3.2 Å, with six molecules in the AU in space group P3221. A reduction of the protein concentration to 4 mg.mL⁻¹ occasionally led to crystals with a different, more wedge like shape and a different space group (P22121), not showing any more the inter-grown clusters seen initially. These crystals diffracted up to 1.5 Å resolution and contained only two molecules in the AU. For both proteins, two molecules are present in the AU, yet they are monomeric in solution as found by size exclusion chromatography. The HMA6 and HMA8 nucleotide binding and phosphorylation domains (N/P domain, HMA6N/P, residues 576–853 and HMA8N/P, residues 526–806) were expressed and purified, yet in both cases a fraction of the protein was found to be C-terminally truncated. While the full N/P domains were dimeric
<table><caption>Table 2. Data collection and refinement statistics of HMA8<sup>N</sup>.</caption><thead><tr><th></th><th></th></tr></thead><tbody><tr><td><strong>Data collection</strong></td><td></td></tr><tr><td>Diffraction source</td><td>BM14 (ESRF)</td></tr><tr><td>Wavelength (Å)</td><td>0.953</td></tr><tr><td>Temperature (K)</td><td>100</td></tr><tr><td>Detector</td><td>MARCCD 225</td></tr><tr><td>Crystal-detector distance (mm)</td><td>160.85</td></tr><tr><td>Rotation range per image (°)</td><td>1.0</td></tr><tr><td>Total rotation range (°)</td><td>150</td></tr><tr><td>Space group</td><td>P2<sub>1</sub>2<sub>1</sub>2<sub>1</sub></td></tr><tr><td>a, b, c</td><td>44.157 48.257 108.295</td></tr><tr><td>α, β, γ</td><td>90, 90, 90</td></tr><tr><td>Mosaicity (°)</td><td>0.134</td></tr><tr><td>Resolution range</td><td>34.22–1.756 (1.86–1.756)</td></tr><tr><td>Total No. of reflections</td><td>141918 (22757)</td></tr><tr><td>No. of unique reflections</td><td>23737 (3786)</td></tr><tr><td>Completeness (%)</td><td>99.4</td></tr><tr><td>Redundancy</td><td>6.0 (6.0)</td></tr><tr><td>&lt;I/σ(I)&gt;</td><td>16.5 (2.64)</td></tr><tr><td><strong>Refinement statistics</strong></td><td></td></tr><tr><td>R factor</td><td>19.81</td></tr><tr><td>R<sub>free</sub></td><td>23.76</td></tr><tr><td>No. of non-H atoms</td><td></td></tr><tr><td>Protein</td><td>1787</td></tr><tr><td>Water</td><td>125</td></tr><tr><td>Mean B factors (Å<sup>2</sup>)</td><td></td></tr><tr><td>Protein</td><td>43</td></tr><tr><td>Water</td><td>43.6</td></tr><tr><td>R.m.s. deviation from ideal</td><td></td></tr><tr><td>Bond length (Å)</td><td>0.013</td></tr><tr><td>Bond angles (°)</td><td>1.103</td></tr><tr><td>Ramachandran plot (%)</td><td></td></tr><tr><td>Residues in favoured region</td><td>97.4</td></tr><tr><td>Residues in allowed region</td><td>2.6</td></tr><tr><td>Outliers</td><td>0</td></tr><tr><td>Molprobity Score</td><td>1.35</td></tr><tr><td>PDB code</td><td>5LBK</td></tr></tbody></table>

doi:10.1371/journal.pone.0165666.t002

in solution the truncated fragments were monomeric. Upon removal of the truncated fraction, the full-length constructs were stable, yet did not crystallize.

## Structure of HMA6 and HMA8 in comparison to others

The overall structure of HMA6<sup>N</sup>-AA and HMA8<sup>N</sup> (Fig 1A and 1B, respectively) exhibits a twisted, antiparallel β-sheet made of six strands being flanked by two pairs of short α-helices, one on the concave (α1 and α2) and one on the convex (α3 and α4) side of the sheet (Fig 1C). The overall fold is similar to those found in other P-type ATPases [17,32] (Tables 3 and 4). While the central β-sheet is highly conserved within the subfamily, the main differences, as in the Menkes and Wilson proteins [33,34], arise in the loop length and the number of helices
<figure><img src="image_4.png" /><figcaption>**Fig 1.** 3D structures of HMA6 and HMA8 N-domains. **A** The structure of HMA6<sup>N</sup>-AA with β-strands shown in blue and α-helices in red. Missing loops are outlined by dashed lines. **B** The structure of HMA8<sup>N</sup>; colouring is as in **A**. **C** Superimposition of the two structures with HMA6<sup>N</sup>-AA in red and HMA8<sup>N</sup> in blue.</figcaption></figure>

doi:10.1371/journal.pone.0165666.g001

(Figs 2 and 3). In HMA6<sup>N</sup>-AA two loops are disordered—L1 between β1 and α1 (res. 617–625) and L5 between α4 and β5 (res. 712–716, with α4 partly and completely disordered in chain A and B, respectively), while in HMA8<sup>N</sup> only one loop is flexible (res. 656–668). The first missing loop in HMA6 corresponds to an insertion not present in most other P<sub>IB</sub>-type ATPases (Fig 2). The position of the second missing loop is identical in both structures, differing only in the length, as in HMA8 an insertion is present, not typically found in the subfamily (Fig 2). The position of helix α4, which precedes the loop and is also partially disordered, relative to a helix
<table><caption>Table 3. A matrix of sequence identity values among the P<sub>1B</sub>-type ATPase N-domains compared to HMA (1, 6 and 8) N-domains.</caption><thead><tr><th></th><th>HMA1<sup>N</sup></th><th>HMA6<sup>N</sup></th><th>HMA8<sup>N</sup></th><th>LpCopA</th><th>AfCopA</th><th>AfCopB</th><th>SsCopB</th><th>ATP7A</th><th>ATP7B</th></tr></thead><tbody><tr><td>HMA1N</td><td>-</td><td>29</td><td>23</td><td>22.6</td><td>20.3</td><td>25.4</td><td>18.3</td><td>15.6</td><td>18.9</td></tr><tr><td>HMA6N</td><td>29</td><td>-</td><td>34.2</td><td>28.5</td><td>27.7</td><td>28.7</td><td>26.1</td><td>27.5</td><td>23.3</td></tr><tr><td>HMA8N</td><td>23</td><td>34.2</td><td>-</td><td>23.1</td><td>31.3</td><td>32.8</td><td>31.2</td><td>23.3</td><td>24.2</td></tr><tr><td>LpCopA</td><td>22.6</td><td>28.5</td><td>23.1</td><td>-</td><td>32.7</td><td>34.2</td><td>28.6</td><td>26.5</td><td>24.8</td></tr><tr><td>AfCopA</td><td>20.3</td><td>27.7</td><td>31.3</td><td>32.7</td><td>-</td><td>40.5</td><td>35.2</td><td>31.3</td><td>29.6</td></tr><tr><td>AfCopB</td><td>25.4</td><td>28.7</td><td>32.8</td><td>34.2</td><td>40.5</td><td>-</td><td>34.7</td><td>26.3</td><td>23.9</td></tr><tr><td>SsCopB</td><td>18.3</td><td>26.1</td><td>31.2</td><td>28.6</td><td>35.2</td><td>34.7</td><td>-</td><td>26</td><td>28.4</td></tr><tr><td>ATP7A</td><td>15.6</td><td>27.5</td><td>23.3</td><td>26.5</td><td>31.3</td><td>26.3</td><td>26</td><td>-</td><td>54</td></tr><tr><td>ATP7B</td><td>18.9</td><td>23.3</td><td>24.2</td><td>24.8</td><td>29.6</td><td>23.9</td><td>28.4</td><td>54</td><td>-</td></tr></tbody></table>

doi:10.1371/journal.pone.0165666.t003

in the P-domain is indicative of the conformation, open or closed, of the two domains. While not specific contacts are present in structures where the N- and P-domain are present, the P-domain could aid in stabilizing this disordered region. Furthermore, the insertion found in HMA8 could interact tighter with the P-domain and report on the N-domain inclination. In the two HMA structures helix, α4 is more tightly packed against the β-sheet than in other members of the subfamily, which doesn't seem to be caused by the absence of the P-domain, as other isolated N-domains show a less tight packing [33,34]. Within the AU, both structures are highly similar with an rmsd of 0.43 Å and 0.33 Å for HMA6<sup>N</sup>-AA and HMA8<sup>N</sup>, respectively. When comparing the two proteins they also show a high similarity with an rmsd of about 1 Å (Fig 1).

## Effects of surface mutants

The effect of the double alanine mutant in HMA6, proposed by the SERp server, seems to have two aspects. Glu709A, in the partially disordered helix α4, points to the solvent and the glutamate of the wild type sequence, when modelled in the structure, clashes with the side chain of Leu705 of the symmetry related chain A* (Fig 4). Likewise, the same clash of Glu709 of chain A* with Leu705A would occur. On the other hand, the side chain of Glu710 points towards the centre of the domain and would clash with Val730 of β6. If not mutated glutamate 709 and 710 would be 1.2–1.5 Å from the next side chain, with the clash preventing the crystal formation. Several crystal contacts are present between the chains in two axis, while these mutations allowed for a crystal contact in the third axis, between α4 and α4'. The same region in chain B, including helix α4 (704–716), is disordered, potentially caused by the lack, unlike for chain A, of a stabilizing interaction with a symmetry related chain.

<table><caption>Table 4. A matrix of RMSD values (in Å) among the N-domains of structural homologs of HMA6<sup>N</sup> and HMA8<sup>N</sup> compared in Fig 3 using the CA positions.</caption><thead><tr><th></th><th>HMA6<sup>N</sup> 5LBD</th><th>HMA8<sup>N</sup> 5LBK</th><th>LpCopA 3RFU</th><th>AfCopA 2B8E</th><th>AfCopB 3SKX</th><th>SsCopB 2IYE</th><th>ATP7B 2ARF</th></tr></thead><tbody><tr><td>HMA6<sup>N</sup> 5LBD</td><td>-</td><td>0.73</td><td>1.04</td><td>0.77</td><td>1.35</td><td>1.14</td><td>4.02</td></tr><tr><td>HMA8<sup>N</sup> 5LBK</td><td>0.73</td><td>-</td><td>1.19</td><td>2.9</td><td>0.92</td><td>0.87</td><td>5.86</td></tr><tr><td>LpCopA 3RFU</td><td>1.04</td><td>1.19</td><td>-</td><td>1.26</td><td>1.32</td><td>1.22</td><td>3.6</td></tr><tr><td>AfCopA 2B8E</td><td>0.77</td><td>2.9</td><td>1.26</td><td>-</td><td>1.51</td><td>1.95</td><td>3.44</td></tr><tr><td>AfCopB 3SKX</td><td>1.35</td><td>0.92</td><td>1.32</td><td>1.51</td><td>-</td><td>1.05</td><td>4.39</td></tr><tr><td>SsCopB 2IYE</td><td>1.14</td><td>0.87</td><td>1.22</td><td>1.95</td><td>1.05</td><td>-</td><td>4.01</td></tr><tr><td>ATP7B 2ARF</td><td>4.02</td><td>5.86</td><td>3.6</td><td>3.44</td><td>4.39</td><td>4.01</td><td>-</td></tr></tbody></table>

doi:10.1371/journal.pone.0165666.t004
<figure><img src="image_5.png" /><figcaption>**Fig 2. Sequence alignment of P<sub>1B</sub>-type ATPase N-domains.** Alignment of the region covering the N-domains of HMA1, HMA6 and HMA8 from *A. thaliana*, CopA and CopB from *A. fulgidus*, CopA from *L. pneumophila*, CopB from *S. solfataricus*, the human Menkes (ATP7A) and Wilson proteins (ATP7B) and ZntA from *S. sonnei*. Disordered regions in HMA6 and HMA8 are marked by a red rectangle. The mutated residues in HMA6 are marked in orange. The secondary structures of HMA6 and HMA8 are indicated above and below the alignment, respectively. Multiple sequence alignment was done using MUSCLE (http://www.ebi.ac.uk/Tools/msa/muscle) and visualized using ESPript 3.0 [35].</figcaption></figure>

doi:10.1371/journal.pone.0165666.g002

ATP binding properties

ATP binding to the isolated N-domain of HMA6 and HMA8 was assessed by ITC in different experimental conditions (various proteins and ATP concentrations, presence or absence of magnesium). None of them showed a signal of interaction between HMA6<sup>N</sup> or HMA8<sup>N</sup> and ATP. In parallel, several attempts to obtain the structure of the HMA N-domains in complex with ATP were made by soaking or co-crystallization on HMA8<sup>N</sup>. They led to diffracting crystals but no density of the nucleotide could be observed. An accurate review of the literature shows that nucleotide binding studies of isolated N-domains of P-type ATPases are scarce. A few data have been reported such as on the rat Na<sup>+</sup>K<sup>+</sup>-ATPase from rat [36], the rabbit skeletal muscle Ca<sup>2+</sup>-ATPase SERCA1a [37] and the K<sup>+</sup>-ATPase KdpB from *E. coli* [38]. Apparent ATP affinities measured on these N-domains are low (Kd of 5–25 mM for the Na<sup>+</sup>K<sup>+</sup>-ATPase; 0.7–2.4 mM for SERCA1a; 1.4 mM for KdpB), up to 10 times lower than the apparent affinity determined with the full-length proteins. For a given protein the range of Kd values is sometimes large since measured by different methods (Chemical shift, ITC, intrinsic fluorescence). These values do not depend on the presence of magnesium except for the Na<sup>+</sup>K<sup>+</sup>-ATPase N-domain that unexpectedly exhibits a lower affinity for ATP in presence of magnesium.
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Comparison of P<sub>1B</sub>-type ATPase N-domains of structural homologs. Superimposition of HMA6<sup>N</sup>-AA (gray), HMA8<sup>N</sup> (red), LpCopA (3RFU, green), AfCopA (2B8E, cyan) and AfCopB (3SKX, orange), SsCopB (2IYE, yellow) and the Wilson disease protein (2ARF, blue).</figcaption></figure>

doi:10.1371/journal.pone.0165666.g003

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Effect of the double alanine mutant in HMA6<sup>N</sup>. Interface of helix 4 from chain A and the symmetry related chain A* with glutamate modelled at position 709 and 710, matching the WT sequence instead of the two alanine, as found in the structure. Helix 4 is shown in orange in the symmetry related side chain marked with an asterisk.</figcaption></figure>

doi:10.1371/journal.pone.0165666.g004
Similar studies have also been carried out on the highly homologous human Cu<sup>+</sup>-ATPases
ATP7B [34,39–41] and ATP7A [33]. The N-domains of these two proteins were shown to dis-
play similar apparent Kd for ATP ranging from 10 to 110 µM. The data available on AfCopA
[17,42–44] and ssCopB [32] were obtained on the [N + P]-domains and therefore cannot be
strictly compared with the present work in terms of ATP binding.

A striking difference between the N-domains of the human and plant Cu⁺-ATPases is the presence of a large unfolded region found between β3 and β4 only found in ATP7A and ATP7B (Fig 2). Such a region is also missing in the N-domains of prokaryotic Cu⁺-ATPases such as AfCopA for which no ATP binding has been reported. This additional region found in the N-domains of mammalian Cu⁺-ATPases even though not directly involved in nucleotide binding could give some flexibility to the isolated domain allowing a measurable binding of ATP. Indeed, it has been shown in ATP7A that decrease in size of the adenine ring binding cavity caused by the motion of both α1- α2 and β2-β3 loops, allows a tight packing around the ATP molecule [33]. So, even faint, a conformation change of the isolated N-domain is required for ATP binding in ATP7A (and probably ATP7B).

A detailed analysis of the nucleotide binding site among the Cu⁺-ATPases is presented
below.

Nucleotide binding site analysis

The HMA N-domains structures were compared to the [N + P] domains structures of the two prokaryotic Cu⁺-ATPases AfCopA (pdb 3A1D) and SsCopB (pdb 2YJ4) that were solved in the presence of ATP or ADP. No major structural changes are expected upon the binding of ATP. A superimposition shows that the nucleotide binding site and the position of conserved residues is basically identical between the structures (Fig 5A) as well as a good fit of the binding pocket to the nucleotide. The three conserved glycine residues G675, G677 and G690 (HMA6 numbering) lining one side of the binding pocket and the glutamate residue E640 (HMA6 numbering), normally in contact with the N1 and N6 atoms of the nucleotide, are occupying

<figure><img src="image_8.png" /><figcaption>**Fig 5.** A Superimposition of the ATP binding site of HMA6N-AA (red), HMA8N (orange), SsCopB (blue, 2YJ4) and AfCopA (light-blue, 3A1D).
Conserved residues are shown as sticks, coloured as the corresponding chain and numbered according to their position in HMA6. ATP and ADP of SsCopB and AfCopA are shown as sticks with carbon atoms show in grey, nitrogen atoms in blue, oxygen atoms in red and phosphate atoms in orange. **B** Structure of HMA6N-AA and **C** HMA8N, each with the ATP of the SsCopB structure modelled. For the positioning of the nucleotide the N-domains of HMA6N-AA, HMA8N and SsCopB (2YJ4) were superimposed and as the conserved residues are spatially in close proximity a model of HMA6N-AA and HMA8N together with the ATP of SsCopB was generated. Conserved residues important for the binding are shown in sticks and the ATP as spheres. Colouring as in A.</figcaption></figure>

doi:10.1371/journal.pone.0165666.g005
positions ready for nucleotide binding. In contrast, histidine H645 (HMA6 numbering), part of the HP motif and found elsewhere in contact with the α- and β-phosphate [17], is pointing away from the binding site and would have to adopt a different position for proper interaction with the phosphate chain (Fig 5B and 5C). Of all available structures of the P<sub>1B</sub> subfamily, a part of the histidine residues is already in such a position, while the remaining adopt an orientation similar to the one observed in the HMA structures. Therefore, it seems that the proper orientation of the histidine is not a consequence of the nucleotide binding but rather a prerequisite, possibly aided by the P-domain, providing a conceivable explanation for our failure to observe ATP binding to the isolated HMA domains. In summary the modelled position of ATP in the N-domain (Fig 5B and 5C) most likely provides an accurate description, while changes in the side chain orientation of the histidine H645 (HMA6 numbering) are needed for proper interaction.

An interesting observation in the case of HMA6 is the presence of glutamate 672, not commonly found at this position, pointing to the nucleotide binding cavity and guiding the position of the conserved H645. While missing in HMA8, this interaction probably contributes to the inability of the isolated domain to bind ATP. Of interest is also a second salt bridge in the vicinity of the nucleotide binding domain, E673-R697 in HMA6 and E612-R637 in HMA8, linking β2 to α3. A salt bridge is also observed in LpCopA (pdb 3RFU), however the charges are inverted (K493-E514), while in the Menkes protein, only van der Waals contacts are present (pdb 2KMX). This bridge sits either next to, or interacts via, the backbone with the highly conserved glycine in β3 (677 in HMA6, 616 in HMA8 or 494 in LpCopA), part of a highly flat region of the β-sheet in contact with the adenine ring, and probably aides in positioning of this glycine. While the loop between β2 and β3 is not changing its position upon nucleotide binding in these structures, it is moving in the case of SsCopB (pdb 2IYE), which lacks helices 3 and 4 and is not stabilized by the mentioned salt-bridge. The salt bridge could therefore control the position of the two strands β2-β3 relative to the cavity. Interestingly, the charge distributions at the surface delineating the cavity differ from the known N-domain structures (S1 Fig). Electrostatic properties of this surface, might not only guide the entrance of the nucleotide to the cavity and the interactions with the P-domain, but also influence the kinetic properties of these pumps. Marked differences between ATP/ADP ratios in the cytosol and the chloroplast stroma [45], as well as between the pH of these two subcellular compartments (increase of the stromal alkalization on illumination), might also explain differences in electrostatic properties of this surface (HMA8 requires the efficient use of ATP at a low ATP/ADP ratio found in the stroma).

Although too few data are yet available on Cu<sup>+</sup>-ATPases to draw a general conclusion on their ATP binding mechanism, it sounds realistic to correlate the affinity of ATP to the ability of the isolated N-domain to accommodate efficiently the nucleotide and therefore to adapt its conformation. In that sense, our study reveals that the N-domains of AtHMA6 and AtHMA8 exhibit significant differences compared to those of ATP7A and ATP7B, the only Cu<sup>+</sup>-ATPases for which ATP binding was measured on the isolated N-domains. It also suggests that ATP binding mechanism of AtHMA6 and AtHMA8 (and by extension to prokaryotic Cu<sup>+</sup>-ATPases) might resemble more that of KdpB.

**Supporting Information**

**S1 Fig. Electrostatic potential surfaces of P1B-type ATPases.** Red denotes negative and blue positive charges, respectively. All structures are limited to the N-domains and oriented the same way, with the nucleotide binding site in the centre and the nucleotide, if present, shown as sticks. For each structure the pdb code is given. In summary the first structures of plant P1B-type ATPase N-domains show their high similarity to their bacterial counterparts, with
the binding amino acids in place for interaction with ATP. The central histidine of the HP
motif needs to be reoriented in order to interact with the adenine ring, providing a possible
clue why these isolated domains show a very weak affinity for ATP. In addition, a so far not
noted salt bridge linking the central β-sheet to α3 stabilizes the glycine. Finally, modifications
of the surface helped in obtaining crystals of the HMA6 N-domain that were otherwise unat-
tainable.

(TIF)

**Acknowledgments**

HM was supported by two post-doc fellowships, first from the European Community's Seventh Framework Programme FP7/2007-2013 under grant agreement no. HEALTH-F4-2007-201924, EDICT Consortium and then from the Grenoble Alliance for Integrated Structural Cell Biology GRAL (ANR-10-LABX-49-01). EPP was supported by the Institut Universitaire de France. SR was supported by a chair from the French Centre National de la Recherche Scientifique and from the University Grenoble Alpes. These funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We acknowledge the platforms of the Grenoble Instruct Center (ISBG; UMS 3518 CNRS-CEA-UGA-EMBL) supported by the French Infrastructure for Integrated Structural Biology Initiative FRISBI (ANR-10-INSB-05-02) and by the Grenoble Alliance for Integrated Structural Cell Biology GRAL (ANR-10-LABX-49-01) within the Grenoble Partnership for Structural Biology (PSB) The X-ray diffraction data were collected at the European Synchrotron Radiation Facility (ESRF) Grenoble on beamlines BM14 and ID23-1.

**Author Contributions**

**Conceptualization:** HM NR PC DSB EPP SR.

**Formal analysis:** HM ES EPP SR.

**Funding acquisition:** NR EPP SR.

**Supervision:** EPP SR.

**Writing – original draft:** HM SR.

**Writing – review & editing:** HM NR PC DSB EPP SR.

**References**

1. Kühlbrandt W. Biology, structure and mechanism of P-type ATPases. Nat Rev Mol Cell Biol. 2004; 5: 282–295. doi: 10.1038/nrm1354 PMID: 15071553

2. Post RL, Hegyvary C, Kume S. Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphatase. J Biol Chem. 1972; 247: 6530–6540. PMID: 4263199

3. Palmgren MG, Nissen P. P-Type ATPases. Annual Review of Biophysics. 2011; 40: 243–266. doi: 10.1146/annurev.biophys.093008.131331 PMID: 21351879

4. Axelsen KB, Palmgren MG. Evolution of substrate specificities in the P-type ATPase superfamily. J Mol Evol. 1998; 46: 84–101. PMID: 9419228

5. Bramkamp M, Altendorf K. Functional Modules of KdpB, the Catalytic Subunit of the Kdp-ATPase from Escherichia coli†. Biochemistry. 2004; 43: 12289–12296. doi: 10.1021/bi048727d PMID: 15379567

6. Tsivkovskii R, Efremov RG, Lutsenko S. The Role of the Invariant His-1069 in Folding and Function of the Wilson's Disease Protein, the Human Copper-transporting ATPase ATP7B. J Biol Chem. 2003; 278: 13302–13308. doi: 10.1074/jbc.M300034200 PMID: 12551905
7. Abdel-Ghany SE. Contribution of plastocyanin isoforms to photosynthesis and copper homeostasis in Arabidopsis thaliana grown at different copper regimes. Planta. 2009; 229: 767-779. doi: 10.1007/s00425-008-0869-z PMID: 19084994

8. Kieselbach T, Bystedt M, Hynds P, Robinson C, Schröder WP. A peroxidase homologue and novel plastocyanin located by proteomics to the Arabidopsis chloroplast thylakoid lumen. FEBS Letters. 2000; 480: 271-276. doi: 10.1016/S0014-5793(00)01890-1 PMID: 11034343

9. Puig S, Andrés-Colás N, García-Molina A, Peñarrubia L. Copper and iron homeostasis in Arabidopsis: responses to metal deficiencies, interactions and biotechnological applications. Plant, Cell & Environment. 2007; 30: 271-290. doi: 10.1111/j.1365-3040.2007.01642.x PMID: 17263774

10. Catty P, Boutigny S, Miras R, Joyard J, Rolland N, Seigneurin-Berny D. Biochemical characterization of AtHMA6/PAA1, a chloroplast envelope Cu(I)-ATPase. J Biol Chem. 2011; 286: 36188-36197. doi: 10.1074/jbc.M111.241034 PMID: 21878617

11. Shikanai T, Müller-Moulé P, Munekage Y, Niyogi KK, Pilon M. PAA1, a P-type ATPase of Arabidopsis, functions in copper transport in chloroplasts. Plant Cell. 2003; 15: 1333-1346. doi: 10.1105/tpc.011817 PMID: 12782727

12. Seigneurin-Berny D, Gravot A, Auroy P, Mazard C, Kraut A, Finazzi G, et al. HMA1, a new Cu-ATPase of the chloroplast envelope, is essential for growth under adverse light conditions. J Biol Chem. 2006; 281: 2882-2892. doi: 10.1074/jbc.M508333200 PMID: 16282320

13. Tapken W, Ravet K, Pilon M. Plastocyanin Controls the Stabilization of the Thylakoid Cu-transporting P-type ATPase PAA2/HMA8 in Response to Low Copper in Arabidopsis. J Biol Chem. 2012; 287: 18544-18550. doi: 10.1074/jbc.M111.318204 PMID: 22493454

14. Abdel-Ghany SE, Müller-Moulé P, Niyogi KK, Pilon M, Shikanai T. Two P-Type ATPases Are Required for Copper Delivery in Arabidopsis Thaliana Chloroplasts. Plant Cell. 2005; 17: 1233-1251. doi: 10.1105/tpc.104.030452 PMID: 15772282

15. Sautron E, Mayerhofer H, Giustini C, Pro D, Crouzy S, Ravaud S, et al. HMA6 and HMA8 are two chloroplast Cu+-ATPases with different enzymatic properties. Biosci Rep. 2015; doi: 10.1042/BSR20150065 PMID: 26182363

16. Gourdon P, Liu X-Y, Skjørringe T, Morth JP, Møller LB, Pedersen BP, et al. Crystal structure of a copper-transporting PIB-type ATPase. Nature. 2011; 475: 59-64. doi: 10.1038/nature10191 PMID: 21716286

17. Tsuda T, Toyoshima C. Nucleotide recognition by CopA, a Cu+-transporting P-type ATPase. EMBO J. 2009; 28: 1782-1791. doi: 10.1038/emboj.2009.143 PMID: 19478797

18. Wang K, Sitsel O, Meloni G, Autzen HE, Andersson M, Klymchuk T, et al. Structure and mechanism of Zn2+-transporting P-type ATPases. Nature. 2014; advance online publication. doi: 10.1038/nature13618 PMID: 25132545

19. Zhou P, Lugovskoy AA, Wagner G. A solubility-enhancement tag (SET) for NMR studies of poorly behaving proteins. J Biomol NMR. 2001; 20: 11-14. PMID: 11430750

20. Walden H. Selenium incorporation using recombinant techniques. Acta Crystallogr D Biol Crystallogr. 2010; 66: 352-357. doi: 10.1107/S0907444909038207 PMID: 20382987

21. Dimasi N, Flot D, Dupeux F, Márquez JA. Expression, crystallization and X-ray data collection from microcrystals of the extracellular domain of the human inhibitory receptor expressed on myeloid cells IREM-1. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2007; 63: 204-208. doi: 10.1107/S1744309107004903 PMID: 17329815

22. Kabsch W. XDS. Acta Crystallogr D Biol Crystallogr. 2010; 66: 125-132. doi: 10.1107/S0907444909047337 PMID: 20124692

23. Evans PR, Murshudov GN. How good are my data and what is the resolution? Acta Crystallographica Section D Biological Crystallography. 2013; 69: 1204-1214. doi: 10.1107/S0907444913000061 PMID: 23793146

24. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. Phaser crystallographic software. J Appl Crystallogr. 2007; 40: 658-674. doi: 10.1107/S0021889807021206 PMID: 19461840

25. Emsley P, Lohkamp B, Scott WG, Cowtan K. Features and development of Coot. Acta Crystallogr D Biol Crystallogr. 2010; 66: 486-501. doi: 10.1107/S0907444910007493 PMID: 20383002

26. Adams PD, Afonine PV, Bunkóczi G, Chen VB, Davis IW, Echols N, et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr. 2010; 66: 213-221. doi: 10.1107/S0907444909052925 PMID: 20124702

27. Painter J, Merritt EA. TLSMD web server for the generation of multi-group TLS models. Journal of Applied Crystallography. 2006; 39: 109-111. doi: 10.1107/S0021889805038987
28. Sheldrick GM. A short history of SHELX. Acta Crystallographica Section A Foundations of Crystallography. 2007; 64: 112–122. doi: 10.1107/S0108767307043930 PMID: 18156677

29. Cowtan K. The Buccaneer software for automated model building. 1. Tracing protein chains. Acta Crystallogr D Biol Crystallogr. 2006; 62: 1002–1011. doi: 10.1107/S0907444906022116 PMID: 16929101

30. Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, et al. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallographica Section D Biological Crystallography. 2009; 66: 12–21. doi: 10.1107/S0907444909042073 PMID: 20057044

31. Goldschmidt L, Cooper DR, Derewenda ZS, Eisenberg D. Toward rational protein crystallization: A Web server for the design of crystallizable protein variants. Protein Sci. 2007; 16: 1569–1576. doi: 10.1110/ps.072914007 PMID: 17656576

32. Lübben M, Güldenhaupt J, Zoltner M, Deigweiher K, Haebel P, Urbanke C, et al. Sulfate acts as phosphate analog on the monomeric catalytic fragment of the CPx-ATPase CopB from Sulfolobus solfataricus. J Mol Biol. 2007; 369: 368–385. doi: 10.1016/j.jmb.2007.03.029 PMID: 17434529

33. Banci L, Bertini I, Cantini F, Inagaki S, Migliardi M, Rosato A. The binding mode of ATP revealed by the solution structure of the N-domain of human ATP7A. J Biol Chem. 2010; 285: 2537–2544. doi: 10.1074/jbc.M109.054262 PMID: 19917612

34. Dmitriev O, Tsivkovskii R, Abildgaard F, Morgan CT, Markley JL, Lutsenko S. Solution structure of the N-domain of Wilson disease protein: distinct nucleotide-binding environment and effects of disease mutations. Proc Natl Acad Sci USA. 2006; 103: 5302–5307. doi: 10.1073/pnas.0507416103 PMID: 16567646

35. Robert X, Gouet P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014; 42: W320–W324. doi: 10.1093/nar/gku316 PMID: 24753421

36. Hilge M, Siegal G, Vuister GW, Günter P, Gloor SM, Abrahams JP. ATP-induced conformational changes of the nucleotide-binding domain of Na,K-ATPase. Nat Struct Biol. 2003; 10: 468–474. doi: 10.1038/nsb924 PMID: 12730684

37. Abu-Abed M, Millet O, MacLennan DH, Ikura M. Probing nucleotide-binding effects on backbone dynamics and folding of the nucleotide-binding domain of the sarcoplasmic/endoplasmic-reticulum Ca2+-ATPase. Biochem J. 2004; 379: 235–242. doi: 10.1042/BJ20040168 PMID: 14987197

38. Haupt M, Bramkamp M, Heller M, Coles M, Deckers-Hebestreit G, Herkenhoff-Hesselmann B, et al. The holo-form of the nucleotide binding domain of the KdpFABC complex from Escherichia coli reveals a new binding mode. J Biol Chem. 2006; 281: 9641–9649. doi: 10.1074/jbc.M508290200 PMID: 16354672

39. Morgan CT, Tsivkovskii R, Kosinsky YA, Efremov RG, Lutsenko S. The distinct functional properties of the nucleotide-binding domain of ATP7B, the human copper-transporting ATPase: analysis of the Wilson disease mutations E1064A, H1069Q, R1151H, and C1104F. J Biol Chem. 2004; 279: 36363–36371. doi: 10.1074/jbc.M404553200 PMID: 15205462

40. Rodriguez-Granillo A, Sedlak E, Wittung-Stafshede P. Stability and ATP binding of the nucleotide-binding domain of the Wilson disease protein: effect of the common H1069Q mutation. J Mol Biol. 2008; 383: 1097–1111. doi: 10.1016/j.jmb.2008.07.065 PMID: 18692069

41. Liu L, O'Grady C, Dalrymple SA, Prasad L, Dmitriev OY, Delbaere LTJ. Crystallization and preliminary X-ray studies of the N-domain of the Wilson disease associated protein. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2009; 65: 621–624. doi: 10.1107/S1744309109017023 PMID: 19478447

42. Sazinsky MH, Mandal AK, Argüello JM, Rosenzweig AC. Structure of the ATP Binding Domain from the Archaeoglobus fulgidus Cu++-ATPase. J Biol Chem. 2006; 281: 11161–11166. doi: 10.1074/jbc.M510708200 PMID: 16495228

43. Jayakanthan S, Roberts SA, Weichsel A, Arguello JM, McEvoy MM. Conformations of the apo, substrate-bound and phosphate-bound ATP-binding domain of the Cu(II) ATPase CopB illustrate coupling of domain movement to the catalytic cycle. Bioscience reports. 2012; doi: 10.1042/BSR20120048 PMID: 22663904

44. Meng D, Bruschweiler-Li L, Zhang F, Brüschweiler R. NMR backbone resonance assignments of the N, P domains of CopA, a copper-transporting ATPase, in the apo and ligand bound states. Biomol NMR Assign. 2015; 9: 129–133. doi: 10.1007/s12104-014-9558-y PMID: 24706033

45. Stitt M, Lilley RM, Heldt HW. Adenine Nucleotide Levels in the Cytosol, Chloroplasts, and Mitochondria of Wheat Leaf Protoplasts. Plant Physiol. 1982; 70: 971–977. doi: 10.1104/pp.70.4.971 PMID: 16662653