Cotranslational protein folding

Engineering of CYP82Y1, a cytochrome P450 monooxygenase: a key enzyme in noscapine biosynthesis in opium poppy

Protein engineering provides a powerful base for the circumvention of challenges tied with characteristics accountable for enzyme functions. CYP82Y1 introduces a hydroxyl group (-OH) into C1 of N-methylcanadine as the substrate to yield 1-hydroxy-N-methylcanadine. This chemical process has been found to be the gateway to noscapine biosynthesis. Owning to the importance of CYP82Y1 in this biosynthetic pathway, it has been selected as a target for enzyme engineering. The insertion of tags to the N- and C-terminal of CYP82Y1 was assessed for their efficiencies for improvement of the physiological performances of CYP82Y1. Although these attempts achieved some positive results, further strategies are required to dramatically enhance the CYP82Y1 activity. Here methods that have been adopted to achieve a functionally improved CYP82Y1 will be reviewed. In addition, the possibility of recruitment of other techniques having not yet been implemented in CYP82Y1 engineering, including the substitution of the residues located in the substrate recognition site, formation of the synthetic fusion proteins, and construction of the artificial lipid-based scaffold will be discussed. Given the fact that the pace of noscapine synthesis is constrained by the CYP82Y1-catalyzing step, the methods proposed here are capable of accelerating the rate of reaction performed by CYP82Y1 through improving its properties, resulting in the enhancement of noscapine accumulation.

CHIP: A Co-chaperone for Degradation by the Proteasome and Lysosome

Protein homeostasis relies on a balance between protein folding and protein degradation. Molecular chaperones like Hsp70 and Hsp90 fulfill well-defined roles in protein folding and conformational stability via ATP-dependent reaction cycles. These folding cycles are controlled by associations with a cohort of non-client protein co-chaperones, such as Hop, p23, and Aha1. Pro-folding co-chaperones facilitate the transit of the client protein through the chaperone-mediated folding process. However, chaperones are also involved in proteasomal and lysosomal degradation of client proteins. Like folding complexes, the ability of chaperones to mediate protein degradation is regulated by co-chaperones, such as the C-terminal Hsp70-binding protein (CHIP/STUB1). CHIP binds to Hsp70 and Hsp90 chaperones through its tetratricopeptide repeat (TPR) domain and functions as an E3 ubiquitin ligase using a modified RING finger domain (U-box). This unique combination of domains effectively allows CHIP to network chaperone complexes to the ubiquitin-proteasome and autophagosome-lysosome systems. This chapter reviews the current understanding of CHIP as a co-chaperone that switches Hsp70/Hsp90 chaperone complexes from protein folding to protein degradation.

Mechanistic insights into global suppressors of protein folding defects

Most amino acid substitutions in a protein either lead to partial loss-of-function or are near neutral. Several studies have shown the existence of second-site mutations that can rescue defects caused by diverse loss-of-function mutations. Such global suppressor mutations are key drivers of protein evolution. However, the mechanisms responsible for such suppression remain poorly understood. To address this, we characterized multiple suppressor mutations both in isolation and in combination with inactive mutants. We examined six global suppressors of the bacterial toxin CcdB, the known M182T global suppressor of TEM-1 β-lactamase, the N239Y global suppressor of p53-DBD and three suppressors of the SARS-CoV-2 spike Receptor Binding Domain. When coupled to inactive mutants, they promote increased in-vivo solubilities as well as regain-of-function phenotypes. In the case of CcdB, where novel suppressors were isolated, we determined the crystal structures of three such suppressors to obtain insight into the specific molecular interactions responsible for the observed effects. While most individual suppressors result in small stability enhancements relative to wildtype, which can be combined to yield significant stability increments, thermodynamic stabilisation is neither necessary nor sufficient for suppressor action. Instead, in diverse systems, we observe that individual global suppressors greatly enhance the foldability of buried site mutants, primarily through increase in refolding rate parameters measured in vitro. In the crowded intracellular environment, mutations that slow down folding likely facilitate off-pathway aggregation. We suggest that suppressor mutations that accelerate refolding can counteract this, enhancing the yield of properly folded, functional protein in vivo.

Biosynthetic Protein Folding and Molecular Chaperons

The problem of linear polypeptide chain folding into a unique tertiary structure is one of the fundamental scientific challenges. The process of folding cannot be fully understood without its biological context, especially for big multidomain and multisubunit proteins. The principal features of biosynthetic folding are co-translational folding of growing nascent polypeptide chains and involvement of molecular chaperones in the process. The review summarizes available data on the early events of nascent chain folding, as well as on later advanced steps, including formation of elements of native structure. The relationship between the non-uniformity of translation rate and folding of the growing polypeptide is discussed. The results of studies on the effect of biosynthetic folding features on the parameters of folding as a physical process, its kinetics and mechanisms, are presented. Current understanding and hypotheses on the relationship of biosynthetic folding with the fundamental physical parameters and current views on polypeptide folding in the context of energy landscapes are discussed.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Identification and functional characterization of Cofilin-1 as a new member of antimicrobial protein

Cofilin-1 (Cfl1), a member of the ADF/cofilin family, has been identified as one of differentially expressed proteins in human dendritic cells challenged with lipopolysaccharide (LPS), suggesting that it may be involved in immune response. Here we showed that zebrafish cfl1 was markedly up-regulated by LPS and LTA treatment. We also showed that zebrafish recombinant Cfl1 (rCfl1) not only bound to the Gram-negative and positive bacteria A. hydrophila and S. aureus as well as their signature molecules LPS and LTA but also inhibited the growth of the bacteria. Moreover, we found that the heparin-binding motif-containing regions of Cfl1, i.e., Cfl1_9-25, Cfl1_34-51 and Cfl1_108-125, like rCfl1, were also able to bind to LPS and LTA and to inhibit the bacterial growth. rCfl1, Cfl1_9-25, Cfl1_34-51, and Cfl1_108-125 were all able to cause bacterial cell destruction, to induce membrane depolarization, and to stimulate intracellular ROS production. Finally, we showed that zebrafish Cfl1 could protect developing embryos/larvae against attack by the potential pathogen A. hydrophila. These data together indicate that zebrafish Cfl1 plays an immune-relevant role as a newly-characterized antimicrobial protein.

Hydrophobic and electrostatic interactions modulate protein escape at the ribosomal exit tunnel

After translation, nascent proteins must escape the ribosomal exit tunnel to attain complete folding to their native states. This escape process also frees up the ribosome tunnel for a new translation job. In this study, we investigate the impacts of energetic interactions between the ribosomal exit tunnel and nascent proteins on the protein escape process by molecular dynamics simulations using partially coarse-grained models that incorporate hydrophobic and electrostatic interactions of the ribosome tunnel of Haloarcula marismortui with nascent proteins. We find that, in general, attractive interactions slow down the protein escape process, whereas repulsive interactions speed it up. For the small globular proteins considered, the median escape time correlates with both the number of hydrophobic residues, N_h, and the net charge, Q, of a nascent protein. A correlation coefficient exceeding 0.96 is found for the relation between the median escape time and a combined quantity of N_h + 5.9Q, suggesting that it is ∼6 times more efficient to modulate the escape time by changing the total charge than the number of hydrophobic residues. The estimated median escape times are found in the submillisecond-to-millisecond range, indicating that the escape does not delay the ribosome recycling. For various types of the tunnel model, with and without hydrophobic and electrostatic interactions, the escape time distribution always follows a simple diffusion model that describes the escape process as a downhill drift of a Brownian particle, suggesting that nascent proteins escape along barrier-less pathways at the ribosome tunnel.

The ribosome modulates folding inside the ribosomal exit tunnel

Proteins commonly fold co-translationally at the ribosome, while the nascent chain emerges from the ribosomal exit tunnel. Protein domains that are sufficiently small can even fold while still located inside the tunnel. However, the effect of the tunnel on the folding dynamics of these domains is not well understood. Here, we combine optical tweezers with single-molecule FRET and molecular dynamics simulations to investigate folding of the small zinc-finger domain ADR1a inside and at the vestibule of the ribosomal tunnel. The tunnel is found to accelerate folding and stabilize the folded state, reminiscent of the effects of chaperonins. However, a simple mechanism involving stabilization by confinement does not explain the results. Instead, it appears that electrostatic interactions between the protein and ribosome contribute to the observed folding acceleration and stabilization of ADR1a.

Rational design strategies for functional reconstitution of plant cytochrome P450s in microbial systems

Plant natural products (NPs) are of pharmaceutical and agricultural significance, yet the low abundance is largely impeding the broad investigation and utilization. Microbial bioproduction is a promising alternative sourcing to plant NPs. Cytochrome P450s (CYPs) play an essential role in plant secondary metabolism, and functional reconstitution of plant CYPs in the microbial system is one of the major challenges in establishing efficient microbial plant NP bioproduction. In this review, we briefly summarized the recent progress in rational engineering strategies for enhanced activity of plant CYPs in Escherichia coli and Saccharomyces cerevisiae, two commonly used microbial hosts. We believe that in-depth foundational investigations on the native microenvironment of plant CYPs are necessary to adapt the microbial systems for more efficient functional reconstitution of plant CYPs.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Prion Protein Translocation Mechanism Revealed by Pulling Force Studies

The mammalian prion protein (PrP) engages with the ribosome-Sec61 translocation channel complex to generate different topological variants that are either physiological, or involved in neurodegenerative diseases. Here, we describe cotranslational folding and translocation mechanisms of PrP coupled to an Xbp1-based arrest peptide as folding sensor, to measure forces acting on PrP nascent chain. Our data reveal two main pulling events followed by a minor third one exerted on the nascent chains during their translocation. Using those force landscapes, we show that a specific sequence within an intrinsically disordered region, containing a polybasic and glycine-proline rich residues, modulates the second pulling event by interacting with TRAP complex. This work also delineates the sequence of events involved in generation of PrP toxic transmembrane topologies during its synthesis. Our results shed new insight into the folding of such a topological complex protein, where marginal pulling by the signal sequence, together with the flanking downstream sequence in the mature domain, primarily drives an overall inefficient translocation resulting in the nascent chain to adopt alternative topologies.

Chaperone action of a cofactor in protein folding

Previously, we have reported that ATP accelerates the folding and unfolding of Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is a glycolytic enzyme utilizing NAD⁺ as a cofactor. Because ATP and NAD⁺ share the ADP moiety, we hypothesized that NAD⁺ also accelerates the folding of GAPDH and that the common structural motif between ATP and NAD⁺ is responsible for the chaperone activity. After confirming that NAD⁺ indeed accelerates the folding of GAPDH, we examined the chaperone activity of the structural fragments of NAD⁺ (ADP, AMP, adenosine, and nicotinamide monophosphate). Our finding showed that ADP and AMP significantly speed up the folding of GAPDH, while adenosine and nicotinamide monophosphate do not. ADP and AMP also dramatically speed up the unfolding of GAPDH by selectively stabilizing a transition state in which GAPDH has a partially unfolded conformation. Similar to the previously reported effect of ATP on the equilibrium unfolding of GAPDH, a partially unfolded intermediate also accumulates in the presence of ADP and AMP. Based on the effect of the structural fragments of NAD⁺ on the folding of GAPDH, we identified that AMP is the structural determinant of the chaperone activity of ATP and NAD⁺ . Also, we propose a plausible model to explain how NAD⁺ accelerates the folding of GAPDH through a stepwise development of molecular interactions with the protein.

Mitochondria-Associated Proteostasis

Mitochondria are essential organelles in eukaryotes. Most mitochondrial proteins are encoded by the nuclear genome and translated in the cytosol. Nuclear-encoded mitochondrial proteins need to be imported, processed, folded, and assembled into their functional states. To maintain protein homeostasis (proteostasis), mitochondria are equipped with a distinct set of quality control machineries. Deficiencies in such systems lead to mitochondrial dysfunction, which is a hallmark of aging and many human diseases, such as neurodegenerative diseases, cardiovascular diseases, and cancer. In this review, we discuss the unique challenges and solutions of proteostasis in mitochondria. The import machinery coordinates with mitochondrial proteases and chaperones to maintain the mitochondrial proteome. Moreover, mitochondrial proteostasis depends on cytosolic protein quality control mechanisms during crises. In turn, mitochondria facilitate cytosolic proteostasis. Increasing evidence suggests that enhancing mitochondrial proteostasis may hold therapeutic potential to protect against protein aggregation-associated cellular defects.

Unraveling the Role of Silent Mutation in the ω-Subunit of Escherichia coli RNA Polymerase: Structure Transition Inhibits Transcription

The bacterial RNA polymerase is a multi-subunit enzyme complex composed of six subunits, α2ββ'σω. The function of this enzyme is to transcribe the DNA base sequence to the RNA intermediate, which is ultimately translated to protein. Though the contribution of each subunit in RNA synthesis has been clearly elucidated, the role of the smallest ω-subunit is still unclear despite several studies. Recently, a study on a dominant negative mutant of rpoZ has been reported in which the mutant was shown to render the RNA polymerase defective in transcription initiation (ω6, N60D) and gave an insight on the function of ω in RNA polymerase. Serendipitously, we also obtained a silent mutant, and the mutant was found to be lethal during the isolation of toxic mutants. The primary focus of this study is to understand the mechanistic details of this lethality. Isolated ω shows a predominantly unstructured circular dichroism profile and becomes α-helical in the enzyme complex. This structural transition is perhaps the reason for this lack of function. Subsequently, we generated several silent mutants of ω to investigate the role of codon bias and the effect of rare codons with respect to their position in rpoZ. Not all silent mutations affect the structure. RNA polymerase when reconstituted with structurally altered silent mutants of ω is transcriptionally inactive. The CodonPlus strain, which has surplus tRNA, was used to assess for the rescue of the phenotype in lethal silent mutants.

Versatile format of minichaperone-based protein fusion system

Hydrophobic recombinant proteins often tend to aggregate upon expression into inclusion bodies and are difficult to refold. Producing them in soluble forms constitutes a common bottleneck problem. A fusion system for production of insoluble hydrophobic proteins in soluble stable forms with thermophilic minichaperone, GroEL apical domain (GrAD) as a carrier, has recently been developed. To provide the utmost flexibility of the system for interactions between the carrier and various target protein moieties a strategy of making permutated protein variants by gene engineering has been applied: the original N- and C-termini of the minichaperone were linked together by a polypeptide linker and new N- and C-termini were made at desired parts of the protein surface. Two permutated GrAD forms were created and analyzed. Constructs of GrAD and both of its permutated forms fused with the initially insoluble N-terminal fragment of hepatitis C virus' E2 protein were tested. Expressed fusions formed inclusion bodies. After denaturation, all fusions were completely renatured in stable soluble forms. A variety of permutated GrAD variants can be created. The versatile format of the system provides opportunities for choosing an optimal pair between particular target protein moiety and the best-suited original or specific permutated carrier.

Nature and Regulation of Protein Folding on the Ribosome

Co-translational protein folding is an essential process by which cells ensure the safe and efficient production and assembly of new proteins in their functional native states following biosynthesis on the ribosome. In this review, we describe recent progress in probing the changes during protein synthesis of the free energy landscapes that underlie co-translational folding and discuss the critical coupling between these landscapes and the rate of translation that ultimately determines the success or otherwise of the folding process. Recent developments have revealed a variety of mechanisms by which both folding and translation can be modulated or regulated, and we discuss how these effects are utilised by the cell to optimise the outcome of protein biosynthesis.

The Ribosome Cooperates with a Chaperone to Guide Multi-domain Protein Folding

Multi-domain proteins, containing several structural units within a single polypeptide, constitute a large fraction of all proteomes. Co-translational folding is assumed to simplify the conformational search problem for large proteins, but the events leading to correctly folded, functional structures remain poorly characterized. Similarly, how the ribosome and molecular chaperones promote efficient folding remains obscure. Using optical tweezers, we have dissected early folding events of nascent elongation factor G, a multi-domain protein that requires chaperones for folding. The ribosome and the chaperone trigger factor reduce inter-domain misfolding, permitting folding of the N-terminal G-domain. Successful completion of this step is a crucial prerequisite for folding of the next domain. Unexpectedly, co-translational folding does not proceed unidirectionally; emerging unfolded polypeptide can denature an already-folded domain. Trigger factor, but not the ribosome, protects against denaturation. The chaperone thus serves a previously unappreciated function, helping multi-domain proteins overcome inherent challenges during co-translational folding.

Acidic pH Mediates Changes in Antigenic and Oligomeric Conformation of Herpes Simplex Virus gB and Is a Determinant of Cell-Specific Entry

Herpes simplex virus (HSV) is an important human pathogen with a high worldwide seroprevalence. HSV enters epithelial cells, the primary site of infection, by a low-pH pathway. HSV glycoprotein B (gB) undergoes low pH-induced conformational changes, which are thought to drive membrane fusion. When neutralized back to physiological pH, these changes become reversible. Here, HSV-infected cells were subjected to short pulses of radiolabeling, followed by immunoprecipitation with a panel of gB monoclonal antibodies (MAbs), demonstrating that gB folds and oligomerizes rapidly and cotranslationally in the endoplasmic reticulum. Full-length gB from transfected cells underwent low-pH-triggered changes in oligomeric conformation in the absence of other viral proteins. MAbs to gB neutralized HSV entry into cells regardless of the pH dependence of the entry pathway, suggesting a conservation of gB function in distinct fusion mechanisms. The combination of heat and acidic pH triggered irreversible changes in the antigenic conformation of the gB fusion domain, while changes in the gB oligomer remained reversible. An elevated temperature alone was not sufficient to induce gB conformational change. Together, these results shed light on the conformation and function of the HSV-1 gB oligomer, which serves as part of the core fusion machinery during viral entry.IMPORTANCE Herpes simplex virus (HSV) causes infection of the mouth, skin, eyes, and genitals and establishes lifelong latency in humans. gB is conserved among all herpesviruses. HSV gB undergoes reversible conformational changes following exposure to acidic pH which are thought to mediate fusion and entry into epithelial cells. Here, we identified cotranslational folding and oligomerization of newly synthesized gB. A panel of antibodies to gB blocked both low-pH and pH-neutral entry of HSV, suggesting conserved conformational changes in gB regardless of cell entry route. Changes in HSV gB conformation were not triggered by increased temperature alone, in contrast to results with EBV gB. Acid pH-induced changes in the oligomeric conformation of gB are related but distinct from pH-triggered changes in gB antigenic conformation. These results highlight critical aspects of the class III fusion protein, gB, and inform strategies to block HSV infection at the level of fusion and entry.

Protein escape at the ribosomal exit tunnel: Effects of native interactions, tunnel length, and macromolecular crowding

How fast a post-translational nascent protein escapes from the ribosomal exit tunnel is relevant to its folding and protection against aggregation. Here, by using Langevin molecular dynamics, we show that non-local native interactions help decrease the escape time, and foldable proteins generally escape much faster than same-length, self-repulsive homopolymers at low temperatures. The escape process, however, is slowed down by the local interactions that stabilize the α-helices. The escape time is found to increase with both the tunnel length and the concentration of macromolecular crowders outside the tunnel. We show that a simple diffusion model described by the Smoluchowski equation with an effective linear potential can be used to map out the escape time distribution for various tunnel lengths and various crowder concentrations. The consistency between the simulation data and the diffusion model, however, is found only for the tunnel length smaller than a crossover length of 90 Å-110 Å, above which the escape time increases much faster with the tunnel length. It is suggested that the length of ribosomal exit tunnel has been selected by evolution to facilitate both the efficient folding and the efficient escape of single-domain proteins. We show that macromolecular crowders lead to an increase in the escape time, and attractive crowders are unfavorable for the folding of nascent polypeptide.

Protein Sequences Recapitulate Genetic Code Evolution

Several hypotheses predict ranks of amino acid assignments to genetic code's codons. Analyses here show that average positions of amino acid species in proteins correspond to assignment ranks, in particular as predicted by Juke's neutral mutation hypothesis for codon assignments. In all tested protein groups, including co- and post-translationally folding proteins, 'recent' amino acids are on average closer to gene 5' extremities than 'ancient' ones. Analyses of pairwise residue contact energies matrices suggest that early amino acids stereochemically selected late ones that stablilize residue interactions within protein cores, presumably producing 5'-late-to-3'-early amino acid protein sequence gradients. The gradient might reduce protein misfolding, also after mutations, extending principles of neutral mutations to protein folding. Presumably, in self-perpetuating and self-correcting systems like the genetic code, initial conditions produce similarities between evolution of the process (the genetic code) and 'ontogeny' of resulting structures (here proteins), producing apparent teleonomy between process and product.

Co-translational protein targeting facilitates centrosomal recruitment of PCNT during centrosome maturation in vertebrates

As microtubule-organizing centers of animal cells, centrosomes guide the formation of the bipolar spindle that segregates chromosomes during mitosis. At mitosis onset, centrosomes maximize microtubule-organizing activity by rapidly expanding the pericentriolar material (PCM). This process is in part driven by the large PCM protein pericentrin (PCNT), as its level increases at the PCM and helps recruit additional PCM components. However, the mechanism underlying the timely centrosomal enrichment of PCNT remains unclear. Here, we show that PCNT is delivered co-translationally to centrosomes during early mitosis by cytoplasmic dynein, as evidenced by centrosomal enrichment of PCNT mRNA, its translation near centrosomes, and requirement of intact polysomes for PCNT mRNA localization. Additionally, the microtubule minus-end regulator, ASPM, is also targeted co-translationally to mitotic spindle poles. Together, these findings suggest that co-translational targeting of cytoplasmic proteins to specific subcellular destinations may be a generalized protein targeting mechanism.

Before a cell divides, it creates a copy of its genetic material (DNA) and evenly distributes it between the new ‘daughter’ cells with the help of a complex called the mitotic spindle. This complex is made of long cable-like protein chains called microtubules.

To ensure that each daughter cell receives an equal amount of DNA, structures known as centrosomes organize the microtubules during the division process. Centrosomes have two rigid cores, called centrioles, which are surrounded by a matrix of proteins called the pericentriolar material. It is from this material that the microtubules are organized.

The pericentriolar material is a dynamic structure and changes its size by assembling and disassembling its protein components. The larger the pericentriolar material, the more microtubules can form. Before a cell divides, it rapidly expands in a process called centrosome maturation. A protein called pericentrin initiates the maturation by helping to recruit other proteins to the centrosome. Pericentrin molecules are large, and it takes the cell between 10 and 20 minutes to make each one. Nevertheless, the cell can produce and deliver large quantities of pericentrin to the centrosome in a matter of minutes. We do not yet know how this happens.

To investigate this further, Sepulveda, Antkowiak, Brust-Mascher et al. used advanced microscopy to study zebrafish embryos and human cells grown in the laboratory. The results showed that cells build and transport pericentrin at the same time. Cells use messenger RNA molecules as templates to build proteins. These feed into protein factories called ribosomes, which assemble the building blocks in the correct order. Rather than waiting for the pericentrin production to finish, the cell moves the active factories to the centrosome with the help of a molecular motor called dynein. By the time the pericentrin molecules are completely made by ribosomes, they are already at the centrosome, ready to help with the recruitment of other proteins during centrosome maturation.

These findings improve our understanding of centrosome maturation. The next step is to find out how the cell coordinates this process with the recruitment of other proteins to the centrosome. It is also possible that the cell uses similar processes to deliver other proteins to different parts of the cell.

Sequential search leads to faster, more efficient fragment-based de novo protein structure prediction

Motivation

Most current de novo structure prediction methods randomly sample protein conformations and thus require large amounts of computational resource. Here, we consider a sequential sampling strategy, building on ideas from recent experimental work which shows that many proteins fold cotranslationally.

Results

We have investigated whether a pseudo-greedy search approach, which begins sequentially from one of the termini, can improve the performance and accuracy of de novo protein structure prediction. We observed that our sequential approach converges when fewer than 20 000 decoys have been produced, fewer than commonly expected. Using our software, SAINT2, we also compared the run time and quality of models produced in a sequential fashion against a standard, non-sequential approach. Sequential prediction produces an individual decoy 1.5-2.5 times faster than non-sequential prediction. When considering the quality of the best model, sequential prediction led to a better model being produced for 31 out of 41 soluble protein validation cases and for 18 out of 24 transmembrane protein cases. Correct models (TM-Score > 0.5) were produced for 29 of these cases by the sequential mode and for only 22 by the non-sequential mode. Our comparison reveals that a sequential search strategy can be used to drastically reduce computational time of de novo protein structure prediction and improve accuracy.

Availability and implementation

Data are available for download from: http://opig.stats.ox.ac.uk/resources. SAINT2 is available for download from: https://github.com/sauloho/SAINT2.

Contact

saulo.deoliveira@dtc.ox.ac.uk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Does mRNA structure contain genetic information for regulating co-translational protein folding?

Currently many facets of genetic information are illdefined. In particular, how protein folding is genetically regulated has been a long-standing issue for genetics and protein biology. And a generic mechanistic model with supports of genomic data is still lacking. Recent technological advances have enabled much needed genome-wide experiments. While putting the effect of codon optimality on debate, these studies have supplied mounting evidence suggesting a role of mRNA structure in the regulation of protein folding by modulating translational elongation rate. In conjunctions with previous theories, this mechanistic model of protein folding guided by mRNA structure shall expand our understandings of genetic information and offer new insights into various biomedical puzzles.

Genetic Code Optimization for Cotranslational Protein Folding: Codon Directional Asymmetry Correlates with Antiparallel Betasheets, tRNA Synthetase Classes

A new codon property, codon directional asymmetry in nucleotide content (CDA), reveals a biologically meaningful genetic code dimension: palindromic codons (first and last nucleotides identical, codon structure XZX) are symmetric (CDA = 0), codons with structures ZXX/XXZ are 5'/3' asymmetric (CDA = - 1/1; CDA = - 0.5/0.5 if Z and X are both purines or both pyrimidines, assigning negative/positive (-/+) signs is an arbitrary convention). Negative/positive CDAs associate with (a) Fujimoto's tetrahedral codon stereo-table; (b) tRNA synthetase class I/II (aminoacylate the 2'/3' hydroxyl group of the tRNA's last ribose, respectively); and (c) high/low antiparallel (not parallel) betasheet conformation parameters. Preliminary results suggest CDA-whole organism associations (body temperature, developmental stability, lifespan). Presumably, CDA impacts spatial kinetics of codon-anticodon interactions, affecting cotranslational protein folding. Some synonymous codons have opposite CDA sign (alanine, leucine, serine, and valine), putatively explaining how synonymous mutations sometimes affect protein function. Correlations between CDA and tRNA synthetase classes are weaker than between CDA and antiparallel betasheet conformation parameters. This effect is stronger for mitochondrial genetic codes, and potentially drives mitochondrial codon-amino acid reassignments. CDA reveals information ruling nucleotide-protein relations embedded in reversed (not reverse-complement) sequences (5'-ZXX-3'/5'-XXZ-3').

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

Rare ER protein misfolding-mistrafficking disorders: Therapeutic developments

The presence of a functional protein at the appropriate location in the cell is the result of the processes of transcription, translation, folding and trafficking to the correct destination. There are numerous diseases that are caused by protein misfolding, mainly due to mutations in the respective gene. The consequences of this misfolding may be that proteins effectively lose their function, either by being removed by the cellular quality control machinery or by accumulating at the incorrect intracellular or extracellular location. A number of mutations that lead to protein misfolding and affect trafficking to the final destination, e.g. Cystic fibrosis, Wilson's disease, and Progressive Familial Intrahepatic 1 cholestasis, result in proteins that retain partial function if their folding and trafficking is restored either by molecular or pharmacological means. In this review, we discuss several mutant proteins within this class of misfolding diseases and provide an update on the status of molecular and therapeutic developments and potential therapeutic strategies being developed to counter these diseases.

High Yield of Recombinant Protein in Shaken E. coli Cultures with Enzymatic Glucose Release Medium EnPresso B

Expression of recombinant proteins in sufficient quantities is essential for protein structure-function studies. The most commonly used method for recombinant protein production is overexpression in E. coli cultures. However, producing high yields of functional proteins in E. coli can be a challenge in conventional shaken cultures. This is often due to nonoptimal growth conditions, which result in low cell yields and insoluble or incorrectly folded target protein. To overcome the shortcomings of shake flask cultivation, we present a cultivation method based on enzymatic glucose delivery. This system mimics the fed-batch principle used in bioreactor cultivations and provides high yields of biomass and recombinant proteins in shaken cultivations.

Prion-Associated Toxicity is Rescued by Elimination of Cotranslational Chaperones

The nascent polypeptide-associated complex (NAC) is a highly conserved but poorly characterized triad of proteins that bind near the ribosome exit tunnel. The NAC is the first cotranslational factor to bind to polypeptides and assist with their proper folding. Surprisingly, we found that deletion of NAC subunits in Saccharomyces cerevisiae rescues toxicity associated with the strong [PSI+] prion. This counterintuitive finding can be explained by changes in chaperone balance and distribution whereby the folding of the prion protein is improved and the prion is rendered nontoxic. In particular, the ribosome-associated Hsp70 Ssb is redistributed away from Sup35 prion aggregates to the nascent chains, leading to an array of aggregation phenotypes that can mimic both overexpression and deletion of Ssb. This toxicity rescue demonstrates that chaperone modification can block key steps of the prion life cycle and has exciting implications for potential treatment of many human protein conformational disorders.

Misfolded proteins can be toxic to cells, causing pathologies such as Alzheimer’s disease, Parkinson’s disease, prion diseases, and ALS. One mechanism by which organisms combat protein misfolding involves molecular chaperones, proteins that help other proteins fold correctly. Here, we describe a novel role for a family of chaperones called the nascent polypeptide-associated complex (NAC). The NAC is a group of proteins that exist in all multicellular organisms, yet we do not fully understand its function. Using yeast as a model system, we have found that deletion of NAC subunits can reduce the toxicity associated with misfolded proteins. This work has implications for human protein misfolding diseases, as modulation of the NAC may present a viable therapeutic avenue by which to slow the progression of neurodegeneration and other protein conformational disorders.

Folding and escape of nascent proteins at ribosomal exit tunnel

We investigate the interplay between post-translational folding and escape of two small single-domain proteins at the ribosomal exit tunnel by using Langevin dynamics with coarse-grained models. It is shown that at temperatures lower or near the temperature of the fastest folding, folding proceeds concomitantly with the escape process, resulting in vectorial folding and enhancement of foldability of nascent proteins. The concomitance between the two processes, however, deteriorates as temperature increases. Our folding simulations as well as free energy calculation by using umbrella sampling show that, at low temperatures, folding at the tunnel follows one or two specific pathways without kinetic traps. It is shown that the escape time can be mapped to a one-dimensional diffusion model with two different regimes for temperatures above and below the folding transition temperature. Attractive interactions between amino acids and attractive sites on the tunnel wall lead to a free energy barrier along the escape route of the protein. It is suggested that this barrier slows down the escape process and consequently promotes correct folding of the released nascent protein.

Folding and assembly of the large molecular machine Hsp90 studied in single-molecule experiments

Folding of small proteins often occurs in a two-state manner and is well understood both experimentally and theoretically. However, many proteins are much larger and often populate misfolded states, complicating their folding process significantly. Here we study the complete folding and assembly process of the 1,418 amino acid, dimeric chaperone Hsp90 using single-molecule optical tweezers. Although the isolated C-terminal domain shows two-state folding, we find that the isolated N-terminal as well as the middle domain populate ensembles of fast-forming, misfolded states. These intradomain misfolds slow down folding by an order of magnitude. Modeling folding as a competition between productive and misfolding pathways allows us to fully describe the folding kinetics. Beyond intradomain misfolding, folding of the full-length protein is further slowed by the formation of interdomain misfolds, suggesting that with growing chain lengths, such misfolds will dominate folding kinetics. Interestingly, we find that small stretching forces applied to the chain can accelerate folding by preventing the formation of cross-domain misfolding intermediates by leading the protein along productive pathways to the native state. The same effect is achieved by cotranslational folding at the ribosome in vivo.

A minichaperone-based fusion system for producing insoluble proteins in soluble stable forms

We have developed a fusion system for reliable production of insoluble hydrophobic proteins in soluble stable forms. A carrier is thermophilic minichaperone, GroEL apical domain (GrAD), a 15 kDa monomer able to bind diverse protein substrates. The Met-less variant of GrAD has been made for further convenient use of Met-specific CNBr chemical cleavage, if desired. The Met-less GrAD retained stability and solubility of the original protein. Target polypeptides can be fused to either C-terminus or N-terminus of GrAD. The system has been tested with two unrelated insoluble proteins fused to the C-terminus of GrAD. One of the proteins was also fused to GrAD N-terminus. The fusions formed inclusion bodies at 25°C and above and were partly soluble only at lower expression temperatures. Most importantly, however, after denaturation in urea, all fusions without exception were completely renatured in soluble stable forms that safely survived freezing-thawing as well as lyophilization. All fusions for both tested target proteins retained solubility at high concentrations for days. Functional analysis revealed that a target protein may retain functionality in the fusion. Convenience features include potential thermostability of GrAD fusions, capacity for chemical and enzymatic cleavage of a target and His6 tag for purification.

Proteome Analysis Unravels Mechanism Underling the Embryogenesis of the Honeybee Drone and Its Divergence with the Worker (Apis mellifera lingustica)

The worker and drone bees each contain a separate diploid and haploid genetic makeup, respectively. Mechanisms regulating the embryogenesis of the drone and its mechanistic difference with the worker are still poorly understood. The proteomes of the two embryos at three time-points throughout development were analyzed by applying mass spectrometry-based proteomics. We identified 2788 and 2840 proteins in the worker and drone embryos, respectively. The age-dependent proteome driving the drone embryogenesis generally follows the worker's. The two embryos however evolve a distinct proteome setting to prime their respective embryogenesis. The strongly expressed proteins and pathways related to transcriptional-translational machinery and morphogenesis at 24 h drone embryo relative to the worker, illustrating the earlier occurrence of morphogenesis in the drone than worker. These morphogenesis differences remain through to the middle-late stage in the two embryos. The two embryos employ distinct antioxidant mechanisms coinciding with the temporal-difference organogenesis. The drone embryo's strongly expressed cytoskeletal proteins signify key roles to match its large body size. The RNAi induced knockdown of the ribosomal protein offers evidence for the functional investigation of gene regulating of honeybee embryogenesis. The data significantly expand novel regulatory mechanisms governing the embryogenesis, which is potentially important for honeybee and other insects.

How Co-translational Folding of Multi-domain Protein Is Affected by Elongation Schedule: Molecular Simulations

Co-translational folding (CTF) facilitates correct folding in vivo, but its precise mechanism remains elusive. For the CTF of a three-domain protein SufI, it was reported that the translational attenuation is obligatory to acquire the functional state. Here, to gain structural insights on the underlying mechanisms, we performed comparative molecular simulations of SufI that mimic CTF as well as refolding schemes. A CTF scheme that relied on a codon-based prediction of translational rates exhibited folding probability markedly higher than that by the refolding scheme. When the CTF schedule is speeded up, the success rate dropped. These agree with experiments. Structural investigation clarified that misfolding of the middle domain was much more frequent in the refolding scheme than that in the codon-based CTF scheme. The middle domain is less stable and can fold via interactions with the folded N-terminal domain. Folding pathway networks showed the codon-based CTF gives narrower pathways to the native state than the refolding scheme.

Proteins are synthesized in vivo by ribosome from their N-termini. When N-terminal fragments of nascent proteins get out of the ribosome exit, they start folding, which is called co-translational folding. It has been suggested that well-scheduled co-translational folding schemes would facilitate correct acquisition of their native structures for some multi-domain proteins. In particular, an un-ambiguous experiment was recently reported for a model protein, SufI where pauses at certain positions in the translational elongation are obligatory for efficient folding. Here, for the first time to our knowledge, we performed molecular dynamics simulations of SufI with co-translational folding as well as re-folding schemes. We found a co-translational folding shceme with rare codon-based pauses indeed increased the success ratio of folding, which is consistent with recent experiments. On top, molecular simulations provided much of structural insights on the folding routes and misfolding in the case of re-folding scheme. This explains why pauses in the translational elongation rescue SufI from misfolding.

Analysis of transcription factor networks using IVV method

We have developed a simple and totally in vitro selection procedure based on cell-free cotranslation using a highly stable and efficient in vitro virus (IVV). Cell-free cotranslation of tagged bait and prey proteins is advantageous for the formation of protein complexes and allows high-throughput analysis of protein-protein interactions (PPI) as a result of providing in vitro instead of in vivo preparation of bait proteins. The use of plural selection rounds and a two-step purification of the IVV selection, followed by in vitro post-selection, is advantageous for decreasing false positives. This simple IVV selection system based on cell-free cotranslation is applicable to high-throughput and comprehensive analysis of transcription factor networks.

CHIP: a co-chaperone for degradation by the proteasome

Protein homeostasis relies on a balance between protein folding and protein degradation. Molecular chaperones like Hsp70 and Hsp90 fulfil well-defined roles in protein folding and conformational stability via ATP dependent reaction cycles. These folding cycles are controlled by associations with a cohort of non-client protein co-chaperones, such as Hop, p23 and Aha1. Pro-folding co-chaperones facilitate the transit of the client protein through the chaperone mediated folding process. However, chaperones are also involved in ubiquitin-mediated proteasomal degradation of client proteins. Similar to folding complexes, the ability of chaperones to mediate protein degradation is regulated by co-chaperones, such as the C terminal Hsp70 binding protein (CHIP). CHIP binds to Hsp70 and Hsp90 chaperones through its tetratricopeptide repeat (TPR) domain and functions as an E3 ubiquitin ligase using a modified RING finger domain (U-box). This unique combination of domains effectively allows CHIP to network chaperone complexes to the ubiquitin-proteasome system. This chapter reviews the current understanding of CHIP as a co-chaperone that switches Hsp70/Hsp90 chaperone complexes from protein folding to protein degradation.

Ligand binding to a high-energy partially unfolded protein

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP(+) as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP(+) and found that NADP(+) binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP(+) is diminished upon partial unfolding. Based on known crystallographic structures of NADP(+) -bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine-binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP(+) . Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high-energy non-native forms.

The effect of tRNA levels on decoding times of mRNA codons

The possible effect of transfer ribonucleic acid (tRNA) concentrations on codons decoding time is a fundamental biomedical research question; however, due to a large number of variables affecting this process and the non-direct relation between them, a conclusive answer to this question has eluded so far researchers in the field. In this study, we perform a novel analysis of the ribosome profiling data of four organisms which enables ranking the decoding times of different codons while filtering translational phenomena such as experimental biases, extreme ribosomal pauses and ribosome traffic jams. Based on this filtering, we show for the first time that there is a significant correlation between tRNA concentrations and the codons estimated decoding time both in prokaryotes and in eukaryotes in natural conditions (−0.38 to −0.66, all P values <0.006); in addition, we show that when considering tRNA concentrations, codons decoding times are not correlated with aminoacyl-tRNA levels. The reported results support the conjecture that translation efficiency is directly influenced by the tRNA levels in the cell. Thus, they should help to understand the evolution of synonymous aspects of coding sequences via the adaptation of their codons to the tRNA pool.

Nanosecond dynamics of calmodulin and ribosome-bound nascent chains studied by time-resolved fluorescence anisotropy

We report a time-resolved fluorescence anisotropy study of ribosome-bound nascent chains (RNCs) of calmodulin (CaM), a prototypical member of the EF-hand family of calcium-sensing proteins. As shown in numerous studies, in vitro protein refolding can differ substantially from biosynthetic protein folding, which takes place cotranslationally and depends on the rate of polypeptide chain elongation. A challenge in this respect is to characterize the adopted conformations of nascent chains before their release from the ribosome. CaM RNCs (full-length, half-length, and first EF-hand only) were synthesized in vitro. All constructs contained a tetracysteine motif site-specifically incorporated in the first N-terminal helix; this motif is known to react with FlAsH, a biarsenic fluorescein derivative. As the dye is rotationally locked to this helix, we characterized the structural properties and folding states of polypeptide chains tethered to ribosomes and compared these with released chains. Importantly, we observed decelerated tumbling motions of ribosome-tethered and partially folded nascent chains, compared to released chains. This indicates a pronounced interaction between nascent chains and the ribosome surface, and might reflect chaperone activity of the ribosome.

Protein folding on the ribosome studied using NMR spectroscopy

NMR spectroscopy is a powerful tool for the investigation of protein folding and misfolding, providing a characterization of molecular structure, dynamics and exchange processes, across a very wide range of timescales and with near atomic resolution. In recent years NMR methods have also been developed to study protein folding as it might occur within the cell, in a de novo manner, by observing the folding of nascent polypeptides in the process of emerging from the ribosome during synthesis. Despite the 2.3 MDa molecular weight of the bacterial 70S ribosome, many nascent polypeptides, and some ribosomal proteins, have sufficient local flexibility that sharp resonances may be observed in solution-state NMR spectra. In providing information on dynamic regions of the structure, NMR spectroscopy is therefore highly complementary to alternative methods such as X-ray crystallography and cryo-electron microscopy, which have successfully characterized the rigid core of the ribosome particle. However, the low working concentrations and limited sample stability associated with ribosome-nascent chain complexes means that such studies still present significant technical challenges to the NMR spectroscopist. This review will discuss the progress that has been made in this area, surveying all NMR studies that have been published to date, and with a particular focus on strategies for improving experimental sensitivity.

Quantifying elongation rhythm during full-length protein synthesis

Pauses regulate the rhythm of ribosomal protein synthesis. Mutations disrupting even minor pauses can give rise to improperly formed proteins and human disease. Such minor pauses are difficult to characterize by ensemble methods, but can be readily examined by single-molecule (sm) approaches. Here we use smFRET to carry out real-time monitoring of the expression of a full-length protein, the green fluorescent protein variant Emerald GFP. We demonstrate significant correlations between measured elongation rates and codon and isoacceptor tRNA usage, and provide a quantitative estimate of the effect on elongation rate of replacing a codon recognizing an abundant tRNA with a synonymous codon cognate to a rarer tRNA. Our results suggest that tRNA selection plays an important general role in modulating the rates and rhythms of protein synthesis, potentially influencing simultaneous co-translational processes such as folding and chemical modification.

Nutrient signaling in protein homeostasis: an increase in quantity at the expense of quality

The discovery that rapamycin extends the life span of diverse organisms has triggered many studies aimed at identifying the underlying molecular mechanisms. Mammalian target of rapamycin complex 1 (mTORC1) regulates cell growth and may regulate organismal aging by controlling mRNA translation. However, how inhibiting mTORC1 and decreasing protein synthesis can extend life span remains an unresolved issue. We showed that constitutively active mTORC1 signaling increased general protein synthesis but unexpectedly reduced the quality of newly synthesized polypeptides. We demonstrated that constitutively active mTORC1 decreased translation fidelity by increasing the speed of ribosomal elongation. Conversely, rapamycin treatment restored the quality of newly synthesized polypeptides mainly by slowing the rate of ribosomal elongation. We also found distinct roles for mTORC1 downstream targets in maintaining protein homeostasis. Loss of S6 kinases, but not 4E-BP family proteins, which are both involved in regulation of translation, attenuated the effects of rapamycin on the quality of newly translated proteins. Our results reveal a mechanistic connection between mTORC1 and protein quality, highlighting the central role of nutrient signaling in growth and aging.

Coding nucleic acids are chaperons for protein folding: a novel theory of protein folding

The arguments for nucleic acid chaperons are reviewed and three new lines of evidence are added. (1) It was found that amino acids encoded by codons in short nucleic acid loops frequently form turns and helices in the corresponding protein structures. (2) The amino acids encoded by partially complementary (1st and 3rd nucleotides) codons are more frequently co-located in the encoded proteins than expected by chance. (3) There are significant correlations between thermodynamic changes (ddG) caused by codon mutations in nucleic acids and the thermodynamic changes caused by the corresponding amino acid mutations in the encoded proteins. We conclude that the concept of the Proteomic Code and nucleic acid chaperons seems correct from the bioinformatics point of view, and we expect to see direct biochemical experiments and evidence in the near future.

Mechanisms of CFTR Folding at the Endoplasmic Reticulum

In the past decade much has been learned about how Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) folds and misfolds as the etiologic cause of cystic fibrosis (CF). CFTR folding is complex and hierarchical, takes place in multiple cellular compartments and physical environments, and involves several large networks of folding machineries. Insertion of transmembrane (TM) segments into the endoplasmic reticulum (ER) membrane and tertiary folding of cytosolic domains begin cotranslationally as the nascent polypeptide emerges from the ribosome, whereas posttranslational folding establishes critical domain-domain contacts needed to form a physiologically stable structure. Within the membrane, N- and C-terminal TM helices are sorted into bundles that project from the cytosol to form docking sites for nucleotide binding domains, NBD1 and NBD2, which in turn form a sandwich dimer for ATP binding. While tertiary folding is required for domain assembly, proper domain assembly also reciprocally affects folding of individual domains analogous to a jig-saw puzzle wherein the structure of each interlocking piece influences its neighbors. Superimposed on this process is an elaborate proteostatic network of cellular chaperones and folding machineries that facilitate the timing and coordination of specific folding steps in and across the ER membrane. While the details of this process require further refinement, we finally have a useful framework to understand key folding defect(s) caused by ΔF508 that provides a molecular target(s) for the next generation of CFTR small molecule correctors aimed at the specific defect present in the majority of CF patients.

Sensitive measurement of single-nucleotide polymorphism-induced changes of RNA conformation: application to disease studies

Single-nucleotide polymorphisms (SNPs) are often linked to critical phenotypes such as diseases or responses to vaccines, medications and environmental factors. However, the specific molecular mechanisms by which a causal SNP acts is usually not obvious. Changes in RNA secondary structure emerge as a possible explanation necessitating the development of methods to measure the impact of single-nucleotide variation on RNA structure. Despite the recognition of the importance of considering the changes in Boltzmann ensemble of RNA conformers in this context, a formal method to perform directly such comparison was lacking. Here, we solved this problem and designed an efficient method to compute the relative entropy between the Boltzmann ensembles of the native and a mutant structure. On the basis of this theoretical progress, we developed a software tool, remuRNA, and investigated examples of its application. Comparing the impact of common SNPs naturally occurring in populations with the impact of random point mutations, we found that structural changes introduced by common SNPs are smaller than those introduced by random point mutations. This suggests a natural selection against mutations that significantly change RNA structure and demonstrates, surprisingly, that randomly inserted point mutations provide inadequate estimation of random mutations effects. Subsequently, we applied remuRNA to determine which of the disease-associated non-coding SNPs are potentially related to RNA structural changes.