Folding Pathway Research Articles

Kinetic pathway of HIV-1 TAR cotranscriptional folding.

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

The Folding Pathway of ABC Transporter CFTR: Effective and Robust

De novo protein folding into a native three-dimensional structure is indispensable for biological function, is instructed by its amino acid sequence, and occurs along a vectorial trajectory. The human proteome contains thousands of membrane-spanning proteins, whose biosynthesis begins on endoplasmic reticulum-associated ribosomes. Nearly half of all membrane proteins traverse the membrane more than once, including therapeutically important protein families such as solute carriers, G-protein-coupled receptors, and ABC transporters. These mediate a variety of functions like signal transduction and solute transport and are often of vital importance for cell function and tissue homeostasis. Missense mutations in multispan membrane proteins can lead to misfolding and cause disease; an example is the ABC transporter Cystic Fibrosis Transmembrane Conductance Regulator (CFTR). Even though our understanding of multispan membrane-protein folding still is rather rudimental, the cumulative knowledge of 20 years of basic research on CFTR folding has led to development of drugs that modulate the misfolded protein. This has provided the prospect of a life without CF to the vast majority of patients. In this review we describe our understanding of the folding pathway of CFTR in cells, which is modular and tolerates many defects, making it effective and robust. We address how modulator drugs affect folding and function of CFTR, and distinguish protein stability from its folding process. Since the domain architecture of (mammalian) ABC transporters are highly conserved, we anticipate that the insights we discuss here for folding of CFTR may lay the groundwork for understanding the general rules of ABC-transporter folding.

The immune-evasive proline-283 substitution in influenza nucleoprotein increases aggregation propensity without altering the native structure.

Nucleoprotein (NP) is a key structural protein of influenza ribonucleoprotein complexes and is central to viral RNA packing and trafficking. NP also determines the sensitivity of influenza to myxovirus resistance protein 1 (MxA), an innate immunity factor that restricts influenza replication. A few critical MxA-resistant mutations have been identified in NP, including the highly conserved proline-283 substitution. This essential proline-283 substitution impairs influenza growth, a fitness defect that becomes particularly prominent at febrile temperature (39°C) when host chaperones are depleted. Here, we biophysically characterize proline-283 NP and serine-283 NP to test whether the fitness defect is caused by the proline-283 substitution introducing folding defects. We show that the proline-283 substitution changes the folding pathway of NP, making NP more aggregation prone during folding, but does not alter the native structure of the protein. These findings suggest that influenza has evolved to hijack host chaperones to promote the folding of otherwise biophysically incompetent viral proteins that enable innate immune system escape.

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides.

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Quantitative proteomic analysis reveals unique Hsp90 cycle-dependent client interactions.

Hsp90 is an abundant and essential molecular chaperone that mediates the folding and activation of client proteins in a nucleotide-dependent cycle. Hsp90 inhibition directly or indirectly impacts the function of 10-15% of all proteins due to degradation of client proteins or indirect downstream effects. Due to its role in chaperoning oncogenic proteins, Hsp90 is an important drug target. However, compounds that occupy the ATP-binding pocket and broadly inhibit function have not achieved widespread use due to negative effects. More selective inhibitors are needed; however, it is unclear how to achieve selective inhibition. We conducted a quantitative proteomic analysis of soluble proteins in yeast strains expressing wild-type Hsp90 or mutants that disrupt different steps in the client folding pathway. Out of 2,482 proteins in our sample set (approximately 38% of yeast proteins), we observed statistically significant changes in abundance of 350 (14%) of those proteins (log2 fold change ≥ 1.5). Of these, 257/350 (∼73%) with the strongest differences in abundance were previously connected to Hsp90 function. Principal component analysis of the entire dataset revealed that the effects of the mutants could be separated into 3 primary clusters. As evidence that Hsp90 mutants affect different pools of clients, simultaneous co-expression of 2 mutants in different clusters restored wild-type growth. Our data suggest that the ability of Hsp90 to sample a wide range of conformations allows the chaperone to mediate folding of a broad array of clients and that disruption of conformational flexibility results in client defects dependent on those states.

Stimuli-Responsive Control over Self-Assembled Nanostructures in Sequence-Specific Functional Block Copolymers.

The precise sequence of a protein's primary structure is essential in determining its folding pathways. To emulate the complexity of these biomolecules, functional block copolymers consisting of segmented triblocks with distinct functionalities positioned in a sequence-specific manner are designed to control the polymer chain compaction. Triblock polymers P- b -C- b -F and P- b -F- b -C and random diblock copolymer P- b -C- r -F consist of a hydrophilic poly(ethylene oxide) (PEO) block and a hydrophobic block with coumarin (C) and ferrocene (F) moieties that are grafted in a sequence-specific or random manner onto the hydrophilic block. External stimuli such as UVB light, redox, and chemical cues influence the functional hydrophobic block to alter the packing parameters that are monitored with spectroscopic and scattering techniques. Interestingly, the positioning of the stimuli-responsive moiety within the hydrophobic block of P- b -C- b -F, P- b -F- b -C, and P- b -C- r -F affects the extent of the hydrophobic-hydrophilic balance in block copolymers that renders orthogonal control in stimuli-responsive transformation of self-assembled vesicles to micelles.

Integrated analysis of transcriptome and small RNAome reveals regulatory network of rapid and long-term response to heat stress in Rhododendron moulmainense

Main conclusionThe post-transcriptional gene regulatory pathway and small RNA pathway play important roles in regulating the rapid and long-term response of Rhododendron moulmainense to high-temperature stress.The Rhododendron plays an important role in maintaining ecological balance. However, it is difficult to domesticate for use in urban ecosystems due to their strict optimum growth temperature condition, and its evolution and adaptation are little known. Here, we combined transcriptome and small RNAome to reveal the rapid response and long-term adaptability regulation strategies in Rhododendron moulmainense under high-temperature stress. The post-transcriptional gene regulatory pathway plays important roles in stress response, in which the protein folding pathway is rapidly induced at 4 h after heat stress, and alternative splicing plays an important role in regulating gene expression at 7 days after heat stress. The chloroplasts oxidative damage is the main factor inhibiting photosynthesis efficiency. Through WGCNA analysis, we identified gene association patterns and potential key regulatory genes responsible for maintaining the ROS steady-state under heat stress. Finally, we found that the sRNA synthesis pathway is induced under heat stress. Combined with small RNAome, we found that more miRNAs are significantly changed under long-term heat stress. Furthermore, MYBs might play a central role in target gene interaction network of differentially expressed miRNAs in R. moulmainense under heat stress. MYBs are closely related to ABA, consistently, ABA synthesis and signaling pathways are significantly inhibited, and the change in stomatal aperture is not obvious under heat stress. Taken together, we gained valuable insights into the transplantation and long-term conservation domestication of Rhododendron, and provide genetic resources for genetic modification and molecular breeding to improve heat resistance in Rhododendron.

The Mechanism of Folding of Human Frataxin in Comparison to the Yeast Homologue – Broad Energy Barriers and the General Properties of the Transition State

The funneled energy landscape theory suggests that the folding pathway of homologous proteins should converge at the late stages of folding. In this respect, proteins displaying a broad energy landscape for folding are particularly instructive, allowing inferring both the early, intermediate and late stages of folding. In this paper we explore the folding mechanisms of human frataxin, an essential mitochondrial protein linked to the neurodegenerative disorder Friedreich's ataxia. Building upon previous studies on the yeast homologue, the folding pathway of human frataxin is thoroughly examined, revealing a mechanism implying the presence of a broad energy barrier, reminiscent of the yeast counterpart. Through an extensive site-directed mutagenesis, we employed a Φ -value analysis to map native-like contacts in the folding transition state. The presence of a broad energy barrier facilitated the exploration of such contacts in both early and late folding events. We compared results from yeast and human frataxin providing insights into the impact of native topology on the folding mechanism and elucidating the properties of the underlying free energy landscape. The findings are discussed in the context of the funneled energy landscape theory of protein folding.

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.

Multilevel Framework for Analysis of Protein Folding Involving Disulfide Bond Formation.

In this study, a three-layered multicenter ONIOM approach is implemented to characterize the naive folding pathway of bovine pancreatic trypsin inhibitor (BPTI). Each layer represents a distinct level of theory, where the initial layer, encompassing the entire protein, is modeled by a general all-atom force-field GFN-FF. An intermediate electronic structure layer consisting of three multicenter fragments is introduced with the state-of-the-art semiempirical tight-binding method GFN2-xTB. Higher accuracy, specifically addressing the breaking and formation of the three disulfide bonds, is achieved at the innermost layer using the composite DFT method r2SCAN-3c. Our analysis sheds light on the structural stability of BPTI, particularly the significance of interlinking disulfide bonds. The accuracy and efficiency of the multicenter QM/SQM/MM approach are benchmarked using the oxidative formation of cystine. For the folding pathway of BPTI, relative stabilities are investigated through the calculation of free energy contributions for selected intermediates, focusing on the impact of the disulfide bond. Our results highlight the intricate trade-off between accuracy and computational cost, demonstrating that the multicenter ONIOM approach provides a well-balanced and comprehensive solution to describe electronic structure effects in biomolecular systems. We conclude that multiscale energy landscape exploration provides a robust methodology for the study of intriguing biological targets.

Iso-mukaadial acetate and ursolic acid acetate bind to Plasmodium Falciparum heat shock protein 70: towards targeting parasite protein folding pathway.

Plasmodium falciparum is the most lethal malaria parasite. P. falciparum Hsp70 (PfHsp70) is an essential molecular chaperone (facilitates protein folding) and is deemed a prospective antimalarial drug target. The present study investigates the binding capabilities of select plant derivatives, iso-mukaadial acetate (IMA) and ursolic acid acetate (UAA), against P. falciparum using an in silico docking approach. The interaction between the ligands and PfHsp70 was evaluated using plasmon resonance (SPR) analysis. Molecular docking, binding free energy analysis and molecular dynamics simulations were conducted towards understanding the mechanisms by which the compounds bind to PfHsp70. The molecular docking results revealed ligand flexibilities, conformations and positions of key amino acid residues and protein-ligand interactions as crucial factors accounting for selective inhibition of Hsp70. The simulation results also suggest protein-ligand van der Waals forces as the driving force guiding the interaction of these compounds with PfHsp70. Of the two compounds, UAA and IMA bound to PfHsp70 within the micromolar range based on surface plasmon resonance (SPR) based binding assay. Our findings pave way for future rational design of new selective compounds targeting PfHsp70.

The binding selectivity of the C-terminal SH3 domain of Grb2, but not its folding pathway, is dictated by its contiguous SH2 domain

The adaptor protein Grb2, or Growth factor receptor bound protein 2, possesses a pivotal role in the transmission of fundamental molecular signals in the cell. Despite lacking enzymatic activity, Grb2 functions as a dynamic assembly platform, orchestrating intracellular signals through its modular structure. This study delves into the energetic communication of Grb2 domains, focusing on the folding and binding properties of the C-SH3 domain linked to its neighboring SH2 domain. Surprisingly, while the folding and stability of C-SH3 remain robust and unaffected by SH2 presence, significant differences emerge in the binding properties when considered within the tandem context compared to isolated C-SH3. Through a double mutant cycle analysis, we highlighted a subset of residues, located at the interface with the SH2 domain and far from the binding site, finely regulating the binding of a peptide mimicking a physiological ligand of the C-SH3 domain. Our results have mechanistic implications about the mechanisms of specificity of the C-SH3 domain, indicating that the presence of the SH2 domain optimizes binding to its physiological target, and emphasizing the general importance of considering supramodular multidomain protein structures to understand the functional intricacies of protein-protein interaction domains.

AP profiling resolves co-translational folding pathways and chaperone interactions in vivo

AP profiling resolves co-translational folding pathways and chaperone interactions in vivo

From loops to stacks: Early steps in the DNA folding pathway by protamine

From loops to stacks: Early steps in the DNA folding pathway by protamine

Tandem repeats of highly bioluminescent NanoLuc are refolded noncanonically by the Hsp70 machinery.

Chaperones are a large family of proteins crucial for maintaining cellular protein homeostasis. One such chaperone is the 70 kDa heat shock protein (Hsp70), which plays a crucial role in protein (re)folding, stability, functionality, and translocation. While the key events in the Hsp70 chaperone cycle are well established, a relatively small number of distinct substrates were repetitively investigated. This is despite Hsp70 engaging with a plethora of cellular proteins of various structural properties and folding pathways. Here we analyzed novel Hsp70 substrates, based on tandem repeats of NanoLuc (Nluc), a small and highly bioluminescent protein with unique structural characteristics. In previous mechanical unfolding and refolding studies, we have identified interesting misfolding propensities of these Nluc-based tandem repeats. In this study, we further investigate these properties through in vitro bulk experiments. Similar to monomeric Nluc, engineered Nluc dyads and triads proved to be highly bioluminescent. Using the bioluminescence signal as the proxy for their structural integrity, we determined that heat-denatured Nluc dyads and triads can be efficiently refolded by the E. coli Hsp70 chaperone system, which comprises DnaK, DnaJ, and GrpE. In contrast to previous studies with other substrates, we observed that Nluc repeats can be efficiently refolded by DnaK and DnaJ, even in the absence of GrpE co-chaperone. Taken together, our study offers a new powerful substrate for chaperone research and raises intriguing questions about the Hsp70 mechanisms, particularly in the context of structurally diverse proteins.

New insights into the physiological effects of structural environmental enrichment on black sea bream (Acanthopagrus schlegelii) based on a combined analysis of transcriptomics and metabolomics

Environmental enrichment (EE) has been generally accepted as an effective and promising way to enhance the performance and improve the welfare of animals kept in captivity. Evidence for the effects of EE on fish has accumulated rapidly, but little is known about its underlying mechanisms. In this study, black sea bream (Acanthopagrus schlegelii) were reared under enriched and enrichment-free conditions for 30 (post-larval) and 60 (juvenile) days. Physiological, transcriptome, and metabolome data were compared to investigate potential rearing improvements, physiological responses to acute stress, and the underlying mechanisms. The results showed that fish reared under enriched conditions had significantly better growth performance, a higher survival rate, and lower basal cortisol levels. Meanwhile, the cortisol, lactic acid, malondialdehyde (MDA), and total antioxidant capacity (T-AOC) levels in fish subjected to acute stress suggested that EE-reared individuals presented greater stress resistance. However, the effect of EE may decrease with the extension of treatment time for larger fish. The transcriptome analysis showed that energy-related pathways (including oxidative phosphorylation, the citrate cycle, and glycolysis/gluconeogenesis) as well as protein folding, sorting, and degradation-related pathways (protein processing in the endoplasmic reticulum and proteasome pathways); replication and repair related pathways (including pathways involved in nucleotide excision repair, base excision repair, mismatch repair, and homologous recombination); and immune-related pathways (complement and coagulation cascades pathway and the IL-17 signaling pathway) were significantly altered in the EE-reared individuals. In addition, integrated analysis of metabolome and transcriptome data revealed marked enrichment of transcription-associated metabolites involved in cortisol synthesis and secretion, taste transduction, olfactory transduction, and protein digestion and absorption pathways. Our results suggested that EE may influence the growth, fitness, and stress resistance of A. schlegelii by regulating a series of biological processes including energy metabolism, cortisol synthesis, antioxidant enzyme activities, protein folding, DNA replication and repair, and complement and coagulation cascades.

Discovery of a trefoil knot in the RydC RNA: Challenging previous notions of RNA topology

Knots are very common in polymers, including DNA and protein molecules. Yet, no genuine knot has been identified in natural RNA molecules to date. Upon re-examining experimentally determined RNA 3D structures, we discovered a trefoil knot 31, the most basic non-trivial knot, in the RydC RNA. This knotted RNA is a member of a small family of short bacterial RNAs, whose secondary structure is characterized by an H-type pseudoknot. Molecular dynamics simulations suggest a folding pathway of the RydC RNA that starts with a native twisted loop. Based on sequence analyses and computational RNA 3D structure predictions, we postulate that this trefoil knot is a conserved feature of all RydC-related RNAs. The first discovery of a knot in a natural RNA molecule introduces a novel perspective on RNA 3D structure formation and on fundamental research on the relationship between function and spatial structure of biopolymers.

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.

Discordant Effects of Polyamine Depletion by DENSpm and DFMO on β-cell Cytokine Stress and Diabetes Outcomes in Mice.

Type 1 diabetes (T1D) is an autoimmune disease leading to dysfunction and loss of insulin-secreting β cells. In β cells, polyamines have been implicated in causing cellular stress and dysfunction. An inhibitor of polyamine biosynthesis, difluoromethylornithine (DFMO), has been shown to delay T1D in mouse models and preserve β-cell function in humans with recent-onset T1D. Another small molecule, N1,N11-diethylnorspermine (DENSpm), both inhibits polyamine biosynthesis and accelerates polyamine metabolism and is being tested for efficacy in cancer clinical trials. In this study, we show that DENSpm depletes intracellular polyamines as effectively as DFMO in mouse β cells. RNA-sequencing analysis, however, suggests that the cellular responses to DENSpm and DFMO differ, with both showing effects on cellular proliferation but the latter showing additional effects on mRNA translation and protein-folding pathways. In the low-dose streptozotocin-induced mouse model of T1D, DENSpm, unlike DFMO, did not prevent or delay diabetes outcomes but did result in improvements in glucose tolerance and reductions in islet oxidative stress. In nonobese diabetic (NOD) mice, short-term DENSpm administration resulted in a slight reduction in insulitis and proinflammatory Th1 cells in the pancreatic lymph nodes. Longer term treatment resulted in a dose-dependent increase in mortality. Notwithstanding the efficacy of both DFMO and DENSpm in reducing potentially toxic polyamine levels in β cells, our results highlight the discordant T1D outcomes that result from differing mechanisms of polyamine depletion and, more importantly, that toxic effects of DENSpm may limit its utility in T1D treatment.

Distinct patterns of proteostasis network gene expression are associated with different prognoses in melanoma patients

The proteostasis network (PN) is a collection of protein folding and degradation pathways that spans cellular compartments and acts to preserve the integrity of the proteome. The differential expression of PN genes is a hallmark of many cancers, and the inhibition of protein quality control factors is an effective way to slow cancer cell growth. However, little is known about how the expression of PN genes differs between patients and how this impacts survival outcomes. To address this, we applied unbiased hierarchical clustering to gene expression data obtained from primary and metastatic cutaneous melanoma (CM) samples and found that two distinct groups of individuals emerge across each sample type. These patient groups are distinguished by the differential expression of genes encoding ATP-dependent and ATP-independent chaperones, and proteasomal subunits. Differences in PN gene expression were associated with increased levels of the transcription factors, MEF2A, SP4, ZFX, CREB1 and ATF2, as well as markedly different survival outcomes. However, surprisingly, similar PN alterations in primary and metastatic samples were associated with discordant survival outcomes in patients. Our findings reveal that the expression of PN genes demarcates CM patients and highlights several new proteostasis sub-networks that could be targeted for more effective suppression of CM within specific individuals.