|Home | About | Journals | Submit | Contact Us | Français|
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.
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.
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.
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.
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).
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).
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 and , 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 indeed accelerates hemifusion, and the reverse mutation retards it (Figure 4A). Moreover, like WT Udorn virions, virions do not show flow-induced motion at pH drop, while 76% virions do (Videos 5 and 6, Table 1). Kinetic parameters for arrest, derived from the tlag(pHdrop–arrest) distribution for the 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.
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 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 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.
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 and 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).
Glass coverslips were cleaned extensively by two rounds of sequential 20-min sonication in detergent, 1 M potassium hydroxide, acetone then ethanol. They were dried at 100°C for 30 min to 1 hr then plasma cleaned for 3 min at 500 mTorr O2 and 75 W (March Plasmod Plasma Etcher; March Instruments, Inc., Concord, CA). The final sonication in ethanol and plasma cleaning steps were done immediately before each experiment.
PDMS flow-cell fabrication was described previously (Ivanovic et al., 2012). Flow cells were affixed to clean-glass coverslips prior to planar membrane bilayer preparation as described (Floyd et al., 2008). Liposomes used to generate planar bilayers consisted of 4:4:2:0.1:2 × 10−4 ratio of 1, 2, dioleoyl-sn-glycero-3-phosphocholine (DOPC) (Avanti Polar Lipids, Alabaster, AL), 1-oleoyl-2-palmitoyl-sn-glycero-3-phosphocholine (POPC; Avanti Polar Lipids), cholesterol (Avanti Polar Lipids), bovine brain disialoganglioside GD1a (Sigma, St. Louis, MO), and N-((6-(biotinoyl)amino)hexanoyl)-1,2-dihexadecanoyl-sn-glycero-3-phosphoethanolamine(biotin-X DHPE) (Molecular Probes, Life Technologies, Grand Island, NY), where sialic acid residues on GD1a served as receptors for influenza virions. We extruded liposomes through a polycarbonate membrane filter with a 200-nm pore size, loaded liposomes into flow cells at 0.03 ml/min for 1 min, then stopped the flow and allowed 15 min for bilayers to form by vesicle-spreading method (Nollert et al., 1995). The bilayers were washed with neutral pH buffer (50 mM HepesNaOH pH 7.4, 137 mM NaCl and 0.2 mM EDTA) flowing at 0.05 ml/min for 2 min. They were labeled with fluorescent-labeled streptavidin (30 µg/ml; Invitrogen) for 5 min then washed again with neutral pH buffer.
We used the protocol for influenza membrane fusion described previously (Floyd et al., 2008) with modifications. We labeled virions in 5- to 20-µl reactions with a lipophilic dye Octadecyl Rhodamine B Chloride (R18; Invitrogen, Life Technologies, Grand Island, NY) at a quenching concentration, 40 µM, for 2 hr at room temperature. Concentration of total viral protein in virion labeling reactions was between 0.25 and 1 mg/ml. Unincorporated dye was removed by gel filtration using PD-10 desalting column (GE Healthcare) according to manufacturer's recommendations using the neutral pH buffer, and three most concentrated, 250-µl fractions were pooled. Labeled virions were used within 48 hr of labeling. Virions were loaded into the flow cell at 0.05 ml/min for 30 s, and unattached virions were washed with neutral pH buffer within 2 min of introducing the virus into the flow cell. The pH in the flow cell was dropped at 7 min post virus introduction using citrate buffer (10 mM citric acid, 140 mM NaCl and 0.1 mM EDTA) having indicated pH. The flow rate of the low-pH buffer was kept constant at 0.075 ml/min for the duration of imaging, except in a single experiment (Video 2), where the flow was stopped to allow switching from low-pH back to neutral-pH buffer flow. All experiments were performed at 23°C.
We used the microscope configuration described previously (Floyd et al., 2008) except that 488-nm laser (Obis, Coherent, Santa Clara, CA) used to excite fluorescein and 552-nm laser (Obis, Coherent) used to excite R18 were used simultaneously and each was projected onto an entire 512 × 512 pixel EM-CCD sensor (Model C9100-13; Hamamatsu, Bridgewater, NJ). The power of the 488-nm laser was 1–2 µW and that of the 552-nm laser was 0.5–1.5 µW over the approximately 160 µm in diameter area. Acquisition times were between 0.5 and 2 s, depending on pH.
We used MATLAB (MathWorks, Natick, MA) for all data analysis.
The reference time point (t = 0) for our kinetic measurements was the time at which most of the fluorescein signal had dissipated, determined as follows. We integrated the total optical fluorescence power in each imaged field from the onset of imaging to a time when at most a small subset of virions had hemifused. We generated fluorescence-vs-time plots of these data. We fitted the transition step resulting from fluorescein dissipation with the integral of a Gaussian function with a step-down shape, as previously described for dissipation of internal virion fluorescein (Ivanovic et al., 2012), using the following model:
where h is the height of the intensity drop, w is the half-width of the transition (in s), tc is the time at which the transition is half complete, yofs is the residual intensity after dissipation, and m is a decay rate (slower than 1/w); erfc(), the complementary error function, is the integral of a Gaussian, defined as
(see below for the definition of erf(t)). We define t0 as (tc+w).
Individual virion positions after their arrest or after the pH drop (for static virions) were detected automatically and their fluorescence trajectories extracted as previously described (Floyd et al., 2008). We fitted the local intensity jump due to virion arrival and arrest at this location (see Figure 1B) with the integral of a Gaussian with a step-up shape:
where erf(), the error function, is defined as
and the remaining parameters are as described above for the model used to fit fluorescein dissipation. We define ta as (tc + w) − t0.
Hemifusion time for individual virions was determined as previously described (Floyd et al., 2008) using t0 described above as the reference point.
We fitted tlag(pHdrop–arrest), tlag(arrest–hemifusion), tlag(pHdrop–hemifusion) with a gamma distribution, , to derive k and/or N (Figures 1E,F, 2B, 3B,E and Table 1), or to derive k while keeping N fixed (see below) (Figures 1G, 2C). The 95% confidence intervals for the fits were automatically generated.
Our current estimates for N describing either arrest to hemifusion for X31-HA virions or pH drop to hemifusion for Udorn virions was between 3 and 4 (Figures 1F, 2B and Table 1), similar to previously published value of 3 for X31 virions (Floyd et al., 2008). All hemifusion k values in Figures 1G and 2C were derived from gamma-distribution fitting of tlag(arrest–hemifusion) (X31-HA virions) or tlag(pHdrop–hemifusion) (Udorn virions) distributions with fixed N = 3 and free-fitted k.
Because the time for the external pH to drop from neutral to a target value in our flow cell was significant relative to the average tlag(pHdrop–arrest), especially at high proton concentrations (the duration of the pH transition was equal to the mean tlag(pHdrop–arrest) at pH 4 and 4.5 and about one-tenth of that interval at pH 5.65), the true value of N for pH drop to arrest is approached by experiments at low proton concentrations. We derive a value for N of about four from the tlag(pHdrop–arrest) distributions at low proton concentrations (2 to 6 µM or pH 5.65–5.2) (Figure 1F). We get a more accurate estimate for N after pooling tlag(pHdrop–arrest) data from three independent experiments at pH 5.5 (n = 1215) to obtain a value of 3.4 ± 0.2 (Figure 1C, left, and Table 1). All arrest k values in Figure 1G were derived from gamma-distribution fitting of tlag(pHdrop–arrest) distributions with fixed N = 3 and free-fitted k, and thus represent a lower limit estimate.
The efficiency of arrest for all virions at all analyzed pHs is close to 100%. For approximately 20% of X31-HA WT and virions and nearly 100% of Udorn WT and all destabilized mutant virions (, and ), we failed to detect significant movement following pH drop, even in the higher pH range. In some cases, we cannot unambiguously distinguish between (1) virion arrest due to virion or bilayer defects, such as presence of virion-associated denatured HA or parts of the target bilayer that are not continuous with the rest of the membrane surface (e.g., remaining unfused liposomes), and (2) genuine HA triggering and membrane insertion that precede completion of the pH drop and might even occur with low frequency at neutral pH. One reason for this limitation is that virion movement under hydrodynamic force is somewhat retarded at neutral pH even in the context of X31-HA WT virions (Video 9); upon pH drop, X31-HA WT virions initially accelerate and then arrest (Video 10). HA triggering and membrane insertion at neutral pH is probably a rare event for all virions analyzed in this study but one that might occur more frequently for Udorn WT and other virions with a destabilized HA. Indeed, we occasionally observe hemifusion at neutral pH even with X31-HA WT virions (see Video 9).
We believe that the fraction (approximately 20%) of X31-HA WT and virions that never move under hydrodynamic force at high pH mainly results from virion or bilayer defects and not from HA triggering at neutral pH for three reasons: (1) the observed fraction of immobile virions saturates above pH 5.5, while a modest increase in pH to 5.65 has a twofold to threefold effect on mean lag times for both arrest and hemifusion; (2) hemifusion kinetics are not significantly different for pre-arrested and mobile populations, and a majority of virions (65% to 80%) from both populations undergo hemifusion (see Table 1); (3) the observed fraction of immobile virions varies to some extent (15–30%) even when mean arrest and hemifusion times are unchanged from experiment to experiment.
For Udorn WT and destabilized mutants, a combination of factors might explain the nearly complete absence of motion at pH drop: some fraction of HAs might release fusion peptides even at neutral pH (although probably not very frequently, as neutral pH fusion is still a rare event), and they might arrest before low pH-induced motion takes effect (see Video 10 for low pH-induced motion of X31-HA WT virions).
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, , 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.
We thank L'Oreal and AAAS for the 2011 Women in Science Fellowship to TI that enabled construction of a custom TIRF microscope used in the current experiments, Branislav Ivanovic for 3D drawings of HA monomer (Figure 3) and the model of influenza membrane fusion (Figure 6), and Milos Popovic for valuable input and discussion. Robert Garcea generously made his laboratory and equipment available to TI for this work.
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
This paper was supported by the following grants:
The authors declare that no competing interests exist.
TI, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article.
JLC, Acquisition of data.
SPW, Conception and design, Analysis and interpretation of data.
AMvO, Conception and design, Analysis and interpretation of data, Drafting or revising the article.
SCH, Conception and design, Analysis and interpretation of data, Drafting or revising the article.
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for choosing to send your work entitled “Influenza-virus membrane fusion: cooperative fold-back of stochastically induced hemagglutinin intermediates” for consideration at eLife. Your article has been evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.
The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments based on the reviewers' reports.
This paper is an important contribution towards understanding the molecular mechanism of HA fusion. Harrison, van Oijen, and collaborators developed a single-virion fusion assay four years ago to observe the sequence of events that lead to hemifusion (lipid mixing) and complete fusion (content mixing) induced by influenza virus hemagglutinin (HA) upon a controlled pH drop. In the present work they further developed the assay to also monitor the time it takes for individual virions to stably dock to a supported bilayer (referred to as “arrest”). Furthermore, the assay was streamlined in order to test a number of strategic mutations of HA.
The chosen mutations destabilize the interaction between the fusion peptide and the rest of HA in the pre-fusion state (the structure of which has been determined by X-ray crystallography). The single-virion fusion data suggest that the kinetics of both virion arrest and hemifusion are accelerated by the destabilizing mutations, while the apparent number of events (~3 trimers) remains unchanged. This result illustrates that release of the fusion peptide from its pocket is rate-limiting for virion arrest and hemifusion, and it is consistent with a model of an extended HA intermediate that bridges membranes prior to a further conformational change that leads to hemifusion.
The key insight in this paper is to make a link between the kinetic intermediates that lead to virion arrest and hemifusion with the structural requirement of fusion peptide release.
1. A subset of the virons for X31-HA is already arrested prior to the pH drop. Some speculation about the origin of this phenomenon might be warranted. Perhaps some fusion peptides have been released for this variant even at higher pH?
2. What is the efficiency of arrest upon the pH drop for the X31-HA variant? In other words, how many virions transiently dock to the bilayer but do arrest and hemifuse (Movie 1)? Some rough estimation would be sufficient.
3. Likewise, what is the estimated efficiency of hemifusion for those virions that are arrested?
4. As far as we understand it, the final model posits that each trimer releases fusion peptide stochastically, and if 3–4 of them out of about 50 or more trimers do so (regardless of proximity), the virion arrests. The remaining available trimers continue to release fusion peptide at the same individual rate and if this results in peptide release of 3–4 neighboring ones, then the hemi-fusion occurs.
In this model, however, it is not intuitively clear that the kinetics from arrest to hemi-fusion would show a gamma distribution with N ~ 3–4. The reaction kinetics behind the gamma distribution assumes multiple reaction steps to occur sequentially or parallel. The steps leading to virion arrest may well be described by a gamma distribution because one starts with a large number of trimers and fusion peptide release of the first few would occur with essentially the same global rate. However, the situation is different for the kinetics from arrest to hemi-fusion because the model proposes the release of peptides from 3–4 neighboring trimers.
Could the observed lag-time distributions be reproduced with a simple simulation? For example, one could start with a triangular lattice of a certain size and simulate the entire process based on the proposed model. Each trimer could release peptide with a rate of k, and if N of them release peptide, one calls it an arrest; and then if N of them accumulate next to each other, one calls it hemi-fusion. The simulation could then be repeated with different values for N. Such a simulation would be more directly related to the model that is proposed here.
5. It would be nice (but not required) to investigate the effect of the same mutations on content mixing using the content marker that was used in the previous PNAS study. In that previous work, it was concluded that the transition from hemifusion to full fusion is a one step process. One would then expect that the step from hemifusion to full fusion would not be influenced by the HA mutations studied here unless an additional HA trimer is actively engaged to achieve complete fusion.
1. A subset of the virons for X31-HA is already arrested prior to the pH drop. Some speculation about the origin of this phenomenon might be warranted. Perhaps some fusion peptides have been released for this variant even at higher pH?
We have added to the Results and Materials and methods statements to address this observation and speculate about its cause. We have also added two additional movies (Video 9 and Video 10), which speak to this point. In summary, while it is possible and likely that some fusion peptides are released and inserted into the target membrane even at neutral pH (we occasionally observe hemifusion events at neutral pH – see Video 9), we believe that this is a rare occurrence and that most of the pre-arrested events in the context of X31-HA virions at pH 5.5 and above (~20%) result from imperfections in either the bilayer or the virions themselves.
2. What is the efficiency of arrest upon the pH drop for the X31-HA variant? In other words, how many virions transiently dock to the bilayer but do arrest and hemifuse (Movie 1)? Some rough estimation would be sufficient.
3. Likewise, what is the estimated efficiency of hemifusion for those virions that are arrested?
Response to points 2 & 3: Virion arrest preceded hemifusion in all cases and hemifusion followed virion arrest in 75% (X31-HA virions) and 65% (UdornS4G virions) of the cases at pH 5.5. We had originally included this measurement in Table 1, and we have now added a sentence in the Materials and methods further discussing these percentages. We also added a sentence about the relative efficiency of arrest, which was nearly 100% up to the highest pH analyzed (pH 5.65). At pH 5.65 we observe a small subset of virions that continue to move under hydrodynamic force well after the pH drop, but in this case also nearly 100% seem to arrest within a few minutes of pH drop. We cannot estimate the exact percentage of those virions that never arrest or that arrest very late, because we do not have a way to automatically detect and track the virions that never arrest, but careful observation of the movies in slow motion allowed us to conclude that ~100% of virions do arrest at all the pHs we analyzed.
From the subset of virions that were arrested preceding pH drop, 82% (X31-HA virions) and 81% (UdornS4G virions) hemifused at pH 5.5. We have now added these percentages to Table 1.
4. […] Could the observed lag-time distributions be reproduced with a simple simulation? For example, one could start with a triangular lattice of a certain size and simulate the entire process based on the proposed model. Each trimer could release peptide with a rate of k, and if N of them release peptide, one calls it an arrest; and then if N of them accumulate next to each other, one calls it hemi-fusion. The simulation could then be repeated with different values for N. Such a simulation would be more directly related to the model that is proposed here.
This was an astute point and an excellent suggestion. We have done the simulation, added a new figure (now Figure 5), and we have added a new section to the Results (titled ‘Simulation of molecular events in the contact patch between virion and target membrane'). Briefly, we have generated a computer simulation of our predicted model that virion arrest results from independent triggering and irreversible membrane insertion from several HA trimers anywhere at the target membrane interface, and that hemifusion is limited by the same molecular rearrangements, fusion peptide release, and membrane insertion, but requires participation of neighboring HA trimers, and is thus related to a different probability. This exercise not only supported our original contention but narrowed down our interpretation of results, allowing us to state that virion arrest results from 3–4 trimers inserting into the target membrane anywhere in the contact patch and hemifusion from 3–4 neighboring HA trimers undergoing the same molecular rearrangement. In other words, our observed data are inconsistent with participation of as few as 2 or as many as 5 neighboring HA trimers in hemifusion under our current experimental conditions for a range of relevant contact–patch sizes.
Optional: 5. It would be nice (but not required) to investigate the effect of the same mutations on content mixing using the content marker that was used in the previous PNAS study. In that previous work, it was concluded that the transition from hemifusion to full fusion is a one step process. One would then expect that the step from hemifusion to full fusion would not be influenced by the HA mutations studied here unless an additional HA trimer is actively engaged to achieve complete fusion.
This is a good suggestion and one that should be addressed with future work, but we believe it is beyond the scope of the current manuscript.