Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

HIV and influenza fusion peptide interactions with (dis)ordered lipid bilayers: Understanding mechanisms and implications for antimicrobial and antiviral approaches

The interactions of viral fusion peptides from influenza (E4K and Ac-E4K) and human immunodeficiency virus (gp41 and Ac-gp41) with planar lipid bilayers and monolayers was investigated herein. A combination of surface-sensitive techniques, including quartz crystal microbalance with dissipation (QCM-D), Langmuir-Blodgett area-pressure isotherms with Micro-Brewster angle microscopy, and neutron reflectometry, was employed. Differences in the interactions of the viral fusion peptides with lipid bilayers featuring ordered and disordered phases, as well as lipid rafts, were revealed. The HIV fusion peptide (gp41) exhibited strong binding to the DOPC/DOPS bilayer, comprising a liquid disordered phase, with neutron reflectometry (NR) showing interaction with the bilayer's headgroup area. Conversely, negligible binding was observed with lipid bilayers in a liquid ordered phase. Notably, the influenza peptide (E4K) demonstrated slower binding kinetics with DOPC/DOPS bilayers and distinct interactions compared to gp41, as observed through QCM-D. This suggests different mechanisms of interaction with the lipid bilayers: one peptide interacts more within the headgroup region, while the other is more involved in transmembrane interactions. These findings hold implications for understanding viral fusion mechanisms and developing antimicrobials and antivirals targeting membrane interactions. The differential binding behaviours of the viral fusion peptides underscore the importance of considering membrane composition and properties in therapeutic strategy design.

Two modes of fusogenic action for influenza virus fusion peptide

The entry of influenza virus into the host cell requires fusion of its lipid envelope with the host membrane. It is catalysed by viral hemagglutinin protein, whose fragments called fusion peptides become inserted into the target bilayer and initiate its merging with the viral membrane. Isolated fusion peptides are already capable of inducing lipid mixing between liposomes. Years of studies indicate that upon membrane binding they form bend helical structure whose degree of opening fluctuates between tightly closed hairpin and an extended boomerang. The actual way in which they initiate fusion remains elusive. In this work we employ atomistic simulations of wild type and fusion inactive W14A mutant of influenza fusion peptides confined between two closely apposed lipid bilayers. We characterise peptide induced membrane perturbation and determine the potential of mean force for the formation of the first fusion intermediate, an interbilayer lipid bridge called stalk. Our results demonstrate two routes through which the peptides can lower free energy barrier towards fusion. The first one assumes peptides capability to adopt transmembrane configuration which subsequently promotes the creation of a stalk-hole complex. The second involves surface bound peptide configuration and proceeds owing to its ability to stabilise stalk by fitting into the region of extreme negative membrane curvature resulting from its formation. In both cases, the active peptide conformation corresponds to tight helical hairpin, whereas extended boomerang geometry appears to be unable to provide favourable thermodynamic effect. The latter observation offers plausible explanation for long known inactivity of boomerang-stabilising W14A mutation.

Computational methods to study enveloped viral entry

The protein-membrane interactions that mediate viral infection occur via loosely ordered, transient assemblies, creating challenges for high-resolution structure determination. Computational methods and in particular molecular dynamics simulation have thus become important adjuncts for integrating experimental data, developing mechanistic models, and suggesting testable hypotheses regarding viral function. However, the large molecular scales of virus-host interaction also create challenges for detailed molecular simulation. For this reason, continuum membrane models have played a large historical role, although they have become less favored for high-resolution models of protein assemblies and lipid organization. Here, we review recent progress in the field, with an emphasis on the insight that has been gained using a mixture of coarse-grained and atomic-resolution molecular dynamics simulations. Based on successes and challenges to date, we suggest a multiresolution strategy that should yield the best mixture of computational efficiency and physical fidelity. This strategy may facilitate further simulations of viral entry by a broader range of viruses, helping illuminate the diversity of viral entry strategies and the essential common elements that can be targeted for antiviral therapies.

Molecular mechanisms of the influenza fusion peptide: insights from experimental and simulation studies

A key step in infections by enveloped viruses, such as influenza, is the fusion between the viral envelope and the host cell membrane, which allows the virus to insert its genetic material into the host cell and replicate. The influenza virus fusion process is promoted by hemagglutinin (HA), a glycoprotein that contains three identical monomers composed of two polypeptide chains (HA1 and HA2). Early studies on this protein revealed that HA-mediated fusion involves the insertion of the HA2 N-terminal segment into the host membrane and that this segment, known as the fusion peptide, is a key player in the fusion process. This mini-review highlights the main findings that have been obtained by experimental and computational studies on the HA fusion peptide, which give us a glimpse of its mode of action.

Transient Excursions to Membrane Core as Determinants of Influenza Virus Fusion Peptide Activity

Fusion of viral and host cell membranes is a critical step in the life cycle of enveloped viruses. In the case of influenza virus, it is mediated by subunit 2 of hemagglutinin (HA) glycoprotein whose N-terminal fragments insert into the target membrane and initiate lipid exchange. These isolated fragments, known as fusion peptides (HAfp), already possess own fusogenic activity towards liposomes. Although they have long been studied with the hope to uncover the details of HA-mediated fusion, their actual mechanism of action remains elusive. Here, we use extensive molecular dynamics simulations combined with experimental studies of three HAfp variants to fully characterize their free energy landscape and interaction with lipid bilayer. In addition to customary assumed peptides localization at lipid-water interface, we characterize membrane-spanning configurations, which turn out to be metastable for active HAfps and unstable for the fusion inactive W14A mutant. We show that, while the degree of membrane perturbation by surface peptide configurations is relatively low and does not show any mutation-related differences, the effect of deeply inserted configurations is significant and correlates with insertion depth of the N-terminal amino group which is the highest for the wild type HAfp. Finally, we demonstrate the feasibility of spontaneous peptide transition to intramembrane location and the critical role of strictly conserved tryptofan residue 14 in this process.

Computational Insight Into the Mechanism of SARS-CoV-2 Membrane Fusion

Membrane fusion, a key step in the early stages of virus propagation, allows the release of the viral genome in the host cell cytoplasm. The process is initiated by fusion peptides that are small, hydrophobic components of viral membrane-embedded glycoproteins and are typically conserved within virus families. Here, we attempted to identify the correct fusion peptide region in the Spike protein of SARS-CoV-2 by all-atom molecular dynamics simulations of dual membrane systems with varied oligomeric units of putative candidate peptides. Of all of the systems tested, only a trimeric unit of a 40-amino-acid region (residues 816-855 of SARS-CoV-2 Spike) was effective in triggering the initial stages of membrane fusion, within 200 ns of simulation time. Association of this trimeric unit with dual membranes resulted in the migration of lipids from the upper leaflet of the lower bilayer toward the lower leaflet of the upper bilayer to create a structural unit reminiscent of a fusion bridge. We submit that residues 816-855 of Spike represent the bona fide fusion peptide of SARS-CoV-2 and that computational methods represent an effective way to identify fusion peptides in viral glycoproteins.

Intracellular Vesicle Fusion Requires a Membrane-Destabilizing Peptide Located at the Juxtamembrane Region of the v-SNARE

Intracellular vesicle fusion is mediated by soluble N-ethylmaleimide sensitive factor attachment protein receptors (SNAREs) and Sec1/Munc18 (SM) proteins. It is generally accepted that membrane fusion occurs when the vesicle and target membranes are brought into close proximity by SNAREs and SM proteins. In this work, we demonstrate that, for fusion to occur, membrane bilayers must be destabilized by a conserved membrane-embedded motif located at the juxtamembrane region of the vesicle-anchored v-SNARE. Comprised of basic and hydrophobic residues, the juxtamembrane motif perturbs the lipid bilayer structure and promotes SNARE-SM-mediated membrane fusion. The juxtamembrane motif can be functionally substituted with an unrelated membrane-disrupting peptide in the membrane fusion reaction. These findings establish the juxtamembrane motif of the v-SNARE as a membrane-destabilizing peptide. Requirement of membrane-destabilizing peptides is likely a common feature of biological membrane fusion.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

Charged N-terminus of Influenza Fusion Peptide Facilitates Membrane Fusion

Cleavage of hemagglutinin precursor (HA0) by cellular proteases results in the formation of two subunits, HA1 and HA2. The N-terminal fragment of HA2, named a fusion peptide (HAfp), possess a charged, amine N-terminus. It has been shown that the N-terminus of HAfp stabilizes the structure of a helical hairpin observed for a 23-amino acid long peptide (HAfp1-23), whose larger activity than HAfp1-20 has been demonstrated recently. In this paper, we analyze the effect of N-terminal charge on peptide-mediated fusion efficiency and conformation changes at the membrane interface by comparison with the corresponding N-acetylated peptides of 20- and 23-amino acid lengths. We found that higher fusogenic activities of peptides with unmodified amino termini correlates with their ability to form helical hairpin structures oriented perpendicularly to the membrane plane. Molecular dynamics simulations showed that acetylated peptides adopt open and surface-bound conformation more often, which induced less disorder of the phospholipid chains, as compared to species with unmodified amino termini.

Three conserved C-terminal residues of influenza fusion peptide alter its behavior at the membrane interface

The N-terminal fragment of the viral hemagglutinin HA2 subunit is termed a fusion peptide (HAfp). The 23-amino acid peptide (HAfp1-23) contains three C-terminal W21-Y22-G23 residues which are highly conserved among serotypes of influenza A and has been shown to form a tight helical hairpin very distinct from the boomerang structure of HAfp1-20. We studied the effect of peptide length on fusion properties, structural dynamics, and binding to the membrane interface. We developed a novel fusion visualization assay based on FLIM microscopy on giant unilamellar vesicles (GUV). By means of molecular dynamics simulations and spectroscopic measurements, we show that the presence of the three C-terminal W21-Y22-G23 residues promotes the hairpin formation, which orients perpendicularly to the membrane plane and induces more disorder in the surrounding lipids than the less structured HAfp1-20. Moreover, we report cholesterol-enriched domain formation induced exclusively by the longer fusion peptide.

Fusing simulation and experiment: The effect of mutations on the structure and activity of the influenza fusion peptide

During the infection process, the influenza fusion peptide (FP) inserts into the host membrane, playing a crucial role in the fusion process between the viral and host membranes. In this work we used a combination of simulation and experimental techniques to analyse the molecular details of this process, which are largely unknown. Although the FP structure has been obtained by NMR in detergent micelles, there is no atomic structure information in membranes. To answer this question, we performed bias-exchange metadynamics (BE-META) simulations, which showed that the lowest energy states of the membrane-inserted FP correspond to helical-hairpin conformations similar to that observed in micelles. BE-META simulations of the G1V, W14A, G12A/G13A and G4A/G8A/G16A/G20A mutants revealed that all the mutations affect the peptide’s free energy landscape. A FRET-based analysis showed that all the mutants had a reduced fusogenic activity relative to the WT, in particular the mutants G12A/G13A and G4A/G8A/G16A/G20A. According to our results, one of the major causes of the lower activity of these mutants is their lower membrane affinity, which results in a lower concentration of peptide in the bilayer. These findings contribute to a better understanding of the influenza fusion process and open new routes for future studies.

Capturing Spontaneous Membrane Insertion of the Influenza Virus Hemagglutinin Fusion Peptide

Hemagglutinin (HA) is a protein located on the surface of the influenza virus that mediates viral fusion to the host cellular membrane. During the fusion process the HA fusion peptide (HAfp), formed by the first 23 N-terminal residues of HA and structurally characterized by two alpha helices (Helix A and Helix B) tightly packed in a hairpin-like arrangement, is the only part of the virus in direct contact with the host membrane. After encountering the host cell, HAfp is believed to insert into the membrane, thereby initiating the fusion of the viral and host membranes. Detailed characterization of the interactions between the HAfp and cellular membrane is therefore of high relevance to the mechanism of viral entry into the host cell. Employing HMMM membrane representation with enhanced lipid mobility, we have performed a large set of independent simulations of unbiased membrane binding of HAfp. We have been able to capture spontaneous binding and insertion of HAfp consistently in nearly all the simulations. A reproducible membrane-bound configuration emerges from these simulations, despite employing a diverse set of initial configurations. Extension of several of the simulations into full membrane systems confirms the stability of the membrane-bound form obtained from HMMM binding simulations. The resulting model allows for the characterization of important interactions between the peptide and the membrane and the details of the binding process of the peptide for the first time. Upon membrane binding, Helix A inserts much deeper into the membrane than Helix B, suggesting that the former is responsible for hydrophobic anchoring of the peptide into the membrane. Helix B, in contrast, is found to establish major amphipathic interactions at the interfacial region thereby contributing to binding strength of HAfp.

Self-assembly molecular dynamics simulations shed light into the interaction of the influenza fusion Peptide with a membrane bilayer

Influenza virus is one of the most devastating human pathogens. In order to infect host cells, this virus fuses its membrane with the host membrane in a process mediated by the glycoprotein hemagglutinin. During fusion, the N-terminal region of hemagglutinin, which is known as the fusion peptide (FP), inserts into the host membrane, promoting lipid mixing between the viral and host membranes. Therefore, this peptide plays a key role in the fusion process, but the exact mechanism by which it promotes lipid mixing is still unclear. To shed light into this matter, we performed molecular dynamics (MD) simulations of the influenza FP in different environments (water, dodecylphosphocholine (DPC) micelles, and a dimyristoylphosphatidylcholine (DMPC) membrane). While in pure water the peptide lost its initial secondary structure, in simulations performed in the presence of DPC micelles it remained stable, in agreement with previous experimental observations. In simulations performed in the presence of a preassembled DMPC bilayer, the peptide became unstructured and was unable to insert into the membrane as a result of technical limitations of the method used. To overcome this problem, we used a self-assembly strategy, assembling the membrane together with the peptide. These simulations revealed that the peptide can adopt a membrane-spanning conformation, which had not been predicted by previous MD simulation studies. The peptide insertion had a strong effect on the membrane, lowering the bilayer thickness, disordering nearby lipids, and promoting lipid tail protrusion. These results contribute to a better understanding of the role of the FP in the fusion process.

Molecular dynamics, monte carlo simulations, and langevin dynamics: a computational review

Macromolecular structures, such as neuraminidases, hemagglutinins, and monoclonal antibodies, are not rigid entities. Rather, they are characterised by their flexibility, which is the result of the interaction and collective motion of their constituent atoms. This conformational diversity has a significant impact on their physicochemical and biological properties. Among these are their structural stability, the transport of ions through the M2 channel, drug resistance, macromolecular docking, binding energy, and rational epitope design. To assess these properties and to calculate the associated thermodynamical observables, the conformational space must be efficiently sampled and the dynamic of the constituent atoms must be simulated. This paper presents algorithms and techniques that address the abovementioned issues. To this end, a computational review of molecular dynamics, Monte Carlo simulations, Langevin dynamics, and free energy calculation is presented. The exposition is made from first principles to promote a better understanding of the potentialities, limitations, applications, and interrelations of these computational methods.

Computer-Aided Approaches for Targeting HIVgp41

Virus-cell fusion is the primary means by which the human immunodeficiency virus-1 (HIV) delivers its genetic material into the human T-cell host. Fusion is mediated in large part by the viral glycoprotein 41 (gp41) which advances through four distinct conformational states: (i) native, (ii) pre-hairpin intermediate, (iii) fusion active (fusogenic), and (iv) post-fusion. The pre-hairpin intermediate is a particularly attractive step for therapeutic intervention given that gp41 N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) domains are transiently exposed prior to the formation of a six-helix bundle required for fusion. Most peptide-based inhibitors, including the FDA-approved drug T20, target the intermediate and there are significant efforts to develop small molecule alternatives. Here, we review current approaches to studying interactions of inhibitors with gp41 with an emphasis on atomic-level computer modeling methods including molecular dynamics, free energy analysis, and docking. Atomistic modeling yields a unique level of structural and energetic detail, complementary to experimental approaches, which will be important for the design of improved next generation anti-HIV drugs.

Structure and dynamics of a fusion peptide helical hairpin on the membrane surface: comparison of molecular simulations and NMR

The conserved N-terminal residues of the HA2 subunit of influenza hemagglutinin (fusion peptide) are essential for membrane fusion and viral entry. Recent NMR studies showed that the 23-residue fusion peptide forms a helical hairpin that undergoes rocking motion relative to the membrane surface on a nanosecond time scale. To compare with NMR and to obtain a detailed molecular picture of the peptide-membrane interaction, we performed molecular dynamics simulations of the fusion peptide in explicit dimyristoylphosphatidylcholine and with the IMM1 implicit membrane model. To account for low and neutral pH conditions, simulations were performed with acidic groups (E11 and D19) protonated and unprotonated, respectively. The hairpin structure was stable in the simulations, with the N-terminal helix buried more deeply into the hydrophobic membrane interior than the C-terminal helix. Interactions between the tryptophans in the fusion peptide and phospholipid residues contribute to peptide orientation. Higher flexibility of the hairpin was observed in the implicit membrane simulations. Internal correlation functions of backbone N-H vectors were fit to the extended Lipari-Szabo model-free approach to obtain order parameters and correlation times. Good agreement with the NMR results was obtained for orientational fluctuations around the hairpin axis (rotation), but those around the perpendicular axis (tilting) were more limited in the simulations than inferred from the NMR experiments.

The influenza fusion peptide promotes lipid polar head intrusion through hydrogen bonding with phosphates and N-terminal membrane insertion depth

Influenza infection requires fusion between the virus envelope and a host cell endosomal membrane. The influenza hemagglutinin fusion peptide (FP) is essential to viral membrane fusion. It was recently proposed that FPs would fuse membranes by increasing lipid tail protrusion, a membrane fusion transition state. The details of how FPs induce lipid tail protrusion, however, remain to be elucidated. To decipher the molecular mechanism by which FPs promote lipid tail protrusion, we performed molecular dynamics simulations of the wild-type (WT) FP, fusogenic mutant F9A, and nonfusogenic mutant W14A in model bilayers. This article presents the peptide-lipid interaction responsible for lipid tail protrusion and a related lipid perturbation, polar head intrusion, where polar heads are sunk under the membrane surface. The backbone amides from the four N-terminal peptide residues, deeply inserted in the membrane, promoted both perturbations through H bonding with lipid phosphates. Polar head intrusion correlated with peptides N-terminal insertion depth and activity: the N-termini of WT and F9A were inserted deeper into the membrane than nonfusogenic W14A. Based on these results, we propose that FP-induced polar head intrusion would complement lipid tail protrusion in catalyzing membrane fusion by reducing repulsions between juxtaposed membranes headgroups. The presented model provides a framework for further research on membrane fusion and influenza antivirals.

Plasticity and conformational equilibria of influenza fusion peptides in model lipid bilayers

Membrane fusion is critical to eukaryotic cellular function and crucial to the entry of enveloped viruses such as influenza and human immunodeficiency virus. Influenza viral entry in the host cell is mediated by a 20-23 amino acid long sequence, called the fusion peptide. In the last years, possible structures for the fusion peptide and their implication in the membrane fusion initiation have been proposed; these ranging from an inverted V shaped α-helical structure to an α-helical hairpin, or to a complete α-helix. Here we develop a coarse grained approach to describe effectively the plasticity of the fusion peptide and the explored conformational states. We describe also a trimeric assembly for the fusion peptide and analyse the explored states in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine model membrane. For the single fusion peptide systems the kink angle observed experimentally for the V shaped structure shows a strong correlation with the orientation of the fusion peptide within the lipid bilayer. The trimeric fusion peptide model also experiences different conformational states and represents a more realistic model for the anchoring mechanism of one influenza haemagglutinin molecule. This article is part of a Special Issue entitled: Viral Membrane Proteins - Channels for Cellular Networking.

Lipid tail protrusion in simulations predicts fusogenic activity of influenza fusion peptide mutants and conformational models

Fusion peptides from influenza hemagglutinin act on membranes to promote membrane fusion, but the mechanism by which they do so remains unknown. Recent theoretical work has suggested that contact of protruding lipid tails may be an important feature of the transition state for membrane fusion. If this is so, then influenza fusion peptides would be expected to promote tail protrusion in proportion to the ability of the corresponding full-length hemagglutinin to drive lipid mixing in fusion assays. We have performed molecular dynamics simulations of influenza fusion peptides in lipid bilayers, comparing the X-31 influenza strain against a series of N-terminal mutants. As hypothesized, the probability of lipid tail protrusion correlates well with the lipid mixing rate induced by each mutant. This supports the conclusion that tail protrusion is important to the transition state for fusion. Furthermore, it suggests that tail protrusion can be used to examine how fusion peptides might interact with membranes to promote fusion. Previous models for native influenza fusion peptide structure in membranes include a kinked helix, a straight helix, and a helical hairpin. Our simulations visit each of these conformations. Thus, the free energy differences between each are likely low enough that specifics of the membrane environment and peptide construct may be sufficient to modulate the equilibrium between them. However, the kinked helix promotes lipid tail protrusion in our simulations much more strongly than the other two structures. We therefore predict that the kinked helix is the most fusogenic of these three conformations.

Membrane fusion is a common process critical to both cellular function and infection by enveloped viruses. Influenza is a particularly useful model system for studying fusion because the fusion reaction is accomplished by a single protein, hemagglutinin. Furthermore, mutations to the membrane-inserted portion of hemagglutinin have been identified that do not detectably alter the rest of the protein but can either arrest fusion halfway or block it entirely. For influenza at least, it seems that the membrane-inserted hemagglutinin peptide plays a critical role in promoting fusion, perhaps by increasing the local disorder of lipid bilayers. However, we lack a mechanistic understanding sufficient to predict the activity of fusion peptide mutants from their sequence. Here, we have used lipid tail protrusion as a way to measure how much fusion peptides disorder their surrounding bilayer; we see a strong relationship between lipid tail protrusion and the ability of fusion peptide mutants to promote lipid mixing between membranes. Our simulations also predict that this lipid tail protrusion is much more common when the peptides adopt a kinked helix structure than when they are straight or hairpin-like. We therefore hypothesize that, while all three types of structure likely undergo conformational exchange, the kinked helix structure is most active in promoting fusion.

The influenza fusion peptide adopts a flexible flat V conformation in membranes

Knowledge about the influenza fusion peptide (FP) membrane insertion mode is crucial for understanding its fusogenic mechanism. NMR and electron paramagnetic resonance experiments showed that in micelles, the FP inserted as a fixed-angle inverted V. In membranes, however, it was shown to insert as a straight α-helix (by molecular-dynamics simulations) and to adopt multiple kinked conformations (by solid-state NMR). In this work we performed explicit-solvent molecular-dynamics simulations of the influenza FP, and its F9A and W14A mutants, in POPC membranes. The Hα1 chemical shifts predicted from the molecular-dynamics structures are in excellent agreement with the experimental values obtained for the three peptides. The peptide orientation and conformations observed from the simulations lead to a flexible flat-V model in which the peptide lies almost flat on the membrane surface and alternates between kinked and straight-helix conformations.

Helical hairpin structure of influenza hemagglutinin fusion peptide stabilized by charge-dipole interactions between the N-terminal amino group and the second helix

The fusion domain of the influenza coat protein hemagglutinin HA2, bound to dodecyl phosphocholine micelles, was recently shown to adopt a structure consisting of two antiparallel α-helices, packed in an exceptionally tight hairpin configuration. Four interhelical H(α) to C═O aliphatic H-bonds were identified as factors stabilizing this fold. Here, we report evidence for an additional stabilizing force: a strong charge-dipole interaction between the N-terminal Gly(1) amino group and the dipole moment of helix 2. pH titration of the amino-terminal (15)N resonance, using a methylene-TROSY-based 3D NMR experiment, and observation of Gly(1 13)C' show a strongly elevated pK = 8.8, considerably higher than expected for an N-terminal amino group in a lipophilic environment. Chemical shifts of three C-terminal carbonyl carbons of helix 2 titrate with the protonation state of Gly(1)-N, indicative of a close proximity between the N-terminal amino group and the axis of helix 2, providing an optimal charge-dipole stabilization of the antiparallel hairpin fold. pK values of the side-chain carboxylate groups of Glu(11) and Asp(19) are higher by about 1 and 0.5 unit, respectively, than commonly seen for solvent-exposed side chains in water-soluble proteins, indicative of dielectric constants of ε = ∼30 (Glu(11)) and ∼60 (Asp(19)), placing these groups in the headgroup region of the phospholipid micelle.

Single mutation effects on conformational change and membrane deformation of influenza hemagglutinin fusion peptides

The single mutation effect on the conformational change and membrane permeation of influenza hemagglutinin fusion peptides has been studied with molecular dynamics simulations. A total of seven peptides, including wild-type fusion peptide and its six single point mutants (G1E, G1S, G1V, G4V, E11A, and W14A, all with no fusion or hemifusion activity) are examined systematically, which covers a wide range of mutation sites as well as mutant residue types (both hydrophobic and hydrophilic). The wild-type shows a kink structure (inversed V-shape), which facilitates the interaction between the fusion peptide and the lipid bilayer, as well as the interaction between the two arms of the fusion peptide. All mutants show a strong tendency toward a linear alpha-helix conformation, with the initial kink structure in the wild-type broken. More interestingly, one of the key hydrophobic residues around the initial kink region, Phe-9, is found to flip away from the membrane surface in most of these mutants. This conformational change causes a loss of key interactions between the original two arms of the inversed V-shape of the wild-type, thus disabling the kink structure, which results in the stabilization of the linear alpha-helix structure. The fusion peptides also display significant impact on the membrane structure deformation. The thickness of the lipid bilayer surrounding the wild-type fusion peptide decreases significantly, which induces a positive curvature of lipid bilayer. All the single mutations examined here reduce this membrane structural deformation, supporting the fusion activity data from experiments.

Atomic-resolution simulations predict a transition state for vesicle fusion defined by contact of a few lipid tails

Membrane fusion is essential to both cellular vesicle trafficking and infection by enveloped viruses. While the fusion protein assemblies that catalyze fusion are readily identifiable, the specific activities of the proteins involved and nature of the membrane changes they induce remain unknown. Here, we use many atomic-resolution simulations of vesicle fusion to examine the molecular mechanisms for fusion in detail. We employ committor analysis for these million-atom vesicle fusion simulations to identify a transition state for fusion stalk formation. In our simulations, this transition state occurs when the bulk properties of each lipid bilayer remain in a lamellar state but a few hydrophobic tails bulge into the hydrophilic interface layer and make contact to nucleate a stalk. Additional simulations of influenza fusion peptides in lipid bilayers show that the peptides promote similar local protrusion of lipid tails. Comparing these two sets of simulations, we obtain a common set of structural changes between the transition state for stalk formation and the local environment of peptides known to catalyze fusion. Our results thus suggest that the specific molecular properties of individual lipids are highly important to vesicle fusion and yield an explicit structural model that could help explain the mechanism of catalysis by fusion proteins.

Membrane fusion is a common underlying process critical to neurotransmitter release, cellular trafficking, and infection by many viruses. Proteins have been identified that catalyze fusion, and mutations to these proteins have yielded important information on how fusion occurs. However, the precise mechanism by which membrane fusion begins is the subject of active investigation. We have used atomic-resolution simulations to model the process of vesicle fusion and to identify a transition state for the formation of an initial fusion stalk. Doing so required substantial technical advances in combining high-performance simulation and distributed computing to analyze the transition state of a complex reaction in a large system. The transition state we identify in our simulations involves specific structural changes by a few lipid molecules. We also simulate fusion peptides from influenza hemagglutinin and show that they promote the same structural changes as are required for fusion in our model. We therefore hypothesize that these changes to individual lipid molecules may explain a portion of the catalytic activity of fusion proteins such as influenza hemagglutinin.

Conformational sampling of influenza fusion peptide in membrane bilayers as a function of termini and protonation states

Influenza fusion peptide is critical for mediating the fusion of viral and host cell membranes during viral entry. The interaction of monomeric influenza fusion peptide with membranes is studied with replica exchange molecular dynamics simulations using a new implicit membrane model to effectively reach microsecond to millisecond time scales. The conformational sampling of the fusion peptide was studied as a function of different N- and C-termini, including an experimental construct with an additional C-terminal tag, as well as a function of protonation of acidic residues. It is found that the influenza fusion peptide mostly adopts helical structures with a pronounced kink at residues 11-13 with both N-terminal and C-terminal helices oriented mostly parallel to the membrane surface. A charged C-terminus and the presence of a charge C-terminal tag significantly alters the conformational sampling of the fusion peptide and results in more diverse conformational ensembles that include obliquely inserted N-terminal peptide structures. Protonation of acidic residues also affects the conformational sampling, however, based on pK(a) shift estimates the overall effect of pH = 5 on the conformational sampling of the influenza fusion peptide appears to be only minor.

13C-13C correlation spectroscopy of membrane-associated influenza virus fusion peptide strongly supports a helix-turn-helix motif and two turn conformations

The influenza virus fusion peptide (IFP) is the N-terminal domain of the viral hemagglutinin protein, binds to the endosomal membrane, and plays a critical role in fusion between the viral and endosomal membranes which is a primary step in infection. The IFP is also an important system for testing simulation methods for membrane-associated peptides. In detergent, the IFP forms helix-turn-helix and helix-turn-strand structures at pH 5.0 and 7.4, respectively, while simulations in membranes by different groups have yielded conflicting results with some reports of a continuous helix without a turn. In this study, (13)C-(13)C NMR correlation spectra were obtained for the membrane-associated IFP and the (13)C chemical shifts supported a helix-turn-helix motif at both pH 5.0 and 7.4 with an alternate turn conformation at pH 5.0 that was absent at pH 7.4. The alternate conformation was correlated with protonation of the side chain of Glu-11 in the turn and with greater fusion at pH 5.0. The structures are overall consistent with the hypothesis of "inverted V" membrane location of the IFP with insertion of the N-terminal region into the membrane and contact of the turn with the lipid/water interface. The positions of hydrophobic residues in the pH 5.0 structure may favor membrane insertion with resultant increased membrane perturbation and fusion rate. In addition to their functional relevance, these IFP structures are important reference data for simulations of the membrane-associated IFP which can in principle detect the full conformational distribution of the IFP.

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.

Evaluation of Glu11 and Gly8 of the H5N1 influenza hemagglutinin fusion peptide in membrane fusion using pseudotype virus and reverse genetics

Influenza viruses gain entry into host cells by binding to cellular receptors and promoting the fusion of the viral envelope with the host cell membrane. The fusion peptide of influenza hemagglutinin (HA) is crucial for fusion. To examine the structural and functional roles of amino acids E11 and G8 of the H5 HA fusion peptide, a series of fusion mutants was generated. We determined the effect of each mutation on fusion activity and infection of rescued recombinant virus by polykaryon formation, cell-cell fusion assay, HA pseudovirus transduction and reverse genetics. Our findings indicate that E11V and E11A mutants dramatically inhibit fusion and that at position 11 a polar residue such as glutamic acid or serine may be desirable for preserving the fusion activity. More interestingly, one mutation (G8E) raised the threshold pH of polykaryon formation. Our results suggest that G8 as well as E11 play an important functional and structural role in membrane fusion and that the polarity of E11 is crucial for fusion activity. Finally, we developed an assay based on a reporter gene plus pseudotyped virus that could sensitively detect fusion activity.

The Membrane-proximal Fusion Domain of HIV-1 GP41 Reveals Sequence-specific and Fine-tuning Mechanism of Membrane Binding

The membrane interface-partitioning region preceding the transmembrane anchor of the human immunodeficiency virus type 1 (HIV-1) gp41 envelope protein is one of the sites responsible for virus binding to its host cell membrane and subsequent fusion events. Here, we used molecular modeling techniques to assess membrane interactions, structure, and hydrophobic properties of the fusion-active peptide representing this region, several of its homologs from different HIV-1 strains, as well as a peptide - defective gp41 phenotype - unable to mediate cell-cell fusion and virus entry. It is shown that the wild-type peptides bind to the water-membrane interface in alpha-helical conformation, while the mutant adopts partly destabilized helix-break-helix structure on the membrane surface. The wild-type peptides reveal specific "tilted oblique-oriented" pattern of hydrophobicity on their surfaces - the property specific for fusion regions of other viruses. Fusion peptides penetrate into the membrane with their N-termini and reveal "fine-tuning" interactions with membrane and water environments: the shift of this balance (e.g., due to point mutations) may dramatically change the mode of membrane binding, and therefore, may cause loss of fusion activity. The modeling results agree well with experimental data and provide a strategy to delineate fusogenic regions in amino acid sequences of viral proteins.

[The current trend of avian influenza viruses in bioinformatics research]

OBJECTIVES

Since the first human infection from avian influenza was reported in Hong Kong in 1997, many Asian countries have confirmed outbreaks of highly pathogenic H5N1 avian influenza viruses. In addition to Asian countries, the EU authorities also held an urgent meeting in February 2006 at which it was agreed that Europe could also become the next target for H5N1 avian influenza in the near future. In this paper, we provide the general and applicable information on the avian influenza in the bioinformatics field to assist future studies in preventive medicine.

METHODS

We introduced some up-to-date analytical tools in bioinformatics research, and discussed the current trends of avian influenza outbreaks. Among the bioinformatics methods, we focused our interests on two topics: pattern analysis using the secondary database of avian influenza, and structural analysis using the molecular dynamics simulations in vaccine design.

RESULTS

Use of the public genome databases available in the bioinformatics field enabled intensive analysis of the genetic patterns. Moreover, molecular dynamic simulations have also undergone remarkable development on the basis of the high performance supercomputing infrastructure these days.

CONCLUSIONS

The bioinformatics techniques we introduced in this study may be useful in preventive medicine, especially in vaccine and drug discovery.

Configuration of influenza hemagglutinin fusion peptide monomers and oligomers in membranes

The 20 N-terminal residues of the HA2 subunit of influenza hemagglutinin (HA), known as the fusion peptide, play a crucial role in membrane fusion. Molecular dynamics simulations with implicit solvation are employed here to study the structure and orientation of the fusion peptide in membranes. As a monomer the alpha-helical peptide adopts a shallow, slightly tilted orientation along the lipid tail-head group interface. The average angle of the peptide with respect to membrane plane is 12.4 degrees . We find that the kinked structure proposed on the basis of NMR data is not stable in our model because of the high energy cost related to the membrane insertion of polar groups. Because hemagglutinin-mediated membrane fusion is promoted by low pH, we examined the effect of protonation of the Glu and Asp residues. The configurations of the protonated peptides were slightly deeper in the membrane but at similar angles. Finally, because HA is a trimer, we modeled helical fusion peptide trimers. We find that oligomerization affects the insertion depth of the peptide and its orientation with respect to the membrane: a trimer exhibits equally favorable configurations in which some or all of the helices in the bundle insert obliquely deep into the membrane.

Autocatalytic Activation of Influenza Hemagglutinin

Enveloped viruses contain surface proteins that mediate fusion between the viral and target cell membranes following an activating stimulus. Acidic pH induces the influenza virus fusion protein hemagglutinin (HA) via irreversible refolding of a trimeric conformational state leading to exposure of hydrophobic fusion peptides on each trimer subunit. Herein, we show that cells expressing fowl plague virus HA demonstrate discrete switching behavior with respect to the HA conformational change. Partially activated states do not exist at the scale of the cell, activation of HA leads to aggregation of cell surface trimers, and newly synthesized HA refold spontaneously in the presence of previously activated HA. These observations imply a feedback mechanism involving self-catalyzed refolding of HA and thus suggest a mechanism similar to the autocatalytic refolding and aggregation of prions.

Modeling a spin-labeled fusion peptide in a membrane: implications for the interpretation of EPR experiments

Site-directed spin-labeling and electron paramagnetic resonance are powerful tools for studying structure and conformational dynamics of proteins, especially in membranes. The position of the spin label is used as an indicator of the position of the site to which it is attached. The interpretation of these experiments is based on the assumptions that the spin label does not affect the peptide configuration and that it has a fixed orientation and distance with respect to the protein backbone. Here, the validity of these assumptions is examined through implicit membrane molecular dynamics simulations of the influenza hemagglutinin fusion peptide that has been labeled with methanethiosulfonate spin label. We find that the methanethiosulfonate spin label can occasionally induce peptide orientations that differ from those adopted by the wild-type peptide. Furthermore, the spin-label resides, on average, several Angstroms deeper in the membrane than the corresponding backbone C(alpha)-atom even at sites pointing toward the solvent. The nitroxide spin label exhibits flexibility and adopts various configurations depending on the surrounding residues.

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in micelles and bilayers: conformational analysis of peptide and lipids

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in two differently sized dodecylphosphocholine micelles and a palmitoyl oleoyl phosphatidylcholine bilayer were generated to analyze the influence of the environment. Four independent trajectories (5 ns each for the bilayer, and 2 ns each for the micelles) were generated for each system. The peptide lies at the surface of the micelles, while its N-terminal region inserts deeply in the bilayer. This leads to a substantial increase of the solvation and rigidity of the peptide in micelles as compared to the bilayer. The average structures, nevertheless, are similar in all three systems and agree reasonably with micelle-based NMR structures. When in the bilayer, the peptide increases the chain gauche population and area of adjacent lipids in the same binding leaflet, while it has the opposite effect for the nearby lipids of the other leaflet. These changes, which occur spontaneously to fill voids and defects, cause a decrease in the thickness of the membrane in the neighborhood of the peptide. They would be expected to promote positive curvature, as consistent with the formation of the convex bulge, or "nipple", in the initial stage of membrane fusion. An extension of the classical surfactant theory of Israelachvili based on shapes is proposed to introduce the concept of a "dynamically induced shape" of the membrane lipids by the peptide.

Effect of Lipid Composition on the “Membrane Response” Induced by a Fusion Peptide

To understand the initial stages of membrane destabilization induced by viral proteins, the factors important for binding of fusion peptides to cell membranes must be identified. In this study, effects of lipid composition on the mode of peptides' binding to membranes are explored via molecular dynamics (MD) simulations of the peptide E5, a water-soluble analogue of influenza hemagglutinin fusion peptide, in two full-atom hydrated lipid bilayers composed of dimyristoyl- and dipalmitoylphosphatidylcholine (DMPC and DPPC, respectively). The results show that, although the peptide has a common folding motif in both systems, it possesses different modes of binding. The peptide inserts obliquely into the DMPC membrane mainly with its N-terminal alpha helix, while in DPPC, the helix lies on the lipid/water interface, almost parallel to the membrane surface. The peptide seriously affects structural and dynamical parameters of surrounding lipids. Thus, it induces local thinning of both bilayers and disordering of acyl chains of lipids in close proximity to the binding site. The "membrane response" significantly depends upon lipid composition: distortions of DMPC bilayer are more pronounced than those in DPPC. Implications of the observed effects to molecular events on initial stages of membrane destabilization induced by fusion peptides are discussed.

Plasticity of influenza haemagglutinin fusion peptides and their interaction with lipid bilayers

A detailed molecular dynamics study of the haemagglutinin fusion peptide (N-terminal 20 residues of the HA2 subunits) in a model bilayer has yielded useful information about the molecular interactions leading to insertion into the lipids. Simulations were performed on the native sequence, as well as a number of mutant sequences, which are either fusogenic or nonfusogenic. For the native sequence and fusogenic mutants, the N-terminal 11 residues of the fusion peptides are helical and insert with a tilt angle of approximately 30 degrees with respect to the membrane normal, in very good agreement with experimental data. The tilted insertion of the native sequence peptide leads to membrane bilayer thinning and the calculated order parameters show larger disorder of the alkyl chains. These results indicate that the lipid packing is perturbed by the fusion peptide and could be used to explain membrane fusion. For the nonfusogenic sequences investigated, it was found that most of them equilibrate parallel to the interface plane and do not adopt a tilted conformation. The presence of a charged residue at the beginning of the sequence (G1E mutant) resulted in a more difficult case, and the outcomes do not fall straightforwardly into the general picture. Sequence searches have revealed similarities of the fusion peptide of influenza haemagglutinin with peptide sequences such as segments of porin, amyloid alpha eta peptide, and a peptide from the prion sequence. These results confirm that the sequence can adopt different folds in different environments. The plasticity and the conformational dependence on the local environment could be used to better understand the function of fusion peptides.