Folding Kinetics Research Articles

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.

A General Framework to Interpret Hydrogen-Deuterium Exchange Native Mass Spectrometry of G-Quadruplex DNA.

G-quadruplexes (G4s) are secondary structures formed by guanine-rich oligonucleotides involved in various biological processes. However, characterizing G4s is challenging, because of their structural polymorphism. Here, we establish how hydrogen-deuterium exchange native mass spectrometry (HDX/MS) can help to characterize the G4 structures and dynamics in solution. We correlated the time range of G4 exchange to the number of guanines involved in the inner and outer tetrads. We also established relationships among exchange rates, numbers of tetrads and bound cations, and stability. The use of HDX/native MS allows for the determination of tetrads formed and assessment of G4 stability at a constant temperature. A key finding is that stable G4s exchange through local fluctuations (EX2 exchange), whereas less stable G4s also undergo exchange through partial or complete unfolding (EX1 exchange). Deconvolution of the bimodal isotope distributions resulting from EX1 exchange provides valuable insight into the kinetics of folding and unfolding processes and allows one to detect and characterize transiently unfolded intermediates, even if scarcely populated. HDX/native MS thus represents a powerful tool for a more comprehensive exploration of the folding landscapes of G4s.

Molecular bases for strong phenotypic effects of single synonymous codon substitutions in the E. coli ccdB toxin gene

BackgroundSingle synonymous codon mutations typically have only minor or no effects on gene function. Here, we estimate the effects on cell growth of ~ 200 single synonymous codon mutations in an operonic context by mutating almost all positions of ccdB, the 101-residue long cytotoxin of the ccdAB Toxin-Antitoxin (TA) operon to most degenerate codons. Phenotypes were assayed by transforming the mutant library into CcdB sensitive and resistant E. coli strains, isolating plasmid pools, and subjecting them to deep sequencing. Since autoregulation is a hallmark of TA operons, phenotypes obtained for ccdB synonymous mutants after transformation in a RelE toxin reporter strain followed by deep sequencing provided information on the amount of CcdAB complex formed.ResultsSynonymous mutations in the N-terminal region involved in translation initiation showed the strongest non-neutral phenotypic effects. We observe an interplay of numerous factors, namely, location of the codon, codon usage, t-RNA abundance, formation of anti-Shine Dalgarno sequences, predicted transcript secondary structure, and evolutionary conservation in determining phenotypic effects of ccdB synonymous mutations. Incorporation of an N-terminal, hyperactive synonymous mutation, in the background of the single synonymous codon mutant library sufficiently increased translation initiation, such that mutational effects on either folding or termination of translation became more apparent. Introduction of putative pause sites not only affects the translational rate, but might also alter the folding kinetics of the protein in vivo.ConclusionIn summary, the study provides novel insights into diverse mechanisms by which synonymous mutations modulate gene function. This information is useful in optimizing heterologous gene expression in E. coli and understanding the molecular bases for alteration in gene expression that arise due to synonymous mutations.

Modeling the Effect on a Novel Fungal Peptaibol Placed in an All-Atom Bacterial Membrane Mimicking System via Accelerated Molecular Dynamics Simulations.

We previously reported on a novel peptaibol, named Tripleurin XIIc (TPN), an 18-residue long sequence produced by the fungus Trichoderma pleuroti. We elucidated its 3D structure via classical and accelerated molecular dynamics simulation (aMD) methods and reported the folding dynamics of TPN in water and chloroform solvents. Peptaibols, in general, are insoluble in water, as they are amphipathic and may prefer hydrophobic environments like transmembrane regions. In this study, we attempted to use aMD simulations to model an all-atom bacterial membrane system while placing a TPN molecule in its vicinity. The results highlighted that TPN was able to introduce some disorder into the membrane and caused lipid clustering. It could also enter the transmembrane region from the water-bilayer interface. The structural dynamics of TPN in the transmembrane region revealed a single energetically stable conformation similar to the one obtained from water and chloroform solvent simulations reported by us previously. However, this linear structure was found to be at the local energy minimum (stable) in water but at a metastable intermediate state (higher energy) in chloroform. Therefore, it could be said that the water solvent can be successfully used for folding simulations of peptaibols.

Line Shape Analysis of 19F NMR-Monitored Chemical Denaturation of a Fold-Switching Protein RfaH Reveals Its Slow Folding Dynamics.

The recent discovery of metamorphic proteins, which can switch between multiple conformations under native conditions, has challenged the well-established one sequence-one structure paradigm of protein folding. This is exemplified in the C-terminal domain of the multidomain transcription factor RfaH, which converts from an α-helical coiled-coil conformation in its autoinhibited state to a β-barrel conformation upon activation. Here, we use multisite line shape analysis of 19F NMR-monitored equilibrium chemical denaturation measurements of two 19F-labeled variants of full-length RfaH, to show that it folds/unfolds slowly on the NMR time scale, in an apparent all-or-none fashion at physiological pH and room temperature in the 3.3-4.8 M urea concentration range. The significant thermodynamic stability and slow unfolding rate (kinetic stability) are likely responsible for maintaining the closed autoinhibited state of RfaH, preventing spurious interactions with RNA polymerase (RNAP) when not functional. Our results provide a quantitative understanding of the folding-function relationship in the model fold-switching protein RfaH.