Nothing Special   »   [go: up one dir, main page]
Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

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

Abstract

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.

Author Summary

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.

Citation: Tanaka T, Hori N, Takada S (2015) How Co-translational Folding of Multi-domain Protein Is Affected by Elongation Schedule: Molecular Simulations. PLoS Comput Biol 11(7): e1004356. https://doi.org/10.1371/journal.pcbi.1004356

Editor: Alexander MacKerell, University of Maryland, UNITED STATES

Received: January 16, 2015; Accepted: May 22, 2015; Published: July 9, 2015

Copyright: © 2015 Tanaka et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: ST was supported by Grant-in-Aid for Scientific Research (A) of the Japan Society of Promotion of Science (JSPS), the Grant-in-Aid for Scientific Research on Innovative Area of the Ministry of Education Culture, Sports, Science, and Technology (MEXT), and the Strategic Programs for Innovative Research "Supercomputational Life Science" of MEXT. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

While in vitro folding dynamics of single-domain proteins has been relatively well understood by now[1,2], several additional factors make in vivo protein folding much more difficult to characterize. About 70% of proteins have multiple domains and inter-domain interactions often cause many metastable intermediates and can hamper folding to the native states [3,4]. Cellular environment is highly crowded by macromolecules, which affects folding kinetics and could cause aggregation [5–7]. To circumvent some of these difficulties, several types of molecular chaperones facilitate folding [8]. During protein synthesis in ribosome, nascent polypeptides start folding co-translationally [9].

Co-translational folding (CTF) has been suggested for in vivo folding mechanism since 1960’s [10] and there is no room to doubt its relevance both in bacteria and in eukaryotic cells [11]. Many elements in the CTF have been characterized [12]. First of all, many proteins, once denatured in a test tube, do not refold with high probability, whereas they fold in the CTF condition. Thus, as a rule of thumb, the CTF condition facilitates correct folding of many proteins [13,14]. Ribosome is not just a machine for synthesis, but also helps folding of nascent chains at the exit tunnel and on the surface [15]. The translation elongation is not at uniform rate, but there are some regions on mRNA where the elongation is markedly slowed down [16–18].

This so-called elongation attenuation can be realized by a few mechanisms. Most notably, for a given codon, the elongation rate is affected by its cognate tRNA binding kinetics, thus depending on the concentration of the cognate tRNA [19]. The concentration of cognate tRNAs are highly correlated with the frequencies of codon usage for each of species. There are some codons, of which the cognate tRNAs have markedly low concentration[19]. These rare codons sometime appeared in mRNA as a cluster, which often leads to translational attenuation. On top, some portions of mRNA form partial secondary structures, which may slow down the elongation contributing to the elongation attenuation as well[20]. It was anticipated that the locations of the attenuation might have evolved to facilitate the CTF. Some of them appear near domain boundaries of multi-domain proteins [21]. By synonymous substitution of rare codons, one can speed up the translation elongation at a certain position without changing amino acid sequence, which led to reduce or impair functions and/or protease resistance for some proteins, such as an acetyl-transferase [22] and SufI [16].

SufI in E. coli was recently used to test the role of translational attenuation in the CTF [16]. SufI, an about 450-residue protein, is made of three domains; N- (blue in Fig 1A), M- (green), and C- domains (red), in order. Zhang et al. first identified three clusters of rare codons, two of which indeed exhibited elongation attenuation [16]. Synonymous substitutions of some rare codons in these regions led to reduction or impair of protease resistance. Separately, using a cell-free system, they also increased the concentrations of the corresponding tRNAs, which showed the similar results to the above synonymous substitution experiment. It should also be noted that, they found no interactions of SufI with molecular chaperones. Thus, these experiments provide us an unambiguous evidence of biological importance of the elongation attenuation for efficient folding in the CTF condition.

thumbnail
Download:
Fig 1. Simulation setup.

A) Strucuture of SufI. N-, M-, and C-domains are depicted in blue, green, and red. Linkers that connet two domains are depicted in yellow. B) The codon-based elongation rate by Spencer et al’s algorithm. A threshold is introduced. The region where the elongation rate is slower than the threshold is drawn in red. C) The elongation schedule used in the CTFcodon simulations. Regions marked in red in B take long elongation time. D) Three CTF schems. The CTFfast (dashed), the CTFslow (dotted), and the CTFcodon (solid) lines. E) A schematic view of the system including the wall-and-tunnel potential.

https://doi.org/10.1371/journal.pcbi.1004356.g001

These experimental data can be complemented with theoretical and computational analysis to deepen our understanding on the CTF mechanisms. Previously, lattice Monte Carlo simulations [23,24] and statistical theories [25,26] addressed physical aspects of CTF mechanisms. Coarse-grained molecular dynamics (CG MD) was used to investigate interaction with ribosome in the CTF [15,27–29]. These works helped understanding general and conceptual aspects of the CTF, but they were not specific enough to compare with experimental data of specific substrate proteins. It is time to start computational study of CTF for a specific protein, of which clear experimental data are available. This enables us to address structural aspects of CTF mechanisms, which is indeed the purpose of this work and we chose SufI for it.

Since the CTF becomes non-trivial primarily for relatively large and multi-domain proteins (SufI has three domains and is about 450 residue long (Fig 1A)), all-atom MD simulations are not feasible for this problem at the moment. By now, no all-atom MD simulation for folding to the native structure of multi-domain proteins was reported. To overcome size and time scale limit in all-atom MD simulations, protein folding simulations have commonly performed by coarse-grained (CG) models that are based on the energy landscape theory [30,31]. In particular, these simulations include medium-to-large proteins, such as multi-domain proteins[32–34].

Yet, to address mechanisms of the CTF and, in particular, an impact of elongation attenuation by CGMD simulations, technically, there are two major issues. First, we need to realize misfolding as well as correct folding in a well-balanced manner. Thus, the CG model needs to be calibrated so that an energy landscape is globally funneled in one hand and modestly rugged in the other hand. There have been a considerable number of studies towards hybrid modeling of structure-based potentials for globally funneled landscape with sequence-dependent terms for modestly rugged surfaces [35]. Yet, it should be noted that, currently, there is no established manner to balance the two aspects. Thus, here we unavoidably take a heuristic and empirical approach. Second, we need to design a scheme that mimics co-translational folding in silico. Quantitative kinetic measurements and detailed mechanisms of translation attenuation are not available at the moment, which led us to take a rather simplistic modeling of CTF scheme. Albeit these limitations, with the current CGMD, we can simulate complete folding and misfolding events of full-length SufI hundreds of times in scheme that mimics the CTF.

In this paper, we first describe computational modeling of CGMD for the CTF. Then, we performed the CTF and, as a control, the refolding simulations of SufI, comparing these results. Characteristics of misfolded structures are then analyzed. Next, folding networks for these simulations clarify impacts of CTF and the elongation attenuation on folding reaction mechanisms. Finally, the correlation between the degree of folding and the translation elongation time was investigated.

Results and Discussions

Computational modeling

In the current CG modeling, each amino acid is represented by one bead located at the Cα position. For folding simulations by CGMD, the so-called perfect-funnel model, or often called Go model, has been widely used giving many insightful lessons for folding dynamics [36–39]. However, the perfect-funnel approximation may not be sufficient to study CTF dynamics where successful folding competes with misfolding, or non-native traps. The latters are, by definition, not realized by the perfect funnel approximation. To this end, here we developed a hybrid CG model where we added a generic hydrophobic (HP) interaction potential VHP to the Go model potential VGo; the latter is responsible for globally funnel-like shape of the landscape, while the former makes the landscape modestly rugged leading to many metastable non-native traps. Concretely, the entire potential function of a protein is aVGo+bVHP. The Go potential was parameterized based on the atomic interaction at the native structure, called the AICG model developed by Li et al [32]. The HP interaction is a generic many-body potential that estimates how a hydrophobic residue is buried by other residues [40]. The HP interactions were applied not only natively interacting pairs, but also any residues. Detailed potential functions are described in Materials and Methods, Coarse-grained model.

As is well-known, proteins in vivo are gradually synthesized by ribosome from their N-termini and released from the ribosome exit tunnel, which we try to mimic in a simple manner. In the protocol, amino acids are added one by one to the C-terminus of the nascent polypeptide chain with certain “translation” rates (Fig 1B and 1C. See Materials and Methods, Coarse-grained model for details). To investigate effects of elongation attenuation, we employed the following three translation rate schemes (Fig 1D): 1) The uniformly fast translation scheme (a dashed line, designated as CTFfast), 2) the uniformly slow translation scheme (a dotted line, CTFslow), and 3) the non-uniform codon-based translation scheme (a solid line in Fig 1D and 1C, CTFcodon) that is dependent on the cognate tRNA concentration. We note that, in our scheme, the in silico translation rate is not proportional to the translation rate predicted from cognate tRNA concentrations. The translation attenuation was linked to a cluster of rare codons, which implies that the attenuation is a collective phenomenon and possesses distinct phases. Thus, using a threshold of the predicted translation rate, we introduced a two-phase approximation where the in silico translation is either "normal" or "slowed". The slowed translation phase, of which translation speed is 100 times slower than the normal case, corresponds to the translation attenuation. Since there is no quantitative kinetic measurement on the attenuation, this two-phase approximation and use of slowing factor 100 are rather simple, possibly over-simplified, schemes. Yet, we consider it qualitatively captures some of the major features of the translation attenuation. As far as the slowed phase is sufficiently slower that the normal phase, we expect qualitatively similar results. The relation between the translational time scales and the inherent folding time scales is of crucial importance, which will be discussed at the end of the results. Detailed in silico elongation scheme is described in Materials and Methods, Translational elongation scheme. Additionally, Vtunnel was introduced to mimic the ribosome steric effect that is realized by a combination of a wall and a tunnel (Fig 1E). Note that we did not include any molecular representation of ribosome and thus the tunnel is merely to restrict the nascent chain in a confined geometry. During elongation, a polypeptide chain is tethered to the base of the tunnel. On average, about 28 residues resided in the tunnel (S1 Fig). After completing the elongation, the chain is released from the base. We note that the exit tunnel was included to account for the gap between the residue at the catalytic center and the segment that can fold. The codon-based translation rate is based on the codon (sequence) at the catalytic center. In principle, some alpha helical structures can be formed in the exit tunnel depending of the sequence (although, retrospectively, we did not find it).

For comparison, we also performed folding simulations of SufI in a refolding scheme, where a full-length polypeptide chain started folding from denatured conformations obtained by high temperature simulations. No wall-and-tunnel potential Vtunnel was utilized in this scheme.

MD simulations were performed at 0.82TF*, where TF* is an upper limit of denaturation temperature in our CG model. To determine the temperature, starting from the native state of SufI, we performed unfolding simulations for 1 x 108 time steps at many temperatures. The lowest temperature at which we observed unfolding was defined as the upper limit of denaturation temperature TF*. (S2 Fig). We note that, even with the CG modeling, accurately calculating the denaturation temperature is a formidable task for this size of proteins; using the standard replica-exchange method or multi-canonical ensemble method, we did not succeeded to obtain the reversible folding/unfolding trajectories.

Co-translational folding and refolding simulations

First we compare a representative folding trajectory via the codon-based co-translational folding (CTFcodon) scheme with that via the refolding scheme. Fig 2 illustrates folding time courses quantified as the so-called Q-score defined as the fraction of formed contacts that exist at the native structure, together with some representative snapshots.

thumbnail
Download:
Fig 2. Representative time courses of folding simulations.

A) A time course of a refolding trajectory. B) that of the CTFcodon. Some snapshots were drawn with the same color code as Fig 1A.

https://doi.org/10.1371/journal.pcbi.1004356.g002

In the refolding trajectory shown in Fig 2A, the protein first acquired one globular region, which roughly corresponds to the N-domain. After a while, another globular region was formed, which contains, roughly, the C-domain and a half of M-domain. They gradually coalesced and made a single globular structure, which was a deep misfolded trap; the protein stayed in this trap until the end of the simulation.

On the other hand, the CTFcodon trajectory in Fig 2B showed markedly different time course. A cooperative folding of the N-domain at ~ 0.2 × 108 time step is followed by the folding of M-domain at ~ 1.2 × 108 time step. Subsequently, at ~ 1.7 × 108 time step, the protein folded to near native structure in which the C-domain is partly misfolded. Finally, at around 1.9 × 108 time step, it quickly transited into the native-like conformation.

More quantitatively, we repeated folding simulations of SufI 100 times both in the refolding and the CTFcodon schemes. In each trajectory, we judged whether the protein is folded or not by a set of native-ness scores, Q-scores, at the final 100 structures of the simulation (0 ≤ Q ≤ 1. Q = 1 at the native structure. We have both a generous and a stringent criteria for the judgment of folding. See Materials and Methods, Criteria for folding for more detail). Using a stringent criterion of folding, of 100 trajectories we found 18 successful folding cases in the refolding scheme (Table 1). Whereas, the CTFcodon resulted in 35 cases of correct folding. To clarify the statistical significance of the difference, we computed the histograms of Q-scores of the final structures in each scheme (S3 Fig). The difference in Q-score probability distributions was tested by the Kolmogorov-Smirnov test, which gave p-value of 0.000174 (Table 2, See also Table 3 for pairwise Mann-Whitney U tests). Thus, we conclude that the codon-based CTF simulation can fold SufI with significantly higher probability than the refolding can do.

thumbnail
Download:
Table 1. Number of successful folding casesa out of 100 trajectories.

https://doi.org/10.1371/journal.pcbi.1004356.t001

thumbnail
Download:
Table 2. Pairwise Kolmogorov-Smirnov tests to check the difference of the histograms of Qtotal-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t002

thumbnail
Download:
Table 3. Pairwise Mann-Whitney U tests to check the difference of the histograms of Qtotal-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t003

Effects of translational attenuations in co-translational folding

We then investigate effects of translational attenuation regions in SufI sequence, that was studied in experiments [16]. Experimentally, accelerating translation at certain slow translating regions, either by synonimous substitutions or by increasing concentrations of the rare tRNAs, inpaired SufI functions, most likely, due to misfolding. To test this idea in simulations, we conducted folding simulations by the CTF scheme, in which the chain is elongated with a uniform and fast rate across the entire chain (CTFfast).

In the same way as the CTFcodon case, we repeated the CTFfast simulations 100 times. Using the same criteria for the judgment of folding, i.e., Q-scores, we found only 20 cases of successful folding, which is much fewer than the CTFcodon scheme. The statistical analysis of the distribution suggested that the difference is significant (p = 0.00822). Actually, the result by the CTFfast scheme is statistically indistinguishable to that by the refolding scheme (p = 0.556). This is consistent with the experiment of Zhang et al [16].

Experimentally, lowering temperature could rescue the low-folding yield of the impaired folding scheme, which we now test in simulations. For the purpose, we performed folding simulations of SufI by the CTF where the elongation is slow and is in a uniform rate entirely (CTFslow). Of 100 simulations, we found 25 successful foldings by the same criteria as above. The statistics test resulted in no significance between th CTFslow and the CTFcodon schemes, while a subtle p value, p = 0.14 for the comparison between the slow and the fast CTF schemes.

To understand the CTF, comparison between the translation time scale and the folding time scale is of central importance. To estimate relevant folding time scales, for individual domains, we performed kinetic folding simulations. Time required to reach structures that have Q > 0.5 was computed for each of domains (S4 Fig). First, the M-domain is rather unstable and we could not observe successful folding of the standalone M-domain. The time scale for rough folding of N-domain τN−fold was 1.4 × 107 time steps, which is longer than that of the C-domain, τC−fold = 3.6 x 106 time steps. Interestingly, τN−fold is longer than the time to complete translation by the CTFfast scheme, τtranslation−fast ~ 4.4 × 106, but is comparable to that by the CTFslow scheme, τtranslation−slow ~ 1.3×107. Importantly, when the time for completion of the translation of N-domain is comparable to or longer than the average folding time of N-domain, the success ratio of SufI is high.

Misfolded structures

To understand why the codon-based CTF can facilitate folding of SufI, we now look into misfolded structures. For each of the four folding schemes, we analyzed probabilities of misfolding of individual domains at the ends of simulations (Fig 3A. Statistical test given in Tables 4–9). Here, the misfolded state was judged by the Q-scores of individual domains (To help understanding of typical Q-scores in SufI, we tabulated Q-scores of individual domains as well as those of interface for every snapshots in Fig 2 as S1 Table). Clearly, misfolding in the N-domain and the M-domain occurred with the highest probability by the refolding scheme, which is followed by the CTFfast scheme. The CTFcodon showed the smallest probabilities of misfolding for these domains. Of the four schemes, the rank order in misfolding of N- and M-domains is well (anti-)correlated with the probability of successful folding of the full-length SufI. (Table 1) In particular, probabilities of misfolding of the M-domain are markedly different between the refolding and the codon-based CTF. We note that the M-domain is not very stable and cannot fold as an isolated domain (S4 Fig). Folding of M-domain is achieved by structural support of the N-domain. In CTF schemes, when M-domain is synthesized and released from the exit tunnel, the N-domain has large chance to be folded. The folding of the C-domain is not much different among the four schemes.

thumbnail
Download:
Fig 3. Misfolding in SufI domains.

A) Fractions of misfolded domains at the end of simulations in four different schemes; the refolding (black), the CTFfast (red), the CTFslow (green), and the CTFcodon (blue). B) Representative final structures of misfolding. i) structure that is misfolded in N-domain. ii) Misfolded in M-domain. iii) (right) Misfolded in C-domain. (left) Native structure for comparison. See text for the explanation of the block arrows.

https://doi.org/10.1371/journal.pcbi.1004356.g003

thumbnail
Download:
Table 4. Pairwise Kolmogorov-Smirnov tests to check the difference of the histograms of N-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t004

thumbnail
Download:
Table 5. Pairwise Kolmogorov-Smirnov tests to check the difference of the histograms of M-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t005

thumbnail
Download:
Table 6. Pairwise Kolmogorov-Smirnov tests to check the difference of the histograms of C-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t006

thumbnail
Download:
Table 7. Pairwise Mann-Whitney U tests tests to check the difference of the histograms of N-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t007

thumbnail
Download:
Table 8. Pairwise Mann-Whitney U tests to check the difference of the histograms of M-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t008

thumbnail
Download:
Table 9. Pairwise Mann-Whitney U tests to check the difference of the histograms of C-domain’s Q-scores of the final structures.

https://doi.org/10.1371/journal.pcbi.1004356.t009

We now show some representative misfolded structures (Fig 3B). A conformation in Fig 3B (i) taken from a refolding trajectory, is misfolded in the N-domain, while the M- and C-domains are well-folded (The non-native contact map is given in S6 Fig (i)). In this structure, C-terminal end of the N-domain is unfolded and is flipped out to the left side in the figure (See the block arrow. Also, see S6 Fig (i) for many non-native contacts in C-terminus of the N-domain). With this flipped-out segment, three domains coalesced to form near-native domain-domain interfaces. Once the interfaces are firmly formed, the protein is topologically trapped and an escape event from this trap was not realized. Fig 3B(ii) illustrates a case where C-terminus segment of the M-domain was entangled with the C-domain (the block arrow. See also the non-native contact map S6 Fig (ii)). Again, the domain-domain interfaces are near-native like, which makes an escape from this trap difficult. The right cartoon of Fig 3B(iii) shows the case that a N-terminal segment of the C-domain, 314–340 residues (shown in red-and-gray striped pattern with block arrow) goes through different paths from the native structure (the left cartoon of Fig 3B (iii)).

Folding network

Next, we investigated the ensemble of folding pathways for the CTF and the refolding schemes. To clarify folding pathways, we drew folding networks where nodes represent discretized conformational states and links represent transitions between the states[41,42].

Conformational states were discretized by the native-ness scores (Q-scores) and by the non-native contact scores (N-scores) (See Materials and Methods, Discretization of states by Q-scores of parts for more details). For each domain and each interface between domains, we defined Q-score and N-score (we have six Q- and N-scores, in total). As usual, the Q-score measures fraction of formed native contacts. The N-score is defined as the number of non-native contacts normalized by its maximal number observed. Each Q-score is categorized into 5 classes, while each N-score is divided into 3 classes. Together, we have as many as 56 × 36 ~ 1.1 × 107 states (nodes). To simplify the network, we removed any loops that go from a node and return to the same node later. All 100 trajectories were used to draw a network for each folding scheme.

We depict folding networks of SufI for four different folding schemes (Figs 4 and S5). Comparing the folding networks of the refolding (Fig 4A) and the CTFcodon (Fig 4B) schemes, we found, first of all, that the network for the refolding has much more nodes (3284 nodes) than the CTFcodon has (820 nodes). By refolding, the protein exhibited much more divergent conformational states, many of which are characterized by low Q-scores and high N-scores. Second, while the refolding scheme did not show any dominant pathways, the CTFcodon has a clear folding route from the top in the figure to the bottom. Obviously, the CTF enforced SufI to fold vectorially from N-terminal, which provided constraints to the order of domain folding events. In contrast, the refolding scheme made a protein fold freely from any segments resulting into diverse transitions. The CTF restricts kinetics of proteins and reduces conformational ensemble being observed, and are consistent with earlier theoretical works[23].

thumbnail
Download:
Fig 4. SufI folding networks for refolding (A) and for the codon-based CTF (B).

The refolding network possesses 3284 nodes, while the codon-based CTF has only 820 nodes. The size of nodes represent their probabilities. The darkness of the node represents native-ness. The darker one is closer to the native. Diamonds, triangle, and stars indicate that N-, M-, and C-domains are pre-dominantly unfolded, respectively. When pre-dominantly unfolded domains are no uniquely decided, circles are used.

https://doi.org/10.1371/journal.pcbi.1004356.g004

The folding network for the CTFfast scheme (S5A Fig) apparently looks similar to that of the refolding. The number of nodes found was 3108, which is only slightly fewer than that of the refolding network, i.e. 3284. The transitions are diverse with no dominant pathway to the native state.

On the other hand, the slow CTF scheme showed the folding network (S5B Fig) rather similar to that in the CTFcodon scheme. The number of nodes found in the slow CTF was 1096, which is slightly larger than that found in the CTFcodon, i.e. 820. We see a single and nearly identical folding route in these two schemes.

Correlation between translation rate and folding

It is interesting to ask to what extent the translation rate is designed (optimized), via codon usage, to facilitate folding. To this end, here we investigate the correlation, if any, between a putative translation rate and the degree of folding. For the former, we simply use the translation rate, in arbitrary unit, predicted by an algorithm proposed in Spencer et al (Fig 5B) [43]. This translation rate is encoded in the codon usage as well as tRNA concentrations and other factors, but not apparently dependent on the physical chemistry of folding. For the degree of folding acquisition, we defined the progress of native-ness ΔQi in a nascent chain of the length i as where 〈Q〉L=i is the average Q-score when the nascent chain has the length i and 〈 〉100 trajectories means the average over 100 trajectories in the slow CTF scheme. If ΔQi is high at i-th residue, a nascent chain gains Q-score without disturbance from more C-terminal region of the chain. Note that if we used the codon-based CTF scheme in calculation of ΔQi, it would naturally correlate with the translation rate. Importantly, however, we did not bias the CTF by the codon usage. Instead, we used a uniform and slow CTF scheme. Thus, ΔQi is not directly related to the difference in the translation rate, but is a purely physicochemical quantity determined by the amino acid sequence. We note that ΔQi was smoothed by a window average of the 5-residue windows to reduce the noise.

thumbnail
Download:
Fig 5. The correlation between the inverse of elongation rate and the degree of folding in SufI.

A) The degree of folding acquisition ΔQi after averaging over the window size 5. B) One over the translation rate computed from the Spencer et al.’s algorithm [43]. Experimentally-detected translational attenuation regions, 33-40kDa (281-326th residues) and 25-28kDa (214-240th residues), are shaded in grey [16]. C) The scattered plot of the translation time and the degree of folding. Here, residues 200–350 are used. The correlation coefficient was 0.51.

https://doi.org/10.1371/journal.pcbi.1004356.g005

The ΔQi profile shown in Fig 5A exhibits several peaks. First, we focus on the peaks that correspond to folding of M-domain because it is the most difficult event. We find a high ΔQi region around 280–310, which well correlates with a translational attenuation region, 33-40kDa region (281–326 residues, grey shaded in Fig 5). Experimentally, synonymous substitutions of rare codons in this region reduced resistance to a protease [16]. The other translational attenuation experimentally tested is 25-28kDa (214–240 residues), in which synonymous substitution of two leucine codons impaired the protease resistance of SufI. In Fig 5A, we see peak in the ΔQi profile at ~245. More quantitatively, by using 200th-350th residues, we computed the correlation between the ΔQi profile and the translation rate profile (Fig 5C) finding the correlation coefficient 0.51. Thus, they are indeed, albeit modestly, correlated.

The highest peak of the ΔQi profile in Fig 5A is located at 166-th residue, which corresponds to the situation that the N-domain (1–143 residues) is mostly released from the ribosome exit tunnel. (Remember that the average number of residues in the exit tunnel is 28 as in S1 Fig). However, the translation profile in Fig 5B does not indicate any attenuation in this region. It seems that misfolding in the N-domain is not very probable in any CTFs and thus translational attenuation at this point is not required for successful folding.

Conclusions

Comprehensively performing molecular simulations of co-translational folding (CTF) and refolding of SufI, we elucidated mechanisms of how translational attenuation can facilitate correct folding from structural perspectives. First, coarse-grained simulations showed that the codon-based CTF, CTFcodon, exhibited higher probability of correct folding than the refolding did. When the translational attenuation is removed, the CTFfast simulations resulted in the success rate similar to that by the refolding scheme. When the elongation was uniformly slowed down, the CTFslow simulation gave essentially the same results as those of CTFcodon. These are all consistent with recent experiments. On top, the simulations provided much of structural and mechanistic insights. Specifically for SufI, we found that the M-domain is least stable and can fold only when it is supported by the pre-folded N-domain. Once a segment of the M-domain is entangled with either N- or C-domain, an escape from the trap was difficult. Combining molecular simulations with biochemical experiments provided detailed mechanistic understanding of CTFs.

A recent theoretical study suggested that, under certain situations, fast translation can coordinate folding to the native structure [44]. Apparently, this is not the case in our SufI simulations. Whether slower or faster translation facilitates the correct folding depends on the folding kinetic network as was shown in [44]. We need some more investigations for specific proteins, through which we know which scenarios are more common.

We note that the current CG modeling has some limitations. One of the major limitations is on the time scales. Using the CG modeling, one cannot easily estimate the absolute time scales of folding and translation. Using a low viscosity in Langevin dynamics and structure-based potentials, we speeded up the folding kinetics some orders of magnitude. Translation kinetic parameters in the normal and slowed phases are not accurately known. This makes quantitative comparison difficult. Another limitation is the balance between the structure-based potential and the sequence-dependent terms, which was determined empirically here. Accurate modeling of these balances is highly desired in future work.

Materials and Methods

Model protein

In this study, we studied folding of a three-domain protein SufI [45] (Fig 1, PDB code: 2UXT). Starting with the PDB structure 2UXT, we removed the His-tag and modeled missing residues by MODELLER [46], resulting in the 443-residue long protein model. The model structured was refined by the energy minimization with AMBER [47]. Using Pfam’s [48], we defined three domains; N-terminal domain as 1–143, the middle (M-) domain as 160–300, and the C-terminal domain as region 314–443. Segments between two domains are termed linkers. The linker between M- and C- domains are rather long and extended.

Coarse-grained model

In the simulation, one residue is represented by one CG particle which locates at Cα position. We used our in-house developing software CafeMol for all the simulations [49].

The potential energy function consists of the native-based AICG2+ potential (VGo) and non-local many body hydrophobic interaction potential (VHP). The total energy Vtotal for the refolding simulation is given as where a and b are coefficients to control the balance between two terms. The potential for the CTF simulations is written as

The native-based potential VGo is defined as [32]:

The first term keeps virtual bonds between consecutive amino acids, the second and the third terms represent statistical potential for virtual bond-angles and virtual dihedral-angles [50]. The fourth and the fifth terms define native-based local interactions [32]. The sixth term is non-local contact interaction for natively contacting pairs. The last term is a generic excluded volume interaction (See [32] for more details).

For the hydrophobic interaction, we take the function developed in[40], which is written in the form: where is a parameter that reflects the hydrophobicity of amino acids for the amino acid type A(i). SHP represents the buried-ness of the amino acid i and is defined as: where clinear and ρmin are constants and ρi represents local density and is calculated by: where nA(i) is the number of heavy atoms that defines the amino acid A(i) represents and nmax,A(i) is the maximum coordination number for particle type A(i). The function uHP represents the degree of the contact between particle i and particle j and is defined as below a sigmoidal function:

We note that the described hydrophobic interaction potential was first developed for a CG model that uses different resolution from the current work. Thus we need to re-parameterize the function. We estimated parameters rmin,A(i),A(j) rmax,A(i),A(j), , nmax,A(i), and nA(i) for each amino acid types in the following way. Using Dunbrack’s culled PDB set [51], we analyzed radius distributions of twenty types of amino acids. For details, if a distance between heavy atoms of two amino acids is less than RvdW,i + RvdW,j + RvdW,H2O, where RvdW,i is the van der Walls radius of the atom i, we defined the distance between Cα's as an effective distance, obtaining a set of radial distribution of 20x20 amino acid combinations. Then, we defined 95% confidence coefficient of their histograms as rmax,A(i),A(j) and we set rmin,A(i),A(j) = rmax,A(i),A(j) −4. is taken from hydrophobic indices of Fauchere & Pliska [52]. All these parameters are included in the latest CafeMol and publicly available.

In the total potential energy Vtotal, VGo(R) is responsible to make a globally funnel like energy landscape, while VHP(R) makes the landscape modestly rugged via physicochemical interactions. Thus, the balance of the two potentials is of central importance in the simulation. Since it is not straightforward to decide the coefficients in ab initio manner, in this work, instead, we took an empirical approach. At the native structure of SufI, we first calculated the potential energy VGo(Rnat). We assumed that this value is a reasonable energy at the native structure and fixed this value at the native structure. We then express it as a linear combination of VGo(R) and VHP(R). Formally, it can be written as where a free parameter a was decided fully empirically. With several values of a, we performed some preliminary simulations by the refolding scheme, estimating the probability of the successful folding, of which the criterion is defined below. We ended up with a = 0.8, by which about 20% of runs could reach the native structure. The coefficient b = 2.13 was derived from this procedure.

To reproduce a steric effect of ribosome exit tunnel and surface, we added a pure repulsive wall-and-tunnel potential Vtunnel defined as: where di is the distance between the particle i and the wall-and-tunnel. The default parameters in CafeMol were used for εex and C. The radius and length of the tunnel were set as rT = 15 Å, lT = 90 Å, respectively.

We note that all the interaction potentials here are temperature independent. Since hydrophobic interactions are effective interaction that itself depends on temperature, one can include temperature dependence as in Chan et al for more accurate modeling [53].

Langevin dynamics

Molecular dynamics was simulated by the Langevin equation at the constant temperature T, where γi is a friction constant and ξi is a random force. This random force satisfies 〈ξi(t)〉 = 0, and 〈ξi(t)ξj(t')〉 = 2γ(kT / m)δ(t – t')δi,j. The stationary distribution generated by this Langevin equation is the Boltzmann distribution for a given temperature T. The force Fi is derived from partial differentiation of the potential energy function. For numerical integration, we used the scheme in [54,55]. In the simulation, γ is 0.02, and a finite time step Δt = 0.1.

Translation elongation scheme

In simulations that mimic CTF, we increased the chain length of the nascent polypeptide one by one residue and used a wall-and-tunnel potential that represents the rough geometry of the ribosome exit tunnel (Fig 1E). The C-terminal residue of the nascent chain was fixed to the base of the ribosome tunnel during the elongation and is released when the final residue was “synthesized”. We assumed that the time scale for the covalent bond formation (the synthesis) is much shorter than time scale to wait the cognate tRNA and that for folding, and thus the synthesis is treated as the instantaneous change in the chain length. We also ignored any mechanistic factors possibly involved in the synthesis. Simply, we shifted the nascent chain toward the exit direction and added one residue at the base of the tunnel at the one-step elongation.

In a scheme that mimics the CTF rate that depends on codon (CTFcodon), we used the translation elongation profile derived from Spencer’s algorithm [43]. Spencer’s algorithm generates relative translation profile for each organism (Fig 1B), by distinguishing Watson-Crick interactions from non-Watson-Crick (wobble) interactions. We note that we took this relatively simple algorithm, although there can be other algorithms. The number of tRNA genes for every codon was referred from Genomic tRNA database (gtrnadb.ucsc.edu) [56]. The mechanistic detail of the translational attenuation is unknown at the moment, and so, when an elongation rate is under a threshold, we defined the codon as a rare codon and the elongation was attenuated for 106steps per one residue. For other residues, we used the elongation rate as 104steps per one residue. The scheme was termed the CTFcodon (See Fig 1C).

To test the effect of synonymous substitution that remove the translational attenuation, we set the CTF in which the elongation rate is fast and uniform. A protein is elongated at the rate of 104steps per one residue. This is termed as the CTFfast.

To test the effect of lowering temperature, we set the CTF in which the elongation rate is uniform and is slow (CTFslow). The elongation speed is 3×105 steps per one residue.

As a control, we also set up the refolding scheme. In this scheme, a wall-and-tunnel potential was not used and the full-length SufI was present from the beginning. The initial unfolded conformation was prepared by constant temperature simulation at a high temperature for 107 time steps from the native state. This was sufficient to prepare a fully unfolded structure.

In all four schemes, we ran 100 trajectories, and each trajectory is simulated for 3×108 time steps, including the time for translation in the cases of CTF schemes. The comparison of three elongation schemes is given in Fig 1D.

Criteria for folding

To judge whether SufI is folded or not, we introduced multiple native-ness scores, i.e., Q-scores.

In general, the widely used Q-score is defined as the number of formed native contacts relative to that presents in the native structures. First, an amino acid pair ij is defined as the native contact when one atom, except hydrogen, in i-th residue is within 6.5Å from at least one atom in j-th residue in the native structure. For natively contacting pairs, we check the Cα-Cα distance in a snapshot of folding simulations. If it is within 1.2 times the corresponding distance at the native structure, we assign the contact is formed in the snapshot.

The Q-score can be defined either for the full-length protein Qtotal, or for any part of the protein, both of which were used in this work. Qtotal is convenient to quantify an overall native-ness by one value. When Qtotal is above 0.95, SufI takes native state with high probability (this is called as a generous criterion for the native state). During the analysis, however, we noticed that, for multi-domain proteins such as SufI, the completion of folding cannot easily be assessed by Qtotal alone. For example, we found that individual domains are all correctly folded, while some domain-domain interfaces are not. Since the number of native contacts for the domain interface is much less than those within the domains, these structures often take Qtotal values close to one. (Even worse is that these values can be within the thermal fluctuation range of the Qtotal at the true native state.) To distinguish these misfolded structures, we need to check Q-scores for every domain-domain interfaces, separately. Specifically for SufI, we introduced Q-scores for individual domains (three in total) as well as Q-scores for domain-domain interfaces (three in total), leading to six Q-scores of parts. When all the six Q-scores of parts are above their thresholds, we stringently assigned the structure well-folded (a stringent criterion for the native state).

Discretization of states by Q-scores of parts

Q-scores for N-, M-, and C-domains and for N-M, N-C, and M-C domain-domain interfaces are classified by four thresholds. We located those thresholds at the local minima of statistical weight distributions. Specifically, thresholds of N-domain’s Q-score is [0.50, 0.63, 0.88, 0.94], those of M-domain: [0.30, 0.64, 0.78, 0.90], C-domain: [0.39, 0.59, 0.90, 0.95], N-M interface: [0.10, 0.30, 0.64, 0.83], N-C interface: [0.19, 0.44, 0.76, 0.90], M-C interface: [0.31, 0.56, 0.81, 0.90].

Drawing folding network

To draw a folding network, we used a physical model of network, which is called a spring-electrical model [57]. In this model, each node is represented as a mass point and possesses a positive charge. If two nodes linked each other, the pair of nodes has an elastic energy. We seek locations of nodes that minimize the total “energy” function. We obtained an optimized network structure by a simulated annealing.

To discretize structural conformations, we classified six Q-scores of parts and six N-scores of parts. Here, N-score represents degree of formed non-native contacts and was defined as the number of formed non-native contacts relatively to the maximal number of the same contacts. Based on the thresholds, we can assign conformations to one of 56×36 nodes and represent a trajectory by a polygonal line that transits from a node to another. For simplicity, we removed any loops. Here, a loop is a sequence of transitions that start from and return to one node.

Supporting Information

S1 Table. Q-scores of intra-domains and inter-domains for snapshots depicted in Fig 2.

https://doi.org/10.1371/journal.pcbi.1004356.s001

(DOCX)

S1 Fig. Observed numbers of residues that resided in the ribosome tunnel upon the elongation.

The average is about 28 residues.

https://doi.org/10.1371/journal.pcbi.1004356.s002

(TIF)

S2 Fig. Temperature dependence of SufI unfolding.

Starting from a denatured state, we performed folding simulations for 108 time steps. Temperatur is given in CafeMol unit. The sudden drop in average Q-score was found at the temperatur 440, which corresponds to TF*. Folding simulations were conducted at 360, which corresponds to 0.82 TF.

https://doi.org/10.1371/journal.pcbi.1004356.s003

(TIF)

S3 Fig. Histograms of Qtotal-score in the final conformations.

In each folding scheme, the last 100 snapshots (corresponding to 105 time steps) are used.

https://doi.org/10.1371/journal.pcbi.1004356.s004

(TIF)

S4 Fig. Folding time course of standalone domains of SufI in normal (A) and in logarithmic (B) scales.

(C) The linear fitting is used to obtain folding times of individual domain. Blue, green, and red curves correspond to folding of N-, M-, and C-domains.

https://doi.org/10.1371/journal.pcbi.1004356.s005

(TIF)

S5 Fig. Protein folding networks drawn from CTFfast and CTFslow folding schemes.

The meaning of symbols are identical to those in Fig 4.

https://doi.org/10.1371/journal.pcbi.1004356.s006

(TIF)

S6 Fig. Non-native contact maps of the representative misfolded structures.

The upper right triangle part shows the probability map of non-native map formed in the last 100 snapshots (corresponding to 105 time steps) in representative trajectories. The lower triangl part shows the native contact map obtained from the native structure. The (i),(ii) and (iii) are three representative misfolded structures corresponding to the same symbols in Fig 3B.

https://doi.org/10.1371/journal.pcbi.1004356.s007

(TIF)

Author Contributions

Conceived and designed the experiments: TT NH ST. Performed the experiments: TT. Analyzed the data: TT NH. Contributed reagents/materials/analysis tools: TT. Wrote the paper: TT NH ST.

References

  1. 1. Onuchic JN, Wolynes PG. Theory of protein folding. Curr Opin Struct Biol. Elsevier; 2004;14: 70–75. pmid:15102452
  2. 2. Kubelka J, Hofrichter J, Eaton W a. The protein folding “speed limit”. Curr Opin Struct Biol. 2004;14: 76–88. pmid:15102453
  3. 3. Han J-H, Batey S, Nickson A a, Teichmann S a, Clarke J. The folding and evolution of multidomain proteins. Nat Rev Mol Cell Biol. 2007;8: 319–330. pmid:17356578
  4. 4. Fitter J. The perspectives of studying multi-domain protein folding. Cell Mol Life Sci. 2009;66: 1672–1681. pmid:19183848
  5. 5. Gershenson A, Gierasch LM. Protein folding in the cell: challenges and progress. Curr Opin Struct Biol. Elsevier Ltd; 2011;21: 32–41. pmid:21112769
  6. 6. Sarkar M, Smith AE, Pielak GJ. Impact of reconstituted cytosol on protein stability. Proc Natl Acad Sci U S A. 2013;110: 19342–19347. pmid:24218610
  7. 7. Ebbinghaus S, Dhar A, McDonald JD, Gruebele M. Protein folding stability and dynamics imaged in a living cell. Nat Methods. Nature Publishing Group; 2010;7: 319–323. pmid:20190760
  8. 8. Hartl FU, Bracher A, Hayer-Hartl M. Molecular chaperones in protein folding and proteostasis. Nature. Nature Publishing Group; 2011;475: 324–332. pmid:21776078
  9. 9. Zhang G, Ignatova Z. Folding at the birth of the nascent chain: coordinating translation with co-translational folding. Curr Opin Struct Biol. Elsevier Ltd; 2011;21: 25–31. pmid:21111607
  10. 10. Zipser D, Perrin D. Complementation on ribosomes. Cold Spring Harbor Symposia on Quantitative Biology. 1963. pp. 533–537.
    • 11. Fedorov AN, Baldwin TO. Cotranslational Protein Folding. J Biol Chem. 1997;272: 32715–32718. pmid:9407040
    • 12. Gloge F, Becker AH, Kramer G, Bukau B. Co-translational mechanisms of protein maturation. Curr Opin Struct Biol. 2014;24: 24–33. pmid:24721450
    • 13. Frydman J, Erdjument-Bromage H, Tempst P, Hartl FU. Co-translational domain folding as the structural basis for the rapid de novo folding of firefly luciferase. Nat Struct Biol. 1999;6: 697–705. pmid:10404229
    • 14. Ugrinov KG, Clark PL. Cotranslational folding increases GFP folding yield. Biophys J. Biophysical Society; 2010;98: 1312–1320.
    • 15. O’Brien EP, Christodoulou J, Vendruscolo M, Dobson CM. New scenarios of protein folding can occur on the ribosome. J Am Chem Soc. 2011;133: 513–526. Available: http://pubs.acs.org/doi/abs/10.1021/ja107863z pmid:21204555
    • 16. Zhang G, Hubalewska M, Ignatova Z. Transient ribosomal attenuation coordinates protein synthesis and co-translational folding. Nat Struct Mol Biol. 2009;16: 274–280. pmid:19198590
    • 17. Komar AA. A pause for thought along the co-translational folding pathway. Trends Biochem Sci. 2009;34: 16–24. pmid:18996013
    • 18. O’Brien EP, Ciryam P, Vendruscolo M, Dobson CM. Understanding the Influence of Codon Translation Rates on Cotranslational Protein Folding. Acc Chem Res. 2014;47: 1536–1544. pmid:24784899
    • 19. Spencer P, Barral J. Genetic code redundancy and its influence on the encoded polypeptides. Comput Struct Biotechnol J. 2012;1. Available: http://journals.sfu.ca/rncsb/index.php/csbj/article/view/13
    • 20. Gloge F, Becker AH, Kramer G, Bukau B. Co-translational mechanisms of protein maturation. Curr Opin Struct Biol. 2014;24: 24–33. pmid:24721450
    • 21. Purvis IJ, Bettany a J, Santiago TC, Coggins JR, Duncan K, Eason R, et al. The efficiency of folding of some proteins is increased by controlled rates of translation in vivo. A hypothesis. J Mol Biol. 1987;193: 413–417. Available: http://www.ncbi.nlm.nih.gov/pubmed/3298659 pmid:3298659
    • 22. Komar AA, Lesnik T, Reiss C. Synonymous codon substitutions affect ribosome traffic and protein folding during in vitro translation. FEBS Lett. 1999;462: 387–391. Available: http://www.ncbi.nlm.nih.gov/pubmed/10622731 pmid:10622731
    • 23. Morrissey MP, Ahmed Z, Shakhnovich EI. The role of cotranslation in protein folding: a lattice model study. Polymer. Elsevier; 2004;45: 557–571.
    • 24. Krobath H, Shakhnovich EI, Faísca PFN. Structural and energetic determinants of co-translational folding. J Chem Phys. 2013;138: 215101. pmid:23758397
    • 25. Deane CM, Dong M, Huard FPE, Lance BK, Wood GR. Cotranslational protein folding—fact or fiction? Bioinformatics. 2007;23: i142–148. pmid:17646290
    • 26. Saunders R, Deane CM. Synonymous codon usage influences the local protein structure observed. Nucleic Acids Res. 2010;38: 6719–6728. pmid:20530529
    • 27. Elcock AH. Molecular Simulations of Cotranslational Protein Folding: Fragment Stabilities, Folding Cooperativity, and Trapping in the Ribosome. PLoS Comput Biol. 2006;2: e98. pmid:16789821
    • 28. O’Brien EP, Hsu S-TD, Christodoulou J, Vendruscolo M, Dobson CM. Transient tertiary structure formation within the ribosome exit port. J Am Chem Soc. 2010;132: 16928–16937. pmid:21062068
    • 29. O’Brien EP, Vendruscolo M, Dobson CM. Prediction of variable translation rate effects on cotranslational protein folding. Nat Commun. Nature Publishing Group; 2012;3: 868. pmid:22643895
    • 30. Bryngelson JD, Wolynes PG. Spin glasses and the statistical mechanics of protein folding. Proc Natl Acad Sci U S A. 1987;84: 7524–7528. pmid:3478708
    • 31. Onuchic JN, Luthey-Schulten Z, Wolynes PG. Theory of protein folding: the energy landscape perspective. Annu Rev Phys Chem. 1997;48: 545–600. pmid:9348663
    • 32. Li W, Terakawa T, Wang W, Takada S. Energy landscape and multiroute folding of topologically complex proteins adenylate kinase and 2ouf-knot. Proc Natl Acad Sci U S A. 2012;109: 17789–17794. pmid:22753508
    • 33. Wang Y, Chu X, Suo Z, Wang E, Wang J. Multidomain Protein Solves the Folding Problem by Multifunnel Combined Landscape: Theoretical Investigation of a Y-Family DNA Polymerase. J Am Chem Soc. 2012; https://doi.org/10.1021/ja3045663
      • 34. Ito M, Ozawa T, Takada S. Folding Coupled with Assembly in Split Green Fluorescent Proteins Studied by Structure-Based Molecular Simulations. J Phys Chem B. 2013; https://doi.org/10.1021/jp4032817
        • 35. Chen T, Song J, Chan HS. Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding. Curr Opin Struct Biol. Elsevier; 2015;30: 32–42. pmid:25544254
        • 36. Taketomi H, Ueda Y, Go N. Studies on protein folding, unfolding and fluctuations by computer simulation. Int J Pept Protein Res. Wiley Online Library; 1975;7: 445–459. pmid:1201909
        • 37. Clementi C, Nymeyer H, Onuchic JN. Topological and energetic factors: what determines the structural details of the transition state ensemble and “en-route” intermediates for protein folding? An investigation for small globular proteins. J Mol Biol. 2000;298: 937–953. pmid:10801360
        • 38. Koga N, Takada S. Roles of native topology and chain-length scaling in protein folding: a simulation study with a Go-like model. J Mol Biol. 2001;313: 171–180. pmid:11601854
        • 39. Hills RD, Brooks CL. Insights from coarse-grained gō models for protein folding and dynamics. Int J Mol Sci. 2009;10: 889–905. pmid:19399227
        • 40. Fujitsuka Y, Takada S, Luthey-Schulten ZA, Wolynes PG. Optimizing physical energy functions for protein folding. Proteins. 2004;54: 88–103. pmid:14705026
        • 41. Hori N, Chikenji G, Berry RS, Takada S. Folding energy landscape and network dynamics of small globular proteins. Proc Natl Acad Sci U S A. 2009;106: 73–78. pmid:19114654
        • 42. Rao F, Caflisch A. The protein folding network. J Mol Biol. 2004;342: 299–306. pmid:15313625
        • 43. Spencer PS, Siller E, Anderson JF, Barral JM. Silent substitutions predictably alter translation elongation rates and protein folding efficiencies. J Mol Biol. Elsevier Ltd; 2012;422: 328–335. pmid:22705285
        • 44. O’Brien EP, Vendruscolo M, Dobson CM. Kinetic modelling indicates that fast-translating codons can coordinate cotranslational protein folding by avoiding misfolded intermediates. Nat Commun. Nature Publishing Group; 2014;5.
        • 45. Tarry M, Arends SJR, Roversi P, Piette E, Sargent F, Berks BC, et al. The Escherichia coli cell division protein and model Tat substrate SufI (FtsP) localizes to the septal ring and has a multicopper oxidase-like structure. J Mol Biol. Elsevier Ltd; 2009;386: 504–19. pmid:19135451
        • 46. Fiser A, Do RKG, Sali A. Modeling of loops in protein structures. Protein Sci. 2000;9: 1753–1773. pmid:11045621
        • 47. Case DA, Darden TA, T.E. Cheatham I, Simmerling CL, Wang J, Duke RE, et al. AMBER 11. University of California, San Francisco; 2010.
          • 48. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40: D290–301. pmid:22127870
          • 49. Kenzaki H, Koga N, Hori N, Kanada R, Li W, Okazaki K, et al. CafeMol: A Coarse-Grained Biomolecular Simulator for Simulating Proteins at Work. J Chem Theory Comput. 2011;7: 1979–1989.
          • 50. Terakawa T, Takada S. Multiscale ensemble modeling of intrinsically disordered proteins: p53 N-terminal domain. Biophys J. Biophysical Society; 2011;101: 1450–1458.
          • 51. Wang G, Dunbrack RL. PISCES: a protein sequence culling server. Bioinformatics. 2003;19: 1589–1591. pmid:12912846
          • 52. Fauchere JL, Pliska V. Hydrophobic parameters-pi of amino-acid side-chains from the partitioning of N-acetyl-amino-acid amides. Eur J Med Chem. 1983;18: 369–375.
          • 53. Chan HS, Zhang Z, Wallin S, Liu Z. Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models. Annu Rev Phys Chem. 2011;62: 301–26. pmid:21453060
          • 54. Honeycutt JD, Thirumalai D. The nature of folded states of globular proteins. Biopolymers. 1992;32: 695–709. pmid:1643270
          • 55. Guo Z, Thirumalai D. Kinetics of protein folding: nucleation mechanism, time scales, and pathways. Biopolymers. 1995;36: 83–102.
          • 56. Chan PP, Lowe TM. GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res. 2009;37: D93–97. pmid:18984615
          • 57. Hu Y. Efficient, High-Quality Force-Directed Graph Drawing. Math J. 2005;10: 37–71.