PCuAC domains from methane-oxidizing bacteria use a histidine brace to bind copper

Copper is critically important for methanotrophic bacteria because their primary metabolic enzyme, particulate methane monooxygenase (pMMO), is copper-dependent. In addition to pMMO, many other copper proteins are encoded in the genomes of methanotrophs, including proteins that contain periplasmic copper-A chaperone (PCuAC) domains. Using bioinformatics analyses, we identified three distinct classes of PCuAC domain-containing proteins in methanotrophs, termed PmoF1, PmoF2, and PmoF3. PCuAC domains from other types of bacteria bind a single Cu(I) ion via an HXnMX21/22HXM motif, which is also present in PmoF3, but PmoF1 and PmoF2 lack this motif entirely. Instead, the PCuAC domains of PmoF1 and PmoF2 bind only Cu(II), and PmoF1 binds additional Cu(II) ions in a His-rich extension to its PCuAC domain. Crystal structures of the PmoF1 and PmoF2 PCuAC domains reveal that Cu(II) is coordinated by an N-terminal histidine brace HX10H motif. This binding site is distinct from those of previously characterized PCuAC domains but resembles copper centers in CopC proteins and lytic polysaccharide monooxygenase (LPMO) enzymes. Bioinformatics analysis of the entire PCuAC family reveals previously unappreciated diversity, including sequences that contain both the HXnMX21/22HXM and HX10H motifs, and sequences that lack either set of copper-binding ligands. These findings provide the first characterization of an additional class of copper proteins from methanotrophs, further expand the PCuAC family, and afford new insight into the biological significance of histidine brace-mediated copper coordination.

methanotrophs, bacteria that utilize methane as their sole source of carbon and energy (15), pCu A C genes are frequently contained within the same operon that encodes the copperdependent pMMO enzyme. pMMO, the central metabolic enzyme in methanotrophs, comprises three subunits, PmoC, PmoA, and PmoB, arranged in a trimer of PmoCAB protomers (16 -19). In most ␣-proteobacterial methanotrophs, the pmo operons encoding these three subunits contain up to four additional genes encoding the proteins PmoD, CopC, and CopD and a PCu A C domain (9,20,21). In the methanotroph Methylosinus trichosporium OB3b, these four genes are coregulated with the pmoC, pmoA, and pmoB genes, exhibiting mild up-regulation in response to copper (21). Previous studies have shown that PmoD (21) and CopC (9) are both copper-binding proteins and that PmoD is needed for copper-dependent growth (21). Notably, the periplasmic domain of PmoD has been demonstrated to form an unusual Cu A site in vitro (21,22), but how it is loaded with copper in vivo remains unknown.
In addition to the pCu A C genes in some pmo operons, many methanotrophs also encode other PCu A C proteins that are differentially regulated in response to copper (21,23). A number of methanotroph pCu A C genes do not have a sco1 neighbor so their functions remain unknown and are likely distinct from those of PCu A C proteins that have been studied previously. Also of interest, methanotrophic PCu A Cs encoded within pmo operons frequently have an additional C-terminal extension that includes a cluster of histidine residues.
To further understand the diversity of the PCu A C family, we have performed detailed bioinformatics analyses and have biochemically and structurally characterized two selected methanotroph PCu A C proteins, PmoF1 and PmoF2. These two homologs lack the HX n MX 21/22 HXM Cu(I)-binding motif present in other family members (3,4,6) and in the third methanotrophic PCu A C homolog, PmoF3. Instead, both PmoF1 and PmoF2 bind Cu(II) using a ligand set similar to the histidine brace coordination observed in CopCs (9,24) and lytic polysaccharide monooxygenase (LPMO) enzymes (25,26). Furthermore, PmoF1 contains a His-rich extension that binds additional copper. Bioinformatics analyses indicate that divergent copper-binding sites in PCu A C domains are more prevalent in, but not unique to, methanotrophs. Thus, the PCu A C family is more diverse than appreciated previously, and the different types of copper-binding sites coupled with distinct genomic neighborhoods may correlate with distinct biological functions, including, in the cases of PmoF1 and PmoF2, ones specific to methanotroph copper homeostasis.

Bioinformatics analysis of methanotrophic PCu A C proteins
To identify PCu A C domain-containing proteins in methanotrophs, we searched all methanotroph genomes deposited in the IMG-JGI database for genes encoding proteins containing a PCu A C (pfam04314) domain. PCu A C domains are found primarily in the ␣-proteobacterial methanotroph strains, with only the Methylococcus genus represented among ␥-proteobacteria (Fig. 1A). Many ␣-proteobacterial methanotroph genomes contain up to three PCu A C-encoding genes. The majority of these are found in one of three different genome neighborhoods: within a pmo operon immediately downstream of pmo-CAB, pmoD, copC, and copD (which we denote pmoF1); neighboring a lone pmoD gene (pmoF2); or ϳ9 kb away from the operon encoding components of the soluble methane monooxygenase and the direct neighbor of a sco1 gene (pmoF3) (Fig. 1, B and C). All three pmoF genes are present and exhibit copper-dependent regulation in M. trichosporium OB3b. Specifically, pmoF1 exhibits slight up-regulation in response to copper similar to pmoCAB, and pmoF2 is also up-regulated in response to copper, whereas pmoF3 is markedly down-regulated (21). Both the PmoF1 and PmoF2 homologs are exclu- Figure 1. PCu A C domain-containing proteins in methanotrophs. A, sequence similarity network of all PCu A C proteins in methanotrophs colored by taxonomy. Each colored node represents a single PCu A C domain sequence, and the gray edges connecting the nodes represent the degree of sequence identity between nodes, with shorter edges denoting more closely related proteins. B, same sequence similarity network as in A colored by PmoF type as defined by genome neighborhood. C, typical gene neighborhoods in which pmoF1, pmoF2, and pmoF3 are found. D, domain organization of PmoF1, PmoF2, and PmoF3. The locations of signal peptides (SP), PCu A C domains, and the His-rich region are indicated. E, sequence logos (generated by Skylign (56)) for PmoF1, PmoF2, and PmoF3 homologs in the region of the canonical PCu A C Cu(I)-binding motif. The locations of the canonical HX n MX 21/22 HXM residues are indicated by boxes.
PCu A C domains bind copper using a histidine brace sively found in ␣-proteobacterial strains, suggesting that these proteins may be linked to pMMO function in these organisms, but not in other methanotrophs.
All three PmoF proteins contain signal peptides at their N termini preceding their PCu A C domains (Fig. 1D). In addition, PmoF1 and PmoF2 homologs include C-terminal extensions, whereas PmoF3 homologs do not ( Fig. 1D and Fig. S1A). Previously studied PCu A C proteins contain a Cu(I)-binding HX n MX 21/22 HXM motif (3,4,6,7) and/or have been implicated in copper transfer to Sco1 (3,7,8,10,14). Whereas the PmoF3 proteins share the HX n MX 21/22 HXM motif and also frequently neighbor Sco1 proteins, the conserved metal-binding motif is notably absent in the PmoF1 and PmoF2 proteins ( Fig. 1D and Fig. S1). The C-terminal extension of PmoF1 proteins (ϳ60 residues), which has no predicted secondary structure, includes a cluster of 5-11 histidines within its C-terminal half ( Fig. 1E and Fig. S1). By contrast, the ϳ60 residue C-terminal extensions of PmoF2 proteins lack the concentrated cluster of histidine residues ( Fig. 1E and Fig. S1). Despite the absence of the canonical PCu A C copper-binding motif in PmoF1 and PmoF2, their genomic neighborhoods and previously observed copperdependent regulation patterns suggest that they may also bind copper.

PmoF1 and PmoF2 are Cu(II)-binding proteins
To test for copper binding experimentally, we expressed fulllength PmoF1 (locus tag Met49242_1449, residues 33-202, PmoF1-FL) from Methylocystis sp. Rockwell and full-length PmoF2 (locus tag MettrDRAFT_3935, residues 29 -193, PmoF2-FL) from M. trichosporium OB3b in Escherichia coli as N-terminally SUMO-tagged constructs. These proteins were purified to homogeneity, and the SUMO tags were cleaved using SUMO protease to yield proteins with the authentic N terminus that would be produced upon signal peptide cleavage. Given the preference for Cu(I) exhibited by other members of the PCu A C family, we first investigated Cu(I) binding by anaerobically loading PmoF1-FL and PmoF2-FL with 2-5 molar eq of tetrakis(acetonitrile) Cu(I) hexafluorophosphate. After removing unbound copper using a desalting column, no bound copper could be detected by inductively coupled plasma-optical emission spectroscopy (ICP-OES). Thus, neither PmoF1-FL nor PmoF2-FL binds Cu(I), unlike previously characterized PCu A C domain-containing proteins that all reportedly bind 1 eq of Cu(I) (3)(4)(5).
We then assessed Cu(II) binding in a similar manner, adding 0 -10 Cu(II) equivalents. PmoF1-FL binds 2.2 Ϯ 0.6 molar eq of Cu(II) ( Fig. 2A). To localize the copper-binding sites, we generated an additional construct lacking the C-terminal extension (residues 33-153, PmoF1-⌬C). The addition of 0 -5 Cu(II) equivalents indicates PmoF1-⌬C binds 0.75 Ϯ 0.3 Cu(II) ions ( Fig. 2A). Together, these data suggest that the His-rich extension, which in Methylocystis sp. Rockwell PmoF1 includes 7 histidine residues (Fig. S2), contains additional copper-binding sites. The variability in copper-binding stoichiometry of PmoF1-FL may be related to changing accessibility of the presumably unstructured C-terminal extension. Other members of the PCu A C family, including the AccA protein from Neisseria species (12) and PcuC from Bradyrhizobium diazoefficiens (14), also contain C-terminal His-or Met-rich extensions. In the case of B. diazoefficiens PcuC, this extension delivers copper to Sco1 and to cytochrome c oxidase (14). In other types of proteins, such extensions have also been implicated in copper binding. For example, some Cu/Zn superoxide dismutases include an N-terminal His-rich region that binds copper and has been proposed to either directly insert the metal into a neighboring protein molecule or increase the local copper concentration near the active site (27). Similarly, the N-terminal metal-binding domain of the Cu(I)-specific CopB P 1B -ATPase contains a large number of histidine residues and binds up to eight additional copper ions (28). It is possible that the His-rich extension in PmoF1-FL binds periplasmic copper and delivers it to the PCu A C domain or to other protein targets. Alternatively, it may serve a regulatory function.
In contrast to PmoF1-FL, the addition of excess Cu(II) to PmoF2-FL results in 0.90 Ϯ 0.4 Cu(II) ions bound per protein molecule (Fig. 2B). PmoF2 homologs also contain a C-terminal extension to the PCu A C domain. Many of these extensions include two histidine residues in the equivalent part of the sequence to the PmoF1 His-rich region (Fig. S1D). M. trichosporium OB3b PmoF2, however, does not have these two histidine residues, but does have three other histidine residues within its C-terminal region, His-152, His-154, and His-175 (Fig. S2). The PmoF2 PCu A C domain alone (residues 29 -159, PmoF2-⌬C), which includes His-152 and His-154, but not His-175, binds as many Cu(II) equivalents as PmoF1-FL, 0.94 Ϯ 0.2 Cu(II) per domain (Fig. 2B), indicating that unlike PmoF1,

PCu A C domains bind copper using a histidine brace
PmoF2-FL does not bind significantly more copper than the PCu A C domain alone.

Crystal structures reveal that PmoF proteins bind copper using a histidine brace
Given the absence of the HX n MX 21/22 HXM motif in both PmoF1 and PmoF2 and the unusual preference for Cu(II) rather than Cu(I), we pursued structural studies of both apo-and copper-loaded PCu A C domains from each protein. We obtained crystals of the apo-SM PmoF1-⌬C (a construct identical to PmoF1-⌬C but with two additional N-terminal vector-derived residues, a serine and a methionine) and PmoF2-⌬C, and solved the structures to 1.85 and 1.65 Å resolutions, respectively (Table 1). Both structures comprise primarily two ␤-sheets in a Greek key topology (Fig. 3, A and B), similar to the overall fold of other structurally characterized members of the PCu A C family. Superpositions of the SM PmoF1-⌬C and PmoF2-⌬C structures with that of the ECuC protein from S. lividans yield r.m.s.d. values of 1.42 Å over 106 C␣ atoms and 1.82 Å over 108 C␣ atoms, respectively (Fig. 4A). The SM PmoF1-⌬C and PmoF2-⌬C structures are also highly similar to one another with an r.m.s.d. of 1.70 Å over 119 C␣ atoms (Fig. 3). One notable divergence between previously characterized PCu A C domains and the current structures is the presence of an additional ␤ strand at the N termini of PmoF1 and PmoF2 (termed ␤0) that extends the ␤-sheet formed by strands ␤1, ␤3, and ␤5 (Fig. 4). We also crystallized Cu(II)-loaded forms of each protein and determined the structures of PmoF1-⌬C and SM PmoF2-⌬C to 1.91 and 2.05 Å resolutions, respectively ( Table 1). The overall fold is similar in both cases (Fig. 3), with r.m.s.d. values of 0.63 and 1.36 Å for superposition of 120 C␣ atoms from the apo-and copper-loaded forms of PmoF1 and PmoF2, respectively. Thus, copper binding does not induce significant conformational changes.
The copper-loaded PmoF1-⌬C structure contains a single copper ion in square pyramidal geometry (Fig. 3A). This copper ion is coordinated by the ␦-nitrogen atoms of two histidine residues at the N terminus of the protein (bond lengths ϳ2 Å). These two histidine ligands derive from a HX 10 H motif. The N terminus itself provides a third nitrogen ligand in the same plane, forming a T-shaped histidine brace coordination motif ( Fig. 3A and Table 2). An axial water molecule and a glutamate side chain oxygen from the second molecule in the asymmetric unit are located 2.7/3.1 and 2.2 Å from the copper ion, respectively. Additionally, size-exclusion chromatography with multiangle light scattering (SEC-MALS) analysis of the copperloaded PmoF1-⌬C revealed a mix of species with an average molar mass of 17,930 Da Ϯ 6%, which is between the expected molar masses of a monomer and a dimer (Fig. S3A). However, SEC-MALS analysis of copper-loaded PmoF1-FL indicates a single species with the expected mass of a monomer (Fig. S3B). Because the glutamate coordination observed in the crystal structure can only occur in a dimeric species, this coordination most likely represents a crystallographic artifact. Copper coordination by a glutamate from a neighboring molecule in the crystal lattice is also observed in one of the LPMO structures (29).
The copper-loaded SM PmoF2-⌬C structure similarly reveals a single copper ion coordinated by an HX 10 H motif (Fig. 3B). Although the construct used to obtain this structure included two additional vector-derived residues upstream of the authen-

PCu A C domains bind copper using a histidine brace
tic N-terminal histidine (we did not obtain crystals of the PmoF2-⌬C construct), the copper coordination, similar to that in the PmoF1-⌬C structure, includes the ␦-nitrogen atoms from two histidine side chains and the N terminus of the peptide. The peptide backbone from the vector-derived residues provides two additional deprotonated amide ligands, oriented similarly to the water molecule and glutamate side chain observed in the PmoF1-⌬C structure (Fig. 3B). Consistent with the crystallographic observations, copper-loaded PmoF2-⌬C (as well as copper-loaded PmoF2-FL) is a single, monomeric species by SEC-MALS (Fig. S3, C and D). For both PmoF1-⌬C and SM PmoF2-⌬C, there is a large conformational change in the side chains of the two coordinating histidine residues between the apo-and copper-loaded forms (Fig. 3).
Despite the similar tertiary structure (Fig. 4A), the copperbinding sites observed in the PmoF1-⌬C and SM PmoF2-⌬C structures are completely different from those in the other structurally characterized PCu A C domains. In the other PCu A C domains, the copper ion is coordinated by one histidine residue and one methionine residue from strand ␤5 (e.g. His-127 and Met-129 in ECuC), a second histidine from the ␤4 strand antiparallel to it (e.g. His-98 in ECuC), and a second methionine in a short extension of ␤4 termed ␤4Ј (e.g. Met-105 in ECuC). Not only do PmoF1 and PmoF2 lack both methionine ligands (Methylocystis sp. Rockwell PmoF1 has an arginine and an alanine, whereas M. trichosporium OB3b PmoF2 has a leucine and a glutamate), but in both structures, ␤4 is significantly shortened, and ␤4Ј is absent (Fig. 4B). PmoF1 also has a tyrosine residue in the equivalent position to ECuC His-98. Likewise, the N-terminal regions of the proteins diverge significantly. In other previously characterized PCu A C domains, neither of the coordinating histidine residues found in PmoF1 and PmoF2 is conserved, and the ␤0 strand is absent, precluding formation of a histidine brace copper coordination site (Fig. 4).

PCu A C domains bind copper using a histidine brace
In the copper-loaded SM PmoF2-⌬C structure, we also observed strong electron density that could not be modeled as a water molecule in a similar position to the Cu(I) site observed in other PCu A C proteins (7) (Fig. S4, A and B). However, anomalous difference Fourier maps calculated using data collected at the copper absorption edge do not exhibit peaks at this position, indicating that it is not a copper ion. Due to the presence of low levels of zinc in this particular sample (0.2 eq per protein molecule as determined by ICP-OES), we have modeled a Zn(II) ion in this position, coordinated by His-111, His-83, and Glu-113 (Fig. S4C). In support of this model, Zn(II)-binding experiments indicate that PmoF2-⌬C and PmoF2-FL bind 0.78 Ϯ 0.4 and 1.16 Ϯ 0.1 eq, respectively (Fig. S5). In one of the two molecules in the asymmetric unit, His-87 provides an additional nitrogen ligand, but it is disordered in the second molecule (Fig. S4D). Whereas His-83 is completely conserved, His-87 is 33% conserved, and both His-111 and Glu-113 are only 22% conserved among PmoF2 sequences. Because some samples of PmoF2-FL contained Ͼ1 eq of Cu(II) (Fig. 2B), we cannot rule out the possibility that there is a second metalbinding site in PmoF2-FL (presumably featuring histidine ligation), although the lack of conservation of these additional histidine residues suggests that it is an artifact.

Electron paramagnetic resonance (EPR) characterization of the PmoF1 and PmoF2 Cu(II) sites
We further characterized the Cu(II) sites of PmoF1 and PmoF2 by X-band continuous-wave (CW) EPR spectroscopy. PmoF1-⌬C loaded with 1.0 Cu(II) equivalent exhibits a mononuclear Cu(II) EPR spectrum (g ϭ [2.25, 2.07, 2.02], 63 Cu A 1 ϭ 560 MHz or 187 ϫ 10 Ϫ4 cm Ϫ1 ) (Fig. 5A). This Cu(II) signal is assigned to the only copper site observed in the PmoF1-⌬C crystal structure, the square pyramidal N-terminal histidine brace site (Fig. 3A). Consistent with this assignment, the ratio of g 1 /A 1 , which corresponds to the tetrahedral distortion of a Cu(II) center (in which a larger number indicates greater distortion), is 120 cm for the PmoF1-⌬C Cu(II) signal, a value indicative of square planar equatorial ligation (30,31). Indeed, all PmoF1 and PmoF2 Cu(II) species reported herein fall within the range for square planar equatorially ligated Cu(II) sites (31), g 1 /A 1 ϭ 105-135 cm (except for one Cu(II) species with g 1 /A 1 ϭ 136 cm).
PmoF1-FL loaded with 0.9 eq of Cu(II) exhibits an EPR spectrum (g ϭ [2.24, 2.07, 2.02], 63 Cu A 1 ϭ 535 MHz or 178 ϫ 10 Ϫ4 cm Ϫ1 , g 1 /A 1 ϭ 126 cm) resembling that of PmoF1-⌬C (Fig. 5A). One notable difference between the two spectra is the wellresolved 14 N hyperfine splitting along g 2 in the PmoF1-FL, but not PmoF1-⌬C, spectrum. Unlike PmoF1-⌬C, PmoF1-FL contains the ligands for both the N-terminal histidine brace Cu(II) site and a His-rich region ( Fig. 1D and Fig. S1D) which binds Cu(II) (Fig. 2A). The simplest interpretation of these results is that in PmoF1-FL, the His-rich region binds Cu(II) before the histidine brace site and binds with a similar ligation and geometry to the histidine brace site, resulting in a similar but not identical EPR spectrum. Consequently, the EPR spectrum of PmoF1-FL loaded with 0.9 eq of Cu(II) corresponds to Cu(II) bound in the His-rich region.
The addition of 3 eq of Cu(II) to PmoF1-FL results in additional Cu(II) binding ( Fig. 2A) but minimal changes to the EPR spectrum (g ϭ [2.24, 2.07, 2.02], 63 Cu A 1 ϭ 545 MHz or 182 ϫ 10 Ϫ4 cm Ϫ1 , g 1 /A 1 ϭ 123 cm) (Fig. 5A), the two most evident being 1) the broadened lowest-field Cu A 1 hyperfine line and 2) decreased resolution of the 14 N hyperfine splitting along g 2 . These results are consistent with the superstoichiometrically Cu(II)-loaded PmoF1-FL spectra being the sum of overlapping signals observed for PmoF1-⌬C and the signal observed for PmoF1-FL loaded with 0.9 eq of Cu(II). Therefore, in the superstoichiometrically Cu(II)-loaded PmoF1-FL sample, both the N-terminal histidine brace site and the His-rich region have Cu(II) bound. These combined results indicate that the additional Cu(II) sites in PmoF1-FL involve ligation and geometry similar, but not identical, to that of the N-terminal histidine brace Cu(II) site.
PmoF2-⌬C loaded with 1.0 eq of Cu(II) exhibits a Cu(II) EPR spectrum (g ϭ [2.25, 2.07, 2.03], 63 Cu A 1 ϭ 525 MHz or 175 ϫ 10 Ϫ4 cm Ϫ1 , g 1 /A 1 ϭ 129 cm) with well-resolved 14 N hyperfine splitting along g 2 (Fig. 5B). This Cu(II) signal is assigned to the N-terminal histidine brace Cu(II)-binding site, which is the only Cu(II) site observed in the crystal structure (Fig. 3B). The  (Fig. 5B). This suggests that both PmoF2-⌬C and PmoF2-FL bind Cu(II) in the N-terminal histidine brace Cu(II) site, but in PmoF2-FL, the site is not as rigidly structured as in PmoF2-⌬C, and therefore the 14 N hyperfine splitting along g 2 is not as wellresolved in PmoF2-FL as it is in PmoF2-⌬C.

Bioinformatics analyses of the entire PCu A C family
To assess the prevalence of the PmoF type of copper-binding site within the larger PCu A C family, we extended our bioinformatics analysis to all bacterial genomes. We constructed a sequence similarity network of all PCu A C domain-containing proteins, 22,649 proteins in total (Fig. 6). The vast majority of PCu A C proteins are encoded by bacterial genomes, particularly those of proteobacteria and actinobacteria. Although the sequence similarity network provides information regarding .0] MHz, anisotropic g-strain (EasySpin g-strain parameters 0.03, 0.03, and 0 for g 3 , g 2 , and g 1 , respectively) with additional 60 G Gaussian isotropic line-width broadening.

PCu A C domains bind copper using a histidine brace
global sequence similarity, it does not necessarily correspond to conservation of the metal-binding motif(s). Therefore, we aligned all of the bacterial sequences to the hidden Markov model (HMM) for the PCu A C family and analyzed each sequence for the presence of either the canonical HX n MX 21/22 HXM motif, the N-terminal HX 10 H-binding motif identified in this study, or both. Whereas the canonical binding motif is present in 93.0% of the total sequences, it is absent in the remainder, including 0.5% of the total PCu A C domain-containing proteins, such as PmoF1 and PmoF2, that have the N-terminal HX 10 H motif (122 sequences in total, 22 of which are from methanotroph genomes) (Fig.  S6A). Similar to the majority of the methanotroph genomes represented among this group, the majority of the other species represented are also ␣-proteobacteria. More functional data will be required to ascertain whether these proteins represent a functionally distinct group. Of the total sequences (2836 sequences), 12.5% contain both the canonical Cu(I)-binding motif and the histidine brace Cu(II)-bind-ing motif, whereas 6.5% (1468 sequences) lack both Cu-binding motifs. The latter are more frequently found in actinobacterial genomes (Fig. S6B). Interestingly, a protein from Neisseria meningitidis that has both sites has been reported to bind two copper ions, although the specific coordinating residues have not been identified experimentally (12).

Discussion
Previously characterized PCu A C domains bind Cu(I) with two histidines and two methionines from an HX n MX 21/22 HXM motif (3,4,6), whereas PmoF1 and PmoF2 bind Cu(II) with a histidine brace provided by an HX 10 H motif. Histidine brace coordination of copper is fairly rare but has been observed in several other proteins. Most relevant to PmoF1 and PmoF2 (Fig. 7, A and B) is the nearly identical coordination observed for the Cu(II) site in the periplasmic copper-binding protein CopC, which contains an aspartate in or near the fourth equatorial ligation position (9, 24, 32) (Fig. 7C and Table 2). The Figure 6. Sequence similarity network for the PCu A C family. Shown is a sequence similarity network for all bacterial PCu A C domain-containing proteins, colored by putative copper-binding motif. Each colored node represents a set of sequences that share Ն90% sequence identity, with the length of the gray edges connecting the nodes indicative of the relative sequence similarity between them. The clusters that contain a sequence corresponding to proteins that have been structurally characterized are indicated by colored boxes. PCu A C domains bind copper using a histidine brace majority of CopCs, denoted C 0-1 CopCs, have only this site, although one subset, the C 1-1 CopCs, also contain a distinct Cu(I) site coordinated by primarily methionine residues (9). The C 0-1 CopCs are usually found in genomic proximity or as fusions to CopDs, which are proposed to import copper, and CopC has been proposed to deliver copper to CopD, a function consistent with the solvent-exposed location of its binding site (20,33,34).
Similarly, PmoF1 and PmoF2 may play a role in copper delivery to various target proteins. The gene encoding PmoF1 is part of the pmo operon (Fig. 1C), and PmoF1 could interact with pMMO itself, PmoD, CopC, or CopD. Interestingly, the gene encoding PmoF2 is not found in a pmo operon but is paired with a pmoD gene (Fig. 1C). PmoD proteins are only found in methane-and ammonia-oxidizing bacteria (21), an observation that, combined with the genomic proximity to PmoF-encoding genes and the prevalence of the HX 10 H motif-containing PCu A C homologs in ␣-proteobacterial methanotrophs, suggests a function specific to methanotroph copper homeostasis. One possibility is that PmoF supplies PmoD with copper, whereas the canonical HX n MX 21/22 HXM motif-containing PCu A C proteins function in the assembly of the cytochrome c oxidase copper centers. Studies to address this hypothesis are under way.
Besides expanding the diversity of the PCu A C family, characterization of PmoF1 and PmoF2 has implications for the role of histidine brace-ligated copper sites in biology. Whereas PmoFs and CopCs likely use the solvent-exposed histidine brace site for copper transfer, similar sites are found in some copper enzymes. In particular, the defining feature of the copper active site of LPMOs is histidine brace coordination (26,35,36) (Fig. 7D). These enzymes hydroxylate C-H bonds of glycosidic linkages, which then are cleaved by elimination and are also referred to as PMOs, as the lytic elimination reaction is not actually enzymatic (37). The histidine brace coordination geometry has been proposed to facilitate the formation of highvalent copper-oxygen intermediates for oxidation of an unactivated C-H bond in the substrate (26). Although the same effect has been proposed for the Cu B site in pMMO (26), which has similar coordination, albeit with a third histidine (18,19,38), recent work suggests that methane oxidation does not occur at this site (38).
Given that the CopC and PmoF proteins are not involved in cleavage of strong C-H bonds, the main advantage of their histidine brace may be to provide a versatile solvent-exposed site for copper transfer or other functions. Importantly, a solventexposed site is also absolutely necessary for LPMOs, given that their substrate is lignocellulosic biomass (39). This requirement, rather than an ability to confer special oxidative power, may explain the presence of a histidine brace in LPMOs. Additionally, variations in the coordination, such as ligation by the ⑀-nitrogen atom of the second histidine in LPMOs rather than the ␦-nitrogen atom in PmoFs and CopCs, could affect reactivity. Steric, electronic, and chemical (pK a ) (40) differences may govern whether a histidine brace copper site is catalytic; ligation by the ⑀ nitrogen provides more flexibility for oxygen and substrate binding (41). These considerations are important to understanding the linkage between histidine brace coordination geometry and biological function.

Construct design and cloning
The full-length PmoF1 ( SM PmoF1-FL), excluding the signal peptide and comprising residues 33-202 from Methylocystis sp. Rockwell (locus tag Met49242_1449) was subcloned into the pSGC-His vector with an N-terminal TEV-cleavable His 6 tag, a gift from Anthony Gizzi and Steven Almo (Albert Einstein College of Medicine), using primer set 1 (Table S1). Due to the presence of two vector-derived residues (Ser-Met) at the N terminus in SM PmoF1-FL, we also generated a construct (PmoF1-FL) in which we inserted a SUMO tag via Gibson assembly using primer sets 5 and 6 immediately N-terminal to His-33 to generate the authentic N terminus by cleavage with SUMO protease. To produce constructs containing only the PCu A C domain (residues 33-153), residue 154 was mutated to a stop codon using primer set 2 and the QuikChange XL mutagenesis kit (Agilent) to generate His 6 -tagged constructs SM PmoF1-⌬C and PmoF1-⌬C.
The DNA for PmoF2 from M. trichosporium OB3b (locus tag MettrDRAFT_3935), encoding residues 29 -193, was synthesized by Genscript with an EcoRI restriction site at the 3Ј end and a BamHI restriction site at the 5Ј end. This sequence was subcloned into the pPR-IBA2 vector (IBA Lifesciences) using EcoRI and BamHI restriction sites to produce Strep-tagged fulllength PmoF2 ( SM PmoF2-FL). A SUMO tag was then inserted immediately upstream of His-29 by Gibson assembly using primer sets 5 and 6 (Table S1) to generate the PmoF2-FL construct. To produce PmoF2 constructs containing only the PCu A C domain (residues 29 -159), residue 160 was mutated to a stop codon using primer set 7 (Table S1) with the QuikChange XL mutagenesis kit (Agilent) to generate Strep-tagged constructs SM PmoF2-⌬C and PmoF2-⌬C.

Protein expression and purification
All PmoF1 and PmoF2 constructs described above were transformed and expressed following the protocol described below. Each plasmid was transformed into BL21* (DE3) cells (Novagen). Overnight cultures were grown in lysogeny broth medium, inoculated into autoinduction medium (42), and grown with shaking at 180 rpm at 37°C. When cultures reached

PCu A C domains bind copper using a histidine brace
A 600 ϭ 0.6, the temperature was lowered to 22°C overnight. Cells were harvested by centrifugation at 6000 ϫ g for 20 min and resuspended in lysis buffer (25 mM PIPES, pH 7.0, 250 mM NaCl) supplemented with 1 mM DTT, 1 mg/ml DNase I, and 1 mM phenylmethylsulfonyl fluoride. Cells were lysed by sonication for 20 min with 1-s pulses, and the cell debris was removed by centrifugation at 24,000 ϫ g for 1 h.
For the His-tagged PmoF1 proteins, the clarified lysate was loaded onto Ni-NTA resin. The column was then washed with 25 mM PIPES, pH 7.0, 250 mM NaCl, 20 mM imidazole, and proteins were eluted with 25 mM PIPES, pH 7.0, 250 mM NaCl, 500 mM imidazole. Proteins were then incubated with His 6tagged TEV (in the case of the SM PmoF1-⌬C construct) or His 6tagged SUMO (in the case of the PmoF1-FL and PmoF1-⌬C constructs) proteases overnight at 4°C with nutation (43,44). Cleavage products of the PmoF1-⌬C construct (lacking the His-rich C-terminal extension) were separated by loading the reaction onto Ni-NTA resin; washing with 2 column volumes of buffer containing 25 mM PIPES, pH 7.0, 250 mM NaCl, 20 mM imidazole; and collecting the flow-through. For the full-length PmoF1 constructs, cleavage products were dialyzed into saltfree buffer (25 mM PIPES, pH 7.0) using 10-kDa MWCO Snake-Skin dialysis tubing. They were then separated using a HiTrap Q FF anion exchange column (GE Healthcare Life Sciences), from which PmoF1 eluted in the flow-through.
For the Strep-tagged PmoF2 proteins, clarified lysate was applied to a StrepTactin column and washed with 5-10 column volumes of 25 mM PIPES, pH 7.0, 250 mM NaCl. The Streptagged protein was then eluted in 25 mM PIPES, pH 7.0, 250 mM NaCl, 2.5 mM D-desthiobiotin. Proteins were then incubated with His-tagged TEV (in the case of SM PmoF1-⌬C) or SUMO (in the case of PmoF2-FL and PmoF2-⌬C) proteases overnight at 4°C with nutation (43,44). The proteins were then loaded onto Ni-NTA resin to remove the His-tagged protease, and the flow-through containing the cleavage products was reapplied to the StrepTactin column to remove the Strep-tagged SUMO. The flow-through from this column, containing the untagged PmoF2 protein, was concentrated by centrifugation to ϳ1-2 ml in 10-kDa MWCO concentrators (Millipore).
The untagged PmoF1 and PmoF2 proteins were each subjected to a final purification step on a HiLoad 16/600 Superdex 75 column (GE Healthcare). Peak fractions were pooled and concentrated by centrifugation. The protein concentration was measured by absorbance at A 280 using a Nanodrop. The following extinction coefficients were used as determined using the ExPASy ProtParam tool: PmoF1-⌬C, 15 (Fig. S7).

Copper-and zinc-loading experiments
To test for Cu(I) binding, a 10 mM tetrakis(acetonitrile)Cu(I) hexafluorophosphate solution was prepared in 100% acetonitrile in a Coy anaerobic chamber. The buffer solution was deoxygenated on a Schlenk line by completion of three purge cycles with high-purity argon gas. After cycling, the buffer was left under argon for 30 -45 min before it was taken into the anaerobic chamber. To ensure removal of all dissolved oxygen, the buffer was opened and stirred overnight in the chamber. The buffer was supplemented with 5 mM 2-mercaptoethanol prior to copper-loading experiments. Then 1-5 molar eq of the Cu(I) solution were added to 20 M protein in deoxygenated buffer in the anaerobic chamber. The solution was incubated at room temperature for 1-2 h and desalted on a PD-10 column equilibrated with deoxygenated buffer containing 2-mercaptoethanol. To test for Cu(II) binding, 1-10 molar eq of CuSO 4 were slowly added to 10 -50 M protein in 25 mM PIPES, pH 7.0, 250 mM NaCl. This solution was incubated on ice for 1-2 h prior to removing unbound metal using a PD-10 desalting column (GE Healthcare). Proteins were eluted in 3.5 ml of 25 mM PIPES, pH 7.0, 250 mM NaCl and subsequently concentrated to ϳ500 l using 10-kDa MWCO concentrators (Millipore). Concentrated and desalted proteins loaded with 0 -10 eq of Cu(II) were used to prepare samples for ICP-OES.
To test for Zn(II) binding to PmoF2 constructs, 1 and 5 molar eq of ZnSO 4 ⅐7H 2 O were slowly added to 20 M protein in 25 mM PIPES, pH 7.0, 250 mM NaCl. This solution was incubated on ice for 2-3 h prior to removing unbound metal using a PD-10 desalting column. Proteins were eluted and concentrated as described above, and concentrated and desalted proteins loaded with 0, 1, and 5 eq of Zn(II) were used to prepare samples for ICP-OES.

ICP-OES
To determine the copper content of all proteins and zinc content of PmoF2-⌬C and PmoF2-FL, the metal-loaded, desalted samples were digested in a 3% nitric acid solution in metal-free tubes (VWR). A standard curve was generated from a serial dilution of a custom multielement standard (Inorganic Ventures) that was also digested in 3% nitric acid. The metal contents of the standards and experimental samples were determined using a Thermo iCAP 7600 ICP-OES instrument (Quantitative Bio-element Imaging Center core facility at Northwestern University). Metal concentrations and stoichiometry were calculated based on the standard curve.

EPR spectroscopy
Prior to EPR measurements, each sample was Cu(II)-loaded and desalted as described above and then concentrated using 10-kDa MWCO concentrators (Millipore) via centrifugation. ϳ180 l of protein solution (protein concentrations ranging from 200 to 450 M in buffer containing 25 mM PIPES, pH 7.0, and 250 mM NaCl) was aliquoted into a Wilmad quartz X-band EPR tube (Sigma-Aldrich). This tube was then flash-frozen in liquid nitrogen; the sample was stored in liquid nitrogen until EPR measurements were conducted. CW X-band EPR measurements were conducted at 77 K utilizing a modified Varian E-4 spectrometer, equipped with a finger-Dewar filled with liquid nitrogen. Simulation of EPR spectra was performed using EasySpin (45).

Crystallization and structure determination
An initial crystal structure was obtained using copper-loaded SM PmoF1-⌬C. The protein crystallized at 3.5 mg/ml in 1-l PCu A C domains bind copper using a histidine brace drops at a 1:1 ratio with a precipitant composed of 0.2 M potassium thiocyanate and 20% PEG 3350 (PEGs Suite, Qiagen) and was cryo-cooled in liquid nitrogen prior to data collection. A data set was collected at 1.362 Å on this crystal at beamline 21-ID-D (Life Sciences Collaborative Access Team (LS-CAT), Advanced Photon Source). The data were processed using HKL2000 (46) to 1.6 Å resolution in space group P2 1 2 1 2 1 . The structure was solved using the program CRANK2 (47) as implemented in CCP4i2 (48). This program identified two Cu sites and built an initial model comprising 235 residues with R work / R free ϭ 26.0%/29.2%. The structure was further refined using iterative rounds of refinement in phenix.refine (49) and modelbuilding in Coot (50). The final model has R work /R free ϭ 17.3%/ 19.2%. In this structure, two copper ions are coordinated at a crystallographic interface between two protein chains with an additional ligand provided by what we have modeled as a thiocyanate ion from the crystallization condition ( Fig. S8 and Table S2). This crystal structure, although not biologically relevant due to the thiocyanate-mediated dimer, was used as a search model to solve all subsequent crystal structures by molecular replacement.
Apo-SM PmoF1-⌬C initially crystallized at 5 mg/ml in 0.1 M sodium acetate, pH 4.6, 0.2 M ammonium acetate, and 30% PEG 4000 (Classics Suite I, Qiagen) using a 96-well screen and 1-l drops in a 1:1 ratio of protein/precipitant. Conditions were optimized to 0.1 M sodium acetate, pH 4.6, 0.2 M ammonium acetate, and 34% PEG 4000. A data set was collected at LS-CAT beamline 21-ID-D and was processed to 1.85 Å resolution using XDS (51) in space group P2 1 2 1 2 1 . An initial solution was obtained using molecular replacement in phenix.phaser (49) with chain A of the copper-bound SM PmoF1-⌬C structure as a search model.
For the apo-PmoF2-⌬C structure, an initial crystallization hit was obtained in 0.1 M sodium acetate, pH 4.6, and 2.0 M ammonium sulfate in the Classics Suite (Qiagen) at 10 mg/ml protein concentration using a 96-well screen and 1-l drops in a 1:1 ratio of protein/precipitant. These conditions were further optimized to obtain well-diffracting single crystals in a precipitant solution of 0.1 M sodium acetate, pH 4.6, and 1.727 M ammonium sulfate using protein at 2.5 mg/ml. The crystal used for data collection was briefly transferred to a 3.4 M sodium malonate solution for cryoprotection prior to cryocooling in liquid nitrogen. Data collected at LS-CAT beamline 21-ID-D were processed using XDS (51) in space group P4 1 2 1 2. The structure was solved using phenix.phaser (49) with chain A of the Cu-bound SM PmoF1-⌬C structure as a search model, followed by autobuilding resulting in a model including 102 residues and R work /R free of 24.7%/27.5%. Further rounds of refine-ment and model building were performed using phenix.refine (49) and Coot (50), respectively.
For the copper-bound structure of SM PmoF2-⌬C, the protein crystallized in 0.2 M magnesium acetate, 0.1 M sodium cacodylate, pH 6.5, and 20% PEG 8000, in 1-l drops with a 1:1 protein/precipitant ratio and protein concentration of 10 mg/ml. A single crystal was cryocooled in liquid nitrogen prior to data collection at LS-CAT beamline 21-ID-D. A data set was collected and processed to 2.05 Å resolution in space group P3 2 21 using XDS (51). A partial model was obtained using phenix.phaser (49) with the same search model as described above, which was then input into autosol to obtain an initial solution by MR-SAD that resulted in building 234 residues with R work / R free of 22.0%/26.4%. Further rounds of refinement and model building were performed using phenix.refine (49) and Coot (50), respectively.

SEC-MALS
PmoF1-⌬C and PmoF2-⌬C were prepared for analysis by incubating each protein with 1 molar eq CuSO 4 on ice for 1-2 h prior to desalting on a PD10 column. The desalted, copperloaded proteins were then concentrated to 3.5 mg/ml in 10-kDa MWCO concentrators. The protein samples were filtered, and 300 l was applied to a Superdex 75 10/300 GL column (GE Healthcare) column in 250 mM NaCl, 25 mM PIPES, pH 7.0, at room temperature. Analysis was performed using an Agilent 1260 series HPLC system equipped with diode-array detection absorbance in-line with a DAWN HELEOS II multiangle static light-scattering detector (Wyatt Technology), a QELS (quasielastic light scattering) dynamic light-scattering detector (Wyatt Technology), and a T-rEx differential refractive index detector (Wyatt Technology) Each sample was run at 0.4 ml/min for 60 min, and data were processed using Astra software version 5.3.4 (Wyatt Technology).

Bioinformatics analyses
Methanotroph genomes were searched for genes encoding PCu A C domain-containing proteins. These sequences and their associated data were identified in the IMG-JGI database by searching all methanotroph genuses (as identified in Ref. 52) for pfam04314 (File S1). The sequences were then used to generate a sequence similarity network by submitting the FASTA file generated above (Files S2 and S3) to EFI-EST using an alignment score of 40 and an E-value of 5 (53). The resulting networks were visualized using Cytoscape (54). To identify metalbinding sites, the sequences were aligned to the PCu A C HMM using Clustal Omega (55) (File S4).
To generate the sequence similarity network for the full complement of PCu A C domain-containing proteins, sequences were identified in IMG-JGI by searching all bacterial genomes using the HMM for pfam04314, yielding a total of 22,649 sequences in total (File S5). A sequence similarity network was generated using the user-supplied FASTA file generated in this manner using EFI-EST as described above for the methanotroph sequences. Additional parameters included using an alignment score cut-off of 30, an E-value of 5, restricting the sequence length to 100 -350 amino acids, and representing sequences with Ͼ90% sequence identity by a single node (Files PCu A C domains bind copper using a histidine brace