- Open Access
Myristoylation drives dimerization of matrix protein from mouse mammary tumor virus
Retrovirology volume 13, Article number: 2 (2016)
Myristoylation of the matrix (MA) domain mediates the transport and binding of Gag polyproteins to the plasma membrane (PM) and is required for the assembly of most retroviruses. In betaretroviruses, which assemble immature particles in the cytoplasm, myristoylation is dispensable for assembly but is crucial for particle transport to the PM. Oligomerization of HIV-1 MA stimulates the transition of the myristoyl group from a sequestered to an exposed conformation, which is more accessible for membrane binding. However, for other retroviruses, the effect of MA oligomerization on myristoyl group exposure has not been thoroughly investigated.
Here, we demonstrate that MA from the betaretrovirus mouse mammary tumor virus (MMTV) forms dimers in solution and that this process is stimulated by its myristoylation. The crystal structure of N-myristoylated MMTV MA, determined at 1.57 Å resolution, revealed that the myristoyl groups are buried in a hydrophobic pocket at the dimer interface and contribute to dimer formation. Interestingly, the myristoyl groups in the dimer are mutually swapped to achieve energetically stable binding, as documented by molecular dynamics modeling. Mutations within the myristoyl binding site resulted in reduced MA dimerization and extracellular particle release.
Based on our experimental, structural, and computational data, we propose a model for dimerization of MMTV MA in which myristoyl groups stimulate the interaction between MA molecules. Moreover, dimer-forming MA molecules adopt a sequestered conformation with their myristoyl groups entirely buried within the interaction interface. Although this differs from the current model proposed for lentiviruses, in which oligomerization of MA triggers exposure of myristoyl group, it appears convenient for intracellular assembly, which involves no apparent membrane interaction and allows the myristoyl group to be sequestered during oligomerization.
Mouse mammary tumor virus (MMTV), a causative agent of mammary breast cancer and T cell lymphomas in mice (for review see ), is a representative of the Betaretrovirus genus. In contrast to HIV-1 and other retroviruses that assemble immature particles during budding at the plasma membrane (PM), betaretroviruses assemble immature particles within the cytoplasm prior to transport to the PM for budding. In general, the hexameric lattice of immature retroviral particles is formed by association of 1500 to 2500 molecules of the multidomain structural polyprotein Gag. The N-terminal domain of Gag, matrix (MA) protein, is responsible for Gag targeting and PM binding. During or shortly after budding, the viral protease-dependent cleavage of Gag initiates rearrangement of individual mature proteins, after which MA molecules remain associated with the enveloping membrane and form the outer shell of the viral particle .
In most retroviruses, the N-terminal glycine residue of the MA domain of Gag is co-translationally modified with myristic acid (reviewed in ). Myristoylation of MA and the highly basic patch of amino acids located near its N-terminus mediate the interaction of Gag with the PM. The myristoyl group serves as a hydrophobic anchor associated with the lipidic acyl chains in the PM, while the positively charged surface patch in MA interacts electrostatically with the polar heads of phospholipids. Both the myristoyl group and the basic patch function in concert to facilitate efficient MA membrane anchoring (for review see ). In HIV-1 and other viruses that assemble particles at the PM, prevention of myristoylation by mutation of the N-terminal glycine to alanine (G2A) dramatically reduces binding of Gag to PM and inhibits viral particle formation [5–9]. However, in betaretroviruses such as Mason-Pfizer monkey virus (M-PMV) and MMTV, the analogous G2A mutation does not interfere with the process of intracytoplasmic assembly but completely blocks the transport of viral particles to the PM [10, 11]. The MA domain not only targets the PM and mediates the association of viral particles, but it also dictates the lipid-binding specificity (for review see [12, 13]).
In addition to its role in promoting the association with the PM, some other functions have been attributed to the MA domain of Gag. In M-PMV, a highly conserved short stretch of amino acids within MA domain, termed the cytoplasmic targeting/retention signal (CTRS), interacts with Tctex-1, a component of the dynein motor machinery , and is responsible for intracellular targeting of Gag to the pericentriolar region for assembly. Similarly, CTRS in the MMTV MA domain was recently found to mediate pericentriolar targeting of Gag . Several mutations within the M-PMV MA domain can abrogate transport of assembled particles to the PM in a manner similar to the G2A mutation in Gag, and these changes in MA might negatively affect Gag/Env interaction at the recycling endosome and prohibit usage of the vesicular transport system [16–18]. The lentiviral MA domain may be directly involved in incorporation of Env into virions, at least for HIV-1 .
Structures of MA protein from ten retroviruses have been resolved by NMR or X-ray crystallography [20–29]. Despite the low sequence homology of these proteins, their overall three-dimensional organization is remarkably similar, consisting of a globular core composed of four or five α-helices (reviewed in ). However, structures of myristoylated forms [myr(+)] have been reported only for MA proteins from the lentiviruses HIV-1 , HIV-2 , feline immunodeficiency virus (FIV) , and the betaretrovirus M-PMV . The solution structure of myr(+) HIV-1 MA revealed an equilibrium between a monomeric form with the myristoyl group sequestered inside the protein core and a trimeric form with exposed myristoyl groups . Transition between the monomeric and trimeric states of myr(+) HIV-1 MA is entropically regulated, and higher protein concentration stimulated and stabilized formation of the trimeric form . In contrast, nonmyristoylated [myr(−)] HIV-1 MA appeared mainly in the monomeric state in solution  but in the trimeric state in crystals . Despite the structural similarity between the HIV-1 and HIV-2 MA proteins, myristoylation does not influence the trimerization ability of HIV-2 MA . However, betaretroviral MA from M-PMV readily forms oligomers (dimers and trimers) in its myr(−) form , while myr(+) MA is monomeric in solution . Based on extensive studies of HIV-1 MA, it is generally accepted that the myristoyl group can adopt a sequestered conformation within the globular head of MA as well as an exposed conformation where it is more accessible and can facilitate membrane association. The process of transition between these two states, termed the myristoyl switch, can be influenced by several factors, including protein concentration, oligomerization status, and pH [30, 35].
Here, we investigated the role of myristoylation on oligomerization of MMTV MA and determined the crystal structure of myr(+) MMTV MA, which represents the first X-ray structure of a myr(+) retroviral MA protein. Our structural data reveal the sequestered conformation of the myristoyl group bound into the MMTV MA dimer interface, and analytical ultracentrifugation and molecular dynamics analyses show its crucial role for MA dimerization. Based on the structural and molecular modeling data, we constructed MMTV MA mutant with a reduced ability to form dimers in vitro. Identical mutations were also tested in the context of MMTV virus and their impact upon assembly and transport of immature particles was analyzed.
Myristoylation stimulates dimerization of MMTV MA in solution
Typically, myristoyl group participates in protein subcellular localization by facilitating protein–membrane and protein–protein interactions (reviewed in ). To determine the influence of myristoylation on self-association of the MA domain of MMTV Gag, we performed analytical ultracentrifugation using bacterially expressed and purified myr(+) and myr(−) MA, containing five C-terminally attached amino acids from the pp21 domain followed by a His-tag . The purified myr(+) MA contained less than 0.1 % of myr(−) MA as documented by mass spectrometry analysis (Additional file 1: Figure S1). Sedimentation velocity experiments were performed with recombinant myr(+) and myr(−) MA proteins at two different concentrations (0.5 and 5.0 mg/mL), and we analyzed their tendency to form oligomers by fitting the sedimentation data using a continuous size distribution model. The comparison of all resultant size distributions is shown in Fig. 1. At the lower concentration, myr(−) MA formed one discrete particle with sedimentation coefficient s20,w = 1.57 S. This value corresponds well with the mass of the MA monomer, pointing to a moderately elongated particle with approximate dimensions of 2–3 × 4–6 nm. At the higher concentration, myr(−) MA formed two particles, corresponding to a monomer (s20,w = 1.69 S; 74 %) and a dimer (s20,w = 2.12 S; 26 %). In contrast, at the lower concentration, myr(+) MA formed only dimers with s20,w = 2.12 S. At the higher concentration, myr(+) MA predominantly formed dimers (s20,w = 2.27 S; 98 %) but also formed a particle corresponding to a tetramer (s20,w = 3.69 S; 2 %). Thus, we conclude that the ability of MMTV MA to form dimers is stimulated by its myristoylation.
The dimerization of MMTV MA driven by its myristoylation also provides a likely explanation for the separation of His-tagged myr(+) and myr(−) MMTV MAs on an IMAC column that we described previously . A dimer of MA has two His-tags which results in higher affinity to the IMAC column. The efficiency of the separation is influenced by the difference in the ability of myr(+) and myr(−) MAs to form dimers. The formation of myr(+):myr(−) mixed dimers decreases the separation efficiency.
Structure of myristoylated MMTV MA
To provide structural evidence for the role of the myristoyl group in MA dimerization, we determined the structure of myr(+) MA prepared by the procedure described by Dolezal et al. . The protein crystallized in the P3 1 21 space group with two molecules in the asymmetric unit. The structure was solved by single isomorphous replacement with anomalous scattering (SIRAS) using potassium iodide. The final dataset was collected to 1.57 Å resolution, and the final model was refined to R = 22.6 % (Rfree = 26.3 %). Crystal parameters, data collection statistics, and refinement statistics are summarized in Table 1. In the refined structure, myristoyl groups and their connections to the N-terminal glycines were well-defined in the electron density map, as were the first 91 of 111 amino acid residues. The last 20 residues (the C-terminus of the protein and the His-tag) lacked electron density, likely due to a disordered conformation. The structure was deposited in the Protein Data Bank under PDB ID 4ZV5.
The asymmetric unit contains a dimer formed by two MA molecules. The structures of both MA subunits are virtually identical except for the position of the myristoyl group and the first five amino acids; the residual mean square deviation of Cα atoms for residues 7–92 is 0.306 Å (the N-terminal glycine is designated as residue 2). The monomeric subunit (Fig. 2) is a helical bundle comprising five alpha helices and a single-turn 310-helix (helix 5). The CTRS region (spanning residues 45–62) is accessible on the surface of the molecule in an orientation similar to that observed for M-PMV MA [14, 32]. Despite the low sequence similarity of retroviral MA proteins (Table 2), which precluded the use of the molecular replacement approach in the structure determination, the position and orientation of the helices follow the structural pattern of other retroviral MA proteins.
Interestingly, the myristoyl group of one of the dimer-forming MA molecules is not buried in the hydrophobic core of the protein chain to which it is covalently attached, as is common for the sequestered form of myr(+) MA proteins [28, 30–32]. Instead, the myristoyl group is sequestered in the hydrophobic pocket formed by the interface of both MA molecules in the asymmetric dimer unit (Fig. 3). This interface is formed predominantly by helices 1 and 3. However, the myristoyl group of the second MA molecule in the dimer extrudes into the neighboring dimer unit. Therefore, the two nearest dimer units in the crystal are interconnected by two myristoyl groups, forming a tetramer of MA molecules around the two-fold axis of the P3 1 21 space group (Fig. 3a, b). Regardless of the dimer unit to which the myristoyl groups are covalently attached, their positions in the hydrophobic interface are virtually symmetrical (Fig. 3c).
To analyze the relevance of the oligomeric organization in the protein crystal, we used PDBePISA , an interactive tool for the exploration of macromolecular interfaces. The PISA software detected several interfaces, of which the largest was 945 Å2 per MA monomer. The dominant interface lies between the two MA molecules in the single asymmetric unit. We did not identify any significant interface between MA dimers, with the exception of the interface represented solely by interconnecting myristoyl groups.
Because myr(+) MMTV MA predominantly forms dimers in solution, as confirmed by analytical ultracentrifugation analysis, the myristoyl swapping between the neighboring asymmetric units appears to be a specific result of crystal packing. In solution, the extruded myristoyl group is likely inserted into the hydrophobic pocket of its own dimer, forming a symmetric dimer with swapped myristoyl groups.
Molecular dynamics simulation and calculation of interaction energies support the myristoyl swapping model in MMTV MA dimer formation
Based on our experimental results, we propose that the swapping of myristoyl groups within the dimer is important for dimerization of MMTV MA. To support this proposal, we performed molecular dynamics (MD) simulations and calculated interaction energies for four theoretically possible MMTV MA dimers (Fig. 4). All the MD trajectories became stable after 6 ns of simulation. The root-mean-square deviations (RMSD) of the protein backbone fluctuated around 1 Å, while those of the myristoyl groups reached up to 1.7 Å in the case of dimer 3 (D3). Interaction energy calculations using the MM-PBSA approach showed the order of binding as follows: D2 (−130 kcal/mol, swapped myristoyls) >D1 (−108 kcal/mol, the asymmetric unit) >D3 (−89 kcal/mol, unswapped myristoyls) >D4 (−68 kcal/mol, nonmyristoylated). The results of the MD simulation indicate that myristoyl groups are energetically more favorable when swapped between the two dimer-forming MA molecules (Fig. 4, D2) than when unswapped (Fig. 4, D3). Furthermore, the absence of a single myristoyl group within the dimer interface causes the interaction energy to drop (Fig. 4, D1). To evaluate the relative contribution of individual residues to the dimerization, we performed a per-residue decomposition of the MM-PBSA interaction energies. The results, summarized in Table 3, show that the summed contribution of the myristoyl groups is three-fold higher for swapped myristoyl groups (−35.3 kcal/mol) than for unswapped myristoyl groups (−11.7 kcal/mol), which is comparable to the summed contribution of the most strongly contributing amino acid residues (−10 to −18 kcal/mol). These results support myristoyl swapping within the dimer as a driver for dimerization of MMTV MA.
Mutations within the dimerization interface impair MMTV MA dimerization
To verify contribution of residues to MA dimerization, we mutated selected residues oriented to the interior of dimerization interface (Fig. 5). Because many residues from dimerization interface are located in the regions important for proper biological function of MMTV MA (myristoylation signal at the N-terminus, and cytoplasmic targeting/retention signal spanning residues 45–62 ), our choice of the residues available for mutation was restricted. We selected residues L11, F12, and V15 and prepared a triple mutant L11N/F12N/V15N (MA-3N) in both myristoylated and nonmyristoylated form. Analytical ultracentrifugation analyses of purified mutant myr(+) and myr(−) MA mutants were performed at two concentrations (0.5 and 5.0 mg/mL) under conditions identical to those used for WT MA. Figure 6 shows that mutations L11N, F12N, and V15N impaired MA dimerization. Contrary to the WT myr(−) MA, which partially formed dimers at the higher concentration, myr(−) MA-3N formed only monomers at both concentrations (s20,w = 1.53 S). At the higher concentration, myr(+) MA-3N formed only dimers (s20,w = 2.20 S) whereas wild-type myr(+) MA also formed a small amount of tetramers (2 %). At the lower concentration, myr(+) MA-3N formed monomers (s20,w = 1.81 S; 52 %) and dimers (s20,w = 2.39 S; 48 %) in contrast to wild-type myr(+) MA, which formed only dimers. The effect of L11N/F12N/V15N mutations on dimerization of myr(+) MA was lesser compared to the contribution of myristoyl group to this process. These results show that myristoylation itself is required, but not sufficient for efficient MMTV MA dimerization.
Mutations affecting MMTV MA dimerization also influence extracellular particle production
To analyze impact of the mutations affecting MA dimerization on MMTV particle production, we performed pulse-chase experiments in transiently transfected human embryonic kidney 293T cells. For MMTV particle production, we used the chimeric proviral vector pSMt-HYB/D26A, which contains the heterologous M-PMV LTR promoter and a mutation inactivating the viral protease . This construct facilitates detection of intracytoplasmic particles and it will be denoted hereafter as wild-type (WT).
Following transfection into 293T cells, WT and MA-3N mutant viruses were analyzed by metabolic labeling for their ability to express stable Gag polyproteins and produce extracellular immature particles (Fig. 7). Particle-producing cells were pulse-labeled for 1 h and chased for 12 h. Similar levels of stable Gag precursors were synthesized in all transfected cells and extracellular particles were observed for WT and mutant virus as well. As shown in Fig. 7a, the amount of pelletable structures in media produced by MA-3N mutant was lower than those produced by WT virus. To measure the relative efficiency of particle formation, the amount of extracellular assembled Gag as a fraction of total Gag expressed in cells was determined for each construct and then compared with that of WT (Fig. 7b). MA-3N virus produced particles with approximately half the efficiency of WT.
To verify the results obtained in human 293T cells, we used the rat mammary cell line RBA, which is suitable for production of MMTV upon transient transfection . Because the expression level of MMTV proteins is significantly lower in RBA cells than in 293T cells, we examined the accumulation of assembled Gag in culture media by Western blot analysis after 48 h of production in transfected cells (Fig. 8). Proteins were separated on 15 % SDS PAGE gels and detected with a rabbit anti-MMTV capsid polyclonal antibody. Under these conditions, we were able to detect extracellular pelletable structures for both WT and MA mutant and as expected, the particle release of MMTV MA-3N from RBA cells was reduced. These results suggest that MMTV bearing mutations L11N, F12N and V15N in the MA domain exhibits reduced efficiency of extracellular particle production.
To further investigate the effect of impaired MA dimerization on extracellular particle production, we analyzed transiently transfected 293T cells expressing WT and MA-3N virus by transmission electron microscopy (TEM) (Fig. 9). In contrast to the WT producing cells, only a minimum of particles budding from the plasma membranes and accumulation of both isolated and clusters of assembled particles were observed in cells expressing the MA-3N mutant virus. Aberrantly formed particles were detected with a low frequency similar for both WT and MA-3N MMTV. These observations suggest a reduced efficiency of the intracellular transport as an explanation for the partial defect of extracellular particle production in MA-3N mutant virus.
Here, we demonstrated that myristoylation stimulates dimerization of MMTV MA and identified structural requirements for this interaction. Molecular dynamics modeling based on our experimentally determined structure of myr(+) MA confirmed the key role of the myristoyl group for dimer formation and revealed energetically favorable binding of the myristoyl group within the dimer interface. Myristoylated MMTV MA shares structural similarity with other retroviral MAs, but it differs in the orientation of the sequestered myristoyl group. In other myr(+) MAs characterized to date, the myristoyl group is sequestered in the hydrophobic core of the molecule to which it is covalently attached [28, 30–32]. In contrast, the myristoyl group of MMTV MA is entirely buried in the dimerization interface between two MA molecules. Moreover, dimerization of MMTV MA supports the sequestered conformation of the myristoyl, which contrasts with the myristoyl switch model based on HIV-1 MA, in which oligomerization supports its extrusion .
Despite the fact that extrusion of one myristoyl group from one MMTV MA dimer unit to the neighboring asymmetric unit appears to be the result of crystal packing, it demonstrates an important aspect related to MA function. The process by which one MMTV MA molecule adopts an exposed myristoyl group conformation is analogous to the myristoyl switch model, in which the MA domain of a retroviral Gag polyprotein interacts with phospholipid membranes  and stimulates Gag multimerization. Thus, our results indicate that myr(+) MMTV MA can adopt both sequestered and exposed conformations analogous to those described for myr(+) HIV-1 MA.
The MMTV MA dimer interface forms a deep hydrophobic pocket that is energetically favorable for myristoyl binding. Our data from analytical ultracentrifugation clearly show that the myristoyl group stimulates MA dimer formation at high and low protein concentrations. The absence of the myristoyl group reduces the amount of MA dimer in solution even at high protein concentration (5 mg/mL). At a ten-fold lower concentration, myr(−) MA does not dimerize. It appears that the myristoyl group actively participates in formation of the dimer interface by interconnecting both monomeric MMTV MA units and that hydrophobic interactions are primarily responsible for dimer stability. Accordingly, no homotypic interactions were detected when myr(−) MMTV MA protein fused at the N-terminus with a reporter domain was tested for self-association in the yeast two-hybrid system . Mutations of residues from dimerization interface (L11N/F12N/V15N), designed according to the MA crystal structure and molecular dynamic calculations, resulted in decreased MA dimerization. Interestingly, these mutations influenced also transport of assembled particles to PM and resulted in decreased production of extracellular particles. The current model describing the relocation of betaretroviral particles to PM involves anterograde transport mediated by microtubule-associated Env-containing vesicles [18, 41–43]. However, a detailed mechanism of this transport is not yet fully understood. These questions remain to be investigated.
Myristoylation has been demonstrated to influence oligomerization of HIV-1 MA in solution; myr(−) MA is monomeric and myr(+) MA resides in monomer-trimer equilibrium . Despite a known tendency of retroviral MAs to form dimers or trimers at high protein concentrations, they are predominantly monomeric in solution. To date, the only oligomeric state of a retroviral MA protein proven to be biologically relevant is the trimer initially discovered in the crystal structures of myr(−) simian immunodeficiency virus (SIV) MA  and myr(−) HIV-1 MA . NMR spectroscopy showed that in solution, myr(+) HIV-1 MA also has a tendency to form trimers analogous to those observed in the crystal structures . Although it was shown that myristoylation positively influences trimerization , the contribution of the myristoyl group to the trimerization of HIV-1 MA is of a different nature than our model proposes for dimerization of MMTV MA. In the HIV-1 trimer, the myristoyl groups are exposed and positioned outside of the trimer interface and thus cannot participate in its formation (Fig. 10).
In betaretrovirus assembly, there is no apparent requirement for MA interaction with the PM, and intracytoplasmic formation of immature particles instead relies on the internal scaffolding functions of domains located N-terminal to capsid protein [39, 40, 44–47]. In contrast to HIV-1 MA, the oligomerization capacity of betaretroviral M-PMV MA, which is capable of forming dimers and trimers in its myr(−) form [34, 48], is not stimulated by myristoylation . Instead, the structural changes caused by the sequestered myristoyl negatively affect the oligomerization of the MA molecule. It was suggested that the oligomerization capacity of M-PMV MA is restored by conformational changes induced by the binding of MA to the membrane . These differences from myristoyl-stimulated dimerization of MMTV MA suggest that betaretroviruses may have developed different strategies to sequester the MA myristoyl group until it needs to adopt an exposed conformation for membrane binding.
We provide experimental data demonstrating that dimerization of MMTV MA is stimulated by its myristoylation. We clarified the molecular mechanism of this stimulation by determining the crystal structure of the myr(+) MMTV MA protein. Based on structural data and molecular modeling, we propose a model for symmetric dimer formation of myr(+) MMTV MA in solution. Our mechanism describing how dimerization of MMTV MA relies on swapped myristoyl groups represents a novelty not only among retroviral MA proteins but also for dimerization of any myristoylated protein. Our model suggests that oligomerization of betaretroviral MA protein is accompanied by sequestration of the myristoyl group that participates in the formation of the interaction interface. Although it contrasts with the model proposed for HIV-1 and other retroviruses that assemble with membrane support, our model is well-suited for intracytoplasmic betaretroviral assembly. Biological experiments further indicated that proper MA dimerization is important for efficient MMTV particle release.
Plasmid pET22-MA-His used for the bacterial expression of MMTV MA protein was described previously . Mutations of residues 11L, 12F and 15V to asparagines were introduced by PCR into pSMtv  plasmid containing MMTV gag gene sequence using forward mutagenic primer 5′-TGGGGGTCTCGGGCTCAAAAGGGCAGAAAAACAATGTTTCTAATTTACAAAGGCTCCTCTCAGAGAGGGGTC-3′ in combination with reverse primer 5′-TTTCAATGGCAGCGGTTCCC-3′ and pSMtv as template. Resulting PCR product was cleaved with AvaI–PflMI and ligated into pSMtv vector digested in the same way to produce pSMtv-MA-3N shuttle vector. Mutated gag sequence was used as a template for PCR using primers MA-Nde-F and MA-6H-Xho-R and the resulting mutant MA sequence was cleaved with NdeI–XhoI enzymes and ligated into identicaly prepared pET22b vector for the bacterial expression. The primer sequences used for this subcloning will be available upon request. In parallel, the complete gag sequence carrying the desired mutations was subcloned from pSmtv-MA-3N plasmid using KasI–XbaI enzymes and subcloned into identicaly treated MMTV proviral construct pSMt-HYB/D26A . All DNA segments resulting from PCR were verified by sequencing.
Expression and purification of recombinant MMTV MA protein
Recombinant myr(−) and myr(+) MMTV MA proteins were expressed and purified as previously described . The molecular masses of the proteins were analyzed by mass spectrometry (MALDI TOF/TOF) as previously described . To approximately determine the minimal detectable amount of myr(−) MA in the sample of myr(+) MA, we mixed a sample of myr(+) MA which we considered “100 %” myristoylated with a sample of myr(−) in the ratio 1000:1 (2:0.002 mg/mL) and analyzed the mixture by mass spectrometry (Additional file 1: Figure S1B). The 0.1 % concentration of myr(−) MA was detectable in the mixture and therefore we can state that the sample of myr(+) contained less than 0.1 % of myr(−) MA, which we suppose to not interfere with the subsequent analyses. After purification, the proteins were dialyzed against 50 mM Na2HPO4, 300 mM NaCl, 2.5 mM TCEP, pH 6.0 (buffer A), and concentrated to an appropriate concentration by ultrafiltration using an Amicon Ultra 3K device (Millipore).
Sedimentation velocity analyses were performed using a ProteomeLab XL-I analytical ultracentrifuge equipped with an An50Ti rotor (Beckman Coulter). Protein samples were concentrated to 0.5 or 5.0 mg/mL in buffer A, which was also used as a reference buffer. The experiments were carried out at 50,000 rpm and 4 °C; 200 absorbance scans were recorded at 6 min intervals with 30 µm spatial resolution at 280 nm for the diluted samples and at 250 nm for the concentrated samples with 12 and 3 mm Epon centerpieces, respectively. Buffer density and protein partial specific volume were estimated in SEDNTERP 1.09 (http://www.jphilo.mailway.com). Data were analyzed with SEDFIT 14.1  using a c(s) continuous size distribution model. The approximate ratio of monomeric and dimeric forms was calculated from the peak areas.
Myristoylated MMTV MA was concentrated to a final concentration of 10 mg/mL. Initial crystallization trials were performed with the help of a Crystal Gryphon crystallization workstation (Art Robbins Instruments) by sitting drop vapor diffusion method at 19 °C in 96-well plates; 0.2 μL protein solution was mixed with 0.2 μL reservoir solution and the mixture was equilibrated over 200 μL reservoir solution. The PEGs Suite and JSCG Core I Suite (QIAGEN) were used for the initial crystallization condition screen. Initial microcrystals appeared in several days under the following conditions: 0.2 M potassium chloride and 20 % PEG 3350. Further optimization involved changing to the hanging drop mode in 24-well crystallization plates (EasyXtal DG-Tool, QIAGEN). Final crystals were obtained by mixing 3 μL MMTV MA complex solution with 1 μL reservoir solution composed of 0.2 M potassium chloride and 20 % PEG 3350 as precipitant, and these were directly cryocooled in liquid nitrogen. For the phasing experiment, crystals were soaked for 10 min in a solution of 0.2 M potassium iodide and 20 % PEG 3350 and cryocooled in liquid nitrogen.
Data collection and structure determination
Diffraction data for phasing were collected using a 1.5 Å wavelength for native and iodide-soaked crystals. Data were collected to 1.9 Å resolution at 100 K for both datasets. The native dataset for the final structure was collected using a 0.9184 Å wavelength at 100 K. All datasets were collected at the MX14.2 beamline at BESSY, Berlin, Germany , and processed using the XDS program . The structure was solved with SHELXC/D/E programs  using HKL2MAP GUI  or GUI for SHELX programs . Macromolecular phasing was performed according to the SIRAS method, using datasets from the native crystal and crystal soaked in potassium iodide (both measured using a 1.5 Å wavelength). The initial model was improved and rebuilt with the program Buccaneer . This was followed by manual rebuilding with COOT  and refinement with Refmac5  using the native dataset measured at 1.5 Å wavelength. At the end of this stage, the refinement statistics were as follows: R = 26.0 % (Rfree = 32.4 %). The final model, calculated from the dataset measured at 0.9184 Å wavelength, was refined to R = 22.6 % (Rfree = 26.3 %). The crystal parameters, data collection statistics, and refinement statistics are summarized in Table 1.
Molecular dynamics simulation and calculation of interaction energies
Four possible dimers were derived from the crystal structure of myr(+) MMTV MA (Fig. 4). The first dimer (D1) was the asymmetric unit without any modifications. The second dimer (D2) was derived from the asymmetric unit by connecting the myristoyl group from monomer 2B to monomer 1B (Fig. 3a), which yielded a pseudo-symmetric dimer with swapped myristoyl groups. The third dimer (D3) was derived from D2 by swapping the connections of the proteins to the myristoyl groups, which yielded a dimer with unswapped myristoyl groups. The fourth dimer (D4) was derived from D2 by removing the myristoyl groups.
The systems were prepared for molecular dynamics (MD) simulations essentially as previously described , with a few modifications. The two histidine residues (both located on the protein surface far from the interface) were modeled in the Nε monoprotonated state. This setup was created in the LEaP module of AMBER14 . The amino acid parameters were taken from the ff14SB AMBER force field , and the myristoyl parameters were obtained using RESP charge fitting  and GAFF parameters . Chloride counterions were added to the solvated dimers to neutralize the total charge (four for D1–D3 and six for D4).
The initial relaxation included optimization of the myristoyl and Gly2 residues in 250 cycles of steepest descent. The other minimization steps, warming, and equilibration MD followed a previously described protocol . The length of production MD was 10 ns. Monomer–monomer interaction energies were calculated using the MM-PBSA module of AMBER14 on 100 snapshots of the last 4 ns of the simulation. The Poisson-Boltzmann equation was solved, and Debye–Hückel screening with 150 mM ionic strength was used. The solvent-accessible surface-area-dependent term was employed for nonelectrostatic solvation free energies. The interaction energies were decomposed on a per-residue basis so that 1–4 interactions were added to either electrostatic or van der Waals contributions.
Cell lines and transfections
Human kidney embryonic cells 293T (ATCC) and rat mammary cells RBA (ATCC) were maintained in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10 % fetal bovine serum and 20 mM l-glutamine (Sigma Aldrich). 293T cells were transfected by FuGENE HD (Roche Molecular Biochemicals) reagent and RBA cells were transfected using Lipofectamine LTX reagent (Invitrogen), both according to the manufacturer’s instructions.
Metabolic labeling and immunoprecipitation
Cells were grown on 60 mm culture dishes, 24 h post-transfection were 3 times washed with PBS buffer and then pulse-labeled for 60 min at 37 °C in 1 mL of DMEM without methionine and cysteine supplemented with 100 µCi of [35S]methionine-[35S]cysteine protein-labeling mix (MP Biochemicals). The label was chased for the desired period of time by replacing the labeling medium with complete DMEM. For analysis of total intracellular viral proteins, the cells were lysed in buffer containing 0.15 M NaCl, 50 mM Tris–HCl (pH 7.5), 1 % Triton X-100, and 1 % deoxycholate (Lysis Buffer). Cellular debris and nuclei were removed by centrifugation for 1 min at 14,000 rcf. The lysate was then adjusted to 0.1 % SDS and the cell-associated viral proteins were immunoprecipitated by rabbit polyclonal anti-CA antiserum for 1 h and after addition of protein A immobilized on Sepharose beads (Invitrogen) the suspension was incubated for 2 h. For analysis of released virus particles, tissue culture medium was first filtered through a 0.45 μm filter. Particles were then collected by centrifugation through a 20 % w/v sucrose cushion at 35,000 rpm for 1 h in a SW 41 Ti rotor (Beckman Coulter). Pelleted virions were dissolved in Lysis Buffer containing additional 0.1 % SDS and the virus-associated proteins were immunoprecipitated as described for cell lysates. Immunoprecipitates were examined by separation on 10 % SDS-PAGE gels followed by phosphorimager analysis. Band intensities for Gag were acquired on a Typhoon system using ImageQuant software (Amersham).
Transiently transfected RBA cells were grown on 100 mm culture dishes and harvested 48 h post-transfection. Trypsinized cells were pelleted for 1 min at 14,000 rcf and resuspended in 0.1 ml PBS and 0.2 ml of SDS protein loading buffer (PLB 2×) Virions from the culture supernatants were collected as described above for metabolically labeled proteins and resulting pellets were resuspended in 30 μl of PLB (2×). Analyzed proteins were separated on 10 % SDS-PAGE gel, blotted onto a nitrocellulose membrane, and detected by using rabbit polyclonal anti-CA antibody.
Examined cells were washed in PBS and fixed by 3 % glutaraldehyde in 0.1 M cacodylate buffer, pH 7.4. The cells were postfixed with 1 % osmium tetroxide, dehydrated in graded ethanol solutions and embedded in epoxy resin AGAR 100. Ultrathin sections (70 nm) were stained with uranyl acetate and lead citrate. The samples were analyzed with a JEOL JEM-1200EX electron microscope operated at 60 kV.
Availability of supporting data
The structure of myristoylated MMTV MA protein was deposited in the Protein Data Bank under PDB ID 4ZV5.
cytoplasmic targeting/retention signal
human immunodeficiency virus type 1
mouse mammary tumor virus
Mason-Pfizer monkey virus
Ross SR. Mouse mammary tumor virus molecular biology and oncogenesis. Viruses. 2010;2:2000–12.
Swanstrom R, Wills JW: Synthesis, assembly, and processing of viral proteins. In: Coffin JM, Hughes SH, Varmus. HE, editors. Retroviruses. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 1997. p. 263–334.
Maurer-Stroh S, Eisenhaber F. Myristoylation of viral and bacterial proteins. Trends Microbiol. 2004;12:178–85.
Dick RA, Vogt VM. Membrane interaction of retroviral Gag proteins. Front Microbiol. 2014;5:1–11.
Zhou W, Parent LJ, Wills JW, Resh MD. Identification of a membrane-binding domain within the amino-terminal region of human immunodeficiency virus type 1 Gag protein which interacts with acidic phospholipids. J Virol. 1994;68:2556–69.
Bryant M, Ratner L. Myristoylation-dependent replication and assembly of human immunodeficiency virus 1. Proc Natl Acad Sci USA. 1990;87:523–7.
Göttlinger HG, Sodroski JG, Haseltine WA. Role of capsid precursor processing and myristoylation in morphogenesis and infectivity of human immunodeficiency virus type 1. Proc Natl Acad Sci USA. 1989;86:5781–5.
Freed EO, Orenstein JM, Buckler-White AJ, Martin MA. Single amino acid changes in the human immunodeficiency virus type 1 matrix protein block virus particle production. J Virol. 1994;68:5311–20.
Ono A, Freed EO. Binding of human immunodeficiency virus type 1 Gag to membrane: role of the matrix amino terminus. J Virol. 1999;73:4136–44.
Rhee SS, Hunter E. Myristylation is required for intracellular transport but not for assembly of D-type retrovirus capsids. J Virol. 1987;61:1045–53.
Zábranský A, Hadravová R, Stokrová J, Sakalian M, Pichová I. Premature processing of mouse mammary tumor virus Gag polyprotein impairs intracellular capsid assembly. Virology. 2009;384:33–7.
Hamard-Peron E, Muriaux D. Retroviral matrix and lipids, the intimate interaction. Retrovirology. 2011;8:15.
Chukkapalli V, Ono A. Molecular determinants that regulate plasma membrane association of HIV-1 Gag. J Mol Biol. 2011;410:512–24.
Vlach J, Lipov J, Rumlová M, Veverka V, Lang J, Srb P, Knejzlík Z, Pichová I, Hunter E, Hrabal R, Ruml T. D-retrovirus morphogenetic switch driven by the targeting signal accessibility to Tctex-1 of dynein. Proc Natl Acad Sci USA. 2008;105:10565–70.
Zhang G, Sharon D, Jovel J, Liu L, Wine E, Tahbaz N, Indik S, Mason A. Pericentriolar targeting of the mouse mammary tumor virus GAG protein. PLoS One. 2015;10:e0131515.
Rhee SS, Hunter E: Amino acid substitutions within the matrix protein of type D retroviruses affect assembly, transport and membrane association of a capsid. EMBO J 1991;10:535–546.
Stansell E, Tytler E, Walter MR, Hunter E: An early stage of Mason-Pfizer monkey virus budding is regulated by the hydrophobicity of the Gag matrix domain core. J Virol 2004;78:5023–5031.
Stansell E, Apkarian R, Haubova S, Diehl WE, Tytler EM, Hunter E: Basic residues in the Mason-Pfizer monkey virus gag matrix domain regulate intracellular trafficking and capsid-membrane interactions. J Virol 2007;81:8977–88.
Murakami T, Freed EO. Genetic evidence for an interaction between human immunodeficiency virus type 1 matrix and alpha-helix 2 of the gp41 cytoplasmic tail. J Virol. 2000;74:3548–54.
Massiah MA, Starich MR, Paschall C, Summers MF, Christensen AM, Sundquist WI. Three-dimensional structure of the human immunodeficiency virus type 1 matrix protein. J Mol Biol. 1994;244:198–223.
Rao Z, Belyaev AS, Fry E, Roy P, Jones IM, Stuart DI. Crystal structure of SIV matrix antigen and implications for virus assembly. Nature. 1995;378:743–7.
Matthews S, Mikhailov M, Burny A, Roy P. The solution structure of the bovine leukaemia virus matrix protein and similarity with lentiviral matrix proteins. EMBO J. 1996;15:3267–74.
Christensen AM, Massiah MA, Turner BG, Sundquist WI, Summers MF. Three-dimensional structure of the HTLV-II matrix protein and comparative analysis of matrix proteins from the different classes of pathogenic human retroviruses. J Mol Biol. 1996;264:1117–31.
Conte MR, Klikova M, Hunter E, Ruml T, Matthews S. The three-dimensional solution structure of the matrix protein from the type D retrovirus, the Mason-Pfizer monkey virus, and implications for the morphology of retroviral assembly. EMBO J. 1997;16:5819–26.
McDonnell JM, Fushman D, Cahill SM, Zhou W, Wolven A, Wilson CB, Nelle TD, Resh MD, Wills J, Cowburn D. Solution structure and dynamics of the bioactive retroviral M domain from Rous sarcoma virus. J Mol Biol. 1998;279:921–8.
Hatanaka H, Iourin O, Rao Z, Fry E, Kingsman A, Stuart DI. Structure of equine infectious anemia virus matrix protein. J Virol. 2002;76:1876–83.
Riffel N, Harlos K, Iourin O, Rao Z, Kingsman A, Stuart D, Fry E. Atomic resolution structure of Moloney murine leukemia virus matrix protein and its relationship to other retroviral matrix proteins. Structure. 2002;10:1627–36.
Saad JS, Ablan SD, Ghanam RH, Kim A, Andrews K, Nagashima K, Soheilian F, Freed EO, Summers MF. Structure of the myristylated human immunodeficiency virus type 2 matrix protein and the role of phosphatidylinositol-(4,5)-bisphosphate in membrane targeting. J Mol Biol. 2008;382:434–47.
Serrière J, Robert X, Perez M, Gouet P, Guillon C. Biophysical characterization and crystal structure of the Feline Immunodeficiency Virus p15 matrix protein. Retrovirology. 2013;10:64.
Tang C, Loeliger E, Luncsford P, Kinde I, Beckett D, Summers MF. Entropic switch regulates myristate exposure in the HIV-1 matrix protein. Proc Natl Acad Sci USA. 2004;101:517–22.
Brown L, Cox C, Baptiste J, Summers H, Button R, Bahlow K, Spurrier V, Kyser J, Luttge B, Kuo L, Freed E, Summers M. NMR Structure of the myristylated feline immunodeficiency virus matrix protein. Viruses. 2015;7:2210–29.
Prchal J, Srb P, Hunter E, Ruml T, Hrabal R. The structure of myristoylated Mason-Pfizer monkey virus matrix protein and the role of phosphatidylinositol-(4,5)-bisphosphate in its membrane binding. J Mol Biol. 2012;423:427–38.
Hill CP, Worthylake DK, Bancroft DP, Christensen AM, Sundquist WI. Crystal structures of the trimeric human immunodeficiency virus type 1 matrix protein: implications for membrane association and assembly. Proc Natl Acad Sci USA. 1996;93:3099–104.
Vlach J, Srb P, Prchal J, Grocký M, Lang J, Ruml T, Hrabal R. Nonmyristoylated matrix protein from the Mason-Pfizer monkey virus forms oligomers. J Mol Biol. 2009;390:967–80.
Fledderman EL, Fujii K, Ghanam RH, Waki K, Prevelige PE, Freed EO, Saad JS. Myristate exposure in the human immunodeficiency virus type 1 matrix protein is modulated by pH. Biochemistry. 2010;49:9551–62.
Martin DDO, Beauchamp E, Berthiaume LG. Post-translational myristoylation: fat matters in cellular life and death. Biochimie. 2011;93:18–31.
Doležal M, Zábranský A, Hrabal R, Ruml T, Pichová I, Rumlová M. One-step separation of myristoylated and nonmyristoylated retroviral matrix proteins. Protein Expr Purif. 2013;92:94–9.
Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372:774–97.
Zábranský A, Hoboth P, Hadravová R, Stokrová J, Sakalian M, Pichová I. The noncanonical Gag domains p8 and n are critical for assembly and release of mouse mammary tumor virus. J Virol. 2010;84:11555–9.
Zábranský A, Sakalian M, Pichová I. Localization of self-interacting domains within betaretrovirus Gag polyproteins. Virology. 2005;332:659–66.
Sfakianos JN, Hunter E. M-PMV capsid transport is mediated by Env/Gag interactions at the pericentriolar recycling endosome. Traffic. 2003;4:671–80.
Clark J, Grznarova P, Stansell E, Diehl W, Lipov J, Spearman P, Ruml T, Hunter E. A Mason-Pfizer Monkey virus Gag-GFP fusion vector allows visualization of capsid transport in live cells and demonstrates a role for microtubules. PLoS One. 2013;8:e83863.
Pereira LE, Clark J, Grznarova P, Wen X, LaCasse R, Ruml T, Spearman P, Hunter E. Direct evidence for intracellular anterograde co-transport of M-PMV Gag and Env on microtubules. Virology. 2014;449:109–19.
Sakalian M, Hunter E. Separate assembly and transport domains within the Gag precursor of Mason-Pfizer monkey virus. J Virol. 1999;73:8073–82.
Knejzlík Z, Smékalová Z, Ruml T, Sakalian M. Multimerization of the p12 domain is necessary for Mason-Pfizer monkey virus Gag assembly in vitro. Virology. 2007;365:260–70.
Sakalian M, Dittmer SS, Gandy AD, Rapp ND, Zábranský A, Hunter E. The Mason-Pfizer monkey virus internal scaffold domain enables in vitro assembly of human immunodeficiency virus type 1 Gag. J Virol. 2002;76:10811–20.
Sommerfelt MA, Rhee SS, Hunter E. Importance of p12 protein in Mason-Pfizer monkey virus assembly and infectivity. J Virol. 1992;66:7005–11.
Srb P, Vlach J, Prchal J, Grocký M, Ruml T, Lang J, Hrabal R. Oligomerization of a retroviral matrix protein is facilitated by backbone flexibility on nanosecond time scale. J Phys Chem B. 2011;115:2634–44.
Schuck P. Size-distribution analysis of macromolecules by sedimentation velocity ultracentrifugation and Lamm equation modeling. Biophys J. 2000;78:1606–19.
Mueller U, Darowski N, Fuchs MR, Förster R, Hellmig M, Paithankar KS, Pühringer S, Steffien M, Zocher G, Weiss MS. Facilities for macromolecular crystallography at the Helmholtz-Zentrum Berlin. J Synchrotron Radiat. 2012;19:442–9.
Kabsch W. XDS. Acta Crystallogr Sect D Biol Crystallogr 2010; 66:125–132.
Sheldrick GM. Experimental phasing with SHELXC/D/E: combining chain tracing with density modification. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 4):479–85.
Pape T, Schneider TR. HKL2MAP: a graphical user interface for macromolecular phasing with SHELX programs. J Appl Crystallogr. 2004;37:843–4.
Schneider TR, Sheldrick GM. Substructure solution with SHELXD. Acta Crystallogr D Biol Crystallogr. 2002;58:1772–9.
Cowtan K. The Buccaneer software for automated model building. 1. Tracing protein chains. Acta Crystallogr D Biol Crystallogr. 2006;62:1002–11.
Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr D Biol Crystallogr. 2004;60:2126–32.
Pannu NS, Murshudov GN, Dodson EJ, Read RJ. Incorporation of prior phase information strengthens maximum-likelihood structure refinement. Acta Crystallogr D Biol Crystallogr. 1998;54:1285–94.
Žáková L, Kletvíková E, Veverka V, Lepšík M, Watson CJ, Turkenburg JP, Jiráček J, Brzozowski AM. Structural integrity of the B24 site in human insulin is important for hormone functionality. J Biol Chem. 2013;288:10230–40.
Case DA, Babin V, Berryman JT, Betz RM, Cai Q, Cerutti DS, T.E. Cheatham I, Darden TA, Duke RE, Gohlke H, Goetz AW, Gusarov S, Homeyer N, Janowski P, Kaus J, Kolossváry I, Kovalenko A, Lee TS, LeGrand S, Luchko T, Luo R, Madej B, Merz KM, Paesani F, Roe DR, Roitberg A, Sagui C, Salomon-Ferrer R, Seabra G, Simmerling CL, et al.: AMBER 14. 2014.
Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C: Comparison of multiple amber force fields and development of improved protein backbone parameters. Proteins Struct Funct Genet 2006;712–725.
Bayly CCI, Cieplak P, Cornell WD, Kollman PA. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. J Phys… 1993; 97:10269–80.
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA. Development and testing of a general Amber force field. J Comput Chem. 2004;25:1157–74.
MD, AZ, and IP designed the study and wrote the manuscript. MD purified the MMTV MA proteins. AZ prepared the plasmid constructs and performed the tissue culture experiments. JD and JB crystallized and determined the structure of myristoylated MMTV MA protein. OV performed the analytical ultracentrifugation analysis. ML performed the molecular dynamics simulation. RH performed the TEM analysis. IP supervised the project. All authors read and approved the final manuscript.
This work was supported by research projects RVO 61388963 and RVO 68378050; by Charles University (UNCE 204025/2012, SVV 260079/2014); by Ministry of Education, Czech Republic (NPU projects LO1302 and LO1304); and by the Gilead Sciences and IOCB Research Center. ML acknowledges the financial support of the Czech Science Foundation (grant number P208/12/G016). Diffraction data were collected on beamline MX14.2 at BESSY, Berlin, Germany. The authors would like to thank Hillary Hoffman for proofreading of the manuscript.
The authors declare that they have no competing interests.
Michal Doležal and Aleš Zábranský contributed equally
Additional file 1: Figure S1. MS MALDI TOF/TOF analysis of myristoylated, myr(+), and nonmyristoylated, myr(−), MMTV MA. (A) The sample of myristoylated MA contained myr(+) MA (m/z = 13006). The sample of nonmyristoylated MA contained myr(−) MA (m/z = 12795) and MA with uncleaved initial methionine (m/z = 12927). (B) To approximately determine the minimal detectable amount of myr(−) MA in the sample of myr(+) MA, the sample of myr(+) MA was mixed with the sample of myr(−) in the ratio 1000:1 (2:0.002 mg/mL). This amount of myr(−) MA was detectable whereas the sample of myr(+) contained no detectable amount of myr(−) MA.
About this article
Cite this article
Doležal, M., Zábranský, A., Dostál, J. et al. Myristoylation drives dimerization of matrix protein from mouse mammary tumor virus. Retrovirology 13, 2 (2016). https://doi.org/10.1186/s12977-015-0235-8
- Matrix protein
- Molecular dynamics
- Mouse mammary tumor virus
- Myristoyl switch
- Analytical ultracentrifugation
- X-ray crystallography