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.

Structure and function in rhodopsin: correct folding and misfolding in point mutants at and in proximity to the site of the retinitis pigmentosa mutation Leu-125-->Arg in the transmembrane helix C.

L125R is a mutation in the transmembrane helix C of rhodopsin that is associated with autosomal dominant retinitis pigmentosa. To probe the orientation of the helix and its packing in the transmembrane domain, we have prepared and studied the mutations E122R, I123R, A124R, S127R, L125F, and L125A at, and in proximity to, the above mutation site. Like L125R, the opsin expressed in COS-1 cells from E122R did not bind 11-cis-retinal, whereas those from I123R and S127R formed the rhodopsin chromophore partially. A124R opsin formed the rhodopsin chromophore (lambda max 495 nm) in the dark, but the metarhodopsin II formed on illumination decayed about 6.5 times faster than that of the wild type and was defective in transducin activation. The mutant opsins from L125F and L125A bound 11-cis-retinal only partially, and in both cases, the mixtures of the proteins produced were separated into retinal-binding and non-retinal-binding (misfolded) fractions. The purified mutant rhodopsin from L125F showed lambda max at 500 nm, whereas that from L125A showed lambda max at 503 nm. The mutant rhodopsin L125F showed abnormal bleaching behavior and both mutants on illumination showed destabilized metarhodopsin II species and reduced transducin activation. Because previous results have indicated that misfolding in rhodopsin is due to the formation of a disulfide bond other than the normal disulfide bond between Cys-110 and Cys-187 in the intradiscal domain, we conclude from the misfolding in mutants L125F and L125A that the folding in vivo in the transmembrane domain is coupled to that in the intradiscal domain.

Photolabeling reveals the proximity of the alpha-neurotoxin binding site to the M2 helix of the ion channel in the nicotinic acetylcholine receptor.

A photoactivatable derivative of neurotoxin II from Naja naja oxiana containing a 125I-labeled p-azidosalicylamidoethyl-1,3'-dithiopropyl label at Lys-25 forms a photo-induced cross-link with the delta subunit of the membrane-bound Torpedo californica nicotinic acetylcholine receptor (AChR). The cross-linked radioactive receptor peptide was isolated by reverse-phase HPLC after tryptic digestion of the labeled delta subunit. The sequence of this peptide, delta-(260-277), and the position of the label at Ala-268 were established by matrix-assisted laser-desorption-ionization mass spectrometry based on the molecular mass and on post-source decay fragment analysis. With the known dimensions of the AChR molecule, of the photolabel, and of alpha-neurotoxin, finding the cross-link at delta Ala-268 (located in the upper part of the channel-forming transmembrane helix M2) means that the center of the alpha-neurotoxin binding site is situated at least approximately 40 A from the extracellular surface of the AChR, proximal to the channel axis.

Two nuclear localization signals present in the basic-helix 1 domains of MyoD promote its active nuclear translocation and can function independently.

MyoD, a member of the family of helix-loop-helix myogenic factors that plays a crucial role in skeletal muscle differentiation, is a nuclear phosphoprotein. Using microinjection of purified MyoD protein into rat fibroblasts, we show that the nuclear import of MyoD is a rapid and active process, being ATP and temperature dependent. Two nuclear localization signals (NLSs), one present in the basic region and the other in the helix 1 domain of MyoD protein, are demonstrated to be functional in promoting the active nuclear transport of MyoD. Synthetic peptides spanning these two NLSs and biochemically coupled to IgGs can promote the nuclear import of microinjected IgG conjugates in muscle and nonmuscle cells. Deletion analysis reveals that each sequence can function independently within the MyoD protein since concomittant deletion of both sequences is required to alter the nuclear import of this myogenic factor. In addition, the complete cytoplasmic retention of a beta-galactosidase-MyoD fusion mutant protein, double deleted at these two NLSs, argues against the existence of another functional NLS motif in MyoD.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

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.

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.

Alanine is helix-stabilizing in both template-nucleated and standard peptide helices

Alanine-based peptides of defined sequence and length show measurable helix contents, allowing them to be used as a model system both for analyzing the mechanism of helix formation and for investigating the contributions of side-chain interactions to protein stability. Extensive characterization of many peptide sequences with varying amino acid contents indicates that the favorable helicity of alanine-based peptides can be attributed to the large helix-stabilizing propensity of alanine. Based on their analysis of alanine-rich sequences N-terminally linked to a synthetic helix-inducing template, Kemp and coworkers [Kemp, D. S., Boyd, J. G. & Muendel, C. C. (1991) Nature (London) 352, 451–454; Kemp, D. S., Oslick, S. L. & Allen, T. J. (1996) J. Am. Chem. Soc. 118, 4249–4255] argue that alanine is helix-indifferent, however, and that the favorable helix contents of alanine-based peptides must have some other explanation. Here, we show that the helix contents of template-nucleated sequences are influenced strongly by properties of the template–helix junction. A model in which the helix propensities of residues at the template–peptide junction are treated separately brings the results from alanine-based peptides and template-nucleated helices into agreement. The resulting model provides a physically plausible resolution of the discrepancies between the two systems and allows the helix contents of both template-nucleated and standard peptide helices to be predicted by using a single set of helix propensities. Helix formation in both standard peptides and template–peptide conjugates can be attributed to the large intrinsic helix-forming tendency of alanine.

Interacting helical faces of subunits a and c in the F1Fo ATP synthase of Escherichia coli defined by disulfide cross-linking

Subunits a and c of Fo are thought to cooperatively catalyze proton translocation during ATP synthesis by the Escherichia coli F1Fo ATP synthase. Optimizing mutations in subunit a at residues A217, I221, and L224 improves the partial function of the cA24D/cD61G double mutant and, on this basis, these three residues were proposed to lie on one face of a transmembrane helix of subunit a, which then interacted with the transmembrane helix of subunit c anchoring the essential aspartyl group. To test this model, in the present work Cys residues were introduced into the second transmembrane helix of subunit c and the predicted fourth transmembrane helix of subunit a. After treating the membrane vesicles of these mutants with Cu(1,10-phenanthroline)2SO4 at 0°, 10°, or 20°C, strong a–c dimer formation was observed at all three temperatures in membranes of 7 of the 65 double mutants constructed, i.e., in the aS207C/cI55C, aN214C/cA62C, aN214C/cM65C, aI221C/cG69C, aI223C/cL72C, aL224C/cY73C, and aI225C/cY73C double mutant proteins. The pattern of cross-linking aligns the helices in a parallel fashion over a span of 19 residues with the aN214C residue lying close to the cA62C and cM65C residues in the middle of the membrane. Lesser a–c dimer formation was observed in nine other double mutants after treatment at 20°C in a pattern generally supporting that indicated by the seven landmark residues cited above. Cross-link formation was not observed between helix-1 of subunit c and helix-4 of subunit a in 19 additional combinations of doubly Cys-substituted proteins. These results provide direct chemical evidence that helix-2 of subunit c and helix-4 of subunit a pack close enough to each other in the membrane to interact during function. The proximity of helices supports the possibility of an interaction between Arg210 in helix-4 of subunit a and Asp61 in helix-2 of subunit c during proton translocation, as has been suggested previously.

The efficiency of propulsion by a rotating flagellum

[At very low Reynolds number, the regime in which fluid dynamics is governed by Stokes equations, a helix that translates along its axis under an external force but without an external torque will necessarily rotate. By the linearity of the Stokes equations, the same helix that is caused to rotate due to an external torque will necessarily translate. This is the physics that underlies the mechanism of flagellar propulsion employed by many microorganisms. Here, I examine the linear relationships between forces and torques and translational and angular velocities of helical objects to understand the nature of flagellar propulsion.]

Oncogenic transformation induced by the Qin protein is correlated with transcriptional repression

The retroviral oncogene qin codes for a protein that belongs to the family of the winged helix transcription factors. The viral Qin protein, v-Qin, differs from its cellular counterpart, c-Qin, by functioning as a stronger transcriptional repressor and a more efficient inducer of tumors. This observation suggests that repression may be important in tumorigenesis. To test this possibility, chimeric proteins were constructed in which the Qin DNA-binding domain was fused to either a strong repressor domain (derived from the Drosophila Engrailed protein) or a strong activator domain (from the herpes simplex virus VP16 protein). The chimeric transcriptional repressor, Qin–Engrailed, transformed chicken embryo fibroblasts in culture and induced sarcomas in young chickens. The chimeric activator, Qin–VP16, failed to transform cells in vitro or in vivo and caused cellular resistance to oncogenic transformation by Qin. These data support the conclusion that the Qin protein induces oncogenic transformation by repressing the transcription of genes which function as negative growth regulators or tumor suppressors.

Mechanical separation of the complementary strands of DNA

We describe the mechanical separation of the two complementary strands of a single molecule of bacteriophage λ DNA. The 3′ and 5′ extremities on one end of the molecule are pulled progressively apart, and this leads to the opening of the double helix. The typical forces along the opening are in the range of 10–15 pN. The separation force signal is shown to be related to the local GC vs. AT content along the molecule. Variations of this content on a typical scale of 100–500 bases are presently detected.

A curved RNA helix incorporating an internal loop with G·A and A·A non-Watson–Crick base pairing

The crystal structure of the RNA dodecamer 5′-GGCC(GAAA)GGCC-3′ has been determined from x-ray diffraction data to 2.3-Å resolution. In the crystal, these oligomers form double helices around twofold symmetry axes. Four consecutive non-Watson–Crick base pairs make up an internal loop in the middle of the duplex, including sheared G·A pairs and novel asymmetric A·A pairs. This internal loop sequence produces a significant curvature and narrowing of the double helix. The helix is curved by 34° from end to end and the diameter is narrowed by 24% in the internal loop. A Mn2+ ion is bound directly to the N7 of the first guanine in the Watson–Crick region following the internal loop and the phosphate of the preceding residue. This Mn2+ location corresponds to a metal binding site observed in the hammerhead catalytic RNA.

Functional analysis of DNA bending and unwinding by the high mobility group domain of LEF-1

LEF-1 (lymphoid enhancer-binding factor 1) is a cell type-specific member of the family of high mobility group (HMG) domain proteins that recognizes a specific nucleotide sequence in the T cell receptor (TCR) α enhancer. In this study, we extend the analysis of the DNA-binding properties of LEF-1 and examine their contributions to the regulation of gene expression. We find that LEF-1, like nonspecific HMG-domain proteins, can interact with irregular DNA structures such as four-way junctions, albeit with lower efficiency than with specific duplex DNA. We also show by a phasing analysis that the LEF-induced DNA bend is directed toward the major groove. In addition, we find that the interaction of LEF-1 with a specific binding site in circular DNA changes the linking number of DNA and unwinds the double helix. Finally, we identified two nucleotides in the LEF-1-binding site that are important for protein-induced DNA bending. Mutations of these nucleotides decrease both the extent of DNA bending and the transactivation of the TCRα enhancer by LEF-1, suggesting a contribution of protein-induced DNA bending to the function of TCRα enhancer.

Ol-Prx 3, a member of an additional class of homeobox genes, is unimodally expressed in several domains of the developing and adult central nervous system of the medaka (Oryzias latipes)

Large-scale genetic screens for mutations affecting early neurogenesis of vertebrates have recently been performed with an aquarium fish, the zebrafish. Later stages of neural morphogenesis have attracted less attention in small fish species, partly because of the lack of molecular markers of developing structures that may facilitate the detection of discrete structural alterations. In this context, we report the characterization of Ol-Prx 3 (Oryzias latipes-Prx 3). This gene was isolated in the course of a large-scale screen for brain cDNAs containing a highly conserved DNA binding region, the homeobox helix-three. Sequence analysis revealed that this gene belongs to another class of homeobox genes, together with a previously isolated mouse ortholog, called OG-12 [Rovescalli, A. C., Asoh, S. & Nirenberg, M. (1996) Proc. Natl. Acad. Sci. USA 93, 10691–10696] and with the human SHOX gene [Rao, E., Weiss, B., Fukami, M., Rump, A., Niesler, B., et al. (1997) Nat. Genet. 16, 54–62], thought to be involved in the short-stature phenotype of Turner syndrome patients. These three genes exhibit a moderate level of identity in the homeobox with the other genes of the paired-related (PRX) gene family. Ol-Prx 3, as well as the PRX genes, are expressed in various cartilaginous structures of head and limbs. These genes might thus be involved in common regulatory pathways during the morphogenesis of these structures. Moreover, this paper reports a complex and monophasic pattern of Ol-Prx 3 expression in the central nervous system, which differs markedly from the patterns reported for the PRX genes, Prx 3 excluded: this gene begins to be expressed in a variety of central nervous system territories at late neurula stage. Strikingly, it remains turned on in some of the derivatives of each territory during the entire life of the fish. We hope this work will thus help identify common features for the PRX 3 family of homeobox genes.