Folding Kinetics Research Articles

Vitamin K Epoxide Reductase Complex-Protein Disulphide Isomerase Assemblies in the Thiol-Disulphide Exchange Reactions: Portrayal of Precursor-to-Successor Complexes.

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme that converts vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents supplied by its redox partner through thiol-disulphide exchange reactions. The functionally related molecular complexes assembled during this process have never been described, except for a proposed de novo model of a 'precursor' complex of hVKORC1 associated with protein disulphide isomerase (PDI). Using numerical approaches (in silico modelling and molecular dynamics simulation), we generated alternative 3D models for each molecular complex bonded either covalently or non-covalently. These models differ in the orientation of the PDI relative to hVKORC1 and in the cysteine residue involved in forming protein-protein disulphide bonds. Based on a comparative analysis of these models' shape, folding, and conformational dynamics, the most probable putative complexes, mimicking the 'precursor', 'intermediate', and 'successor' states, were suggested. In addition, we propose using these complexes to develop the 'allo-network drugs' necessary for treating blood diseases.

Mechanism of DNA origami folding elucidated by mesoscopic simulations

Many experimental and computational efforts have sought to understand DNA origami folding, but the time and length scales of this process pose significant challenges. Here, we present a mesoscopic model that uses a switchable force field to capture the behavior of single- and double-stranded DNA motifs and transitions between them, allowing us to simulate the folding of DNA origami up to several kilobases in size. Brownian dynamics simulations of small structures reveal a hierarchical folding process involving zipping into a partially folded precursor followed by crystallization into the final structure. We elucidate the effects of various design choices on folding order and kinetics. Larger structures are found to exhibit heterogeneous staple incorporation kinetics and frequent trapping in metastable states, as opposed to more accessible structures which exhibit first-order kinetics and virtually defect-free folding. This model opens an avenue to better understand and design DNA nanostructures for improved yield and folding performance.

Impact of Hydrodynamic Interactions on the Kinetic Pathway of Protein Folding.

Protein folding is a fundamental process critical to cellular function and human health, but it remains a grand challenge in biophysics. Hydrodynamic interaction (HI) plays a vital role in the self-organization of soft and biological materials, yet its role in protein folding is not fully understood despite folding occurring in a fluid environment. Here, we use the fluid particle dynamics method to investigate many-body hydrodynamic couplings between amino acid residues and fluid motion in the folding kinetics of a coarse-grained four-α-helices bundle protein. Our results reveal that HI helps select fast folding pathways to the native state without being kinetically trapped, significantly speeding up the folding kinetics compared to its absence. First, the directional flow along the protein backbone expedites protein collapse. Then, the incompressibility-induced squeezing flow effects retard the accumulation of non-native hydrophobic contacts, thus preventing the protein from being trapped in local energy minima during the conformational search of the native structure. We also find that the significance of HI in folding kinetics depends on temperature, with a pronounced effect under biologically relevant conditions. Our findings suggest that HI, particularly the short-range squeezing effect, may be crucial in avoiding protein misfolding.

Identical sequences, different behaviors: Protein diversity captured at the single-molecule level

The classical “one sequence, one structure, one function” paradigm has shaped much of our intuition of how proteins work inside the cell. Partially due to the insight provided by bulk biochemical assays, individual biomolecules are often assumed to behave as identical entities, and their characterization relies on ensemble averages that flatten any conformational diversity into a unique phenotype. While the emergence of single-molecule techniques opened the gates to interrogating individual molecules, technical shortcomings typically limit the duration of these measurements, which precludes a complete characterization of an individual protein and, hence, capturing the heterogeneity among molecular populations. Here, we introduce an ultrastable magnetic tweezers design, which enables us to measure the folding dynamics of a single protein during several uninterrupted days with high temporal and spatial resolution. Thanks to this instrumental development, we fully characterize the nanomechanics of two proteins with a very distinct force response, the talin R3IVVI domain and protein L. Days-long recordings on the same protein individual accumulate thousands of folding transitions with submicrosecond resolution, allowing us to reconstruct their free energy landscapes and describe how they evolve with force. By mapping the nanomechanical identity of many different protein individuals, we directly capture their molecular diversity as a quantifiable dispersion on their force response and folding kinetics. By significantly expanding the measurable timescales, our instrumental development offers a tool for profiling individual molecules, opening the gates to directly characterizing biomolecular heterogeneity.

Diffusive dynamics of a model protein chain in solution.

A Markov state model is a powerful tool that can be used to track the evolution of populations of configurations in an atomistic representation of a protein. For a coarse-grained linear chain model with discontinuous interactions, the transition rates among states that appear in the Markov model when the monomer dynamics is diffusive can be determined by computing the relative entropy of states and their mean first passage times, quantities that are unchanged by the specification of the energies of the relevant states. In this paper, we verify the folding dynamics described by a diffusive linear chain model of the crambin protein in three distinct solvent systems, each differing in complexity: a hard-sphere solvent, a solvent undergoing multi-particle collision dynamics, and an implicit solvent model. The predicted transition rates among configurations agree quantitatively with those observed in explicit molecular dynamics simulations for all three solvent models. These results suggest that the local monomer-monomer interactions provide sufficient friction for the monomer dynamics to be diffusive on timescales relevant to changes in conformation. Factors such as structural ordering and dynamic hydrodynamic effects appear to have minimal influence on transition rates within the studied solvent densities.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics.

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

The Biochemical Impact of Extracting an Embedded Adenylate Kinase Domain Using Circular Permutation.

Adenylate kinases (AKs) have evolved AMP-binding and lid domains that are encoded as continuous polypeptides embedded at different locations within the discontinuous polypeptide encoding the core domain. A prior study showed that AK homologues of different stabilities consistently retain cellular activity following circular permutation that splits a region with high energetic frustration within the AMP-binding domain into discontinuous fragments. Herein, we show that mesophilic and thermophilic AKs having this topological restructuring retain activity and substrate-binding characteristics of the parental AK. While permutation decreased the activity of both AK homologues at physiological temperatures, the catalytic activity of the thermophilic AK increased upon permutation when assayed >30 °C below the melting temperature of the native AK. The thermostabilities of the permuted AKs were uniformly lower than those of native AKs, and they exhibited multiphasic unfolding transitions, unlike the native AKs, which presented cooperative thermal unfolding. In addition, proteolytic digestion revealed that permutation destabilized each AK in differing manners, and mass spectrometry suggested that the new termini within the AMP-binding domain were responsible for the increased proteolysis sensitivity. These findings illustrate how changes in contact order can be used to tune enzyme activity and alter folding dynamics in multidomain enzymes.

Design of Functional Globular β-Sheet Miniproteins.

The design of discrete β-sheet peptides is far less advanced than e. g. the design of α-helical peptides. The reputation of β-sheet peptides as being poorly soluble and aggregation-prone often hinders active design efforts. Here, we show that this reputation is unfounded. We demonstrate this by looking at the β-hairpin and WW domain. Their structure and folding have been extensively studied and they have long served as model systems to investigate protein folding and folding kinetics. The resulting fundamental understanding has led to the development of hyperstable β-sheet scaffolds that fold at temperatures of 100 °C or high concentrations of denaturants. These have been used to design functional miniproteins with protein or nucleic acid binding properties, in some cases with such success that medical applications are conceivable. The β-sheet scaffolds are not always completely rigid, but can be specifically designed to respond to changes in pH, redox potential or presence of metal ions. Some engineered β-sheet peptides also exhibit catalytic properties, although not comparable to those of natural proteins. Previous reviews have focused on the design of stably folded and non-aggregating β-sheet sequences. In our review, we now also address design strategies to obtain functional miniproteins from β-sheet folding motifs.

Sequence-based prediction of the effects of histones H1 post-translational modifications: impact on the features related to the function

Post-translational modifications modulate histones H1 activity but their impact on proteins features was not studied so far. Therefore, this work was intended to answer how the most common modifications, i.e. acetylation, methylation, phosphorylation and ubiquitination, can influence on histones H1 to alter their physicochemical and molecular properties. Investigations were done with the use of sequence-based predictors trained on various protein features. Because a full set of histones H1 modifications is not included in the databases of histone proteins, the survey was performed on the human, animals, plants, fungi and protist sequences selected from UniProtKB/Swiss-Prot database. Quantitative proportions of modifications were similar between the groups of organisms (CV = 0.11) but different within the group (p < 0.05). The effects of modifications were evaluated with the use of mutated sequences obtained through the substitution of modified residue of Lys, Ser and Thr by a neutral residue of the Ala. An advantage of deleterious mutations at the sites of acetylation, methylation and ubiquitination over the sites of phosphorylation (p < 0.05) indicate that this modification have more redundant character. Modifications evoke an increase of protein solubility and stability as well as acceleration of folding kinetics and a weaken of binding affinity. Besides, they also maintain a higher extent of intrinsic structural disorder. The obtained results prove that modifications should be perceived as relevant factors influencing physicochemical features determining molecular properties. Thus, histones H1 functioning is strictly correlated with the status of modifications. Communicated by Ramaswamy H. Sarma

Folding dynamics of histone replacement by protamine

Folding dynamics of histone replacement by protamine

4-Thiaproline accelerates the slow folding phase of proteins containing cis prolines in the native state by two orders of magnitude.

The cis/trans isomerization of peptidyl-prolyl peptide bonds is often the bottleneck of the refolding reaction for proteins containing cis proline residues in the native state. Proline (Pro) analogues, especially C4-substituted fluoroprolines, have been widely used in protein engineering to enhance the thermodynamic stability of peptides and proteins and to investigate folding kinetics. 4-thiaproline (Thp) has been shown to bias the ring pucker of Pro, to increase the cis population percentage of model peptides in comparison to Pro, and to diminish the activation energy barrier for the cis/trans isomerization reaction. Despite its intriguing properties, Thp has been seldom incorporated into proteins. Moreover, the impact of Thp on the folding kinetics of globular proteins has never been reported. In this study, we show that upon incorporation of Thp at cisPro76 into the thioredoxin variant Trx1P the half-life of the refolding reaction decreased from ~2 h to ~35 s. A dramatic acceleration of the refolding rate could be observed also for the protein pseudo wild-type barstar upon replacement of cisPro48 with Thp. Quantum chemical calculations suggested that the replacement of the Cγ H2 group by a sulfur atom in the pyrrolidine ring, might lower the barrier for cis/trans rotation due to a weakened peptide bond. The protein variants retained their thermodynamic stability upon incorporation of Thp, while the catalytic and enzymatic activities of the modified Trx1P remained unchanged. Our results show that the Pro isostere Thp might accelerate the rate of the slow refolding reaction for proteins containing cis proline residues in the native state, independent from the local structural environment.

Structural effects of inosine substitution in telomeric DNA quadruplex.

The telomeric DNA, a distal region of eukaryotic chromosome containing guanine-rich repetitive sequence of (TTAGGG)n, has been shown to adopt higher-order structures, specifically G-quadruplexes (G4s). Previous studies have demonstrated the implication of G4 in tumor inhibition through chromosome maintenance and manipulation of oncogene expression featuring their G-rich promoter regions. Besides higher order structures, several regulatory roles are attributed to DNA epigenetic markers. In this work, we investigated how the structural dynamics of a G-quadruplex, formed by the telomeric sequence, is affected by inosine, a prevalent modified nucleotide. We used the standard (TTAGGG)n telomere repeats with guanosine mutated to inosine at each G position. Sequences (GGG)4, (IGG)4, (GIG)4, (GGI)4, (IGI)4, (IIG)4, (GII)4, and (III)4, bridged by TTA linker, are studied using biophysical experiments and molecular modeling. The effects of metal cations in quadruplex folding were explored in both Na+ and K+ containing buffers using CD and UV-melting studies. Our results show that antiparallel quadruplex topology forms with the native sequence (GGG)4 and the terminal modified DNAs (IGG)4 and (GGI)4 in both Na+ and K+ containing buffers. Specifically, quadruplex hybrid was observed for (GGG)4 in K+ buffer. Among the other modified sequences, (GIG)4, (IGI)4 and (GII)4 show parallel features, while (IIG)4 and (III)4 show no detectable conformation in the presence of either Na+ or K+. Our studies indicate that terminal lesions (IGG)4 and (GGI)4 may induce certain unknown conformations. The folding dynamics become undetectable in the presence of more than one inosine substitution except (IGI)4 in both buffer ions. In addition, both UV melting and CD melting studies implied that in most cases the K+ cation confers more thermodynamic stability compared to Na+. Collectively, our conformational studies revealed the diverse structural polymorphisms of G4 with position dependent G-to-I mutations in different ion conditions.

A TeZla micromixer for interrogating the early and broad folding landscape of G-quadruplex via multistage velocity descending

Biomacromolecular folding kinetics involves fast folding events and broad timescales. Current techniques face limitations in either the required time resolution or the observation window. In this study, we developed the TeZla micromixer, integrating Tesla and Zigzag microstructures with a multistage velocity descending strategy. TeZla achieves a significant short mixing dead time (40 µs) and a wide time window covering four orders of magnitude (up to 300 ms). Using this unique micromixer, we explored the folding landscape of c-Myc G4 and its noncanonical-G4 derivatives with different loop lengths or G-vacancy sites. Our findings revealed that c-Myc can bypass folding intermediates and directly adopt a G4 structure in the cation-deficient buffer. Moreover, we found that the loop length and specific G-vacancy site could affect the folding pathway and significantly slow down the folding rates. These results were also cross-validated with real-time NMR and circular dichroism. In conclusion, TeZla represents a versatile tool for studying biomolecular folding kinetics, and our findings may ultimately contribute to the design of drugs targeting G4 structures.

Structural Antagonism-Aided Conformational Regulation Enables an Aptamer-Loop G-Quadruplex Modular Sensor of β-Lactoglobulin.

A simple, reliable method for identifying β-lactoglobulin (β-LG) in dairy products is needed to protect those with β-LG allergies. A common, practical strategy for target detection is designing simplified nucleic acid nanodevices by integrating functional components. This work presents a label-free modular β-LG aptasensor consisting of an aptamer-loop G-quadruplex (G4), the working conformation of which is regulated by conformational antagonism to ensure respective module functionality and the related signal transduction. The polymorphic conformations of the module-fused sequence are systematically characterized, and the cause is revealed as shifting antagonistic equilibrium. Combined with conformational folding dynamics, this helped regulate functional conformations by fine-tuning the sequences. Furthermore, the principle of specific β-LG detection by parallel G4 topology is examined as binding on the G4 aptamer loop by β-LG to reinforce the G4 topology and fluorescence. Finally, a label-free, assembly-free, succinct, and turn-on fluorescent aptasensor is established, achieving excellent sensitivity across five orders of magnitude, rapidly detecting β-LG within 22-min. This study provides a generalizable approach for the conformational regulation of module-fused G4 sequences and a reference model for creating simplified sensing devices for a variety of targets.

Quantifying the energy landscape in weakly and strongly disordered frictional media.

We investigate the "roughness" of the energy landscape of a system that diffuses in a heterogeneous medium with a random position-dependent friction coefficient α(x). This random friction acting on the system stems from spatial inhomogeneity in the surrounding medium and is modeled using the generalized Caldira-Leggett model. For a weakly disordered medium exhibiting a Gaussian random diffusivity D(x) = kBT/α(x) characterized by its average value ⟨D(x)⟩ and a pair-correlation function ⟨D(x1)D(x2)⟩, we find that the renormalized intrinsic diffusion coefficient is lower than the average one due to the fluctuations in diffusivity. The induced weak internal friction leads to increased roughness in the energy landscape. When applying this idea to diffusive motion in liquid water, the dissociation energy for a hydrogen bond gradually approaches experimental findings as fluctuation parameters increase. Conversely, for a strongly disordered medium (i.e., ultrafast-folding proteins), the energy landscape ranges from a few to a few kcal/mol, depending on the strength of the disorder. By fitting protein folding dynamics to the escape process from a metastable potential, the decreased escape rate conceptualizes the role of strong internal friction. Studying the energy landscape in complex systems is helpful because it has implications for the dynamics of biological, soft, and active matter systems.

Locus folding mechanisms determine modes of antigen receptor gene assembly.

The dynamic folding of genomes regulates numerous biological processes, including antigen receptor (AgR) gene assembly. We show that, unlike other AgR loci, homotypic chromatin interactions and bidirectional chromosome looping both contribute to structuring Tcrb for efficient long-range V(D)J recombination. Inactivation of the CTCF binding element (CBE) or promoter at the most 5'Vβ segment (Trbv1) impaired loop extrusion originating locally and extending to DβJβ CBEs at the opposite end of Tcrb. Promoter or CBE mutation nearly eliminated Trbv1 contacts and decreased RAG endonuclease-mediated Trbv1 recombination. Importantly, Trbv1 rearrangement can proceed independent of substrate orientation, ruling out scanning by DβJβ-bound RAG as the sole mechanism of Vβ recombination, distinguishing it from Igh. Our data indicate that CBE-dependent generation of loops cooperates with promoter-mediated activation of chromatin to juxtapose Vβ and DβJβ segments for recombination through diffusion-based synapsis. Thus, the mechanisms that fold a genomic region can influence molecular processes occurring in that space, which may include recombination, repair, and transcriptional programming.

Ligand-dependent folding and unfolding dynamics and free energy landscapes of acylphosphatase.

Acylphosphatase (AcP) is an enzyme which catalyses the hydrolysis of acylphosphate. The binding with the phosphate ion (Pi) assumes significance in preserving both the stability and enzymatic activity of AcP. While previous studies using single molecule force spectroscopy explored the mechanical properties of AcP, the influence of Pi on its folding and unfolding dynamic behaviors remains unexplored. In this work, using stable magnetic tweezers, we measured and compared the force-dependent folding and unfolding rates of AcP in the Tris buffer and phosphate buffer within a force range from 2 pN to 40 pN. We found that Pi exerts no discernible effect on the folding dynamics but consistently decreases the force-dependent unfolding rate of AcP by a constant ratio across the entire force spectrum. The free energy landscapes of AcP in the absence and presence of Pi are constructed. Our results reveal that Pi selectively binds to the native state of AcP, stabilizing it and suggesting the general properties of specific ligand-receptor interactions.

Effects of altered backbone composition on the folding kinetics and mechanism of an ultrafast-folding protein.

Sequence-encoded protein folding is a ubiquitous biological process that has been successfully engineered in a range of oligomeric molecules with artificial backbone chemical connectivity. A remarkable aspect of protein folding is the contrast between the rapid rates at which most sequences in nature fold and the vast number of conformational states possible in an unfolded chain with hundreds of rotatable bonds. Research efforts spanning several decades have sought to elucidate the fundamental chemical principles that dictate the speed and mechanism of natural protein folding. In contrast, little is known about how protein mimetic entities transition between an unfolded and folded state. Here, we report effects of altered backbone connectivity on the folding kinetics and mechanism of the B domain of Staphylococcal protein A (BdpA), an ultrafast-folding sequence. A combination of experimental biophysical analysis and atomistic molecular dynamics simulations performed on the prototype protein and several heterogeneous-backbone variants reveal the interplay among backbone flexibility, folding rates, and structural details of the transition state ensemble. Collectively, these findings suggest a significant degree of plasticity in the mechanisms that can give rise to ultrafast folding in the BdpA sequence and provide atomic level insights into how protein mimetic chains adopt an ordered folded state.

Deciphering the influence of Y12L and N17H substitutions on the conformation and oligomerization of human calcitonin.

The abnormal aggregation of human calcitonin (hCT) hormone peptides impairs their physiological function, leading to harmful immune responses and cytotoxicity, which limits their clinical utility. Interestingly, a representative hCT analog incorporating Y12L and N17H substitutions (DM-hCT) has shown reduced aggregation tendencies while maintaining bioactivity. But the molecular mechanism of Y12L and N17H substitutions on the conformational dynamics of hCT remains unclear. Here, we systematically investigated the folding and self-assembly dynamics of hCT and DM-hCT using atomistic discrete molecular dynamics (DMD) simulations. Our findings revealed that hCT monomers predominantly adopted unstructured conformations with dynamic helices. Oligomerization of hCT resulted in the formation of β-sheet-rich aggregates and β-barrel intermediates. The Y12L and N17H substitutions enhanced helical conformations and suppressed β-sheet formation in both monomers and oligomers. These substitutions stabilized the dynamic helices and disrupted aromatic interactions responsible for β-sheet formation at residue 12. Notably, DM-hCT assemblies still exhibited β-sheets in phenylalanine-rich and C-terminal hydrophobic regions, suggesting that future optimizations should focus on these areas. Our simulations provide insights into the molecular mechanisms underlying hCT aggregation and the amyloid-resistant effects of Y12L and N17H substitutions. These findings have valuable implications for the development of clinical hCT analogs.

Kinetic FRET Assay to Measure Binding-Induced Conformational Changes of Nucleic Acids.

The interaction of small molecules or proteins with RNA or DNA often involves changes in the nucleic acid (NA) folding and structure. A biophysical characterization of these processes helps us to understand the underlying molecular mechanisms. Here, we propose kinFRET (kinetics Förster resonance energy transfer), a real-time ensemble FRET methodology to measure binding and folding kinetics. With kinFRET, the kinetics of conformational changes of NAs (DNA or RNA) upon analyte binding can be directly followed via a FRET signal using a chip-based biosensor. We demonstrate the utility of this approach with two representative examples. First, we monitored the conformational changes of different formats of an aptamer (MN19) upon interaction with small-molecule analytes. Second, we characterized the binding kinetics of RNA recognition by tandem K homology (KH) domains of the human insulin-like growth factor II mRNA-binding protein 3 (IMP3), which reveals distinct kinetic contributions of the two KH domains. Our data demonstrate that kinFRET is well suited to study the kinetics and conformational changes of NA-analyte interactions.