reSeArCH Letter
26. Gottlieb, H. E., Kotlyar, V. & Nudelman, A. NMR chemical shifts of common
laboratory solvents as trace impurities. J. Org. Chem. 62, 7512–7515
(1997).
27. Holland, H. L., Andreana, P. R. & Brown, F. M. Biocatalytic and chemical routes
to all the stereoisomers of methionine and ethionine sulfoxides. Tetrahedron
Asymm. 10, 2833–2843 (1999).
28. Lavine, T. F. The formation, resolution, and optical properties of the
diastereoisomeric sulfoxides derived from l-methionine. J. Biol. Chem. 169,
477–491 (1947).
29. Koc, A., Gasch, A. P., Rutherford, J. C., Kim, H. Y. & Gladyshev, V. N. Methionine
sulfoxide reductase regulation of yeast lifespan reveals reactive oxygen
species-dependent and -independent components of aging. Proc. Natl Acad.
Sci. USA 101, 7999–8004 (2004).
30. Lherbet, C., Gravel, C. & Keillor, J. W. Synthesis of S-alkyl l-homocysteine
analogues of glutathione and their kinetic studies with γ-glutamyl
transpeptidase. Bioorg. Med. Chem. Lett. 14, 3451–3455 (2004).
31. Vertommen, D. et al. The disulphide isomerase DsbC cooperates with the
oxidase DsbA in a DsbD-independent manner. Mol. Microbiol. 67, 336–349
(2008).
32. Arts, I. S. et al. Dissecting the machinery that introduces disulfide bonds in
Pseudomonas aeruginosa. MBio 4, e00912–13 (2013).
33. Vizcaíno, J. A. et al. ProteomeXchange provides globally co-ordinated
proteomics data submission and dissemination. Nature Biotechnol. 30,
223–226 (2004).
colony-forming units (c.f.u.) per millilitre. The absolute c.f.u. at time-point 0 (used
as 100%) was ~108 cells per millilitre in all experiments. For strains CH194, CH196
and CH197, the same protocol was used with chloramphenicol (25μgml−1) and
arabinose (0.2%) added to the cultures.
SDS survival assays. Cells (MG1655 and BE107) were grown at 37°C with shaking
in 10ml of LB (in 100ml flasks). When cells reached A600nm≈0.8, 5ml samples
were transferred to conical polypropylene centrifuge tubes (50ml, Sarstedt) and
HOCl (2mM) was added. After 5min of incubation, samples were taken and
diluted in PBS buffer to ~2×103 cells per millilitre. Aliquots (100μl) were then
spread on LB agar plates containing SDS (1%). Colonies were counted the next day.
Data set construction and phylogenetic analyses. A non-redundant local pro-
tein database containing 1,342 complete prokaryotic proteomes available in NCBI
ried with the BlastP program (default parameters)39, using YedY (NP_416480) and
YedZ (NP_416481) of E. coli strain K-12 substrate MG1655 as a seed. Distinction
between homologous and non-homologous sequences was assessed by visual
inspection of each BlastP output (no arbitrary cut-off on the E value or score). To
ensure that we did not overlook divergent YedY or YedZ proteins, iterative BlastP
queries were performed using homologues identified at each step as new seeds.
The list of YedY and YedZ homologues is provided in Supplementary Data 1. The
retrieved sequences were aligned using MAFFT version 7 (default parameters40
;
34. Roberts, D. M. et al. Chemical synthesis and expression of a calmodulin
gene designed for site-specific mutagenesis. Biochemistry 24, 5090–5098
(1985).
35. Grimaud, R. et al. Repair of oxidized proteins. Identification of a new
methionine sulfoxide reductase. J. Biol. Chem. 276, 48915–48920
(2001).
Supplementary Data 2 and 3). Each alignment was visually inspected and manually
refined when necessary using the ED program from the MUST package41. Regions
where the homology between amino-acid positions was doubtful were removed
by using BMGE software (BLOSUM30 similarity matrix42).
For each homologue, the genomic context was investigated using MGcV
(Microbial Genomic context Viewer43). The domain composition and protein
and a gamma distribution with four categories47. On the basis of the resulting tree,
the subfamily containing the sequence from E. coli was identified and selected for
further phylogenetic investigations. The corresponding sequences were realigned
using MAFFT version 7. The resulting alignment was trimmed with BMGE as
previously described.
Maximum likelihood trees were computed using PHYML version 3.1 (ref. 48)
with the Le and Gascuel model (amino-acid frequencies estimated from the data
set) and a gamma distribution (four discrete categories of sites and an estimated
alpha parameter) to take into account variations in evolutionary rate across sites.
Branch robustness was estimated by the non-parametric bootstrap procedure
implemented in PhyML (100 replicates of the original data set with the same
parameters). Bayesian inferences were performed using MrBayes 3.2 (ref. 49) with
a mixed model of amino-acid substitution including a gamma distribution (four
discrete categories). MrBayes was run with four chains for one million generations
and trees were sampled every 100 generations. To construct the consensus tree, the
first 2,000 trees were discarded as ‘burn in’.
36. Tsvetkov, P. O. et al. Calorimetry and mass spectrometry study of oxidized
calmodulin interaction with target and differential repair by methionine
sulfoxide reductases. Biochimie 87, 473–480 (2005).
37. Cascales, E., Bernadac, A., Gavioli, M., Lazzaroni, J. C. & Lloubes, R. Pal
lipoprotein of Escherichia coli plays a major role in outer membrane integrity.
J. Bacteriol. 184, 754–759 (2002).
38. Miller, J. A Short Course in Bacterial Genetics Unit 3, 72–74 (Cold Spring Harbor
Laboratory Press, 1992).
39. Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein
database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).
40. Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software
version 7: improvements in performance and usability. Mol. Biol. Evol. 30,
772–780 (2013).
41. Philippe, H. MUST, a computer package of Management Utilities for Sequences
and Trees. Nucleic Acids Res. 21, 5264–5272 (1993).
42. Criscuolo, A. & Gribaldo, S. BMGE (Block Mapping and Gathering with
Entropy): a new software for selection of phylogenetic informative regions from
multiple sequence alignments. BMC Evol. Biol. 10, 210 (2010).
43. Overmars, L., Kerkhoven, R., Siezen, R. J. & Francke, C. MGcV: the microbial
genomic context viewer for comparative genome analysis. BMC Genomics 14,
209 (2013).
44. Finn, R. D. et al. Pfam: the protein families database. Nucleic Acids Res. 42,
D222–D230 (2014).
45. Petersen, T. N., Brunak, S., von Heijne, G. & Nielsen, H. SignalP 4.0:
discriminating signal peptides from transmembrane regions. Nature Methods
8, 785–786 (2011).
22. Baba, T. et al. Construction of Escherichia coli K-12 in-frame, single-gene
knockout mutants: the Keio collection. Mol Syst Biol 2, 2006.0008 (2006).
23. Bremer, E., Silhavy, T. J., Weisemann, J. M. & Weinstock, G. M. Lambda placMu:
a transposable derivative of bacteriophage lambda for creating lacZ protein
fusions in a single step. J. Bacteriol. 158, 1084–1093 (1984).
24. Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in
Escherichia coli K-12 using PCR products. Proc. Natl Acad. Sci. USA 97,
6640–6645 (2000).
46. Sonnhammer, E. L., von Heijne, G. & Krogh, A. A hidden Markov model for
predicting transmembrane helices in protein sequences. Proc. Int. Conf. Intell.
Syst. Mol. Biol. 6, 175–182 (1998).
47. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2 – approximately maximum-
likelihood trees for large alignments. PLoS One 5, e9490 (2010).
48. Guindon, S. et al. New algorithms and methods to estimate maximum-
likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59,
307–321 (2010).
25. Mandin, P. & Gottesman, S. A genetic approach for finding small RNAs
regulators of genes of interest identifies RybC as regulating the DpiA/DpiB
two-component system. Mol. Microbiol. 72, 551–565 (2009).
49. Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference
and model choice across a large model space. Syst. Biol. 61, 539–542
(2012).
© 2015 Macmillan Publishers Limited. All rights reserved