
Conformational Plasticity of the Essential Membrane-associated Mannosyltransferase PimA from Mycobacteria*
David Giganti
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
¶Department of Biological Sciences, Columbia University, New York, New York 10027,
Jorge Alegre-Cebollada
¶Department of Biological Sciences, Columbia University, New York, New York 10027,
Saioa Urresti
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
David Albesa-Jové
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
Ane Rodrigo-Unzueta
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
Natalia Comino
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
Michael Kachala
‖European Molecular Biology Laboratory, Hamburg Outstation, c/o Deutsches Elektronen-Synchrotron (DESY), Notkestrasse 85, D-22603 Hamburg, Germany, and
Sonia López-Fernández
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
Dmitri I. Svergun
‖European Molecular Biology Laboratory, Hamburg Outstation, c/o Deutsches Elektronen-Synchrotron (DESY), Notkestrasse 85, D-22603 Hamburg, Germany, and
Julio M. Fernández
¶Department of Biological Sciences, Columbia University, New York, New York 10027,
Marcelo E. Guerin
From the ‡Unidad de Biofísica, Consejo Superior de Investigaciones Científicas-Universidad del País Vasco/Euskal Herriko Unibertsitatea, Barrio Sarriena s/n, Leioa, Bizkaia 48940, Spain,
§Departamento de Bioquímica, Universidad del País Vasco, Leioa, Bizkaia 48940, Spain,
**IKERBASQUE, Basque Foundation for Science, Bilbao 48011, Spain
Abstract
Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential glycosyltransferase (GT) that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides, lipomannan, and lipoarabinomannan, which are key glycolipids/lipoglycans of the mycobacterial cell envelope. PimA belongs to a large family of peripheral membrane-associated GTs for which the understanding of the molecular mechanism and conformational changes that govern substrate/membrane recognition and catalysis remains a major challenge. Here we used single molecule force spectroscopy techniques to study the mechanical and conformational properties of PimA. In our studies, we engineered a polyprotein containing PimA flanked by four copies of the well characterized I27 protein, which provides an unambiguous mechanical fingerprint. We found that PimA exhibits weak mechanical stability albeit displaying β-sheet topology expected to unfold at much higher forces. Notably, PimA unfolds following heterogeneous multiple step mechanical unfolding pathways at low force akin to molten globule states. Interestingly, the ab initio low resolution envelopes obtained from small angle x-ray scattering of the unliganded PimA and the PimA·GDP complexed forms clearly demonstrate that not only the “open” and “closed” conformations of the GT-B enzyme are largely present in solution, but in addition, PimA experiences remarkable flexibility that undoubtedly corresponds to the N-terminal “Rossmann fold” domain, which has been proved to participate in protein-membrane interactions. Based on these results and on our previous experimental data, we propose a model wherein the conformational transitions are important for the mannosyltransferase to interact with the donor and acceptor substrates/membrane.
Introduction
Glycosyltransferases (GTs)6 play a central role in nature. GTs catalyze the stereo- and regiospecific transfer of a sugar moiety from nucleotide-sugar or lipid-phospho-sugar donors to a wide range of acceptor substrates, including mono-, oligo-, and polysaccharides; lipids; proteins; small organic molecules; and deoxyribonucleic acids (1). As a consequence, GTs generate a significant amount of structural diversity in biological systems. This structural information is particularly apparent not only in the maintenance of the structural integrity of the cell but also in the modulation of molecular recognition events, including cell-signaling, cell-cell, and cell-pathogen interactions (2). Glycosyl transfer can proceed with either “inversion” or “retention” of the anomeric configuration with respect to the sugar donor. “Inverting” GTs seem to follow a single displacement mechanism with an oxocarbenium ion-like transition state and an asynchronous SN2 mechanism analogous to that observed for inverting glycosyl hydrolases (1). In contrast, the catalytic mechanism for “retaining” enzymes is actually a matter of strong debate. Retaining GTs are believed to follow the “internal return-like” mechanism in which leaving group departure and nucleophilic attack occur on the same face of the sugar (3) involving either a short lived oxocarbenium ion intermediate (SNi-like; Ref. 4) or an oxocarbenium ion transition state (SNi; Ref. 5). However, recent reports suggest the presence of a covalent intermediate in the α3-glycosyltransferase and the human blood group-synthesizing α-(1→3)-N-acetylgalactosaminyltransferase and α-(1→3)-galactosyltransferase (6, 7). Whether or not retaining GTs could proceed via different catalytic mechanisms is a notion that will certainly need further experimental support. Importantly, to achieve such enzyme·transition state complexes, a spatial rearrangement of the active site or more distant regions of the GTs would often be required. However, to date, the impact of substrate and enzyme dynamics and conformational changes in GT-mediated catalysis is largely unknown (1).
Phosphatidyl-myo-inositol mannosyltransferase A (PimA) is a retaining membrane-associated GT that initiates the biosynthetic pathway of essential phosphatidyl-myo-inositol mannosides, lipomannan, and lipoarabinomannan in mycobacteria (8–10). Phosphatidyl-myo-inositol mannoside, lipomannan, and lipoarabinomannan are considered not only essential structural components of the cell but also important molecules implicated in host-pathogen interactions (10–13). PimA catalyzes the transfer of a mannose residue from guanosine 5′-diphosphate mannose (GDP-Manp) to the 2-position of phosphatidyl-myo-inositol (PI) to form PI monomannoside on the cytoplasmic side of the plasma membrane. The crystal structure of PimA from Mycobacterium smegmatis has been recently solved in complex with both the donor substrate GDP-Manp and GDP, one of the reaction products. Both structures superimpose well (root mean square deviation of 0.3 Å for 361 identical residues) and reveal that the enzyme displays the typical GT-B fold of GTs, one of the two major structural folds described for the nucleotide-sugar-dependent enzymes (Ref. 14 and see Fig. 1A). This fold was first described for the 351-amino acid DNA-modifying β-glucosyltransferase from family GT63, an inverting GT from bacteriophage T4, and was found to be structurally related to the catalytic core of glycogen phosphorylase (15–18). The GT-B fold displays two “Rossmann fold” domains separated by a deep cleft that includes the catalytic center. Therefore, an important interdomain movement has been predicted/demonstrated in some members of this superfamily, including MurG (19), glycogen synthase (20–22), and MshA (23), during substrate binding and catalysis. It is generally accepted that in GT-B enzymes the nucleotide-sugar donors mainly bind to the C-terminal domain of the protein, whereas the N-terminal domain is involved in acceptor substrate recognition. As acceptors exhibit a marked diversity of chemical structures compared with nucleotide-sugar donors, the N-terminal domains reflect this variability by showing different rearrangements of secondary structural elements (24). On the basis of primary sequence homology analysis, a glycogen phosphorylase/glycosyltransferase family motif has been suggested to be present in many GT-B enzymes (18, 25). However, GT-B enzymes do not seem to share any strictly conserved residue (26).
Schematic representation of I272-PimA-I272 polyprotein. A, the crystal structure of PimA in complex with GDP displays the typical GT-B fold of GTs. The N-terminal domain contains the residues 1–168 (orange) and the helical C-terminal extremity 348–373 (yellow), whereas the C-terminal domain includes residues 169–310 (red) and the segment 311–347 (yellow). The last two C-terminal helices (yellow) participate in the hinge. B, configuration of the atomic force microscope experiments. The polyprotein is tethered between the tip of a cantilever and a gold slide on a piezo positioner. PimA (red and orange) is flanked by four immunoglobulin I27 domains (blue).
Structural, biochemical, and biophysical evidence indicates that PimA undergoes significant conformational changes upon substrate binding and catalysis involving loop flexibility and domain motions (14, 27). In particular, the binding of the sugar donor substrate GDP-Manp triggers an important interdomain rearrangement from an “open” to a “closed” state that stabilizes the enzyme and markedly increases its affinity for the acceptor substrate, PI. According to the crystal structures of the PimA·GDP and PimA·GDP-Manp complexes, the interaction of the enzyme with the β-phosphate of GDP is essential for this conformational change to occur. The open-to-closed motion brings together critical residues from the N- and C-terminal domains, allowing the formation of a functionally competent active site. In contrast, the binding of PI to the enzyme had the opposite effect, inducing the formation of a more relaxed complex with PimA (27). Membrane association seems to be mediated by a combination of positively charged residues on the amphipathic α-helix 2 and surface-exposed hydrophobic residues on the N-terminal domain of the protein (14, 27).
This study describes a detailed investigation of the mechanical properties of PimA by using single molecule force spectroscopy techniques. In combination with the solution structure of PimA by small-angle X-ray scattering, we propose a model for substrate binding wherein protein flexibility and conformational transitions play a prominent role.
EXPERIMENTAL PROCEDURES
Methods
Recombinant PimA from M. smegmatis (MsPimA; MSMEG_2935) was produced in Escherichia coli and purified to apparent homogeneity as described previously (27).
I272-PimA-I272 Polyprotein Production
The polyprotein I272-PimA-I272 was generated by using recombinant DNA technology based on an iterative process of digestion and ligation of DNA fragments as described previously (28). The pimA gene from M. smegmatis was amplified by standard PCR using oligonucleotide primers pimA_BamHI_Fwd (5′-CGCGGATCCGGCTCTGGCGCGATGCGTATCGGGATGGTCTGCC-3′) and pimA_SacI_Rev (5′-ATAATAGAGCTCTTAGCAACAAGATCTGACCGATTCTCCGGCCGTCTG-3′), Phusion DNA polymerase (Finnzymes, Thermo Fisher Scientific), and the plasmid pET29a-MspimA as the DNA template (27). The PCR fragment was digested with BamHI and SacI and ligated to the vector pT7-I272 harboring two consecutive copies of the immunoglobulin (Ig) domain of the giant muscle protein titin, I27. The resulting plasmid, pT7-I272-PimA, was further digested with BglII and SacI, and the fragment LF-pT7-I272-PimA was purified. The pT7-I272 was digested with BamHI and SacI, and the fragment SF-I272 was purified. Both fragments LF-pT7-I272-PimA and SF-I272 were ligated, and the resulting construct was named pT7-I272-PimA-I272. Digestion of pT7-I272-PimA-I272 with BamHI and SacI allowed in-frame cloning of the I272-PimA-I272 into the expression vector pQE80L (Qiagen, Chatsworth, CA). As depicted in supplemental Fig. 1S, the full-length construct resulted in the following amino acid additions. (i) The N-terminal sequence is Met-Arg-Gly-Ser-His6-Gly-Ser-I27. (ii) The junction between the I27 domains (BamHI-BglII hybrid site) is Arg-Ser. (iii) The junction between the second I27 domain and PimA is Arg-Ser-Gly-Ser-Gly-Ala. (iv) The junction between PimA and the third I27 domain is Arg-Ser. In addition, the resulting construct contained two Cys residues at the C terminus, allowing the covalent attachment of the polyprotein to the gold-coated coverslips. E. coli BL3 cells transformed with pQE80L-I272-PimA-I272 were grown in 400 ml of LB medium supplemented with 30 μg/ml kanamycin and 50 μg/ml carbenicillin at 37 °C. When the culture reached an A600 value of 0.6, I272-PimA-I272 expression was induced by adding 1 mm isopropyl β-d-1-thiogalactopyranoside (MIP). After overnight incubation at 25 °C, cells were harvested and resuspended in 16 ml of 50 mm Tris-HCl, pH 8.0 (solution A) containing protease inhibitors (Complete EDTA-free, Roche Applied Science). Cells were then disrupted by sonication (five cycles of 1 min each), and the suspension was centrifuged for 20 min at 10,000 × g. The supernatant was applied to a Talon resin column (Clontech) followed by a gel filtration step in a Superdex 200 column (GE Healthcare) equilibrated in 10 mm Hepes, pH 7.2, 1 mm EDTA, 150 mm NaCl (solution B). The resulting I272-PimA-I272 preparation displayed a single protein band when run on a 12% NuPAGE® Bis-Tris precast gel stained with SimplyBlueTM SafeStain (Invitrogen). The purified I272-PimA-I272 polyprotein was stored at −80 °C for further use in single molecule experiments.
I27-PimA Protein Production
The pimA gene from M. smegmatis was amplified as described for the I272-PimA-I272 construct. The PCR fragment was digested with BamHI and SacI and then ligated to the vector pT7-I27 harboring one copy of the I27 module. The resulting plasmid, pT7-I27-PimA, was further digested with BamHI and SacI and ligated to the expression vector pQE80L (Qiagen). As depicted in supplemental Fig. 2S, the full-length construct resulted in the following amino acid additions. (i) The N-terminal sequence is Met-Arg-Gly-Ser-His6-Gly-Ser-I27. (ii) The junction between the I27 domain and PimA is Arg-Ser-Gly-Ser-Gly-Ala. In addition, the final construct contained two Cys residues at the C terminus. E. coli DH5α cells (Invitrogen) transformed with pQE80L-I27-PimA were grown in 3000 ml of 2×YT medium supplemented with 50 μg/ml carbenicillin at 37 °C. When the culture reached an A600 value of 0.6, I27-PimA expression was induced by adding 1 mm isopropyl β-d-1-thiogalactopyranoside (MIP). After overnight incubation at 18 °C, cells were harvested and resuspended in 40 ml of solution A containing protease inhibitors (Complete EDTA-free, Roche Applied Science). Cells were then disrupted by sonication (five cycles of 1 min each), and the suspension was centrifuged for 20 min at 10,000 × g. The supernatant was applied onto a HisTrap column (GE Healthcare) followed by a gel filtration step in a Superdex 75 column (GE Healthcare) equilibrated in buffer 50 mm Tris-HCl, pH 7.5, 150 mm NaCl (solution C). The resulting I27-PimA preparation displayed a single protein band when run on a 12% NuPAGE Bis-Tris precast gel stained with SimplyBlue SafeStain (Invitrogen). The purified I27-PimA fusion protein was stored at 4 °C and then concentrated for small angle x-ray scattering (SAXS) measurements by using a Vivaspin 20 spin concentrator (Vivascience) with a 10-kDa-molecular mass cutoff. The homogeneity of I27-PimA was confirmed by dynamic light scattering (data not shown).
Single Molecule Force Spectroscopy Measurements
5–15 μl of I272-PimA-I272 from the most concentrated fraction obtained after gel filtration were deposited onto an evaporated gold coverslip. Veeco silicon nitride MLCT (Bruker, Camarillo, CA) cantilevers were mounted in a fluid cell carrier. The laser beam was centered at the tip of the cantilever and carefully calibrated. Spring constants around 15 piconewtons/nm were obtained using the equipartition theorem (29). The cantilever was kept in contact with the gold coverslip for 1–2 s before retraction. We used a custom-built atomic force microscopy setup that allows a strict control of the force based on an electronic feedback. Experiments were performed at room temperature in 10 mm Hepes, pH 7.2, 1 mm EDTA, 150 mm NaCl. Only traces displaying four steps corresponding to I27 unfolding were considered. Because of the low stability of the construct I272-PimA-I272, single molecule experiments were performed within 3–4 days after purification. During that time, the purified I272-PimA-I272 polyprotein construct was demonstrated to conserve full enzymatic activity in vitro (Ref. 8 and data not shown). To calculate the contour length of the pulled polypeptide, we applied the wormlike chain model with a persistent length of 0.4 nm. Hence, the wormlike chain curves in Fig. 3 correspond to the theoretical extension of the N- and C-terminal Rossmann fold domains, 163 and 119 amino acids, respectively, and I27. For a more accurate estimation of the contour length, the length of the folded N- and C-terminal Rossmann fold domains (1.8 and 2.3 nm, respectively) was subtracted. The wormlike chain was fitted in force-extension curves as well to estimate the initial extension before the rupture of the first I27 domain (30).
Extension of I272-PimA-I272 during a linear increase in force. A and B, typical force-ramp traces of I272-PimA-I272 in the absence and presence of GDP, respectively. Step size events of PimA are annotated in red. C and D, distribution of step sizes and unfolding forces of I27 and PimA unfolding events from 50 unambiguous traces. Theoretical extensions following the wormlike chain of I27 and the N-terminal and C-terminal β-sheets of PimA are plotted in blue, red dashed, and red dotted lines, respectively. pN, piconewtons.
Small Angle X-ray Scattering Measurements
Synchrotron x-ray diffraction data for recombinant purified MsPimA apo and MsPimA in the presence of GDP (27) were collected on a pixel Pilatus 1M detector at the European Molecular Biology Laboratory X33 beamline (Deutsches Elektronen-Synchrotron, Hamburg, Germany; Refs. 31 and 32), whereas data for the fusion I27-PimA protein were collected at the European Molecular Biology Laboratory P12 beamline (Deutsches Elektronen-Synchrotron) on a pixel Pilatus 2M detector. The scattering patterns were measured with a 2-min exposure time for protein samples at a minimum of three different protein concentrations ranging from 1 to 10 mg/ml at the X33 and with a 1-s exposure time at the P12. Samples were incubated with 2 mm dithiothreitol before data collection. To check for radiation damage, four 30-s exposures were compared at the X33, and 20 50-millisecond exposures were compared at the P12; no radiation damage was observed. Using the sample-to-detector distance of 2.7 m (3.1 m at the P12), the range of momentum transfer values is 0.009 < q < 0.5 Å−1 at the X33 and 0.008 < q < 0.47 Å−1 at the P12 (q = 4π sin(θ)/λ where 2θ is the scattering angle and λ = 1.5 Å (1.24 Å for the P12 beamline) is the x-ray wavelength). Data were processed using standard procedures and extrapolated to infinite dilution by the program package PRIMUS (33). The forward scattering (I(0)) was evaluated using the Guinier approximation (34) assuming the intensity is represented as I(q) = I(0)exp(−(qRg)2/3) for a very small range of momentum transfer values (q < 1.3/Rg). The maximum dimensions (Dmax), the interatomic distance distribution functions (P(r)), and the radii of gyration (Rg) were computed using GNOM (35). The molecular mass of the protein was evaluated by comparison of the forward scattering with that from a reference solution of bovine serum albumin (66 kDa).
Ab Initio Shape Determination
The low resolution structures of MsPimA apo and MsPimA·GDP complex were calculated ab initio by using the program DAMMIF (36). Structure clustering and averaging were carried out with the program DAMCLUST (37). The results and statistics are summarized in Table 1. For each data set, pairwise alignment of structures was carried out, and the normalized spatial discrepancy (NSD) was calculated (38). The clustering process groups structures with lower NSD. In general, NSD values close to 1 indicate that the two structures are similar. Structures that deviate by more than 2 are taken as outliers and removed. For all calculations, three outlier models were rejected. Average NSD values within clusters are between 0.5–7, indicating a very good structural agreement for each cluster. Furthermore, average NSD values for structures of different clusters are between 0.7 and 1.0, indicating similarity between clusters.
TABLE 1
Parameters of the unliganded and GDP-bound forms of PimA and the I27-PimA fusion protein calculated from the SAXS data

a Values in parentheses are estimated errors approximated to the last decimal place.
b q-range used for calculation of P(r) function.
c Rg value extrapolated to 0 mg/ml concentration.
d Molecular mass estimated by calculation of I(0) and comparison against bovine serum albumin.
e Values in parentheses are standard deviations calculated for 20 ab initio calculations.
f MONSA R-factors are reported for each experimental curve employed in the calculation (i.e. SAXS data of I27-PimA and PimA apo). NSD values are reported for the individual phases reconstructed (i.e. I27 and PimA apo).
The ab initio reconstruction of the fusion protein I27-PimA was carried out by MONSA. MONSA is based on DAMMIN and allows bead modeling of macromolecular complexes by simultaneously fitting multiple experimental data (39, 40). The search volume was defined as a sphere of radius equal to half the maximum dimension of I27-PimA (Dmax = 120 Å) as computed by GNOM. Within the search volume, a grid of points is defined with a constant distance of 2 Å between points, resulting in dense packing of spheres. Each sphere within the search volume might correspond to one of the following components: solvent, I27, or PimA, represented by an integer number, 0, 1, or 2. To facilitate computation and to keep physical sense with respect to the experimental data (i.e. low resolution of data and connectivity of each component), a number of conditions are applied during structural calculations accounting for the likely connectivity of two arbitrarily selected atoms belonging to a given component and its compactness. The minimization procedure fits the two experimental data sets, PimA and I27-PimA, to the corresponding calculated data obtained from computed structures. It is carried out by simulated annealing (41) starting from a random configuration X0 where the assignment of a single bead to a given component is changed randomly (a move from X to X′). In the case that the change produces an improvement of fitness (Δ = f(X′) − f(X) < 0), it is accepted. Otherwise, it can still be accepted with a probability of exp(−Δ/T) where T is the annealing temperature, which decreases as the calculation progresses toward a minimum. The described approach allowed us to unambiguously differentiate the relative location of I27 and PimA within the fusion protein and to correctly orientate the N- and C-terminal domains of PimA as observed in the crystal structure in the solution structures obtained by SAXS techniques.
Molecular Modeling
SITUS (42) was used to perform rigid body docking of high resolution crystal structures into solution structures. To this end, the bead models obtained by SAXS techniques were transformed to three-dimensional volumes using the pdb2vol program (part of the SITUS package) followed by exhaustive search docking, including all six degrees of freedom (three rotation angles and three translation coordinates). The fitting improvement is assessed by a cross-correlation coefficient function (CC). Cross-correlation coefficient values range between 0 and 1 where 1 indicates excellent agreement between the three-dimensional volume (obtained by SAXS) and the crystal structures (42).
RESULTS AND DISCUSSION
PimA Displays an Intrinsic Low Mechanical Stability
Single molecule force spectroscopy has been proven to be a powerful technique not only to understand the role of force in mechanical stability but also to study protein conformational dynamics and enzyme catalysis (43, 44) For instance, applying a calibrated force to the substrate of an enzyme can probe the dynamics of the active site during catalysis with subangstrom resolution (45). To study the mechanical properties of PimA under force, we engineered a polyprotein composed of a single PimA copy flanked on either side by two copies of an Ig domain of the giant muscle protein titin, I27 (I272-PimA-I272) (see “Experimental Procedures” and Ref. 46).
We applied force to this polyprotein with an atomic force microscope as shown in Fig. 1B. I27 modules provide an unambiguous mechanical fingerprint to identify single molecule events. Two cysteines at the C terminus allow a covalent binding of I272-PimA-I272 to a gold surface on top of a piezoelectric positioner (Fig. 1B and Ref. 28). To calculate the expected initial extension, we considered the sequences linking the I27 domains and PimA in the polyprotein. The N terminus before the first I27 domain contains 11 amino acids, whereas the C terminus after the last I27 contains four amino acids. Twelve extra residues separate the I27 modules in the polyprotein (supplemental Fig. 1S). Therefore, the initial extension of I272-PimA-I272, including the unfolding of the 386 amino acids of PimA, involves 386 + 12 + 12 + 4 = 414 amino acids and four folded I27 domains (4.4 nm each). Considering a contour length per amino acid of 0.4 nm (47), this corresponds to a maximal contour length of 183 nm. The polyprotein was first stretched at a constant rate of 400 nm/s. The opposing force exerted by the polyprotein was measured from the deflection of the atomic force microscope cantilever. The resulting force extension traces displayed an initial extension region characterized by the presence of marginal unfolding peaks. The initial extension was followed by a sawtooth pattern of three or four consecutive unfolding events and a final detachment peak (Fig. 2). The observed unfolding peaks of ∼200 piconewtons spaced by ∼28–29 nm (number of events = 138) fit well to the wormlike chain model of polymer elasticity and correspond to the characteristic fingerprint of the I27 module (48). Based on the contour length of the first I27 peak, we determined the length of the initial extension as ∼180 nm (supplemental Fig. 1S and Fig. 2), which is in agreement with a fully extended PimA and sequences linking the I27 domains (182 nm; see above). The addition of 5 mm GDP gives similar curves without any distinct peak during the extension of PimA (data not shown). Altogether, these results demonstrate that I272-PimA-I272 is tethered from both ends, and PimA is subject to force. It is worth noting that the addition of protein domains with known mechanical features has been successfully used to study new protein folds under force (46, 49). In the case of proteins without a uniform and discernible mechanical response to force, this indispensable fingerprint guaranties that the initial extension is relevant to the pulled polyprotein (50). PimA seems to unfold at low force without significant energy barriers limiting its extensibility, indicating that the mannosyltransferase displays an intrinsic weak mechanical resistance to force.

Mechanical properties of I272-PimA-I272 as visualized by the force-extension mode of atomic force microscope. A, typical force-extension profile of I272-PimA-I272 essentially showing the length of the initial extension of PimA as ∼180 nm. The GDP-bound form of PimA displayed a similar behavior (data not shown). B, the histogram includes experimental results from both PimA apo and PimA-GDP conditions. pN, piconewtons.
PimA Unfolds following Heterogeneous Multiple Step Mechanical Unfolding Pathways
To further characterize PimA unfolding, we used a force ramp protocol with a linearly increasing force (51). This protocol explores a wide range of forces and is therefore more suitable to examine and segregate different unfolding pathways. We increased the force from 0 to 250 piconewtons over a time period of 5 s, obtaining 50 well defined traces in the presence and absence of GDP. Extension of I272-PimA-I272 was identified by the presence of four uniform I27 unfolding steps (Fig. 3, A and B, and Ref. 50). For each distinguishable step appearing in the force ramp, the step size and the force at which the event occurs are reported in the scatter plots in Fig. 3. Remarkably, PimA unfolds following multiple steps at low force. The steps observed during the ramp are well defined and can be characterized by its size and force. However, the large dispersion in step sizes and unfolding forces for PimA unfolding reflects the heterogeneity of mechanical unfolding pathways, analogous to the mechanical unfolding of molten globule states (52, 53). Similar results were obtained in the presence of 5 mm GDP (dissociation constant Kd of 0.03 μm (ΔH = −14.0 kcal·mol−1; Fig. 2), which is a well known thermal stabilizer of the enzyme. GDP binds tightly to the C-terminal domain through a network of hydrogen bonds and stabilizes the enzyme by coupling the two N- and C-terminal lobes in a closed state (27).
Fig. 4 reveals a particularly well resolved trace illustrating a possible pathway in which the theoretical extension coincides with the one-step rupture of the 163 amino acids forming the N-terminal β-sheet and the 119 amino acids forming the C-terminal β-sheet (Fig. 4A). We propose that the first four shorter steps bear out the detachment of the C-terminal helices (Fig. 4B, yellow). Then the linker and the N-terminal domain unfold in a two-step manner. The last step is compatible with the full extension of the C-terminal β-sheet. This single trace depicts one of the multiple routes that are available to fully extend PimA under force. Step sizes larger than the theoretical unfolding length of the two β-sheets would suggest the co-rupture of larger structural elements.
A single unfolding pathway of I272-PimA-I272. A, a well resolved trace illustrating a PimA unfolding pathway. The force applied increases from 0 to 250 piconewtons (pN) in 5 s. The first steps describe extensions of 6, 4, 6, 6, 22, 35, and 42 nm as the force increases, revealing intermediates along the unfolding pathway of PimA. The 42-nm step, colored in red, correlates with the unfolding of the C-terminal Rossmann fold. The four 24-nm steps occurring at higher force are colored in blue and correspond to I27 unfolding. Gray arrows display the direction of the force applied to PimA and I27 in their native states. B, topology of the proposed intermediates after unraveling the hinge (α7) and the last three α-helices (α12, α13, and α14) of PimA. The N-terminal domain unfolds first followed by the C-terminal domain. In contrast to PimA, I27 unfolds in an all-or-none unfolding pathway.
The mechanical properties of proteins are strongly influenced by their structural elements (28). To date, the strongest mechanical element corresponds to the β-sandwich topology characteristic of Ig and Ig-like proteins (54). As depicted in Figs. 3 and and4,4, I27 unfolds in an all-or-none unfolding pathway identically to proteins with similar topology such as ubiquitin and protein L (51, 55). In all these proteins, two parallel β-strands linking the N and C termini firmly protect the majority of native contacts before unraveling the rest of the protein. As soon as the parallel β-strand interactions break, the remaining structural elements unfold without any force resistance. Indeed, mutations in this structural motif in Ig or Ig-like proteins markedly reduce mechanostability (56). Interestingly, different unfolding forces have been reported for homologous Ig domains. The presence of parallel β-strands motifs does not fully explain the mechanostability of I27, and other properties have to be considered (48).
PimA topology is remarkably more complex than that of I27. PimA shares with mechanically stable proteins the presence of β-sheets in both Rossmann fold domains (14) with numerous parallel β-strands (Figs. 1 and and4),4), suggesting a strong mechanical stability. The N-terminal end folds into a β-strand (β1; residues 1–7) and makes strong interactions with β2 and β4 strands of the central β-sheet core of the first domain. In contrast, the C-terminal region terminates with the classical GT-B pattern of three long α-helices, α12 (residues 315–328), α13 (residues 330–347), and α14 (residues 349–363), with α14 crossing over to the N-terminal domain (Figs. 1 and and4).4). This hinge region is completed by residues 165–190, which include α7. Because α-helices are mechanically weak elements (57, 58), all helices are expected to unfold first at low force, allowing PimA to reach a conformation where the two Rossmann fold domains are separated. In this scenario, the GDP binding site is dislocated very early in the unfolding pathway. Therefore, when the force is applied to a PimA·GDP complex, the affinity for GDP is expected to drop following the mechanical disruption of the binding site. This model agrees with the observation that mechanical unfolding of PimA in the presence or absence of GDP gives rise to similar distributions of events (Fig. 3B).
Once the hinge has extended, parallel β-strands in both Rossmann fold domains are exposed to force (Fig. 4B). This situation induced by force is equivalent to native mechanostable proteins such as I27. However, contrary to I27, PimA unfolds with a remarkably wide distribution of step sizes at low force where no common structure dominates. Our results suggest that the presence of parallel β-strands is not sufficient to explain the mechanical properties of PimA where the presence of marginal stable structural elements within the protein allows the Rossmann fold domains to unfold along a variety of pathways. We propose that the partial unfolding of helical regions and the hinge at low force leads to a marked perturbation of the PimA energy landscape where a significant number of native contacts are disrupted before the two β-sheets are under force.
The Open and Closed States of PimA as Visualized by SAXS
To gain insight into the conformational changes that occur in PimA, the solution structures of the enzyme in its unliganded and GDP-bound forms were solved by SAXS. SAXS proved to be a powerful technique capable of providing structural information of flexible and dynamic proteins in solution (59). The interatomic distance distribution function (P(r)) computed for PimA apo and PimA in the presence of GDP is shown in Fig. 5. The molecular mass determined from the scattering data confirmed that the protein is monomeric in solution in agreement with size exclusion chromatography and analytical ultracentrifugation experiments (see “Experimental Procedures,” Table 1, and Refs. 14 and 27). It is worth noting that the protein remained in a monomeric state even in the presence of GDP.

Solution SAXS data statistics. A, scattering curves of PimA apo, PimA·GDP complex, and I27-PimA fusion polyprotein. B, P(r) function distributions of PimA apo, PimA·GDP complex, and I27-PimA fusion polyprotein. C, plot of Rg as a function of protein concentration. SAXS measurements were performed at concentrations of 2.2, 3.5, and 5.5 mg/ml PimA. Rg values at a concentration of 0 mg/ml of 28.5 and 27.5 Å were calculated for PimA apo and PimA·GDP complex, respectively.
The Rg values obtained for PimA apo and the PimA·GDP complex revealed a reduction in Rg (ΔRg) of −1.0 Å, indicating that GDP binding leads to a compaction of PimA, clearly supporting the existence of the open and closed states of the enzyme (Table 1 and Ref. 27). Several lines of experimental evidence provide strong confidence in this model. The relevance of the β-PO4 in the stabilization of the closed conformation of PimA was first studied by limited proteolysis (27). PimA was rapidly degraded after incubation with elastase. N-terminal sequencing of the two predominant species of 23 and 15 kDa revealed two exposed sites located in α9 and the connecting loop β7-β8 at the junction between N- and C-terminal domains. The α9 contains two important residues involved in GDP-Manp recognition: Asp-253 and Lys-256. Interestingly, when PimA was incubated in the presence of GDP, the enzyme was protected from the action of elastase even after 90 min, suggesting a conformational rearrangement. A close inspection of the crystal structure of the PimA·GDP complex revealed that the enzyme crystallizes in a closed conformation with the active site buried between both Rossmann fold domains (14). Two residues of the N-terminal domain were found to interact with GDP: Pro-14 in the connecting loop β1-α1 stabilizes the guanine heterocycle by a van der Waals stacking interaction, whereas Gly-16 at the top of α1 hydrogen bonds the β-PO4 of GDP (14). Moreover, Arg-196 and Lys-202, largely conserved in PimA orthologues but located in the C-terminal domain, make electrostatic interactions with the β-PO4, thus restricting its position into the catalytic site. Importantly, the α-PO4 of GDP does not interact with any particular residue from the enzyme (14). Structural comparisons of the unliganded and nucleotide-diphospho-sugar-bound forms of MshA from Corynebacterium glutamicum (CgMshA) and the glycogen synthase from Agrobacterium tumefaciens (AtGS) and E. coli (EcGS) revealed that a substantial subdomain rotation is required to achieve the closed state (20, 21, 23). The residues Gly-23, Arg-231, and Lys-236 in CgMshA; Gly-18, Arg-299, and Lys-304 in AtGS; and Gly-18, Arg-300, and Lys-305 in EcGS interact with the β-PO4 of uridine 5-diphosphate (UDP) and adenosine 5-diphosphate (ADP), respectively, and are structurally equivalent to Gly-16, Arg-196, and Lys-202 in the closed conformation of MsPimA.
The “open-to-closed” motion has also been monitored in the absence and presence of GDP and guanosine by sedimentation velocity analytical ultracentrifugation studies on pure MsPimA (27). The unliganded form of the enzyme sedimented as a single homogeneous species with an average sedimentation coefficient of 3.22 S, which is consistent with an asymmetric monomeric protein. Upon the addition of equimolar GDP, the sedimentation coefficient increased to 3.53 S, indicating the formation of a more symmetrical and compact structure. The addition of guanosine in which the α-PO4 and β-PO4 are missing did not significantly affect the sedimentation coefficient value of the mannosyltransferase, supporting the requirement of the β-PO4 to trigger the closed conformation. Moreover, isothermal titration calorimetry measurements revealed that guanosine bound to PimA with a binding constant ∼103-fold smaller than that of GDP. The GDP binding also produced a stabilizing effect on MsPimA characterized by an increment of 3.5 °C in the melting temperature value as observed by differential scanning calorimetry and circular dichroism (27). In summary, GDP induces the closing movement of PimA with a clear increment of the sedimentation coefficient value that is correlated with a reduction of the radii of gyration and accompanied by a marked stabilization of the enzyme.
Conformational Flexibility of PimA
The ab initio low resolution calculated envelopes for PimA apo and PimA·GDP complex displayed the characteristic two-domain organization of the GT-B enzymes with the deep fissure between the N- and C-terminal lobes including the catalytic center (see “Experimental Procedures,” Table 1, and Figs. 6 and and7).7). However, a comparison of the crystal structure of PimA·GDP complex with the solution structures obtained by SAXS revealed significant differences in the conformation of one of the Rossmann fold domains of the protein. To unambiguously identify the N- and C-terminal domains in the ab initio models, SAXS data were collected on a I27-PimA polyprotein where the 10-kDa I27 module was fused to the N-terminal end of MsPimA (see Table 1 and “Experimental Procedures” for details). The experimental data clearly demonstrate the position of the I27 protein in close proximity to a protruding lobe on the N-terminal domain (Table 1 and Fig. 6). It is worth noting that the different structures depicted in Fig. 7 do not directly prove the flexibility of PimA but rather stem from the potential ambiguity of ab initio reconstructions in SAXS. Nevertheless, they provide insight into the possible structural flexibility by depicting possible protein conformations compatible with the scattering data.

Low resolution models of PimA in solution. Left, average low resolution structure of PimA apo with the high resolution crystal structure of PimA·GDP complex (Protein Data Bank code 2GEK) fitted by rigid body docking. Middle, average low resolution structure of PimA·GDP complex with the high resolution crystal structure of PimA·GDP complex fitted by rigid body docking. Right, average low resolution structure of I27-PimA fusion polyprotein with the high resolution crystal structures of I27 and PimA·GDP complex fitted by rigid body docking.

Gallery of selected PimA ab initio SAXS models. Individual ab initio reconstruction runs of PimA apo (left) and PimA·GDP complex (right) showing structural variability are shown.
The N-terminal domain of PimA has been proved to participate in acceptor recognition involving important conformational changes in the protein (14, 27). Specifically, when incubated with PI, both the unliganded and the GDP-bound forms of the enzyme became highly sensitive to elastase treatment, indicating that PI triggers a significant conformational change able to modify the closed GDP-induced conformation. Furthermore, analytical ultracentrifugation experiments demonstrated that the addition of PI to the enzyme resulted in a significant change in the sedimentation coefficient values of both PimA apo and PimA·GDP complex consistent with the formation of less compact structures, which correlates with a reduction of the melting temperature by 1.5 and 0.4 °C, respectively. The crystal structure of the PimA·GDP complex displayed two regions on the N-terminal domain comprising residues 118–140 and 149–165 with particularly high B-factor values suggesting structural flexibility. These residues are located in three α-helices, α4, α5, and α6, and primarily define a hydrophobic tunnel, which is compatible with the binding of fatty acyl chains of PI or phospholipids of the mycobacterial plasma membrane. Thus, this particular region of PimA could certainly experience conformational flexibility and might explain the presence of the extra lobe in the SAXS structure.
The conformation of membrane-associated proteins can certainly be affected by binding to the lipid bilayer, thus modulating protein function. These conformational changes can be quite different in nature, ranging from subtle modifications such as the reorientation of lateral chain residues to more complex structural rearrangements (60). For example, the orientation of transmembrane domains in polytopic integral membrane-associated proteins has been proved to be dependent on membrane phospholipid composition during initial assembly as well as on changes in lipid composition postassembly (61). A dramatic conformational change has also been reported for the pore-forming toxin pneumolysin where two of the four domains refold along with the deformation of the membrane (62). Interestingly, recent studies suggest that the molten globule states play an important role in cellular processes that involve membrane interaction. The translocation of the N-terminal domain of the anthrax lethal factor across the membrane requires mechanical unfolding of a molten globule-like state (63). Moreover, the interaction of the peripheral membrane-associated α-lactalbumin to phospholipid bilayers of various compositions revealed the membrane-bound state to resemble a molten globule-like state (64). It is worth noting that the amphipathic helix C is a critical determinant for the adsorption and further integration of the α-lactalbumin into the membrane. In that sense, PimA displays an amphipathic helix on its N-terminal domain that has been proved to be essential for membrane binding and activity in vitro and in vivo (27). The occurrence of metastable states of PimA might be reflective of conformational variability of the enzyme restricted to the N-terminal domain and the association to the acceptor substrate and/or lipid bilayer. However, further studies will be required to define in detail the residues responsible for the observed structural variability.
Concluding Remarks
Remarkable progress has been made in recent years in our understanding of the catalytic mechanism and structural basis of glycosyl transfer. However, the study of the conformational changes and dynamics that govern substrate recognition and catalysis remains a major challenge in the field of GTs. The combination of single molecule force spectroscopy and SAXS techniques outlines the high variability of the PimA apo and GDP-bound states. The structural data indicate that the open and closed conformations of the enzyme are largely present in solution, but in addition, PimA experiences remarkable flexibility that undoubtedly corresponds to the N-terminal Rossmann fold domain of the protein. This domain has been proved to participate in acceptor-membrane interaction. Unfolding of PimA under force leads to heterogeneous unfolding trajectories that seem to be reflective of its intrinsic dynamic nature. According to our previously reported and current experimental data, we propose a model wherein the flexibility and conformational transitions confer adaptability of PimA to the substrates and the membrane, which seems to be of importance during the catalytic cycle.
Finally, protein dynamics and conformational changes are known to be critically important in drug discovery and development strategies (65). Although GTs play a central role in a broad range of human diseases, there is a clear lack of potent and specific GT inhibitors (66). The information presented herein highlights the conformational plasticity of PimA, information that should be considered in the rational design of GT inhibitors.
Acknowledgments
We acknowledge the European Molecular Biology Laboratory (EMBL) for synchrotron beam time allocation at Double Ring Storage (DORIS) and Positron-Electron Tandem Ring Accelerator (PETRA) storage rings, EMBL/Deutsches Elektronen-Synchrotron (DESY), Hamburg, Germany and the X33 and P12 beamline staff for technical assistance during data collection. We especially thank the BioStruct-X project for supporting access to structural biology facilities. We gratefully acknowledge all members of the Structural Glycobiology Group (Unidad de Biofísica) and the Single Protein Mechanics and Engineering Laboratory (Columbia University) for valuable scientific discussions.
*This work was supported, in whole or in part, by National Institutes of Health Grants HL61228/HL66030 (to J. M. F.). This work was also supported by European Commission Contracts LSHP-CT-2005-018923 and HEALTH-F3-2011-260872, Spanish Ministry of Science and Innovation Contract SAF2010-19096, and IKERBASQUE and the Basque Government (to M. E. G.).
This article contains supplemental Figs. 1S and 2S.
6The abbreviations used are:
- GT
- glycosyltransferase
- GDP-Manp
- guanosine 5′-diphosphate mannose
- PI
- phosphatidyl-myo-inositol
- PimA
- phosphatidyl-myo-inositol mannosyltransferase A
- SAXS
- small angle x-ray scattering
- Bis-Tris
- 2-[bis(2-hydroxyethyl)amino]-2-(hydroxymethyl)propane-1,3-diol
- NSD
- normalized spatial discrepancy
- GS
- glycogen synthase
- Rg
- radius of gyration.
REFERENCES
Articles from The Journal of Biological Chemistry are provided here courtesy of American Society for Biochemistry and Molecular Biology


