Molecular computation: RNA solutions to chess problems

We have expanded the field of “DNA computers” to RNA and present a general approach for the solution of satisfiability problems. As an example, we consider a variant of the “Knight problem,” which asks generally what configurations of knights can one place on an n × n chess board such that no knight is attacking any other knight on the board. Using specific ribonuclease digestion to manipulate strands of a 10-bit binary RNA library, we developed a molecular algorithm and applied it to a 3 × 3 chessboard as a 9-bit instance of this problem. Here, the nine spaces on the board correspond to nine “bits” or placeholders in a combinatorial RNA library. We recovered a set of “winning” molecules that describe solutions to this problem.

Folding and activity of circularly permuted forms of a polytopic membrane protein

The transmembrane subunit of the Glc transporter (IICBGlc), which mediates uptake and concomitant phosphorylation of glucose, spans the membrane eight times. Variants of IICBGlc with the native N and C termini joined and new N and C termini in the periplasmic and cytoplasmic surface loops were expressed in Escherichia coli. In vivo transport/in vitro phosphotransferase activities of the circularly permuted variants with the termini in the periplasmic loops 1 to 4 were 35/58, 32/37, 0/3, and 0/0% of wild type, respectively. The activities of the variants with the termini in the cytoplasmic loops 1 to 3 were 0/25, 0/4 and 24/70, respectively. Fusion of alkaline phosphatase to the periplasmic C termini stabilized membrane integration and increased uptake and/or phosphorylation activities. These results suggest that internal signal anchor and stop transfer sequences can function as N-terminal signal sequences in a circularly permuted α-helical bundle protein and that the orientation of transmembrane segments is determined by the amino acid sequence and not by the sequential appearance during translation. Of the four IICBGlc variants with new termini in periplasmic loops, only the one with the discontinuity in loop 4 is inactive. The sequences of loop 4 and of the adjacent TM7 and TM8 are conserved in all phosphoenolpyruvate-dependent carbohydrate:phosphotransferase system transporters of the glucose family.

Methionine-independent initiation of translation in the capsid protein of an insect RNA virus

Protein synthesis is believed to be initiated with the amino acid methionine because the AUG translation initiation codon of mRNAs is recognized by the anticodon of initiator methionine transfer RNA. A group of positive-stranded RNA viruses of insects, however, lacks an AUG translation initiation codon for their capsid protein gene, which is located at the downstream part of the genome. The capsid protein of one of these viruses, Plautia stali intestine virus, is synthesized by internal ribosome entry site-mediated translation. Here we report that methionine is not the initiating amino acid in the translation of the capsid protein in this virus. Its translation is initiated with glutamine encoded by a CAA codon that is the first codon of the capsid-coding region. The nucleotide sequence immediately upstream of the capsid-coding region interacts with a loop segment in the stem–loop structure located 15–43 nt upstream of the 5′ end of the capsid-coding region. The pseudoknot structure formed by this base pair interaction is essential for translation of the capsid protein. This mechanism for translation initiation differs from the conventional one in that the initiation step controlled by the initiator methionine transfer RNA is not necessary.

Coupling between protein folding and allostery in the GroE chaperonin system

GroEL is an allosteric protein that facilitates protein folding in an ATP-dependent manner. Herein, the relationship between cooperative ATP binding by GroEL and the kinetics of GroE-assisted folding of two substrates with different GroES dependence, mouse dihydrofolate reductase (mDHFR) and mitochondrial malate dehydrogenase, is examined by using cooperativity mutants of GroEL. Strong intra-ring positive cooperativity in ATP binding by GroEL decreases the rate of GroEL-assisted mDHFR folding owing to a slow rate of the ATP-induced transition from the protein-acceptor state to the protein-release state. Inter-ring negative cooperativity in ATP binding by GroEL is found to affect the kinetic partitioning of mDHFR, but not of mitochondrial malate dehydrogenase, between folding in solution and folding in the cavity underneath GroES. Our results show that protein folding by this “two-stroke motor” is coupled to cooperative ATP binding.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Denaturant-induced movement of the transition state of protein folding revealed by high-pressure stopped-flow measurements

The small all-β protein tendamistat folds and unfolds with two-state kinetics. We determined the volume changes associated with the folding process by performing kinetic and equilibrium measurements at variable pressure between 0.1 and 100 MPa (1 to 1,000 bar). GdmCl-induced equilibrium unfolding transitions reveal that the volume of the native state is increased by 41.4 ± 2.0 cm3/mol relative to the unfolded state. This value is virtually independent of denaturant concentration. The use of a high-pressure stopped-flow instrument enabled us to measure the activation volumes for the refolding (ΔVf0‡) and unfolding reaction (ΔVu0‡) over a broad range of GdmCl concentrations. The volume of the transition state is 60% native-like (ΔVf0‡ = 25.0 ± 1.2 cm3/mol) in the absence of denaturant, indicating partial solvent accessibility of the core residues. The volume of the transition state increases linearly with denaturant concentration and exceeds the volume of the native state above 6 M GdmCl. This result argues for a largely desolvated transition state with packing deficiencies at high denaturant concentrations and shows that the structure of the transition state depends strongly on the experimental conditions.

Structure and function in rhodopsin: Further elucidation of the role of the intradiscal cysteines, Cys-110, -185, and -187, in rhodopsin folding and function

The disulfide bond between Cys-110 and Cys-187 in the intradiscal domain is required for correct folding in vivo and function of mammalian rhodopsin. Misfolding in rhodopsin, characterized by the loss of ability to bind 11-cis-retinal, has been shown to be caused by an intradiscal disulfide bond different from the above native disulfide bond. Further, naturally occurring single mutations of the intradiscal cysteines (C110F, C110Y, and C187Y) are associated with retinitis pigmentosa (RP). To elucidate further the role of every one of the three intradiscal cysteines, mutants containing single-cysteine replacements by alanine residues and the above three RP mutants have been studied. We find that C110A, C110F, and C110Y all form a disulfide bond between C185 and C187 and cause loss of retinal binding. C185A allows the formation of a C110–C187 disulfide bond, with wild-type-like rhodopsin phenotype. C187A forms a disulfide bond between C110 and C185 and binds retinal, and the pigment formed has markedly altered bleaching behavior. However, the opsin from the RP mutant C187Y forms no rhodopsin chromophore.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Stretching lattice models of protein folding

A new class of experiments that probe folding of individual protein domains uses mechanical stretching to cause the transition. We show how stretching forces can be incorporated in lattice models of folding. For fast folding proteins, the analysis suggests a complex relation between the force dependence and the reaction coordinate for folding.

Observation of strange kinetics in protein folding

Highly nonexponential folding kinetics in aqueous solution have been observed during temperature jump-induced refolding of two proteins, yeast phosphoglycerate kinase and a ubiquitin mutant. The observations are most easily interpreted in terms of downhill folding, which posits a heterogeneous ensemble of structures en route to the folded state. The data are also reconciled with exponential kinetics measured under different experimental conditions and with titration experiments indicating cooperative folding.

Long RNA hairpins that contain inosine are present in Caenorhabditis elegans poly(A)+ RNA

Adenosine deaminases that act on RNA (ADARs) are RNA-editing enzymes that convert adenosine to inosine within double-stranded RNA. In the 12 years since the discovery of ADARs only a few natural substrates have been identified. These substrates were found by chance, when genomically encoded adenosines were identified as guanosines in cDNAs. To advance our understanding of the biological roles of ADARs, we developed a method for systematically identifying ADAR substrates. In our first application of the method, we identified five additional substrates in Caenorhabditis elegans. Four of those substrates are mRNAs edited in untranslated regions, and one is a noncoding RNA edited throughout its length. The edited regions are predicted to form long hairpin structures, and one of the RNAs encodes POP-1, a protein involved in cell fate decisions.

Ubiquitination of RNA polymerase II large subunit signaled by phosphorylation of carboxyl-terminal domain

A sensitive assay using biotinylated ubiquitin revealed extensive ubiquitination of the large subunit of RNA polymerase II during incubations of transcription reactions in vitro. Phosphorylation of the repetitive carboxyl-terminal domain of the large subunit was a signal for ubiquitination. Specific inhibitors of cyclin-dependent kinase (cdk)-type kinases suppress the ubiquitination reaction. These kinases are components of transcription factors and have been shown to phosphorylate the carboxyl-terminal domain. In both regulation of transcription and DNA repair, phosphorylation of the repetitive carboxyl-terminal domain by kinases might signal degradation of the polymerase.

Structural characterization of an engineered tandem repeat contrasts the importance of context and sequence in protein folding

To test a different approach to understanding the relationship between the sequence of part of a protein and its conformation in the overall folded structure, the amino acid sequence corresponding to an α-helix of T4 lysozyme was duplicated in tandem. The presence of such a sequence repeat provides the protein with “choices” during folding. The mutant protein folds with almost wild-type stability, is active, and crystallizes in two different space groups, one isomorphous with wild type and the other with two molecules in the asymmetric unit. The fold of the mutant is essentially the same in all cases, showing that the inserted segment has a well-defined structure. More than half of the inserted residues are themselves helical and extend the helix present in the wild-type protein. Participation of additional duplicated residues in this helix would have required major disruption of the parent structure. The results clearly show that the residues within the duplicated sequence tend to maintain a helical conformation even though the packing interactions with the remainder of the protein are different from those of the original helix. It supports the hypothesis that the structures of individual α-helices are determined predominantly by the nature of the amino acids within the helix, rather than the structural environment provided by the rest of the protein.

Cell cycle-dependent regulation of RNA polymerase I transcription: The nucleolar transcription factor UBF is inactive in mitosis and early G1

Transcription of ribosomal RNA genes by RNA polymerase (pol) I oscillates during the cell cycle, being maximal in S and G2 phase, repressed during mitosis, and gradually recovering during G1 progression. We have shown that transcription initiation factor (TIF)-IB/SL1 is inactivated during mitosis by cdc2/cyclin B-directed phosphorylation of TAFI110. In this study, we have monitored reactivation of transcription after exit from mitosis. We demonstrate that the pol I factor UBF is also inactivated by phosphorylation but recovers with different kinetics than TIF-IB/SL1. Whereas TIF-IB/SL1 activity is rapidly regained on entry into G1, UBF is reactivated later in G1, concomitant with the onset of pol I transcription. Repression of pol I transcription in mitosis and early G1 can be reproduced with either extracts from cells synchronized in M or G1 phase or with purified TIF-IB/SL1 and UBF isolated in the presence of phosphatase inhibitors. The results suggest that two basal transcription factors, e.g., TIF-IB/SL1 and UBF, are inactivated at mitosis and reactivated by dephosphorylation at the exit from mitosis and during G1 progression, respectively.

Laser-mediated, site-specific inactivation of RNA transcripts

The biological function of specific gene products often is determined experimentally by blocking their expression in an organism and observing the resulting phenotype. Chromophore-assisted laser inactivation using malachite green (MG)-tagged antibodies makes it possible to inactivate target proteins in a highly restricted manner, probing their temporally and spatially resolved functions. In this report, we describe the isolation and in vitro characterization of a MG-binding RNA motif that may enable the same high-resolution analysis of gene function specifically at the RNA level (RNA-chromophore-assisted laser inactivation). A well-defined asymmetric internal bulge within an RNA duplex allows high affinity and high specificity binding by MG. Laser irradiation in the presence of low concentrations of MG induces destruction of the MG-binding RNA but not of coincubated control RNA. Laser-induced hydrolysis of the MG-binding RNA is restricted predominantly to a single nucleotide within the bulge. By appropriately incorporating this motif into a target gene, transcripts generated by the gene may be effectively tagged for laser-mediated destruction.