Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Evidence in support of a four transmembrane-pore-transmembrane topology model for the Arabidopsis thaliana Na+/K+ translocating AtHKT1 protein, a member of the superfamily of K+ transporters

The Arabidopsis thaliana AtHKT1 protein, a Na+/K+ transporter, is capable of mediating inward Na+ currents in Xenopus laevis oocytes and K+ uptake in Escherichia coli. HKT1 proteins are members of a superfamily of K+ transporters. These proteins have been proposed to contain eight transmembrane segments and four pore-forming regions arranged in a mode similar to that of a K+ channel tetramer. However, computer analysis of the AtHKT1 sequence identified eleven potential transmembrane segments. We have investigated the membrane topology of AtHKT1 with three different techniques. First, a gene fusion alkaline phosphatase study in E. coli clearly defined the topology of the N-terminal and middle region of AtHKT1, but the model for membrane folding of the C-terminal region had to be refined. Second, with a reticulocyte-lysate supplemented with dog-pancreas microsomes, we demonstrated that N-glycosylation occurs at position 429 of AtHKT1. An engineered unglycosylated protein variant, N429Q, mediated Na+ currents in X. laevis oocytes with the same characteristics as the wild-type protein, indicating that N-glycosylation is not essential for the functional expression and membrane targeting of AtHKT1. Five potential glycosylation sites were introduced into the N429Q. Their pattern of glycosylation supported the model based on the E. coli-alkaline phosphatase data. Third, immunocytochemical experiments with FLAG-tagged AtHKT1 in HEK293 cells revealed that the N and C termini of AtHKT1, and the regions containing residues 135–142 and 377–384, face the cytosol, whereas the region of residues 55–62 is exposed to the outside. Taken together, our results show that AtHKT1 contains eight transmembrane-spanning segments.

Novel polymers: Molecular to nanoscale order in three dimensions

The assembly of polymer chains in solution is a powerful method that is leading to the preparation of interesting and unique macromolecular-based synthetic nanostructures. Specific control over the intramolecular and intermolecular physical interactions dictates either the folding of single chains or the aggregation and ordering of multiple chains. This control is provided through the selective placement of functional groups along the polymer backbone and the relative strengths of their attractive and repulsive interactions.

Kinetic stability as a mechanism for protease longevity

The folding of the extracellular serine protease, α-lytic protease (αLP; EC 3.4.21.12) reveals a novel mechanism for stability that appears to lead to a longer functional lifetime for the protease. For αLP, stability is based not on thermodynamics, but on kinetics. Whereas this has required the coevolution of a pro region to facilitate folding, the result has been the optimization of native-state properties independent of their consequences on thermodynamic stability. Structural and mutational data lead to a model for catalysis of folding in which the pro region binds to a conserved β-hairpin in the αLP C-terminal domain, stabilizing the folding transition state and the native state. The pro region is then proteolytically degraded, leaving the active αLP trapped in a metastable conformation. This metastability appears to be a consequence of pressure to evolve properties of the native state, including a large, highly cooperative barrier to unfolding, and extreme rigidity, that reduce susceptibility to proteolytic degradation. In a test of survival under highly proteolytic conditions, homologous mammalian proteases that have not evolved kinetic stability are much more rapidly degraded than αLP. Kinetic stability as a means to longevity is likely to be a mechanism conserved among the majority of extracellular bacterial pro-proteases and may emerge as a general strategy for intracellular eukaryotic proteases subject to harsh conditions as well.

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.

Symmetry and the energy landscapes of biomolecules

The role of symmetry in the folding of proteins is discussed using energy landscape theory. An analytical argument shows it is much easier to find sequences with funneled energy landscape capable of fast folding if the structure is symmetric. The analogy with phase transitions of small clusters with magic numbers is discussed.

A statistical mechanical model for β-hairpin kinetics

Understanding the mechanism of protein secondary structure formation is an essential part of the protein-folding puzzle. Here, we describe a simple statistical mechanical model for the formation of a β-hairpin, the minimal structural element of the antiparallel β-pleated sheet. The model accurately describes the thermodynamic and kinetic behavior of a 16-residue, β-hairpin-forming peptide, successfully explaining its two-state behavior and apparent negative activation energy for folding. The model classifies structures according to their backbone conformation, defined by 15 pairs of dihedral angles, and is further simplified by considering only the 120 structures with contiguous stretches of native pairs of backbone dihedral angles. This single sequence approximation is tested by comparison with a more complete model that includes the 215 possible conformations and 15 × 215 possible kinetic transitions. Finally, we use the model to predict the equilibrium unfolding curves and kinetics for several variants of the β-hairpin peptide.

Optimizing the stability of single-chain proteins by linker length and composition mutagenesis

Linker length and composition were varied in libraries of single-chain Arc repressor, resulting in proteins with effective concentrations ranging over six orders of magnitude (10 μM–10 M). Linkers of 11 residues or more were required for biological activity. Equilibrium stability varied substantially with linker length, reaching a maximum for glycine-rich linkers containing 19 residues. The effects of linker length on equilibrium stability arise from significant and sometimes opposing changes in folding and unfolding kinetics. By fixing the linker length at 19 residues and varying the ratio of Ala/Gly or Ser/Gly in a 16-residue-randomized region, the effects of linker flexibility were examined. In these libraries, composition rather than sequence appears to determine stability. Maximum stability in the Ala/Gly library was observed for a protein containing 11 alanines and five glycines in the randomized region of the linker. In the Ser/Gly library, the most stable protein had seven serines and nine glycines in this region. Analysis of folding and unfolding rates suggests that alanine acts largely by accelerating folding, whereas serine acts predominantly to slow unfolding. These results demonstrate an important role for linker design in determining the stability and folding kinetics of single-chain proteins and suggest strategies for optimizing these parameters.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

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.

Polypeptides of the Maize Amyloplast Stroma1 : Stromal Localization of Starch-Biosynthetic Enzymes and Identification of an 81-Kilodalton Amyloplast Stromal Heat-Shock Cognate

In the developing endosperm of monocotyledonous plants, starch granules are synthesized and deposited within the amyloplast. A soluble stromal fraction was isolated from amyloplasts of immature maize (Zea mays L.) endosperm and analyzed for enzyme activities and polypeptide content. Specific activities of starch synthase and starch-branching enzyme (SBE), but not the cytosolic marker alcohol dehydrogenase, were strongly enhanced in soluble amyloplast stromal fractions relative to soluble extracts obtained from homogenized kernels or endosperms. Immunoblot analysis demonstrated that starch synthase I, SBEIIb, and sugary1, the putative starch-debranching enzyme, were each highly enriched in the amyloplast stroma, providing direct evidence for the localization of starch-biosynthetic enzymes within this compartment. Analysis of maize mutants shows the deficiency of the 85-kD SBEIIb polypeptide in the stroma of amylose extender cultivars and that the dull mutant lacks a >220-kD stromal polypeptide. The stromal fraction is distinguished by differential enrichment of a characteristic group of previously undocumented polypeptides. N-terminal sequence analysis revealed that an abundant 81-kD stromal polypeptide is a member of the Hsp70 family of stress-related proteins. Moreover, the 81-kD stromal polypeptide is strongly recognized by antibodies specific for an Hsp70 of the chloroplast stroma. These findings are discussed in light of implications for the correct folding and assembly of soluble, partially soluble, and granule-bound starch-biosynthetic enzymes during import into the amyloplast.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

Modular construction for function of a ribonucleoprotein enzyme: the catalytic domain of Bacillus subtilis RNase P complexed with B.subtilis RNase P protein

The bacterial RNase P holoenzyme catalyzes the formation of the mature 5′-end of tRNAs and is composed of an RNA and a protein subunit. Among the two folding domains of the RNase P RNA, the catalytic domain (C-domain) contains the active site of this ribozyme. We investigated specific binding of the Bacillus subtilis C-domain with the B.subtilis RNase P protein and examined the catalytic activity of this C-domain–P protein complex. The C-domain forms a specific complex with the P protein with a binding constant of ∼0.1 µM. The C-domain–P protein complex and the holoenzyme are equally efficient in cleaving single-stranded RNA (∼0.9 min–1 at pH 7.8) and substrates with a hairpin–loop 3′ to the cleavage site (∼40 min–1). The holoenzyme reaction is much more efficient with a pre-tRNA substrate, binding at least 100-fold better and cleaving 10–500 times more efficiently. These results demonstrate that the RNase P holoenzyme is functionally constructed in three parts. The catalytic domain alone contains the active site, but has little specificity and affinity for most substrates. The specificity and affinity for the substrate is generated by either the specificity domain of RNase P RNA binding to a T stem–loop-like hairpin or RNase P protein binding to a single-stranded RNA. This modular construction may be exploited to obtain RNase P-based ribonucleoprotein complexes with altered substrate specificity.

Topological challenges to DNA replication: Conformations at the fork

The unwinding of the parental DNA duplex during replication causes a positive linking number difference, or superhelical strain, to build up around the elongating replication fork. The branching at the fork and this strain bring about different conformations from that of (−) supercoiled DNA that is not being replicated. The replicating DNA can form (+) precatenanes, in which the daughter DNAs are intertwined, and (+) supercoils. Topoisomerases have the essential role of relieving the superhelical strain by removing these structures. Stalled replication forks of molecules with a (+) superhelical strain have the additional option of regressing, forming a four-way junction at the replication fork. This four-way junction can be acted on by recombination enzymes to restart replication. Replication and chromosome folding are made easier by topological domain barriers, which sequester the substrates for topoisomerases into defined and concentrated regions. Domain barriers also allow replicated DNA to be (−) supercoiled. We discuss the importance of replicating DNA conformations and the roles of topoisomerases, focusing on recent work from our laboratory.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins.

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.