article: Influenza virus penetrates cells by fusion of viral and endosomal membranes catalyzed by the viral hemagglutinin (HA). Structures of the initial and final states of the HA trimer define the fusion endpoints, but do not specify intermediates. We have characterized these transitions by analyzing low-pH-induced fusion kinetics of individual virions and validated the analysis by computer simulation. We detect initial engagement with the target membrane of fusion peptides from independently triggered HAs within the larger virus-target contact patch; fusion then requires engagement of three or four neighboring HA trimers. Effects of mutations in HA indicate that withdrawal of the fusion peptide from a pocket in the pre-fusion trimer is rate-limiting for both events, but the requirement for cooperative action of several HAs to bring the fusing membranes together leads to a long-lived intermediate state for single, extended HA trimers. This intermediate is thus a fundamental aspect of the fusion mechanism.
  Fusion of two lipid-bilayer membranes is a thermodynamically favorable process, but it crosses a high kinetic barrier as the two bilayers approach each other. Efficient fusion therefore requires a catalyst, a role served in living cells by a fusion protein or protein complex. The influenza virus hemagglutinin (HA) has become an important paradigm of a fusion catalyst, in part because of early crystallographic and mechanistic studies and in part because of continued concern about a virus that caused tens of millions of deaths during the 20th century. HA facilitates fusion by undergoing a large-scale conformational change, coupled to the two fusing membranes (virus and target). Our current picture of HA-mediated membrane fusion, illustrated in Figure 1A (Harrison, 2008), comes from HA structures in both pre- and post-fusion conformations and from inferences about transient intermediate states. HA is a homotrimer, synthesized as an inactive precursor, [HA0]3, and activated for fusion by proteolytic cleavage of each chain into HA1 and HA2, yielding [HA1-HA2]3. At the N-terminus of HA2 is a hydrophobic ‘fusion peptide' , which following cleavage inserts firmly into a pocket near the axis of the trimer (Chen et al. , 1998). Exposure to low pH, which during infection occurs in an endosome, causes the HA1 ‘head' to separate from the HA2 ‘stem' and enables a set of HA2 conformational transformations: (1) release of the fusion peptide from its pre-fusion pocket; (2) HA2 extension; (3) insertion of the fusion peptide into the target membrane; (4) fold-back of the extended HA2 intermediate (Figure 1A) (Skehel and Wiley, 2000). This last step brings together the fusion peptide and the C-terminal, transmembrane segment of each HA2, anchored respectively in the viral and target membranes, which thus approach each other, either as apposed protrusions or as a single, target-membrane protrusion (Kuzmin et al. , 2001; Lee, 2010). Fusion then ensues, initially as hemifusion (merger of only the proximal membrane leaflets) and then as formation of a continuous aqueous channel. 10. 7554/eLife. 00333. 003Figure 1. Single-virion analysis of fusion-promoting conformational change in influenza virus HA. (A) Hydrophobic fusion peptide (red asterisk) is initially inserted into a pocket near the trimer threefold axis. HA assumes an ‘open' conformation upon proton binding allowing fusion-peptide release and membrane insertion. The fold-back of the extended intermediate causes hemifusion. The known pre-fusion and post-fusion HA structures are colored, and the inferred structural transitions are showed in gray. (B) Left: Tile view from pH drop (t0) of a virion initially displaying directed motion (white arrow) followed by arrest (ta) then hemifusion (th). Red circle marks the final virion location. Arrowheads mark two virions that were arrested at pH drop and hemifused at or just before t40 s and t68 s respectively. Right: Fluorescence trace of the virion circled in (A) (red line), line fitting the timing of virion arrest at its target location (black line) and parameters t0 (green vertical line), ta (blue horizontal line) and th (orange horizontal line) and arrest to hemifusion, ta–h (dark orange horizontal line). (C) tlag (pHdrop–arrest) for all virions for which arrest values could be derived (left) and tlag (arrest–hemifusion) for all virions for which both arrest and hemifusion values could be derived (right) with gamma-distribution fit (black line). Data are pooled from three independent experiments. (D) Mean tlag (pHdrop–arrest) and tlag (arrest–hemifusion). Error bars represent the standard error of the mean. (E) and (F) N derived from fitting tlag (pHdrop-arrest) (E) and tlag (arrest–hemifusion) (F) with gamma probability density. (G) Rate constants derived from fitting tlag (arrest–hemifusion) and tlag (pHdrop–arrest) with gamma probability density and keeping N fixed (N = 3) for both processes. (B–G) X31-HA virions have X31 HA in otherwise Udorn genetic background. Data shown are from representative experiments performed on the same day (n = 50 to 150) unless indicated that multiple experiments were pooled. Error bars represent 95% confidence interval for gamma fit-derived values unless otherwise indicated. Please refer to Figure 1—figure supplement 1 for all histogram and gamma-distribution fit data plotted in (D–G). DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 00310. 7554/eLife. 00333. 004Figure 1—figure supplement 1. Histograms of tlag (pHdrop–arrest) (left) and tlag (arrest–hemifusion) (right) for X31-HA virions. Data are fitted with gamma-probability density where both N and k were fitted (black line) or N was fixed (N = 3) (orange line). Fit-derived values are indicated below each histogram in corresponding black and orange font colors. Parentheses indicate 95% confidence interval for each fit-derived value. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 004 We report here a series of experiments that probe the relationship between HA structural properties and kinetic intermediates in the fusion pathway. Floyd et al. (2008) devised a method to monitor in real time the fusion of individual influenza virus particles with planar bilayers. Their results led to the conclusion that fusion requires on average three HA trimers, each of which independently undergoes the same, rate-limiting rearrangement, but they left undetermined the relationship between this conclusion and the inferred intermediates in Figure 1A. In the experiments reported here, we correlate HA structure with observed variations in fusion kinetics by comparing rates for appropriate HA mutants. We conclude that irreversible engagement of fusion peptides from 3–4 neighboring [HA1-HA2]3 trimers, within a much larger virus-target-membrane interface, leads to subsequent rearrangements that rapidly and cooperatively induce membrane merger. Release of the fusion peptide from its pocket is rate-limiting for membrane engagement. A long-lived membrane-inserted extended intermediate is a fundamental aspect of the fusion mechanism.
  We recorded in real time a large number of individual influenza virions (approximately 1000 virions per field of view), labeled with a lipophilic fluorophore (R18), as they fused with a supported planar bilayer (see ‘Materials and methods' ). We followed hemifusion as a spike in individual virion fluorescence resulting from R18 fluorescence dequenching upon dilution into the target membrane (Figure 1B). We also observed a previously undetected fusion intermediate. Upon initiating the flow of buffer for pH exchange, most virions started to move under the hydrodynamic force while retaining contacts with the target bilayer. The particles arrested at various times after the pH drop, but invariably preceding hemifusion (Figure 1B; Video 1). The arrest was irreversible (Videos 2 and 3). A subset (approximately 20% at pH 5. 5 and 5. 65) of virions had already arrested at the onset of imaging or had done so before the pH drop. For lower final pH, a larger fraction had arrested by the time the pH drop was complete (Table 1). We determined lag times between pH drop and individual virion arrest (tlag (pHdrop–arrest) ) and between arrest and hemifusion (tlag (arrest–hemifusion) ) whenever both values could be extracted from the data (Figure 1B, C), for a range of proton concentrations between 2 and 100 µM (pH 5. 65–4) (Figure 1D). Mean values for tlag (pHdrop–arrest) and tlag (arrest–hemifusion) show the same pH dependence, with mean tlag (pHdrop–arrest) being about an order of magnitude shorter throughout the tested pH range. 10. 7554/eLife. 00333. 005Video 1. X31-HA WT virion hemifusion at pH 5. 5 from t0 to t230 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 00510. 7554/eLife. 00333. 006Video 2. Virion arrest is an irreversible intermediate of fusion. X31-HA WT virions were imaged. pH was dropped from 7. 4 to 5. 65 then brought back up to 7. 4. Buffer flow was kept constant except when it was stopped to allow for the source buffer exchange back to neutral (between approximately t45 s and t95 s). DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 00610. 7554/eLife. 00333. 007Video 3. Virion arrest is an irreversible intermediate of fusion. A different field of view of the same experimental lane used in Video 2 after the events imaged in Video 2. pH was dropped to pH 5. 65. There is a marked reduction in the mobile fraction in Video 3 relative to early times shown in Video 2. Furthermore, prearrested virions proceeded to hemifusion despite the intermediate reneutralization step. The videos are shown at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 00710. 7554/eLife. 00333. 008Table 1. Arrest and hemifusion statistics for X31-HA WT and arrest statistics for S4GHA2Udorn virionsDOI: http: //dx. doi. org/10. 7554/eLife. 00333. 008Virion arrest to hemifusion*pH drop to virion arrest†pH 5. 5X31-HA WT‡X31-HA WT‡S4GHA2Udorn§Number of virions3801215401Mean lag time (s) #105. 1 ± 2. 710. 3 ± 0. 29 ± 0. 4N¶3. 7 ± 0. 63. 4 ± 0. 23. 4 ± 0. 3Rate constant (s−1) ¶0. 032 ± 0. 0060. 33 ± 0. 020. 37 ± 0. 03Mobile at pH drop**81%76%Mobile that hemifused**75%65%Static that hemifused**82%81%pH 4. 5Number of virions155392385Mean lag time (s) #13. 3 ± 0. 61. 4 ± 0. 0041. 4 ± 0. 004Mobile at pH drop**46%38%*tlag (arrest–hemifusion) for all virions for which both arrest and hemifusion values could be derived. †tlag (pHdrop–arrest) for all mobile virions for which arrest values could be derived. ‡X31-HA WT data are pooled from three independent experiments. §S4GHA2Udorn data are pooled from two independent experiments. #Errors represent the standard error of the mean. ¶Errors represent the 95% confidence interval for the values derived from gamma-probability fits shown in Figure 1C (X31-HA virions) or Figure 4C (S4GHA2Udorn virions). **Percentages are derived from entire data sets. We can fit the distributions for both tlag (pHdrop–arrest) and tlag (arrest–hemifusion) with one describing a requirement for N independent events (either parallel or sequential; Feller, 1968), each with rate constant, k (Figure 1C and Table 1). For any given particle, the two lag times are uncorrelated, ruling out any mutual dependence on overall particle properties such as length (Figure 2). We have further controlled for virion length by using only shorter-virion fractions in our experiments (see ‘Materials and methods' ). From the distributions, we derive N = 3–4 both for the number of events required to arrest virions following pH drop (Figure 1E and Table 1) and for the number of events required for hemifusion of the arrested particle (Figure 1F and Table 1) (see ‘Materials and methods' ); the latter is similar to the value obtained previously from tlag (pHdrop–hemifusion) distributions (Floyd et al. , 2008). 10. 7554/eLife. 00333. 009Figure 2. Virion arrest and hemifusion lag times are uncorrelated. tlag (arrest–hemifusion) vs tlag (pHdrop–arrest) for X31-HA virions at pH 5. 5 (n = 380). Data are pooled from three independent experiments also shown in Figure 1C, Table 1 and Figure 4A, B (X31–HA). DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 009 A simple interpretation of virion arrest is that insertion of fusion peptides from a number of independent HA trimers (3 or 4) into the target membrane anchors the particle. Bulk experiments have shown that short incubations of virions with target membranes at pH 5 and 0°C lead to insertion of a small subset of fusion peptides and stable virus anchoring not associated with membrane fusion (Tsurudome et al. , 1992). We attribute the majority of pre-arrested events during high pH experiments (pH 5. 5 and above) to imperfections in the bilayer or defective virions; the observed increase in the immobile fraction at lower pH might result from genuine triggering and membrane insertion of HA fusion peptides upon proton binding during the pH transition (see ‘Materials and methods' ). To probe the molecular mechanism of arrest and its relationship to the mechanism of hemifusion (for which membrane insertion of the fusion peptide is clearly critical), we generated recombinant virus particles with site-directed mutations in HA. We used a set of plasmids derived from the A/Udorn/72 H3 influenza strain. The HA of Udorn is 97% identical in amino-acid sequence to that of X-31, the virus used in previous experiments (Floyd et al. , 2008) and also the source of HA in otherwise Udorn genetic background in the experiments in Figure 1. Nonetheless, recombinant virions with Udorn HA had shorter hemifusion lag times in the physiologically relevant pH regime (pH > 4. 5) than did those with X-31 HA (Figure 3A, Videos 1 and 4). Moreover, the Udorn-HA particles did not show directed motion under flow at any observable time point—that is, they had arrested by the time the low-pH transition was complete (see ‘Materials and methods' ). The Udorn hemifusion lag time (tlag (pHdrop–hemifusion) ) distribution between pH 5. 65 and 4 gave a pH-independent value of N close to 3 (Figure 3B), showing that differences between Udorn and X31 HAs do not affect the number of rate-limiting steps between virion arrest and hemifusion, but only their individual rate constants (Figure 3C). 10. 7554/eLife. 00333. 010Figure 3. Udorn virions have accelerated hemifusion kinetics. (A) Mean tlag (pHdrop–hemifusion) for Udorn and mean tlag (arrest–hemifusion) for X31-HA virions. Error bars represent the standard error of the mean. (B) N derived from fitting tlag (pHdrop–hemifusion) with gamma probability density. (C) Rate constants derived from fitting tlag (pHdrop-arrest) for Udorn and tlag (arrest–hemifusion) for X31-HA virions with gamma probability density and keeping each N fixed (N = 3). (B–C) Error bars represent 95% confidence interval for gamma fit-derived values. (A–C) Data shown are from representative experiments performed on the same day (n = 50 to 350). Please refer to Figure 3—figure supplement 1 for all histogram and gamma-distribution fit data plotted in (A–C). (D) Top left: Ribbon representation of X31-HA trimer (Weis et al. , 1990) showing positions of all residues that differ in Udorn-HA (arrowhead) including Gly4HA2 (asterisk). Top right: Cartoon of an HA monomer emphasizing Asp112HA2-fusion peptide hydrogen bond network and showing positions of residues along the HA2 chain. Bottom: Close-up of the Asp112HA2-fusion peptide network (region marked with an orange square on top left). DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01010. 7554/eLife. 00333. 011Figure 3—figure supplement 1. Histograms of tlag (arrest–hemifusion) for X31-HA virions (left) and of tlag (pHdrop–hemifusion) for Udorn WT virions (right). Data are fitted with gamma-probability density where both N and k were fitted (black line) or N was fixed (N = 3) (orange line). Fit-derived values are indicated below each histogram in corresponding black and orange font colors. Parentheses indicate 95% confidence interval for each fit-derived value. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 011 Comparison of the Udorn and X-31 HA amino-acid sequences suggests that one particular difference might account for the accelerated fusion kinetics of the former—a substitution of serine for glycine at position four in HA2. Udorn is the only strain in the database with this substitution, glycine being otherwise universally conserved (Nobusawa et al. , 1991; Cross et al. , 2009). In the pre-fusion conformation of HA, Gly4HA2 participates in a network of polar hydrogen-bond interactions with the carboxylate of Asp112 HA2, which is also strictly conserved (Figure 3D) (Wilson et al. , 1981; Weis et al. , 1990; Russell et al. , 2004). A serine substitution would weaken or interrupt this interaction, because the glycine has a backbone conformation not allowed for other residues. Previous studies have reported an elevated pH threshold of fusion for cell-surface expressed HA harboring substitutions at Gly4 or Asp112 (Gething et al. , 1986; Steinhauer et al. , 1995). Gly4HA2 might also favor a tight conformation of the membrane-inserted fusion-peptide helices (Lorieau et al. , 2010), but its conservation does not seem critical for free fusion-peptide insertion into membranes (Wharton et al. , 1988). We generated a series of HA mutants in the context of intact, recombinant virions, to test the consequences for fusion kinetics of the Gly4HA2–Asp112HA2 interaction and of the fusion-peptide sequestration in activated HA. These included G4SHA2X31 and S4GHA2Udorn, as well as D112AHA2 in both backgrounds. Figure 4A, B and Videos 1 and 4–8 summarize our analysis of hemifusion rates for these variants at pH 5. 5. The mutation G4SHA2X31 indeed accelerates hemifusion, and the reverse mutation S4GHA2Udorn retards it (Figure 4A). Moreover, like WT Udorn virions, G4SHA2X31 virions do not show flow-induced motion at pH drop, while 76% S4GHA2Udorn virions do (Videos 5 and 6, Table 1). Kinetic parameters for arrest, derived from the tlag (pHdrop–arrest) distribution for the S4GHA2Udorn particles are indistinguishable from those for X31-HA WT (Figure 4C and Table 1). These observations show that a single mutation, and thus probably the same molecular process, affects both the rate of arrest and the rate of hemifusion. 10. 7554/eLife. 00333. 012Figure 4. Fusion-peptide release from its pre-fusion pocket is a rate-limiting fusion-inducing molecular rearrangement in the physiologically relevant pH regime. (A) pH 5. 5, mean tlag (arrest–hemifusion) for X31-HA WT, S4GHA2Udorn and mean tlag (pHdrop–hemifusion) for Udorn WT, G4SHA2X31, D112AHA2Udorn and D112AHA2X31. (B) N derived from gamma-probability-density fits of the data analyzed also in (A). (C) Histogram of tlag (arrest–hemifusion) distribution for S4GHA2Udorn virions at pH 5. 5 with the gamma-probability-density fit (black line). (D) pH 4. 5, mean tlag (pHdrop–hemifusion) for indicated virions. (E) N derived from gamma-probability-density fits of the data analyzed also in (D). (A–E) Data shown are from pooled independent experiments (3, for Udorn and X31-HA WT virions, and 2, for Udorn and X31-HA mutant virions; n = 150 to 900). Error bars represent the standard error of the mean (A and D) or the 95% confidence interval for gamma fit-derived values (B and E). Please refer to Figure 4—figure supplement 1 for all histogram and gamma-distribution fit data plotted in (A, B, D and E). DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01210. 7554/eLife. 00333. 013Figure 4—figure supplement 1. Histograms of tlag (arrest–hemifusion) distributions for X31-HA WT, S4GHA2Udorn virions and of tlag (pHdrop–hemifusion) for Udorn WT, G4SHA2X31, D112AHA2Udorn and D112AHA2X31 at pH 5. 5 (top). Histograms of tlag (pHdrop–hemifusion) for all virions at pH 4. 5 (bottom). Data are fitted with gamma-probability density where both N and k were fitted (black line) or N was fixed (N = 3) (orange line). Fit-derived values are indicated below each histogram in corresponding black and orange font colors. Parentheses indicate 95% confidence interval for each fit-derived value. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01310. 7554/eLife. 00333. 014Video 4. Udorn WT virion hemifusion at pH 5. 5 from t0 to t90 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01410. 7554/eLife. 00333. 015Video 5. G4SHA2X31 virion hemifusion at pH 5. 5 from t0 to t90 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01510. 7554/eLife. 00333. 016Video 6. S4GHA2Udorn virion hemifusion at pH 5. 5 from t0 to t180 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01610. 7554/eLife. 00333. 017Video 7. D112AHA2X31 virion hemifusion at pH 5. 5 from t0 to t40 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01710. 7554/eLife. 00333. 018Video 8. D112AHA2Udorn virion hemifusion at pH 5. 5 from t0 to t35 s at 20× the actual rate. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 018 Data from D112AHA2 mutants confirm that the Gly4HA2–Asp112HA2 interaction has an important influence on hemifusion rate at pH 5. 5. D112AHA2 mutations accelerate hemifusion of both Udorn and X31-HA virions, but the effect on the latter is more pronounced (Figure 4A). We infer from these observations that the G4SHA2 mutation partly destabilizes contacts between the terminal fusion-peptide residues and Asp112HA2 and that the D112AHA2 change then eliminates this source of stability completely. None of the mutations change the values of N derived from the pH 5. 5 tlag (arrest–hemifusion) distributions (or for static virions, tlag (pHdrop–hemifusion) distributions) (Figure 4B), showing that these changes do not affect the number of steps. We conclude that for these viruses, release of the fusion peptide from its pre-fusion pocket is a rate-limiting step in the transition from arrest to hemifusion. The consistent changes in fusion kinetics produced by the mutations, including the good fit obtained from a probability distribution that assumes concomitant changes in the rates of all N events (see ‘Materials and methods' ), further support our proposal that N is related to the number of HA trimers that must independently expose their fusion peptides before a hemifusion-inducing transition can follow. The rate of arrest for either of the destabilized D112AHA2 variants was too rapid for us to detect flow-induced motion (Videos 7 and 8). We therefore could not determine N from arrest distributions for those viruses, but we could nonetheless infer that fusion-peptide release from its pocket is rate-limiting for virion arrest as well as for hemifusion. At pH 4. 5 and below, the mutants described above have indistinguishable mean values for tlag (pHdrop–hemifusion) (Figure 4D). Moreover, the lag times have become nearly pH independent, presumably because most of the relevant titrating groups (carboxylates) have become protonated. A subset (smaller than at pH 5. 5) of X31-HA WT and S4GHA2Udorn virions exhibit directed motion under the force of buffer flow at the pH drop, both with mean times of just under 1. 5 s (Table 1). Fitting the tlag (pHdrop–hemifusion) distributions with a gamma probability density yields unchanged values of N (between 3 and 6) for all the variants, confirming that what changes with pH is not the number of rate-limiting rearrangements leading to hemifusion but only their rate (Figure 4E). Protonation of Asp112 probably explains why mutations that affect stability of the fusion peptide in its pre-fusion pocket and influence the fusion rate at pH 5. 5 have little or no effect on the fusion rate at pH 4. 5 and below. The unchanged value of N and the observation that a fraction of the X31-HA WT and S4GHA2Udorn particles show directed motion after the pH drop suggest that even at pH 4. 5, some rate-limiting rearrangement precedes membrane insertion of the fusion peptide. The influence of mutations at positions 4 and 112 in HA2 on hemifusion at the physiologically relevant pH of 5. 5 demonstrates that the steps determining the rate of arrest (from the time of pH drop) and the rate of transition from arrest to hemifusion are the same, that is, both processes are limited by the rate of fusion-peptide exposure. Both also result from the action of several HA trimers each undergoing this rate-limiting rearrangement as indicated by the unchanged shape of their individual lag time distributions (unchanged value of N) under vastly different proton concentrations (Figure 1). The rate constants derived from gamma distribution fitting of the tlag (pHdrop–arrest) data are, however, about five to tenfold greater than those for tlag (arrest-hemifusion) (Figure 1), suggesting that the measured outcome in each case, arrest and hemifusion, is governed by different requirement (s). A simple interpretation of these observations is that virion arrest results from fusion-peptide release and irreversible target-membrane insertion from several HA trimers anywhere in the area of virus-target membrane contact and that hemifusion is limited by the same HA rearrangements, fusion-peptide release and membrane insertion, but requires a specific geometry for the participating HAs, that is, their proximity to each other, and is thus related to a different probability. To test the predictions of our model directly and to probe the relationship between the value of N derived from gamma-distribution fitting of hemifusion lag times and the actual number of participating HA trimers in the hemifusion reaction, we developed the following computer simulation (Figure 5). We defined a contact patch of various sizes between 50 and 300 HAs as a circular area of hexagonally arranged trimers. Our lower-limit estimate of a contact area that can accommodate up to about 50 HA trimers, shown in Figure 6, is based on dimensions of a spherical virion with 55 nm in membrane-to-membrane diameter (Calder et al. , 2010). Virions used in our current experiments are on average approximately 130 nm long and 55 nm wide (Ivanovic et al. , 2012; see ‘Materials and methods' ), and thus the actual contact area can on average accommodate up to about 120 HA trimers. For each hypothetical virion we obtained lag times for individual HA triggering and membrane insertion events by random selection from an exponentially decaying function with a fixed rate, kHApre-fusion–HAextended. We defined tlag (pHdrop–arrest) as the time when the Nath trimer within the contact area extended, and tlag (arrest–hemifusion) as the time from virion arrest to the time when Nhth neighbor of the previously extended trimers extended. Figure 5A shows definitions of HA neighbors for Nh values between 2 and 5 and accompanying simulation-derived arrest and hemifusion lag time distributions that assume Na = 3 and a contact patch of 121 HA trimers. 10. 7554/eLife. 00333. 019Figure 5. Computer simulation of the proposed model for HA-mediated membrane fusion. (A) Left: A contact patch containing 121 HA trimers and examples of individual virions and their contact-area HAs that have extended and inserted their fusion peptides into the target membrane by the time of hemifusion, for simulated Nh values of 2–5 and assuming either a contact patch consisting solely of active HAs (ii–v) or containing 20% inactive HAs randomly distributed throughout the contact area (vi). The neighboring HA clusters inducing hemifusion in each case are circled in green. Right: Na = 3 and patch size = 121 containing either 0% or 20% inactive HAs as indicated in the corresponding diagrams; simulation-derived tlag (pHdrop–arrest) distributions (i) and tlag (arrest–hemifusion) distributions for Nh = 2–5 (ii–vi) along with gamma-distribution fits (black curves) and fit-derived values for N and k. Parentheses indicate 95% confidence interval for each fit-derived value. We collected data for 500 virions in each example. (B) (i) Na = 3; effect of patch size (PS) and Nh either on gamma-fit-derived values for N from hemifusion lag-time distributions (top) or on ratios of hemifusion and virion arrest mean lag times (bottom). (ii) Patch size = 121; effect of Na and Nh either on gamma-fit-derived values for N from hemifusion lag-time distributions (top) or on ratios of hemifusion and virion arrest mean lag times (bottom). (iii) Na = 3, Nh = 3, patch size = 121; gamma-fit-derived N values (top) and mean lag times (bottom) for arrest and hemifusion at a range of inactive HA concentrations. Data shown are averages, and error bars, when included, represent standard deviation for five independent simulations each performed for 500 virions. For MATLAB script and accompanying functions please refer to Source code 1. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 01910. 7554/eLife. 00333. 020Figure 6. Stages of HA-mediated membrane fusion. Virion dimensions in the diagram are based on a spherical particle, 55 nm in membrane-to-membrane diameter (Calder et al. , 2010). There are approximately 50 trimers in the contact area between the virus and the target cell, shown in green (HA1) and gray (HA2); the trimers outside of the contact area are shown in white. The contact area illustrated represents a lower limit, and the contact area for the somewhat elongated particles in our current experiments can likely accommodate up to approximately 120 HA trimers (see ‘Results' and ‘Materials and methods' ). For clarity we present a virion surface containing only HA, but the conclusions we draw are independent of specific HA arrangement and potential irregularities in this region of contact (see Figure 5). We have also not included the target membrane in the perspective views (top), nor have we depicted the sialic-acid receptors in the cross-section views (bottom). The t0, ta and th labels are included to correlate the molecular transitions depicted in the model with the events defined in Figure 1B. (A) Virion engages a target membrane by HA-receptor interactions but can still move freely relative to the target-membrane area of contact. (B) Fusion peptides of several independent HA trimers have inserted into the target-membrane surface, immobilizing the virion on the target-cell surface. (C) Engagement of fusion peptides from three neighboring HA trimers allows subsequent rearrangements leading to hemifusion. HA1s are left out for clarity from the bottom panel, which shows three neighboring trimers about to fold back. DOI: http: //dx. doi. org/10. 7554/eLife. 00333. 020 Simulation-derived distributions agree remarkably well with our experimentally obtained data and predict Na = 3–4, Nh = 3–4 and are most consistent with our independently predicted contact patch area of approximately 120 HAs (Figure 5B, i and ii). The gamma-distribution-derived value of N for virion arrest closely approximates the number of trimers leading to arrest, Na, but the apparent rate constant is substantially higher than the rate constant governing HA extension and membrane insertion (Figure 5A, i). The gamma-distribution-derived value of N for hemifusion is related to the number of neighboring HA trimers participating in hemifusion, but it only approximates this value for smaller contact patch sizes (close to 50 HA trimers) and is greater than the true Nh for larger patch sizes (Figure 5B, i). The value of Nh = 2 yields a distribution that cannot be approximated by a gamma distribution for any analyzed patch size and value of Na (see Figure 5A, ii, for one such example). Likewise, the value of Nh = 5 yields a gamma distribution value of N > 5 for all analyzed patch sizes and N close to 7 for the predicted patch size of approximately 120 HA trimers (Figure 5B, i), inconsistent with our experimentally derived values. Our experimental data—the gamma-distribution-derived N and the ratio of mean hemifusion and arrest times—are consistent with Nh = 3–4 for patch sizes between 50 and 150 HA trimers. The apparent N-value-range between 3 and 6 (as observed for various mutants in Figure 4) corresponds to the actual Nh range between 3 and 4 and a patch size of about 120. We have also addressed the dependence of the predicted data on any irregularities that might be present within the area of virus-target membrane contact, such as divergence from ideal hexagonal arrangement of HA trimers, presence of inactive HAs or intercalation of the viral neuraminidase (NA). We found that the key parameters—the gamma-distribution-derived N value and the ratio of mean hemifusion and arrest lag times—are independent of the percentage of active HAs within the contact patch, down to a lattice-point occupancy of 40% (Figure 5B, iii). Further reduction in active HA density to 20% yielded a sudden drop in the apparent N derived from gamma-distribution fitting of simulation-derived tlag (arrest–hemifusion) data, consistent with previously published experimental evidence at pH 3, a proton concentration suggested to induce substantial HA denaturation (Floyd et al. , 2008).
  Many early studies of HA-mediated fusion kinetics relied on fusion of HA-expressing cells with red blood cells (Spruce et al. , 1989; Zimmerberg et al. , 1994) or planar bilayers (Melikyan et al. , 1993a, 1993b). In the latter case, measurements of conductance and capacitance showed reversible opening and closing of flickering pores, following initial acidification and commitment to secure pore opening only after a number of seconds. Good evidence that several HAs participate in pore formation came from red-blood cell fusion experiments in which the surface concentration of HA in the expressing cells varied over a roughly thirteenfold range (Danieli et al. , 1996). Hill-plot treatment of the fusion times led to an estimate of 3–4 HAs per fusion event. A limitation of any cell–cell fusion experiment is the large contact zone and hence some ambiguity concerning the number of independent, local fusion events that occur in any one observation. Single-virion fusion experiments confirmed the earlier estimate of 3–4 HAs and showed that lipid-mixing precedes content mixing, as expected for an obligatory hemifusion intermediate (Floyd et al. , 2008). In that work, the single, rate-limiting step for content mixing had a rate constant of about 0. 05/s. Calculations based on continuum models and molecular simulations have suggested various structures for hemifusion intermediates. The pathway with the lowest energy barrier appears to be through a so-called ‘hemifusion stalk' , derived from the apposed leaflets (Kuzmin et al. , 2001). Some proposals for the subsequent transition to a fusion pore postulate expansion of the stalk into a ‘hemifusion diaphragm' , in which the distal leaflets form a continuous bilayer between the two membrane-bounded compartments; others suggest that a stalk can open directly into a pore. Hemifusion diaphragms can be seen by electron microscopy (e. g. , Hernandez et al. , 2012); evidence for stalk-like structures comes from x-ray analysis of crystalline lipid phases (Yang and Huang, 2002). Both calculations and observations suggest that hemifusion diaphragms are stable, long-lived intermediates and that their formation might be a dead-end side reaction. Recent work on calcium-triggered SNARE/synaptotagmin/complexin-mediated fusion of lipid vesicles on a physiological time scale shows that the reaction proceeds too quickly (<100 ms) to detect even a hemifusion intermediate and that fusion from a hemifusion diaphragm is slow (tens of seconds) and inefficient (Diao et al. , 2012). In influenza-virus fusion, tight packing of proteins on the viral surface probably makes a transition from stalk or other point-contact intermediate to diaphragm even less likely than in SNARE-mediated fusion, in which protein coverage of the fusing membranes is sparser. The experimental results we describe above apply the previously developed single-virion experimental design (Floyd et al. , 2008) to virus particles with site-directed mutations in the HA protein, enabling us to relate the kinetics of the lipid-mixing step to molecular events in the HA trimer inferred from the pre-fusion and post-fusion structures. The experiments report on conformational rearrangements in the catalyst (HA); their interpretation does not depend on the details of membrane organization in the substrates (the hemifusing membranes). They provide strong evidence for a long-lived, extended HA intermediate and for a rapid and cooperative transition to lipid mixing when a critical number of these intermediates has accumulated within a contact neighborhood. Thus, they begin to elucidate features of the ‘gray zone' in Figure 1A. Figure 6 illustrates the model derived from our current experiments (Figures 1–4), as supported by our computer simulation (Figure 5). Binding to sialic-acid-containing glycans is weak, and local forces, such as flow, can displace the particle (either by ‘rolling' , with dissociation at the rear edge accompanied by new attachments at the front, or by translation across the membrane surface together with any associated lipids). The contact region between a virus particle and the target membrane is relatively large (probably 50–150 HA trimers, depending on the particle morphology: see comments above and Figure 6 caption). Any fusion peptide released within this area of contact and inserted into the target membrane will contribute to resisting buffer flow; the probability of arrest is then the probability for fusion-peptide exposure and insertion from 3–4 trimers within the area of membrane contact. Because hemifusion requires that fusion peptides of 3–4 adjacent trimers insert into the target membrane, the apparent rate constant for hemifusion relates to a different probability—that of inserting fusion peptides from several neighboring HAs. The data suggest that HA trimers independently undergo the first part of their low-pH-induced conformational change and pause as extended intermediates with their fusion peptides inserted into the target membrane. Capture by the target membrane restrains a trimer from collapsing to induce hemifusion until 2–3 neighbors help generate sufficient pull to deform the membrane. This interpretation also argues that when these conditions have been met, the final HA rearrangements are fast and cooperative. Cooperativity of the final HA conformational rearrangements can derive either from contacts between neighboring HAs or simply from the linkage that is present among all participating HAs by virtue of their common insertion in the two hemifusing membranes. We favor the latter explanation. There is no evidence for defined, lateral contacts between HA trimers, and any such contacts, were they present, could not persist during the substantial changes in HA conformation and orientation that accompany the progression to hemifusion. The accumulation of N extended HA2 neighbors, each fluctuating toward a more stable, fully collapsed conformation while pulling against the elastic restoring force of the two membranes, will lead to an abrupt and cooperative event, when the critical Nth molecule joins the cluster (rather like the abrupt consequence of adding a new member to the team on one side of a previously balanced tug-of-war). Evidence for the rapidity of the final HA rearrangements comes from the results of our simulation and gamma distribution analysis and the sensitivity of the derived rate constant to the mutations we have introduced. Any additional steps, not directly sensitive to those mutations, must be faster than any one of the steps of fusion-peptide engagement with the target membrane. Engagement of the fusion peptide with a target membrane requires formation of an extended intermediate to project it beyond the outer margin of adjacent HA1 heads (Figure 6B). For an intermediate to form that resembles the inner core of the ultimate, low-pH triggered structure (as drawn for the extended structure in Figure 1A), the HA1 heads must dissociate from their contact with an HA2 loop (HA2 residues 58–74: see Figure 3D) between the short first α-helix (HA2, residues 38–57) and the long central α-helix (HA2, residues 75–127) (Godley et al. , 1992), and that loop must in turn reverse direction and become itself a helix. Extension of the central trimeric coiled-coil probably drives the latter process, as the amino-acid sequence of the loop strongly favors such a conformation (Carr and Kim, 1993). Further necessary structural changes, in addition to release of the fusion peptide from its pocket, are dissociation of two strands (HA2, residues 21–38) from the edge of the small β-sheet near the base of the HA trimer and displacement outward of this segment of HA2 along with the short first helix, which could unfold and refold as it adds to the extending, central helical coiled-coil. The structure thus pictured would project the fusion peptide no more than about 25 Å beyond the palisade of un-triggered HAs (Figures 1A and 6), even if the HA2 segment that will fold back against the central core remains roughly in its pre-fusion state (again, as drawn for the extended structure in Figure 1A). In short, the intermediate that mediates arrest and hemifusion must indeed be almost completely extended, at least for part of its lifetime. The structures of the initial and final states do not dictate a unique ordering of the events just listed, but some must clearly precede others. Our results show that breaking the interactions that restrain the fusion peptide is rate-limiting both for virion arrest and for hemifusion at the pH of endosomes from which influenza virus penetrates (approximately 5 to 5. 5); this interpretation is consistent with the dominance hierarchy of mutations that influence the threshold pH for fusion (Steinhauer et al. , 1996). Mutations that affect the pH dependence of fusion map to a number of sites distributed across the HA2 trimer interface (Daniels et al. , 1985). Changes in the stability of this interface will influence fluctuations from the equilibrium structure and hence the ‘window of opportunity' for the fusion peptide to withdraw productively from its pocket. Emergence of the fusion peptide can in principle precede the other events, which then must follow rapidly (or the peptide will reinsert); perhaps more frequently, extension of the central coiled-coil by residues 20–74 may exert a ‘tug' to extract the fusion peptide. The fusion-peptide contacts, conserved in all field-isolated influenza HA virions, form upon cleavage of HA0 (Chen et al. , 1998; Russell et al. , 2004). Their stability is critical for viral transmission, not only in the passage from one host or cell to another, but also during entry, as fusion will be relatively unproductive if it occurs before the pH of the surrounding compartment has dropped enough to allow M2-mediated acidification of the particle interior (Ivanovic et al. , 2012), which induces dissociation of the eight ribonucleoproteins (RNPs) of the viral genome from the matrix protein, M1 (Martin and Helenius, 1991). Fusion too rapid for adequate acidification would release an inactive M1-RNP complex. Evidence for extended intermediates for other viral fusion proteins is largely indirect. The most extensive data are for HIV gp41, in which binding of peptides such as C34 (Chan et al. , 1998) and T-20 (Wild et al. , 1994), mimics for the outer layer of the final, hairpin structure, indicates that the inner-core coiled-coil must be present before the outer layer folds back against it. Outer-layer peptides from paramyxovirus F proteins similarly inhibit fusion (Lambert et al. , 1996). Electron microscopy of human parainfluenza virus 5 (PIV5) fusing with target vesicles shows a distance between viral and target membranes that corresponds to the gap expected for the putative extended intermediate of F, which (unlike the HA intermediate) should project much farther from the viral membrane than does the pre-fusion form (Kim et al. , 2011). The relative kinetics of HA-mediated hemifusion and arrest show that because fusion requires several neighboring HA trimers, the lifetime of an isolated extended intermediate is substantial. Release of the fusion peptide is rate limiting for generating an extended HA intermediate, but subsequent collapse is rapid only when three or four neighboring extended intermediates have appeared, and the mean lifetime of the ‘pioneer' extended HA in any patch can be as long as a minute or more, depending on pH. Thus, a relatively long-lived extended intermediate could prove to be a clinically useful target for an influenza-virus entry inhibitor: the inherent lower limit to the stability of the pre-fusion HA combined with the requirement for the coordinated action of multiple, independently-triggered HAs to induce fusion, curbs the ability of a virus to escape such treatment by mutations that accelerate fusion-inducing rearrangements. The only clinically approved inhibitor of virus fusion targets the extended intermediate of the HIV fusion protein, gp41 (Wild et al. , 1994; Kilby and Enron, 2003); analogous inhibitors have been designed for paramyxoviruses (Lambert et al. , 1996), influenza virus (Lee et al. , 2011) and flaviviruses (Schmidt et al. , 2010a, 2010b). For HIV fusion, the results of timing-of-inhibitor-addition experiments (Gallo et al. , 2001; Miyauchi et al. , 2009), the low surface density of the HIV-envelope protein (Chertova et al. , 2002; Zhu et al. , 2003), and the indication, from experiments with virions bearing mixtures of active and inactive spikes, that one HIV-envelope trimer may be sufficient to mediate fusion (Yang et al. , 2005), have led to the inference that a single, trimeric gp41 mediates fusion through a very slow fold-back transition. The data presented here for influenza HA show that there can be a long-lived species, analogous to the extended form of gp41, even when the fold-back transition itself is fast. We suggest that the long delay between HIV binding and membrane fusion may reflect a sequence of events similar to those in our model for influenza virus fusion. In other words, we propose that HIV fusion may require recruitment of (at least) a second envelope trimer, following attachment of the first. The delay time would then be a combination of the time needed to recruit one or more additional envelope proteins and the time needed to activate them with CD4 and a coreceptor. Extensions of the approaches taken here to HIV and other viruses will be necessary to probe the true similarities and differences among them.
  We based the protocol for recombinant influenza rescue on that developed by Neumann et al. (1999) and later adapted for A/Udorn/72 by Takeda et al. (2003) and which we described previously (Ivanovic et al. , 2012). All virions in the current study contained either Udorn or X31 HA, and the remaining genome segments were derived from A/Udorn/72. Virus plaque purification, passaging and purification of the third passage virus were performed as previously described (Ivanovic et al. , 2012) with the following modifications. We passaged both D112AHA2 mutants at reduced temperature (34°C instead of 37°C) and increased multiplicity of infection (0. 05 PFU/ml instead of 0. 001 PFU/ml). All final (third passage) infections were performed at 34°C. For the current experiments we used only the fractions enriched in ‘short' virions (average membrane-to-membrane length approximately 130 nm) (we previously published details about fractionating virions by size on sucrose gradients; Ivanovic et al. , 2012). We introduced specific mutations in Udorn or X31 HAs by site-directed mutagenesis with QuickChange (Stratagene, La Jolla, CA) based on manufacturer' s recommendations. We prepared two independent clones (derived from different plaques after initial virion rescue) for each WT and mutant viruses. We show data for only a single set of clones, but we reached the same basic conclusions with the second set: that G4SHA2X31 and S4GHA2Udorn mutations reverse the arrest and hemifusion rate properties and that D112AHA2 mutations accelerate hemifusion in both Udorn and X31-HA backgrounds, but the effect on X31-HA virions is more pronounced. We verified the presence of engineered mutations and the absence of spurious mutations in the open reading frame of the HA segment RNA of the purified viral products as follows. RNA was extracted using RNeasy mini kit (Qiagen, Valencia, CA), full-length HA segment RNA was reverse transcribed and amplified with universal primers specific for HA genome segment ends (Hoffmann et al. , 2001) using One-step RT-PCR kit (Qiagen) and sequenced using the same set of primers. The G4SHA2 mutation in Udorn HA2 (GenBank: ABD79032. 1), present also in the DNA used to derive Udorn HA by reverse genetics, is probably an anomaly of laboratory propagation at some point in the history of deriving the strain. Indeed, G4SHA2 was one of two mutations that arose in one isolate after several passages of a recombinant virus in mammalian cells (Lin et al. , 1997). Completeness of HA cleavage/activation in the purified virus preps was verified by Western Blotting using anti-HA-tag mouse monoclonal antibody (Cell Signaling, Danvers, MA), anti-mouse HRP secondary antibody and ECL Plus Western Blotting Detection System (GE Healthcare, Little Chalfont, UK). We used MATLAB (MathWorks, Natick, MA) for all data analysis. We developed a computer simulation of our proposed model using MATLAB (MathWorks). We defined a hexagonal lattice of HA trimers, in which each HA has exactly six neighbors each at a distance a. We then defined a circle of diameter, D=2×nπ×a, where n is an approximate number of HA trimers in the contact patch, and then identified the positions and the total number of HAs within the circle, nactual. For simulations involving a reduced fraction of active HAs, a given fraction of nactual HAs at random positions within a contact patch were flagged as inactive. We assumed that for each HA trimer, its conversion to an extended, membrane-inserted form is a single-step process, governed by a rate constant, kHApre-fusion-HAextended. We determined an arbitrary value for the extension rate constant of 0. 0025/s and kept this value unchanged in all our simulations. We first obtained lag times for extension of each HA trimer in a given patch by random selection from an exponentially decaying function with rate kHApre-fusion–HAextended. We sorted lag times for all HAs in ascending order. tlag (pHdrop–arrest) was defined as the lag time of Nath active HA trimer in a contact patch (where Na is a predetermined number of HAs that results in virion arrest). tlag (arrest–hemifusion) was defined as the lag time of that HA trimer that extended next to a defined number (Nh − 1) of active neighbors with shorter lag times minus the tlag (pHdrop–arrest). Nh is a value between 2 and 5 and is a predetermined number of extended active neighbors that results in hemifusion. Definitions of neighbors for Nh values between 2 and 5 that were used in our simulations are outlined in Figure 5A. The process was repeated for 500 virions in a typical simulation experiment. Simulation data for tlag (pHdrop–arrest) and tlag (arrest–hemifusion) were fitted as described under ‘Data analysis' for experiment-derived data.
summary: Influenza is caused by viruses that infect birds and mammals. These viruses enter cells when two lipid bilayers—one surrounding the virus, the other enclosing the cellular compartment into which the virus has been engulfed—merge to form a single unified membrane. This process, known as membrane fusion, allows the RNA of the virus to gain access to the host cell' s molecular machinery, which it commandeers to produce multiple copies of itself and to direct the assembly of new virus particles. The process of membrane fusion generally includes an intermediate hemifused state in which only the adjacent monolayers from each bilayer have merged. In addition to its role in virology, membrane fusion is critical for many other biological processes, including exocytosis, protein trafficking and the fertilization of eggs by sperm. Efficient membrane fusion requires a catalyst, and a glycoprotein known as the influenza hemagglutinin performs this role for the influenza virus. The hemagglutinin is found on the surface of the virus, and a typical influenza virus particle can have a few hundred such molecules on its surface. When an influenza virus particle binds to the surface of a cell (as a result of these hemagglutinin molecules interacting with cellular receptor molecules), the cell engulfs the virus into an internal compartment called an endosome. Acidification of the endosome, part of the cell' s normal activity, triggers a sequence of conformational changes in the hemagglutinin molecules on the surface of the virus. One part of the hemagglutinin inserts itself into the endosomal membrane, and further conformational changes draw the endosomal and viral membranes together into an intermediate, hemifused state; the process then continues until fusion of the two membranes is complete. Previous work has suggested that an average of three hemagglutinin molecules are required to fuse the endosomal and viral membranes. Ivanovic et al. have now investigated the molecular details of this process and described the time course of conformational changes undergone by the hemagglutinin molecules from the moment the pH is lowered within the endosome until the time when hemifusion of the endosomal and viral membranes is complete. They find, among other things, that hemifusion proceeds rapidly only when three or four immediately adjacent hemagglutinin molecules have inserted into the endosomal membrane. Since membrane fusion is a very general cellular process, the findings of Ivanovic et al. are relevant to many areas of cell biology, in addition to having potential applications in virology.
keywords: structural biology and molecular biophysics
  microbiology and infectious disease
title: Influenza-virus membrane fusion by cooperative fold-back of stochastically induced hemagglutinin intermediates