Characterizing a Histidine Switch Controlling pH-Dependent Conformational Changes of the Influenza Virus Hemagglutinin

During the fusion of the influenza virus to the host cell, bending of the HA2 chain of hemagglutinin into a hairpin-shaped structure in a pH-dependent manner facilitates the fusion of the viral envelope and the endosomal membrane. To characterize the structural and dynamical responses of the hinge region of HA2 to pH changes and examine the role of a conserved histidine in this region (the hinge histidine), we have performed an extensive set of molecular dynamics (MD) simulations of 26-residue peptides encompassing the hinge regions of several hemagglutinin subtypes under both neutral and low pH conditions, modeled by the change of the protonation state of the hinge histidine. More than 70 sets of MD simulations (collectively amounting to 25.1 μs) were performed in both implicit and explicit solvents to study the effect of histidine protonation on structural dynamics of the hinge region. In both explicit and implicit solvent simulations, hinge bending was consistently observed upon the protonation of the histidine in all the simulations starting with an initial straight helical conformation, whereas the systems with a neutral histidine retained their primarily straight conformation throughout the simulations. Conversely, the MD simulations starting from an initially bent conformation resulted in the formation of a straight helical structure upon the neutralization of the hinge histidine, whereas the bent structure was maintained when the hinge histidine remained protonated. Finally, mutation of the hinge histidine to alanine abolishes the bending response of the peptide altogether. A molecular mechanism based on the interaction of the hinge histidine with neighboring acidic residues is proposed to be responsible for its role in controlling the conformation of the hinge. We propose that this might present a common mechanism for pH-controlled structural changes in helical structures when histidines act as the pH sensor.     

Early steps of the conformational change of influenza virus hemagglutinin to a fusion active state: stability and energetics of the hemagglutinin

A conformational change of the homotrimeric glycoprotein hemagglutinin (HA) of influenza virus mediates fusion between the viral envelope and the endosome membrane. The conformational change of the HA ectodomain is triggered by the acidic pH of the endosome lumen. An essential step of the conformational change is the formation of an extended coiled-coil motif exposing the hydrophobic fusion peptide toward the target membrane. The structures of the neutral-pH, non-fusion active conformation of the HA ectodomain and of a fragment of the ectodomain containing the coiled-coil motif are known. However, it is not known by which mechanism protonation triggers the conformational change of the stable neutral-pH conformation of the ectodomain. Here, recent studies on the stability of the HA ectodomain at neutral pH, the energetics of the conformational change toward the fusion-active state and of the unfolding of the HA ectodomain are summarised. A model for the early steps of the conformational change of the HA ectodomain is presented. The model implicates that protonation leads to a partial dissociation of the distal domains of the HA monomers that is driven by electrostatic repulsion. The opening of the ectodomain enables water to enter the ectodomain. The interaction of water with respective sequences originally shielded from contact with water drives the formation of the coiled-coil structure.     

Early steps of the conformational change of influenza virus hemagglutinin to a fusion active state: Stability and energetics of the hemagglutinin

A conformational change of the homotrimeric glycoprotein hemagglutinin (HA) of influenza virus mediates fusion between the viral envelope and the endosome membrane. The conformational change of the HA ectodomain is triggered by the acidic pH of the endosome lumen. An essential step of the conformational change is the formation of an extended coiled-coil motif exposing the hydrophobic fusion peptide toward the target membrane. The structures of the neutral-pH, non-fusion active conformation of the HA ectodomain and of a fragment of the ectodomain containing the coiled-coil motif are known. However, it is not known by which mechanism protonation triggers the conformational change of the stable neutral-pH conformation of the ectodomain. Here, recent studies on the stability of the HA ectodomain at neutral pH, the energetics of the conformational change toward the fusion-active state and of the unfolding of the HA ectodomain are summarised. A model for the early steps of the conformational change of the HA ectodomain is presented. The model implicates that protonation leads to a partial dissociation of the distal domains of the HA monomers that is driven by electrostatic repulsion. The opening of the ectodomain enables water to enter the ectodomain. The interaction of water with respective sequences originally shielded from contact with water drives the formation of the coiled-coil structure.     

Helix-coil transition of PrP106-126: molecular dynamic study.

A set of 34 molecular dynamic (MD) simulations totaling 305 ns of simulation time of the prion protein-derived peptide PrP106-126 was performed with both explicit and implicit solvent models. The objective of these simulations is to investigate the relative stability of the alpha-helical conformation of the peptide and the mechanism for conversion from the helix to a random-coil structure. At neutral pH, the wild-type peptide was found to lose its initial helical structure very fast, within a few nanoseconds (ns) from the beginning of the simulations. The helix breaks up in the middle and then unwinds to the termini. The spontaneous transition into the random coil structure is governed by the hydrophobic interaction between His(111) and Val(122). The A117V mutation, which is linked to GSS disease, was found to destabilize the helix conformation of the peptide significantly, leading to a complete loss of helicity approximately 1 ns faster than in the wild-type. Furthermore, the A117V mutant exhibits a different mechanism for helix-coil conversion, wherein the helix begins to break up at the C-terminus and then gradually to unwind towards the N-terminus. In most simulations, the mutation was found to speed up the conversion through an additional hydrophobic interaction between Met(112) and the mutated residue Val(117), an interaction that did not exist in the wild-type peptide. Finally, the beta-sheet conformation of the wild-type peptide was found to be less stable at acidic pH due to a destabilization of the His(111)-Val(122), since at acidic pH this histidine is protonated and is unlikely to participate in hydrophobic interaction.     

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.     

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 α-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 °. 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.     

The 1-127 HA2 construct of influenza virus hemagglutinin induces cell-cell hemifusion.

Conformational changes in the HA2 subunit of influenza hemagglutinin (HA) are coupled to membrane fusion. We investigated the fusogenic activity of the polypeptide FHA2 representing 127 amino-terminal residues of the ectodomain of HA2. While the conformation of FHA2 both at neutral and at low pH is nearly identical to the final low-pH conformation of HA2, FHA2 still induces lipid mixing between liposomes in a low-pH-dependent manner. Here, we found that FHA2 induces lipid mixing between bound cells, indicating that the "spring-loaded" energy is not required for FHA2-mediated membrane merger. Although, unlike HA, FHA2 did not form an expanding fusion pore, both acidic pH and membrane concentrations of FHA2, required for lipid mixing, have been close to those required for HA-mediated fusion. Similar to what is observed for HA, FHA2-induced lipid mixing was reversibly blocked by lysophosphatidylcholine and low temperature, 4 degrees C. The same genetic modification of the fusion peptide inhibits both HA- and FHA2-fusogenic activities. The kink region of FHA2, critical for FHA2-mediated lipid mixing, was exposed in the low-pH conformation of the whole HA prior to fusion. The ability of FHA2 to mediate lipid mixing very similar to HA-mediated lipid mixing is consistent with the hypothesis that hemifusion requires just a portion of the energy released in the conformational change of HA at acidic pH.     

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.     

pH-dependent membrane fusion activity of a synthetic twenty amino acid peptide with the same sequence as that of the hydrophobic segment of influenza virus hemagglutinin

A twenty amino acid hydrophobic peptide with the same sequence as that of the HA2 N-terminal segment of influenza virus hemagglutinin was synthesized and studied as to its fusion activity. The peptide caused rapid and efficient fusion of egg yolk phosphatidylcholine sonicated vesicles at acidic pH but not at neutral pH. The threshold pH was ca. 6.2 and the maximum fusion occurred at pH 4.8, the half-maximal pH for fusion being 5.6. The pH dependence was similar to that of the parent virus. The fusion efficiency was dependent on the ration of lipid to peptide, increasing with decreasing ratio. The fusion can be rapidly switched on and off by adjusting the pH, to the acidic side and neutral, respectively. The peptide with an acetylated or succinylated N-terminus also showed low pH-induced fusion activity but the pH range was shifted by ca. 1 unit to the acidic side. The results indicate that the HA2 hydrophobic segment in the virus fusion protein is directly involved in the fusion reaction and protonation of the acidic residues in the segment is required for the activity.     

Studies on the conformation of protonated DNA

Experiments were made to demonstrate the predominant protonation effects and structural changes of the ordered double helical DNA structure and denatured state of DNA. Spectrophotometric titrations performed at different wavelengths indicate that cytosine can be protonated in the DNA double helical molecule to a high extent without breakdown of the secondary structure. With DNA heat-denatured under severe conditions the protonation of cytosine can be measured at 280, 295, and 300 mμ: the apparent pK value obtained was ～4.6. The protonated double helical conformation of the DNA molecule differs from the unprotonated state, which follows from the decrease of the thermal stability and from changes in the ORD curves. The ORD of a GC-rich DNA indicates a novel Cotton effect with positive rotations at ～260 mμ in 0.02M KCl below pH 4.0 to pH 3.3. The occurrence of the new peak parallels the extent of protonated cytosine measured by the spectrophotometric titrations. It is concluded that the protonated cytosine in the double helical structure is responsible for the difference between the protonated DNA conformation and the native state at neutral pH.     

An Induced Pocket for the Binding of Potent Fusion Inhibitor CL-385319 with H5N1 Influenza Virus Hemagglutinin

The influenza glycoprotein hemagglutinin (HA) plays crucial roles in the early stage of virus infection, including receptor binding and membrane fusion. Therefore, HA is a potential target for developing anti-influenza drugs. Recently, we characterized a novel inhibitor of highly pathogenic H5N1 influenza virus, CL-385319, which specifically inhibits HA-mediated viral entry. Studies presented here identified the critical binding residues for CL-385319, which clustered in the stem region of the HA trimer by site-directed mutagenesis. Extensive computational simulations, including molecular docking, molecular dynamics simulations, molecular mechanics generalized Born surface area (MM_GBSA) calculations, charge density and Laplacian calculations, have been carried out to uncover the detailed molecular mechanism that underlies the binding of CL-385319 to H5N1 influenza virus HA. It was found that the recognition and binding of CL-385319 to HA proceeds by a process of "induced fit" whereby the binding pocket is formed during their interaction. Occupation of this pocket by CL-385319 stabilizes the neutral pH structure of hemagglutinin, thus inhibiting the conformational rearrangements required for membrane fusion. This "induced fit" pocket may be a target for structure-based design of more potent influenza fusion inhibitors.     

Structural intermediates in the low pH-induced transition of influenza hemagglutinin

The hemagglutinin (HA) glycoproteins of influenza viruses play a key role in binding host cell receptors and in mediating virus-host cell membrane fusion during virus infection. Upon virus entry, HA is triggered by low pH and undergoes large structural rearrangements from a prefusion state to a postfusion state. While structures of prefusion state and postfusion state of HA have been reported, the intermediate structures remain elusive. Here, we report two distinct low pH intermediate conformations of the influenza virus HA using cryo-electron microscopy (cryo-EM). Our results show that a decrease in pH from 7.8 to 5.2 triggers the release of fusion peptides from the binding pockets and then causes a dramatic conformational change in the central helices, in which the membrane-proximal ends of the central helices unwind to an extended form. Accompanying the conformational changes of the central helices, the stem region of the HA undergoes an anticlockwise rotation of 9.5 degrees and a shift of 15 Å. The HA head, after being stabilized by an antibody, remains unchanged compared to the neutral pH state. Thus, the conformational change of the HA stem region observed in our research is likely to be independent of the HA head. These results provide new insights into the structural transition of HA during virus entry.     

Conformation and interaction with the membrane models of the amino-terminal peptide of influenza virus hemagglutinin HA2 at fusion pH.

Conformations of a 48-mer peptide corresponding to the amino-terminal region of influenza HA2 in aqueous and membranous environments were studied. In aqueous solution the peptide was found to be oligomeric and its helicity was enhanced at higher concentrations. The conformation in phospholipid bilayer and insertion depth into the sodium dodecyl sulfate (SDS) micelle for the fusion peptide were in line with those determined for the amino-terminal 25-mer analog. The turn of residues 28-31 found in the crystal structure of hemagglutinin at neutral pH persisted in the presence of SDS at pH 5.0. Except for the turn, conformational lability of the amino portion of HA2 is suggested by comparison of the secondary structure determined herein with that obtained with the influenza fusion protein crystallized in the aqueous phase at neutral pH. The backbone amide proton exchange experiment suggested an interaction with the micellar surface for the segment carboxy-terminal to the fusion peptide domain.     

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.     

The intrinsic helical propensities of the helical fragments in prion protein under neutral and low pH conditions: a replica exchange molecular dynamics study

Replica exchange molecular dynamics simulations in neutral and acidic aqueous solutions were employed to study the intrinsic helical propensities of three helices in both Syrian hamster (syPrP) and human (huPrP) prion proteins. The helical propensities of syPrP HA and huPrP HA are very high under both pH conditions, which implies that HA is barely involved in the helix-to-β transition. The SyPrP HB chain has a strong tendency to adopt an extended conformation, which is possibly involved in the mechanism of infectious prion diseases in Syrian hamster. HuPrP HC has more of a preference for the extended conformation than huPrP HA and huPrP HB do, which leads to the conjecture that it is more likely to be the source of β-rich structure for human prion protein. We also noticed that the presence of salt bridges is not correlated with helical propensity, indicating that salt bridges do not stabilize helices.     

Conformational aspects of the acid-induced fusion mechanism of influenza virus hemagglutinin. Circular dichroism and fluorescence studies

Circular dichroism and tryptophan fluorescence spectroscopy have been used to investigate the structures of the influenza virus membrane glycoprotein hemagglutinin, acid-treated hemagglutinin, and fragments of hemagglutinin derived by proteolysis. The conformational change in hemagglutinin which occurs at the pH of membrane fusion (pH 5-6) was associated with a significant change of the environment of tyrosine residues, a change in the environment of tryptophan residues, but no changes in secondary structure. Tryptic digestion of the hemagglutinin in its low pH conformation which releases one of the subunit polypeptides (HA1) caused minimal changes in tyrosine and tryptophan environments but a small secondary structural change in HA1. The secondary structure of the remainder of the molecule (HA2) was very similar to that predicted from the known x-ray crystallographic structure of the native molecule. However, fluorescence spectroscopy indicated a tertiary change in structure in the coiled coil of alpha-helices which form the fibrous central stem of the molecule. These results are consistent with a conformational change required for membrane fusion which involves a decrease of HA1/HA1, HA1/HA2 interactions and changes in tertiary structure not accompanied by changes in secondary structure.     

The pH dependence of flavivirus envelope protein structure: insights from molecular dynamics simulations

The flavivirus membrane fusion is triggered by the acid pH of the endosomes after virus endocytosis. The proposed mechanism involves changes in the protonation state of conserved histidine residues of the E protein present in the viral surface that undergoes a series of structural rearrangements that result in the fusion between the endosome and viral bilayers. We studied the pH dependence of E protein rearrangements of dengue virus type 2, used as a model, in the pH range experimented by the virus along the fusion process. We employed a low computational cost scheme to explore the behavior of the E protein by molecular dynamics (MD) simulations of complete systems that include the protein, the solvent, and ions. The procedure alternates cyclically the update of the ionization states of the protein residues with common MD steps applied to the new ionization configuration. Important pH-dependent protein structure rearrangements consistent with the changes of the protonation states of conserved histidine residues were observed. The involvement of other conserved residues in the flavivirus in the rearrangements was also identified. The results show interesting correlations with a proposed model for the fusion mechanism, as well as the experimentally identified key residues, contributing to a better understanding of the structural changes in protein E that lead to the fusion process.     

Molecular dynamics simulations of Ac-3Aib-Cage-3Aib-NHMe

Solvents play a stabilising role with the more stable conformations obtained in polar solvents than in vacuo. We investigate to what extent the structural propensities of the pentacyclo-undecane (PCU) cage polypeptide chain of the type Ac-3Aib-Cage-3Aib-NHMe are influenced in implicit water and in explicit solvents: methanol (MEOH), dimethyl sulphoxide (DMSO) and TIP3P water. The sampling of the α-helical conformations of the PCU cage polypeptide was investigated using the in-house modified PARM94 force-field parameters. Analysis of 50 ns molecular dynamics (MD) simulations revealed a tendency of the PCU cage polypeptide to assume bent structures, especially in polar solvents. The choice of solvents was designed to relate the simulations to physiological conditions. The individual amino-isobutyric acid residues predominantly sampled the right-handed and left-handed 3₁₀-helical conformations, indicating that the helical conformations are preferred in all four environments (in vacuo, MEOH, water and DMSO). Additionally, the 100 ns replica exchange MD (REMD) simulations of the PCU cage polypeptide in implicit water revealed more conformational variety present than in explicit solvents, and is more consistent with previous theoretical studies on the PCU cage residue. The present theoretical results may help in rationalising experimental results on these PCU cage polypeptides, and definitely show the importance of a dynamical approach for a correct interpretation and prediction of the conformational behaviour of the PCU cage molecules in different environments.     

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.