The folding of a family of three-helix bundle proteins: spectrin R15 has a robust folding nucleus, unlike its homologous neighbours

Memory Effects Explain the Fractional Viscosity Dependence of Rates Associated with Internal Friction: Simple Models and Applications to Butane Dihedral Rotation

Barrier-crossing rates of biophysical processes, ranging from simple conformational changes to protein folding, often deviate from the Kramers prediction of an inverse viscosity dependence. In many recent studies, this has been attributed to the presence of internal friction within the system. Previously, we showed that memory-dependent friction arising from the nonequilibrium solvation of a single particle crossing a smooth one-dimensional barrier can also cause such a deviation and be misinterpreted as internal friction. Here we introduce a simple diatom model and show that even in the absence of explicit solvent, internal memory effects arise due to coupling of the reaction coordinate motion with frictionally orthogonal degrees of freedom. This results in a fractional viscosity dependence and a deviation from Kramers' theory, typically attributed to the presence of internal friction. This model therefore mimics several biological processes where a local conformational change of a biomolecule is often influenced by its surroundings. This gives rise to an apparent "internal friction" commonly measured in terms of empirical fitting parameters α and σ. We propose a microscopic measure of this internal friction using Grote-Hynes theory which employs memory-dependent friction. We use butane to demonstrate the effect of coupling strength on the internal friction in realistic systems. This model therefore can serve the purpose of understanding internal friction in biological systems in terms of such coupling.

Probing the Energy Landscape of Spectrin R15 and R16 and the Effects of Non-native Interactions

Understanding the details of a protein folding mechanism can be a challenging and complex task. One system with an interesting folding behavior is the α-spectrin domain, where the R15 folds three-orders of magnitude faster than its homologues R16 and R17, despite having similar structures. The molecular origins that explain these folding rate differences remain unclear, but our previous work revealed that a combined effect produced by non-native interactions could be a reasonable cause for these differences. In this study, we explore further the folding process by identifying the molecular paths, metastable states, and the collective motions that lead these unfolded proteins to their native state conformation. Our results uncovered the differences between the folding pathways for the wild-type R15 and R16 and an R16 mutant. The metastable ensembles that speed down the folding were identified using an energy landscape visualization method (ELViM). These ensembles correspond to similar experimentally reported configurations. Our observations indicate that the non-native interactions are also associated with secondary structure misdocking. This computational methodology can be used as a fast, straightforward protocol for shedding light on systems with unclear folding or conformational traps.

Temperature Dependence of Internal Friction of Peptides

Internal friction is a valuable concept to describe the kinetics of proteins. As is well known, internal friction can be modulated by solvent features (such as viscosity). How can internal friction be affected by environmental temperature? The answer to this question is not evident. In the present work, we approach this problem with simulations on two model peptides. The thermodynamics and relaxation kinetics are characterized through long molecular dynamics simulations, with the viscosity modulated by varying the mass of solvent molecules. Based on the extrapolation to zero viscosity together with scaling of the relaxation time scales, we discover that internal friction is almost invariant at various temperatures. Controlled simulations further support the idea that internal friction is independent of environmental temperature. Comparisons between the two model peptides help us to understand the diverse phenomena in experiments.

Cotranslational folding cooperativity of contiguous domains of α-spectrin

Proteins synthesized in the cell can begin to fold during translation before the entire polypeptide has been produced, which may be particularly relevant to the folding of multidomain proteins. Here, we study the cotranslational folding of adjacent domains from the cytoskeletal protein α-spectrin using force profile analysis (FPA). Specifically, we investigate how the cotranslational folding behavior of the R15 and R16 domains are affected by their neighboring R14 and R16, and R15 and R17 domains, respectively. Our results show that the domains impact each other's folding in distinct ways that may be important for the efficient assembly of α-spectrin, and may reduce its dependence on chaperones. Furthermore, we directly relate the experimentally observed yield of full-length protein in the FPA assay to the force exerted by the folding protein in piconewtons. By combining pulse-chase experiments to measure the rate at which the arrested protein is converted into full-length protein with a Bell model of force-induced rupture, we estimate that the R16 domain exerts a maximal force on the nascent chain of ∼15 pN during cotranslational folding.

The burst-phase folding intermediate of ribonuclease H changes conformation over evolutionary history

The amino acid sequence encodes the energy landscape of a protein. Therefore, we expect evolutionary mutations to change features of the protein energy landscape, including the conformations adopted by a polypeptide as it folds to its native state. Ribonucleases H (RNase H) from Escherichia coli and Thermus thermophilus both fold via a partially folded intermediate in which the core region of the protein (helices A-D and strands 4-5) is structured. Strand 1, however, uniquely contributes to the T. thermophilus RNase H folding intermediate (Icore+1 ), but not the E. coli RNase H intermediate (Icore ) (Rosen & Marqusee, PLoS One 2015). We explore the origin of this difference by characterizing the folding intermediate of seven ancestral RNases H spanning the evolutionary history of these two homologs. Using fragment models with or without strand 1 and FRET probes to characterize the folding intermediate of each ancestor, we find a distinct evolutionary trend across the family-the involvement of strand 1 in the folding intermediate is an ancestral feature that is maintained in the thermophilic lineage and is gradually lost in the mesophilic lineage. Evolutionary sequence changes indeed modulate the conformations present on the folding landscape and altered the folding trajectory of RNase H.

Comparative mechanical unfolding studies of spectrin domains R15, R16 and R17

Spectrins belong to repetitive three-helix bundle proteins that have vital functions in multicellular organisms and are of potential value in nanotechnology. To reveal the unique physical features of repeat proteins we have studied the structural and mechanical properties of three repeats of chicken brain α-spectrin (R15, R16 and R17) at the atomic level under stretching at constant velocities (0.01, 0.05 and 0.1 Å·ps^-1) and constant forces (700 and 900 pN) using molecular dynamics (MD) simulations at T = 300 K. 114 independent MD simulations were performed and their analysis has been done. Despite structural similarity of these domains we have found that R15 is less mechanically stable than R16, which is less stable than R17. This result is in agreement with the thermal unfolding rates. Moreover, we have observed the relationship between mechanical stability, flexibility of the domains and the number of aromatic residues involved in aromatic clusters.

Dystrophin and Spectrin, Two Highly Dissimilar Sisters of the Same Family

Dystrophin and Spectrin are two proteins essential for the organization of the cytoskeleton and for the stabilization of membrane cells. The comparison of these two sister proteins, and with the dystrophin homologue utrophin, enables us to emphasise that, despite a similar topology with common subdomains and a common structural basis of a three-helix coiled-coil, they show a large range of dissimilarities in terms of genetics, cell expression and higher level structural organisation. Interactions with cellular partners, including proteins and membrane phospholipids, also show both strikingly similar and very different behaviours. The differences between dystrophin and spectrin are also illustrated by the large variety of pathological anomalies emerging from the dysfunction or the absence of these proteins, showing that they are keystones in their function of providing a scaffold that sustains cell structure.

Evolutionary trend toward kinetic stability in the folding trajectory of RNases H

Proper folding of proteins is critical to producing the biological machinery essential for cellular function. The rates and energetics of a protein's folding process, which is described by its energy landscape, are encoded in the amino acid sequence. Over the course of evolution, this landscape must be maintained such that the protein folds and remains folded over a biologically relevant time scale. How exactly a protein's energy landscape is maintained or altered throughout evolution is unclear. To study how a protein's energy landscape changed over time, we characterized the folding trajectories of ancestral proteins of the ribonuclease H (RNase H) family using ancestral sequence reconstruction to access the evolutionary history between RNases H from mesophilic and thermophilic bacteria. We found that despite large sequence divergence, the overall folding pathway is conserved over billions of years of evolution. There are robust trends in the rates of protein folding and unfolding; both modern RNases H evolved to be more kinetically stable than their most recent common ancestor. Finally, our study demonstrates how a partially folded intermediate provides a readily adaptable folding landscape by allowing the independent tuning of kinetics and thermodynamics.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.