If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
During the late phase of HIV-1 infection, viral Gag polyproteins are targeted to the plasma membrane (PM) for assembly. Gag localization at the PM is a prerequisite for the incorporation of the envelope protein (Env) into budding particles. Gag assembly and Env incorporation are mediated by the N-terminal myristoylated matrix (MA) domain of Gag. Nonconservative mutations in the trimer interface of MA (A45E, T70R, and L75G) were found to impair Env incorporation and infectivity, leading to the hypothesis that MA trimerization is an obligatory step for Env incorporation. Conversely, Env incorporation can be rescued by a compensatory mutation in the MA trimer interface (Q63R). The impact of these MA mutations on the structure and trimerization properties of MA is not known. In this study, we employed NMR spectroscopy, X-ray crystallography, and sedimentation techniques to characterize the structure and trimerization properties of HIV-1 MA A45E, Q63R, T70R, and L75G mutant proteins. NMR data revealed that these point mutations did not alter the overall structure and folding of MA but caused minor structural perturbations in the trimer interface. Analytical ultracentrifugation data indicated that mutations had a minimal effect on the MA monomer–trimer equilibrium. The high-resolution X-ray structure of the unmyristoylated MA Q63R protein revealed hydrogen bonding between the side chains of adjacent Arg-63 and Ser-67 on neighboring MA molecules, providing the first structural evidence for an additional intermolecular interaction in the trimer interface. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.
). Binding of HIV-1 Gag polyproteins to the inner leaflet of the PM is mediated by the matrix (MA) domain, which harbors a bipartite signal represented by an N-terminal myristoyl (myr) group and a highly basic region (HBR). It is demonstrated that HIV-1 Gag association with membranes is regulated by multiple factors including electrostatic and hydrophobic interactions, protein multimerization, cellular and viral RNA, and specific phospholipids such as phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) (
Rous sarcoma virus gag has No specific requirement for phosphatidylinositol-(4,5)-bisphosphate for plasma membrane association in vivo or for liposome interaction in vitro.
). On the other hand, sedimentation equilibrium data confirmed that the MA protein resides in monomer–trimer equilibrium. NMR–based structural studies have shown that the myr group of MA can adopt sequestered and exposed conformations (
). Extrusion of the myr group can be increased by factors that promote protein self-association, such as increasing protein concentration or inclusion of the capsid (CA) domain (
). In an attempt to characterize the trimeric form of MA, we recently engineered a stable recombinant protein construct by fusing a foldon domain (FD) of phage T4 fibritin to the MA C terminus (
Foldon, the natural trimerization domain of T4 fibritin, dissociates into a monomeric A-state form containing a stable beta-hairpin: Atomic details of trimer dissociation and local beta-hairpin stability from residual dipolar couplings.
). Hydrogen–deuterium exchange mass spectrometry (HDX-MS) data supported an MA–MA interface that is consistent with that observed in the crystal structure of the myr(–)MA trimer (
). A plethora of structural and biophysical studies have provided invaluable insights into retroviral MA proteins binding to phospholipids and membrane mimetics (
The structure of myristoylated Mason-Pfizer monkey virus matrix protein and the role of phosphatidylinositol-(4,5)-bisphosphate in its membrane binding.
). These studies identified key molecular determinants of MA–membrane interactions and advanced our knowledge of retroviral assembly.
Studies by Barklis and coworkers revealed that the HIV-1 MA and MACA proteins assemble as hexagonal cage lattices on PS/cholesterol membrane monolayers and as hexamers of trimers in the presence of PI(4,5)P2 (
). The hexamer of trimers assembly is increasingly used as a model to explain the mechanism by which the envelope (Env) protein is incorporated into virus particles (
). In its processed form, Env consists of two noncovalently associated subunits, a surface glycoprotein (gp120) and a transmembrane domain (gp41), which are products of proteolytic cleavage of the viral gp160 precursor protein (reviewed in [
]). There is strong evidence that the cytoplasmic tail of gp41 (gp41CT) plays a functional role in Env incorporation in physiologically relevant cell types (
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
). NMR studies of the HIV-1 gp41CT protein revealed that the N-terminal 45 residues lack secondary structure and are not associated with the membrane. The C-terminal 105 residues, however, form three membrane-bound amphipathic α-helices with distinctive structural features such as variable degree of membrane penetration, hydrophobic and basic surfaces, clusters of aromatic residues, and a network of cation–π interactions (
Genetic and biochemical studies suggested that the MA domain of Gag and the gp41CT play distinct, yet complementary, roles in Env incorporation into budding particles (
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
Virion incorporation of envelope glycoproteins with long but not short cytoplasmic tails is blocked by specific, single amino acid substitutions in the human immunodeficiency virus type 1 matrix.
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
Virion incorporation of envelope glycoproteins with long but not short cytoplasmic tails is blocked by specific, single amino acid substitutions in the human immunodeficiency virus type 1 matrix.
). [Note: The N-terminal Met, which is absent in the myristoylated protein, is designated as residue 1. In contrast, other studies considered the N-terminal Gly of the myristoylated protein as residue 1]. When mapped on the putative hexamer of trimer model of MA, these residues were found to point toward the hexamer centers (Fig. 1) (
). It was also shown that substitution of residue Gln63 with Arg suppressed Env incorporation defects of the L13E, E17K, L31E, V35E, or E99V MA mutations and of a gp41CT mutation that has the same phenotype (
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
). Freed and coworkers provided biochemical evidence that MA trimerization is an obligatory step in the assembly of infectious HIV-1 and demonstrated a correlation between loss of MA trimerization and loss of Env incorporation (
). It was suggested that the Q63R mutation may stabilize the trimer structure such that MA lattices, which form large hexamer holes, are favored over those that feature small hexamer holes (Fig. 1) (
). In another study, biochemical and pull-down experiments suggested that the Q63R mutation promoted interaction with gp41CT; this mutation did not alter the organization of MA on a membrane monolayer (
). Nonconservative mutations at positions Ala45, Ser72, and Leu75 located in the trimer interface (Fig. 1) also impaired Env incorporation and infectivity, and biochemical cross-linking experiments showed that these mutations reduced trimerization of MA (
). Although genetic and biochemical studies provided evidence for a correlation between Env incorporation and formation of MA trimers, the structural or biophysical evidence for such a correlation is lacking. In particular, it is not known how single amino acid mutations that impair or rescue Env incorporation impact the structure and oligomerization properties of MA.
Figure 1MA trimer model and residues implicated in Env incorporation.A, structure of the myr(–)MA trimer (PDB ID: 1HIW) showing residues implicated in Env incorporation (green, blue, and red spheres). B, schematic representation of the hexamer of trimer arrangement of MA based on studies by Alfadhli et al. (
). In this model, a ∼45 nm central aperture is formed by residues at the center of the hexamer (green). It is suggested that the gp41CT protein is accommodated in the central aperture. Mutations of residues in the trimer interface (red spheres) significantly disrupted Env incorporation. C, perturbations of the putative hexamer or trimer interface in the MA lattice are thought to create smaller central aperture (∼30 nm) that may cause a steric exclusion of gp41CT.
In this report, we employed NMR spectroscopy, X-ray crystallography, and analytical ultracentrifugation (AUC) techniques to characterize the structure and trimerization properties of HIV-1 MA mutants that have been shown to impair Env incorporation (L75G, A45E, and T70R) or suppress Env incorporation defects (Q63R). We show that these mutations had no adverse effect on the structure or folding of MA and had only a minimal effect on the monomer–trimer equilibrium. NMR data revealed that these amino acid substitutions caused chemical shift perturbations (CSPs) for residues located in the trimer interface. The X-ray structure of myr(–)MA Q63R protein revealed hydrogen bonding between the side chains of adjacent Arg63 and Ser67 on neihbouring MA molecules located in the trimer interface. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.
Results
Effect of point mutations on the structure of MA
Point mutations in the HIV-1 MA protein (A45E, Q63R, T70R, and L75G) were made by site-directed mutagenesis of the MA gene embedded in a coexpression vector harboring yeast N-myristoyltransferase (
). Point mutations had no detectable effect on protein expression or the efficiency of myristoylation, which was judged by the elution profile of the protein from the hydrophobic interaction column. Two-dimensional (2D) 1H-15N HSQC NMR spectrum provides a structure-sensitive signal for each amide group in a protein, and thus, it can be used to report on changes to protein structure caused by amino acid substitutions. As shown in Figures 2 and S1, the HSQC spectra for the four MA mutant proteins are very similar to the spectrum of the wild-type (WT) protein, indicating that amino acid substitutions had no adverse effect on the structure and folding of the protein.
Effect of amino acid substitutions on the myr switch
Previous NMR studies of the WT MA protein have shown that a subset of 1H and 15N resonances for amino acid residues 3–18, Val35, Trp36, Arg39, Gly49, Glu52, His89, and Gln90 shift progressively toward the corresponding frequencies observed for myr(–)MA upon increasing protein concentration. These chemical shifts were attributed to exposure of the myr group and a concomitant shift in the monomer–trimer equilibrium toward the trimeric species (
). Subsequent NMR studies revealed that point mutations in the N-terminal region of MA (e.g., V7R, L8A, or L8I) that impaired membrane targeting of Gag and inhibited virus assembly and release (
Characterization of human immunodeficiency virus type 1 matrix revertants: Effects on virus assembly, Gag processing, and Env incorporation into virions.
). Structural data revealed conformational changes that appeared to be responsible for stabilizing the myristate-sequestered species and inhibiting exposure (
). Herein, we assessed whether the A45E, Q63R, T70R and L75G mutations in MA had any effect on the positioning of the myr group and/or the concentration-dependent myr-exposure. We obtained 1H-15N HSQC spectra for WT and A45E, Q63R, T70R, and L75G mutant MA proteins at three protein concentrations (50, 150, and 450 μM) (Figs. 2 and S2). Like WT MA, for all mutants the 1H and 15N resonances for residues 3–18, Val35, Trp36, Arg39, Gly49, Glu52, His89, and Gln90 were sensitive to protein concentration and shifted progressively toward the corresponding resonances observed for myr(–)MA upon increasing protein concentration (Figs. 2 and S2). Of note, although these signals shift toward the corresponding resonances of the myr(–)MA protein, the chemical shifts for these residues in the fully myr-exposed form of MA may not be at identical positions. Therefore, it is not possible to determine the population of the monomer and trimer states from the NMR data. However, based on the monomer–trimer association constants described below, it is estimated that the trimer population is ∼10%, 35%, and 63% at 50, 150, and 450 μM, respectively. Altogether, these results indicate that none of the four amino acid substitutions had any detectable effect on the concentration-dependent myr switch.
Figure 2NMR spectra of WT and mutant MA at different protein concentrations. Overlay of 2D 1H-15N HSQC spectra for WT, myr(–), and mutant MA proteins collected at different concentrations (50 μM [black], 150 μM [red], and 450 μM [green]). Like WT MA, for all four mutant proteins a subset of 1H-15N resonances of MA shifted toward the corresponding signals of myr(–)MA (blue), indicating a shift in equilibrium toward the myr-exposed state.
Chemical shift comparison between the WT and mutant MA proteins
To compare the backbone 1H, 15N and 13C chemical shifts with those of the WT MA protein, triple-resonance experiments were collected for all mutants. We mapped the 1H and 15N chemical shift changes caused by amino acid substitutions in MA and assessed whether they induced structural and/or conformational changes in distant regions of the MA protein, especially within the trimer interface. Typically, only a few signals corresponding to amino acids in the vicinity of the mutation site exhibit major chemical shift changes in the 1H-15N HSQC spectrum. Chemical shift differences between the WT and mutant MA proteins are provided in Figure 3, and a cartoon representation of the X-ray structure of the WT MA protein onto which the locations of residues that exhibited significant amide 1H-15N chemical shift change are mapped (Fig. 4). As expected, the largest differences are observed in the vicinity of mutation site. Interestingly, other regions that are distant from the mutation site were also perturbed. For the A45E mutant, in addition to the CSPs observed for signals of residues 39–47, which are in the vicinity of the mutation site, resonances corresponding to Glu74 and Leu75 also exhibited significant CSPs (Fig. 3). Interestingly, these two residues are quite distant from the mutation site within the same molecule but are intermolecularly adjacent (∼4.5 Å) in the trimer structure of myr(–)MA (Fig. 4). This result suggests that substitution of Ala45 to Glu may have induced structural perturbations in the trimer interface. Next, we examined the Q63R mutant. Gln63 is located in α-helix IV, with its side chain projecting toward the loop comprised of residues 42–46 (∼4.1 Å to residue Val46) within the same molecule. As shown in Figure 3, in addition to the mutation site, two other regions exhibited significant CSPs. These comprise residues 42–47 and 74–76. Mapping out the CSPs on the trimer structure of myr(–)MA revealed that this mutation can cause structural perturbations for residues 42–46 within the same molecule (Fig. 4). Therefore, it is likely that the CSPs corresponding to residues 42–46 are caused by intramolecular structural adjustments or perturbations as a result of this substitution. However, whereas amino acid residues 74–75 are intramolecularly distant from Gln63, these residues are located in the trimer interface (Fig. 4). In this case, the effect of Q63R mutation could be direct on residues 74–76 or indirect through a local structural perturbation of residues 42–47, which are also close to residues 74–76 (Fig. 4).
Figure 3Chemical shift differences between the WT and mutant MA proteins. Normalized 1H-15N chemical shift differences are plotted versus residue number. Δδ > 0.05 ppm indicates significant chemical shift changes. NMR spectra were collected at 150 μM.
Figure 4Chemical shift mapping on the structure of the myr(–)MA trimer. Structural mapping of the chemical shift differences onto the X-ray structure of WT MA trimer (PDB ID: 1HIW). The location of the amino acid mutation is marked with blue spheres. Positions of residues whose amide resonances exhibited Δδ > 0.05 ppm are highlighted in red (intermolecular) and green (intramolecular).
Among the four mutants, the largest CSPs have been observed in the HSQC spectra of the MA T70R and L75G proteins (Fig. 3). For T70R, the most significant CSPs are for residues localized in the vicinity of the mutation site (amino acid residues 69–76). Structural analysis of the MA trimer revealed that the side chain of Thr70, which resides in the loop connecting α-helices IV and V, is sandwiched between the side chains of Gln59 and Gln63 of a neighboring molecule (Fig. 4). Interestingly, the amide 1H-15N resonances of these two residues did not exhibit significant chemical shift changes. On the other hand, mutation of Leu75 to Gly induced substantial CSPs that extended beyond the mutation site, including residues in the trimer interface (Figs. 3 and 4). Of note, residue Leu75, located in helix V, is not involved in intermolecular interactions with a neighboring molecule. However, within the same molecule, the side chain of Leu75 is relatively close to the side chain of Leu64 located in helix V (4 Å) and is spatially close to the loop made of residues 67–72. As discussed above, residues 67–72 are located in the trimer interface and are in close proximity to the loop made of residues 43–47. Therefore, it is conceivable that Leu75 contributes to the orientation of the loop.
Taken together, NMR data provided insights into structural perturbations within the trimer interface induced by the amino acid substitutions. However, we were unable to precisely determine these structural changes in the trimer interface by NMR methods because, like the WT MA protein (
), detection of intermolecular contacts (e.g., NOE) proved to be technically challenging. In summary, we provided structural evidence that substitution of amino acid residues in the trimer interface can directly or indirectly contribute to the (de)stabilization of the MA trimer form.
Crystal structure of the myr(–)MA Q63R protein
As discussed above, to explain the role of MA trimerization in Env incorporation, it was suggested that the Q63R mutation in MA may stabilize the trimer structure such that MA lattices that form large hexamer holes are favored over those that feature small hexamer holes (Fig. 1) (
). However, it is not clear how this substitution can stabilize the trimer form. We have shown that, similar to the WT MA protein, MA Q63R is in monomer–trimer equilibrium. As indicated by the NMR data, substitution of Gln63 to Arg did not induce significant CSPs of residues in the trimeric interface. To be able to assess the structural changes in the MA trimer, we determined the high-resolution structure of myr(–)MA Q63R protein by X-ray crystallography. Atomic coordinates and structure factors for two crystal forms have been deposited in the PDB (codes 7JXR and 7JXS; Table 1). Interestingly, the two crystal forms of myr(–)MA Q63R yielded structures that are essentially identical with slightly different parameters (Fig. S3 and Table 1). Each asymmetric unit of the crystal structure contains two trimers. The structure of the myr(–)MA Q63R trimer in both crystal forms is virtually identical to that of the WT myr(–)MA protein with slight orientation differences of monomers in the trimeric assembly (Fig. 5) (
). The structures in the two crystal forms revealed only minor conformational changes in the loops connecting helices I and II and helices V and VI (Fig. S3). Intriguingly, the guanidinium group of Arg63 is 2.5 Å (averaged across 12 examples in the asymmetric unit of the two crystal forms) from the hydroxyl group of Ser67 of an adjacent MA molecule, which forms new hydrogen bonding in the trimer interface (Fig. 5). This result demonstrates that substitution of Gln63 with Arg stabilizes the trimer structure of MA in the crystal lattice through H-bonding between the side chains of Arg63 and Ser67 located on two neighboring MA molecules.
Table 1X-ray diffraction and refinement statistics
Parameter
Crystal 1
Crystal 2
Wavelength
1.000
1.000
Resolution range
46.21–2.04 (2.113–2.04)
36.82–2.35 (2.434–2.35)
Space group
P 1 21 1
P 1 21 1
Unit cell
63.13 97.92 84.81 90 99.94 90
64.54 90.72 73.16 90 102.38 90
Total reflections
171742
100393
Unique reflections
61,490 (6078)
32,573 (3286)
Multiplicity
2.8
3.1
Completeness (%)
95.03 (94.19)
94.51 (95.36)
Mean I/sigma(I)
28.86 (2.50)
25.95 (2.38)
Wilson B-factor
42.84
49.51
R-merge
0.067 (0.408)
0.072 (0.431)
R-meas
0.083 (0.508)
0.088 (0.515)
R-pim
0.048 (0.298)
0.050 (0.277)
CC1/2
0.987 (0.726)
0.987 (0.807)
Reflections used in refinement
61,481 (6078)
32,566 (3285)
Reflections used for R-free
1995 (197)
1999 (202)
R-work
0.2032 (0.2525)
0.2170 (0.2481)
R-free
0.2340 (0.3166)
0.2560 (0.2723)
Number of nonhydrogen atoms
5507
5286
macromolecules
5141
5111
Ligands
149
108
Solvent
217
67
Protein residues
638
632
RMS(bonds)
0.008
0.008
RMS(angles)
1.2
1.07
Ramachandran favored (%)
100
98.06
Ramachandran allowed (%)
0.00
1.61
Ramachandran outliers (%)
0.00
0.32
Rotamer outliers (%)
1.07
3.22
Clashscore
5.99
8.92
Average B-factor
51.98
65.77
macromolecules
51.67
65.83
Ligands
59.99
65.32
Solvent
53.77
61.45
Number of TLS groups
40
29
The values in parentheses correspond to the highest resolution data shell.
Figure 5High-resolution crystal structure of the HIV-1 myr(–)MA Q63R protein.A, superposition of WT myr(–)MA trimer protein (PDB code 1HIW; blue) and the two trimers found in the asymmetric unit of the crystal structure of myr(–)MA Q63R (yellow and gray). The structure of the myr(–)MA Q63R mutant is essentially identical to that of the WT myr(–)MA protein. B, surface representation of the trimer arrangement of the myr(–)MA Q63R protein. Amino acid residues Arg63 and Ser67 are shown as sticks. C, cartoon representation of myr(–)MA Q63R mutant showing the H-bond formed between the side chains of Arg63 and Ser67. Residue Thr70 is also shown to highlight the close proximity to Arg63.
As mentioned earlier, previous sedimentation equilibrium (SE) studies have shown that, while myr(–)MA is monomeric at all tested protein concentrations and pH values (
). These results were explained by a (de)protonation process of the His89 imidazole ring. Deprotonation of His89 destabilizes the salt bridge formed between His89(Hδ2) and Glu12(COO-), leading to tight sequestration of the myr group and a shift in the equilibrium from trimer to monomer. To determine how the various MA mutations impacted the oligomerization properties of the MA protein in solution, we collected sedimentation velocity (SV) and SE data on the WT and mutant MA proteins at pH 5.5. Consistent with previous data, the SV profile for myr(–)MA exhibits a sharp sedimentation boundary indicating a monomeric state (Fig. 6). The MA protein, however, exhibits a broad sedimentation boundary with a larger S value consistent with a monomer–trimer equilibrium. The MA mutants exhibited slightly variable sedimentation boundaries with an overall sedimentation coefficient (c(s)) distribution similar to that of the WT MA protein. Interestingly, the sedimentation boundary for the L75G mutant is sharper than the WT and the other mutants and shifted to a smaller S value, suggesting that the monomer–trimer equilibrium is shifted toward the monomeric form (Fig. 6). Next, we collected SE data for the WT and mutant MA proteins at pH 5.5 (Fig. S4). Like the WT protein, SE data for all mutants best fit a monomer–trimer model. Interestingly, while the Ka values for the A45E, T70R, and Q63R mutant proteins are similar to the WT protein, the corresponding value for the L75G mutant is approximately twofold lower (Table 2), consistent with a lower trimer population. Taken together, the SV and SE data indicate that whereas none of the four mutations in MA significantly altered the monomer–trimer equilibrium in solution, the L75G mutation appears to shift the equilibrium toward the monomer form.
Figure 6Sedimentation velocity data. Sedimentation coefficient distributions, c(s), obtained from the sedimentation profiles for WT and mutant MA proteins at pH 5.5.
Point mutations of MA residues Leu13, Glu17, Leu31, Val35, and Glu99, located at the tip of the hexamer centers in the hexamer of trimer model (Fig. 1), were found to impair Env incorporation without affecting virus particle formation (
Virion incorporation of envelope glycoproteins with long but not short cytoplasmic tails is blocked by specific, single amino acid substitutions in the human immunodeficiency virus type 1 matrix.
). Freed and coworkers have provided biochemical evidence that MA trimerization is an obligatory step in the assembly of infectious HIV-1 and demonstrated a correlation between loss of MA trimerization and loss of Env incorporation (
). Nonconservative mutations at, or near, the trimer interface in the crystal structure (Fig. 1) inhibited MA trimerization and yielded particles with impaired Env incorporation and infectivity (
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
). The current hypothesis suggests that Q63R mutation may stabilize the trimer structure such that MA lattices, which form large hexamer holes, are favored over those that feature small hexamer holes (Fig. 1) (
). This hypothesis was further supported by the finding that the Q63R mutation promoted interaction with gp41CT without altering the organization of MA on a membrane layer (
). A correlation between MA trimerization and gp41CT binding was also suggested in a biochemical study involving MA mutants and MACA proteins, supplemented with inositol polyphosphates (
In this study, we have shown that: (i) mutations that either impaired (A45E, T70R, and L75G) or rescued (Q63R) Env incorporation had no adverse effects on the structures and folding of the MA protein. This result rules out the possibility that such mutations produce a nonfunctional protein that indirectly impacts its ability to form a trimer or participate in other functions involved in Env incorporation. (ii) None of the mutations had any detectable effect on the efficiency of myristoylation or concentration-dependent myr switch. NMR data obtained at different protein concentrations indicate that all four mutants display a concentration-dependent myr switch mechanism, similar to the WT MA protein. (iii) Mutations appear to induce minor structural perturbations in the trimer interface, which may directly or indirectly contribute to the (de)stabilization of the MA trimer form. As detected by NMR, the largest CSPs were observed for the L75G mutant. This result is not surprising given that Leu75 is located in an α-helical motif. Substitution with a glycine reside, which is sometimes known as a helix breaker, may have induced a larger structural change than the other three mutants. As mentioned earlier, the effect of Leu75 mutation on MA trimerization is probably indirect as this residue is not involved in intermolecular interactions with a neighboring MA molecule. Leu75 is in close proximity to Leu64, which is in juxtaposition to the loop made of amino acid residues 67–72 that interacts with amino acid residues 42–46 of a neighboring molecule. Therefore, it is conceivable that Leu75 contributes to the orientation of this loop to enable interaction with an adjacent molecule. (iv) A new hydrogen bond is formed between the side chains of adjacent Arg63 and Ser67 residues on neighbouring MA molecules in the center of the trimer interface of myr(–)MA Q63R, providing the first structural evidence for an additional intermolecular interaction in the trimer interface. (v) SV and SE data show that the A45E, T70R, and Q63R mutations had minimal effect on the monomer–trimer equilibrium, whereas the L75G mutation enhanced the population of the monomeric form by twofold.
As discussed above, the structures of monomeric MA and myr(–)MA are nearly identical (
). The MA protein, however, is in monomer–trimer equilibrium. To determine whether the MA–MA interface in solution resembles that in the X-ray structure of myr(–)MA (
). Hydrogen–deuterium exchange MS data supported a MA–MA interface that is consistent with that observed in the crystal structure of the myr(–)MA trimer (
). Of note, the region that exhibited a significant increase in protection from deuteration in the HDX-MS assay is the same region identified here as a putative trimer interface. Our data presented here demonstrate that neither of the three mutants that were proposed to inhibit MA trimerization (A45E, T70R, and L75G) completely shifted the equilibrium to a monomer, nor did the Q63R mutation, which rescues Env incorporation, form a stable trimer in solution even though a new H-bond is observed between the side chains of Arg63 and Ser67 from neighboring subunits (Fig. 5). This discrepancy between the crystal structure and solution data is still not well understood. It is possible that other forces such as the position of the myr group and the conformation of the N-terminal helix may also influence the oligomerization properties. In the absence of detailed structural data of the myristoylated MA trimer, it is not possible to determine how all factors contribute to stabilization or destabilization of the trimer form. Furthermore, it is worth mentioning that the mutations that impair or rescue Env incorporations were studied in the context of full-length Gag. Other factors that were found to influence Gag oligomerization include the CA and SP1 domains, NC-RNA interaction, and membrane binding ([
It is important to highlight the link between the structural findings and the previous biological findings. It is worth noting that previous studies indicated that residue Gln63 is not crucial for Env incorporation in the context of otherwise-WT MA and that Q63R is unique for its ability to fully rescue Env incorporation defects (
). It appears that the type of amino acid substitution at position Gln63 is key for its function in restoring Env incorporation. When introduced as a single mutation in MA, Q63R only minimally enhanced Env incorporation (
). However, the effect was more dramatic in restoring Env incorporation to a WT level when this mutation was introduced along with the Env-defective L13E mutant (
). By substituting Gln63 with E, G, K, L, N, or W, it was shown that all mutants replicated with WT kinetics in Jurkat cells and that none of the single mutations severely impaired virus release, infectivity, or Env incorporation. When the L13E/Q63[E/G/K/L/N/W] double mutants were subjected to the same analysis, none of the Gln63 mutants was able to fully rescue the virus replication, infectivity, and Env incorporation defects imposed by the L13E mutant (
). Interestingly, even the L13E/Q63K mutant exhibited a partial rescue, as infectivity in TZM-bl cells was comparable with that of L13E/Q63R, no rescue of Env incorporation was apparent. Taken together, Q63R appears to be unique for its ability to fully rescue Env incorporation defects. Our findings complemented the biological data by showing that the H-bond formed between the side chains of Arg63 and Ser67 is key for the stability of the MA trimer and hence to its function in rescuing Env incorporation (
). It is possible that slight changes to trimerization may be amplified in a biological context, due to the 2D constraint imposed by the membrane or the presence of other oligomerization promoting domains of Gag such as CA.
Role of Ser67 and Thr70 in Env incorporation in the context of WT, L13E, Q63R, and L13E/Q63R clones was also investigated (
). Whereas the S67A mutation had no effect on the phenotypes of the four viral clones, T70A blocked the ability of Q63R to rescue L13E infectivity and Env incorporation, although Q63R/T70A was as infectious as Q63R alone. The single T70A mutant was also impaired for infectivity and Env incorporation, suggesting that Thr70 may be involved in MA function. Another interesting result is the finding that S67R behaved like Q63R in its ability to rescue the defect imposed by L13E, indicating that an Arg at either residue 63 or 67 could rescue the L13E defect (
). Based on the structure of the myr(–)MA Q63R protein (Fig. 5), it is clear that Ser67 plays an important role in forming a H-bond with Arg63. It is expected that an Arg at position 67 is also capable of forming an H-bond with the side chain of Gln63 on a neighboring molecule, therefore stabilizing the MA trimer structure. The close proximity of Arg63 and Ser67 and their role in the stabilization of the MA trimer are supported by the finding that MA forms dimer and trimer in the glutaraldehyde cross-linking assay conducted with the double mutant (Q63K/S67K) in Jurkat and MT4 cell lines (
Combining previous studies with the data presented here, it appears that the correlation between trimerization of the MA domain of Gag and Env incorporation is not straightforward. Mutation of residues located in the trimer interface, on the tip of the hexamer ring, or even distant residues can impact Env incorporation. In a recent study, new compensatory mutations that rescue MA trimer interface mutants with severely impaired Env incorporation in MT-4 T-cells were identified (
). Viruses with MA L75G mutation acquired mutations G11S, V35I/F44I, V35I/E52K, V46I, and E52K. For viruses with MA L75E mutation, compensatory mutations D14N/V45I, Q28K, V35I/F44L, and E52K were identified. Of note, among all these secondary mutations only Phe44 and Val46 are located in the trimer interface. It was suggested that mutations that are distant from the trimer interface may indirectly influence the trimer interface by inducing potential structural or conformational changes (
). Structural studies to discern conformational changes induced by the compensatory mutations are warranted.
Our data support the hypothesis that mutations in MA that block Env incorporation may do so by disrupting an otherwise-stable MA lattice, rather than by disrupting a specific MA–gp41 interaction. The importance of the trimer interface in rescuing Env incorporation suggests that a stable MA trimer is key to trap Env into the Gag lattice. This model is also supported by recent high-resolution single-molecule tracking data (
), which demonstrated that Env trimers are confined to subviral regions of a budding Gag lattice, supporting a model where direct interactions and steric corralling between the gp41CT and MA lattice promote Env trapping. In summary, the new findings advance our understanding of the complementary, yet distinct, roles the MA domain of Gag and gp41CT play in Env incorporation.
Experimental procedures
Plasmid construction and protein expression
A plasmid encoding for HIV-1 MA (pNL4-3 strain) and yeast N-terminal myristoyl transferase was provided by Dr Michael Summers (Howard Hughes Medical Institute, University of Maryland). The A45E, T70R, Q63R and L75G MA mutant constructs were generated using a QuickChange Lightning site-directed mutagenesis kit (Agilent Technologies). Forward and reverse primers (Integrated DNA Technologies) extended 15 base pairs on either side of the mutation codon. Mutations were verified by plasmid sequencing at the Heflin Genomics Core at the University of Alabama at Birmingham. WT and mutant MA and myr(–)MA proteins were prepared as described (
NMR data were collected at 35 °C on a Bruker Avance II (700 MHz 1H) or Avance III (600 or 850 MHz 1H) spectrometers equipped with cryogenic triple-resonance probes, processed with NMRPIPE (
). The backbone resonances were assigned using standard triple resonance data (HNCA, HN(CO)CA, HN(CO)CACB, and HNCACB) collected at 35 °C on 100–300 μM samples in a buffer containing 50 mM sodium phosphates (pH 5.5), 50 mM NaCl, and 1 mM tris(2-carboxyethyl)phosphine hydrochloride (TCEP). The triple-resonance experiments were collected as nonuniformly sampled (NUS) sparse data (20% sampling density in indirect dimensions) according to schemes generated using hmsIST (
). Chemical shift perturbations were calculated as , where ΔδH and ΔδN are 1H and 15N chemical shift changes, respectively. Histograms for the chemical shift changes were generated using gnuplot software (http://www.gnuplot.info) (
The HIV-1 myr(–)MA Q63R protein used in the crystallization trials was at 20 mg/ml in a buffer containing 10 mM Tris.HCl (pH 8), 1 mM EDTA, 2 mM β-mercaptoethanol, and 50 mM NaCl. Diffraction quality crystals were obtained using hanging-drop vapor diffusion in a solution of 40% PEG (molecular weight 400), 0.1 M sodium acetate (pH 4), and 50 mM lithium sulfate. Crystals were cryo-cooled in the same conditions. X-ray diffraction data were collected at the Advance Photon Source, Southeast Regional Collaborative Access Team (SER-CAT) Beamline 22-ID. Raw intensity data were processed with the HKL2000 software package (
). Visualization of structures was performed using Pymol.
Analytical ultracentrifugation
Sedimentation velocity (SV) and sedimentation equilibrium (SE) experiments were performed on a Beckman XL-I Optima ultracentrifuge equipped with a four-hole An-60 rotor (Beckman Coulter). Cells were equipped with double-sector, charcoal-filled epon centerpieces with 12 mm path lengths and sapphire windows. Prior to AUC data collection, protein samples were run on a size-exclusion chromatography column (Superdex 75, 10/300GL, cytiva) in a buffer containing 50 mM sodium phosphates (pH 5.5), 100 mM NaCl, and 2 mM TCEP. Loading concentrations of MA samples ranged from 30 to 100 μM. For the SV experiments, rotor speed was set at 40,000 rpm. SE experiments were performed at 22,000, 26,000, and 30,000 rpm. All experiments were carried out at 20 °C. Scans were collected at a wavelength of 280 and 250 nm for the SV and SE experiments, respectively. Partial specific volumes (v-bar) and molar extinction coefficients were calculated by using the program SENDTERP, and buffer densities were measured pycnometrically. SV data were analyzed using SEDFIT (
). Sedimentation coefficients obtained from SV experiments were corrected to 20 °C and infinite dilution in water (s20,w). Equilibrium association constants were obtained by global fitting of SE scans collected at a single concentration and at three rotor speeds using a monomer–trimer equilibrium model. The molecular weight of the monomer (15.6 kDa) and the complex stoichiometry (n = 3) were treated as fixed parameters. Reference concentrations and baselines for each scan were floated parameters, and the Ka for the monomer–trimer equilibrium was fit as a global parameter.
Data availability
The atomic coordinates and structure factors (codes 7JXR and 7JXS) have been deposited in the Protein Data Bank (http://wwpdb.org/). The NMR chemical shift data for HIV-1 MA A45E, Q63R, L75G, and T70R are available from the Biological Magnetic Resonance Data Bank under BMRB accession numbers 50640, 50641, 50642, and 50643, respectively. The raw data described in the article can be shared upon request by directly contacting [email protected]
Conflict of interest
The authors declare that they have no conflicts of interest with the contents of this article.
Acknowledgments
We thank Dr Aaron Lucius (University of Alabama at Birmingham) for helping with the AUC experiments. We also thank the O'Neal Comprehensive Cancer Center at the University of Alabama at Birmingham (funded by grant P30 CA013148 from the National Institutes of Health, NIH) for supporting the High-Field NMR facility. X-ray data were collected at the SER-CAT 22-ID beamline at the Advanced Photon Source, Argonne National Laboratory. SER-CAT is supported by its member institutions and equipment grants (S10_RR25528 and S10_RR028976) from the NIH. Use of the Advanced Photon Source was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. W-31-109-Eng-38.
Author contributions
G. N. E., R. H. G., T. J. G., and J. S. S. designed the experiments. G. N. E. and R. H. G. expressed, purified, and characterized the proteins. G. N. E. and J. S. S. performed the NMR experiments and analyzed the results. T. J. G. collected and analyzed the X-ray data. G. N. E. and J. S. S. wrote the paper. G. N. E., R. H. G., T. J. G., and J. S. S. edited the paper.
Funding and additional information
This work was supported by grants 5R01GM117837 and 9R01AI150901 from the NIH to J. S. S. The High-Field NMR facility at the University of Alabama at Birmingham was established through the NIH (1S10RR026478) and is currently supported by the O’Neal comprehensive cancer center (NIH grant P30 CA013148). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Rous sarcoma virus gag has No specific requirement for phosphatidylinositol-(4,5)-bisphosphate for plasma membrane association in vivo or for liposome interaction in vitro.
Foldon, the natural trimerization domain of T4 fibritin, dissociates into a monomeric A-state form containing a stable beta-hairpin: Atomic details of trimer dissociation and local beta-hairpin stability from residual dipolar couplings.
The structure of myristoylated Mason-Pfizer monkey virus matrix protein and the role of phosphatidylinositol-(4,5)-bisphosphate in its membrane binding.
Elucidating the mechanism by which compensatory mutations rescue an HIV-1 matrix mutant defective for gag membrane targeting and envelope glycoprotein incorporation.
Virion incorporation of envelope glycoproteins with long but not short cytoplasmic tails is blocked by specific, single amino acid substitutions in the human immunodeficiency virus type 1 matrix.
Characterization of human immunodeficiency virus type 1 matrix revertants: Effects on virus assembly, Gag processing, and Env incorporation into virions.