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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Inhibition of Influenza H7 Hemagglutinin-Mediated Entry

The recent outbreak of H7N9 influenza in China is of high concern to public health. H7 hemagglutinin (HA) plays a critical role in influenza entry and thus HA presents an attractive target for antivirals. Previous studies have suggested that the small molecule tert-butyl hydroquinone (TBHQ) inhibits the entry of influenza H3 HA by binding to the stem loop of HA and stabilizing the neutral pH conformation of HA, thereby disrupting the membrane fusion step. Based on amino acid sequence, structure and immunogenicity, H7 is a related Group 2 HA. In this work we show, using a pseudovirus entry assay, that TBHQ inhibits H7 HA-mediated entry, as well as H3 HA-mediated entry, with an IC50∼6 µM. Using NMR, we show that TBHQ binds to the H7 stem loop region. STD NMR experiments indicate that the aromatic ring of TBHQ makes extensive contact with the H7 HA surface. Limited proteolysis experiments indicate that TBHQ inhibits influenza entry by stabilizing the H7 HA neutral pH conformation. Together, this work suggests that the stem loop region of H7 HA is an attractive target for therapeutic intervention and that TBHQ, which is a widely used food preservative, is a promising lead compound.     

Hydrophobic Effect Drives Oxygen Uptake in Myoglobin via Histidine E7

Since the elucidation of the myoglobin (Mb) structure, a histidine residue on the E helix (His-E7) has been proposed to act as a gate with an open or closed conformation controlling access to the active site. Although it is believed that at low pH, the His-E7 gate is in its open conformation, the full relationship between the His-E7 protonation state, its conformation, and ligand migration in Mb is hotly debated. We used molecular dynamics simulations to first address the effect of His-E7 protonation on its conformation. We observed the expected shift from the closed to the open conformation upon protonation, but more importantly, noted a significant difference between the conformations of the two neutral histidine tautomers. We further computed free energy profiles for oxygen migration in each of the possible His-E7 states as well as in two instructive Mb mutants: Ala-E7 and Trp-E7. Our results show that even in the closed conformation, the His-E7 gate does not create a large barrier to oxygen migration and permits oxygen entry with only a small rotation of the imidazole side chain and movement of the E helix. We identify, instead, a hydrophobic site in the E7 channel that can accommodate an apolar diatomic ligand and enhances ligand uptake particularly in the open His-E7 conformation. This rate enhancement is diminished in the closed conformation. Taken together, our results provide a new conceptual framework for the histidine gate hypothesis.     

High-resolution structural profile of hylaseptin-4: Aggregation,membrane topology and pH dependence of overall membrane binding process

Hylaseptin-4 (HSP-4, GIGDILKNLAKAAGKAALHAVGESL-NH₂) is an antimicrobial peptide originally isolated from Hypsiboas punctatus tree frog. The peptide has been chemically synthetized for structural investigations by CD and NMR spectroscopies. CD experiments reveal the high helical content of HSP-4 in biomimetic media. Interestingly, the aggregation process seems to occur at high peptide concentrations either in aqueous solution or in presence of biomimetic membranes, indicating an increase in the propensity of the peptide for adopting a helical conformation. High-resolution NMR structures determined in presence of DPC-d₃₈ micelles show a highly ordered α-helix from amino acid residues I2 to S24 and a smooth bend near G14. A large separation between hydrophobic and hydrophilic residues occurs up to the A16 residue, from which a shift in the amphipathicity is noticed. Oriented solid-state NMR spectroscopy show a roughly parallel orientation of the helical structure along the POPC lipid bilayer surface, with an insertion of the hydrophobic N-terminus into the bilayer core. Moreover, a noticeable pH dependence of the aggregation process in both aqueous and in biomimetic membrane environments is attributed to a single histidine residue (H19). The protonation degree of the imidazole side-chain might help in modulating the peptide-peptide or peptide-lipid interactions. Finally, molecular dynamics simulations confirm the orientation and preferential helical conformation and in addition, show that HSP-4 tends to self-aggregate in order to stabilize its active conformation in aqueous or phospholipid bilayer environments.     

How hormone receptor-DNA binding affects nucleosomal DNA: the role of symmetry.

Molecular dynamics simulations have been employed to determine the optimal conformation of an estrogen receptor DNA binding domain dimer bound to a consensus response element, ds(AGGTCACAGTGACCT), and to a nonconsensus response element, ds(AGAACACAGTGACCT). The structures simulated were derived from a crystallographic structure and solvated by a sphere (45-A radius) of explicit water and counterions. Long-range electrostatic interactions were accounted for during 100-ps simulations by means of a fast multipole expansion algorithm combined with a multiple time-step scheme in the molecular dynamics package NAMD. The simulations demonstrate that the dimer induces a bent and underwound (10.7 bp/turn) conformation in the DNA. The bending reflects the dyad symmetry of the receptor dimer and can be described as an S-shaped curve in the helical axis of DNA when projected onto a plane. A similar bent and underwound conformation is observed for nucleosomal DNA near the nucleosome's dyad axis that reflects the symmetry of the histone octamer. We propose that when a receptor dimer binds to a nucleosome, the most favorable dimer-DNA and histone-DNA interactions are achieved if the respective symmetry axes are aligned. Such positioning of a receptor dimer over the dyad of nucleosome B in the mouse mammary tumor virus promoter is in agreement with experiment.     

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.     

Structural characterizations of fusion peptide analogs of influenza virus hemagglutinin. Implication of the necessity of a helix-hinge-helix motif in fusion activity

Infection by enveloped viruses initially involves membrane fusion between viral and host cell membranes. The fusion peptide plays a crucial role in triggering this reaction. To clarify how the fusion peptide exerts this specific function, we carried out biophysical studies of three fusion peptide analogs of influenza virus hemagglutinin HA2, namely E5, G13L, and L17A. E5 exhibits an activity similar to the native fusion peptide, whereas G13L and L17A, which are two point mutants of the E5 analog, possess much less fusion activity. Our CD data showed that the conformations of these three analogs in SDS micelles are pH-dependent, with higher alpha-helical contents at acidic pH. Tryptophan fluorescence emission experiments indicated that these three analogs insert deeper into lipid bilayers at acidic pH. The three-dimensional structure of the E5 analog in SDS micelles at pH 4.0 revealed that two segments, Leu(2)-Glu(11) and Trp(14)-Ile(18), form amphipathic helical conformations, with Gly(12)-Gly(13) forming a hinge. The hydrophobic residues in the N- and C-terminal helices form a hydrophobic cluster. At neutral pH, however, the C-terminal helix of Trp(14)-Ile(18) reduces dramatically, and the hydrophobic core observed at acidic pH is severely disrupted. We suggest that the disruption of the C-terminal helix renders the E5 analog fusion-inactive at neutral pH. Furthermore, the decrease of the hinge and the reduction of fusion activity in G13L reveal the importance of the hinge in fusion activity. Also, the decrease in the C-terminal helix and the reduction of fusion activity in L17A demonstrates the importance of the C-terminal helix in fusion activity. Based on these biophysical studies, we propose a model that illustrates the structural change of the HA2 fusion peptide analog and explains how the analog interacts with the lipid bilayer at different pH values.     

Conformation of the Ras-binding domain of Raf studied by molecular dynamics and free energy simulations

Recognition of Ras by its downstream target Raf is mediated by a Ras-recognition region in the Ras-binding domain (RBD) of Raf. Residues 78–89 in this region occupy two different conformations in the ensemble of NMR solution structures of the RBD: a fully α-helical one, and one where 87–90 form a type IV β-turn. Molecular dynamics simulations of the RBD in solution were performed to explore the stability of these and other possible conformations of both the wild-type RBD and the R89K mutant, which does not bind Ras. The simulations sample a fully helical conformation for residues 78–89 similar to the NMR helical structures, a conformation where 85–89 form a 3₁₀-helical turn, and a conformation where 87–90 form a type I |iB-turn, whose free energies are all within 0.3 kcal/mol of each other. NOE patterns and H^α chemical shifts from the simulations are in reasonable agreement with experiment. The NMR turn structure is calculated to be 3 kcal/mol higher than the three above conformations. In a simulation with the same implicit solvent model used in the NMR structure generation, the turn conformation relaxes into the fully helical conformation, illustrating possible structural artifacts introduced by the implicit solvent model. With the Raf R89K mutant, simulations sample a fully helical and a turn conformation, the turn being 0.9 kcal/mol more stable. Thus, the mutation affects the population of RBD conformations, and this is expected to affect Ras binding. For example, if the fully helical conformation of residues 78–89 is required for binding, its free energy increase in R89K will increase the binding free energy by about 0.6 kcal/mol. Proteins 31:186–200, 1998. © 1998 Wiley-Liss, Inc.     

A computational study of nucleosomal DNA flexibility

Molecular dynamics simulations of the nucleosome core particle and its isolated DNA free in solution are reported. The simulations are based on the implicit solvent methodology and provide insights into the nature of large-scale structural fluctuations and flexibility of the nucleosomal DNA. In addition to the kinked regions previously identified in the x-ray structure of the nucleosome, the simulations support the existence of a biochemically identified distorted region of the DNA. Comparison of computed relative free energies shows that formation of the kinks is associated with little, if any, energy cost relative to a smooth, ideal conformation of the DNA superhelix. Isolated nucleosomal DNA is found to be considerably more flexible than expected for a 147 bp stretch of DNA based on its canonical persistence length of 500 A. Notably, the significant bending of the DNA observed in our simulations occurs without breaking of Watson-Crick bonds. The computed relative stability of bent conformations is sensitive to the ionic strength of the solution in the physiological range; the sensitivity suggests possible experiments that might provide further insights into the structural origins of the unusual flexibility of the DNA.