PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosbiolPLoS BiologySubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Biol. Jan 2014; 12(1): e1001776.
Published online Jan 28, 2014. doi:  10.1371/journal.pbio.1001776
PMCID: PMC3904827
Segmental Helical Motions and Dynamical Asymmetry Modulate Histidine Kinase Autophosphorylation
Ariel E. Mechaly,¤ Nathalie Sassoon, Jean-Michel Betton, and Pedro M. Alzari*
Institut Pasteur, Unité de Microbiologie Structurale and CNRS UMR 3528, Paris, France
Frederick Hughson, Academic Editor
Princeton University, United States of America
* E-mail: pedro.alzari/at/pasteur.fr
The authors have declared that no competing interests exist.
The author(s) have made the following declarations about their contributions: Conceived and designed the experiments: AEM JMB PMA. Performed the experiments: AEM NS JMB. Analyzed the data: AEM NS JMB PMA. Wrote the paper: AEM JMB PMA.
¤Current address: Unit of Protein Crystallography, Institut Pasteur de Montevideo, Montevideo, Uruguay.
Received June 14, 2013; Accepted December 13, 2013.
Histidine kinases (HKs) are dimeric receptors that participate in most adaptive responses to environmental changes in prokaryotes. Although it is well established that stimulus perception triggers autophosphorylation in many HKs, little is known on how the input signal propagates through the HAMP domain to control the transient interaction between the histidine-containing and ATP-binding domains during the catalytic reaction. Here we report crystal structures of the full cytoplasmic region of CpxA, a prototypical HK involved in Escherichia coli response to envelope stress. The structural ensemble, which includes the Michaelis complex, unveils HK activation as a highly dynamic process, in which HAMP modulates the segmental mobility of the central HK α-helices to promote a strong conformational and dynamical asymmetry that characterizes the kinase-active state. A mechanical model based on our structural and biochemical data provides insights into HAMP-mediated signal transduction, the autophosphorylation reaction mechanism, and the symmetry-dependent control of HK kinase/phosphatase functional states.
Bacteria use two-component signal transduction systems (TCSs), composed of a histidine kinase (HK) and a downstream response regulator (RR), to perceive a broad range of external stimuli and trigger an appropriate response. Prototypical HKs are transmembrane homodimeric receptors with an external sensor region coupled to a cytoplasmic catalytic core through an adaptor HAMP domain. By a mechanism that is still poorly understood, the external input signal propagates through the HAMP domain to control the activity of the catalytic domain, which adds a phosphate tag, first to the HK itself (“autophosphorylation”), ready for subsequent transfer to the RR. In this study, we report the first detailed structures of the full cytoplasmic domain of a prototypical HK, CpxA, which is involved in the detection and response to damage of the cell envelope of Escherichia coli. The structures provide snapshots of different kinase states along the autophosphorylation reaction. Our combined genetic, biochemical, and structural data point to a mechanical model for HK activation in which the input signal induces coordinated helical motions that regulate catalysis by asymmetrically repositioning the catalytic domains.
Bacteria commonly use two-component signal transduction pathways to couple environmental stimuli to adaptive responses [1]. Typically composed of a sensor histidine kinase (HK) and a response regulator (RR), two-component systems (TCSs) are abundant in bacterial genomes, reflecting the diversity of stimuli that bacteria are capable of perceiving and responding to [2][4]. Independently of the nature of input signals, TCSs share a common signal transduction mechanism [3]. Stimulus perception triggers a two-step phosphorylation cascade consisting in the autophosphorylation of the HK at a conserved histidine residue followed by the transference of the phosphoryl group to an aspartate in the cognate RR, usually a transcription factor that initiates the bacterial response by modifying target gene(s) expression.
HKs are homodimeric receptors with a modular architecture, which in its simplest form consists of a sensor domain and a cytoplasmic transmitter core. Whereas the signal-dependent sensor domains are highly variable, the transmitter core is more conserved. It contains a dimerization and histidine phosphotransfer (DHp) domain coupled to a catalytic ATP-binding (CA) domain. Mainly depending on the cellular location of the sensor domain (either extracytoplasmic, membrane-anchored, or cytosolic), a wide variety of HK architectures exists [4],[5]. However, the prevalent form of HKs is a transmembrane receptor with an extracytoplasmic sensor domain coupled to the cytoplasmic catalytic core through an adaptor HAMP domain [6], a small protein module that is commonly found in Histidine kinases, Adenyl cyclases, Methyl-accepting chemotaxis proteins and Phosphatases. During the last few years, a number of structural studies have shed light into various aspects of the TCS mechanism of signal transduction. Crystal structures of isolated HK sensor domains have revealed the conformational consequences of stimulus perception [7][9], whereas those of HK transmitter cores illustrated the conversion between different functional states through DHp interhelical rearrangements [10],[11], and two structures of cytoplasmic HK–RR complexes provided the structural basis underlying the RR dephosphosphorylation reaction [12],[13]. Nevertheless, some important questions on TCS signal transduction and autophosphorylation mechanisms remain largely unanswered. For instance, it is poorly understood how a stimuli-induced symmetry/asymmetry switch regulates the HK catalytic activities. Moreover, a full structural characterization of the autophosphorylation reaction remains elusive, possibly hampered by the intrinsic flexibility and highly dynamic nature of the catalytic core [10], as well as by the transient nature of the reaction intermediate in HKs [14]. There is also a lack of structural data on prototypic (HAMP-containing) HKs. So far, none of the HKs for which the transmitter core has been crystallized [10][13],[15] belong to the predominant HK subfamily containing HAMP domains directly linked to the catalytic core [4].
To gain insight into the mode of action of prototypical HKs, we chose to focus on Escherichia coli CpxA, the HK component of the Cpx signaling system that regulates an envelope stress response in E. coli [16],[17]. The sensor kinase CpxA senses a variety of envelope stresses, including misfolded proteins, and transduces this signal to the RR CpxR through conserved phosphotransfer reactions. Phosphorylated CpxR then activates transcription of genes whose products are involved in protein physiology in the envelope [18]. Cpx regulon includes a large number of members [19], mainly periplasmic chaperones and proteases [20][23]. In addition to CpxA and CpxR, the Cpx system includes the auxiliary periplamic protein CpxP, a dimeric α-helical protein [24],[25] that could interact with the CpxA periplasmic sensor domain and was shown to tightly regulate autokinase activity in vivo [26] and in vitro [27]. Here we report the 3D structures of the full cytoplasmic region of CpxA in several different crystal forms. These structures show CpxA adopting different asymmetric conformations, including a state that corresponds to the autophosphorylating Michaelis complex and provides insights into the reaction mechanism. The ensemble of structures illustrates the highly dynamic and strongly asymmetric nature of the HK kinase-active state. Based on our mutagenesis, structural, and biochemical data, we propose a mechanical model in which propagation of the conformational signal through the HAMP domain can regulate autokinase activity by inducing segmental movements and pronounced bending of the central HK α-helices, which in turn control the mobility of the ATP-binding catalytic domains.
Overall Structure of the CpxA Transmitter Core
The full cytoplasmic region of E. coli CpxA (CpxAHDC), encompassing the HAMP signaling domain (residues 188–237) and the DHp and CA domains of the catalytic core (residues 238–457), has been produced as recombinant protein for structural studies (Figure 1A). The 3D structures of wild-type CpxAHDC (or the point mutant CpxAHDC_M228V) in complex with nucleotides have been determined in five distinct crystalline environments (Table 1) using either single-wavelength anomalous diffraction (SAD) or molecular replacement methods (see Materials and Methods). The CpxAHDC structure consists of a homodimer, in which the HAMP and DHp domains from both subunits assemble into an elongated central coiled-coil region flanked on either side by the two CA domains (Figure 1B).
Figure 1
Figure 1
Overall structure of CpxA.
Table 1
Table 1
Crystallographic data collection and refinement statistics.
The polypeptide chains of the HAMP domain were fully traced into SAD-phased experimental electron density maps of the CpxAHDC orthorhombic and trigonal crystal forms (Figure S1). The domains are formed by two parallel helices (α1 and N-terminal part of α2) linked by a long connector (Figure 1B), which despite a low sequence identity adopt a dimeric four-helix bundle similar to the Af1503 and VicK HAMP structures (Figure S2) [15],[28][30]. In some crystal forms, however, only the C-terminal α2 helix from both monomers could be traced in the electron density maps, because CpxAHDC has a tendency to establish interactions via its central-core helices, leading to the formation of dimers-of-dimers that sterically interfere with the positioning of helix α1 (Figure S3A). Indeed, previous biophysical studies of full-length CpxA reconstituted in nanodiscs had indicated that the protein could exist as a mixture of dimers and tetramers (Subrini, PhD thesis, 2011). We carried out SAXS studies of CpxAHDC at different protein concentrations, which also demonstrate that CpxAHDC is prone to associate into dimers-of-dimers, as scattering data suggest a concentration-dependent equilibrium between dimers and tetramers in solution (Figure S3). Guinier analysis shows no significant nonspecific protein aggregation, albeit the Rg dependence with protein concentration points to interparticle interactions that would likely favor the tetrameric state at protein concentrations similar to those used for the crystallization trials (>10 mg/ml). However, extrapolation of these data to very low concentration suggests that CpxA homodimers are the functionally relevant form of the protein, given the low protein abundance at the E. coli inner membrane.
The DHp antiparallel two-helix hairpins (helices α2 and α3) from both protomers also assemble into a four-helix bundle (Figure 1B) as observed in other HKs structures [10],[11],[31]. Both the HAMP and DHp coiled coils display the characteristic heptad repeats periodicity, in which the inner heptad positions a and d are occupied by hydrophobic residues. Besides the presence of a heptad repeat stutter at Leu232 [32], HAMP and DHp domains merge into a long continuous helix (α2) (Figure 1B), in a very similar way as previously reported for the chimeric fusion between AF1053 HAMP and DHp EnvZ [28]. At the C-terminus of each monomer, DHp helix α3 connects through a short flexible linker (residues 300–304) to the CA domain (Figure 1B). These catalytic domains adopt the ATP-binding Bergerat fold [33], a two-layer 3α/5β sandwich fold typical of the GHKL (bacterial gyrase, HSP90, histidine kinase, MutL) superfamily. In HKs, this topology is characterized by four conserved sequence motifs termed G1, G2, F, and N boxes (Figure S4) that include residues engaged in direct contacts with the bound ATP or confer flexibility to a flexible lid loop that surround the nucleotide binding cleft.
Conformational and Dynamical Asymmetry of the CpxA Catalytic Core
The 3D structure of CpxAHDC is highly asymmetric, not only by a different orientation of the CA domains, but also as a consequence of two potential bending points in the central helix α2 (Figure 1B). The first kink at Ser238 is located at the boundary between the HAMP and DHp domains, where there is a discontinuity between the parallel (HAMP) and the antiparallel (DHp) four-helix bundles. Helical bending at this position is probably required to avoid steric clashes at a layer of polar residues (Gln239–Gln240). The second kink is caused by the presence of a highly conserved proline residue (Pro253) in the proximity of the phosphorylatable His248, and has been also observed in other HK structures, such as HK853 [11], KinB [34], and VicK [15]. As a consequence of different bending angles in the two central α3 helices, the overall structure of the homodimer displays a curved asymmetric conformation (Figure S5), similar to that reported for the wild-type Af1503 HAMP–EnvZ DHp chimera [28]. It is interesting to note that the presence of two-helix–disrupting centers within α2 implies segmental helical movements in response to any torsional or orientational stress signal transmitted through the HAMP domain, which could easily result in asymmetric helical bending of the DHp domain.
When considering together all CpxAHDC crystal structures, it transpires that the DHp helical bending promotes a strikingly different dynamic behavior of the two CA domains in the homodimer, primarily due to distinct interactions with the central four-helix bundle. Thus, whereas one CA domain (from subunit A) is retained in a fixed—inactive—conformation close to the DHp helices (α2 from subunit B and α3 from subunit A) through a relatively large contact interface (~1,200 Å2), the other CA domain with apparent unrestrained mobility adopts disparate orientations (Figure 2A) and displays significantly smaller CA–DHp surface contact areas, the largest of which (~800 Å2) corresponds to the autophosphorylating Michaelis complex (see below). This strongly suggests a highly dynamic autophosphorylation process, leading primarily to a single phosphorylation event per homodimer, a feature already described for various HKs [10],[35],[36]. In fact, in vitro autokinase activity assays show that CpxAHDC phosphorylation is never complete (Figure 2B), with only ~70% of the protein phosphorylated at steady state (Figure 2C). These results suggest that CpxA could be mainly hemiphosphorylated in vivo under stress conditions, as the higher in vitro autophosphorylation efficiency might be attributed to subunit exchange during the assay, as observed for EnvZ [37].
Figure 2
Figure 2
Conformational and dynamical asymmetry of the homodimer.
The Autophosphorylating Michaelis Complex
In all homodimers, the fixed CA domain is sequestered by the DHp domain in an inactive kinase conformation, with the ATP-γP positioned far away from the phosphorylatable His248. Instead, the active site formed by the second, mobile CA domain is poised for catalysis in three independent structures crystallized in two distinct crystal forms (trigonal and hexagonal space groups, Table 1), with the ATP γ-P and the acceptor His248 Nε atom well positioned for the phosphotransfer reaction to proceed (Figure 3A). The catalytic cores (DHp + CA domains) are very similar to each other in the three independent crystal structures (Figure S6); they can be superimposed with r.m.s.d. of 0.9–1.4 Å and display an identical network of hydrogen bonding interactions involving key functional residues for ATP binding and catalysis. This highly conserved active site geometry in distinct crystalline environments emphasizes the stability of the CA–DHp interface in the Michaelis complex and lends further support to its functional relevance.
Figure 3
Figure 3
Autophosphorylating Michaelis complex.
Most residues interacting with the nucleotide belong to the CA domain (Figure S7), except for the phosphoacceptor His248 and the conserved Arg251, both of which make strong hydrogen bonding interactions with oxygen atoms from the ATP γ-P. The ATP phosphate groups are further stabilized in the binding pocket at the correct position for catalysis through additional contacts made by their oxygen atoms with the backbone amide nitrogen atoms of Leu419, Gly420, and Leu421, the hydroxyl groups of Thr417 and Tyr364, the guanidinium group of Arg363, and the carboxamide group of Asn360 (Figure 3A and Figure S7). Most of these residues belong to conserved sequence motifs in the HK family [38]; in particular, Asn360 is a strictly conserved residue from the N-box in the CA domain that is essential for the autophosphorylation reaction [39],[40].
The geometry of the active site points to Glu249, the amino acid adjacent to the phosphoacceptor His248 in the DHp domain, as a key functional residue. Its carboxylic side chain makes hydrogen-bonding interactions with both the Nδ atom of His248 and the carboxamide group of Asn356 from the CA domain (Figure 3A). These interactions favor the deprotonation of the His248 Nε atom, which performs the nucleophilic attack on the ATP γ-P (Figure 3B). Consistent with such a role, mutations of this glutamate were found to abolish kinase activity in other HKs such as NtrB [41] and CdrS [42], and an acidic residue at this position (immediately following the phosphoacceptor histidine) is highly conserved in major HK subfamilies HisKA (Pfam00512) and HisKA_3 (Pfam07730), which represents about 90% of all known HK sequences. Further evidence supporting a key functional role of this residue is provided by the structural comparison of the CpxA active site with that of the chemotaxis HK CheA [43]. Although the two proteins have a different topology, a CheA glutamic acid (Glu67) is spatially equivalent to CpxA Glu249 and interacts in an analogous way with the phosphoacceptor His45 (Figure S8), strongly suggesting a similar functional role for the acidic residue in the autophosphorylation mechanism.
Also important is the interaction between Glu249 and Asn356 across the DHp–CA interface (Figure 3A), which could help to achieve the correct disposition and polarization of Glu249 for catalysis. In agreement with this hypothesis, we found that the substitution Asn356–Tyr presented a strong kinase-deficient phenotype (Table 2) in a genetic screen designed to search for mutations that affect kinase activity in vivo (see below). It is interesting to note that a conserved glutamate at the equivalent position of CpxA Asn356 acts as a general base in the catalytic reaction of ATPases from the GHKL superfamily, which share a similar ATP-binding domain with HKs [44],[45].
Table 2
Table 2
Identification of intragenic suppressors.
In the Michaelis complex, the CA carrying the ATP molecule and the phosphorylatable His belong to different subunits of the homodimer, which structurally confirms trans-autophosphorylation for CpxA (Figure 3C). However, it is worth noting that trans-phosphorylation does not seem to be a generic feature of HKs; recent work by Ashenberg and coworkers [46] has shown that the cis/trans character of the reaction depends not on the local structure of the active site, but on the handedness of the DHp domain, which is determined by the loop at the base of the DHp four-helix bundle (i.e., between helices α2 and α3 in CpxA; Figure 1). This model is now further confirmed by the structural comparison of CpxA (trans-autotophosphorylation) and VicK (cis-autophosphorylation, [15]). Although both proteins have opposite handedness of the central four-helix bundle, their structures revealed a very similar active CA–DHp interaction (Figure S9).
In addition to the polar interactions described above, two hydrophobic residues from the CA domain (Phe403 in the F box and Leu419 in the G2 box) play a key dual role in stabilizing the DHp–CA interface for both the active and inactive CA conformations (Figure 3D). In the active conformation, the side chain of Leu419 is in contact with the phosphorylatable His248, while Phe403 interacts with Leu292 and the aliphatic moiety of Arg296 in DHp helix α3 from the opposite monomer. In the second, inactive, subunit the same two residues participate in anchoring the CA domain to a hydrophobic binding pocket mainly defined by residues Met287, Leu291, and Met 294 from helix α3 (in the same monomer) and residues Leu242, Ile246, and Leu250 from helix α2 (in the opposite monomer), accounting for a larger DHp–CA domain interface. Interestingly, some of these residues correspond to the a and d core positions of the DHp coiled-coil, which become partially accessible in one of the monomers due to the asymmetric conformation (see below).
These differences in the hydrophobic anchoring of Phe403 and Leu419 to the DHp domain (Figure 3D) might also contribute to modulate ATP binding affinity during the catalytic cycle. We have observed that the CpxAHDC homodimer exhibited half-occupation of the active site in some crystal forms, and in these cases, only the less mobile CA domains (i.e., those sequestered by stronger hydrophobic interactions with the DHp domain) consistently displayed a bound nucleotide covered by the lid loop (Figure S10), whereas the CA domains with higher mobility have an empty nucleotide-binding cleft. These observations suggest that proper positioning of the lid loop through specific interactions with the DHp domain might trigger ATP binding. This hypothesis is consistent with recent observations from Mycobacterium tuberculosis DosS [47], which is unable to bind ATP until the CA domain is properly positioned for catalysis, and is reminiscent of previous findings on CheA dimers [48], where the two sites have largely different affinities for ATP. In our case, a lower nucleotide substrate affinity for the floppy CA domains might facilitate nucleotide turnover during the HK activation cycle.
The HAMP Domain and CpxA Signal Transduction
Despite the overall asymmetry of the CpxA homodimer, the HAMP four-helix bundle is perfectly symmetric as was also observed in the Af1503, Aer2, and VicK structures [15],[28][30],[49],[50]. A strict 2-fold axis passing along the length of the coiled-coil relates the two HAMP subunits (Figure 4A), which adopts a packing mode closer to a canonical knobs-into-holes model [51] than the unusual complementary x-da observed in the wild-type Af1503 structures [28][30]. However, a negative slope in the Crick angle deviation plot [29],[52] suggests a departure from an ideal knobs-into-holes packing mode due to a weak heptad repeat periodicity at α1 (Figure 4B). Overall, the observed CpxA–HAMP conformation resembles that observed in several Af1503 HAMP mutants associated with a kinase-active state [28][30].
Figure 4
Figure 4
HAMP structure and function.
The similarities in helical packing modes and overall helical bending between CpxA and the Af1503 HAMP–EnvZ DHp fusion chimera (Figure S5) prompted us to investigate whether CpxA follows an analogous two-state mechanism of signal transduction. Based on studies of Af1503, Hulko et al. [30] proposed a model for HAMP-mediated signal transduction, in which the HAMP coiled-coil can adopt two distinct conformations depending on the input signal. In particular, they found that the single point mutation (Ala291–Val) in Af1503 displaced the equilibrium between the two forms [30] and partially reversed the helical bending asymmetry of the wild-type Af1503–EnvZ crystal structure [28]. We therefore decided to introduce the equivalent mutation in CpxA (Ala197–Val; see Figure S2). However, the responses to extracytoplasmic stress due to the overproduction of a folding-defective mutant of MalE in the periplasm (MalE31, [53]), or to the presence of phenethyl alcohol in the inner membrane [54], were similar in bacteria expressing either CpxAA197V or wild-type CpxA (Figure 4C).
To gain further insight into the signaling state(s) of the HAMP domain, we set up a genetic screening to search for mutations that affect the CpxA kinase activity in vivo. For this, a plasmid encoding a CpxA variant, which lacks the periplasmic sensor domain (CpxAΔP) and therefore displays constitutive kinase activity, was randomly mutagenized and transformed into a cpxA-null strain background. In these cells, the high β-galactosidase activity of the transcriptional fusion (correlated to the kinase activity) provided a means to isolate intragenic suppressors that confer a kinase-deficient/phosphatase-dominant (K/P+; i.e., Lac) phenotype, easily monitored on X-gal containing plates. From 48 Lac colonies examined, the identified mutations were reconstructed into the wild-type cpxA gene and their in vivo activity was evaluated (Table 2). The strongest phenotype corresponded to the point mutation Met228–Val in the HAMP, which completely abolished the Cpx-system response (Figure 4C). In vitro studies revealed that the point mutant CpxAHDC_M228V is less efficient than wild-type CpxA in phosphorylating CpxR (Figure 4D), at least in part as a consequence of a reduced autophosphorylation rate (Figure 4E). Because Met228 occupies a heptad-repeat inner position in the four-helix bundle, its substitution by a smaller valine side-chain would promote a rearrangement of the interhelical packing to avoid the creation of a cavity within the hydrophobic core (Figure 4A). Consistent with a more stable four-helix bundle packing, our SAXS data suggest that the Met228–Val substitution decreases the tendency of CpxAHDC to assemble into tetramers (dimers of dimers), which partially disrupts the overall HAMP structure. Thus, in contrast to wild-type CpxAHDC, CpxAHDC_M228V is mainly dimeric in solution, with a relatively constant Rg value at different protein concentrations (Table S1 and Figure S3D).
A Mechanical Model to Account for Signal Transduction and Autophosphorylation
As described above, our biochemical and structural data indicate that HAMP functioning is consistent with a two-state model for signal transduction [49],[55], in which different types of stress (torsional, rotational, or translational) coming from the membrane would break the conformational symmetry of the homodimer by inducing pronounced helical bending of the full CpxA transmitter core and a strong dynamical asymmetry in CA mobility (Figure 5). An important issue here is that, due to the presence of the helix-disrupting points at Ser238 and Pro253 in the long α2 helix, the overall bending movement would arise from segmental helical mobility. Thus, the movements of the central part of helix α2 (between Ser238 and Pro253, roughly corresponding to the H-box) in each monomer differ from those of the HAMP domain (i.e., α1 and N-terminus of α2 in both monomers) and the membrane-distal part of the DHp domain (i.e., C-terminus of α2 and α3 in both monomers), both of which rotate as separate rigid bodies. As a consequence, in the asymmetric (kinase-active) state, the C-terminus of α3 in one monomer moves away from the central part of α2 in the other monomer, thus partially exposing hydrophobic core residues (Figure 5, insert). This exposure creates a binding interface that retains one CA domain in a fixed inactive conformation (Figure 3D). Instead, on the opposite face of the DHp domain, the same bending movements release the second CA domain, which becomes available to form a transient active site configuration poised for catalysis (Figure 3A). This mechanical model highlights the importance of sequestering the CA domain (through interactions with the DHp domain) as a major mechanism to switch off kinase activity [28], and possibly to bring the catalytic core into phosphatase-active (symmetric) or phosphotransferase-active (asymmetric) states. Taken together, our results point to the importance of conformational and dynamical asymmetry in modulating the autophosphorylation reaction in HKs, in much the same way as ligand-induced asymmetry of homodimeric HK sensor domains was found to play a major role in two-component signal transduction [8],[9].
Figure 5
Figure 5
Mechanical model for HK autophosphorylation control.
Cloning, Expression, and Purification
CpxAHDC construct was cloned into the expression vector pET28a+ (Novagen) and expressed in E. coli BLI5 cells. His-tagged CpxAHDC was purified by immobilized metal affinity chromatography followed by gel filtration using standard procedures. For the production of seleno-methionyl (SeMet) CpxAHDC, the same plasmid was transformed into the autotrophic E. coli B834 (DE3) strain. Differing only in the growth media, SelenoMet Medium (MDL), and L-selenomethionine (50 mg/L) instead of LB, SetMet proteins were purified exactly as the unlabeled ones. Mass spectrometry was used to confirm SeMet incorporation. Full-length CpxA was expressed, solubilized, and purified as described previously [56].
Crystallization and Data Collection
Purified proteins were concentrated using ultrafiltration VivaSpin devices (30 kDa MWCO, Millipore) to ~15 mg/ml prior crystallization trials. Crystals were obtained by vapor diffusion at 18°C. All crystal forms grew under similar crystallization conditions [100 mM Tris-HCl (pH 8.5), 1.75 M Ammonium sulphate, and 25% (v/v) glycerol], with the exception of the two monoclinic forms, which were obtained using a solution containing 100 mM Tris-HCl (pH 8.5), 25%(w/v) PEG3350, and 200 mM lithium sulphate as precipitant. In most of the cases, crystals appeared in a few days and reached their maximum size within 2 wk.
Data were collected at ESRF and Soleil synchrotron beamlines (Table 1). Data were indexed and integrated with program XDS [57], followed by space group determination with Pointless, and scaling with Aimless from the CCP4 software package [58].
Structure Determination and Refinement
Diffraction data from Se-labeled protein crystals were used to solve the structures of the trigonal (P3121) and hexagonal (P6122) crystal forms. Se atom substructure and phases were estimated using the AutoSol wizard of PHENIX [59]. After density modification, the resulting electron density maps were readily interpretable for most of the protein chain. Initial models were manually built with Coot [60] by fitting into the experimental map template helices and a homology model of the CA domain. The orthorhombic crystal form (C2221) was solved by SAD from a Pt derivative dataset at 4.1 Å resolution. Heavy atom positions and phases were estimated using truncated data at 6.5 Å resolution (the anomalous measurability limit) with the AutoSol wizard. Despite the poor quality of the experimental map, it was possible to locate the CA domains and build a main-chain trace for the rest of the protein. Given the nonisomorphism between the Pt-derivative and the native dataset at 3.65 Å resolution, an extensive rigid body refinement was necessary to reorient the three dimers present in the native ASU prior to further refinement. Monoclinic crystal forms were determined by molecular replacement with Phaser [61] using the previously solved structures as search models. All models were refined with Buster [62] including as restrains NCS, TLS, a reference model, and experimental phases when available. Final models were validated with MolProbity [63]. Data collection and refinement statistics are shown in Table 1. Figures were generated and rendered with PyMOL [64].
Small Angle X-Ray Scattering
SAXS data were collected at 20°C at beamline BM29 at the ESRF. Data reduction and analysis following the standard procedures were performed using the ATSAS program package [65]. All samples were in 25 mM Tris-HCl buffer, pH 8.5, containing 300 mM NaCl, 5 mM MgCl2, 5%(v/v) glycerol, and 0.5 mM TECP.
Screening for Phosphatase-Dominant (KP+) Mutations
Plasmid pLCBA was constructed by subcloning cpxA from pIV3cpxA [49] into the low-copy pLCB vector [66]. This plasmid expressed cpxA about 1.5-fold in strains grown in LB at 37°C in the absence of arabinose (Betton, unpublished results). The deletion of the periplasmic sensor domain of CpxA was made by the one-step site-directed deletion method [67] using pLCBA as DNA template, and the primers 5′-ATCTAGATTCACCGCTATTACTGCTGATTGTCACCATG-3′ and 5′-TAGATCTAACCCATTGTAGTTTTGGTTGTAGTTGTG-3′. The PCR fragment was first treated with DpnI, and then with XbaI before ligation and transformation. The resulting plasmid, pLCBAΔP, encoded a CpxA variant in which the residues 29–163 are substituted by NLDS. Error-prone PCR was performed using the GeneMorph II Random Mutagenesis kit (Agilent Technologies) according to the manufacturer's instructions. Amplification using about 325 ng of cpxAΔP gene as DNA template and 20 thermal cycles was performed to aim for 1~3 mutations/kb. The purified PCR fragment was digested, ligated into pLCB, and transformed into NS54 E. coli strain. Approximately 20,000 transformants were platted on LB-agar containing X-gal and 48 white colonies were picked for confirming the linkage of the Lac phenotype to the plasmid. After DNA sequencing, single suppressor mutations were moved by site-directed mutagenesis on the wild-type cpxA gene cloned into pLCBA. To validate the genetic analysis, β-galactosidase activities conferred by the single mutant cpxA alleles were assayed in vivo in NS54 strain.
β-Galactosidase Assays
NS54 cells harboring the various pLCB derivatives were grown overnight in LB medium containing chloramphenicol at 37°C, subcultured (1[ratio]100) into 5 ml of the same medium, and then grown to midlog phase at 37°C. The β-galactosidase activity of permeabilized cells was assayed as described previously [53]. A minimum of four independent determinations was averaged to obtain the indicated values.
Autophosphorylation and Phosphotransfer Assays
Autokinase activity on soluble CpxAHDC and phosphotransfer from full-length CpxA to CpxRN (receiver domain) were assayed by Phos-tag acrylamide gel electrophoresis. Autokinase activity assays on full-length CpxA (as Brij35 complexes) were performed using radioactive ATP as described [56]. For phosphotransfer assays, full-length CpxA (10 µM) in buffer H (25 mM HEPES, pH 7.8, containing 100 mM NaCl, 50 mM KCl, 0.05% Brij35, and 5 mM MgCl2) was phosphorylated with 100 µM ATP at 25°C for 20 min. Then, purified CpxRN was added to the reaction to a final concentration of 20 µM. Samples were withdrawn at various time intervals, mixed with SDS sample buffer, and kept on ice until all reactions were completed. Phospho-proteins were separated by Phos-tag acrylamide gel electrophoresis. Polyacrylamide gels were polymerized with 25 µM Phos-tag acrylamide and 50 µM MnCl2. Following electrophoresis, gels were stained with Coomassie blue and analyzed by densitometry using the software Quantity One (Biorad). Data were fitted using a nonlinear regression with program Kaleidagraph (Synergy Software).
Accession Numbers
The crystal structures of CpxA-nucleotide complexes described in this article have been deposited in the Protein Data Bank under accession numbers 4BIU (CpxAHDC–ADP, orthorhombic C2221 form), 4BIV (CpxAHDC–ATP, trigonal P3121 form), 4BIW (CpxAHDC_M228V–AMPPNP, hexagonal P6122 form), 4BC0 (CpxAHDC–ATP, hexagonal P6122 form), 4BIX (CpxAHDC_M228V–ADP, monoclinic C2 form), and 4BIY (CpxAHDC_M228V–ADP, monoclinic I2 form).
Figure S1
(A) The trigonal crystal form was solved by SAD using Se-labeled protein crystals. Anomalous difference Fourier electron density map (yellow mesh) is contoured at the 5 σ level. Selenomethionine residues are shown as sticks. (B) MR-SAD density modified electron density map (blue mesh) calculated using structure factor anomalous differences and a partial model omitting the HAMP (residues 188–237). Final model is superposed to this unbiased electron density map contoured at the 1.5 σ level.
(TIF)
Figure S2
Structural superposition of the CpxA HAMP domain with those of Af1503 [29] and VicK [15]. The two mutated residues in CpxA, Ala197 (equivalent to Af1503 Ala291) and Met228, are indicated.
(TIF)
Figure S3
(A) Tetrameric form of CpxA (dimer of dimers) as observed in the hexagonal and monoclinic crystal forms. Formation of these tetramers was found to sterically interfere with the positioning of HAMP helix α1. (B) Experimental SAXS curves for wild-type CpxAHDC at different protein concentrations in the presence of 5 mM ADP. (C) Same for the point mutant CpxAHDC_M228V. (D) Linear dependence of the radius of gyration (Rg) as a function of protein concentration for CpxAHDC (green circles) and CpxAHDC_M228V (blue circles).
(TIF)
Figure S4
Cartoon representation of the ATP-binding CA domain. Highlighted in colors are the highly conserved sequence motifs (G1, G2, N, and F boxes) in the GHKL superfamily.
(TIF)
Figure S5
A similar helical core bending was observed in the crystal structures of CpxAHDC (left) and wild-type Af1503 HAMP–EnvZ DHp chimera (right).
(TIF)
Figure S6
Structural superposition reveals that the conformation of the catalytic core of wild-type CpxAHDC in complex with ATP (yellow) crystallized in the trigonal space group P3121 is very similar to that observed in the hexagonal (P6122) crystal form for either the wild-type CpxAHDC–ATP (pink) or the point mutant CpxAHDC_M228V–AMPPNP (green) complexes. The catalytic core (residues 238–455) of CpxAHDC–ATP (trigonal form) and CpxAHDC_M228V–AMPPNP (hexagonal form) crystal structures can be superposed with an RMSD of 1.36 Å. More significant changes are observed, however, for the N-terminal helix of the HAMP domain due to different crystal packing environments.
(TIF)
Figure S7
Protein–nucleotide hydrogen-bonding interactions as observed in the crystal structure of CpxAHDC in complex with AMPPNP (hexagonal crystal form).
(TIF)
Figure S8
The structural superposition of CpxA DHp domain (in green) and CheA P1 domain (PDB ID 1TQG, in yellow) reveals a spatially equivalent glutamic acid that makes hydrogen-bonding interactions with Nδ of the catalytic histidine and serves to enhance nucleophilic catalysis. However, the two proteins display a different topology: whereas the two residues in CpxA are part of the conserved H-box in helix α2, they belong to different helices in CheA.
(TIF)
Figure S9
Cartoon and schematic representations of the DHp four-helix bundle assembly and CA domain positioning in (A) trans-autophosphorylating CpxAHDC (hexagonal crystal form) and (B) cis-autophosphorylating VicK (4I5S) crystal structures.
(TIF)
Figure S10
Half occupation of the ATP binding sites. Cartoon representation of the single occupied ADP-binding site in the CpxAHDC_M228V homodimer crystallized in the monoclinic C2 crystal form (the second monomer in this crystal structure has no bound nucleotide). The green mesh corresponds to a 2 Å resolution σA weighted difference electron density map (mFo-DFc) contoured at the 3 σ level. The ADP molecule and the Mg2+ ion (shown as ball-and-sticks) were omitted from the model before map calculation. Also shown are residues Phe403 and Leu419 engaged in the CA–DHp interface. It can be argued that the presence of ADP (instead of ATP or AMPPNP) in the crystals might account for the observed half-occupancy. Against this hypothesis, however, the structure of CpxAHDC in complex with AMPPNP (determined at 2.7 Å resolution in the same monoclinic C2 crystal form; unpublished data) showed the same half-occupancy pattern.
(TIF)
Table S1
CpxAHDC and CpxAHDC_M228V SAXS-derived parameters.
(DOCX)
Acknowledgments
We acknowledge A. Haouz and P. Weber (PF6, Institut Pasteur) for performing robot-driven crystallization trials, the ESRF and Synchrotron Soleil for provision of synchrotron radiation facilities, and their staffs for the helpful assistance.
Abbreviations
AMP-PNPAdenosine 5′-(β,γ-imido)triphosphate
CAcatalytic and ATP binding domain; DHp, dimerization and histidine phosphorylation
HAMPHistidine kinase, Adenyl cyclase, Methyl-accepting proteins, Phosphatase
HKhistidine kinase
RRresponse regulator
TCStwo-component system

Funding Statement
This work has been partly supported by grants from the Institut Pasteur and the Agence Nationale de la Recherche (ANR, France). AEM was recipient of a DK fellowship from the Basque Government. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Stock AM, Robinson VL, Goudreau PN (2000) Two-component signal transduction. Annu Rev Biochem 69: 183–215. [PubMed]
2. Capra EJ, Laub MT (2012) Evolution of two-component signal transduction systems. Annu Rev Microbiol 66: 325–347. [PubMed]
3. Gao R, Stock AM (2009) Biological insights from structures of two-component proteins. Annu Rev Microbiol 63: 133–154. [PMC free article] [PubMed]
4. Szurmant H, White RA, Hoch JA (2007) Sensor complexes regulating two-component signal transduction. Curr Opin Struct Biol 17: 706–715. [PMC free article] [PubMed]
5. Cheung J, Hendrickson WA (2010) Sensor domains of two-component regulatory systems. Curr Opin Microbiol 13: 116–123. [PMC free article] [PubMed]
6. Aravind L, Ponting CP (1999) The cytoplasmic helical linker domain of receptor histidine kinase and methyl-accepting proteins is common to many prokaryotic signalling proteins. FEMS Microbiol Lett 176: 111–116. [PubMed]
7. Moore JO, Hendrickson WA (2009) Structural analysis of sensor domains from the TMAO-responsive histidine kinase receptor TorS. Structure 17: 1195–1204. [PubMed]
8. Moore JO, Hendrickson WA (2012) An asymmetry-to-symmetry switch in signal transmission by the histidine kinase receptor for TMAO. Structure 20: 729–741. [PMC free article] [PubMed]
9. Neiditch MB, Federle MJ, Pompeani AJ, Kelly RC, Swem DL, et al. (2006) Ligand-induced asymmetry in histidine sensor kinase complex regulates quorum sensing. Cell 126: 1095–1108. [PMC free article] [PubMed]
10. Albanesi D, Martín M, Trajtenberg F, Mansilla MC, Haouz A, et al. (2009) Structural plasticity and catalysis regulation of a thermosensor histidine kinase. Proc Natl Acad Sci U S A 106: 16185–16190. [PubMed]
11. Marina A, Waldburger CD, Hendrickson WA (2005) Structure of the entire cytoplasmic portion of a sensor histidine-kinase protein. EMBO J 24: 4247–4259. [PubMed]
12. Casino P, Rubio V, Marina A (2009) Structural insight into partner specificity and phosphoryl transfer in two-component signal transduction. Cell 139: 325–336. [PubMed]
13. Yamada S, Sugimoto H, Kobayashi M, Ohno A, Nakamura H, et al. (2009) Structure of PAS-linked histidine kinase and the response regulator complex. Structure 17: 1333–1344. [PubMed]
14. Dago AE, Schug A, Procaccini A, Hoch JA, Weigt M, et al. (2012) Structural basis of histidine kinase autophosphorylation deduced by integrating genomics, molecular dynamics, and mutagenesis. Proc Natl Acad Sci U S A 109: E1733–E1742. [PubMed]
15. Wang C, Sang JY, Wang JW, Su MY, Downey JS, et al. (2013) Mechanistic insights revealed by the crystal structure of a histidine kinase with signal transducer and sensor domains. PLoS Biol 11: e1001493. [PMC free article] [PubMed]
16. Hunke S, Keller R, Müller VS (2011) Signal integration by the Cpx-envelope stress system. FEMS Microbiol Lett 326: 12–22. [PubMed]
17. Vogt SL, Raivio TL (2011) Just scratching the surface: an expanding view of the Cpx envelope stress response. FEMS Microbiol Lett 326: 2–11. [PubMed]
18. Raivio TL, Silhavy TJ (1999) The sigmaE and Cpx regulatory pathways: overlapping but distinct envelope stress responses. Curr Opin Microbiol 2: 159–165. [PubMed]
19. De Wulf P, McGuire AM, Liu X, Lin EC (2002) Genome-wide profiling of promoter recognition by the two-component response regulator CpxR-P in Escherichia coli. J Biol Chem 277: 26652–26661. [PubMed]
20. Danese PN, Snyder WB, Cosma CL, Davis LJ, Silhavy TJ (1995) The Cpx two-component signal transduction pathway of Escherichia coli regulates transcription of the gene specifying the stress-inducible periplasmic protease, DegP. Genes Dev 9: 387–398. [PubMed]
21. Pogliano J, Lynch AS, Belin D, Lin EC, Beckwith J (1997) Regulation of Escherichia coli cell envelope proteins involved in protein folding and degradation by the Cpx two-component system. Genes Dev 11: 1169–1182. [PubMed]
22. Price NL, Raivio TL (2009) Characterization of the Cpx regulon in Escherichia coli strain MC4100. J Bacteriol 191: 1798–1815. [PMC free article] [PubMed]
23. Raivio TL, Laird MW, Joly JC, Silhavy TJ (2000) Tethering of CpxP to the inner membrane prevents spheroplast induction of the cpx envelope stress response. Mol Microbiol 37: 1186–1197. [PubMed]
24. Thede GL, Arthur DC, Edwards RA, Buelow DR, Wong JL, et al. (2011) Structure of the periplasmic stress response protein CpxP. J Bacteriol 193: 2149–2157. [PMC free article] [PubMed]
25. Zhou X, Keller R, Volkmer R, Krauss N, Scheerer P, et al. (2011) Structural basis for two-component system inhibition and pilus sensing by the auxiliary CpxP protein. J Biol Chem 286: 9805–9814. [PubMed]
26. Raivio TL, Popkin DL, Silhavy TJ (1999) The Cpx envelope stress response is controlled by amplification and feedback inhibition. J Bacteriol 181: 5263–5272. [PMC free article] [PubMed]
27. Fleischer R, Heermann R, Jung K, Hunke S (2007) Purification, reconstitution, and characterization of the CpxRAP envelope stress system of Escherichia coli. J Biol Chem 282: 8583–8593. [PubMed]
28. Ferris HU, Dunin-Horkawicz S, Hornig N, Hulko M, Martin J, et al. (2012) Mechanism of regulation of receptor histidine kinases. Structure 20: 56–66. [PubMed]
29. Ferris HU, Dunin-Horkawicz S, Mondéjar LG, Hulko M, Hantke K, et al. (2011) The mechanisms of HAMP-mediated signaling in transmembrane receptors. Structure 19: 378–385. [PubMed]
30. Hulko M, Berndt F, Gruber M, Linder JU, Truffault V, et al. (2006) The HAMP domain structure implies helix rotation in transmembrane signaling. Cell 126: 929–940. [PubMed]
31. Inouye M, Ikura M, Tomomori C, Tanaka T, Dutta R, et al. (1999) Solution structure of the homodimeric core domain of Escherichia coli histidine kinase EnvZ. Nat Struct Biol 6: 729–734. [PubMed]
32. Stewart V, Chen LL (2010) The S helix mediates signal transmission as a HAMP domain coiled-coil extension in the NarX nitrate sensor from Escherichia coli K-12. J Bacteriol 192: 734–745. [PMC free article] [PubMed]
33. Dutta R, Qin L, Inouye M (1999) Histidine kinases: diversity of domain organization. Mol Microbiol 34: 633–640. [PubMed]
34. Bick MJ, Lamour V, Rajashankar KR, Gordiyenko Y, Robinson CV, et al. (2009) How to switch off a histidine kinase: crystal structure of Geobacillus stearothermophilus KinB with the inhibitor Sda. J Mol Biol 386: 163–177. [PMC free article] [PubMed]
35. Jiang P, Peliska JA, Ninfa AJ (2000) Asymmetry in the autophosphorylation of the two-component regulatory system transmitter protein Nitrogen Regulator II of Escherichia coli. Biochemistry 39: 5057–5065. [PubMed]
36. Trajtenberg F, Grana M, Ruetalo N, Botti H, Buschiazzo A (2010) Structural and enzymatic insights into the ATP binding and autophosphorylation mechanism of a sensor histidine kinase. J Biol Chem 285: 24892–24903. [PubMed]
37. Cai S-J, Inouye M (2003) Spontaneous subunit exchange and biochemical evidence for trans-autophosphorylation in a dimer of Escherichia coli histidine kinase (EnvZ). J Mol Biol 329: 495–503. [PubMed]
38. Wolanin PM, Thomason PA, Stock JB (2002) Histidine protein kinases: key signal transducers outside the animal kingdom. Genome Biol 3: REVIEWS3013. [PMC free article] [PubMed]
39. Hsing W, Russo FD, Bernd KK, Silhavy TJ (1998) Mutations that alter the kinase and phosphatase activities of the two-component sensor EnvZ. J Bacteriol 180: 4538–4546. [PMC free article] [PubMed]
40. Bilwes AM, Quezada CM, Croal LR, Crane BR, Simon MI (2001) Nucleotide binding by the histidine kinase CheA. Nat Struct Biol 8: 353–360. [PubMed]
41. Atkinson MR, Ninfa AJ (1993) Mutational analysis of the bacterial signal-transducing protein kinase/phosphatase nitrogen regulator II (NRII or NtrB). J Bacteriol 175: 7016–7023. [PMC free article] [PubMed]
42. Willett JW, Kirby JR (2012) Genetic and biochemical dissection of a HisKA domain identifies residues required exclusively for kinase and phosphatase activities. PLoS Genet 8: e1003084. [PMC free article] [PubMed]
43. Quezada CM, Hamel DJ, Gradinaru C, Bilwes AM, Dahlquist FW, Crane BR, Simon MI (2005) Structural and chemical requirements for histidine phosphorylation by the chemotaxis kinase CheA. J Biol Chem 280: 30581–30585. [PubMed]
44. Dutta R, Inouye M (2000) GHKL, an emergent ATPase/kinase superfamily. TIBS Biochem Sci 251: 24–28. [PubMed]
45. Corbett KD, Berger JM (2005) Structural dissection of ATP turnover in the prototypical GHL ATPase TopoVI. Structure 13: 873–882. [PubMed]
46. Ashenberg O, Keating AE, Laub MT (2013) Helix bundle loops determine whether histidine kinases autophosphorylate in cis or in trans. J Mol Biol 425: 1198–1209. [PMC free article] [PubMed]
47. Cho HY, Lee YH, Bae YS, Kim E, Kang BS (2013) Activation of ATP binding for the autophosphorylation of DosS, a Mycobacterium tuberculosis histidine kinase lacking an ATP lid motif. J Biol Chem 288: 12437–12447. [PubMed]
48. Eaton AK, Stewart RC (2009) The two active sites of Thermotoga maritima CheA dimers bind ATP with dramatically different affinities. Biochemistry 48: 6412–6422. [PMC free article] [PubMed]
49. Airola MV, Sukomon N, Samanta D, Borbat PP, Freed JH, et al. (2013) HAMP domain conformers that propagate opposite signals in bacterial chemoreceptors. PLoS Biol 11: e1001479. [PMC free article] [PubMed]
50. Airola MV, Watts KJ, Bilwes AM, Crane BR (2010) Structure of concatenated HAMP domains provides a mechanism for signal transduction. Structure 18: 436–448. [PMC free article] [PubMed]
51. Crick FHC (1953) The packing of alpha-helices: simple coiled-coils. Acta Crystallogr 6: 689–697.
52. Dunin-Horkawicz S, Lupas AN (2010) Measuring the conformational space of square four-helical bundles with the program samCC. J Struct Biol 170: 226–235. [PubMed]
53. Hunke S, Betton JM (2003) Temperature effect on inclusion body formation and stress response in the periplasm of Escherichia coli. Mol Microbiol 50: 1579–1589. [PubMed]
54. Clarke EJ, Voigt CA (2011) Characterization of combinatorial patterns generated by multiple two-component sensors in E. coli that respond to many stimuli. Biotechnol Bioeng 108: 666–675. [PMC free article] [PubMed]
55. Parkinson JS (2010) Signaling mechanisms of HAMP domains in chemoreceptors and sensor kinases. Annu Rev Microbiol 64: 101–122. [PubMed]
56. Miot M, Betton J-M (2010) Reconstitution of the Cpx signaling system from cell-free synthesized proteins. N Biotechnol 283: 277–281. [PubMed]
57. Kabsch W (2010) Xds. Acta Crystallogr D Biol Crystallogr 66: 125–132. [PMC free article] [PubMed]
58. Collaborative Computational Project N (1994) The CCP4 suite: programs for protein crystallography. Acta Crystallogr D Biol Crystallogr 50: 760–763. [PubMed]
59. Terwilliger TC, Adams PD, Read RJ, McCoy AJ, Moriarty NW, et al. (2009) Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard. Acta Crystallogr D Biol Crystallogr 65: 582–601. [PMC free article] [PubMed]
60. Emsley P, Cowtan K (2004) Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr 60: 2126–2132. [PubMed]
61. McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, et al. (2007) Phaser crystallographic software. J Appl Crystallogr 40: 658–674. [PubMed]
62. Bricogne G, Blanc E, Brandl M, Flensburg C, Keller P, et al. . (2011) BUSTER version 2.11.1. Cambridge, UK: Global Phasing Ltd.
63. Davis I, Leaver-Fay A, Chen V, Block J, Kapral G, et al. (2007) MolProbity: all-atom contacts and structure validation for proteins and nucleic acids. Nucleic Acids Res 35: W375–383. [PMC free article] [PubMed]
64. DeLano W (2002) The PyMOL Molecular Graphics System. San Carlos, CA: DeLano Scientific.
65. Konarev PV, Petoukhov MV, Volkov VV, Svergun DI (2006) ATSAS 2.1, a program package for small-angle scattering data analysis. J Appl Crystallogr 39: 277–286.
66. Puertas JM, Betton JM (2009) Engineering an efficient secretion of leech carboxypeptidase inhibitor in Escherichia coli. Microb Cell Fact 8: 57. [PMC free article] [PubMed]
67. Liu HT, Naismith JH (2009) A simple and efficient expression and purification system using two newly constructed vectors. Prot Expr Purif 63: 102–111. [PMC free article] [PubMed]
Articles from PLoS Biology are provided here courtesy of
Public Library of Science