Folding Kinetics Research Articles

Interconnections between m6A RNA modification, RNA structure, and protein-RNA complex assembly.

Protein-RNA complexes exist in many forms within the cell, from stable machines such as the ribosome to transient assemblies like the spliceosome. All protein-RNA assemblies rely on spatially and temporally coordinated interactions between specific proteins and RNAs to achieve a functional form. RNA folding and structure are often critical for successful protein binding and protein-RNA complex formation. RNA modifications change the chemical nature of a given RNA and often alter its folding kinetics. Both these alterations can affect how and if proteins or other RNAs can interact with the modified RNA and assemble into complexes. N6-methyladenosine (m6A) is the most common base modification on mRNAs and regulatory noncoding RNAs and has been shown to impact RNA structure and directly modulate protein-RNA interactions. In this review, focusing on the mechanisms and available quantitative information, we discuss first how the METTL3/14 m6A writer complex is specifically targeted to RNA assisted by protein-RNA and other interactions to enable site-specific and co-transcriptional RNA modification and, once introduced, how the m6A modification affects RNA folding and protein-RNA interactions.

Protein Folding Stability and Kinetics in Alginate Hydrogels.

Proteins are commonly encapsulated in alginate gels for drug delivery and tissue-engineering applications. However, there is limited knowledge of how encapsulation impacts intrinsic protein properties such as folding stability or unfolding kinetics. Here, we use fast relaxation imaging (FReI) to image protein unfolding in situ in alginate hydrogels after applying a temperature jump. Based on changes in the Förster resonance energy transfer (FRET) response of FRET-labeled phosphoglycerate kinase (PGK), we report the quantitative impact of multiple alginate hydrogel concentrations on protein stability and folding dynamics. The gels stabilize PGK by increasing its melting temperature up to 18.4 °C, and the stabilization follows a nonmonotonic dependence on the alginate density. In situ kinetic measurements also reveal that PGK deviates more from two-state folding behavior in denser gels and that the gel decreases the unfolding rate and accelerates the folding rate of PGK, compared to buffer. Phi-value analysis suggests that the folding transition state of an encapsulated protein is structurally similar to that of folded protein. This work reveals both beneficial and negative impacts of gel encapsulation on protein folding, as well as potential mechanisms contributing to altered stability.

Single-molecule FRET observes opposing effects of urea and TMAO on structurally similar meso- and thermophilic riboswitch RNAs.

Bacteria live in a broad range of environmental temperatures that require adaptations of their RNA sequences to maintain function. Riboswitches are regulatory RNAs that change conformation upon typically binding metabolite ligands to control bacterial gene expression. The paradigmatic small class-I preQ1 riboswitches from the mesophile Bacillus subtilis (Bsu) and the thermophile Thermoanaerobacter tengcongensis (Tte) adopt similar pseudoknot structures when bound to preQ1. Here, we use UV-melting analysis combined with single-molecule detected chemical denaturation by urea to compare the thermodynamic and kinetic folding properties of the two riboswitches, and the urea-countering effects of trimethylamine N-oxide (TMAO). Our results show that, first, the Tte riboswitch is more thermotolerant than the Bsu riboswitch, despite only subtle sequence differences. Second, using single-molecule FRET, we find that urea destabilizes the folded pseudoknot structure of both riboswitches, yet has a lower impact on the unfolding kinetics of the thermodynamically less stable Bsu riboswitch. Third, our analysis shows that TMAO counteracts urea denaturation and promotes folding of both the riboswitches, albeit with a smaller effect on the more stable Tte riboswitch. Together, these findings elucidate how subtle sequence adaptations in a thermophilic bacterium can stabilize a common RNA structure when a new ecological niche is conquered.

High-Accuracy Prediction of Stabilizing Surface Mutations to the Three-Helix Bundle, UBA(1), with EmCAST.

The accurate modeling of energetic contributions to protein structure is a fundamental challenge in computational approaches to protein analysis and design. We describe a general computational method, EmCAST (empirical Cα stabilization), to score and optimize the sequence to the structure in proteins. The method relies on an empirical potential derived from the database of the Cα dihedral angle preferences for all possible four-residue sequences, using the data available in the Protein Data Bank. Our method produces stability predictions that naturally correlate one-to-one with the experimental results for solvent-exposed mutation sites. EmCAST predicted four mutations that increased the stability of a three-helix bundle, UBA(1), from 2.4 to 4.8 kcal/mol by optimizing residues in both helices and turns. For a set of eight variants, the predicted and experimental stabilizations correlate very well (R2 = 0.97) with a slope near 1 and with a 0.16 kcal/mol standard error for EmCAST predictions. Tests against literature data for the stability effects of surface-exposed mutations show that EmCAST outperforms the existing stability prediction methods. UBA(1) variants were crystallized to verify and analyze their structures at an atomic resolution. Thermodynamic and kinetic folding experiments were performed to determine the magnitude and mechanism of stabilization. Our method has the potential to enable the rapid, rational optimization of natural proteins, expand the analysis of the sequence/structure relationship, and supplement the existing protein design strategies.

Appraisal of blob-Based Approaches in the Prediction of Protein Folding Times.

A series of reports published in the last 3 years has illustrated that a blob-based model (BBM) can predict the folding time of proteins from their primary amino acid (aa) sequence based on three simple rules established to characterize the long-range backbone dynamics (LRBD) of racemic polypeptides. The sole use of LRBD to predict protein folding times with the BBM represents a radical departure from all other prediction methods currently applied to determine protein folding times, which rely instead on parameters such as the structure content, folding kinetics, chain length, amino acid properties, or contact topography of proteins. Furthermore, the built-in modularity of the BBM enables the parametrization and inclusion of new phenomena affecting the LRBD of polypeptides, while its conceptual simplicity makes it an interesting new mathematical tool for studying protein folding. However, its novelty implies that its relationship with many other methods used to predict protein folding times has not been well researched. Consequently, the purpose of this report is to uncover the physical phenomena encountered during protein folding that are best described by the BBM through the identification of parameters that have been recognized over the years as being strong predictors for protein folding, such as protein size, topology, structural class, and folding kinetics. This was accomplished by determining the parameters most strongly correlated with the folding times predicted by the BBM. While the BBM in its present form appears to be a good indicator of the folding times of the vast majority of the 195 proteins considered so far, this report finds that it excels for moderately large proteins that are primarily composed of locally formed structural motifs such as α-helices or for proteins that fold in multiple steps. Altogether, these observations based on the use of the BBM support the notion that proteins fold the way they do because the LRBD of polypeptides is mostly driven by the local interactions experienced between aa's within reach of one another.

Unveiling Medin Folding and Dimerization Dynamics and Conformations via Atomistic Discrete Molecular Dynamics Simulations.

Medin is a principal component of localized amyloid found in the vasculature of individuals over 50 years old. Its amyloid aggregation has been linked to endothelial dysfunction and vascular inflammation, contributing to the pathogenesis of various vascular diseases. Despite its significance, the structures of the medin monomer, oligomer, and fibril remain elusive, and the dynamic processes of medin aggregation are not fully understood. In this study, we comprehensively investigated the medin folding and dimerization dynamics and conformations using atomistic discrete molecular dynamics simulations. Our simulation results suggested that the folding initiation of the medin involved the formation of β-sheets around medin30-41 and medin42-50, with subsequent capping of other segments to their β-sheet edges. Medin monomers typically consisted of three or four β-strands, along with a dynamic N-terminal helix. Two isolated medin peptides readily aggregated into a β-sheet-rich dimer, displaying a strong aggregation propensity. Dimerization of medin not only enhanced the β-sheet conformations but also led to the formation of β-barrel oligomers. The aggregation tendencies of medin1-18 and medin19-29 were relatively weak. However, the segments of medin30-41 and medin42-50 played a crucial role as they primarily formed a β-sheet core and facilitated medin1-18 and medin19-29 to form intra- and interpeptide β-sheets. The findings highlight the critical role of the medin30-41 and medin42-50 regions in stabilizing the monomer structure and driving the medin amyloid aggregation. These regions could potentially serve as promising targets for designing antiamyloid inhibitors against amyloid aggregation of medin. Additionally, our study provides a full picture of the monomer conformations and dimerization dynamics for medin, which will help better understand the pathology of medin aggregation.

335 Cryogenic electron microscopy studies of human CFTR folding and conformational dynamics elucidate a novel potentiator binding site in the nucleotide binding domains that supports an alternative channelgating mechanism

335 Cryogenic electron microscopy studies of human CFTR folding and conformational dynamics elucidate a novel potentiator binding site in the nucleotide binding domains that supports an alternative channelgating mechanism

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time.

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

Einstein Model of a Graph to Characterize Protein Folded/Unfolded States.

The folded structures of proteins can be accurately predicted by deep learning algorithms from their amino-acid sequences. By contrast, in spite of decades of research studies, the prediction of folding pathways and the unfolded and misfolded states of proteins, which are intimately related to diseases, remains challenging. A two-state (folded/unfolded) description of protein folding dynamics hides the complexity of the unfolded and misfolded microstates. Here, we focus on the development of simplified order parameters to decipher the complexity of disordered protein structures. First, we show that any connected, undirected, and simple graph can be associated with a linear chain of atoms in thermal equilibrium. This analogy provides an interpretation of the usual topological descriptors of a graph, namely the Kirchhoff index and Randić resistance, in terms of effective force constants of a linear chain. We derive an exact relation between the Kirchhoff index and the average shortest path length for a linear graph and define the free energies of a graph using an Einstein model. Second, we represent the three-dimensional protein structures by connected, undirected, and simple graphs. As a proof of concept, we compute the topological descriptors and the graph free energies for an all-atom molecular dynamics trajectory of folding/unfolding events of the proteins Trp-cage and HP-36 and for the ensemble of experimental NMR models of Trp-cage. The present work shows that the local, nonlocal, and global force constants and free energies of a graph are promising tools to quantify unfolded/disordered protein states and folding/unfolding dynamics. In particular, they allow the detection of transient misfolded rigid states.

Protein–lipid charge interactions control the folding of outer membrane proteins into asymmetric membranes

Biological membranes consist of two leaflets of phospholipid molecules that form a bilayer, each leaflet comprising a distinct lipid composition. This asymmetry is created and maintained in vivo by dedicated biochemical pathways, but difficulties in creating stable asymmetric membranes in vitro have restricted our understanding of how bilayer asymmetry modulates the folding, stability and function of membrane proteins. In this study, we used cyclodextrin-mediated lipid exchange to generate liposomes with asymmetric bilayers and characterize the stability and folding kinetics of two bacterial outer membrane proteins (OMPs), OmpA and BamA. We found that excess negative charge in the outer leaflet of a liposome impedes their insertion and folding, while excess negative charge in the inner leaflet accelerates their folding relative to symmetric liposomes with the same membrane composition. Using molecular dynamics, mutational analysis and bioinformatics, we identified a positively charged patch critical for folding and stability. These results rationalize the well-known ‘positive-outside’ rule of OMPs and suggest insights into the mechanisms that drive OMP folding and assembly in vitro and in vivo.

Structural Modifications and Novel Protein-Binding Sites in Pre-miR-675-Explaining Its Regulatory Mechanism in Carcinogenesis.

Pre-miR-675 is a microRNA expressed from the exon 1 of H19 long noncoding RNA, and the atypical expression of pre-miR-675 has been linked with several diseases and disorders including cancer. To execute its function inside the cell, pre-miR-675 is folded into a particular conformation, which aids in its interaction with several other biological molecules. However, the exact folding dynamics of pre-miR-675 and its protein-binding motifs are currently unknown. Moreover, how H19 lncRNA and pre-miR-675 crosstalk and modulate each other's activities is also unclear. The detailed structural analysis of pre-miR-675 in this study determines its earlier unknown conformation and identifies novel protein-binding sites on pre-miR-675, thus making it an excellent therapeutic target against cancer. Co-folding analysis between H19 lncRNA and pre-miR-675 determine structural transformations in pre-miR-675, thus describing the earlier unknown mechanism of interaction between these two molecules. Comprehensively, this study details the conformation of pre-miR-675 and its protein-binding sites and explains its relationship with H19 lncRNA, which can be interpreted to understand the role of pre-miR-675 in the development and progression of tumorigenesis and designing new therapeutics against cancers.

Effects of sequence-dependent non-native interactions in equilibrium and kinetic folding properties of knotted proteins.

Determining the role of non-native interactions in folding dynamics, kinetics, and mechanisms is a classic problem in protein folding. More recently, this question has witnessed a renewed interest in light of the hypothesis that knotted proteins require the assistance of non-native interactions to fold efficiently. Here, we conduct extensive equilibrium and kinetic Monte Carlo simulations of a simple off-lattice C-alpha model to explore the role of non-native interactions in the thermodynamics and kinetics of three proteins embedding a trefoil knot in their native structure. We find that equilibrium knotted conformations are stabilized by non-native interactions that are non-local, and proximal to native ones, thus enhancing them. Additionally, non-native interactions increase the knotting frequency at high temperatures, and in partially folded conformations below the transition temperatures. Although non-native interactions clearly enhance the efficiency of transition from an unfolded conformation to a partially folded knotted one, they are not required to efficiently fold a knotted protein. Indeed, a native-centric interaction potential drives the most efficient folding transition, provided that the simulation temperature is well below the transition temperature of the considered model system.

Programmable Dual‐Responsive Actuation of Single‐Hydrogel‐Based Bilayer Actuators by Photothermal and Skin Layer Effects with Graphene Oxides

AbstractShape morphing of stimuli‐responsive composite hydrogels has received considerable attention in different research fields. Although various multilayer structures with dissimilar materials are studied to achieve shape morphing, combining swellable hydrogel layers with non‐swellable layers results in issues with interface adhesion and structural integrity. In this study, single‐hydrogel‐based bilayer actuators comprising poly(N‐isopropylacrylamide) (PNIPAM) matrices and graphene oxide (GO)–PNIPAM hinges are presented. Upon temperature rising, the PNIPAM hydrogel acts as the passive layer due to the formation of dense microstructures near the surface (i.e., the skin layer effect), whereas the GO‐PNIPAM hydrogel functions as the active layer, maintaining porous due to structural modification by the presence of GO. Under light exposure, the GO‐PNIPAM hinges experience selective heating due to the photothermal effect of GO. Consequently, the resulting bilayer structures exhibit programmable dual‐responsive 3D shape morphing. Additionally, the folding kinetics of these actuators can be adjusted based on the applied stimulus (temperature changes or light), as they are driven by different mechanisms, the skin layer, or photothermal effects, respectively. Furthermore, the hinge‐based bilayer structures demonstrate walking and steering locomotion by light exposure. This approach can lead to advances in soft robotics, biomimetic systems, and autonomous soft actuators in hydrogel‐based systems.

Characterization of the folding and binding properties of the PTB domain of FRS2 with phosphorylated and unphosphorylated ligands

PTB (PhosphoTyrosine Binding) domains are protein domains that exert their function by binding phosphotyrosine residues on other proteins. They are commonly found in a variety of signaling proteins and are important for mediating protein-protein interactions in numerous cellular processes. PTB domains can also exhibit binding to unphosphorylated ligands, suggesting that they have additional binding specificities beyond phosphotyrosine recognition. Structural studies have reported that the PTB domain from FRS2 possesses this peculiar feature, allowing it to interact with both phosphorylated and unphosphorylated ligands, such as TrkB and FGFR1, through different topologies and orientations. In an effort to elucidate the dynamic and functional properties of these protein-protein interactions, we provide a complete characterization of the folding mechanism of the PTB domain of FRS2 and the binding process to peptides mimicking specific regions of TrkB and FGFR1. By analyzing the equilibrium and kinetics of PTB folding, we propose a mechanism implying the presence of an intermediate along the folding pathway. Kinetic binding experiments performed at different ionic strengths highlighted the electrostatic nature of the interaction with both peptides. The specific role of single amino acids in early and late events of binding was pinpointed by site-directed mutagenesis. These results are discussed in light of previous experimental works on these protein systems.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level.

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibilityand the role of co-transcriptional folding in these processes.

Constitutive strength model for spot joints of thin-walled energy-absorbing elements

The goal of this paper was to determine the strength criterion that causes the failure of point-to-point joints of thin-walled profiles under impact load, as well as to determine the strength parameters that guarantee joint continuity during dynamic folding of crash boxes. For this purpose, extensive quasi-static and dynamic materials tests were conducted for the five most popular steel grades used in the automotive industry: DP600, DP800, DP1000, HSLA and DC01. The strength of the spot joints was investigated in acomplex load state using the original instrumentation. This allowed authors to build and validate mathematical models of the load-bearing capacity of point-to-point joints and to conduct simple and efficient numerical simulations of these elements without the need to account for their geometry. Following the validation of mathematical models and load-bearing models of the joints, extensive experiments of quasi-static and dynamic crushing of thin-walled energy-absorbing top-hat (TH) structures were conducted. The data collected made it possible to calibrate the FEM model correctly and conduct numerical simulations of the dynamic crushing of energy-absorbing structures. As a result, the minimum shear strength (Fs) and cross-tension strength (Fn) of the joints were determined as a function of the sheet thickness and the parameters of Re and Rm. The creation of any spot joint enabling the simultaneous carrying of forces Fs and Fn ensures the preservation of joint continuity throughout the entire crushing process and makes it possible for manufacturers to select the type of spot joint based only on its strength.

OligoR: A Native HDX/MS Data Processing Application Dedicated to Oligonucleotides.

Hydrogen-deuterium exchange mass spectrometry (HDX/MS) is increasingly used to study the dynamics of protein conformation. Coupled to native MS, HDX can also characterize the conformations of oligonucleotides and their binding to cations, small molecules, and proteins. Data processing and visualization of native HDX/MS of oligonucleotides requires dedicated software solutions. OligoR is a web-browser-based application that addresses the specific needs of DNA HDX/MS and native MS experiments from raw data in an open format to visualization and export of results. Whole experiments spanning many time points can be processed in minutes for several mass-separated species. To access valuable folding dynamics information, we have developed a simple and robust approach to deconvolute bimodal isotope distributions, even when they are highly overlapping. This approach is based on modeling physically possible isotope distributions determined from chemical formulae and could be extended to any type of analyte (proteins, peptides, sugars, and small molecules). All results are presented in interactive data tables, and publication-quality figures can be generated, customized, and exported.

Emilio Pérez Piñero: invention through abstraction

Emilio Pérez Piñero represents the paradigm of the architect inventor of the twentieth century. His extraordinary mind developed through a series of special circumstances in the early stages of his life, and further, through his own disciplined dedication, he was able to produce a wealth of work of a richness and complexity that go beyond innovation and problem solving. Invention in Pérez Piñero’s work was the result of a masterly mathematical command and a highly developed three-dimensional vision, which enabled him to see possible solutions that had been considered impossible before him. His geodesic domes are true structural landmarks and his folding reticulated structures, in which dynamic folding mechanisms are designed so as to reach equilibrium when fully stretched at the point at which they become static structures, are the culmination of a body of work that has no comparison in the panorama of the architecture of the time.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function.

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Folding of Staphylococcal Nuclease Induced by Binding of Chemically Modified Substrate Analogues Sheds Light on Mechanisms of Coupled Folding/Binding Reactions.

Several proteins have been shown to undergo a shift in the mechanism of ligand binding-induced folding from conformational selection (CS; folding precedes binding) to induced fit (IF; binding precedes folding) with increasing ligand concentration. In previous studies of the coupled folding/binding reaction of staphylococcal nuclease (SNase) in the presence of a substrate analogue, adenosine-3',5'-diphosphate (prAp), we found that the two phosphate groups make important energetic contributions toward stabilizing its complex with the native protein as well as transient conformational states encountered at high ligand concentrations favoring IF. However, the structural contributions of each phosphate group during the reaction remain unclear. To address this question, we relied on fluorescence, nuclear magnetic resonance (NMR), absorption, and isothermal titration calorimetry to study the effects of deletion of the phosphate groups of prAp on the kinetics of ligand-induced folding, using a strategy analogous to mutational ϕ-value analysis to interpret the results. Kinetic measurements over a wide range of ligand concentrations, together with structural characterization of a transient protein-ligand encounter complex using 2D NMR, indicated that, at high ligand concentrations favoring IF, (i) the 5'-phosphate group interacts weakly with denatured SNase during early stages of the reaction, resulting in loose docking of the two domains of SNase, and (ii) the 3'-phosphate group engages in some specific contacts with the polypeptide in the transition state prior to formation of the native SNase-prAp complex.