Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

Electron Tomography of Metaphase Nucleolar Organizer Regions: Evidence for a Twisted-Loop Organization

Metaphase nucleolar organizer regions (NORs), one of four types of chromosome bands, are located on human acrocentric chromosomes. They contain r-chromatin, i.e., ribosomal genes complexed with proteins such as upstream binding factor and RNA polymerase I, which are argyrophilic NOR proteins. Immunocytochemical and cytochemical labelings of these proteins were used to reveal r-chromatin in situ and to investigate its spatial organization within NORs by confocal microscopy and by electron tomography. For each labeling, confocal microscopy revealed small and large double-spotted NORs and crescent-shaped NORs. Their internal three-dimensional (3D) organization was studied by using electron tomography on specifically silver-stained NORs. The 3D reconstructions allow us to conclude that the argyrophilic NOR proteins are grouped as a fiber of 60–80 nm in diameter that constitutes either one part of a turn or two or three turns of a helix within small and large double-spotted NORs, respectively. Within crescent-shaped NORs, virtual slices reveal that the fiber constitutes several longitudinally twisted loops, grouped as two helical 250- to 300-nm coils, each centered on a nonargyrophilic axis of condensed chromatin. We propose a model of the 3D organization of r-chromatin within elongated NORs, in which loops are twisted and bent to constitute one basic chromatid coil.

Myogenic stem cell function is impaired in mice lacking the forkhead/winged helix protein MNF

Myocyte nuclear factor (MNF) is a winged helix transcription factor that is expressed selectively in myogenic stem cells (satellite cells) of adult animals. Using a gene knockout strategy to generate a functional null allele at the Mnf locus, we observed that mice lacking MNF are viable, but severely runted. Skeletal muscles of Mnf−/− animals are atrophic, and satellite cell function is impaired. Muscle regeneration after injury is delayed and incomplete, and the normal timing of expression of cell cycle regulators and myogenic determination genes is dysregulated. Mnf mutant mice were intercrossed with mdx mice that lack dystrophin and exhibit only a subtle myopathic phenotype. In contrast, mdx mice that also lack MNF die in the first few weeks of life with a severe myopathy. Haploinsufficiency at the Mnf locus (Mnf+/−) also exacerbates the mdx phenotype to more closely resemble Duchenne's muscular dystrophy in humans. We conclude that MNF acts to regulate genes that coordinate the proliferation and differentiation of myogenic stem cells after muscle injury. Animals deficient in MNF may prove useful for evaluation of potential therapeutic interventions to promote muscle regeneration for patients having Duchenne's muscular dystrophy.

A general method for determining helix packing in membrane proteins in situ: Helices I and II are close to helix VII in the lactose permease of Escherichia coli

It was previously shown that coexpression of the lactose permease of Escherichia coli in two contiguous fragments leads to functional complementation. We demonstrate here that site-directed thiol crosslinking of coexpressed permease fragments can be used to determine helix proximity in situ without the necessity of purifying the permease. After coexpression of the six N-terminal (N6) and six C-terminal (C6) transmembrane helices, each with a single Cys residue, crosslinking was carried out in native membranes and assessed by the mobility of anti-C-terminal-reactive polypeptides on immunoblots. A Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 28 or 29 (helix I), but not with a Cys residue at position 27, which is on the opposite face of helix I, thereby indicating that the face of helix I containing Pro-28 and Phe-29 is close to helix VII. Similarly, a Cys residue at position 242 or 245 (helix VII) forms a disulfide with a Cys residue at either position 52 or 53 (helix II), but not with a Cys residue at position 54. Furthermore, low-efficiency crosslinking is observed between a Cys residue at position 52 or 53 and a Cys residue at position 361 (helix XI). The results indicate that helix VII lies in close proximity to both helices I and II and that helix II is also close to helix XI. The method should be applicable to a number of different polytopic membrane proteins.

Inactive conformation of the serpin α1-antichymotrypsin indicates two-stage insertion of the reactive loop: Implications for inhibitory function and conformational disease

The serpins are a family of proteinase inhibitors that play a central role in the control of proteolytic cascades. Their inhibitory mechanism depends on the intramolecular insertion of the reactive loop into β-sheet A after cleavage by the target proteinase. Point mutations within the protein can allow aberrant conformational transitions characterized by β-strand exchange between the reactive loop of one molecule and β-sheet A of another. These loop-sheet polymers result in diseases as varied as cirrhosis, emphysema, angio-oedema, and thrombosis, and we recently have shown that they underlie an early-onset dementia. We report here the biochemical characteristics and crystal structure of a naturally occurring variant (Leu-55–Pro) of the plasma serpin α1-antichymotrypsin trapped as an inactive intermediate. The structure demonstrates a serpin configuration with partial insertion of the reactive loop into β-sheet A. The lower part of the sheet is filled by the last turn of F-helix and the loop that links it to s3A. This conformation matches that of proposed intermediates on the pathway to complex and polymer formation in the serpins. In particular, this intermediate, along with the latent and polymerized conformations, explains the loss of activity of plasma α1-antichymotrypsin associated with chronic obstructive pulmonary disease in patients with the Leu-55–Pro mutation.

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.

Detection of competing DNA structures by thermal gradient gel electrophoresis: from self-association to triple helix formation by (G,A)-containing oligonucleotides

Sequence-specific recognition of DNA can be achieved by triple helix-forming oligonucleotides that bind to the major groove of double-helical DNA. These oligonucleotides have been used as sequence-specific DNA ligands for various purposes, including sequence-specific gene regulation in the so-called ‘antigene strategy’. In particular, (G,A)-containing oligonucleotides can form stable triple helices under physiological conditions. However, triplex formation may be in competition with self-association of these oligonucleotides. For biological applications it would be interesting to identify the conditions under which one structure is favoured as compared to the other(s). Here we have directly studied competition between formation of a parallel (G,A) homoduplex and that of a triple helix by a 13 nt (G,A)-containing oligonucleotide. Temperature gradient gel electrophoresis allows simultaneous detection of competition between the two structures, because of their different temperature dependencies and gel electrophoretic mobilities, and characterisation of this competition.

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.

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.