Protein translocation through a tunnel induces changes in folding kinetics: a lattice model study

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.

Dynamics of Co-translational Membrane Protein Integration and Translocation via the Sec Translocon

An important aspect of cellular function is the correct targeting and delivery of newly synthesized proteins. Central to this task is the machinery of the Sec translocon, a transmembrane channel that is involved in both the translocation of nascent proteins across cell membranes and the integration of proteins into the membrane. Considerable experimental and computational effort has focused on the Sec translocon and its role in nascent protein biosynthesis, including the correct folding and expression of integral membrane proteins. However, the use of molecular simulation methods to explore Sec-facilitated protein biosynthesis is hindered by the large system sizes and long (i.e., minute) time scales involved. In this work, we describe the development and application of a coarse-grained simulation approach that addresses these challenges and allows for direct comparison with both in vivo and in vitro experiments. The method reproduces a wide range of experimental observations, providing new insights into the underlying molecular mechanisms, predictions for new experiments, and a strategy for the rational enhancement of membrane protein expression levels.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

Insights into Cotranslational Nascent Protein Behavior from Computer Simulations

Regulation of protein stability and function in vivo begins during protein synthesis, when the ribosome translates a messenger RNA into a nascent polypeptide. Cotranslational processes involving a nascent protein include folding, binding to other macromolecules, enzymatic modification, and secretion through membranes. Experiments have shown that the rate at which the ribosome adds amino acids to the elongating nascent chain influences the efficiency of these processes, with alterations to these rates possibly contributing to diseases, including some types of cancer. In this review, we discuss recent insights into cotranslational processes gained from molecular simulations, how different computational approaches have been combined to understand cotranslational processes at multiple scales, and the new scenarios illuminated by these simulations. We conclude by suggesting interesting questions that computational approaches in this research area can address over the next few years.

Functional implications of ribosomal RNA methylation in response to environmental stress

The study of post-transcriptional RNA modifications has long been focused on the roles these chemical modifications play in maintaining ribosomal function. The field of ribosomal RNA modification has reached a milestone in recent years with the confirmation of the final unknown ribosomal RNA methyltransferase in Escherichia coli in 2012. Furthermore, the last 10 years have brought numerous discoveries in non-coding RNAs and the roles that post-transcriptional modification play in their functions. These observations indicate the need for a revitalization of this field of research to understand the role modifications play in maintaining cellular health in a dynamic environment. With the advent of high-throughput sequencing technologies, the time is ripe for leaps and bounds forward. This review discusses ribosomal RNA methyltransferases and their role in responding to external stress in Escherichia coli, with a specific focus on knockout studies and on analysis of transcriptome data with respect to rRNA methyltransferases.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

Measuring cotranslational folding of nascent polypeptide chains on ribosomes

Protein folding has been studied extensively in vitro, but much less is known about how folding proceeds in vivo. A particular distinction of folding in vivo is that folding begins while the nascent polypeptide chain is still undergoing synthesis by the ribosome. Studies of cotranslational protein folding are inherently much more complex than classical in vitro protein folding studies, and historically there have been few methods available to produce the quantities of pure material required for biophysical studies of the nascent chain, or assays to specifically interrogate its conformation. However, the past few years have produced dramatic methodological advances, which now place cotranslational folding studies within reach of more biochemists, enabling a detailed comparison of the earliest stages of protein folding on the ribosome to the wealth of information available for the refolding of full-length polypeptide chains in vitro.

Computer simulations and theory of protein translocation

The translocation of proteins through pores is central to many biological phenomena, such as mitochondrial protein import, protein degradation, and delivery of protein toxins to their cytosolic targets. Because proteins typically have to pass through constrictions that are too narrow to accommodate folded structures, translocation must be coupled to protein unfolding. The simplest model that accounts for such co-translocational unfolding assumes that both translocation and unfolding are accomplished by pulling on the end of the polypeptide chain mechanically. In this Account, we describe theoretical studies and computer simulations of this model and discuss how the time scales of translocation depend on the pulling force and on the protein structure. Computationally, this is a difficult problem because biologically or experimentally relevant time scales of translocation are typically orders of magnitude slower than those accessible by fully atomistic simulations. For this reason, we explore one-dimensional free energy landscapes along suitably defined translocation coordinates and discuss various approaches to their computation. We argue that the free energy landscape of translocation is often bumpy because confinement partitions the protein's configuration space into distinct basins of attraction separated by large entropic barriers. Favorable protein-pore interactions and nonnative interactions within the protein further contribute to the complexity. Computer simulations and simple scaling estimates show that forces of just 2-6 pN are often sufficient to ensure transport of unstructured polypeptides, whereas much higher forces are typically needed to translocate folded protein domains. The unfolding mechanisms found from simulations of translocation are different from those observed in the much better understood case of atomic force microscopy (AFM) pulling studies, in which proteins are unraveled by stretching them between their N- and C-termini. In contrast to AFM experiments, single-molecule experimental studies of protein translocation have just started to emerge. We describe one example of a collaborative study, in which dwell times of beta-hairpin-forming peptides inside the alpha-hemolysin pore were both measured experimentally and estimated using computer simulations. Analysis of the simulated trajectories has explained the experimental finding that more stable hairpins take, on the average, longer to traverse the pore. Despite the insight we have gained, the general relationship between the structure of proteins and their resistance to mechanically driven co-translocational unfolding remains poorly understood. Future theoretical progress likely will be made in conjunction with single-molecule experiments and will require realistic models to account for specific protein-pore interactions and for solvent effects.

A pause for thought along the co-translational folding pathway

A unifying concept that combines the basic features governing self-organization of proteins into complex three-dimensional structures in vitro and in vivo is still lacking. Recent experimental results and theoretical in silico modeling studies provide evidence showing that mRNA might contain an additional layer of information, beyond the amino acid sequence, that fine-tunes in vivo protein folding, which is largely believed to start as a co-translational process. These findings indicate that translation kinetics might direct the co-translational folding pathway and that translational pausing at rare codons might provide a time delay to enable independent and sequential folding of the defined portions of the nascent polypeptide emerging from the ribosome.

Excursion of a single polypeptide into a protein pore: simple physics, but complicated biology

Despite its fundamental and critical importance in molecular biology and practical medical biotechnology, how a polypeptide interacts with a transmembrane protein pore is not yet comprehensively understood. Here, we employed single-channel electrical recordings to reveal the interactions of short polypeptides and small folded proteins with a robust beta-barrel protein pore. The short polypeptides were approximately 25 residues in length, resembling positively charged targeting presequences involved in protein import. The proteins were consisted of positively charged pre-cytochrome b2 fragments (pb2) fused to the small ribonuclease barnase (approximately 110 residues, Ba). Single-molecule experiments exploring the interaction of a folded pb2-Ba protein with a single beta-barrel pore, which contained negatively charged electrostatic traps, revealed the complexity of a network of intermolecular forces, including driving and electrostatic ones. In addition, the interaction was dependent on other factors, such as the hydrophobic content of the interacting polypeptide, the location of the electrostatic trap, the length of the pb2 presequence and temperature. This single-molecule approach together with protein design of either the interacting polypeptide or the pore lumen opens new opportunities for the exploration of the polypeptide-pore interaction at high temporal resolution. Such future studies are also expected to unravel the advantages and limitations of the nanopore technique for the detection and exploration of individual polypeptides.

Multi-scale simulations provide supporting evidence for the hypothesis of intramolecular protein translocation in GroEL/GroES complexes

The biological function of chaperone complexes is to assist the folding of non-native proteins. The widely studied GroEL chaperonin is a double-barreled complex that can trap non-native proteins in one of its two barrels. The ATP-driven binding of a GroES cap then results in a major structural change of the chamber where the substrate is trapped and initiates a refolding attempt. The two barrels operate anti-synchronously. The central region between the two barrels contains a high concentration of disordered protein chains, the role of which was thus far unclear. In this work we report a combination of atomistic and coarse-grained simulations that probe the structure and dynamics of the equatorial region of the GroEL/GroES chaperonin complex. Surprisingly, our simulations show that the equatorial region provides a translocation channel that will block the passage of folded proteins but allows the passage of secondary units with the diameter of an alpha-helix. We compute the free-energy barrier that has to be overcome during translocation and find that it can easily be crossed under the influence of thermal fluctuations. Hence, strongly non-native proteins can be squeezed like toothpaste from one barrel to the next where they will refold. Proteins that are already fairly close to the native state will not translocate but can refold in the chamber where they were trapped. Several experimental results are compatible with this scenario, and in the case of the experiments of Martin and Hartl, intra chaperonin translocation could explain why under physiological crowding conditions the chaperonin does not release the substrate protein.

In silico protein fragmentation reveals the importance of critical nuclei on domain reassembly

Protein complementation assays (PCAs) based on split protein fragments have become powerful tools that facilitate the study and engineering of intracellular protein-protein interactions. These assays are based on the observation that a given protein can be split into two inactive fragments and these fragments can reassemble into the original properly folded and functional structure. However, one experimentally observed limitation of PCA systems is that the folding of a protein from its fragments is dramatically slower relative to that of the unsplit parent protein. This is due in part to a poor understanding of how PCA design parameters such as split site position in the primary sequence and size of the resulting fragments contribute to the efficiency of protein reassembly. We used a minimalist on-lattice model to analyze how the dynamics of the reassembly process for two model proteins was affected by the location of the split site. Our results demonstrate that the balanced distribution of the "folding nucleus," a subset of residues that are critical to the formation of the transition state leading to productive folding, between protein fragments is key to their reassembly.

Folding kinetics of a lattice protein via a forward flux sampling approach

We implement a forward flux sampling approach [R. J. Allen et al., J. Chem. Phys. 124, 194111 (2006)] for calculating transition rate constants and for sampling paths of protein folding events. The algorithm generates trajectories for the transition between the unfolded and folded states as chains of partially connected paths, which can be used to obtain the transition-state ensemble and the properties that characterize these intermediates. We apply this approach to Monte Carlo simulations of a model lattice protein in open space and in confined spaces of varying dimensions. We study the effect of confinement on both protein thermodynamic stability and folding kinetics; the former by mapping free-energy landscapes and the latter by the determination of rate constants and mechanistic details of the folding pathway. Our results show that, for the range of temperatures where the native state is stable, confinement of a protein destabilizes the unfolded state by reducing its entropy, resulting in increased thermodynamic stability of the folded state. Relative to the folding in open space, we find that the kinetics can be accelerated at temperatures above the temperature at which the unconfined protein folds fastest and that the rate constant increases with the number of constrained dimensions. By examining the statistical properties of the transition-state ensemble, we detect signs of a classical nucleation folding mechanism for a core of native contacts formed at an early stage of the process. This nucleus acts as folding foci and is composed of those residues that have higher probability to form native contacts in the transition-state intermediates, which can vary depending on the confinement conditions of the system.

Single-molecule electrophoresis of beta-hairpin peptides by electrical recordings and Langevin dynamics simulations

We used single-channel electrical recordings and Langevin molecular dynamics simulations to explore the electrophoretic translocation of various beta-hairpin peptides across the staphylococcal alpha-hemolysin (alphaHL) protein pore at single-molecule resolution. The beta-hairpin peptides, which varied in their folding properties, corresponded to the C terminal residues of the B1 domain of protein G. The translocation time was strongly dependent on the electric force and was correlated with the folding features of the beta-hairpin peptides. Highly unfolded peptides entered the pore in an extended conformation, resulting in fast single-file translocation events. In contrast, the translocation of the folded beta-hairpin peptides occurred more slowly. In this case, the beta-hairpin peptides traversed the alphaHL pore in a misfolded or fully folded conformation. This study demonstrates that the interaction between a polypeptide and a beta-barrel protein pore is dependent on the folding features of the polypeptide.