Geula Davidov1
Frank D. Müller2Jens Baumgartner3Ronit Bitton4
Damien Faivre3Dirk Schüler2
Raz Zarivach1*- 1Department of Life Sciences and the National Institute for Biotechnology in the Negev, Ben-Gurion University of the Negev, Beer Sheva, Israel
- 2Department of Microbiology, University of Bayreuth, Bayreuth, Germany
- 3Department of Biomaterials, Max Planck Institute of Colloids and Interfaces (MPI), Potsdam, Germany
- 4Department of Chemical Engineering, Ilse Katz Institute for Nanoscale Science and Technology, Ben Gurion University of the Negev, Beer-Sheva, Israel
Magnetotactic bacteria (MTB) are a diverse group of aquatic bacteria that have the magnetotaxis ability to align themselves along the geomagnetic field lines and to navigate to a microoxic zone at the bottom of chemically stratified natural water. This special navigation is the result of a unique linear assembly of a specialized organelle, the magnetosome, which contains a biomineralized magnetic nanocrystal enveloped by a cytoplasmic membrane. The Magnetospirillum gryphiswaldense MtxA protein (MGR_0208) was suggested to play a role in bacterial magnetotaxis due to its gene location in an operon together with putative signal transduction genes. Since no homology is found for MtxA, and to better understand the role and function of MtxA in MTBés magnetotaxis, we initiated structural and functional studies of MtxA via X-ray crystallography and deletion mutagenesis. Here, we present the crystal structure of the MtxA C-terminal domain and provide new insights into its sequence-structure relationship.
Introduction
Magnetotactic bacteria (MTB) are a heterogeneous group of aquatic microorganisms that share the ability to orient themselves along geomagnetic field lines. This ability is achieved by unique organelles—magnetosomes—that are synthesized by the bacteria for supposedly passive orientation along the magnetic field. This organelle is characterized by its ability to grow one magnetic nanocrystal (greigite, Fe3S4 or magnetite, Fe3O4) per vesicle under ambient conditions (Faivre and Schüler, 2008).
The magnetosomes are arranged into single or multiple chain-like structures that enable the cell to dynamically align along external magnetic fields, a behavior known as magnetotaxis (Schüler, 2008; Lefevre et al., 2014; Zhu et al., 2014). Magnetotaxis, in combination with aerotaxis and perhaps phototaxis, is thought to direct the swimming of cells toward growth-favoring microoxic zones at the bottom of chemically stratified natural waters (Frankel et al., 1997; Frankel and Bazylinski, 2006; Bennet et al., 2014; Popp et al., 2014). Magnetosome formation, magnetite biomineralization as well as magnetosome organization are controlled by a large set of soluble and integral membrane proteins, most of which are unique, largely encoded by genes clustered within a genomic magnetosome island (MAI) (Schübbe et al., 2003; Komeili, 2007; Richter et al., 2007; Kolinko et al., 2014; Nudelman and Zarivach, 2014).
Previously in the alphaproteobacterium Magnetospirillum gryphiswaldense, the MtxA protein (MGR_0208), which is in or attached to the magnetosome membrane, was identified by a comparative genomic study as “MTB-specific.” This means that it is highly conserved in other MTB, but absent from non-magnetic bacteria (Richter et al., 2007). mtxA is encoded within a gene cluster located near the MAI that is also partly conserved among MTB (Richter et al., 2007). Orthologous proteins of MtxA are present in MTB strains Magnetospirillum sp. AMB-1, Magnetospirillum magnetotacticum strain MS-1 and Magnetococcus marinus MC-1. These proteins exhibit high amino acid sequence identity (60%, 60%, and 50%, respectively) and have identical lengths (Matsunaga et al., 2005). All MtxA proteins contain a predicted signal peptide sequence and some of these proteins have a structural element that forms repeat motifs that are defined by possessing duplications of a basic sequence motif usually involved in the formation of a structural element of the protein (Andrade et al., 2001). One of the protein repeats is a 34-amino acid motif called tetratricopeptide repeats (TPR); each repeat of the TPR domain forms a helix-turn-helix structure that serves as a template for protein-protein interactions that can mediate the assembly of multi-protein complexes. The number of observed TPR units in proteins ranges from 3 to 16, which are generally arranged as tandem arrays (Zeytuni and Zarivach, 2012). TPR structural elements have confirmed the general folding pattern with variations in structure, as well as adaptations that provide the TPR-containing proteins with alternative modes of interaction with different binding partners (D'Andrea and Regan, 2003; Zeytuni and Zarivach, 2012).
In its gene cluster, mtxA is followed by a highly conserved gene with similarities to Pfamés adenylate cyclase protein family model and CHASE2 domain in MTB strains. This domain organization belongs to a widespread group of transmembrane receptors that function as sensors for monitoring environmental changes and regulatory circuit function in catabolite repression in microorganisms. Major types of these sensor proteins share conserved intracellular domains (Zhulin et al., 2003). A few of the proteins that participate in intracellular interactions have an immunoglobulin-like (Ig-like) domain which has a rich, β-strand domain that forms an antiparallel, β-sandwich with a topology analogous to an Ig constant domain (Bazan, 1990). Ig-like domain proteins are also found in bacteria and these bacterial immunoglobulin-like (Big) domains are found in various varied functional proteins (Bateman et al., 1996).
Because of its co-localization with putative signal transduction proteins within the same operon, mtxA was suggested to play a role in magnetotaxis (Richter et al., 2007). Another hypothesis of MtxA function is that it could be involved in the biomineralization of magnetite inside the magnetosome vesicle (Tanaka et al., 2006). However, in vitro assays suggested that, despite the possibility that the protein can bind to magnetite, it will likely not affect the crystallization of the mineral in vivo (Baumgartner et al., 2014). Since no obvious homology to any known protein could be found on the sequence level, we initiated genetic and structural studies of MtxA in order to better understand its role and possible involvement in magnetotaxis via structural identification of active folds/structures and via in vivo studies of an mtxA deletion mutant. Here, we present the in vivo cellular and in vitro structural analysis of the MGR_0208 deletion mutant lacking the signal peptide (MtxAΔ1−24).
Materials and Methods
Expression of the MtxAΔ1−24 Gene in Escherichia coli
The gene construct mtxAΔ1−24 was amplified from M. gryphiswaldense MSR-1 genomic DNA by PCR (oligonucleotides produced by MWG Operon) using KOD polymerase (Novagen). The amplified mtxAΔ1−24 gene was ligated between the KpnI and SacI restriction sites of the pET-51b(+)Ek/LIC vector (Novagen), giving rise to plasmid pET51bMtxAΔ1−24MSR1 (Baumgartner et al., 2014). In this construct, the mtxAΔ1−24gene was fused in-frame between an N-terminal Strep-tag and a C-terminal His-tag. MtxAΔ1−24 has a calculated mass of 34,774 Da.
E. coli BL21 strain cells harboring the MtxAΔ1−24MSR-1 plasmid were grown in LB medium containing ampicillin (50 mg/mL) at 37°C at 190 rpm.Isopropyl-d-thiogalactopyranoside (IPTG) was added to 0.5 mM when an optical density of 0.6 OD was measured at 600 nm, to induce protein expression for an additional 16 h at 20°C at 190 rpm. The cells were harvested by centrifugation at 7438 g for 10 min at 4°C.
Selenomethionine (SeMet) incorporation to MtxAΔ1−24 protein was carried out using the SeMet minimal medium (Guerrero et al., 2001). A 30 mL overnight culture was grown in LB medium at 190 rpm. The cells were harvested at 2629 g for 10 min and the cell pellet was used as an inoculum for 1 L of 1X M9 minimal medium supplemented with ampicillin (50 mg/mL) and other additives (Guerrero et al., 2001). 50 μg/mL l-selenomethionine, 100 μg/mL l-lysine, 100 μg/mL l-threonine, 100 μg/mL phenylalanine, 50 μg/mL l-leucine, 50 μg/mL l-valine, and 50 μg/mL l-proline were added when an optical density of 0.6 was measured at 600 nm. The culture was grown for 30 min before IPTG was added to induce protein expression for an additional 16 h at 20°C at 190 rpm. The cells were harvested by centrifugation at 7438 g for 10 min at 4°C.
Purification of MtxAΔ1−24
MtxAΔ1−24-expressing cells were suspended in lysis buffer (20 mM Tris-HCl pH 8.5, 300 mM NaCl, 20 mM imidazole, 0.02% Triton X-100) and incubated with DNase I (10 mg/mL) and a protease inhibitor cocktail [100 mM phenylmethylsulfonyl fluoride (PMSF), 1.2 mg/mL leupeptin and 1 mM pepstatin A] at a ratio of 1:1000 with binding buffer for 20 min at 4°C. The cells were then disrupted by two cycles in a French press pressure cell at 207 MPa. Cell debris was separated by centrifugation at 19,000 g for 90 min at 4°C and the soluble fraction was applied onto a 5 mL bed-volume homemade gravity Ni-NTA column (2.5 cm diameter, Econo-Column Chromatography Columns, Bio-Rad. Ni-NTA His Bind Resin; Lot M0063428, Novagen) that was pre-equilibrated with lysis buffer.
The protein was washed with four different wash buffers, as it leads to the best purified protein without contaminations: I (20 mM Tris-HCl pH 8.5, 1 M NaCl, 25 mM imidazole); II (20 mM Tris-HCl pH 8.5, 700 mM NaCl, 30 mM imidazole); III (20 mM Tris-HCl pH 8.5, 500 mM NaCl, 40 mM imidazole); IV (20 mM Tris-HCl pH 8.5, 300 mM NaCl, 50 mM imidazole) and was eluted with elution buffer (20 mM Tris-HCl pH 8.5, 200 mM NaCl, and 500 mM imidazole). The eluted protein was dialyzed against ion-exchange (IEX) buffer A (20 mM Tris-HCl pH 8.5, 50 mM NaCl) for 16 h at 4°C, after which the protein was applied onto a column of MonoQ 4.6/100 PE (GE Healthcare Biosciences) equilibrated with IEX buffer A. The protein was eluted with a linear gradient of 50–2000 mM NaCl in IEX buffer B.
The relevant protein peak was collected and concentrated using a Vivaspin-4 10,000 mWCO (Sartorius Stedim Biotech GmbH) and applied onto a column of HiLoad 26/60 Superdex 200 (GE Healthcare Biosciences) equilibrated with size-exclusion chromatography (SEC) buffer (20 mM Tris-HCl pH 8.5, 200 mM NaCl). Purified MtxAΔ1−24 was then concentrated to 47 mg/mL for crystallization, flash frozen in liquid nitrogen and stored at −80°C. Sample purity at this stage was analyzed by SDS-PAGE and protein identification was confirmed by tandem mass spectroscopy.
A cell harvesting and purification protocol similar to that used for native MtxAΔ1−24 protein was carried out with the SeMet MtxAΔ1−24 protein. Purified SeMet MtxAΔ1−24 was then concentrated to 52 mg/mL for crystallization, flash frozen in liquid nitrogen and stored at −80°C. Sample purity at this stage was analyzed by SDS-PAGE and the mass spectra obtained from the native MtxAΔ1−24 protein and the SeMet MtxAΔ1−24 protein were compared for differences (data not shown).
Circular Dichroism of MtxAΔ1−24
Circular dichroism (CD) measurements were conducted with a J750 Spectropolarimeter (Jasco Inc, Maryés Court, Easton, USA) equipped with a Pelletier device. MtxAΔ1−24 protein sample was prediluted to 0.2 mg/mL in buffer containing 20 mM Tris-HCl pH 8.5, 200 mM NaCl and measured with a 0.1 cm optical path Suprasil quartz cuvette (Hellma GmbH & Co., Müllheim, Germany). Spectra profiles of the samples were measured at a wavelength range of 190–240 nm at ambient temperature with bandwidth set to 1 nm, scan speed set to 10 nm/min and a time constant of 4 s. The thermal denaturation experiments of MtxAΔ1−24 were conducted by monitoring the dichroic absorption at a wavelength of 222 nm as a function of increased temperature varying from 25°C to 75°C at a heating rate of 1.0°C × min−1. The thermodynamic parameters associated with the temperature-induced denaturation were obtained by nonlinear, least-squares fitting analysis of the temperature dependence of CD, and a two-state denaturation process was assumed for the curve fit.
Small Angle X-ray Scattering (SAXS) Data Collection of MtxAΔ1−24
All protein samples for SAXS measurements were taken from the previously mentioned SEC purification to eliminate products of complex formation or aggregation. MtxAΔ1−24 was diluted with 20 mM Tris-HCl pH 8.5, 200 mM NaCl. SAXS measurements were performed at the French national synchrotron facility SOLEIL, on the SWING beamline. The incident beam energy was 12 keV. The sample-to-detector (Aviex CCD) distance was set to 1892 mm, covering a q-range of 0.004–0.7 Å−1. All experiments were temperature controlled at 25°C. Typically, 55 successive frames of 0.5 s each were recorded for the protein solution and its corresponding buffer. Each frame was first angularly averaged and the final spectrum and experimental error were obtained by averaging over all frames and subtracting the pure solvent spectrum from the sample spectrum. Intensities were scaled using the scattering of water (Carn et al., 2012).
SAXS Data Analysis and Envelope Model of MtxAΔ1−24
The radius of gyration (Rg) was evaluated using the Guinier approximation (Guinier and Fournet, 1955). The GNOM program was used to obtain the pair-distance distribution functions, the corresponding maximum dimension of protein complexes (Dmax) and to determine the value for Rg from the entire scattering profile (Svergun, 1992). Ab initio envelopes were generated by the program DAMMIN (Svergun, 1999) using atomic radii set to the dummy atom packing radius determined by DAMMIN without imposing symmetry operation. Five DAMMIN runs were performed for every sample and an averaged dummy ball model (DBM) was generated by DAMAVER (Volkov and Svergun, 2003). The generated envelope models (DBMs) were fit using the Coot software (Emsley and Cowtan, 2004) and visualized using PyMOL (DeLano, 2002).
Crystallization of MtxAΔ1−24
MtxAΔ1−24 was crystallized using the sitting-drop vapor-diffusion method at 20°C. 0.5 μL of MtxAΔ1−24 and 0.5 μL reservoir solutions were mixed to form the drop. The first crystallization trials were performed with MtxAΔ1−24 concentrated to 12, 15, 20, and 30 mg/mL at 20°C and 13°C in 20 mM Tris-HCl pH 8.5, 200 mM NaCl (Zeytuni et al., 2012). We traced these plates for up to 3 months without any positive crystal hit. The lack of protein crystals could be the result of the long, flexible Strep-tag and, as such, several options were devised to overcome the crystallization problems, such as the removal of the tag by cloning or the use of protease to remove the tag while crystallizing the protein. While trying this approach, we performed crystallization trials with 15 mg/mL MtxAΔ1−24 incubated with trypsin protease (trypsin from bovine pancreas, T8003 SIGMA) at 20°C in 20 mM Tris-HCl pH 8.5, 200 mM NaCl, at a molar ratio of 1:1000 or 1:4000 (trypsin protease:MtxAΔ1−24) in the drop. The initial crystallization conditions were examined using commercial screening kits from Hampton Research (Index, Crystal Screen I + II Screens), Molecular Dimensions (Structure screen I + II HT) and Rigaku (Wizard I + II screens). Crystal hits were observed in different conditions in the Index screen (Hampton Research) and Wizard screen (Rigaku). The same protocol of crystallization trial with trypsin protease was performed for the SeMet MtxAΔ1−24 protein. Crystallization of the SeMet MtxAΔ1−24 protein was with the Index commercial screening kit (Hampton Research). Crystal hits were noticed in different conditions in the Index screen; condition No. 87 produced the same clusters of tiny needle crystals that were obtained using a reservoir containing 0.2 M sodium malonate pH 7.0 and 20% PEG 3350 after 24 h incubation. We then harvested the SeMet MtxAΔ1−24 crystals directly from the screen and flash froze them in liquid nitrogen without any cryo-soaking, due to the high concentration of the PEG 3350 that acted as cryo-protectant. One of the SeMet MtxAΔ1−24 crystals in the condition No. 87 diffracted to 2 Å resolution.
Diffraction Data Collection and Structure Determination of MtxAΔ1−24
The native and SeMet Crystals were harvested and flash-cooled in liquid nitrogen without addition of a cryoprotecting solution. Diffraction data for native MtxAΔ1−24 were collected on beamline ID23-2 at the European Synchrotron Radiation Facility (ESRF, Grenoble, France), resulted in low resolution diffraction (data not showed). Diffraction data for the SeMet MtxAΔ1−24 protein were collected on beamline ID14-4 at the ESRF, which is equipped with an ADSC Q315r mosaic CCD detector. Data collection was performed at 100 K. For the SeMet MtxAΔ1−24 crystal form we scanned the crystal for existence of Selenium signal and to obtained f′ and f″ values (–8.60, 6.92 respectively). A peak data set was measured at a wavelength of 0.98 Å for a total of 267 frames with an oscillation range of 1° and an exposure time of 0.224 s per image. The crystal-to-detector distance was 347.46 mm. Data reduction and scaling were performed with HKL-2000 (Otwinowski and Minor, 1997). This data was used for Se-SAD phase calculation. Phases were obtained by the automatic data processing which has been available and in use on all Joint Structural Biology Group (JSBG) MX beamlines (Monaco et al., 2013). Phase statistics are: correlation coefficient CC(E) = 0.399, figure of merit 0.7551acentric/0.09483centric and phasing power is 1.514. The output map and protein sequence were entered into the auto-build in ARP/wARP web service; http://www.embl-hamburg.de/ARP/ (Langer et al., 2008). The program placed correctly 230 out of the 650 residues and water molecules were added manually. The final model was built by Coot and refined in REFMAC (Vagin et al., 2004; Emsley et al., 2010). Structural figures were prepared with PyMOL (DeLano, 2002).
mtxA Deletion
A markerless in-frame deletion of mtxA (Mgr_0208) was carried out essentially as described (Raschdorf et al., 2014). Briefly, 1 kbp up- and downstream (including the first and the last four codons) of mtxA were amplified from genomic DNA, purified and fused by overlap-extension PCR. The resulting 2 kbp DNA fragment was purified, digested with BamHI and PstI, ligated into pORFM GalK blu and used to transform competent E. coli DH5α cells. Kanamycin resistant, white colonies were selected to isolate the resulting vector pIM01 and DNA sequencing was performed to verify absence of mutations within the homologous regions. Subsequently, the vector was transformed into E. coli BW29427. Transformation of M. gryphiswaldense with pIM01 was carried out by biparental conjugation as described (Schultheiss and Schüler, 2003; Ullrich and Schüler, 2010). Recombinant, merodiploid clones were isolated from kanamycin-supplemented modified flask standard medium (FSM) plates and screened for up or downstream integration of the plasmid. Counterselection was then performed on FSM plates supplemented with 0.5% (w/v) galactose. Successful deletion of mtxA was confirmed by PCR-screening with oligonucleotide primers specific to regions adjacent to the cloned homologous sections.
Phenotypic Characterization of the mtxA Mutant
Swarm assays were carried out on swim agar plates as described (Popp et al., 2014). For transmission electron microscopy (TEM) imaging cells were grown at 25°C under microaerobic conditions to an OD565 of 0.1, fixed in formaldehyde (1%), concentrated, adsorbed onto carbon-coated copper mesh grids, and washed three times with particle-free water. Samples were viewed and recorded with a Morgagni 268 microscope (FEI, Eindhoven, the Netherlands) at an 80-kV accelerating tension. Magnetic responses (Cmag) of exponentially growing M. gryphiswaldense cultures were measured photometrically at 565 nm as reported previously (Schüler et al., 1995).
Results
Protein Sequence Analysis of the Putative Magnetotaxis Protein MtxA
In our first attempt to characterize MtxA, we ran a multiple sequence alignment with different servers. In Figure 1, we present the output from the ClustalW2 server; http://www.ebi.ac.uk/Tools/msa/clustalw2 (Larkin et al., 2007) and ESPript 3.0; http://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi (Gouet et al., 2003). Similarity search analyses revealed high homology at the amino acid sequence level in five different organisms, of which two are not MTB. All mtxA genes in all identified organisms represent uncharacterized proteins with highly conserved residues. The sequence comparison highlights that MtxA shares 61% identity with the M. magneticum AMB-1 amb2230 protein, 58% identity with Magnetospirillum sp. SO-1 H261_01492 protein and Phaeospirillum molischianum DSM 120 PHAMO_340145 protein, 57% identity with Phaeospirillum fulvum MGU-K5 K678_00055 protein, and 51% identity with M. sp. MC-1 Mmc1_3696 protein (Figure 1). All the other servers used (http://www.uniprot.org/align, http://blast.ncbi.nlm.nih.gov/Blast.cgi, http://www.st-va.ncbi.nlm.nih.gov/tools/cobalt) yielded similar results.
Figure 1. Sequence alignment of magnetotaxis protein MtxA (MSR-1). Amino acid sequence alignment of the Magnetospirillum gryphiswaldense MSR-1 magnetotaxis protein MtxA MGR_0208 (MSR-1) with the Phaeospirillum molischianum DSM 120 PHAMO_340145 (DSM 120), the Phaeospirillum fulvum MGU-K5 K678_00055 (MGU–K5), the Magnetospirillum sp SO-1 H261_01492 (SO-1), the Magnetospirillum magneticum AMB-1 amb2230 (AMB-1) and the Magnetococcus sp. MC-1 Mmc1_3696 (MC-1). Strictly conserved residues are highlighted with a red background and highly homologous residues are boxed. The secondary structural elements of MSR-1 defined by the analysis of the structure by using the PSIPRED server are indicated as red coils for α-helices and green arrows for β-strands. The purple arrow indicates for the end of the signal peptide.
To further analyze the protein sequence, we ran secondary-structure prediction on the PSIPRED server; http://bioinf.cs.ucl.ac.uk/psipred (Buchan et al., 2010) (Figure 1). The secondary structure prediction indicated that MtxA has four different structural regions. The first 24 a.a. (blue arrow) that were predicted to fold mainly as an α-helix were previously recognized as a signal peptide sequence (Richter et al., 2007). The second region (25-110 a.a., purple arrow) is mainly predicted to be a β-strand region. The third region (110-230 a.a., black arrow) is predicted to fold as alternating α-helices and β-strands. The last region (231-313 a.a., pink arrow) is predicted to fold mainly as α-helices. The prediction did not find any transmembrane domain or long unstructured regions (Figure 1). We also ran a sequence comparison of MtxA against the Protein Data Bank (PDB), which shows only sequence similarities to very small segments within the MtxA sequence, with ~24–35% sequence identity (data not shown).
Phenotype of the mtxA Deletion Mutant
Characterization of the mtxA mutant under standard culture conditions revealed a phenotype virtually indistinguishable from the wild type (wt). TEM analysis indicated that the mutant cells formed regular, wt-sized and -shaped crystals that are aligned in wt-like chains (Supporting Figure 6). Consistently, the mutant strain exhibited a wt-like magnetic response (not shown). Microscopic inspection showed that the mutant cells were motile, and a swarm assay on swim agar plates exposed to a magnetic field revealed that mtxA mutant cells aligned and moved parallel to the field as the wt, as indicated by distortion of the aerotactic swim ring in semisolid agar (Popp et al., 2014) (Supporting Figure 7). However, although we were unable to discern phenotypic differences between the wt and the mtxA mutant under the common standard growth conditions, this does not entirely preclude that the mtxA deletion might cause more subtle effects under some untested conditions.
Characterization of MtxAΔ1−24
To study the structure-function of MtxA and to avoid the secretion of MtxA to the periplasm, we sub-cloned a deletion mutant that includes Strep- and His-tags for detection and purification (MtxAΔ1−24). MtxAΔ1−24 was expressed as a soluble protein. For protein purification, we started with Ni-NTA affinity chromatography followed by ion exchange (MonoQ) and SEC (Superdex 200) (Figure 2A). Comparing MtxAΔ1−24 to the molecular size marker in the size-exclusion chromatogram indicates that MtxAΔ1−24 eluted as a monomer. SDS-PAGE analysis of the corresponding size-exclusion peak showed that the protein was highly purified and stable (Figure 2A).
Figure 2. MtxAΔ1−24 size-exclusion chromatography and SAXS envelope model. (A) Size exclusion (Superdex 200 HiLoad 26/60) chromatography analysis of the MtxAΔ1−24 (red) and chromatographic separation and calibration curve for the standard proteins (black). The standard proteins are: I: Thyroglobulin (Mr = 669,000), II: Aldolase (Mr = 158,000) and III: Ovalbumin (Mr = 44,000) (Gel Filtration Calibration Kit HMW code no. 28-4038-42, GE). The elution volume of MtxAΔ1−24 (M) is at ~210 mL. Comparison with elution volumes of the standard proteins suggests that the protein is present as a monomer in solution (the molecular weight of MtxAΔ1−24 is 34.77 kDa). The eluted protein was separated by 15% SDS-PAGE gel and stained with Coomassie blue. (B) Small-angle X-ray scattering data of MtxAΔ1−24 in solution shown as individual data points. (C) The ab initio bead model of MtxAΔ1−24 calculated using DAMMAVER with the MtxA structure inside. Dummy atoms are shown as blue spheres, surrounded by a molecular surface and MtxA (orange) shown in cartoon (a.u., arbitrary unit).
To compare the predicted secondary structure elements with MtxAΔ1−24, we measured its (CD) curve at room temperature (Supporting Figure 1A). The observed double-minimum at 205 nm and 225 nm suggests contributions of β-strand (40.6%) and α-helical (3.4%) conformations, as estimated by the K2D3 analysis algorithm (Louis-Jeune et al., 2012). This value corresponds with our secondary structure prediction for a higher percentage of β-strand content. We have also analyzed the thermal stability of MtxAΔ1−24 by measuring the CD temperature dependence (25–75°C) at 222 nm (Supporting Figure 1B). The sigmoidal melting curve with an extrapolated melting temperature (Tm) of 50 ± 2°C further supports our assumption that MtxAΔ1−24 is a stable and soluble protein.
To further characterize the monomeric form of MtxAΔ1−24 in solution, we used small angle X-ray scattering (SAXS). The scattering data was collected at three different protein concentrations of 0.5, 1, and 2 mg/mL (Figure 2B). At these concentrations the scatterings are quite similar, indicating that the molecular dimensions are independent of protein concentration over this concentration range. To obtain a three dimensional model of MtxAΔ1−24, we generated a dummy-ball model (DBM) from the 2 mg/mL SAXS data. The model displays an elongated structure divided into large and small connected ellipsoids (Figure 2C). The long ellipsoid axes dimensions of the large and small ellipsoids are 6 and 3 nm, respectively.
To obtain a high resolution structure, we started with crystallization trials that were performed as described in M&M for native and SeMet MtxAΔ1−24(Supporting Table 1 for the crystal hits). Native crystals did not diffracted beyond 2.8 Å so we focused on crystal hits from the Index condition No.87 that had the largest non-multilayer distinct well faceted SeMet crystals. We harvested the SeMet MtxAΔ1−24crystals directly from the screen and flash froze them in liquid nitrogen without any cryo-soaking, due to the high concentration of the PEG 3350 that acted as cryo-protectant.
One of the SeMet MtxAΔ1−24 crystals diffracted to 2 Å resolution. Analysis of the diffraction data fulfilled the systematic absence rules of space group P212121. The unit cell parameters are a = 40.34 Å, b = 88.95 Å, c = 95.40 Å (Table 1). The crystal unit cell and space group parameters are equal to the parameters that had been obtained from the native MtxAΔ1−24 (data not shown). Since protease was used as part of the crystallization, MtxAΔ1−24 crystals are composed of some of the protein domains that result from the cleavage. The cloned MtxAΔ1−24 protein consists of 325 amino acids with an N-terminal Strep-tag and a C-terminal His-tag, whereas the full MtxA has 313 amino acids. Mass spectrometry analysis of the trypsin treated MtxAΔ1−24 exhibits a high intensity peak around 20.52 kDa and a less intense peak around 41.31 kDa that fits the MtxAΔ1−141 20.67 kDa size (Supporting Figure 2). Furthermore, the residues seen in the electron density maps are L136-S311 (L142-Q313, MtxAΔ1−141). The Matthewsé coefficient calculation that was performed based on the unit cell parameters and the molecular weight for the L136-S311 residues is with a probability of 98%, assuming the presence of two monomers per asymmetric unit and a solvent content of 45.3%. Both are within the normal range of values observed for soluble protein crystals (Matthews, 1968; Winn et al., 2011).
Table 1. Data collection and refinement statistics of SeMet MtxAΔ1−142.
The Overall Fold of MtxAΔ1−141
The crystal structure of MtxAΔ1−141 was determined by the automatic data processing server on beamline ID14-4 and auto-built with ARP/wARP and Coot. The structure was refined to an R factor of 18.5 and Rfree of 22.7% (Table 1) and was deposited to the PDB (4Z29) (map quality is shown in Supplementary Figure 8). The overall structure of MtxAΔ1−141 has two distinguished domain architectures (Figure 3A): I- Ig-like fold (MtxA-Big) in orange, and II- TPR fold (MtxA-TPR) in green. The overall structure of MtxAΔ1−141 was back-fitted on the generated DBM from the SAXS data (Figure 2C). The structure displays a good fit to the longest ellipsoid structure.
Figure 3. MtxAΔ1−141 crystal packing and overall structures. (A) The representative crystal packing of MtxAΔ1−141 structure contains a two-domain architecture: MtxA142−226 (MtxA-Big) and MtxA227−313 (MtxA-TPR) motifs. (B) The overall β-sandwich fold of the MtxA-Big structure; the seven β-sheets are labeled (A–F). The conserved amino acids are in stick representation. (C) Superimposition overlay of MtxA-Big with three of the Dali serverés alignments in cartoon representation reveals the high structural similarity: MtxA-Big in orange, 1J8K in cyan, 1OWW in deep purple and 1TDQ in forest. The molecules are shown in a 90° rotation view. (D) The overall fold of the (MtxA-TPR) structure contains two sequential TPR motifs. The conserved amino acids are represented as sticks. (E) Superimposition overlay of MtxA-TPR with some of the Dali server-aligned structures in cartoon representation reveals the high structural similarity: MtxA-TPR in green, 2Q7F in cyan, 3ASD in pink and 2AVP in blue, with enlargement of the high structural similarity area. All the structural images were prepared using PyMOL.
Electrostatic surface characteristics of the MtxAΔ1−141 show several electrostatic charged patches and a main hydrophobic patch (Supporting Figure 4). The protein exhibits a highly negatively electrostatic potential charge at the bottom of the C-terminal MtxA-TPR domain (left, bottom). The hydrophobic patch covers the whole face created by the TPR and the Ig-like domains (top left, Supporting Figure 4).
The MtxAΔ1−141 structure displays two protein monomers in the asymmetric unit (Supporting Figure 5). The two monomers tightly contact each other with a buried surface area of 418 Å2. The interface between the two monomers is stabilized by a polar interaction network that involves hydrogen bonds, and charged and hydrophobic interactions (Supporting Figure 5). The hydrophobic interactions between the two monomers are comprised of highly conserved residues in MtxA, residues 282-287 “MRPLLI” (Figure 1), located at the Surface of the TPR domain. The surface around this position has a negative electrostatic charge (Supporting Figure 4).
The Fold of MtxA-Big and MtxA-TPR in MtxAΔ1−141 Structure
The first domain—MtxA-Big: L142-L226 (MtxA142−226)—has a common Ig-like fold consisting of seven β-strands forming two antiparallel β-sheets, which pack against each other forming a β-sandwich (Figure 3B). In our structure monomer B, MtxA-Big domain has a higher B-factor than monomer A (37.4/20.3 Å2 respectively, calculated using all atoms, Supporting Figure 5B). Based on the crystal structure, monomer B has less surface area at the crystal contacts than monomer A, which may lead to a higher B-factor due to thermal instability. Alternatively, the higher B-factor may indicate for the MtxA-Big domain flexibility, which is part of the protein function. To see whether the MtxA-Big domain has similarities to other structures with similar topology available in the PDB, we used the Dali server; http://ekhidna.biocenter.helsinki.fi/dali_server/ (Dietmann et al., 2001). One of the best-generated fits with the MtxA-Big domain is fibronectin EDA (PDB ID code 1J8K, Z score = 8.4). Fibronectin type III (Fn-3) structures is a ~90-amino acid domain that forms seven β-strands in conserved regions of the chain (Niimi et al., 2001). These strands are folded into antiparallel β-sandwiches with a topology that is similar to immunoglobulin constant domains (Bazan, 1990). The MtxA-Big domain forms a β-sandwich fold with four β-strands on one side and three on the other side with no disulfide bonds, similar to some members of the Fn-3 family (Potts and Campbell, 1994). The first β-sheet (orange in Figure 3B) comprises strands A, B and E, whereas the second β-sheet (blue in Figure 3B) comprises strands C, D, F, and G. We overlapped MtxA-Big domain with three of the Dali matches: Tenascin-R (PDB ID code 1TDQ-A, Z score = 8.3) (Lundell et al., 2004), fibronectin EDA (PDB ID code 1J8K, Z score = 8.4) (Niimi et al., 2001) and Fibronectin first type III module (PDB ID code 1OWW, Z score = 7.4) (Gao et al., 2003) (Figure 3C). The structural superimposition shows that all adopt the same β-sheet arrangements, topology and conformations with RMS deviation of 1.4–1.57 Å, despite their low sequence identity (Supporting Table 2). Since the sequence identity is very low, we ran multiple sequence alignment in the ClustalW2 server to disclose significant highly conserved residues in Fn-3 structures (Supporting Figure 3A). The sequence alignment brought to light that MtxA-Big has the same highly conserved residues as in Fn-3 members: W161 (in β-sheet B), Y170 (in β-sheet C), and L196 (located in a loop between β-sheet E and β-sheet F) (Hoxha and Campion, 2014). We also found that P198 and G199 residues are highly conserved, as in some of the Fn-3 family members that are part of a structural conformation called the “tyrosine corner” (Hemmingsen et al., 1994). The “tyrosine corner” is a consensus sequence in Fn-3members: L-X-X-G-X-X-Y. MtxA-Big has only partial part of the consensus, L196, P198, G199 residues are conserved, but MtxA-Big lacks the sixth amino acid and instead of the tyrosine it has a histidine residue. In Fn-3 members, the Y and the L are close to each other (β-sheet E through β-sheet F) and the Y forms hydrogen bonds with a second a.a. in the consensus sequence above. The conserved H201 in MtxA-Big forms a hydrogen bond with a second a.a. (T197) and a hydrogen bond with the conserved D175 (Figure 1). This interaction probably preserves the β-sheet sandwich scaffold since MtxA-Big has shorter β-sheets C and D. Looking at the sequence alignment of MtxA, it discloses that there are another two residues that are well conserved: G186 and T222. G186 is in a loop between β-sheet D and β-sheet E and T222 at the C-terminal of the domain (β-sheets G).
The second domain, MtxA-TPR S227-N313 (MtxA227−313), contains five anti-parallel α-helices-and-turn motifs folded as tetra-trico-peptide repeat (TPR) motifs (Figure 3D). Again, we used Dali server for a structural superposition comparing this domain with proteins at the PDB. One of the best-generated fits with the MtxA-TPR domain is MamAR50E (PDB ID code 3ASD-A, Z score = 10.4) (Zeytuni et al., 2011). The TPR is a structural motif based on consensus sequence defined by a pattern of small and large hydrophobic amino acids, that adopts a helix-turn-helix fold that creates the repeating antiparallel α-helices (Zeytuni and Zarivach, 2012). Superimposition of MtxA-TPR domain structure with three Dali matched structures—synthetic consensus TPR protein (PDB ID code 2AVP-A, Z score = 10.4) (Kajander et al., 2007), MamAR50E (PDB ID code 3ASD-A, Z score = 10.4) (Zeytuni et al., 2011), and YrrB protein (PDB ID code 2Q7F-B, Z score = 11.3) (Han et al., 2007)—shows that they all adopt exactly the same TPR fold (Figure 3E) but share low sequence identity (Supporting Table 2). The structural superimposition of MtxA-TPR with 2AVP-A, 3ASD-A and 2Q7F-B results with RMS deviation of 1.2–1.39 Å with similar conformations of helix-turn-helix (Figure 3E and Supporting Table 2). We ran multiple sequence alignment in ClustalW2 server (Supporting Figure 3B) to search for the known conserved residues in the MtxA-TPR motifs because the result from the TPRpred server; http://toolkit.tuebingen.mpg.de/tprpred (Biegert et al., 2006) determined that there are no significant repeats in MtxA. The MtxA-TPR sequence alignment indicates that some of the consensuses pattern of conserved residues exists: G253 (position 8 in the N-terminal motif), G260 (position 15 in the N-terminal motif), A265 (position 20 in the N-terminal motif), and Y290 (position 11 in the N-terminal motif of the second repeat). The MtxA-TPR also keeps the pattern of hydrophobic amino acids, although it has variations in the consensus positions that preserve the motif of a helix-turn-helix (D'Andrea and Regan, 2003). Furthermore, one of these changes is a highly conserved residue in all the MTB strains. The difference in the consensus P32 position (between TPR motifs) (Zeytuni and Zarivach, 2012) is K277 in the MtxA-TPR. This K277 keeps the highly conserved sequence of “K-D-D-N” that holds the TPR-turn in MtxA.
MtxA-Big and MtxA-TPR contact each other extensively in a parallel orientation with a buried surface of 558 Å2. The interface between the two folds is stabilized by a polar and hydrophobic interaction network (Figure 4). The polar network involves polar, hydrophobic and charge interactions. The hydrogen bonds are: T151 with G260, T155 with D231, E162 with R296, E187 with R270, and R190 with D267 (Figure 4A). The amino acids that interact with water molecule and form hydrogen bond network are: R145, A147, S148, D149, T151, E162, W225, E230, T262, V263, D292, L293, and K294 (Figure 4B). The hydrophobic interaction network is comprised of three main regions (Figure 4C). The first region is composed of hydrophobic interactions of Y156 with the side chains of V261, A264, and L256. Another hydrophobic interactionsé network is the side chain of I150 with the side chain of V263, V160, and L293. The last network embraces hydrophobic interactions between the side chain of V160, M188, A265, M266, L295, and L298.
Figure 4. Polar and hydrophobic interaction network of MtxA-Big with MtxA-TPR. (A) “Close up” view of the hydrogen-bond interaction network (blue dotted lines) between the residues from each fold backbone: T151, T155, E162, E187, and R190 (MtxA-Big; orange cartoon; side chains shown as sticks) with D231, G260, D267, R270, and R296 (MtxA-TPR; green cartoon; side chains shown as sticks). (B) Overview of water-mediated interactions (red dotted lines) between the side chains (shown as sticks) and water (red spheres): R145 with K294, R145 with E162, A147 with L293, S148 with D292, D149 with D292, D149 with T262, T151 with V263, and W225 with E230. (C) Surface representation of hydrophobic interacting residues (shown as sticks and surface).
Discussion
In this study, we present the analysis of MtxA, a protein that has been suggested to play a role in magnetotaxis (Schüler, 2008; Zhu et al., 2014). However, considering the lack of a phenotype of the mtxA deletion mutant we conclude that MtxA is essential for magneto-aerotaxis in M. gryphiswaldense, although it might have more subtle functions under certain conditions not tested in our standard assays. MtxA shares homology with orthologous proteins from five different organisms, two of which are not MTB. In the five organisms, the proteins are highly conserved, have not been characterized yet and no functional domain has been associated with them. The amino acid sequence alignment suggests that the conservation covers the whole protein, aside from two major exceptions: amb2230 protein (AMB-1) that lacks the N-terminal (1-63 a.a.) and Mmc1_3696 (MC-1) that has lower identity that causes the sequence alignment to be more homologous than strictly conserved. Our structure prediction of MtxAΔ1−24 indicates that the protein has a very diverse fold but has more β-strands than α-helices, as suggested by the CD results. MtxAΔ1−24 was efficiently expressed in E. coli and is highly soluble as monomers in solution. Additionally, the SAXS DBM demonstrates that MtxAΔ1−24 behaves as an extended ellipsoid with two internal domains of which one is significantly larger than the other. This can indicate that MtxAΔ1−24 is composed of two smaller domains connected by a flexible linker, which can explain the lack of protein crystals for the full length MtxA.
It has been shown previously that flexible proteins can be difficult to crystallize and limited proteolysis can be used to obtain crystals (Dong et al., 2007). This further supports our result that MtxAΔ1−24 crystals appeared only when we crystallized the protein in the presence of trypsin protease, most likely due to its cleavage. The proteolysis with trypsin yielded a smaller and stable domain that generated crystals of excellent quality, which diffracted to 2 Å resolution. We determined the crystal structure of MtxAΔ1−141 and compared it to the initial secondary structure prediction. Overall, the secondary structure was divided correctly into two distinct folded domains: Ig-like and TPR-like, these domains are well established and no server had predicted them in MtxA protein. When looking into the specific secondary structure elements, the secondary structure prediction is deviated from our determined structure further indicating for the inaccuracy of the prediction due to limited sequence similarity of MtxA to other proteins. The MtxA-Big domain adopts the same β-sheet arrangements, topology and conformations as Ig-like proteins, despite the low sequence identity. Furthermore, it has the same highly conserved residues as in Fn-3 protein members (Hoxha and Campion, 2014). Comparing MtxA to Tenascin-R (PDB ID code 1TDQ) indicates that they have a similar structure to Fn-3 protein members and are extracellular proteins. Tenascin-R is eukaryotic member of the tenascin family of extracellular matrix glycoproteins and although it is restricted to the nervous system and affects cell migration, adhesion and differentiation, it has no clinical consequences in knock-out animal model (Anlar and Gunel-Ozcan, 2012). The lack of phenotype is similar to the mtxA mutant as we were unable to discern phenotypic differences under the common standard growth conditions.
Another important structural motif of MtxA-Big is the “tyrosine corner” (Hemmingsen et al., 1994): MtxA protein has histidine instead of the tyrosine (in β-sheet F). The H201 keeps the hydrogen bond with the second a.a. in the consensus sequence (L-X-X-G-X-X-Y) and forms another hydrogen bond with the conserved D175 (in β-sheet C), whereas in members of Fn-3 family there is a hydrophobic interaction of the tyrosine with the last a.a. of β-sheet C. Additionally, all MtxA proteins contain the consensus sequence L-X-P-G-X-H (196-201 a.a.), the “histidine corner” and the highly conserved D175 and 176G in the protein. The “histidine corner” is conserved in human SOD (PDB ID code ISPD) (Deng et al., 1993) and this arrangement makes histidine corner more difficult to accommodate than tyrosine in barrel interiors (Hemmingsen et al., 1994). One of the sites that bind copper in human CuZnSOD protein is the H-V-H at the sequence L-X-X-G-X-H-X-X-H-V-H (Hart et al., 1998); since we assume that MtxA protein does not bind copper, this His corner may indicate for a binding site to other ligands or ions.
The MtxA-TPR domain has a well-known TPR motif although the TPRpred server determined that there are no significant repeats in MtxA protein. MtxA-TPR contains five anti-parallel α-helices-and-turn motifs folded as TPR, and contain the hydrophobic amino acids pattern with some of the motifsé conserved residues (D'Andrea and Regan, 2003). The consensus P32 located between the TPR motifs (Zeytuni and Zarivach, 2012) is modified to K277 in all MtxA proteins (apart from MC-1, P277). This K277 starts a highly conserved sequence in all MtxA proteins: “KDDND/E” (277-281 a.a.). This conserved “zone” exhibits a negative electrostatic potential charge at the interface between MtxA-TPR and MtxA-Big domains (Supporting Figure 3, right, top). Followed by this sequence is another highly conserved patch of MtxA, residues 282-287 “MRPLLI,” that are involved in the hydrophobic interactions between the two monomers in the asymmetric unit.
The hydrophobic patch that covers the whole face created by the TPR and the Ig-like domains may indicate for an interface for interacting with other proteins, since TPR and Ig-like folds are known to be part of protein-protein interactions (Remaut and Waksman, 2006; Zeytuni and Zarivach, 2012). In addition, the C-terminal of MtxA-TPR domain has a very negatively electrostatic potential charge caused by the following amino acids: D248, D249, D278, D279, D281, E302, E308, and E311 (a negative amino acid array). This residue composition may reveal a negatively charged concave surface in MtxA similar to the negatively charged concave surface in YrrB (PDB ID code 2Q7F) (Han et al., 2007).
Since the full MtxA structure, based on our SAXS model, displays an elongated structure divided into large and small connected ellipsoids, we believe that the large ellipsoid is the MtxAΔ1−141 and small ellipsoid is probably the missing MtxA25−141 domain. Since our results did not implicate MtxA function and since we missed a full domain, we wanted to obtain its full structure as a means to obtain a predicted function. For that, we submitted MtxA to the CASP competition (Critical Assessment of Techniques for Protein Structure Prediction) (Kryshtafovych et al., 2014). From the CASP results (T0828) we found two “close” predictions out of the full structure list, the first from nns_TS1 group and the second from QUARK_TS2 group. Structural superimpositions of MtxAΔ1−141 with nns_TS1 and QUARK_TS2 resulted in structural superposition with RMS deviation of 1.54 Å and 1.40 Å, with 58 and 72 superposed Cα atoms at the MtxA-Big and MtxA-TPR fold, respectively. This signifies that the MtxA protein has a unique fold and may have an important role in signaling at the periplasm.
In summary, we have applied a wide range of methodologies in order to characterize MtxAΔ1−24 in solution. We found that MtxA is a monomeric, 34.7 kDa protein with an extended ellipsoid shape containing at least two internal domains. We have determined the MtxAΔ1−141 structure and disclosed its unique fold. The MtxA protein probably has an importance in periplasmic or extracellular interactions that could act as a sensor or mediating interaction with other proteins or even activating signals in the bacteria.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank ESRF (Grenoble, France) beamline ID23-2 and ID14-4 for assistance during data collection. We acknowledge SOLEIL for provision of synchrotron radiation facilities and we would like to thank Dr. Javier Pérez for assistance in using beamline SWING. We thank Dr. Anat Shahar from the Macromolecular Crystallography Research Center (MCRC), Ben-Gurion University of the Negev, Beer Sheva, Israel (BGU). We would like to thank Dr. Mark Karpasas from the Ilse Katz Institute for Nanoscale Science and Technology center, BGU, Beer Sheva, Israel, for the MALDI-TOF analysis. We would like to thank Dr. Tamar Kurzion-Zilbermann for the CD analysis from the Department of Chemistry, BGU, Beer Sheva, Israel. We are grateful to Isabelle Mai for performing mutagenesis experiments and Katharina Silbermann for phenotypic characterization of the mtxA mutant. This was supported by Deutsche Forschungsgemeinschaft (Grants Schu 1080/12-1 and 16-1), the Israel Science Foundation, EMBO Young Investigators Grant and CMST COST Action CM0902. Damien Faivre acknowledges support from the European Research Council (Starting Grant MB2 n°256915).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmolb.2015.00025/abstract
References
Andrade, M. A., Perez-Iratxeta, C., and Ponting, C. P. (2001). Protein repeats: structures, functions, and evolution. J. Struct. Biol. 134, 117–131. doi: 10.1006/jsbi.2001.4392
Anlar, B., and Gunel-Ozcan, A. (2012). Tenascin-R: role in the central nervous system. Int. J. Biochem. Cell Biol. 44, 1385–1389. doi: 10.1016/fiocel.2012.05.009
Bateman, A., Eddy, S. R., and Chothia, C. (1996). Members of the immunoglobulin superfamily in bacteria. Protein Sci. 5, 1939–1941. doi: 10.1002/pro.5560050923
Baumgartner, J., Carillo, M. A., Eckes, K. M., Werner, P., and Faivre, D. (2014). Biomimetic magnetite formation: from biocombinatorial approaches to mineralization effects. Langmuir 30, 2129–2136. doi: 10.1021/la404290c
Bazan, J. F. (1990). Structural design and molecular evolution of a cytokine receptor superfamily. Proc. Natl. Acad. Sci. U.S.A. 87, 6934–6938. doi: 10.1073/pnas.87.18.6934
Bennet, M., McCarthy, A., Fix, D., Edwards, M. R., Repp, F., Vach, P., et al. (2014). Influence of magnetic fields on magneto-aerotaxis. PLoS ONE 9:e101150. doi: 10.1371/journal.pone.0101150
Biegert, A., Mayer, C., Remmert, M., Soding, J., and Lupas, A. N. (2006). The MPI Bioinformatics Toolkit for protein sequence analysis. Nucleic Acids Res. 34, W335–W339. doi: 10.1093/nar/gkl217
Buchan, D. W., Ward, S. M., Lobley, A. E., Nugent, T. C., Bryson, K., and Jones, D. T. (2010). Protein annotation and modelling servers at University College London. Nucleic Acids Res. 38, W563–W568. doi: 10.1093/nar/gkq427
Carn, F., Guyot, S., Baron, A., Perez, J., Buhler, E., and Zanchi, D. (2012). Structural properties of colloidal complexes between condensed tannins and polysaccharide hyaluronan. Biomacromolecules 13, 751–759. doi: 10.1021/bm201674n
DéAndrea, L. D., and Regan, L. (2003). TPR proteins: the versatile helix. Trends Biochem. Sci. 28, 655–662. doi: 10.1016/j.tibs.2003.10.007
DeLano, W. L. (2002). Pymol: an open-source molecular graphics tool. CCP4 Newslett. Protein Crystallogr. 40. Available online at: http://www.ccp4.ac.uk/newsletters/newsletter40/11_pymol.html
Deng, H. X., Hentati, A., Tainer, J. A., Iqbal, Z., Cayabyab, A., Hung, W. Y., et al. (1993). Amyotrophic lateral sclerosis and structural defects in Cu,Zn superoxide dismutase. Science 261, 1047–1051. doi: 10.1126/science.8351519
Dietmann, S., Park, J., Notredame, C., Heger, A., Lappe, M., and Holm, L. (2001). A fully automatic evolutionary classification of protein folds: Dali Domain Dictionary version 3. Nucleic Acids Res. 29, 55–57. doi: 10.1093/nar/29.1.55
Dong, A., Xu, X., Edwards, A. M., Chang, C., Chruszcz, M., Cuff, M., et al. (2007). In situ proteolysis for protein crystallization and structure determination. Nat. Methods 4, 1019–1021. doi: 10.1038/nmeth1118
Emsley, P., and Cowtan, K. (2004). Coot: model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 60, 2126–2132. doi: 10.1107/S0907444904019158
Emsley, P., Lohkamp, B., Scott, W. G., and Cowtan, K. (2010). Features and development of Coot. Acta Crystallogr. D Biol. Crystallogr. 66, 486–501. doi: 10.1107/S0907444910007493
Faivre, D., and Schüler, D. (2008). Magnetotactic bacteria and magnetosomes. Chem. Rev. 108, 4875–4898. doi: 10.1021/cr078258w
Frankel, R. B., and Bazylinski, D. A. (2006). How magnetotactic bacteria make magnetosomes queue up. Trends Microbiol. 14, 329–331. doi: 10.1016/j.tim.2006.06.004
Frankel, R. B., Bazylinski, D. A., Johnson, M. S., and Taylor, B. L. (1997). Magneto-aerotaxis in marine coccoid bacteria. Biophys. J. 73, 994–1000. doi: 10.1016/S0006-3495(97)78132-3
Gao, M., Craig, D., Lequin, O., Campbell, I. D., Vogel, V., and Schulten, K. (2003). Structure and functional significance of mechanically unfolded fibronectin type III1 intermediates. Proc. Natl. Acad. Sci. U.S.A. 100, 14784–14789. doi: 10.1073/pnas.2334390100
Gouet, P., Robert, X., and Courcelle, E. (2003). ESPript/ENDscript: extracting and rendering sequence and 3D information from atomic structures of proteins. Nucleic Acids Res. 31, 3320–3323. doi: 10.1093/nar/gkg556
Guerrero, S. A., Hecht, H. J., Hofmann, B., Biebl, H., and Singh, M. (2001). Production of selenomethionine-labelled proteins using simplified culture conditions and generally applicable host/vector systems. Appl. Microbiol. Biotechnol. 56, 718–723. doi: 10.1007/s002530100690
Han, D., Oh, J., Kim, K., Lim, H., and Kim, Y. (2007). Crystal structure of YrrB: a TPR protein with an unusual peptide-binding site. Biochem. Biophys. Res. Commun. 360, 784–790. doi: 10.1016/j.bbrc.2007.06.129
Hart, P. J., Liu, H., Pellegrini, M., Nersissian, A. M., Gralla, E. B., Valentine, J. S., et al. (1998). Subunit asymmetry in the three-dimensional structure of a human CuZnSOD mutant found in familial amyotrophic lateral sclerosis. Protein Sci. 7, 545–555.
Hemmingsen, J. M., Gernert, K. M., Richardson, J. S., and Richardson, D. C. (1994). The tyrosine corner: a feature of most Greek key beta-barrel proteins. Protein Sci. 3, 1927–1937.
Hoxha, E., and Campion, S. R. (2014). Structure-critical distribution of aromatic residues in the fibronectin type III protein family. Protein J. 33, 165–173. doi: 10.1007/s10930-014-9549-y
Kajander, T., Cortajarena, A. L., Mochrie, S., and Regan, L. (2007). Structure and stability of designed TPR protein superhelices: unusual crystal packing and implications for natural TPR proteins. Acta Crystallogr D Biol. Crystallogr. 63, 800–811. doi: 10.1107/S0907444907024353
Kolinko, I., Lohsse, A., Borg, S., Raschdorf, O., Jogler, C., Tu, Q., et al. (2014). Biosynthesis of magnetic nanostructures in a foreign organism by transfer of bacterial magnetosome gene clusters. Nat. Nanotechnol. 9, 193–197. doi: 10.1038/nnano.2014.13
Komeili, A. (2007). Molecular mechanisms of magnetosome formation. Annu. Rev. Biochem. 76, 351–366. doi: 10.1146/annurev.biochem.74.082803.133444
Kryshtafovych, A., Monastyrskyy, B., and Fidelis, K. (2014). CASP prediction center infrastructure and evaluation measures in CASP10 and CASP ROLL. Proteins 82(Suppl. 2), 7–13. doi: 10.1002/prot.24399
Langer, G., Cohen, S. X., Lamzin, V. S., and Perrakis, A. (2008). Automated macromolecular model building for X-ray crystallography using ARP/wARP version 7. Nat. Protoc. 3, 1171–1179. doi: 10.1038/nprot.2008.91
Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi: 10.1093/bioinformatics/btm404
Lefevre, C. T., Bennet, M., Landau, L., Vach, P., Pignol, D., Bazylinski, D. A., et al. (2014). Diversity of magneto-aerotactic behaviors and oxygen sensing mechanisms in cultured magnetotactic bacteria. Biophys. J. 107, 527–538. doi: 10.1016/j.bpj.2014.05.043
Louis-Jeune, C., Andrade-Navarro, M. A., and Perez-Iratxeta, C. (2012). Prediction of protein secondary structure from circular dichroism using theoretically derived spectra. Proteins 80, 374–381. doi: 10.1002/prot.23188
Lundell, A., Olin, A. I., Morgelin, M., al-Karadaghi, S., Aspberg, A., and Logan, D. T. (2004). Structural basis for interactions between tenascins and lectican C-type lectin domains: evidence for a crosslinking role for tenascins. Structure 12, 1495–1506. doi: 10.1016/j.str.2004.05.021
Matsunaga, T., Okamura, Y., Fukuda, Y., Wahyudi, A. T., Murase, Y., and Takeyama, H. (2005). Complete genome sequence of the facultative anaerobic magnetotactic bacterium Magnetospirillum sp. strain AMB-1. DNA Res. 12, 157–166. doi: 10.1093/dnares/dsi002
Matthews, B. W. (1968). Solvent content of protein crystals. J. Mol. Biol. 33, 491–497. doi: 10.1016/0022-2836(68)90205-2
Monaco, S., Gordon, E., Bowler, M. W., Delageniere, S., Guijarro, M., Spruce, D., et al. (2013). Automatic processing of macromolecular crystallography X-ray diffraction data at the ESRF. J. Appl. Crystallogr. 46, 804–810. doi: 10.1107/S0021889813006195
Niimi, T., Osawa, M., Yamaji, N., Yasunaga, K., Sakashita, H., Mase, T., et al. (2001). NMR structure of human fibronectin EDA. J. Biomol. NMR 21, 281–284. doi: 10.1023/A:1012947209393
Nudelman, H., and Zarivach, R. (2014). Structure prediction of magnetosome-associated proteins. Front. Microbiol. 5:9. doi: 10.3389/fmicb.2014.00009
Otwinowski, Z., and Minor, W. (1997). “Processing of X-ray diffraction data collected in oscillation mode,” in Methods in Enzymology, Vol. 276, eds C. W. Jr. Carter and R. M. Sweet (New York, NY: Academic Press), 307–326.
Popp, F., Armitage, J. P., and Schüler, D. (2014). Polarity of bacterial magnetotaxis is controlled by aerotaxis through a common sensory pathway. Nat. Commun. 5, 5398. doi: 10.1038/ncomms6398
Potts, J. R., and Campbell, I. D. (1994). Fibronectin structure and assembly. Curr. Opin. Cell Biol. 6, 648–655. doi: 10.1016/0955-0674(94)90090-6
Raschdorf, O., Plitzko, J. M., Schüler, D., and Müller, F. D. (2014). A tailored galK counterselection system for efficient markerless gene deletion and chromosomal tagging in Magnetospirillum gryphiswaldense. Appl. Environ. Microbiol. 80, 4323–4330. doi: 10.1128/AEM.00588-14
Remaut, H., and Waksman, G. (2006). Protein-protein interaction through beta-strand addition. Trends Biochem. Sci. 31, 436–444. doi: 10.1016/j.tibs.2006.06.007
Richter, M., Kube, M., Bazylinski, D. A., Lombardot, T., Glockner, F. O., Reinhardt, R., et al. (2007). Comparative genome analysis of four magnetotactic bacteria reveals a complex set of group-specific genes implicated in magnetosome biomineralization and function. J. Bacteriol. 189, 4899–4910. doi: 10.1128/JB.00119-07
Schübbe, S., Kube, M., Scheffel, A., Wawer, C., Heyen, U., Meyerdierks, A., et al. (2003). Characterization of a spontaneous nonmagnetic mutant of Magnetospirillum gryphiswaldense reveals a large deletion comprising a putative magnetosome island. J. Bacteriol. 185, 5779–5790. doi: 10.1128/JB.185.19.5779-5790.2003
Schüler, D. (2008). Genetics and cell biology of magnetosome formation in magnetotactic bacteria. FEMS Microbiol. Rev. 32, 654–672. doi: 10.1111/j.1574-6976.2008.00116.x
Schüler, D., Uhl, R., and Bäuerlein, E. (1995). A simple light scattering method to assay magnetism in Magnetospirillum gryphiswaldense. FEMS Microbiol. Lett. 132, 139–145. doi: 10.1111/j.1574-6968.1995.tb07823.x
Schultheiss, D., and Schüler, D. (2003). Development of a genetic system for Magnetospirillum gryphiswaldense. Arch. Microbiol. 179, 89–94. doi: 10.1007/s00203-002-0498-z
Svergun, D. I. (1992). Determination of the regularization parameter in indirect-transform methods using perceptual criteria. J. Appl. Crystallogr. 25, 495–503. doi: 10.1107/S0021889892001663
Svergun, D. I. (1999). Restoring low resolution structure of biological macromolecules from solution scattering using simulated annealing. Biophys. J. 77, 2896–2896.
Tanaka, M., Okamura, Y., Arakaki, A., Tanaka, T., Takeyama, H., and Matsunaga, T. (2006). Origin of magnetosome membrane: proteomic analysis of magnetosome membrane and comparison with cytoplasmic membrane. Proteomics 6, 5234–5247. doi: 10.1002/pmic.200500887
Ullrich, S., and Schüler, D. (2010). Cre-lox-based method for generation of large deletions within the genomic magnetosome island of Magnetospirillum gryphiswaldense. Appl. Environ. Microbiol. 76, 2439–2444. doi: 10.1128/AEM.02805-09
Vagin, A. A., Steiner, R. A., Lebedev, A. A., Potterton, L., McNicholas, S., Long, F., et al. (2004). REFMAC5 dictionary: organization of prior chemical knowledge and guidelines for its use. Acta Crystallogr. D Biol. Crystallogr. 60, 2184–2195. doi: 10.1107/S0907444904023510
Volkov, V. V., and Svergun, D. I. (2003). Uniqueness of ab initio shape determination in small-angle scattering. J. Appl. Crystallogr. 36, 860–864. doi: 10.1107/S0021889803000268
Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., et al. (2011). Overview of the CCP4 suite and current developments. Acta Crystallogr. D Biol. Crystallogr. 67, 235–242. doi: 10.1107/S0907444910045749
Zeytuni, N., Offer, T., Davidov, G., and Zarivach, R. (2012). Crystallization and preliminary crystallographic analysis of the C-terminal domain of MamM, a magnetosome-associated protein from Magnetospirillum gryphiswaldense MSR-1. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 68, 927–930. doi: 10.1107/S1744309112025638
Zeytuni, N., Ozyamak, E., Ben-Harush, K., Davidov, G., Levin, M., Gat, Y., et al. (2011). Self-recognition mechanism of MamA, a magnetosome-associated TPR-containing protein, promotes complex assembly. Proc. Natl. Acad. Sci. U.S.A. 108, E480–E487. doi: 10.1073/pnas.1103367108
Zeytuni, N., and Zarivach, R. (2012). Structural and functional discussion of the tetra-trico-peptide repeat, a protein interaction module. Structure 20, 397–405. doi: 10.1016/j.str.2012.01.006
Zhu, X., Ge, X., Li, N., Wu, L. F., Luo, C., Ouyang, Q., et al. (2014). Angle sensing in magnetotaxis of Magnetospirillum magneticum AMB-1. Integr. Biol. 6, 706–713. doi: 10.1039/c3ib40259b
Keywords: MtxA, magnetotactic bacteria, magnetosome genomic island, magnetotaxis, magnetosome, tetratricopeptide repeats, immunoglobulin-like domain
Citation: Davidov G, Müller FD, Baumgartner J, Bitton R, Faivre D, Schüler D and Zarivach R (2015) Crystal structure of the magnetobacterial protein MtxA C-terminal domain reveals a new sequence-structure relationship. Front. Mol. Biosci. 2:25. doi: 10.3389/fmolb.2015.00025
Received: 18 November 2014; Accepted: 04 May 2015;
Published: 21 May 2015.
Edited by:
Sandra Macedo-Ribeiro, Instituto de Biologia Molecular e Celular, PortugalReviewed by:
Margarida Archer Frazao, Instituto de Tecnologia Química e Biológica (ITQB-Universidade Nova de Lisboa), PortugalMenico Rizzi, University of Piemonte Orientale Amedeo Avogadro, Italy
David Pignol, Commissariat à l'énergie Atomique (Department of Energy), France
Copyright © 2015 Davidov, Müller, Baumgartner, Bitton, Faivre, Schüler and Zarivach. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Raz Zarivach, Department of Life Sciences, Ben Gurion University of the Negev, PO Box. 653, Beer-Sheva 84105, Israel, zarivach@bgu.ac.il