Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

The central structural feature of the membrane fusion protein subunit from the Ebola virus glycoprotein is a long triple-stranded coiled coil

The ectodomain of the Ebola virus Gp2 glycoprotein was solubilized with a trimeric, isoleucine zipper derived from GCN4 (pIIGCN4) in place of the hydrophobic fusion peptide at the N terminus. This chimeric molecule forms a trimeric, highly α-helical, and very thermostable molecule, as determined by chemical crosslinking and circular dichroism. Electron microscopy indicates that Gp2 folds into a rod-like structure like influenza HA2 and HIV-1 gp41, providing further evidence that viral fusion proteins from diverse families such as Orthomyxoviridae (Influenza), Retroviridae (HIV-1), and Filoviridae (Ebola) share common structural features, and suggesting a common membrane fusion mechanism.

HUNT: launch of a full-length cDNA database from the Helix Research Institute

The Helix Research Institute (HRI) in Japan is releasing 4356 HUman Novel Transcripts and related information in the newly established HUNT database. The institute is a joint research project principally funded by the Japanese Ministry of International Trade and Industry, and the clones were sequenced in the governmental New Energy and Industrial Technology Development Organization (NEDO) Human cDNA Sequencing Project. The HUNT database contains an extensive amount of annotation from advanced analysis and represents an essential bioinformatics contribution towards understanding of the gene function. The HRI human cDNA clones were obtained from full-length enriched cDNA libraries constructed with the oligo-capping method and have resulted in novel full-length cDNA sequences. A large fraction has little similarity to any proteins of known function and to obtain clues about possible function we have developed original analysis procedures. Any putative function deduced here can be validated or refuted by complementary analysis results. The user can also extract information from specific categories like PROSITE patterns, PFAM domains, PSORT localization, transmembrane helices and clones with GENIUS structure assignments. The HUNT database can be accessed at http://www.hri.co.jp/HUNT.

Sequences upstream of the branch site are required to form helix II between U2 and U6 snRNA in a trans-splicing reaction

Three different base paired stems form between U2 and U6 snRNA over the course of the mRNA splicing reaction (helices I, II and III). One possible function of U2/U6 helix II is to facilitate subsequent U2/U6 helix I and III interactions, which participate directly in catalysis. Using an in vitro trans-splicing assay, we investigated the function of sequences located just upstream from the branch site (BS). We find that these upstream sequences are essential for stable binding of U2 to the branch region, and for U2/U6 helix II formation, but not for initial U2/BS pairing. We also show that non-functional upstream sequences cause U2 snRNA stem–loop IIa to be exposed to dimethylsulfate modification, perhaps reflecting a U2 snRNA conformational change and/or loss of SF3b proteins. Our data suggest that initial binding of U2 snRNP to the BS region must be stabilized by an interaction with upstream sequences before U2/U6 helix II can form or U2 stem–loop IIa can participate in spliceosome assembly.

Lifetimes of intermediates in the β-sheet to α-helix transition of β-lactoglobulin by using a diffusional IR mixer

The extremely slow α-helix/β-sheet transition of proteins is a crucial step in amylogenic diseases and represents an internal rearrangement of local contacts in an already folded protein. These internal structural rearrangements within an already folded protein are a critical aspect of biological action and are a product of conformational flow along unknown metastable local minima of the energy landscape of the compact protein. We use a diffusional IR mixer with time-resolved Fourier transform IR spectroscopy capable of 400-μs time resolution to show that the trifluoroethanol driven β-sheet to α-helix transition of β-lactoglobulin proceeds via a compact β-sheet intermediate with a lifetime of 7 ms, small compared with the overall folding time of β-lactoglobulin.

Direct conversion of an oligopeptide from a β-sheet to an α-helix: A model for amyloid formation

A 16-amino acid oligopeptide forms a stable β-sheet structure in water. In physiological solutions it is able to self-assemble to form a macroscopic matrix that stains with Congo red. On raising the temperature of the aqueous solution above 70°C, an abrupt structural transition occurs in the CD spectra from a β-sheet to a stable α-helix without a detectable random-coil intermediate. With cooling, it retained the α-helical form and took several weeks at room temperature to partially return to the β-sheet form. Slow formation of the stable β-sheet structure thus shows kinetic irreversibility. Such a formation of very stable β-sheet structures is found in the amyloid of a number of neurological diseases. This oligopeptide could be a model system for studying the protein conformational changes that occurs in scrapie or Alzheimer disease. The abrupt and direct conversion from a β-sheet to an α-helix may also be found in other processes, such as protein folding and protein–protein interaction. Furthermore, such drastic structure changes may also be exploited in biomaterials designed as sensors to detect environmental changes.

Conditional Synergism between Cryptochrome 1 and Phytochrome B Is Shown by the Analysis of phyA, phyB, and hy4 Simple, Double, and Triple Mutants in Arabidopsis

Wild-type or phyA, phyB, or hy4 mutant Arabidopsis seedlings lacking phytochrome A (phyA), phytochrome B (phyB), or cryptochrome 1 (cry1), respectively, and the double and triple mutants were used in combination with blue-light treatments given simultaneously with red or far-red light. We investigated the interaction between phytochromes and cry1 in the control of hypocotyl growth and cotyledon unfolding. Under conditions deficient for cry1 (short exposures to blue light) or phyB (far-red background), these photoreceptors acted synergistically: Under short exposures to blue light (3 h/d) added to a red-light background, cry1 activity required phyB (e.g. the hy4 mutant was taller than the wild type but the phyBhy4 mutant was not taller than the phyB mutant). Under prolonged exposures to blue light (24 h/d) added to a far-red light background, phyB activity required cry1 (e.g. the phyAphyB mutant was taller than the phyA mutant but the phyAphyBhy4 mutant was not taller than the phyAhy4 mutant). Under more favorable light inputs, i.e. prolonged exposures to blue light added to a red-light background, the effects of cry1 and phyB were independent. Thus, the synergism between phyB and cry1 is conditional. The effect of cry1 was not reduced by the phyA mutation under any tested light condition. Under continuous blue light the triple mutant phyAphyBhy4 showed reduced hypocotyl growth inhibition and cotyledon unfolding compared with the phyAphyB mutant. The action of cry1 in the phyAphyB double mutant was higher under the red-light than the far-red-light background, indicating a synergistic interaction between cry1 and phytochromes C, D, or E; however, a residual action of cry1 independent of any phytochrome is likely to occur.

Cytochrome c′ folding triggered by electron transfer: Fast and slow formation of four-helix bundles

Reduced (FeII) Rhodopseudomonas palustris cytochrome c′ (Cyt c′) is more stable toward unfolding ([GuHCl]1/2 = 2.9(1) M) than the oxidized (FeIII) protein ([GuHCl]1/2 = 1.9(1) M). The difference in folding free energies (ΔΔGf° = 70 meV) is less than half of the difference in reduction potentials of the folded protein (100 mV vs. NHE) and a free heme in aqueous solution (≈−150 mV). The spectroscopic features of unfolded FeII–Cyt c′ indicate a low-spin heme that is axially coordinated to methionine sulfur (Met-15 or Met-25). Time-resolved absorption measurements after CO photodissociation from unfolded FeII(CO)–Cyt c′ confirm that methionine can bind to the ferroheme on the microsecond time scale [kobs = 5(2) × 104 s−1]. Protein folding was initiated by photoreduction (two-photon laser excitation of NADH) of unfolded FeIII–Cyt c′ ([GuHCl] = 2.02–2.54 M). Folding kinetics monitored by heme absorption span a wide time range and are highly heterogeneous; there are fast-folding (≈103 s−1), intermediate-folding (102–101 s−1), and slow-folding (10−1 s−1) populations, with the last two likely containing methionine-ligated (Met-15 or Met-25) ferrohemes. Kinetics after photoreduction of unfolded FeIII–Cyt c′ in the presence of CO are attributable to CO binding [1.4(6) × 103 s−1] and FeII(CO)–Cyt c′ folding [2.8(9) s−1] processes; stopped-flow triggered folding of FeIII–Cyt c′ (which does not contain a protein-derived sixth ligand) is adequately described by a single kinetics phase with an estimated folding time constant of ≈4 ms [ΔGf° = −33(3) kJ mol−1] at zero denaturant.

Force and kinetic barriers to unzipping of the DNA double helix

A theory of the unzipping of double-stranded DNA is presented and is compared to recent micromanipulation experiments. It is shown that the interactions that stabilize the double helix and the elastic rigidity of single strands simply determine the sequence-dependent ≈12-pN force threshold for DNA strand separation. Using a semimicroscopic model of the binding between nucleotide strands, we show that the greater rigidity of the strands when formed into double-stranded DNA, relative to that of isolated strands, gives rise to a potential barrier to unzipping. The effects of this barrier are derived analytically. The force to keep the extremities of the molecule at a fixed distance, the kinetic rates for strand unpairing at fixed applied force, and the rupture force as a function of loading rate are calculated. The dependence of the kinetics and of the rupture force on molecule length is also analyzed.

Site-directed spin labeling and chemical crosslinking demonstrate that helix V is close to helices VII and VIII in the lactose permease of Escherichia coli.

Site-directed chemical cleavage of lactose permease indicates that helix V is in close proximity to helices VII and VIII. To test this conclusion further, permease containing a biotin-acceptor domain and paired Cys residues at positions 148 (helix V) and 228 (helix VII), 148 and 226 (helix VII), or 148 and 275 (helix VIII) was affinity purified and labeled with a sulfhydryl-specific nitroxide spin label. Spin-spin interactions are observed with the 148/228 and 148/275 pairs, indicating close proximity between appropriate faces of helix V and helices VII and VIII. Little or no interaction is evident with the 148/226 pair, in all likelihood because position 226 is on the opposite face of helix VII from position 228. Broadening of the electron paramagnetic resonance spectra in the frozen state was used to estimate distance between the 148/228 and the 148/275 pairs. The nitroxides at positions 148 and 228 or 148 and 275 are within approximately 13-15 A. Finally, Cys residues at positions 148 and 228 are crosslinked by dibromobimane, a bifunctional crosslinker that is approximately 5 A. long, while no crosslinking is detected between Cys residues at positions 148 and 275 or 148 and 226. The results provide strong support for a structure in which helix V is in close proximity to both helices VII and VIII and is oriented in such a fashion that Cys-148 is closer to helix VII.

Combinatorial control of muscle development by basic helix-loop-helix and MADS-box transcription factors.

Members of the MyoD family of muscle-specific basic helix-loop-helix (bHLH) proteins function within a genetic pathway to control skeletal muscle development. Mutational analyses of these factors suggested that their DNA binding domains mediated interaction with a coregulator required for activation of muscle-specific transcription. Members of the myocyte enhancer binding factor 2 (MEF2) family of MADS-box proteins are expressed at high levels in muscle and neural cells and at lower levels in several other cell types. MEF2 factors are unable to activate muscle gene expression alone, but they potentiate the transcriptional activity of myogenic bHLH proteins. This potentiation appears to be mediated by direct interactions between the DNA binding domains of these different types of transcription factors. Biochemical and genetic evidence suggests that MEF2 factors are the coregulators for myogenic bHLH proteins. The presence of MEF2 and cell-specific bHLH proteins in other cell types raises the possibility that these proteins may also cooperate to regulate other programs of cell-specific gene expression. We present a model to account for such cooperative interactions.

Molecular mechanics calculations of the riboacetal internucleotide linkage in double and triple helices.

Structures of Watson-Crick base paired 15-nucleobase oligomer strands in A-type or B-type conformation in which one strand [a strand of alternating nucleotide and riboacetal thymidine nucleoside (RT) units, RP] is DNA and the other is composed of alternating nucleotides and riboacetal nucleosides have been studied by molecular mechanics. Analogously, oligomer strands of RNA in place of DNA have been modeled. The calculations indicate that the RP strand is more stable when complexed in an A-type duplex relative to a B-type form and that this conformational preference is presumably due to the more uniform nature of the former. Nearly planar ribose rings were more commonly observed in the minimized structures of the B-type DNA.RP duplexes as compared with A-type duplexes, despite the fact that planar ribofuranose rings are known to be energetically unfavorable in oligonucleotides. Computed relative stabilities of all duplexes containing the RP strand suggest that such heteroduplexes are less stable than the corresponding double-stranded DNA and double-stranded RNA species. These findings are in agreement with experimental results which show, when equivalent sequences were compared, that a DNA.RNA control forms a more stable duplex than RP hound to a complementary single-stranded RNA strand. In contrast, molecular mechanics studies of complementary triple-helical (DNA)2.RP, (DNA)2.DNA, and (DNA)2.RNA structures indicate that the binding of RP as a Hoogsteen strand stabilizes the underlying duplex to a greater extent compared with native oligonucleotides. These calculations suggest that puckering of the ribose ring in the riboacetal linkage leads to a more favorable interaction with a complementary nucleic acid target than the proposed planar geometry and that this puckering may account for the enhanced binding of RP to a double-stranded target.

Protonatable residues at the cytoplasmic end of transmembrane helix-2 in the signal transducer HtrI control photochemistry and function of sensory rhodopsin I.

Neutral residue replacements were made of 21 acidic and basic residues within the N-terminal half of the Halobacterium salinarium signal transducer HtrI [the halobacterial transducer for sensory rhodopsin I (SRI)] by site-specific mutagenesis. The replacements are all within the region of HtrI that we previously concluded from deletion analysis to contain sites of interaction with the phototaxis receptor SRI. Immunoblotting shows plasmid expression of the htrI-sopI operon containing the mutations produces SRI and mutant HtrI in cells at near wild-type levels. Six of the HtrI mutations perturb photochemical kinetics of SRI and one reverses the phototaxis response. Substitution with neutral amino acids of Asp-86, Glu-87, and Glu-108 accelerate, and of Arg-70, Arg-84, and Arg-99 retard, the SRI photocycle. Opposite effects on photocycle rate cancel in double mutants containing one replaced acidic and one replaced basic residue. Laser flash spectroscopy shows the kinetic perturbations are due to alteration of the rate of reprotonation of the retinylidene Schiff base. All of these mutations permit normal attractant and repellent signaling. On the other hand, the substitution of Glu-56 with the isosteric glutamine converts the normally attractant effect of orange light to a repellent signal in vivo at neutral pH (inverted signaling). Low pH corrects the inversion due to Glu-56 -> Gln and the apparent pK of the inversion is increased when arginine is substituted at position 56. The results indicate that the cytoplasmic end of transmembrane helix-2 and the initial part of the cytoplasmic domain contain interaction sites with SRI. To explain these and previous results, we propose a model in which (i) the HtrI region identified here forms part of an electrostatic bonding network that extends through the SRI protein and includes its photoactive site; (ii) alteration of this network by photoisomerization-induced Schiff base deprotonation and reprotonation shifts HtrI between attractant and repellent conformations; and (iii) HtrI mutations and extracellular pH alter the equilibrium ratios of these conformations.