Conformational selection underlies recognition of a
molybdoenzyme by its dedicated chaperone

M. Lorenzi, L. Sylvi, G. Gerbaud, E. Mileo, F. Halgand, A. Walburger, H. Vezin, V. Belle, B. Guigliarelli, A. Magalon

► To cite this version:

M. Lorenzi, L. Sylvi, G. Gerbaud, E. Mileo, F. Halgand, et al. Conformational selection underlies recognition of a molybdoenzyme by its dedicated chaperone. PLoS ONE, 2012, 7 (11), pp.e49523. 10.1371/journal.pone.0049523. hal-00758737

HAL Id: hal-00758737

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

Submitted on 27 Sep 2018

**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
# Conformational Selection Underlies Recognition of a Molybdoenzyme by Its Dedicated Chaperone

Magali Lorenzi¹, Léa Sylvi², Guillaume Gerbaud¹, Elisabetta Mileo¹, Frédéric Halgand¹*,
Anne Walburger², Hervé Vezin³, Valérie Belle¹*, Bruno Guigliarelli¹, Axel Magalon²*

1 Unité de Bioénergétique et Ingénierie des Protéines (UMR7281), Institut de Microbiologie de la Méditerranée, CNRS & Aix-Marseille Univ, Marseille, France, 2 Laboratoire de Chimie Bactérienne (UMR7283), Institut de Microbiologie de la Méditerranée, CNRS & Aix-Marseille Univ, Marseille, France, 3 Laboratoire de Spectrochimie Infrarouge et Raman (UMR8516), Villeneuve d'Ascq, France

## Abstract

Molecular recognition is central to all biological processes. Understanding the key role played by dedicated chaperones in metalloprotein folding and assembly requires the knowledge of their conformational ensembles. In this study, the NarJ chaperone dedicated to the assembly of the membrane-bound respiratory nitrate reductase complex NarGHI, a molybdenum-iron containing metalloprotein, was taken as a model of dedicated chaperone. The combination of two techniques ie site-directed spin labeling followed by EPR spectroscopy and ion mobility mass spectrometry, was used to get information about the structure and conformational dynamics of the NarJ chaperone upon binding the N-terminus of the NarG metalloprotein partner. By the study of singly spin-labeled proteins, the E119 residue present in a conserved elongated hydrophobic groove of NarJ was shown to be part of the interaction site. Moreover, doubly spin-labeled proteins studied by pulsed double electron-electron resonance (DEER) spectroscopy revealed a large and composite distribution of inter-label distances that evolves into a single preexisting one upon complex formation. Additionally, ion mobility mass spectrometry experiments fully support these findings by revealing the existence of several conformers in equilibrium through the distinction of different drift time curves and the selection of one of them upon complex formation. Taken together our work provides a detailed view of the structural flexibility of a dedicated chaperone and suggests that the exquisite recognition and binding of the N-terminus of the metalloprotein is governed by a conformational selection mechanism.

**Citation:** Lorenzi M, Sylvi L, Gerbaud G, Mileo E, Halgand F, et al. (2012) Conformational Selection Underlies Recognition of a Molybdoenzyme by Its Dedicated Chaperone. PLoS ONE 7(11): e49523. doi:10.1371/journal.pone.0049523

**Editor:** Hendrik W. van Veen, University of Cambridge, United Kingdom

**Received:** June 14, 2012; **Accepted:** October 9, 2012; **Published:** November 19, 2012

**Copyright:** © 2012 Lorenzi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by the Agence Nationale de la Recherche (ANR SPINFOLD n°09-BLAN-0100), the Centre National de la Recherche Scientifique (CNRS) and Aix-Marseille University. The authors are also grateful to the EPR facilities available at the national TGE RPE and the Aix-Marseille University EPR center and, to the Conseil Régional de PACA, the city of Marseille and the University of Provence for financial support in the acquisition of instrumentations. The 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.

* E-mail: magalon@imm.cnrs.fr (AM); belle@imm.cnrs.fr (VB)

□ Current address: Plateforme Protéomique Biogenouest, Campus de Beaulieu, Rennes, France

## Introduction

Highly specific and tightly regulated interactions between proteins are essential for any life process. Moreover, proteins are inherently dynamic and often sample a vast ensemble of conformations. It is now widely accepted that the dynamic nature of proteins plays a critical role not only in molecular recognition but also in the evolution of molecular interactions [1,2,3].

One prominent example is molecular chaperones which adopt numerous structurally distinct conformations to recognize and fold a broad number of different substrate proteins [4,5,6,7]. On the opposite, folding and assembly of metal-containing proteins is ensured by chaperones specific of a given substrate [8,9]. These dedicated chaperones literally orchestrate several events ranging from metal centers insertion, folding, membrane targeting and even translocation of their substrates. Current models for the mechanism of coordinated assembly of metalloproteins suggest that dedicated chaperones initially recognize and bind to the N-terminus of unfolded metalloproteins, which then recruit components of metal cofactor biosynthesis machineries to form the

mature and active metalloproteins [10,11,12,13]. One of the most prominent example is the NarJ chaperone ensuring folding and assembly of the membrane-bound respiratory nitrate reductase complex, NarGHI [9,14], member of a large group of molybdenum containing enzymes [9,15]. NarJ coordinates several maturation events through binding to two distinct sites of the NarG catalytic subunit [12]. Interaction with the N-terminus hampers membrane anchoring of an immature complex while interaction at a second site controls sequential insertion of both an iron-sulfur cluster (FSO) and of the molybdenum cofactor [16,17]. Accordingly, absence of NarJ is associated with a premature membrane anchoring of an immature and inactive NarGH complex to NarI [12,17]. The NarJ binding epitope on the N-terminus of NarG is restricted to the first 15 residues without any influence on the binding properties [18]. Calorimetric experiments demonstrated that this interaction is mostly driven by hydrophobic interactions and strongly modulated by protonation of NarJ [18]. Furthermore, NMR showed that the amphiphilic helix adopted by the N-terminus of NarG within the X-ray structure of the NarGHI
complex [19] is conserved in the isolated NarG(1–15) peptide and remains unchanged upon NarJ binding [18].

Despite a low level of sequence identity, resolution by X-ray crystallography of the structure of several dedicated chaperones for molybdenum containing enzymes indicates a conserved all-helical fold [20,21,22,23,24,25] and allows the description of a new family of chaperones (Pfam PF02613) to which NarJ belongs. While X-ray crystallography provided invaluable high-resolution structural information on dedicated chaperones, NMR is, in many cases, a better source to get information about dynamics and flexibility. NMR studies conducted on *E. coli* NarJ suggest the existence of one or several flexible regions by the absence of a number of peaks in the ¹H,¹⁵N-HSQC spectrum which also precluded complete assignment of the residues. Similar experiments conducted in presence of the NarG(1–15) peptide showed a drastic modification of the HSQC spectrum interpreted as the result of a global conformational change as the result of complex formation [18]. Altogether, these NMR studies did not provide a detailed picture of the conformational dynamics. Moreover, no 3D structure of the complex between the chaperone and the N-terminus of its cognate partner for any dedicated chaperones is available so far.

A better understanding of the initial steps of the assembly
process which consist in the exquisite recognition of the
metalloprotein partner through its N-terminus by the dedicated
chaperone called for the use of alternative strategies. Site-directed
spin labeling combined with EPR spectroscopy is a powerful
method for monitoring the structure and dynamics of soluble and
membrane proteins of arbitrary molecular weight [26]. In site-
directed spin labeling, a radical containing nitroxide spin label is
site specifically introduced onto cysteine residues introduced via
site-directed mutagenesis. The resulting EPR spectrum of the
generated paramagnetic side chain can be used to identify
interaction sites within protein complexes or to detect changes in
protein conformations, and the data can be interpreted in terms of
relative domain movement, backbone dynamics and folding or
unfolding events [26,27,28,29,30]. Moreover, pulsed EPR tech-
niques, specifically double electron-electron resonance (DEER)
experiments, allow the measurement of long-range distances and
distance distributions in multi-labeled proteins [31]. Additional
information can be gained through the use of ion mobility (IM)
experiments coupled to mass spectrometry which allow differen-
tiating molecules according to their charge state (CS), collision
cross-section (CCS) and shape with the determination of drift time
curves reflecting gas phase conformation of proteins and peptides
[32,33,34]. By preserving protein-ligand complexes and the
conformational ensemble of proteins in the gas phase, IM has
emerged as a powerful tool for the study of macromolecular
structures [35,36,37,38]. Such device allows inferring insights on
ligand-protein and protein-protein interactions, and unveiling
conformational changes upon ligand binding [39,40,41,42].

In this study, we aimed to provide new information on protein
dynamics of a dedicated chaperone during the binding process to
the N-terminus of its metalloprotein partner by the combination of
the two above-mentioned complementary biophysical approaches
using the NarJ/NarG(1-15) peptide as a working model. Here, we
clearly identified the location of the interaction site between the
two partners and analyzed in detail the structural flexibility
properties exhibited by the dedicated chaperone NarJ. In addition,
we demonstrated that partner binding results from the selection of
an accessible conformation and its further rearrangement induced
upon recognition.

**Materials and Methods**

Protein expression and purification

Mutants were generated by site-directed mutagenesis of the
pDSNarJT-6His and mutations were verified by DNA sequencing.
Overexpression and purification of NarJ truncated of the last 50
amino acid residues (referred as NarJT) and carrying a C-terminal
hexahistidine tag were carried out as described previously [14,18].

N-terminal NarG peptide

The NarG(1–15) peptide was chemically synthesized and
purified by Synprosis (Marseille, France) as previously reported
[18].

Spin labeling of NarJT

Prior to spin labeling, reduction of cysteine residues was carried
out by incubating the NarJT variants in 20 mM of DTT for
30 min at room temperature. DTT was removed by PD10
desalting column (GE Healthcare) using a 50 mM Tris-HCl
pH 8.0, 100 mM NaCl buffer. Spin label 1-oxyl-2,2,5,5-tetra-
methyl-pyrroline-3-methyl methanethiosulfonate, MTSL, (Toron-
to Research Chemicals Inc.) was immediately added to the sample
at a molar excess of 10 for mono or 20 for double labeling. The
reaction was carried out during one hour for mono spin labeling
and four hours for double spin labeling with four consecutive
additions of MTSL (one per hour) in the dark and at 277 K under
gentle stirring and a continuous flow of argon. The excess of
unbound spin label was removed by PD10 desalting column with
the same elution buffer described above.

Circular dichroïsm

CD spectra were recorded on a Jasco 815 CD spectrometer using 1-mm thick quartz cells in 10 mM sodium phosphate pH 7.5 at 296 K. CD spectra were measured from 260 to 190 nm, at 20 nm/min and were averaged from 2 scans. The spectra were corrected for buffer signal. Mean ellipticity values per residue ([θ]) were calculated as described previously [43]. Protein concentrations of 0.1 mg/mL were used.

Tryptophan fluorescence spectroscopy

Intrinsic tryptophan fluorescence emission was measured by
using a Fluorolog FL3-21 spectrofluorimeter and a 1-cm
pathlength cuvette at 296 K. A NarJT concentration of 2 µM in
a buffer containing 50 mM Tris-HCl pH 7.5, 1 mM MgCl₂ was
used. Increasing concentrations of NarG(1-15) peptide at 50 µM
(initial concentration) were then added and the emission fluores-
cence was scanned. Excitation was at 285 nm and emission was
recorded from 310 to 410 nm with excitation and emission band-
passes at 2 nm. Binding of the peptide was monitored by recording
the variation of intrinsic tryptophan fluorescence of NarJT
produced after addition of increasing concentrations of the
peptide. Corrections for both the variation of volume and the
inner-filter effect of the peptide were performed under the same
conditions by using N-acetyl-L-tryptophanamide instead of NarJT.
A 30% decrease of the tryptophan fluorescence was maximally
observed upon peptide addition. Dissociation constants were
calculated by plotting relative fluorescence peak integration
against ligand concentration and curve-fitting using Origin
software.

Continuous Wave EPR

Room temperature experiment (296 K) EPR spectra were
recorded on an ESP 300E Bruker spectrometer equipped with an
ELEXSYS Super High Sensitivity resonator operating at X-band (9.9 GHz). The microwave power was 10 mW, the magnetic field modulation amplitude was optimized to avoid over-modulation of the signal (in the range of 0.1 mT to 0.3 mT) and the frequency modulation was 100 kHz. Spin quantitation was carried out as described previously [44] and gave a labeling yield of ~80% for mono spin labeling (H21C, Q104C, Q149C, and E119C) and 190% for double spin labeling (H21C/Q104C). Labeled NarJT variants were studied in the absence and presence of the NarG(1–15) peptide in a 10-fold molar excess in a 50 mM Tris-HCl pH 7.9, 100 mM NaCl buffer.

## Simulation of the EPR spectra

The EPR spectra recorded at room temperature were simulated using the EPRSIM-C software program. This program, kindly provided by Dr. J. Strancar (University of Ljubljana, Slovenia), is based on the so-called motional-restricted fast-motion approximation described in details in [45]. Such analysis is well adapted to analyze structural changes in proteins [46]. Briefly, the partial averaging of the hyperfine and *g* tensors is described in this model by the following parameters: an effective rotational correlation time τ and two angles (θ₀ and φ₀) corresponding respectively to the amplitude and the anisotropy of the spin label rotational motion within a cone. The normalization of the two last parameters by $\Omega = (\theta₀\phi₀)/(\pi/2)^2$ represents the free rotational space, which varies from zero (totally restricted movement) to 1 (totally unrestricted movement). Two other parameters were required for simulation: a residual width (ω) and a scalar parameter linked to the polarity of the probe environment (pA) for the adjustment of the hyperfine tensor principal values. Fits were performed by minimizing the χ² value.

## Double electron electron resonance (DEER) experiments and Distance Analysis

NarJT samples with and without NarG(1–15) peptide were prepared in a 50 mM Tris-HCl pH 7.9, 100 mM NaCl buffer. In a second set of experiments, to prevent heterogenous protein concentrations that lead to a shortening of the spin label relaxation time and then to a decrease of the DEER signal to noise ratio, the NarJT samples were prepared in the presence of 30% (v/v) glycerol used as cryoprotectant. As glycerol is known to potentially modify binding affinities and in particular hydrophobic interactions, it has been added only after complex formation and just before rapid freezing in liquid nitrogen. For each sample, the final protein concentration was 115 µM in each case. Experiments were conducted on a Bruker ELEXSYS E580 spectrometer at X-band (9.9 GHz) using the standard MD5 dielectric resonator and at Q-band (34 GHz) using the standard EN 5107D2 resonator. The system was equipped with an Oxford helium temperature regulation unit and the data were acquired at 60 K. This temperature has been optimized according to the relaxation times measured at variable temperatures in the range of 20–100 K with 10 K steps. The four-pulse DEER experiment ($(\pi/2) v_1 - \tau_1 - (\pi)v_1 - \tau - (\pi)v_2 - \tau_1 + \tau_2 - \tau - (\pi)v_1 - \tau_2 - \text{echo})$ was used with observe ($v_1$) pulse durations of 16 ns ($\pi/2$) and 32 ns ($\pi$) at X-band and 20 ns ($\pi/2$) and 40 ns ($\pi$) at Q-band. For both X-band and Q-band experiments, interpulse delays $\tau_1$ of 200 ns were used and $\tau_2$ was adjusted depending on the phase memory constant time, $T_m$, of the sample [47]. Initial $\tau$ of 100 ns was used with increment of 4ns. The pump ELDOR pulse ($v_2$, 32ns) was centered at the central resonance and the observed frequency was set at $v_1 - v_2 \approx 72$ MHz away on the low field side. In both X and Q-band DEER experiments, the total acquisition time was between 8–12 hours corresponding to 40–60 averaged scans. Signal

processing was achieved using the DeerAnalysis2011 software package under Matlab [48]. Background echo decay was corrected by using a homogeneous three-dimensional spin distribution and second order polynomial baseline correction. Tikhonov regularization was applied to the corrected dipolar evolution dataset to obtain inter-spin distance distributions using L-curves [49,50]. For each sample, the optimal regularization factor, representing a compromise between smoothness and resolution, was $\alpha = 100$ according to L curve criterion [49,50].

## Electrospray ionization mass spectrometry (ESI-MS) and Ion Mobility (IM) experiments under non denaturing conditions

ESI-MS spectra of NarJT protein were recorded under non denaturing conditions in the positive ion mode using an electrospray ion source fitted on a quadrupole-time of flight hybrid mass spectrometer equipped with a traveling wave ion mobility (TWIM) cell (Synapt HDMS G1, Waters Corp., Manchester). Protein solutions were introduced at a 3 µL/min flow rate. For ESI-MS studies, backing pressure, capillary voltage, sampling cone voltage, extracting cone voltage, source and desolvation temperatures were 4 mbar, 3 kV, 130 V, 4 V, 373 K and 423 K, respectively. For IM experiments performed under the TWIM mode general instrumental parameters were kept as above. Specific parameters related to IM segregation of ions were: gas pressure in ion mobility cell = 0.6 mbar; gas flow rate in the trap cell = 4.5 mL/min; wave height and wave velocity of the trap and transfer cells were 0.2 V; 300 m.sec⁻¹ and 3 V and 248 m.sec⁻¹, respectively. Wave height and wave velocity in the ion mobility cell were 6 V and 150 m.sec⁻¹, respectively. Mass spectra were recorded for 2 min and 15 min for ESI-MS and IM experiments, respectively. Calibration was achieved using a 1 mg/ ml Cesium Iodide solution dissolved in 70% 2-propanol (Sigma Aldrich). Prior to MS analyses NarJT protein was extensively washed using ammonium acetate solutions from 1 M to 25 mM concentration range and used at a working concentration of 30 µM. MS spectra analysis gave experimental molecular weights of 22000.0±0.5 Da for NarJT and 1951.16 Da±0.25 Da for NarG(1–15) peptide, in agreement with the theoretical values (Table S2). NarJT:NarG(1–15) peptide ratios of 1:1 and 1:3 were used. For CCS calculation of NarJT, a myoglobin solution at 10 µM under denaturing condition was used as calibrate. Rigorous identical experimental IM parameters were used to obtain a calibration curve using the Drift Scope CCS calculator™.

## Results

### Generation of surface-exposed cysteine NarJT variants

A prerequisite for site-directed spin labeling studies is the site-specific introduction of cysteine residues into regions of interest. In absence of structural data for *E. coli* NarJ, positions for the introduction of cysteine residues were selected from the structural model of a truncated form of NarJ named NarJT [18]. This 3D model has been built by homology modeling and was truncated from its 50 C-terminal amino acids because of a lack of similarity. The corresponding NarJT protein was proven to be a valuable construct to study the molecular mechanism of recognition and binding of the N-terminus of the metalloprotein partner as its binding parameters towards the NarG(1–15) peptide remained unchanged as compared with full-length NarJ [18]. Furthermore, docking calculation analysis predicted that the NarG(1–15) peptide interacts within a highly conserved elongated hydrophobic groove of NarJT [18]. With respect to this hydrophobic cavity, four
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Model structure of NarJT. The positions of the four residues targeted for cysteine mutations and spin labeling are indicated. The position of the single Trp residue (W19) used for binding affinity measurement is indicated as well as the estimated distances between the Cα of H21 and Q104 or H21 and E119. NarJT is represented in cartoon by Pymol (http://www.pymol.org). doi:10.1371/journal.pone.0049523.g001</figcaption></figure>

positions were selected for mono-spin labeling: position 119 located inside the cavity, 149 on the opposite side, 21 and 104 on both sides of the cavity (Figure 1). Moreover, a double cysteine variant of NarJT (H21C,Q104C) has been produced.

To determine a possible influence of the mutations and of the presence of the spin label on the structural integrity of NarJT, CD spectra were recorded for all NarJT variants before and after labeling. All CD spectra were similar, showing a shape typical of all-helical proteins with two characteristic minima found at 222 and 208 nm (Figure S1). Hence, neither the substitutions nor the spin label influences the global folding of NarJT (Figure S1). Furthermore NarG(1–15) peptide binding activity of the NarJT variants (before and after spin labeling) was evaluated by intrinsic tryptophan fluorescence quenching thanks to the presence of a unique tryptophan residue (W19) located within the above-mentioned hydrophobic cavity (Figure 1). The data clearly showed a nearly unchanged dissociation constant (Table S1). Altogether, CD and fluorescence data demonstrate that neither mutation nor labeling had an effect on both structural integrity and the ability of NarJT to bind its partner.

## Dynamics of nitroxide side chain grafted on NarJT reveals the interaction site with NarG(1–15) peptide

Room temperature EPR spectra of nitroxide spin label grafted on proteins are very sensitive to the motion of the label side chain leading to a partial motional averaging of the anisotropic components of the g- and hyperfine tensors [26,30,51,52]. Figure 2 (upper panels) shows room temperature EPR spectra of the four labeled NarJT proteins recorded in the absence and presence of NarG(1–15) peptide. The visual inspection of the different spectra clearly shows a drastic change of the EPR spectrum after peptide binding only for the label at position 119. All EPR spectra were simulated so as to achieve a more detailed description of the impact of the NarG(1–15) peptide on the spin label mobility. The simulation allows decomposition of an EPR spectrum into possibly different components, each of them being described by two main parameters: the effective rotational correlation time τ and the free rotational space Ω (see Materials

and Methods). Figure 2 (bottom panels) gives a representation of the results as [τ-Ω] plots for the different positions of the spin label either in the absence or in the presence of NarG(1–15) peptide, with the surface of the spheres being proportional to the relative contribution of the various spectral components (Table 1).

For positions 21 and 104 located on both sides of the hydrophobic cavity, the corresponding EPR spectra are well-simulated with two components of unequal weights (Figure 2). Upon addition of the NarG(1–15) peptide, no significant EPR spectral change is observed for the main component of the position 21, indicating that the environment of the spin label remains unchanged after peptide binding. For the main component of the spin label grafted at position 104, peptide binding induces a slight spectral change as revealed by the shift of the sphere in the [τ-Ω] plot. This observation indicates that small structural rearrangements occur in the region surrounding the residue 104 upon peptide binding. For these two positions, the minor component reflects a very high mobility of the spin label (Ω close to 1 and low values of τ) that could result from free radicals in solution or a small fraction of labeled unfolded NarJT (Figure 2, Table 1). For positions 119 and 149, the fact that two components of relatively equal weights are required to obtain a good fit indicates that the label experiences different environments (Figure 2 and Table 1). This likely reflects the occurrence of conformational sub-ensembles, possibly corresponding to different orientations of the spin label. Indeed, it has been demonstrated that two components EPR spectra can result from two rotameric states of the nitroxide side chain, each experiencing a unique environment, a behavior frequently found at helical sites [44,53,54,55,56]. For position 119 located in the hydrophobic cavity, a drastic change in the EPR spectrum is observed consecutive to the addition of NarG(1–15) peptide (Figure 2). These changes are well illustrated in the [τ-Ω] plot where the two components of equal weights in the absence of peptide evolve into one major component (89%) and a minor one (11%). The latter can be attributed to an unbound fraction of NarJT as the Ω and the τ parameters are unchanged as compared to NarJT alone. The major component reflects a highly restricted mobility of the spin label, as indicated by the very low value of Ω, suggesting that the label is in a buried site consecutive to the association of the peptide with NarJT. On the opposite side of the hydrophobic cavity (position 149), no significant variation in the EPR spectra is observed, indicating that the structural characteristics around the probe remain unchanged upon peptide binding. Altogether, these results demonstrate that the region surrounding the residue 119 is part of the interaction site with the N-terminus of the NarG metalloprotein partner.

## Flexibility and conformational dynamics of NarJT upon NarG(1–15) peptide binding as deduced from DEER experiments

Inter-spin distance measurement by DEER techniques is based on the measurement of the dipolar coupling between two interacting spins distant from r, a quantity that is proportional to 1/r³. In the DEER sequence [47], the dipolar field modulates the spin-echo envelope and the oscillations of this envelope visualized in the time traces are used to calculate inter-spin distances. Typically, this technique is efficient to measure, under low temperature conditions, inter-spin distances in a wide range going from ~1.8 nm up to 6.0 nm in proteins [31,44,57], thus being of interest for the study of conformational transitions captured in the frozen state. Through the use of a double-cysteine mutant (H21C/Q104C) with spin labels encompassing the binding site, DEER experiments conducted at different frequencies allowed to measure inter-spin distance distribution and its variation induced by the
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Dynamic studies of the different labeled sites of NarJT. Upper panel: Room temperature EPR spectra of the different labeled NarJT variants in the absence (A) and presence of 10-fold molar excess of the NarG(1–15) peptide (B) in 50 mM Tris-HCl pH 7.9, 100 mM NaCl buffer. Simulated spectra (dotted grey line) are superimposed on the experimental spectra (solid black line). Simulation of the spectra was conducted as described previously [44]. Lower panel: simulated</figcaption></figure>

parameters given as [τ-Ω] plots, with τ the effective rotational correlation time in ns and Ω the normalized free rotational space, in the absence (dark gray spheres) and presence (light gray spheres) of the peptide. The area of the spheres is proportional to the relative contribution of the EPR spectral shape components.
doi:10.1371/journal.pone.0049523.g002

NarG(1–15) peptide binding process. Figure 3 shows the background corrected X-band DEER time traces acquired at low temperature from the bilabeled NarJT in the absence (A) and presence (C) of the NarG(1–15) peptide and the resulting distance distribution obtained by using Tikhonov regularization (B and D) [48] (see Materials and Methods). In the absence of peptide, the distance distribution is broad and composite with the distinction of three maxima (Figure 3B). Such behavior was found to be reproducible using different sample preparations (data not shown). Therefore, this observation suggests the existence of multiple conformations of NarJT. Interestingly, in the presence of the peptide a narrower distance distribution is observed with a unique maximum centered at 2.7 nm (Figure 3D).

Glycerol is often used in DEER experiments prior freezing of the protein sample as a vitrifying agent to prevent excluded volume effects and heterogenous protein concentrations which lead to a drop of the phase memory time ($T_m$) and limit the DEER time range [58]. The immediate consequence of the presence of glycerol is generally a longer time for signal evolution leading to a better reliability and resolution of the calculated distance distribution [59,60]. In this aim, NarJT samples were studied in the presence of glycerol (30% v/v) both in the absence and presence of Nar(1–15) peptide (Figure 4 A,B,C,D). In the absence of the peptide, the different maxima detected previously are less resolved with the observation of shoulders confirming the presence of multiple conformations (Figure 4B). As previously observed, when NarJT is bound to NarG(1–15) peptide, the distribution of distance converges towards a single one centered at 2.9 nm. The difference in inter-spin distance distribution detected without and with glycerol can be attributed to the effect of the cryoprotectant that has been shown to play a role in the compactness of proteins [58,61].

Finally, DEER experiments were performed at higher frequency (Q-band) for the samples containing glycerol (30% v/v) (Figure 4 E,F,G,H). As demonstrated by Ghimire and coworkers [62], Q-band (34 GHz) increases DEER sensitivity as compared to X-band (9.4 GHz), thus leading to a better signal to noise ratio of the DEER time traces and consequently a more reliable distance distribution. Improved sensitivity can also enable the detection out of the noise of low-frequency components (longer inter-spin distance). In our case, the distance distributions deduced from the Q-band DEER time traces are very similar to those calculated from X-band (Figure 4), confirming that peptide binding induces a transition in distance distribution going from multiple maxima to a single one with a global increase of the distance between the labeled protein regions. This observation highlights the structural dynamics of NarJT during the initial steps of assembly of the metalloprotein partner.

## ESI-MS and ESI-TWIM-MS experiments substantiate a conformational selection mechanism for NarG(1–15) peptide binding

To gain further insights into the structural flexibility observed through a large distribution of inter-label distances, mass spectrometry experiments performed under non denaturing conditions with or without ion mobility segregation were carried out on NarJT alone (Figure 5A) and in presence of the NarG(1–15)
<table><caption>**Table 1.** Simulation parameters of room temperature EPR spectra.</caption><thead><tr><th rowspan="2">Sample</th><th rowspan="2">Component</th><th colspan="3">- NarG(1-15) peptide</th><th colspan="3">+ NarG(1-15) peptide</th></tr><tr><th>τ (ns)</th><th>Ω</th><th>%</th><th>τ (ns)</th><th>Ω</th><th>%</th></tr></thead><tbody><tr><td rowspan="2">H21C</td><td>1</td><td>6.0±0.6</td><td>0.97±0.05</td><td>89</td><td>4.7±0.6</td><td>0.95±0.05</td><td>96</td></tr><tr><td>2</td><td>0.4±0.1</td><td>0.97±0.05</td><td>11</td><td>0.10±0.05</td><td>1.00±0.05</td><td>4</td></tr><tr><td rowspan="2">Q104C</td><td>1</td><td>2.8±0.3</td><td>0.75±0.05</td><td>79</td><td>4.4±0.5</td><td>0.80±0.05</td><td>79</td></tr><tr><td>2</td><td>0.24±0.05</td><td>0.96±0.05</td><td>21</td><td>0.26±0.05</td><td>0.94±0.05</td><td>21</td></tr><tr><td rowspan="2">E119C</td><td>1</td><td>0.9±0.2</td><td>0.40±0.05</td><td>53</td><td>0.8±0.2</td><td>0.23±0.05</td><td>89</td></tr><tr><td>2</td><td>3.5±0.4</td><td>0.64±0.05</td><td>47</td><td>3.5±0.4</td><td>0.64±0.05</td><td>11</td></tr><tr><td rowspan="2">Q149C</td><td>1</td><td>0.9±0.2</td><td>0.85±0.05</td><td>53</td><td>0.8±0.2</td><td>0.86±0.05</td><td>52</td></tr><tr><td>2</td><td>1.3±0.2</td><td>0.70±0.05</td><td>47</td><td>1.3±0.2</td><td>0.68±0.05</td><td>48</td></tr></tbody></table>

Parameters (effective rotational correlation time τ, free conformational space Ω and proportion in %) extracted from the simulation of the RT EPR spectra of the labeled NarJT at positions 21, 104, 119 and 149 either in the absence or in the presence of a molar excess of 10 of NarG(1-15). Simulations were performed using the EPRSIM-C software [44].

doi:10.1371/journal.pone.0049523.t001

peptide (Figure 5B). Interestingly the ESI-MS experiment revealed a tri-modal Charge State Distribution (CSD) suggesting that NarJT exists under at least three conformers: namely conformer A corresponding to +12 to +21 charge states (CS), conformer B corresponding to +9 to +11 CSs and conformer C corresponding to +6 to +8 CSs (Figure 5A). Relative to the total ion current, conformers A, B and C account for 25%, 73% and 2%, respectively. Calculation of the average charge state (Z<sub>av</sub>) (Table 2) allows to propose a compactness scale based on protonation propensity [63], with conformer A being the most open form and conformer C the most compact one. Indeed, folded proteins with compact structures give rise to gaseous ions carrying a relatively low number of charges and presenting a narrow CSD centered on high m/z values. Final proof of the existence of different conformers in equilibrium was provided by IM experiments. As emphasized in Figure 5A NarJT drift chart allowed distinguishing 3 different drift time curves supporting the presence of 3 conformers in equilibrium in agreement with distinct CSDs.

To address the impact of NarG(1-15) peptide binding on NarJT conformational ensemble, mass spectrum was recorded under the same experimental conditions (Figure 5B). In addition to the +2 and +1 CS of the free NarG(1-15) peptide, NarJT alone or complexed to one or two peptide molecules were observed. Experimental masses of these respective species were in agreement with predicted ones (Table S2). Relative to the total ion current, ratios of the different NarJT species were 44%, 54% and 2%, for NarJT alone, NarJT complexed to its peptide and NarJT complexed with two peptides, respectively. The very weak abundance of the NarJT protein complexed with two peptides and the concomitant increase of this species using a higher peptide protein ratio suggest a non-specific binding of the second peptide molecule. Interestingly, the tri-modal CSD observed for NarJT alone (Figure 5A) was replaced herein by two new CSDs having close patterns and corresponding to NarJT alone (CSD = +7 to +10 CSs) and NarJT complexed to its partner peptide (CSD = +7 to +11 CSs), respectively (Figure 5B). Noteworthy is the fact that CSDs and drift time curves of these two new components were close to the major conformer B observed for NarJT alone and as a consequence were renamed B* (NarJT free of peptide) and B*P (NarJT in complex with one peptide) (Figure 5). Origin of the B* conformer was given considering that, in previous work, NarJT-peptide binding was shown to be mostly entropy-driven through hydrophobic interactions [18]. Thereof, the complex is most likely

partly dissociated upon collision activation in the gas phase which is detrimental for hydrophobic interactions [64] and resulting in nearly equal amount of bound (B*P) and unbound (B*) forms. In addition, it cannot be inferred from close CSDs values between B*, B*P and B conformers that they all have the same conformation.

To better document the conformational changes occurring when NarG(1-15) peptide was added to NarJT, the Z<sub>av</sub> was calculated. It appears that both B* and B*P undergo a decrease in the Z<sub>av</sub> values as compared to the conformer B, indicating a reduction of the accessible protein surface for protonation (Table 2). In this context the additional charge observed for B*P was assumed to be carried out by NarG(1-15) peptide. Finally, averaged collision cross-section (CCS<sub>av</sub>) values from the +8 to +10

<figure><img src="image_4.png" /><figcaption>**Figure 3. Distance distribution from X-band DEER experiments.**
(A) and (C): Background corrected DEER time traces (solid line) from a four-pulse DEER experiment on the bi-labeled NarJT (H21C/Q104C) in the absence (A) and presence of 10-fold molar excess of the NarG(1-15) peptide (C) in 50 mM Tris-HCl pH 7.9, 100 mM NaCl buffer. Fits are shown as superimposed gray dotted lines. (B) and (D): Inter-label distance distributions calculated from the DEER time traces by using Tikhonov regularization (DeerAnalysis2011 software package [48]) in the absence (B) and presence of the peptide (D).</figcaption></figure>

doi:10.1371/journal.pone.0049523.g003
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Distance distribution from X-band and Q-band DEER experiments in the presence of 30% glycerol. X-band DEER time traces are shown in (A) to (D) while Q-band DEER time traces are shown in (E) to (G). All samples are prepared in 50 mM Tris-HCl pH 7.9, 100 mM NaCl, 30% (v/v) glycerol buffer. (A) and (C): Background corrected DEER time traces (solid line) from a four-pulse DEER experiment on the bi-labeled NarJT (H21C/Q104C) in the absence (A) and presence of 10-fold molar excess of the NarG(1–15) peptide (C). Fits are shown as superimposed gray dotted lines. (B) and (D): Inter-label distance distributions calculated from the DEER time traces by using Tikhonov regularization (DeerAnalysis2011 software package [48]) in the absence (B) and presence of the peptide (D). (E) and (G): Background corrected DEER time traces (solid line) from a four-pulse DEER experiment on the bi-labeled NarJT (H21C/Q104C) in the absence (E) and presence of 10-fold molar excess of the NarG(1–15) peptide (G). Fits are shown as superimposed gray dotted lines. (F) and (H): Inter-label distance distributions calculated from the DEER time traces by using Tikhonov regularization in the absence (F) and presence of the peptide (H).
doi:10.1371/journal.pone.0049523.g004</figcaption></figure>

CS of B, B\* and B\*P were calculated (Table 2). While the CCS<sub>av</sub> of B\* increased as compared to B\*P due to complex dissociation in the gas phase, a small decrease of the CCS<sub>av</sub> occurred as the result of complex formation when comparing B to B\*P. Altogether, these results suggest that binding of the NarG(1–15) peptide led to the selection of conformer B followed by subtle structural changes.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Mass spectra and drift charts recorded under non denaturing conditions of (A) NarJT and (B) NarJT in the presence of the NarG(1–15) peptide. Each panel (A and B) is composed of the ESI-MS spectra (top) and the corresponding drift chart represented in drift time versus m/z (bottom) obtained by Ion Mobility. Ratios intensities according to the TIC of each spectrum are given in drift charts. Curves depicted in black correspond to drift time curves identified for each NarJT conformers. In panel A, respective charge states of NarJT conformers (Conformer A=filled circle circled, Conformer B=filled circles, Conformer C=diamond) are labeled both on mass spectrum and mobility chart. In panel B, open circles correspond to NarJT alone (conformer B\*); filled circles correspond to NarJT complexed with its partner peptide (conformer B\*P); filled squares correspond to NarJT complexed with 2 molecules of its partner peptide (conformer B\*2P). Arrows in panel B denote +1 and +2 CS of NarG(1–15) peptide observed all along the TWIM cell. A drift time value of 1 bin corresponds to 90 µs.
doi:10.1371/journal.pone.0049523.g005</figcaption></figure>

Overall, our ESI-MS and ESI-TWIM-MS results provide additional support to the structural dynamics of NarJT and suggest that complex formation arise from the conformational selection of one of the NarJT conformers *ie* conformer B.
<table><caption>**Table 2.** ESI-MS and IM data analysis depicting conformational species of NarJT either alone or complexed with NarG(1–15) peptide.</caption><thead><tr><th>Sample</th><th>Z<sub>av</sub></th><th>Ratios in %</th><th>CSD</th><th>CCS<sub>av</sub> in Å<sup>2</sup></th></tr></thead><tbody><tr><td>NarJT conformer A</td><td>16.09</td><td>25</td><td>+12–+21</td><td>N.C.</td></tr><tr><td>NarJT conformer B</td><td>9.97</td><td>73</td><td>+9–+11</td><td>2537</td></tr><tr><td>NarJT conformer C</td><td>7.14</td><td>2</td><td>+6–+8</td><td>N.C.</td></tr><tr><td>NarJT (B*)</td><td>8.65</td><td>44</td><td>+7–+10</td><td>2600</td></tr><tr><td>NarJT (B*P)</td><td>9.68</td><td>54</td><td>+7–+11</td><td>2453</td></tr><tr><td>NarJT (B*2P)</td><td>8.75</td><td>2</td><td>+8–+11</td><td>N.C.</td></tr></tbody></table>

Average charge states ($Z_{av}$), ratios intensities in % relative to TIC (indicated values were given with a ±0.5% standard error), charge state distribution (CSD) and averaged collision cross section (CCS$_{av}$) values in Å$^2$ calculated for +8 to +10 CS of NarJT conformers detected when alone or in mixture with the NarG(1–15) peptide. B*2P stands for NarJT in complex with two peptides.
doi:10.1371/journal.pone.0049523.t002

Discussion

An important role for protein flexibility has been widely discussed in the literature [1,2,3,65] and represents a technical challenge when describing protein interactions. Previous NMR studies conducted on *E. coli* NarJ both in absence and in presence of the NarG(1–15) peptide showed the absence of a number of peaks in the ¹H,¹⁵N-HSQC spectrum which precluded complete assignment of the residues [18]. Moreover, global change of the ¹H,¹⁵N-HSQC spectrum of NarJ upon peptide binding has been interpreted as the result of a conformational change upon complex formation without inferring the underlying molecular mechanisms. Altogether, these NMR studies conducted for the first time on a dedicated chaperone of the NarJ family did not provide a detailed picture of the conformational dynamics. Herein the combination of site-directed spin labeling EPR spectroscopy and IM mass spectrometry has provided a unique window on the conformational substates of the dedicated chaperone NarJ and during the partner binding process, revealing distinct molecular species and conformational dynamics. More importantly, our study explains how structural flexibility of the dedicated chaperone can be exploited for binding of the N-terminus of the metalloprotein during the folding and assembly process.

Strong evidences for the existence of a structural flexibility came from two complementary approaches revealing that the NarJT protein was represented by at least three discrete conformations. DEER experiments performed on a doubly-labeled NarJT variant in frozen state show a broad distance distribution with multiple maxima (Figure 3B and Figure 4B) while IM experiments allow distinguishing three conformers (Figure 5A). Importantly, such structural equilibrium was shown to disappear to the profit of a distance distribution centered on one single maximum or a single protein conformation when the NarG(1–15) peptide was added. It is worth to notice that, by the use of singly spin-labeled proteins, the E119 residue present in a funnel-shaped hydrophobic cavity conserved in all members of the dedicated chaperone is identified as being part of the binding site for the N-terminus of NarG. In particular, the singular location of the E119 residue within this cavity makes it a valuable tool to sense NarG(1–15) peptide binding. Simulation of the EPR spectra revealed two components of equal weights for the E119C variant (Figure 2) reflecting two different environments of the spin probe. Remarkably, the drastic change in the EPR spectrum observed upon complex formation revealed the conversion of these two populations into a major one

of more restricted conformational space resulting from the presence of the NarG(1–15) peptide. These observations are reminiscent of the redistribution of the conformational substates of NarJT obtained upon complex formation from IM experiments.

From a mechanistic point of view, it is important to ask when the conformational transition occurs during the binding process. Several models have been proposed to explain the conformational changes observed between the bound and unbound forms of proteins. In the induced fit mechanism the ligand induces a conformational change upon binding while in the preexisting equilibrium model protein samples an ensemble of conformations at equilibrium conditions and a fraction of which is predisposed to recognize and bind selectively a particular ligand.

In our study, two lines of evidence are highly suggestive of a conformational selection mechanism operating during binding of the N-terminus of NarG by the NarJ chaperone. At first, the distance distribution with one single maximum revealed by DEER in the bound form of NarJT was part of the large and composite distance distribution obtained in the unbound form (Figures 3 and 4). Secondly, the new NarJT conformations (B* and B*P) observed by IM after peptide addition are highly similar to the major conformer B observed for NarJT alone since CSDs and drift time curves of B* and B*P conformers were close to B (Figure 5) suggesting that B* conformer would derive from B. Further analysis revealed a decrease of $Z_{av}$ value associated with B* and B*P that is related to a reduced surface accessibility for protonation as compared to B (Table 2). It is noticeable that B*P $Z_{av}$ value was found to carry one more charge in average as compared to B*. We hypothesized that the presence of an additional charge on B*P could be linked to the binding of the NarG(1–15) peptide forming an amphiphilic α-helix [18,19] with the charge retained specifically on the outward polar α-helix face. Second CCS$_{av}$ values that reflect global gas phase conformation were compared to gain better insights upon global conformational changes occurring both when partner peptide binds to B (ie B*P) or dissociate from B*P (ie B*). The observation of a decrease in CCS$_{av}$ value when the peptide binds to B conformer was consistent with previous calorimetric data showing that binding induces conformational changes leading to a more compact structure as deduced from higher excess partial molar heat capacity of the complex [18]. Interestingly, dissociation of the complex in the gas phase gave rise to a higher CCS$_{av}$ value for the unbound B* form which denote a more open structure.

One tentative model that is consistent with our experimental observations is that NarJT protein is represented by at least 3 conformers in equilibrium and that peptide binding results from the selection of an accessible conformation and its further rearrangement induced upon peptide recognition. An additional interesting observation is that the entropy of binding ($\Delta S_{bind}$) obtainable by calorimetric methods is composed of contributions associated with the protein, the ligand and the solvent. Zakian et al showed a large and positive $\Delta S_{bind}$, which accounts for most of the binding process in the NarJT-NarG(1–15) peptide complex [18]. If such positive entropy variation was originally interpreted as the result of hydrophobic contacts or the loss of water-mediated hydrogen bonds, our results likely support the additional contribution of conformational entropy which has been considered important for molecular recognition [66]. At this stage, it is worth mentioning the absence of structural information to date concerning the mode of interaction between the dedicated chaperones and the N-terminus of their cognate metalloprotein partner, in particular no X-ray structure of such a complex is available so far. During the last decade, many debates have concerned the exact location of the N-terminus binding site within
the dedicated chaperone. By revealing structural dynamic at play
during the recognition and binding process, our results provide
herein a framework for understanding the mode of interaction.

Several studies reported that proteins have often evolved to bind
multiple targets, including proteins, peptides, DNA, and small
molecule substrates. The simplest way to achieve multispecificity is
to recognize each target through a distinct binding interface and
multiple binding domains. However, X-ray structures of several
dedicated chaperones revealed the existence of single-domain
proteins sharing a common all-helical fold with only a low level of
sequence identity (Pfam PF02613) [9]. In such cases, conforma-
tional flexibility is the most frequent explanation for multi-
specificity [67,68]. As demonstrated here for the first time, NarJ
samples different conformational substates and we anticipate that
this situation may hold true for other members of the family.
Structure of those dedicated chaperones may have evolved to
allow interaction with a diverse range of partners such as
components of metal center biosynthesis machineries, or con-
versely, binding partners may have taken advantage of preexisting
conformational diversity [3,69]. Both processes may be important
in the evolution of molecular recognition.

Furthermore no information is available concerning the binding
interface for the multiple partners of NarJ involved in metal center
delivery during the assembly process of the cognate metallopro-
tein. In this context, stabilization of a specific conformer of NarJ
with NarG(1-15) peptide binding and thus redistribution of the
protein conformational ensembles recalls allostery, a mechanism
by which binding of the ligand at one site can affect binding of
others through a propagated change in the protein shape [70,71].
Interestingly, Volkman et al showed that the inactive response
regulator NtrC, a single-domain protein, samples the active state

conformation even in the absence of the ligand [72]. Allostery thus
derives from a redistribution of the conformational ensemble [70].
Overall, the shift in population resulting from peptide binding
could be one of the keys to facilitate subsequent binding of
additional partners at yet unidentified sites of NarJ. Such a
situation may also hold true in other metalloproteins of interest
such as hydrogenases or ureases which rely on dedicated
chaperones for their folding and function [73,74].

**Supporting Information**

<caption>**Figure S1** CD spectra of NarJ, NarJT and all cysteine variants of NarJT before (A) and after (B) spin labeling.
(TIFF)</caption>

<caption>**Table S1** NarG(1–15) peptide binding activity as estimated by intrinsic tryptophan fluorescence quenching.
(PDF)</caption>

<caption>**Table S2** Calculated, theoretical molecular weights and mass errors given in Da and ppm of NarJT species detected under non denaturing conditions.
(PDF)</caption>

**Acknowledgments**

We thank Dr S. Zakian, M. Mure and W. Kaaki for technical support.

**Author Contributions**

Conceived and designed the experiments: GG FH AW HV VB BG AM.
Performed the experiments: ML LS GG EM FH. Analyzed the data: ML
GG EM FH AW HV VB BG AM. Wrote the paper: GG FH VB BG AM.

**References**

1. Boehr DD, Nussinov R, Wright PE. (2009) The role of dynamic conformational ensembles in biomolecular recognition. Nat Chem Biol 5: 789-796.

2. Ma B, Shatsky M, Wolfson HJ, Nussinov R (2002) Multiple diverse ligands binding at a single protein site: a matter of pre-existing populations. Protein Sci 11: 184-197.

3. James LC, Tawfik DS (2003) Conformational diversity and protein evolution-a 60-year-old hypothesis revisited. Trends Biochem Sci 28: 361-368.

4. Braig K, Adams PD, Brunger AT (1995) Conformational variability in the refined structure of the chaperonin GroEL at 2.8 Å resolution. Nat Struct Biol 2: 1083-1094.

5. Roseman AM, Chen S, White H, Braig K, Saibil HR (1996) The chaperonin ATPase cycle: mechanism of allosteric switching and movements of substrate-binding domains in GroEL. Cell 87: 241-251.

6. Krukenberg KA, Forster F, Rice LM, Sali A, Agard DA (2008) Multiple conformations of E. coli Hsp90 in solution: insights into the conformational dynamics of Hsp90. Structure 16: 755-765.

7. Street TO, Lavery LA, Agard DA (2011) Substrate binding drives large-scale conformational changes in the Hsp90 molecular chaperone. Mol Cell 42: 96-105.

8. Sargent F (2007) Constructing the wonders of the bacterial world: biosynthesis of complex enzymes. Microbiology 153: 633-651.

9. Magalon A, Fedor J, Walburger A, Weiner JH (2011) Molybdenum enzymes in bacteria and their maturation. Coord Chem Rev 255: 1159-1178.

10. Oresnik IJ, Ladner CL, Turner RJ (2001) Identification of a twin-arginine leader-binding protein. Mol Microbiol 40: 323-331.

11. Jack RL, Buchanan G, Dubini A, Hatzixanthis K, Palmer T, et al. (2004) Coordinating assembly and export of complex bacterial proteins. Embo J 23: 3962-3972.

12. Vergnes A, Pommier J, Toci R, Blasco F, Giordano G, et al. (2006) NarJ chaperone binds on two distinct sites of the aponitrate reductase of Escherichia coli to coordinate molybdenum cofactor insertion and assembly. J Biol Chem 281: 2170-2176.

13. Genest O, Neumann M, Seduk F, Stocklein W, Mejean V, et al. (2008) Dedicated metallochaperone connects apoenzyme and molybdenum cofactor biosynthesis components. J Biol Chem 283: 21433-21440.

14. Blasco F, Dos Santos JP, Magalon A, Frixon C, Guigliarelli B, et al. (1998) NarJ is a specific chaperone required for molybdenum cofactor assembly in nitrate reductase A of Escherichia coli. Mol Microbiol 28: 435-447.

15. Rothery RA, Workun GJ, Weiner JH (2008) The prokaryotic complex iron-sulfur molybdoenzyme family. Biochim Biophys Acta 1778: 1897-1929.

16. Lanciano P, Vergnes A, Grimaldi S, Guigliarelli B, Magalon A (2007) Biogenesis of a respiratory complex is orchestrated by a single accessory protein. J Biol Chem 282: 17468-17474.

17. Lanciano P, Savoyant A, Grimaldi S, Magalon A, Guigliarelli B, et al. (2007) New method for the spin quantitation of [4Fe-4S](+) clusters with S = (3)/(2). Application to the FS0 center of the NarGHI nitrate reductase from Escherichia coli. J Phys Chem B 111: 13632-13637.

18. Zakian S, Lafitte D, Vergnes A, Pimentel C, Sebban-Kreuzer C, et al. (2010) Basis of recognition between the NarJ chaperone and the N-terminus of the NarG subunit from Escherichia coli nitrate reductase. FEBS J 277: 1886-1895.

19. Bertero MG, Rothery RA, Palak M, Hou C, Lim D, et al. (2003) Insights into the respiratory electron transfer pathway from the structure of nitrate reductase A. Nat Struct Biol 10: 681-687.

20. Tranier S, Iobbi-Nivol C, Birck C, Ilbert M, Mortier-Barriere I, et al. (2003) A novel protein fold and extreme domain swapping in the dimeric TorD chaperone from Shewanella massilia. Structure 11: 165-174.

21. Kirillova O, Chruszcz M, Shumilin IA, Skarina T, Gorodichtchenskaia E, et al. (2007) An extremely SAD case: structure of a putative redox-enzyme maturation protein from Archaeoglobus fulgidus at 3.4 Å resolution. Acta Crystallogr D Biol Crystallogr 63: 348-354.

22. Qiu Y, Zhang R, Binkowski TA, Tereshko V, Joachimiak A, et al. (2008) The 1.38 Å crystal structure of DmsD protein from Salmonella typhimurium, a proofreading chaperone on the Tat pathway. Proteins 71:525-533

23. Ramasamy SK, Clemons WM, Jr (2009) Structure of the twin-arginine signal-binding protein DmsD from Escherichia coli. Acta Crystallogr Sect F Struct Biol Cryst Commun 65: 746-750.

24. Stevens CM, Winstone TM, Turner RJ, Paetzel M (2009) Structural analysis of a monomeric form of the twin-arginine leader peptide binding chaperone Escherichia coli DmsD. J Mol Biol 389: 124-133.

25. Coulhurst SJ, Dawson A, Hunter WN, Sargent F (2012) Conserved signal peptide recognition systems across the prokaryotic domains. Biochemistry 51: 1678-1686.

26. Columbus L, Hubbell WL (2002) A new spin on protein dynamics. Trends Biochem Sci 27: 288-295.

27. Cooper DB, Smith VF, Crane JM, Roth HC, Lilly AA, et al. (2008) SecA, the motor of the secretion machine, binds diverse partners on one interactive surface. J Mol Biol 382: 74-87.

28. Klare JP, Steinhoff HJ (2009) Spin labeling EPR. Photosynthesis Research 102: 377-390.
29. Crane JM, Lilly AA, Randall LL (2010) Characterization of interactions between proteins using site-directed spin labeling and electron paramagnetic resonance spectroscopy. Methods Mol Biol 619: 173-190.

30. Longhi S, Belle V, Fournel A, Guigliarelli B, Carriere F (2011) Probing structural transitions in both structured and disordered proteins using site-directed spin-labeling EPR spectroscopy. J Pept Sci 17: 315-328.

31. Jeschke G, Polyhach Y (2007) Distance measurements on spin-labelled biomacromolecules by pulsed electron paramagnetic resonance. Physical Chemistry Chemical Physics 9: 1895-1910.

32. Giles K, Pringle SD, Worthington KR, Little D, Wildgoose JL, et al. (2004) Applications of a travelling wave-based radio-frequency-only stacked ring ion guide. Rapid Commun Mass Spectrom 18: 2401-2414.

33. Smith DP, Knapman TW, Campuzano I, Malham RW, Berryman JT, et al. (2009) Deciphering drift time measurements from travelling wave ion mobility spectrometry-mass spectrometry studies. Eur J Mass Spectrom 15: 113-130.

34. Utrecht C, Rose RJ, van Duijn E, Lorenzen K, Heck AJ (2010) Ion mobility mass spectrometry of proteins and protein assemblies. Chem Soc Rev 39: 1633-1655.

35. Suckau D, Shi Y, Beu SC, Senko MW, Quinn JP, et al. (1993) Coexisting stable conformations of gaseous protein ions. Proc Natl Acad Sci U S A 90: 790-793.

36. Breuker K, McLafferty FW (2008) Stepwise evolution of protein native structure with electrospray into the gas phase, 10(-12) to 10(2) s. Proc Natl Acad Sci U S A 105: 18145-18152.

37. Rand KD, Pringle SD, Murphy JP 3rd, Fadgen KE, Brown J, et al. (2009) Gas-phase hydrogen/deuterium exchange in a traveling wave ion guide for the examination of protein conformations. Anal Chem 81: 10019-10028.

38. Jurneczko E, Cruickshank F, Porrini M, Nikolova P, Campuzano ID, et al. (2012) Intrinsic disorder in proteins: a challenge for (un)structural biology met by ion mobility-mass spectrometry. Biochem Soc Trans 40: 1021-1026.

39. Williams JP, Phillips HI, Campuzano I, Sadler PJ (2010) Shape changes induced by N-terminal plating of ubiquitin by cisplatin. J Am Soc Mass Spectrom 21: 1097-1106.

40. Atmanene C, Wagner-Rousset E, Malissard M, Chol B, Robert A, et al. (2009) Extending mass spectrometry contribution to therapeutic monoclonal antibody lead optimization: characterization of immune complexes using noncovalent ESI-MS. Anal Chem 81: 6364-6373.

41. Hilton GR, Thalassinos K, Grabenauer M, Sanghera N, Slade SE, et al. (2010) Structural analysis of prion proteins by means of drift cell and traveling wave ion mobility mass spectrometry. J Am Soc Mass Spectrom 21: 845-854.

42. Jenner M, Ellis J, Huang WC, Lloyd Raven E, Roberts GC, et al. (2011) Detection of a Protein Conformational Equilibrium by Electrospray Ionisation-Ion Mobility-Mass Spectrometry. Angew Chem Int Ed Engl 50: 8291-8294.

43. Kelly SM, Price NC (2000) The use of circular dichroism in the investigation of protein structure and function. Curr Protein Pept Sci 1: 349-384.

44. Ranaldi S, Belle V, Woudstra M, Bourgeas R, Guigliarelli B, et al. (2010) Amplitude of Pancreatic Lipase Lid Opening in Solution and Identification of Spin Label Conformational Subensembles by Combining Continuous Wave and Pulsed EPR Spectroscopy and Molecular Dynamics. Biochemistry 49: 2140-2149.

45. Strancar J, Koklic T, Arsov Z, Filipic B, Stopar D, et al. (2005) Spin label EPR-based characterization of biosystem complexity. Journal of Chemical Information and Modeling 45: 394-406.

46. Kavalenka A, Urbancic I, Belle V, Rouger S, Costanzo S, et al. (2010) Conformational analysis of the partially disordered measles virus NTAIL-XD complex by SDL EPR spectroscopy. Biophysical Journal 98: 1055-1064.

47. Pannier M, Veit S, Godt A, Jeschke G, Spiess HW (2000) Dead-time free measurement of dipole-dipole interactions between electron spins. Journal of Magnetic Resonance 142: 331-340.

48. Jeschke G, Chechik V, Ionita P, Godt A, Zimmermann H, et al. (2006) DeerAnalysis2006 - a comprehensive software package for analyzing pulsed ELDOR data. Applied Magnetic Resonance 30: 473-498.

49. Jeschke G (2002) Distance measurements in the nanometer range by pulse EPR. Chemphyschem 3: 927-932.

50. Jeschke G, Panek G, Godt A, Bender A, Paulsen H (2004) Data analysis procedures for pulse ELDOR measurements of broad distance distributions. Applied Magnetic Resonance 26: 223-244.

51. McHaourab HS, Lietzow MA, Hideg K, Hubbell WL (1996) Motion of spin-labeled side chains in T4 lysozyme. Correlation with protein structure and dynamics. Biochemistry 35: 7692-7704.

52. Columbus L, Kalai T, Jeko J, Hideg K, Hubbell WL (2001) Molecular motion of spin labeled side chains in alpha-helices: analysis by variation of side chain structure. Biochemistry 40: 3828-3846.

53. Langen R, Oh KJ, Cascio D, Hubbell WL (2000) Crystal structures of spin labeled T4 lysozyme mutants: implications for the interpretation of EPR spectra in terms of structure. Biochemistry 39: 8396-8405.

54. Guo ZF, Cascio D, Hideg K, Hubbell WL (2008) Structural determinants of nitroxide motion in spin-labeled proteins: Solvent-exposed sites in helix B of T4 lysozyme. Protein Science 17: 228-239.

55. Guo ZF, Cascio D, Hideg K, Kalai T, Hubbell WL (2007) Structural determinants of nitroxide motion in spin-labeled proteins: Tertiary contact and solvent-inaccessible sites in helix G of T4 lysozyme. Protein Science 16: 1069-1086.

56. Bridges MD, Hideg K, Hubbell WL (2010) Resolving Conformational and Rotameric Exchange in Spin-Labeled Proteins Using Saturation Recovery EPR. Applied Magnetic Resonance 37: 363-390.

57. Altenbach C, Kusnetzow AK, Ernst OP, Hofmann KP, Hubbell WL (2008) High-resolution distance mapping in rhodopsin reveals the pattern of helix movement due to activation. Proceedings of the National Academy of Sciences of the United States of America 105: 7439-7444.

58. Vagenende V, Yap MG, Trout BL (2009) Mechanisms of protein stabilization and prevention of protein aggregation by glycerol. Biochemistry 48: 11084-11096.

59. Galiano L, Blackburn ME, Veloro AM, Bonora M, Fanucci GE (2009) Solute effects on spin labels at an aqueous-exposed site in the flap region of HIV-1 protease. J Phys Chem B 113: 1673-1680.

60. Georgieva ER, Roy AS, Grigoryants VM, Borbat PP, Earle KA, et al. (2012) Effect of freezing conditions on distances and their distributions derived from Double Electron Electron Resonance (DEER): a study of doubly-spin-labeled T4 lysozyme. J Magn Reson 216: 69-77.

61. Feng S, Yan YB (2008) Effects of glycerol on the compaction and stability of the wild type and mutated rabbit muscle creatine kinase. Proteins 71: 844-854.

62. Ghimire H, McCarrick RM, Budil DE, Lorigan GA (2009) Significantly improved sensitivity of Q-band PELEDOR/DEER experiments relative to X-band is observed in measuring the intercoil distance of a leucine zipper motif peptide (GCN4-LZ). Biochemistry 48: 5782-5784.

63. Halgand F (2001) Mean charge state and charge state distribution of proteins as structural probes. An electrospray ionisation mass spectrometry study of lysozyme and ribonuclease. Eur J Mass Spectrom 7: 433-439.

64. Loo JA (1997) Studying noncovalent protein complexes by electrospray ionization mass spectrometry. Mass Spectrom Rev 16: 1-23.

65. Wright PE, Dyson HJ (2009) Linking folding and binding. Curr Opin Struct Biol 19: 31-38.

66. Frederick KK, Marlow MS, Valentine KG, Wand AJ (2007) Conformational entropy in molecular recognition by proteins. Nature 448: 325-329.

67. Erijman A, Aizner Y, Shifman JM (2011) Multispecific recognition: mechanism, evolution, and design. Biochemistry 50: 602-611.

68. Schreiber G, Keating AE (2011) Protein binding specificity versus promiscuity. Curr Opin Struct Biol 21: 50-61.

69. Tobi D, Bahar I (2005) Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state. Proc Natl Acad Sci U S A 102: 18908-18913.

70. Gunasekaran K, Ma B, Nussinov R (2004) Is allostery an intrinsic property of all dynamic proteins? Proteins 57: 433-443.

71. Tzeng SR, Kalodimos CG (2011) Protein dynamics and allostery: an NMR view. Curr Opin Struct Biol 21: 62-67.

72. Volkman BF, Lipson D, Wemmer DE, Kern D (2001) Two-state allosteric behavior in a single-domain signaling protein. Science 291: 2429-2433.

73. Watanabe S, Matsumi R, Arai T, Atomi H, Imanaka T, et al. (2007) Crystal structures of [NiFe] hydrogenase maturation proteins HypC, HypD, and HypE: insights into cyanation reaction by thiol redox signaling. Mol Cell 27: 29-40.

74. Musiani F, Bellucci M, Ciurli S (2011) Model structures of Helicobacter pyloriURED/H domains: a putative molecular recognition platform. J Chem Inf Model 51: 1513-1520.