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

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.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

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.

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.

The higher level of complexity of K-Ras4B activation at the membrane

Is nucleotide exchange sufficient to activate K-Ras4B? To signal, oncogenic rat sarcoma (Ras) anchors in the membrane and recruits effectors by exposing its effector lobe. With the use of NMR and molecular dynamics (MD) simulations, we observed that in solution, farnesylated guanosine 5'-diphosphate (GDP)-bound K-Ras4B is predominantly autoinhibited by its hypervariable region (HVR), whereas the GTP-bound state favors an activated, HVR-released state. On the anionic membrane, the catalytic domain adopts multiple orientations, including parallel (∼180°) and perpendicular (∼90°) alignments of the allosteric helices, with respect to the membrane surface direction. In the autoinhibited state, the HVR is sandwiched between the effector lobe and the membrane; in the active state, with membrane-anchored farnesyl and unrestrained HVR, the catalytic domain fluctuates reinlessly, exposing its effector-binding site. Dimerization and clustering can reduce the fluctuations. This achieves preorganized, productive conformations. Notably, we also observe HVR-autoinhibited K-Ras4B-GTP states, with GDP-bound-like orientations of the helices. Thus, we propose that the GDP/GTP exchange may not be sufficient for activation; instead, our results suggest that the GDP/GTP exchange, HVR sequestration, farnesyl insertion, and orientation/localization of the catalytic domain at the membrane conjointly determine the active or inactive state of K-Ras4B. Importantly, K-Ras4B-GTP can exist in active and inactive states; on its own, GTP binding may not compel K-Ras4B activation.-Jang, H., Banerjee, A., Chavan, T. S, Lu, S., Zhang, J., Gaponenko, V., Nussinov, R. The higher level of complexity of K-Ras4B activation at the membrane.

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.

Mechanisms of membrane binding of small GTPase K-Ras4B farnesylated hypervariable region

K-Ras4B belongs to a family of small GTPases that regulates cell growth, differentiation and survival. K-ras is frequently mutated in cancer. K-Ras4B association with the plasma membrane through its farnesylated and positively charged C-terminal hypervariable region (HVR) is critical to its oncogenic function. However, the structural mechanisms of membrane association are not fully understood. Here, using confocal microscopy, surface plasmon resonance, and molecular dynamics simulations, we observed that K-Ras4B can be distributed in rigid and loosely packed membrane domains. Its membrane binding domain interaction with phospholipids is driven by membrane fluidity. The farnesyl group spontaneously inserts into the disordered lipid microdomains, whereas the rigid microdomains restrict the farnesyl group penetration. We speculate that the resulting farnesyl protrusion toward the cell interior allows oligomerization of the K-Ras4B membrane binding domain in rigid microdomains. Unlike other Ras isoforms, K-Ras4B HVR contains a single farnesyl modification and positively charged polylysine sequence. The high positive charge not only modulates specific HVR binding to anionic phospholipids but farnesyl membrane orientation. Phosphorylation of Ser-181 prohibits spontaneous farnesyl membrane insertion. The mechanism illuminates the roles of HVR modifications in K-Ras4B targeting microdomains of the plasma membrane and suggests an additional function for HVR in regulation of Ras signaling.

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.

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.

Detection of closed influenza virus hemagglutinin fusion peptide structures in membranes by backbone (13)CO- (15)N rotational-echo double-resonance solid-state NMR

The influenza virus fusion peptide is the N-terminal ~20 residues of the HA2 subunit of the hemagglutinin protein and this peptide plays a key role in the fusion of the viral and endosomal membranes during initial infection of a cell. The fusion peptide adopts N-helix/turn/C-helix structure in both detergent and membranes with reports of both open and closed interhelical topologies. In the present study, backbone (13)CO-(15)N REDOR solid-state NMR was applied to the membrane-associated fusion peptide to detect the distribution of interhelical distances. The data clearly showed a large fraction of closed and semi-closed topologies and were best-fitted to a mixture of two structures that do not exchange. One of the earlier open structural models may have incorrect G13 dihedral angles derived from TALOS analysis of experimentally correct (13)C shifts.

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.

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.

Models of toxic beta-sheet channels of protegrin-1 suggest a common subunit organization motif shared with toxic alzheimer beta-amyloid ion channels

Antimicrobial peptides (AMPs) induce cytotoxicity by altering membrane permeability. The electrical properties of membrane-associated AMPs as well as their cellular effects have been extensively documented; however their three-dimensional structure is poorly understood. Gaining insight into channel structures is important to the understanding of the protegrin-1 (PG-1) and other AMP cytolytic mechanisms, and to antibiotics design. We studied the beta-sheet channels morphology using molecular dynamics simulations. We modeled PG-1 channels as intrinsic barrel-stave and toroidal membrane pores, and simulated them in zwitterionic and anionic lipid bilayers. PG-1 channels consist of eight beta-hairpins in a consecutive NCCN (N and C represent the beta-hairpin's N- and C-termini) packing organization yielding antiparallel and parallel beta-sheet channels. Both channels preserve the toroidal, but not the barrel-stave pores. The two lipid leaflets of the bilayer bend toward each other at the channels' edges, producing a semitoroidal pore with the outward-pointing hydrophobic residues preventing the polar lipid headgroups from moving to the bilayer center. In all simulated lipid environments, PG-1 channels divide into four or five beta-sheet subunits consisting of single or dimeric beta-hairpins. The channel morphology with subunit organization is consistent with the four to five subunits observed by NMR in the POPE/POPG bilayer. Remarkably, a beta-sheet subunit channel motif is in agreement with Alzheimer ion channels modeled using the universal U-shape beta-strand-turn-beta-strand structure, as well as with high resolution atomic force microscopy images of beta-amyloid channels with four to six subunits. Consistent with the toxic beta-amyloid channels that are ion-conducting, the PG-1 channels permeate anions.