Monitoring the structure of Escherichia coli RNase P RNA in the presence of various divalent metal ions

Lead(II)-induced cleavage can be used as a tool to probe conformational changes in RNA. In this report, we have investigated the conformation of M1 RNA, the catalytic subunit of Escherichia coli RNase P, by studying the lead(II)-induced cleavage pattern in the presence of various divalent metal ions. Our data suggest that the overall conformation of M1 RNA is very similar in the presence of Mg2+, Mn2+, Ca2+, Sr2+ and Ba2+, while it is changed compared to the Mg2+-induced conformation in the presence of other divalent metal ions, Cd2+ for example. We also observed that correct folding of some M1 RNA domains is promoted by Pb2+, while folding of other domain(s) requires the additional presence of other divalent metal ions, cobalt(III) hexamine or spermidine. Based on the suppression of Pb2+ cleavage at increasing concentrations of various divalent metal ions, our findings suggest that different divalent metal ions bind with different affinities to M1 RNA as well as to an RNase P hairpin–loop substrate and yeast tRNAPhe. We suggest that this approach can be used to obtain information about the relative binding strength for different divalent metal ions to RNA in general, as well as to specific RNA divalent metal ion binding sites. Of those studied in this report, Mn2+ is generally among the strongest RNA binders.

The thermodynamic origin of the stability of a thermophilic ribozyme

Understanding the mechanism of thermodynamic stability of an RNA structure has significant implications for the function and design of RNA. We investigated the equilibrium folding of a thermophilic ribozyme and its mesophilic homologue by using hydroxyl radical protection, small-angle x-ray scattering, and circular dichroism. Both RNAs require Mg2+ to fold to their native structures that are very similar. The stability is measured as a function of Mg2+ and urea concentrations at different temperatures. The enhanced stability of the thermophilic ribozyme primarily is derived from a tremendous increase in the amount of structure formed in the ultimate folding transition. This increase in structure formation and cooperativity arises because the penultimate and the ultimate folding transitions in the mesophilic ribozyme become linked into a single transition in the folding of the thermophilic ribozyme. Therefore, the starting point, or reference state, for the transition to the native, functional thermophilic ribozyme is significantly less structured. The shift in the reference state, and the resulting increase in folding cooperativity, is likely due to the stabilization of selected native interactions that only form in the ultimate transition. This mechanism of using a less structured intermediate and increased cooperativity to achieve higher functional stability for tertiary RNAs is fundamentally different from that commonly proposed to explain the increased stability of thermophilic proteins.

O-Mannosylation Protects Mutant Alpha-Factor Precursor from Endoplasmic Reticulum-associated Degradation

Secretory proteins that fail to fold in the endoplasmic reticulum (ER) are transported back to the cytosol and degraded by proteasomes. It remains unclear how the cell distinguishes between folding intermediates and misfolded proteins. We asked whether misfolded secretory proteins are covalently modified in the ER before export. We found that a fraction of mutant alpha-factor precursor, but not the wild type, was progressively O-mannosylated in microsomes and in intact yeast cells by protein O-mannosyl transferase 2 (Pmt2p). O-Mannosylation increased significantly in vitro under ER export conditions, i.e., in the presence of ATP and cytosol, and this required export-proficient Sec61p in the ER membrane. Deletion of PMT2, however, did not abrogate mutant alpha-factor precursor degradation but, rather, enhanced its turnover in intact yeast cells. In vitro, O-mannosylated mutant alpha-factor precursor was stable and protease protected, and a fraction was associated with Sec61p in the ER lumen. Thus, prolonged ER residence allows modification of exposed O-mannosyl acceptor sites in misfolded proteins, which abrogates misfolded protein export from the ER at a posttargeting stage. We conclude that there is a limited window of time during which misfolded proteins can be removed from the ER before they acquire inappropriate modifications that can interfere with disposal through the Sec61 channel.

Hydrophobic moments of protein structures: Spatially profiling the distribution

It is generally accepted that globular proteins fold with a hydrophobic core and a hydrophilic exterior. Might the spatial distribution of amino acid hydrophobicity exhibit common features? The hydrophobic profile detailing this distribution from the protein interior to exterior has been examined for 30 relatively diverse structures obtained from the Protein Data Bank, for 3 proteins of the 30S ribosomal subunit, and for a simple set of 14 decoys. A second-order hydrophobic moment has provided a simple measure of the spatial variation. Shapes of the calculated spatial profiles of all native structures have been found to be comparable. Consequently, profile shapes as well as particular profile features should assist in validating predicted protein structures and in discriminating between different protein-folding pathways. The spatial profiles of the 14 decoys are clearly distinguished from the profiles of their native structures.

Tryptophan zippers: Stable, monomeric β-hairpins

A structural motif, the tryptophan zipper (trpzip), greatly stabilizes the β-hairpin conformation in short peptides. Peptides (12 or 16 aa in length) with four different turn sequences are monomeric and fold cooperatively in water, as has been observed previously for some hairpin peptides. However, the folding free energies of the trpzips exceed substantially those of all previously reported β-hairpins and even those of some larger designed proteins. NMR structures of three of the trpzip peptides reveal exceptionally well-defined β-hairpin conformations stabilized by cross-strand pairs of indole rings. The trpzips are the smallest peptides to adopt an unique tertiary fold without requiring metal binding, unusual amino acids, or disulfide crosslinks.

Non-Arrhenius kinetics for the loop closure of a DNA hairpin

Intramolecular chain diffusion is an elementary process in the conformational fluctuations of the DNA hairpin-loop. We have studied the temperature and viscosity dependence of a model DNA hairpin-loop by FRET (fluorescence resonance energy transfer) fluctuation spectroscopy (FRETfs). Apparent thermodynamic parameters were obtained by analyzing the correlation amplitude through a two-state model and are consistent with steady-state fluorescence measurements. The kinetics of closing the loop show non-Arrhenius behavior, in agreement with theoretical prediction and other experimental measurements on peptide folding. The fluctuation rates show a fractional power dependence (β = 0.83) on the solution viscosity. A much slower intrachain diffusion coefficient in comparison to that of polypeptides was derived based on the first passage time theory of SSS [Szabo, A., Schulten, K. & Schulten, Z. (1980) J. Chem. Phys. 72, 4350–4357], suggesting that intrachain interactions, especially stacking interaction in the loop, might increase the roughness of the free energy surface of the DNA hairpin-loop.

Three-dimensional domain swapping in p13suc1 occurs in the unfolded state and is controlled by conserved proline residues

p13suc1 has two native states, a monomer and a domain-swapped dimer. We show that their folding pathways are connected by the denatured state, which introduces a kinetic barrier between monomer and dimer under native conditions. The barrier is lowered under conditions that speed up unfolding, thereby allowing, to our knowledge for the first time, a quantitative dissection of the energetics of domain swapping. The monomer–dimer equilibrium is controlled by two conserved prolines in the hinge loop that connects the exchanging domains. These two residues exploit backbone strain to specifically direct dimer formation while preventing higher-order oligomerization. Thus, the loop acts as a loaded molecular spring that releases tension in the monomer by adopting its alternative conformation in the dimer. There is an excellent correlation between domain swapping and aggregation, suggesting they share a common mechanism. These insights have allowed us to redesign the domain-swapping propensity of suc1 from a fully monomeric to a fully dimeric protein.

Solution structure of DFF40 and DFF45 N-terminal domain complex and mutual chaperone activity of DFF40 and DFF45

Apoptotic DNA fragmentation is mediated by a caspase-activated DNA fragmentation factor (DFF)40. Expression and folding of DFF40 require the presence of DFF45, which also acts as a nuclease inhibitor before DFF40 activation by execution caspases. The N-terminal domains (NTDs) of both proteins are homologous, and their interaction plays a key role in the proper functioning of this two-component system. Here we report that the NTD of DFF45 alone is unstructured in solution, and its folding is induced upon binding to DFF40 NTD. Therefore, folding of both proteins regulates the formation of the DFF40/DFF45 complex. The solution structure of the heterodimeric complex between NTDs of DFF40 and DFF45 reported here shows that the mutual chaperoning includes the formation of an extensive network of intermolecular interactions that bury a hydrophobic cluster inside the interface, surrounded by intermolecular salt bridges.

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.