- Research article
- Open Access
Characterization of the χψ subcomplex of Pseudomonas aeruginosa DNA polymerase III
© El Houry Mignan et al; licensee BioMed Central Ltd. 2011
- Received: 24 May 2011
- Accepted: 28 September 2011
- Published: 28 September 2011
DNA polymerase III, the main enzyme responsible for bacterial DNA replication, is composed of three sub-assemblies: the polymerase core, the β-sliding clamp, and the clamp loader. During replication, single-stranded DNA-binding protein (SSB) coats and protects single-stranded DNA (ssDNA) and also interacts with the χψ heterodimer, a sub-complex of the clamp loader. Whereas the χ subunits of Escherichia coli and Pseudomonas aeruginosa are about 40% homologous, P. aeruginosa ψ is twice as large as its E. coli counterpart, and contains additional sequences. It was shown that P. aeruginosa χψ together with SSB increases the activity of its cognate clamp loader 25-fold at low salt. The E. coli clamp loader, however, is insensitive to the addition of its cognate χψ under similar conditions. In order to find out distinguishing properties within P. aeruginosa χψ which account for this higher stimulatory effect, we characterized P. aeruginosa χψ by a detailed structural and functional comparison with its E. coli counterpart.
Using small-angle X-ray scattering, analytical ultracentrifugation, and homology-based modeling, we found the N-terminus of P. aeruginosa ψ to be unstructured. Under high salt conditions, the affinity of the χψ complexes from both organisms to their cognate SSB was similar. Under low salt conditions, P. aeruginosa χψ, contrary to E. coli χψ, binds to ssDNA via the N-terminus of ψ. Whereas it is also able to bind to double-stranded DNA, the affinity is somewhat reduced.
The binding to DNA, otherwise never reported for any other ψ protein, enhances the affinity of P. aeruginosa χψ towards the SSB/ssDNA complex and very likely contributes to the higher stimulatory effect of P. aeruginosa χψ on the clamp loader. We also observed DNA-binding activity for P. putida χψ, making this activity most probably a characteristic of the ψ proteins from the Pseudomonadaceae.
- Electrophoretic Mobility Shift Assay
- Sedimentation Coefficient
- SAXS Data
- High Salt Condition
- High Salt Buffer
Pseudomonas aeruginosa is a ubiquitous Gram-negative bacterium that causes high rates of hospital-acquired infections, especially in immunocompromised patients  as well as in patients with cystic fibrosis . P. aeruginosa infections are often difficult to treat because the pathogen is capable of very rapidly acquiring a multitude of resistance mechanisms , making multidrug resistant P. aeruginosa strains increasingly common and raising the need for new antipseudomonal drugs. The replication machinery of P. aeruginosa presents a good target for the development of new antimicrobial agents, since DNA replication is essential for the survival of the pathogen.
DNA polymerase III holoenzyme is the main enzyme involved in the replication of the bacterial genome. In Escherichia coli, this multiprotein complex is composed of ten distinct subunits that are arranged into three functional subassemblies . The core subassembly contains the DNA polymerase and the proofreading exonuclease activities . The β-sliding clamp encircles the DNA duplex while binding to the core and thereby tethering it to the DNA template, ensuring high processivity and speed of the holoenzyme . The clamp loader complex, composed of subunits τ/γ, δ, δ', χ and ψ, uses ATP hydrolysis to load the clamp onto DNA .
The τ/γ, δ and δ' subunits together form the "minimal" clamp loader, (τ/γ)3δδ', which is sufficient for loading the clamp. The χ and ψ subunits are not essential for clamp loading . However, χ is the only direct link between DNA polymerase III and SSB, the single-stranded DNA-binding protein [9, 10] which coats single-stranded DNA (ssDNA) at the lagging strand to protect it from degradation and to prevent hairpin formation. The SSB/χ interaction is necessary for the primase-to-polymerase switch, which involves a competition between DnaG primase and χ in binding to SSB . The ψ subunit bridges χ to the (τ/γ)3δδ' complex and has been shown to play an important role in stabilizing the clamp loader  and increasing its affinity for the β-clamp . In E. coli, χ and ψ form a tight 1:1 complex, the crystal structure of which has been solved . The residues at the interface between χ and ψ are highly conserved across several bacterial species , and the same holds true for a hydrophobic surface pocket of χ surrounded by basic residues which was shown to interact with the amphipathic C-terminus of SSB . The disordered N-terminal part of ψ is also well-conserved, and is required for the binding of χψ to the minimal clamp loader . However, in several bacterial genomes containing χ sequences in which the ψ-binding site appears to be conserved, ψ sequences cannot be found, suggesting that ψ has greatly diverged in these species .
Sequence comparisons between E. coli and P. aeruginosa identified the existence of all genes encoding the subunits of DNA polymerase III except for θ, a non-essential subunit of the core, and ψ . The ψ subunit could only be identified as a 32 kDa protein which copurified with the other DNA polymerase III subunits from P. aeruginosa cell lysate . The start codon of the gene encoding this subunit was found to be incorrectly annotated in the P. aeruginosa PAO1 genome, the actual start codon being a non-canonical UUG present 135 nucleotides upstream from the originally annotated one. Examination of the ψ amino acid sequence showed that it is divergent from and nearly twice the length of E. coli ψ. At the functional level, the P. aeruginosa χψ complex appears to play a more significant role than its E. coli counterpart. Jarvis et al.  observed that, under low salt conditions and at subsaturating levels of the τ3δδ' minimal clamp loader, the addition of P. aeruginosa χψ and SSB has a synergistic effect, increasing the activity of P. aeruginosa τ3δδ' 25-fold. In contrast, under similar conditions, E. coli τ3δδ' was insensitive to the addition of E. coli χψ . At high salt concentrations, however, both χψ and SSB are required for efficient DNA synthesis, in the P. aeruginosa as well as in the E. coli systems .
This report describes the biophysical characterization of the χψ complex of P. aeruginosa, revealing properties which distinguish it from its E. coli counterpart and possibly account for its higher stimulatory effect. Using analytical ultracentrifugation, we show that P. aeruginosa ψ, but not the E. coli ψ subunit, binds to ssDNA under low salt conditions, even when the ssDNA is coated with SSB. This binding results in an increased affinity of P. aeruginosa χψ toward the SSB/ssDNA complex when compared to the E. coli system. In addition, we found that P. putida χψ also binds to ssDNA, suggesting that DNA-binding is a property of ψ proteins from the Pseudomonadaceae.
Biophysical characterization of Pae χψ
To determine the composition of P. aeruginosa χψ (Pae χψ), the complex was subjected to sedimentation equilibrium experiments in the analytical ultracentrifuge under high salt conditions. The concentration gradients obtained at three different rotor speeds (see Methods section) and two protein concentrations [5.8 μM (Additional file 1: Figure S1 A) and 23.3 μM] could be globally fitted with a single molar mass of 46 (± 4) kg/mol. There was no indication of multiple species or aggregation. Since the molar mass of a Pae χψ heterodimer calculated from the amino acid composition is 46.2 kg/mol, the protein exists as a heterodimer in solution. In comparison, examination of the χψ complex of E. coli (Eco χψ) yielded a molar mass of 31 (± 4) kg/mol (Additional file 1: Figure S1 B), confirming the formation of Eco χψ heterodimers (calculated molar mass 31.8 kg/mol) as reported previously .
To characterize its hydrodynamic properties, the Pae χψ heterodimer was examined in sedimentation velocity experiments in the analytical ultracentrifuge under high salt conditions. The data measured at three different protein concentrations were analyzed using the program package SEDFIT . For Pae χψ, the c(s) distributions in a concentration range of 1.9 to 15.5 μM revealed a single species, sedimenting with an s20, W of 2.5 S, and gave no indication of significant impurities or aggregation products (Additional file 2: Figure S2 A). Even at the lowest protein concentration used, no change in the c(s) distribution could be observed, showing that the χψ complex is stable under these conditions with no tendency to dissociate. In comparison, the c(s) analysis for Eco χψ resulted in s20, W= 2.5 S and the complex also showed no tendency to aggregate at higher or to dissociate at lower protein concentrations (Additional file 2: Figure S2 B).
From the calculated molar masses of the χψ heterodimers and the measured sedimentation coefficients, the frictional ratios f/f0 and the hydrodynamic radii rH were calculated. For the P. aeruginosa complex this resulted in f/f0= 1.80 and rH= 4.3 nm, whereas the E. coli complex yielded f/f0= 1.42 and rH= 3.0 nm. Since the frictional ratios for hydrated spherical proteins are expected to be in the range of 1.1 to 1.2 , the shape of both protein complexes seems to deviate substantially from a sphere. For Eco χψ, this is in agreement with the results from the crystal structure , where it was shown that the complex is elongated and that the 26 N-terminal amino acids of ψ are disordered. As the frictional ratio of Pae χψ is even larger, a GlobPlot analysis  of the χ and ψ proteins was performed, which predicted the presence of two large disordered regions in the N-terminus of Pae ψ, while Pae χ was predicted to be mostly globular (data not shown).
To check the GlobPlot prediction for Pae ψ, a truncated version of this protein was constructed in which the 85 N-terminal amino acids were deleted. This construct was coexpressed with full-length Pae χ and both proteins could be copurified as a Pae χψ(Δ1-85) complex. Although the molar mass of this complex was reduced by 9.1 kg/mol due to the truncation, it sedimented essentially with the same sedimentation coefficient as observed for full-length Pae χψ (s20, W= 2.5 S, see additional file 2: Figure S2 C). Accordingly, the frictional ratio and the hydrodynamic radius dropped to f/f0= 1.55 and rH= 3.4 nm, respectively. This result makes it very likely that the N-terminal part of Pae ψ is disordered. It is worth noting, however, that the frictional ratio of the truncated complex is still higher than that of Eco χψ.
SAXS measurements reveal that Pae χψ is elongated in solution
To get structural information of Pae χψ in solution, and to examine the lack of structure within the N-terminal part of Pae ψ in particular, we performed small-angle X-ray scattering (SAXS) experiments of the full-length and the truncated protein complexes (Additional file 3). High-throughput crystallization attempts with both Pae χψ and Pae χψ(Δ1-85) unfortunately failed to produce suitable crystals for X-ray crystallography, thus no high resolution structure could be obtained.
Based on HHpred analysis, the core domains of Pae χψ and Eco χψ show high homologies; we therefore tested whether the ab initio shape calculated from Pae χψ(Δ1-85) scattering data is able to harbor the in silico-generated Pae ψ-Eco χ heterodimer that we modeled. As Pae ψ is much larger than Eco ψ, this model accounts for the putative presence of extended loops in Pae ψ as predicted from sequence alignment (Figure 2A). The shape of the final ab initio SAXS model of Pae χψ(Δ1-85) (Figure 1C) is clearly able to harbor both the modeled χψ-complex and the N-terminal residues 86-125 of Pae ψ, which are not present in our model but in the measured construct.
We further analyzed the SAXS data of the Pae χψ full-length complex using an alternative approach: lacking data to calculate reliable ab initio models due to the large and probably flexible N-terminus of Pae ψ, we performed an ensemble optimization method approach (EOM) . This procedure creates a large pool of random configurations and a genetic algorithm is used to select ensembles being in agreement with the experimental scattering data. For this approach we used our in silico-generated construct of Pae ψ and Eco χ as a rigid (and folded) core and defined the absent N-terminus of Pae ψ as a flexible part containing 125 residues. The RG-distribution obtained for the models in the selected ensembles is quite broad (Figure 1B), supporting the idea of high flexibility (e.g. ). The chosen ensemble of structures was then superposed (Figure 1D) to visualize the different conformations of the N-terminus which all fit to the measured scattering curve. Additionally, we calculated the theoretical sedimentation coefficients of all structures in the chosen ensemble using HYDROPRO . The average s-value of s20, W = 2.6 ± 0.1 S is in very good agreement with our experimental data.
Pae χψ binds to the highly conserved C-terminus of SSB
For Pae χψ and Pae SSB, a similar binding affinity of 1.3 × 105 M-1 was obtained under the same experimental conditions (Figure 3). The binding stoichiometry, however, is lower than in the case of E. coli. As both SSB proteins are known to form homotetramers [31, 32], there exist four potential binding sites for χψ, one at the C-terminus of each monomer. Since the binding isotherm yielded a stoichiometry of n= 3.8 for the E. coli system, each C-terminus of Eco SSB binds one Eco χψ complex (Figure 3). Correspondingly, it has been shown that up to four Eco χ proteins can bind to Eco SSB . In the case of P. aeruginosa, however, not all binding sites can be occupied (n= 2.2, Figure 3). This might be due to the fact that Pae ψ is almost twice the size of Eco ψ and that, as was shown above, the N-terminus of Pae ψ is highly elongated. Thus, it is most probably sterical hindrance that prevents more than two Pae χψ complexes from binding to SSB via χ.
To test whether Pae χψ has the same binding site on SSB as the E. coli complex, we used the deletion mutant Eco SSBQ152* , which is truncated of the last 26 amino acids and which has been shown to be unable to interact with Eco χ . Under high salt conditions, we detected no significant interaction of either Pae or Eco χψ with EcoSSBQ152*, whereas Pae χψ interacted with wild-type Eco SSB with a similar affinity than with Pae SSB (data not shown). Since the 120 N-terminal amino acids of Pae SSB and Eco SSB show 79% homology and the last 7 amino acids are identical , this finding makes it very likely that the highly conserved C-terminus of Pae SSB is the only binding site for Pae χψ.
Pae ψ shows DNA-binding activity
Since the SSB proteins involved in DNA replication form a tight complex with ssDNA, we examined whether Pae χψ not only binds to SSB but also directly to nucleic acids. We therefore tested the binding of Pae χψ to fluorescently labeled 50mer-dsDNA and 50mer-ssDNA in electrophoretic mobility shift assays (EMSA) under low salt conditions. Both substrates clearly show shifts with increasing concentrations of Pae χψ, indicating binding of the protein complex to the nucleic acids (Additional file 5: Figure S5 A and B). However, ssDNA is bound with a somewhat higher affinity than dsDNA (Additional file 5: Figure S5 C), suggesting that the complex seems to preferentially bind to ssDNA.
To characterize this binding, the amount of free Pae χψ was determined from the c(s) distributions, and a binding isotherm was constructed (Figure 4B). The binding stoichiometry, defined as the amount of χψ bound per dT residue, was found to be 0.12. Therefore, one Pae χψ complex is able to bind about eight dT residues with a binding affinity of approximately 5 × 105 M-1. It is worth noting, however, that the binding affinity of Pae χψ to ssDNA is salt-dependant, since it was weakened at increasing salt concentrations, and no binding was observed in high salt buffer (data not shown). In this context, it is interesting to note that under low salt conditions, but not at elevated salt, ssDNA increases the affinity of Eco χ to Eco SSB about 20-fold . Under elevated salt conditions, however, in addition to ssDNA, the presence of all other clamp loader subunits, which themselves do not significantly interact with the SSB/ssDNA complex, is required to increase the SSB-binding affinity of Eco χψ about 1000-fold . This might be mediated by a conformational change in the N-terminus of ψ which accompanies the assembly of χψ into the clamp loader complex . A similar conformational change is expected to occur in the conserved ultimate N-terminus of Pae ψ, which might also increase the ssDNA-binding affinity of Pae χψ under elevated salt conditions.
To determine which protein within the Pae χψ heterodimer is responsible for ssDNA-binding, we wanted to investigate the interaction of poly(dT) with only Pae χ or only Pae ψ. However, both proteins were insoluble if expressed on their own in E. coli and could not be purified. Therefore, we tried to produce the chimeric complexes χEco ψ Pae and χPae ψEco, but only Eco χ and Pae ψ formed a soluble complex. Although these two proteins are from different organisms, they formed a stable homogeneous heterodimeric χEco ψ Pae complex as revealed by sedimentation velocity experiments (data not shown). This is most probably due to the fact that the residues at the χψ interface are conserved between E. coli and P. aeruginosa. Similarly to Pae χψ, χEcoψ Pae is able to bind to poly(dT) under low salt conditions (Additional file 6). Since we showed that Eco χψ does not bind ssDNA (Figure 4A), and under similar conditions no direct interaction between Eco χ and ssDNA could be detected [9, 10, 30, 34], it has to be the ψ subunit within Pae χψ which binds to ssDNA.
To check whether the apparently unstructured N-terminal region of Pae χψ is involved in the interaction with ssDNA, we tested the binding of the truncated Pae χψ(Δ1-85) complex to poly(dT). As can be seen by the smaller increase in the sedimentation coefficient of the protein/ssDNA complex compared to full-length Pae χψ (Figure 4A), less Pae χψ(Δ1-85) can bind to ssDNA. Therefore, the N-terminal region of ψ is involved in ssDNA-binding.
P. putida χψ also binds to ssDNA
When the ψ subunit of P. aeruginosa was discovered , a whole family of ψ proteins present within several species of the Pseudomonadaceae was identified based on sequence comparisons. To test whether ssDNA-binding is a general property of the pseudomonal ψ proteins, we cloned and purified the χψ complex from P. putida. In sedimentation velocity experiments with poly(dT) in low salt buffer, P. putida (Ppu) χψ, similarly to Pae χψ, increased the sedimentation coefficient of poly(dT), indicating ssDNA-binding (Additional file 6).
Pae χψ binds to ssDNA covered by Eco SSB+Gly
The ssDNA at the replication fork is entirely covered by SSB [31, 35]. To have an in vitro setting which most resembles the in vivo situation, the binding of Pae χψ to ssDNA complexed with SSB was analyzed. As this experiment is complicated by the fact that χψ also binds to SSB, we used a C-terminal extension mutant of Eco SSB, Eco SSB+Gly. This protein carries an additional glycine residue at its C-terminus which dramatically weakens its interaction with Eco χ . Using sedimentation velocity experiments, we first confirmed that the binding of both Pae and Eco χψ to this mutant is drastically reduced (Additional file 7).
Pae χψ binds with a higher affinity to SSB/ssDNA complexes than Eco χψ
In the absence of poly(dT), the binding was too weak to determine the binding parameters for either P. aeruginosa or E. coli (Figure 6). Due to the aggregation of Pae χψ at higher protein concentrations at low salt, the binding could not be examined at concentrations which would have allowed for the determination of the binding parameters. The binding of Pae χψ to the SSB/ssDNA complex, however, is so strong that even at the lowest concentration that can be applied in AUC experiments using the absorbance optics, stoichiometric binding was observed and the affinity is larger than 5 × 107 M-1 (Figure 6A). Therefore, the affinity of Pae χψ to SSB/poly(dT) is at least eight times higher than the affinity of the respective E. coli proteins (K= 6 × 106 M-1, Figure 6B). This is most probably a result of the ssDNA-binding affinity of Pae ψ which enhances the stability of the ternary complex, when compared to that of the E. coli system.
The χ and ψ subunits are part of the clamp loader of DNA polymerase III holoenzyme. Whereas Eco χ and Pae χ have a similar size and show about 40% sequence homology, Pae ψ is nearly twice as large as its E. coli counterpart. The N-terminal residues of ψ that are needed for the interaction with the minimal clamp loader  and the residues involved in binding χ  are conserved, however Pae ψ contains additional sequences that are absent from the E. coli protein. We showed by homology-based modeling that these additional sequences do not interfere with secondary structure elements of Eco ψ, but can be assigned to loop regions of the protein (Figure 2B).
Analytical ultracentrifugation experiments revealed that Pae χψ has a higher degree of asymmetry than its E. coli counterpart. This asymmetry is reduced by the truncation of the 85 N-terminal amino acids of Pae ψ. SAXS experiments of Pae χψ and of this truncated variant, as well as molecular modeling, strongly support the idea that the N-terminus of Pae ψ is elongated and probably highly flexible. Our results show that this N-terminus is responsible for ssDNA-binding of Pae χψ under low salt conditions, a property which can not be found in Eco χψ.
Additionally, we could show that also P. putida χψ is able to bind to ssDNA. To our knowledge, this is the first report on the DNA-binding of ψ subunits of DNA polymerase III. Homology searches within the Pseudomonadaceae revealed several proteins related to Pae ψ, all containing extra sequences absent from Eco ψ, which were in many cases highly conserved . One conserved region is located in the unstructured N-terminus of ψ, includes several positively charged amino acids and might be responsible for the binding of the pseudomonal ψ proteins to DNA (Additional file 8).
Using an Eco SSB mutant impaired in χψ binding, we showed that Pae χψ is even able to bind to ssDNA that is covered with SSB. Furthermore, we demonstrated that the affinity of Pae χψ to SSB/ssDNA complexes is significantly increased compared to the E. coli system. This fact is most probably due to the ability of Pae ψ to bind to ssDNA. It has been shown that under low salt conditions only Pae χψ has a significant influence on DNA synthesis by its cognate DNA polymerase III [17, 18]. The addition of Pae χψ and SSB shows a synergistic effect, increasing the activity of the minimal clamp loader Pae τ3δδ' 25-fold, whereas Eco τ3δδ' is insensitive to the addition of Eco χψ [17, 18]. Even though it might not be the only mechanism, the binding of Pae ψ to DNA very likely contributes to the higher stimulatory effect of Pae χψ on the minimal clamp loader in the presence of SSB. Since we showed that the χψ complex of P. putida is also able to bind DNA, this stimulatory effect might be a characteristic of pseudomonal χψ complexes.
In spite of the fundamental aspects of this work, the full understanding of the DNA polymerase III holoenzyme of P. aeruginosa, and of its ψ subunit in particular, may reveal new targets for the development of specific inhibitors of pseudomonal DNA replication. The most striking difference between the E. coli and pseudomonal χψ complexes identified in this work is the ssDNA-binding ability mediated by the N-terminus of Pae ψ. Since we showed that this region is unstructured, it can probably not be used directly in antipseudomonal drug design. However, it is tempting to speculate that this region might fold upon ssDNA-binding and one of our future objectives is to solve the structure of a Pae χψ/ssDNA complex.
Buffers and reagents
Poly(dT) (~1400 nt in length) was purchased from GE Healthcare Life Sciences. Its concentration is given in monomer residues throughout the text and was determined spectrophotometrically, using an absorption coefficient of 8600 M-1 cm-1 at maximum . Protein concentrations were determined spectrophotometrically, using the absorption coefficients at 280 nm calculated from amino acid composition : Pae SSB: 95800 M-1 cm-1; Pae χψ: 51450 M-1 cm-1; Pae χψ(Δ1-85): 45950 M-1 cm-1; Ppu χψ: 55920 M-1 cm-1; Eco χψ: 52940 M-1 cm-1; chimeric χEco ψPae: 63940 M-1 cm-1. For wild-type Eco SSB and Eco SSB+Gly an absorption coefficient of 113000 M-1 cm-1 was used .
Experiments were carried out in potassium phosphate buffer (KPi) pH 7.4, containing two different NaCl concentrations: high salt buffer (0.3 M NaCl, 20 mM KPi pH 7.4, 0.5 mM DTT), and low salt buffer (5 mM NaCl, 5 mM KPi pH 7.4, 0.87 M glycerol, 0.5 mM DTT).
Cloning of the χ and ψ genes of P. aeruginosa, E. coli, and P. putida
Primer sequences are given in additional file 9. The holC and holD genes of P. aeruginosa were amplified by PCR from genomic DNA of strain PAO1. For the holD gene, the start codon was at the position identified in . After Bsa I/Bam HI digestion, the holC PCR product was cloned into pCDFDuet-1 (Novagen), cut with Nco I/Bam HI. The holD PCR product was cut with Bgl II/Eco RI and cloned into the pGEX-6P-1 vector (GE Healthcare Life Sciences), digested with Bam HI/Eco RI.
The holC and holD genes of E. coli were amplified by PCR from genomic DNA of E. coli strain LK111λ . The holC PCR product was cut with Nco I/Xho I, and cloned into pET-15b (Novagen), cut with the same enzymes. This construct was then digested with Xho I, and the holD PCR product, cut with Xho I, was ligated in.
For expression of the chimeric χ Eco ψ Pae complex, the holC gene was amplified by PCR from genomic DNA of E. coli LK111λ. After digestion with Nde I/Xho I, the PCR product was cloned into the vector pCDFDuet-1.
The holC and holD genes of P. putida were amplified by PCR from genomic DNA of P. putida strain KT2440  and were cloned into the vectors pRSFDuet-1 (Novagen) and pETGEX-6P, respectively. To construct the latter, the sequence coding for GST and the multiple cloning site were amplified from pGEX-6P-1, and were cloned into pET-15b. In the resulting vector, protein expression is under control of the T7lac promoter and is therefore more tightly regulated than in pGEX-6P-1. For the holD gene, the start codon was at the position identified in . After digestion with Bsa I/Bam HI, the holC PCR product was cloned into pRSFDuet-1, cut with Nco I/Bam HI. The holD PCR product was cut with Bam HI/Xho I and ligated into pETGEX-6P.
All constructs were checked for errors by sequencing the complete genes (GATC Biotech).
Construction of the N-terminally truncated mutant Pae ψ(Δ1-85)
A 582 bp fragment, encoding Pae ψ truncated of its first 85 amino acids (Pae ψ(Δ1-85)), was amplified by PCR from the genomic DNA of P. aeruginosa strain PAO1 (for sequence of primers see additional file 9). After digestion with Bgl II/Eco RI, the PCR fragment was cloned into pGEX-6P-1, cut with Bam HI/Eco RI, allowing for the expression of N-terminally GST-tagged Pae ψ(Δ1-85).
Protein expression and purification
Pae χψ, Pae χψ(Δ1-85), chimeric χ Eco ψ Pae and Ppu χψ, were expressed as χ/GSTψ. During the purification procedure the GST moiety was cleaved off by PreScission™ protease digestion. Eco χψ was expressed untagged. Expression of the χ proteins was done using Duet vectors (Novagen). While Pae ψ and Pae ψ(Δ1-85) were expressed using pGEX-6P-1, Ppu ψ was expressed from pETGEX-6P.
E. coli Rosetta (DE3) pLyS cells (Novagen) were transformed with the respective vectors encoding χ and GSTψ proteins and were grown to E600 nm= 1.2 in a 10 L fermentor flask at 37°C in LB medium. 5 h after addition of 1 mM IPTG cells were harvested, resuspended in an equal volume of PBS (300 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4, 10% (v/v) glycerol, 1 mM DTT, pH 7.3) and frozen in N2 (liq). 40 g of frozen cells were thawed in the presence of two volumes of PBS buffer, 1 mM EDTA, 1.67 mM DTT, 0.1 mM PMSF, and 10 tablets of complete EDTA-free protease inhibitor (Roche). The suspension was incubated for 20 min at 4°C with 0.3 mg/mL lysozyme, followed by sonification. After centrifugation, the supernatant was applied onto a glutathione sepharose 4B (GE Healthcare Life Sciences) column preequilibrated with PBS buffer. After washing with GS buffer (50 mM Tris-HCl pH 7.0, 150 mM NaCl, 10% glycerol (v/v), 1 mM EDTA, 1 mM DTT), the column was incubated overnight at 4°C with 840 μg of PreScission™ protease (GE Healthcare Life Sciences). The χ and ψ proteins were eluted using GS300 buffer (GS buffer with 0.3 M NaCl). Fractions containing the purest complex were pooled and the proteins were precipitated for 1 h with 164 g/L (NH4)2SO4 (except for Pae χψ(Δ1-85) where 300 g/L were used). After centrifugation, the precipitated proteins were resuspended in GS300 buffer and subjected to size exclusion chromatography (Superdex-75-prepgrade, GE Healthcare Life Sciences) in the same buffer. Fractions that contained contamination-free χψ complex were pooled, dialyzed overnight against high salt buffer containing 10% (v/v) glycerol, flash frozen in N2 (liq) and stored at -80°C.
Expression of the Eco χψ complex was done using the pET-15b vector (Novagen). For purification, cells harvested from the fermentor culture were resuspended in 50 mM Tris-HCl pH 7.5, 10% (w/v) sucrose and frozen in N2 (liq). 30 g of frozen cells were thawed in the presence of 0.1 mM PMSF and cell lysis was done as above, but using 0.1 mg/mL lysozyme. After centrifugation, the proteins in the supernatant were precipitated with 250 g/L (NH4)2SO4 at 4°C. The protein pellet after centrifugation was resuspended in buffer A20 (50 mM Tris pH 7.8, 20 mM NaCl, 10% (v/v) glycerol, 1 mM EDTA, 1 mM DTT), and dialyzed against the same buffer. The solution was applied onto a Q-sepharose fast flow (GE Healthcare Life Sciences) column preequilibrated in the same buffer. After washing with buffer A20, the proteins were eluted in a gradient of 20 mM to 250 mM NaCl in buffer A20. The fractions containing Eco χ and ψ were pooled and the protein complex was precipitated with 300 g/L of (NH4)2SO4 at 4°C. Following centrifugation, the protein pellet was resuspended in buffer A50 (A20 with 50 mM NaCl) and applied onto a size exclusion chromatography column (see above) preequilibrated with the same buffer. Fractions containing a contamination-free χψ complex were pooled and dialyzed overnight against 20 mM Hepes pH 8.0, 100 mM NaCl, 2 mM DTT, 1 mM NaN3, and 3% (w/v) sucrose, flash frozen in N2 (liq), and stored at -80°C.
Analytical ultracentrifugation experiments were performed either in a Beckman Optima XL-A ultracentrifuge or a Beckman/Coulter ProteomeLab XL-I ultracentrifuge, using An-50 Ti rotors. Concentration profiles were measured with the UV-absorption scanning optics of the centrifuge.
Sedimentation velocity experiments were carried out at 20°C in either 12 mm or 3 mm standard double-sector centerpieces filled with 400 μL or 100 μL sample, respectively, at the indicated rotor speeds. For the analysis of protein-protein or protein-DNA interactions, the slower sedimenting molecule was titrated to a constant concentration of the faster sedimenting one. Since the reactions governing the interactions examined in this study are fast compared to the timescale of sedimentation, only two sedimenting boundaries are observed . The slow one represents the free slower sedimenting molecule, whereas the fast one contains complexes of both components and free faster sedimenting molecules.
The measured concentration profiles were evaluated using SEDFIT  which transforms them into diffusion-corrected sedimentation coefficient distributions [c(s) distributions]. As the areas under the separate peaks in the c(s) distributions are a measure of the absorbance of the species represented by the peaks , this information can be used to determine binding isotherms , which were evaluated using a simple interaction model as described previously .
In the case of the χψ/SSB/poly(dT) experiments where a wavelength of 230 nm was used, the extinction coefficients of χψ and the SSB/poly(dT) complexes were determined from the observed absorbances in the centrifuge and the known initial concentrations of the components. Thus, it was possible to correct for variations in the exact wavelength provided by the monochromator of the centrifuge.
For hydrodynamic analysis, measured s-values were corrected to s20, W, using the partial specific volumes calculated from amino acid composition . Since the partial specific volume of complexes of different macromolecules with unknown composition can not be calculated, uncorrected sedimentation coefficients are given in these cases.
Sedimentation equilibrium experiments were carried out at 4°C in either 12 mm or 3 mm double-sector centerpieces filled with 150 μL or 40 μL sample, respectively. Absorbance was detected at 280 nm and rotor speeds of 9000, 13000 and 18000 rpm were used. Analysis of the data was performed as described earlier .
SAXS measurements and modeling of Pae ψ
Sample preparation for SAXS included dialysis in 20 mM Hepes pH 7.4, 50 mM NaCl, 3% (w/v) sucrose, 1 mM NaN3, 1 mM DTT, followed by protein concentration using Vivaspin centrifugal concentrators (Sartorius Stedim Biotech). The flowthrough was taken as buffer reference in SAXS data collection. SAXS data were collected at EMBL/DESY X33 beamline (Hamburg, Germany) at 20°C cell temperature. Molar mass of the samples was determined by comparison of scattering intensity at zero angle I(0) obtained from Guinier analysis (s*RG < 1.3 in ln(I) vs. s2-plot) with reference proteins BSA and lysozyme, and was additionally determined from volume (Porod volume, see e.g. ).
SAXS data were processed and analyzed using the ATSAS package . Pae χψ and Pae χψ(Δ1-85) were measured at protein concentrations of 2, 4 and 8 mg/ml and 1.3, 3.2 and 6.4 mg/ml, respectively. All sample scattering curves were corrected by subtraction of the corresponding buffer scattering curves. Scattering curves were merged as described in e.g.  to obtain optimal data quality at lower and higher s-ranges. For full-length Pae χψ, a flexible modeling approach was performed, using the ensemble optimization method (EOM) described in e.g. [26, 27]. For the truncated Pae χψ(Δ1-85) construct, an independent set of ten ab initio bead models was calculated with GASBOR without prior symmetry information or other restrictions. Models were aligned and averaged using DAMAVER. For better representation, an electron density envelope of the averaged bead model was calculated using the SITUS package . SAXS-shape representation and docking into the SAXS shape were performed with the UCSF chimera package  and the superposition of structures was done using PyMOL .
Modeling of the Pae ψ molecule was performed using the Bioinformatic Toolkit server . The Pae ψ sequence was first analyzed by HHpred , identifying Eco ψ (pdb code: 1em8) as one of the best candidates with highest score (probability 96.66, E-value= 0.02, Identities 21%). This structure was then chosen as a template for MODELLER .
Electrophoretic mobility shift assays (EMSA)
EMSA was performed in 5% polyacrylamide gels in Tris-borate buffer (running buffer: 45 mM Tris, 45 mM boric acid; gels: 25 mM Tris, 25 mM boric acid) using fluorescently FAM- labeled oligonucleotides (Thermo Fischer). Gels were analyzed using a Typhoon Scanner (488 nm laser, GE Healthcare Life Sciences), and quantified using the ImageQuant software (GE Healthcare Life Sciences). 1 μl of DNA was added to the protein in 9 μl low salt buffer to yield the indicated final concentrations. The mixture was incubated for 10 min at room temperature prior to gel electrophoresis at 5°C. Double-stranded DNA was prepared by annealing of oligonucleotides A and B; oligonucleotide A was used as the ssDNA sample (A: 5'6FAM-GGATACGTAACAACGCTTATGCATCGCCGCCGCTACATCCCTGAGCTGAC 3'. B: 5'GTCAGCTCAGGGATGTAGCGGCGGCGATGCATAAGCGTTGTTACGTATCC 3').
The authors declare that they have no competing interests.
The authors thank Lidia Litz and Olga Fettscher for excellent technical assistance, Dr. Claus Urbanke for help with the analytical ultracentrifuge, valuable discussions and for critically reading the manuscript, Dr. Lutz Wiehlmann for supplying the PAO1 genomic DNA and the KT2440 strain and for helpful discussions, Drs. Joachim Greipel, Karl-Peter Hopfner and Dietmar Manstein for helpful discussions, and Mathias Flüggen for performing some of the AUC runs. We thank the staff of X33 beamline at the EMBL/DESY Hamburg for help, excellent support during data collection and fruitful discussions. This work was supported by a Georg Christoph Lichtenberg Stipend provided by the Ph.D. program Infection Biology sponsored by the Ministry of Science and Culture of Lower Saxony to SEHM. Article-processing charges were funded by the program 'Open Access Publication' of the Deutsche Forschungsgemeinschaft (DFG).
- Sadikot RT, Blackwell TS, Christman JW, Prince AS: Pathogen-host interactions in Pseudomonas aeruginosa pneumonia. Am J Respir Crit Care Med. 2005, 171 (11): 1209-1223. 10.1164/rccm.200408-1044SOPubMed CentralView ArticlePubMedGoogle Scholar
- Buchanan PJ, Ernst RK, Elborn JS, Schock B: Role of CFTR, Pseudomonas aeruginosa and Toll-like receptors in cystic fibrosis lung inflammation. Biochem Soc Trans. 2009, 37 (Pt 4): 863-867.View ArticlePubMedGoogle Scholar
- Strateva T, Yordanov D: Pseudomonas aeruginosa - a phenomenon of bacterial resistance. J Med Microbiol. 2009, 58 (Pt 9): 1133-1148.View ArticlePubMedGoogle Scholar
- Onrust R, Finkelstein J, Naktinis V, Turner J, Fang L, O'Donnell M: Assembly of a chromosomal replication machine: two DNA polymerases, a clamp loader, and sliding clamps in one holoenzyme particle. I. Organization of the clamp loader. J Biol Chem. 1995, 270 (22): 13348-13357. 10.1074/jbc.270.22.13348View ArticlePubMedGoogle Scholar
- Scheuermann R, Tam S, Burgers PM, Lu C, Echols H: Identification of the epsilon-subunit of Escherichia coli DNA polymerase III holoenzyme as the dnaQ gene product: a fidelity subunit for DNA replication. Proc Natl Acad Sci USA. 1983, 80 (23): 7085-7089. 10.1073/pnas.80.23.7085PubMed CentralView ArticlePubMedGoogle Scholar
- Stukenberg PT, Studwell-Vaughan PS, O'Donnell M: Mechanism of the sliding beta-clamp of DNA polymerase III holoenzyme. J Biol Chem. 1991, 266 (17): 11328-11334.PubMedGoogle Scholar
- Pritchard AE, Dallmann HG, Glover BP, McHenry CS: A novel assembly mechanism for the DNA polymerase III holoenzyme DnaX complex: association of deltadelta' with DnaX(4) forms DnaX(3)deltadelta'. EMBO J. 2000, 19 (23): 6536-6545. 10.1093/emboj/19.23.6536PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao H, Dong Z, O'Donnell M: DNA polymerase III accessory proteins. IV. Characterization of chi and psi. J Biol Chem. 1993, 268 (16): 11779-11784.PubMedGoogle Scholar
- Kelman Z, Yuzhakov A, Andjelkovic J, O'Donnell M: Devoted to the lagging strand-the subunit of DNA polymerase III holoenzyme contacts SSB to promote processive elongation and sliding clamp assembly. EMBO J. 1998, 17 (8): 2436-2449. 10.1093/emboj/17.8.2436PubMed CentralView ArticlePubMedGoogle Scholar
- Witte G, Urbanke C, Curth U: DNA polymerase III chi subunit ties single-stranded DNA binding protein to the bacterial replication machinery. Nucleic Acids Res. 2003, 31 (15): 4434-4440. 10.1093/nar/gkg498PubMed CentralView ArticlePubMedGoogle Scholar
- Yuzhakov A, Kelman Z, O'Donnell M: Trading places on DNA--a three-point switch underlies primer handoff from primase to the replicative DNA polymerase. Cell. 1999, 96 (1): 153-163. 10.1016/S0092-8674(00)80968-XView ArticlePubMedGoogle Scholar
- Olson MW, Dallmann HG, McHenry CS: DnaX complex of Escherichia coli DNA polymerase III holoenzyme. The chi psi complex functions by increasing the affinity of tau and gamma for delta.delta' to a physiologically relevant range. J Biol Chem. 1995, 270 (49): 29570-29577. 10.1074/jbc.270.49.29570View ArticlePubMedGoogle Scholar
- Anderson SG, Williams CR, O'Donnell ME, Bloom LB: A function for the psi subunit in loading the Escherichia coli DNA polymerase sliding clamp. J Biol Chem. 2007, 282 (10): 7035-7045.View ArticlePubMedGoogle Scholar
- Gulbis JM, Kazmirski SL, Finkelstein J, Kelman Z, O'Donnell M, Kuriyan J: Crystal structure of the chi:psi sub-assembly of the Escherichia coli DNA polymerase clamp-loader complex. Eur J Biochem. 2004, 271 (2): 439-449. 10.1046/j.1432-1033.2003.03944.xView ArticlePubMedGoogle Scholar
- Naue N, Fedorov R, Pich A, Manstein DJ, Curth U: Site-directed mutagenesis of the chi subunit of DNA polymerase III and single-stranded DNA-binding protein of E. coli reveals key residues for their interaction. Nucleic Acids Res. 2011, 39 (4): 1398-1407. 10.1093/nar/gkq988PubMed CentralView ArticlePubMedGoogle Scholar
- Simonetta KR, Kazmirski SL, Goedken ER, Cantor AJ, Kelch BA, McNally R, Seyedin SN, Makino DL, O'Donnell M, Kuriyan J: The mechanism of ATP-dependent primer-template recognition by a clamp loader complex. Cell. 2009, 137 (4): 659-671. 10.1016/j.cell.2009.03.044PubMed CentralView ArticlePubMedGoogle Scholar
- Jarvis TC, Beaudry AA, Bullard JM, Janjic N, McHenry CS: Reconstitution of a minimal DNA replicase from Pseudomonas aeruginosa and stimulation by non-cognate auxiliary factors. J Biol Chem. 2005, 280 (9): 7890-7900.View ArticlePubMedGoogle Scholar
- Jarvis TC, Beaudry AA, Bullard JM, Ochsner U, Dallmann HG, McHenry CS: Discovery and characterization of the cryptic psi subunit of the pseudomonad DNA replicase. J Biol Chem. 2005, 280 (49): 40465-40473. 10.1074/jbc.M508310200View ArticlePubMedGoogle Scholar
- Schuck P: Size-distribution analysis of macromolecules by sedimentation velocity ultracentrifugation and lamm equation modeling. Biophys J. 2000, 78 (3): 1606-1619. 10.1016/S0006-3495(00)76713-0PubMed CentralView ArticlePubMedGoogle Scholar
- Lebowitz J, Lewis MS, Schuck P: Modern analytical ultracentrifugation in protein science: a tutorial review. Protein Sci. 2002, 11 (9): 2067-2079.PubMed CentralView ArticlePubMedGoogle Scholar
- Linding R, Russell RB, Neduva V, Gibson TJ: GlobPlot: Exploring protein sequences for globularity and disorder. Nucleic Acids Res. 2003, 31 (13): 3701-3708. 10.1093/nar/gkg519PubMed CentralView ArticlePubMedGoogle Scholar
- Volkov VV, Svergun DI: Uniqueness of ab initio shape determination in small-angle scattering. J Appl Cryst. 2003, 36: 860-864. 10.1107/S0021889803000268.View ArticleGoogle Scholar
- Putnam CD, Hammel M, Hura GL, Tainer JA: X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution. Q Rev Biophys. 2007, 40 (3): 191-285.View ArticlePubMedGoogle Scholar
- Biegert A, Mayer C, Remmert M, Soding J, Lupas AN: The MPI Bioinformatics Toolkit for protein sequence analysis. Nucleic Acids Res. 2006, W335-9. 34 Web ServerGoogle Scholar
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234 (3): 779-815. 10.1006/jmbi.1993.1626View ArticlePubMedGoogle Scholar
- Bernado P, Mylonas E, Petoukhov MV, Blackledge M, Svergun DI: Structural characterization of flexible proteins using small-angle X-ray scattering. J Am Chem Soc. 2007, 129 (17): 5656-5664. 10.1021/ja069124nView ArticlePubMedGoogle Scholar
- Mertens HD, Svergun DI: Structural characterization of proteins and complexes using small-angle X-ray solution scattering. J Struct Biol. 2010, 172 (1): 128-141. 10.1016/j.jsb.2010.06.012View ArticlePubMedGoogle Scholar
- Garcia De La Torre J, Huertas ML, Carrasco B: Calculation of hydrodynamic properties of globular proteins from their atomic-level structure. Biophys J. 2000, 78 (2): 719-730. 10.1016/S0006-3495(00)76630-6PubMed CentralView ArticlePubMedGoogle Scholar
- Witte G, Fedorov R, Curth U: Biophysical analysis of Thermus aquaticus single-stranded DNA binding protein. Biophys J. 2008, 94 (6): 2269-2279. 10.1529/biophysj.107.121533PubMed CentralView ArticlePubMedGoogle Scholar
- Glover BP, McHenry CS: The chi psi subunits of DNA polymerase III holoenzyme bind to single-stranded DNA-binding protein (SSB) and facilitate replication of an SSB-coated template. J Biol Chem. 1998, 273 (36): 23476-23484. 10.1074/jbc.273.36.23476View ArticlePubMedGoogle Scholar
- Weiner JH, Bertsch LL, Kornberg A: The deoxyribonucleic acid unwinding protein of Escherichia coli. Properties and functions in replication. J Biol Chem. 1975, 250 (6): 1972-1980.PubMedGoogle Scholar
- Genschel J, Litz L, Thole H, Roemling U, Urbanke C: Isolation, sequencing and overproduction of the single-stranded DNA binding protein from Pseudomonas aeruginosa PAO. Gene. 1996, 182 (1-2): 137-143. 10.1016/S0378-1119(96)00535-5View ArticlePubMedGoogle Scholar
- Curth U, Genschel J, Urbanke C, Greipel J: In vitro and in vivo function of the C-terminus of Escherichia coli single-stranded DNA binding protein. Nucleic Acids Res. 1996, 24 (14): 2706-2711. 10.1093/nar/24.14.2706PubMed CentralView ArticlePubMedGoogle Scholar
- Kozlov AG, Jezewska MJ, Bujalowski W, Lohman TM: Binding specificity of Escherichia coli single-stranded DNA binding protein for the chi subunit of DNA pol III holoenzyme and PriA helicase. Biochemistry. 2010, 49 (17): 3555-3566. 10.1021/bi100069sPubMed CentralView ArticlePubMedGoogle Scholar
- Molineux IJ, Pauli A, Gefter ML: Physical studies of the interaction between the Escherichia coli DNA binding protein and nucleic acids. Nucleic Acids Res. 1975, 2 (10): 1821-1837. 10.1093/nar/2.10.1821PubMed CentralView ArticlePubMedGoogle Scholar
- Williams KR, Spicer EK, LoPresti MB, Guggenheimer RA, Chase JW: Limited proteolysis studies on the Escherichia coli single-stranded DNA binding protein. Evidence for a functionally homologous domain in both the Escherichia coli and T4 DNA binding proteins. J Biol Chem. 1983, 258 (5): 3346-3355.PubMedGoogle Scholar
- Urbanke C, Schaper A: Kinetics of binding of single-stranded DNA binding protein from Escherichia coli to single-stranded nuclei acids. Biochemistry. 1990, 29 (7): 1744-1749. 10.1021/bi00459a012View ArticlePubMedGoogle Scholar
- Pace CN, Vajdos F, Fee L, Grimsley G, Gray T: How to measure and predict the molar absorption coefficient of a protein. Protein Sci. 1995, 4 (11): 2411-2423. 10.1002/pro.5560041120PubMed CentralView ArticlePubMedGoogle Scholar
- Lohman TM, Overman LB: Two binding modes in Escherichia coli single strand binding protein-single stranded DNA complexes. Modulation by NaCl concentration. J Biol Chem. 1985, 260 (6): 3594-3603.PubMedGoogle Scholar
- Zabeau M, Stanley KK: Enhanced expression of cro-beta-galactosidase fusion proteins under the control of the PR promoter of bacteriophage lambda. EMBO J. 1982, 1 (10): 1217-1224.PubMed CentralPubMedGoogle Scholar
- Regenhardt D, Heuer H, Heim S, Fernandez DU, Strompl C, Moore ER, Timmis KN: Pedigree and taxonomic credentials of Pseudomonas putida strain KT2440. Environ Microbiol. 2002, 4 (12): 912-915. 10.1046/j.1462-2920.2002.00368.xView ArticlePubMedGoogle Scholar
- Dam J, Velikovsky CA, Mariuzza RA, Urbanke C, Schuck P: Sedimentation velocity analysis of heterogeneous protein-protein interactions: Lamm equation modeling and sedimentation coefficient distributions c(s). Biophys J. 2005, 89 (1): 619-634. 10.1529/biophysj.105.059568PubMed CentralView ArticlePubMedGoogle Scholar
- Durchschlag H: Specific volumes of biological macromolecules and some other molecules of biological interest. Thermodynamic data for biochemistry and biotechnology. Edited by: Hinz H. 1986, 45-128. Berlin, Heidelberg: SpringerView ArticleGoogle Scholar
- Witte G, Urbanke C, Curth U: Single-stranded DNA-binding protein of Deinococcus radiodurans: a biophysical characterization. Nucleic Acids Res. 2005, 33 (5): 1662-1670. 10.1093/nar/gki310PubMed CentralView ArticlePubMedGoogle Scholar
- Konarev PV, Petoukhov MV, Volkov VV, Svergun DI: ATSAS 2.1, a program package for small-angle scattering data analysis. J Appl Cryst. 2006, 39: 277-286. 10.1107/S0021889806004699.View ArticleGoogle Scholar
- Wriggers W, Chacón P: Using Situs for the registration of protein structures with low-resolution bead models from X-ray solution scattering. J Appl Cryst. 2001, 34: 773-776. 10.1107/S0021889801012869.View ArticleGoogle Scholar
- Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem. 2004, 25 (13): 1605-1612. 10.1002/jcc.20084View ArticlePubMedGoogle Scholar
- DeLano WL: The PyMOL Molecular Graphics System. Version 1.2r3pre, Schrödinger, LLC. 2002Google Scholar
- Soding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21 (7): 951-960. 10.1093/bioinformatics/bti125View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.