Stress and strain in staphylococcal nuclease.

Protein molecules generally adopt a tertiary structure in which all backbone and side chain conformations are arranged in local energy minima; however, in several well-refined protein structures examples of locally strained geometries, such as cis peptide bonds, have been observed. Staphylococcal nuclease A contains a single cis peptide bond between residues Lys 116 and Pro 117 within a type VIa beta-turn. Alternative native folded forms of nuclease A have been detected by NMR spectroscopy and attributed to a mixture of cis and trans isomers at the Lys 116-Pro 117 peptide bond. Analyses of nuclease variants K116G and K116A by NMR spectroscopy and X-ray crystallography are reported herein. The structure of K116A is indistinguishable from that of nuclease A, including a cis 116-117 peptide bond (92% populated in solution). The overall fold of K116G is also indistinguishable from nuclease A except in the region of the substitution (residues 112-117), which contains a predominantly trans Gly 116-Pro 117 peptide bond (80% populated in solution). Both Lys and Ala would be prohibited from adopting the backbone conformation of Gly 116 due to steric clashes between the beta-carbon and the surrounding residues. One explanation for these results is that the position of the ends of the residue 112-117 loop only allow trans conformations where the local backbone interactions associated with the phi and psi torsion angles are strained. When the 116-117 peptide bond is cis, less strained backbone conformations are available. Thus the relaxation of the backbone strain intrinsic to the trans conformation compensates for the energetically unfavorable cis X-Pro peptide bond. With the removal of the side chain from residue 116 (K116G), the backbone strain of the trans conformation is reduced to the point that the conformation associated with the cis peptide bond is no longer favorable.     

Modeling protein loops using a phi i + 1, psi i dimer database.

We present an automated method for modeling backbones of protein loops. The method samples a database of phi i + 1 and psi i angles constructed from a nonredundant version of the Protein Data Bank (PDB). The dihedral angles phi i + 1 and psi i completely define the backbone conformation of a dimer when standard bond lengths, bond angles, and a trans planar peptide configuration are used. For the 400 possible dimers resulting from 20 natural amino acids, a list of allowed phi i + 1, psi i pairs for each dimer is created by pooling all such pairs from the loop segments of each protein in the nonredundant version of the PDB. Starting from the N-terminus of the loop sequence, conformations are generated by assigning randomly selected pairs of phi i + 1, psi i for each dimer from the respective pool using standard bond lengths, bond angles, and a trans peptide configuration. We use this database to simulate protein loops of lengths varying from 5 to 11 amino acids in five proteins of known three-dimensional structures. Typically, 10,000-50,000 models are simulated for each protein loop and are evaluated for stereochemical consistency. Depending on the length and sequence of a given loop, 50-80% of the models generated have no stereochemical strain in the backbone atoms. We demonstrate that, when simulated loops are extended to include flanking residues from homologous segments, only very few loops from an ensemble of sterically allowed conformations orient the flanking segments consistent with the protein topology. The presence of near-native backbone conformations for loops from five different proteins suggests the completeness of the dimeric database for use in modeling loops of homologous proteins. Here, we take advantage of this observation to design a method that filters near-native loop conformations from an ensemble of sterically allowed conformations. We demonstrate that our method eliminates the need for a loop-closure algorithm and hence allows for the use of topological constraints of the homologous proteins or disulfide constraints to filter near-native loop conformations.     

The 2.1A crystal structure of the far-red fluorescent protein HcRed: inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1A resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu(64)-Tyr(65)-Gly(66), was observed to adopt both a cis coplanar and a trans non-coplanar conformation. As a result of these two conformations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the beta-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Beta-sheet folding of fragment (16-36) of bovine pancreatic trypsin inhibitor as predicted by Monte Carlo simulated annealing.

A tertiary structure prediction is described using Monte Carlo simulated annealing for the peptide fragment corresponding to residues 16-36 of bovine pancreatic trypsin inhibitor (BPTI). The simulation starts with randomly chosen initial conformations and is performed without imposing experimental constraints using energy functions given for generic interatomic interactions. Out of 20 simulation trials, seven conformations show a sheet-like structure--two strands connected by a turn--although this sheet-like structure is not as rigid as that observed in native BPTI. It is also shown that these conformations are mostly looped and exhibit a native-like right-handed twist. Unlike the case with the C-peptide of RNase A, no conspicuous alpha-helical structure is found in any of the final conformations obtained in the simulation. However, the lowest-energy conformation does not resemble exactly the native structure. This indicates that the rigid beta-sheet conformation of native BPTI merely corresponds to a local minimum of the energy function if the fragment with residues 16-36 is isolated from the native protein. A statistical analysis of all 20 final conformations suggests that the tendency for the peptide segments to form extended beta-strands is strong for those with residues 18-24, and moderate for those with residues 30-35. The segment of residues 25-29 does not tend to form any definite structure. In native BPTI, the former segments are involved in the beta-sheet and the latter in the turn. A folding scenario is also speculated from this analysis.     

Analysis of side-chain conformational distributions in neutrophil peptide-5 NMR structures

The side-chain conformations have been analyzed in the antimicrobial peptide, Neutrophil Peptide-5 (NP-5), whose structure was independently generated from nmr-derived distance constraints using a distance geometry algorithm. The side-chain and peptide dihedral angle distributions in the nmr structures were compared with those constructed from a data base of high-resolution protein crystal structures. The side-chain conformational preferences for NP-5 in solution are significantly different from those observed in the crystal structure data base. These results indicate that the side-chain conformations are quite disordered for many of the residues of NP-5. The absence of a correlation between the width of the conformational distribution and surface accessibility suggests that the disorder may be due to limitations in the structural information extracted from the nmr data rather than to molecular motion. However, it is also observed that the degree of conformational disorder is only weakly correlated with the number of nuclear Overhauser enhancements to a given side chain. Possible reasons for this are discussed. Molecular mechanics refinement of these structures did not significantly change the side-chain populations. Anomolously wide distributions are observed for rotations about the peptide bonds and the disulfide bonds in the NP-5 distance geometry structures, which are improved by the refinement. The very high degree of order observed for the central dihedral angle of the disulfide bond in the high-resolution crystal data base suggests that the rotation about this bond in proteins is determined by the local potential.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Structural compromise of disallowed conformations in peptide and protein structures

Using a data set of 454 crystal structures of peptides and 80 crystal structures of non-homologous proteins solved at ultra high resolution of 1.2 A or better we have analyzed the occurrence of disallowed Ramachandran (phi, psi) angles. Out of 1492 and 13508 non-glycyl residues in peptides and proteins respectively 12 and 76 residues in the two datasets adopt clearly disallowed combinations of Ramachandran angles. These examples include a number of conformational points which are far away from any of the allowed regions in the Ramachandran map. According to the Ramachandran map a given (phi, psi) combination is considered disallowed when two non-bonded atoms in a system of two-linked peptide units with ideal geometry are prohibitively proximal in space. However, analysis of the disallowed conformations in peptide and protein structures reveals that none of the observations of disallowed conformations in the crystal structures correspond to a short contact between non-bonded atoms. A further analysis of deviations of bond lengths and angles, from the ideal peptide geometry, at the residue positions of disallowed conformations in the crystal structures suggest that individual bond lengths and angles are all within acceptable limits. Thus, it appears that the rare tolerance of disallowed conformations is possible by gentle and acceptable deviations in a number of bond lengths and angles, from ideal geometry, over a series of bonds resulting in a net gross effect of acceptable non-bonded inter-atomic distances.     

A new force-field program for the calculation of glycopeptides and its application to a heptacosapeptide-decasaccharide of immunoglobulin G1. Importance of 1-6-glycosidic linkages in carbohydrate.peptide interactions

Energetically favored conformations of glycopeptide 1 were calculated using the newly developed force-field program, GEGOP (geometry of glycopeptides). The three-dimensional structure of glycopeptide 1, which is part of the Fc fragment of IgG1, has been calculated. 1 contains 27 amino acid residues from Pro291 to Lys317 and a biantennary decasaccharide N-linked to Asn297. The conformations of the peptide and the carbohydrate parts are shown to be mutually dependent. Single glycosyl residues of 1 exhibit interaction energies of up to -31.8 kJ/mol with the peptide portion. Generally, only a few of the glycosyl residues of the oligosaccharide moiety express significant interaction energies with the peptide part. No easy prediction is possible of glycosyl residues which exhibit favorable interaction energies. However, in all of the calculated structures, the glycosyl residues of the 1-6-linked branches show strong attractive forces for the peptide part. 1-6-glycosidically linked branches can adopt a larger number of conformations than other linkages due to their high flexibility which allows more favorable interactions with proteins. We developed the GEGOP program in order to be able to study the preferred conformations of large glycopeptides. The program is based on the GESA (geometry of saccharides) program and utilizes the HSEA (hard sphere exo anomeric) force field for the carbohydrate part and the ECEPP/2 (empirical conformation energy program for peptides) force field [Némethy, G., Pottle, M. S. & Scheraga, H. A. (1983) J. Phys. Chem. 87, 1883-1887] for the peptide part. The GEGOP program allows the simultaneous relaxation of all rotational degrees of freedom of these glycoconjugates during the energy optimization process. Thus, mutual interactions between glycosyl and amino acid residues can be studied in detail.     

Solution structure of an LNA hybridized to DNA: NMR study of the d(CT(L)GCT(L)T(L)CT(L)GC):d(GCAGAAGCAG) duplex containing four locked nucleotides

We have used two-dimensional (1)H NMR spectroscopy at 750 MHz to determine a high-resolution solution structure of an oligonucleotide containing restricted nucleotides with a 2'-O, 4'-C-methylene bridge (LNA) hybridized to the complementary DNA strand. The LNA:DNA duplex examined contained four thymidine LNA modifications (T(L), d(C1T(L)2G3C4T(L)5T(L)6C7T(L)8G9C10):d( G11C12A13G14A15A16G17C 18A19G20). A total relaxation matrix approach was used to obtain interproton distance bounds from NOESY cross-peak intensities. These distance bounds were used as restraints in molecular dynamics (rMD) calculations. Forty final structures were generated for the duplex from A-form and B-form DNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the 40 structures of the complex was 0.6 A. The sugar puckerings are averaged values of a dynamic interchange between N- and S-type conformation except in case of the locked nucleotides that were found to be fixed in the C3'-endo conformation. Among the other nucleotides in the modified strand, the furanose ring of C7 and G9 is predominantly in the N-type conformation whereas that of G3 is in a mixed conformation. The furanose rings of the nucleotides in the unmodified complementary strand are almost exclusively in the S-type conformation. Due to these different conformations of the sugars in the two strands, there is a structural strain between the A-type modified strand and the B-type unmodified complementary strand. This strain is relaxed by decreasing the value of rise and compensating with tip, buckle, and propeller twist. The values of twist vary along the strand but for a majority of the base pairs a value even lower than that of A-DNA is observed. The average twist over the sequence is 32+/-1 degrees. On the basis of the structure, we conclude that the high stability of LNA:DNA duplexes is caused by a local change of the phosphate backbone geometry that favors a higher degree of stacking.     

Survey of conformational role of ester bonds in a cyclic depsipeptide. Study on cyclo(-L-Ala-L-Hmb-)2 by energy calculation and NMR spectroscopy.

The effect of ester bond on the conformation of peptide molecule was studied by designing and synthesizing a model tetradepsipeptide cyclo(-L-Ala-L-Hmb-)2 and by analyzing the conformation both theoretically and experimentally. Theoretical analysis showed that both ester and peptide bonds in the calculated low-energy conformations within 3 kcal/mol of the global minimum take a trans but distorted configuration. The distortion is larger in ester bonds than in peptide bonds. Further, the four carbonyls project from one side of the plane of the cyclic backbone, whereas the side chains project from the other side. These results are consistent with the experimental results obtained by NMR measurement; first, the coupling constant deduced from 1H-NMR species in DMSO-d6 is consistent with the dihedral angles of the calculated low-energy conformations; second, results of NOE measurement can reproduce the calculated configuration of the carbonyls and side chains. From the consistency between theoretical and experimental results, it is concluded that this model tetradepsipeptide takes an all-trans backbone conformation in solution and this backbone conformation is stabilized by large distortion in the ester bond, which compensates the strain resulted from the 12-membered cyclic backbone structure consisting only of L-residues.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

Conformation of Leu-Arg-Arg-Ala-Ser-Leu-Gly bound in the active site of adenosine cyclic 3',5'-phosphate dependent protein kinase

Studies utilizing NMR spectroscopy have shown that adenosine cyclic 3',5'-phosphate dependent protein kinase (A-kinase) probably binds Leu-Arg-Arg-Ala-Ser-Leu-Gly (peptide 1) in one of two extended coil conformations (A or B). The relative reactivities of a series of N-methylated peptides based on the structure of peptide 1 might, therefore, be related to how well each can assume the A or B conformation. From estimates of the magnitude of steric interactions that would be induced by N-methylation of an amide in peptide 1 that is locked in either conformation, the ability of each peptide to form that conformation was predicted. The ability of A-kinase to catalyze phosphorylation of the N-methylated peptides correlated well with the ability of each peptide to form conformation A, but not conformation B. In accord with these findings, the reactivity of an unreactive N-methylated peptide was partially restored by a second change, which allowed the peptide to assume conformation A. These results suggest that, when bound in the enzymatic active site, peptide 1 has a conformation that resembles structure A much more closely than structure B.     

Mechanism of action of aspartic proteases. III. Conformational characteristics of HIV-1 protease inhibitor JG-365]

A set of conformations was shown to be characteristic of the free-state spatial structure of substrate-like inhibitor JG-365 for aspartic protease from HIV-1. Among them, the lowest-energy conformations have a folded form of the peptide backbone. The inhibitor has a noncleavable hydroxyethylamine group with an additional chiral center in its structure. Our calculations showed that only the S-isomer of the inhibitor displays conformational characteristics that practically coincide with those of the native substrate for HIV-1 protease. One of the calculated conformations with a completely extended main chain and a relative energy of 9.5 kcal/mol very closely mimics the experimentally observed structure of the inhibitor in the enzyme-inhibitor complex. The realization of this structure is unlikely for a free inhibitor, because it has only a small number of interresidual noncovalent interactions in the extended conformation; these are presumably compensated for by intermolecular interactions at the active site of the enzyme.     

The determinants of bond angle variability in protein/peptide backbones: A comprehensive statistical/quantum mechanics analysis

The elucidation of the mutual influence between peptide bond geometry and local conformation has important implications for protein structure refinement, validation, and prediction. To gain insights into the structural determinants and the energetic contributions associated with protein/peptide backbone plasticity, we here report an extensive analysis of the variability of the peptide bond angles by combining statistical analyses of protein structures and quantum mechanics calculations on small model peptide systems. Our analyses demonstrate that all the backbone bond angles strongly depend on the peptide conformation and unveil the existence of regular trends as function of ψ and/or φ. The excellent agreement of the quantum mechanics calculations with the statistical surveys of protein structures validates the computational scheme here employed and demonstrates that the valence geometry of protein/peptide backbone is primarily dictated by local interactions. Notably, for the first time we show that the position of the H^α hydrogen atom, which is an important parameter in NMR structural studies, is also dependent on the local conformation. Most of the trends observed may be satisfactorily explained by invoking steric repulsive interactions; in some specific cases the valence bond variability is also influenced by hydrogen‐bond like interactions. Moreover, we can provide a reliable estimate of the energies involved in the interplay between geometry and conformations. Proteins 2015; 83:1973–1986. © 2015 Wiley Periodicals, Inc.     

Conformation in solution of porcine brain natriuretic peptide determined by combined use of nuclear magnetic resonance and distance geometry

The conformation in solution of porcine brain natriuretic peptide was determined by combined use of NMR spectroscopy and distance geometry. A set of 157 inter-proton-distance constraints was derived from the two-dimensional NOE spectra, and further a set of three hydrogen bond constraints was obtained from analysis of the temperature dependence of labile protons. The five structures with minimal violations were selected after performing distance-geometry calculations starting from 40 random initial conformations. The distance-geometry structures were further refined by the use of restrained energy minimization and restrained molecular dynamics. This structure shows a compact conformation with the carboxy-terminal region, Asn21-Tyr26, folded back to the disulfide-linked loop region, Cys4-Cys20. The characteristics of the conformation determined are as follows: conformations of the three segments interposed by glycine residues, which are Arg7-Ile12, Ser14-Leu18 and Cys20-Arg25, were well defined and the segments Arg7-Ile12 and Cys20-Arg25 are rather close to each other and nearly parallel. The biological significance of these local conformations is discussed on the basis of comparisons with those of atrial natriuretic peptide reported by Kobayashi et al.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

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.