Non-native Folding Research Articles

Stabilization of Non-Native Folds and Programmable Protein Gelation in Compositionally Designed Deep Eutectic Solvents.

Proteins are adjustable units from which biomaterials with designed properties can be developed. However, non-native folded states with controlled topologies are hardly accessible in aqueous environments, limiting their prospects as building blocks. Here, we demonstrate the ability of a series of anhydrous deep eutectic solvents (DESs) to precisely control the conformational landscape of proteins. We reveal that systematic variations in the chemical composition of binary and ternary DESs dictate the stabilization of a wide range of conformations, that is, compact globular folds, intermediate folding states, or unfolded chains, as well as controlling their collective behavior. Besides, different conformational states can be visited by simply adjusting the composition of ternary DESs, allowing for the refolding of unfolded states and vice versa. Notably, we show that these intermediates can trigger the formation of supramolecular gels, also known as eutectogels, where their mechanical properties correlate to the folding state of the protein. Given the inherent vulnerability of proteins outside the native fold in aqueous environments, our findings highlight DESs as tailorable solvents capable of stabilizing various non-native conformations on demand through solvent design.

Intrinsically disordered proteins/regions and insight into their biomolecular interactions

Proteins may vary from being rigid to having flexible regions to being completely disordered, either as an intrinsically disordered protein (IDP) or having specific intrinsically disordered regions (IDRs). IDPs/IDRs can form complexes otherwise impossible, such as wrapping around the binding partner, hence providing the plasticity needed for achieving assemblies with specific functions. IDRs can exhibit promiscuity, using the same region in the sequence to bind multiple partners, and act as hubs in protein-protein interaction network (an essential part of the cell signalling network). Disorder-to-order transition on binding provides specificity with affinity, optimum for reversibility of the binding, thus offering suitability for regulation and signalling processes. IDRs interactions may be modulated by the environment or covalent modifications; mis-signalling or their unnatural or non-native folding may lead to diseases. This article aims to provide an overview of structural heterogeneity, as seen in IDPs/IDRs, and their role in biological recognition, binding and function.

GrpEL1 regulates mitochondrial unfolded protein response after experimental subarachnoid hemorrhage in vivo and in vitro

Subarachnoid hemorrhage (SAH) is a hemorrhagic stroke disease with high mortality and disability rates. Neurological recovery in early brain injury (EBI) after SAH is a crucial stage to reduce complications and improve the prognosis of patients. The mitochondrial unfolded protein response (UPRmt) is an essential mitochondrial damage repair process, that degrades aggresomes formed by misfolded proteins. UPRmt is a response to cellular stress and enhances mitochondrial homeostasis. GrpEL1 is a nucleotide exchange factor that assists mtHSP70 in nonnative folding proteins in mitochondria. However, the role of UPRmt and GrpEL1 after SAH is unclear. Western blot, Immunofluorescence, Aggresome staining, JC-1 staining were conducted to detect UPRmt after SAH in vivo and in vitro. The results showed that the UPRmt-related proteins HSP60 and mtHSP70 did not change in the EBI after SAH in vivo and in vitro but increased in the isolated mitochondria. In vitro primary neurons treated with oxyhemoglobin (OxyHb) achieved the same result as MG132 induction, increasing neuron protein aggresomes. The expression of GRPEL1 was unchanged in total protein and mitochondrial protein by Western blot. Co-immunoprecipitation (Co-IP) experiments showed that the GRPEL1-mtHSP70 complex decreased after OxyHb treatment. After GRPEL1 overexpression, the GRPEL1-mtHSP70 complex increased, while aggresome in neurons decreased. JC-1 showed an increased mitochondrial membrane potential, ATP content increased, and Western blot analysis revealed decreased cleaved-Caspase 9, suggesting improved mitochondrial function. In conclusion, the reduced GrpEL1-mtHSP70 complex is an essential factor affecting UPRmt in EBI after SAH. Increasing GrpEL1 promotes GrpEL1 and mtHSP70 binding, promoting the neuronal mitochondrial homeostasis, and might be an essential clinical intervention target for EBI after SAH.

Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates

The series of RNA folding events that occur during transcription can critically influence cellular RNA function. Here, we present reconstructing RNA dynamics from data (R2D2), a method to uncover details of cotranscriptional RNA folding. We model the folding of the Escherichia coli signal recognition particle (SRP) RNA and show that it requires specific local structural fluctuations within a key hairpin to engender efficient cotranscriptional conformational rearrangement into the functional structure. All-atom molecular dynamics simulations suggest that this rearrangement proceeds through an internal toehold-mediated strand-displacementmechanism, which can be disrupted with a point mutation that limits local structural fluctuations and rescued with compensating mutations that restore these fluctuations. Moreover, a cotranscriptional folding intermediate could be cleaved invitro by recombinant E.coli RNase P, suggesting potential cotranscriptional processing. These results from experiment-guided multi-scale modeling demonstrate that even an RNA with a simple functional structure can undergo complex folding and processing during synthesis.

The CFTR P67L variant reveals a key role for N-terminal lasso helices in channel folding, maturation, and pharmacologic rescue

Patients with cystic fibrosis (CF) harboring the P67L variant in the cystic fibrosis transmembrane conductance regulator (CFTR) often exhibit a typical CF phenotype, including severe respiratory compromise. This rare mutation (reported in <300 patients worldwide) responds robustly to CFTR correctors, such as lumacaftor and tezacaftor, with rescue in model systems that far exceed what can be achieved for the archetypical CFTR mutant F508del. However, the specific molecular consequences of the P67L mutation are poorly characterized. In this study, we conducted biochemical measurements following low-temperature growth and/or intragenic suppression, which suggest a mechanism underlying P67L that (1) shares key pathogenic features with F508del, including off-pathway (non-native) folding intermediates, (2) is linked to folding stability of nucleotide-binding domains 1 and 2, and (3) demonstrates pharmacologic rescue that requires domains in the carboxyl half of the protein. We also investigated the “lasso” helices 1 and 2, which occur immediately upstream of P67. Based on limited proteolysis, pulse chase, and molecular dynamics analysis of full-length CFTR and a series of deletion constructs, we argue that P67L and other maturational processing (class 2) defects impair the integrity of the lasso motif and confer misfolding of downstream domains. Thus, amino-terminal missense variants elicit a conformational change throughout CFTR that abrogates maturation while providing a robust substrate for pharmacologic repair.

Non-native fold of the putative VPS39 zinc finger domain.

Background: The multi-subunit homotypic fusion and vacuole protein sorting (HOPS) membrane-tethering complex is involved in regulating the fusion of late endosomes and autophagosomes with lysosomes in eukaryotes. The C-terminal regions of several HOPS components have been shown to be required for correct complex assembly, including the C-terminal really interesting new gene (RING) zinc finger domains of HOPS components VPS18 and VPS41. We sought to structurally characterise the putative C-terminal zinc finger domain of VPS39, which we hypothesised may be important for binding of VPS39 to cellular partners or to other HOPS components. Methods: We recombinantly expressed, purified and solved the crystal structure of the proposed zinc-binding region of VPS39. Results: In the structure, this region forms an anti-parallel β-hairpin that is incorporated into a homotetrameric eight-stranded β-barrel. However, the fold is stabilised by coordination of zinc ions by residues from the purification tag and an intramolecular disulphide bond between two predicted zinc ligands. Conclusions: We solved the structure of the VPS39 C-terminal domain adopting a non-native fold. Our work highlights the risk of non-native folds when purifying small zinc-containing domains with hexahistidine tags. However, the non-native structure we observe may have implications for rational protein design.

Non-native fold of the putative VPS39 zinc finger domain

Background: The multi-subunit homotypic fusion and vacuole protein sorting (HOPS) membrane-tethering complex is involved in regulating the fusion of late endosomes and autophagosomes with lysosomes in eukaryotes. The C-terminal regions of several HOPS components have been shown to be required for correct complex assembly, including the C-terminal really interesting new gene (RING) zinc finger domains of HOPS components VPS18 and VPS41. We sought to structurally characterise the putative C-terminal zinc finger domain of VPS39, which we hypothesised may be important for binding of VPS39 to cellular partners or to other HOPS components. Methods: We recombinantly expressed, purified and solved the crystal structure of the proposed zinc-binding region of VPS39. Results: In the structure, this region forms an anti-parallel β-hairpin that is incorporated into a homotetrameric eight-stranded β-barrel. However, the fold is stabilised by coordination of zinc ions by residues from the purification tag and an intramolecular disulphide bond between two predicted zinc ligands. Conclusions: We solved the structure of the VPS39 C-terminal domain adopting a non-native fold. Our work highlights the risk of non-native folds when purifying small zinc-containing domains with hexahistidine tags. However, the non-native structure we observe may have implications for rational protein design.

Non-native fold of the putative VPS39 zinc finger domain

Background: The multi-subunit homotypic fusion and vacuole protein sorting (HOPS) membrane-tethering complex is involved in regulating the fusion of late endosomes and autophagosomes with lysosomes in eukaryotes. The C-terminal regions of several HOPS components have been shown to be required for correct complex assembly, including the C-terminal really interesting new gene (RING) zinc finger domains of HOPS components VPS18 and VPS41. We sought to structurally characterise the putative C-terminal zinc finger domain of VPS39, which we hypothesised may be important for binding of VPS39 to cellular partners or to other HOPS components. Methods: We recombinantly expressed, purified and solved the crystal structure of the proposed zinc-binding region of VPS39. Results: In the structure, this region forms an anti-parallel β-hairpin that is incorporated into a homotetrameric eight-stranded β-barrel. However, the fold is stabilised by coordination of zinc ions by residues from the purification tag and an intramolecular disulphide bond between two predicted zinc ligands. Conclusions: We solved the structure of the VPS39 C-terminal domain adopting a non-native fold. Our work highlights the risk of non-native folds when purifying small zinc-containing domains with hexahistidine tags. However, the non-native structure we observe may have implications for rational protein design.

Transcription Increases the Cooperativity of Ribonucleoprotein Assembly

The synthesis of new ribosomes begins during transcription of the rRNA and is widely assumed to follow an orderly 5' to 3' gradient. To visualize co-transcriptional assembly of ribosomal protein-RNA complexes in real time, we developed a single-molecule platform that simultaneously monitors transcription and protein association with the elongating transcript. Unexpectedly, the early assembly protein uS4 binds newly made pre-16S rRNA only transiently, likely due to non-native folding of the rRNA during transcription. Stable uS4 binding became more probable only in the presence of additional ribosomal proteins that bind upstream and downstream of protein uS4 by allowing productive assembly intermediates to form earlier. We propose that dynamic sampling of elongating RNA by multiple proteins overcomes heterogeneous RNA folding, preventing assembly bottlenecks and initiating assembly within the transcription time window. This may be a common feature of transcription-coupled RNP assembly.

Quantitative Evaluation of Native Protein Folds and Assemblies by Hydrogen Deuterium Exchange Mass Spectrometry (HDX-MS)

Hydrogen deuterium exchange mass spectrometry (HDX-MS) has significant potential for protein structure initiatives but its relationship with protein conformations is unclear. We report on the efficacy of HDX-MS to distinguish between native and non-native proteins using a popular approach to calculate HDX protection factors (PFs) from protein structures. The ability of HDX-MS to identify native protein conformations is quantified by binary structural classification such that merits of the approach for protein modelling can be quantified and better understood. We show that highly accurate PF calculations are not a prerequisite for HDX-MS simulations that are capable of effectively discriminating between native and non-native protein folds. The simulations can also be performed directly on unique structures facilitating high-throughput evaluation of many alternate conformations. The ability of HDX-MS to classify the conformations of homo-protein assemblies is also investigated. In contrast to protein monomers, we show a significant lack of correspondence between the simulated and experimental HDX-MS data for these systems with a subsequent decrease in the ability of HDX-MS to identify native states. However, we demonstrate surprisingly high diagnostic ability of the simulated data for assemblies in which a significant proportion of the individual chains occupy protein-protein interfaces. We relate this to the number of peptides that can sample alternate subunit orientations and discuss these observations within the larger context of applying HDX-MS to evaluate protein structures.Graphical

All-atom four-body knowledge-based statistical potential to distinguish native tertiary RNA structures from nonnative folds

Scientific breakthroughs in recent decades have uncovered the capability of RNA molecules to fulfill a wide array of structural, functional, and regulatory roles in living cells, leading to a concomitantly significant increase in both the number and diversity of experimentally determined RNA three-dimensional (3D) structures. Atomic coordinates from a representative training set of solved RNA structures, displaying low sequence and structure similarity, facilitate derivation of knowledge-based energy functions. Here we develop an all-atom four-body statistical potential and evaluate its capacity to distinguish native RNA 3D structures from nonnative folds based on calculated free energy scores. Atomic four-body nearest-neighbors are objectively identified by their occurrence as tetrahedral vertices in the Delaunay tessellations of RNA structures, and rates of atomic quadruplet interactions expected by chance are obtained from a multinomial reference distribution. Our four-body energy function, referred to as RAMP (ribonucleic acids multibody potential), is subsequently derived by applying the inverted Boltzmann principle to the frequency data, yielding an energy score for each type of atomic quadruplet interaction. Several well-known benchmark datasets reveal that RAMP is comparable with, and often outperforms, existing knowledge- and physics-based energy functions. To the best of our knowledge, this is the first study detailing an RNA tertiary structure-based multibody statistical potential and its comparative evaluation.

An Atomic Four-Body Statistical Potential to Distinguish Native RNA Structures from Nonnative Folds

An Atomic Four-Body Statistical Potential to Distinguish Native RNA Structures from Nonnative Folds

Corrigendum to "All-Atom Four-Body Knowledge-Based Statistical Potentials to Distinguish Native Protein Structures from Nonnative Folds".

[This corrects the article DOI: 10.1155/2017/5760612.].

Partially native intermediates mediate misfolding of SOD1 in single-molecule folding trajectories

Prion-like misfolding of superoxide dismutase 1 (SOD1) is associated with the disease ALS, but the mechanism of misfolding remains unclear, partly because misfolding is difficult to observe directly. Here we study the most misfolding-prone form of SOD1, reduced un-metallated monomers, using optical tweezers to measure unfolding and refolding of single molecules. We find that the folding is more complex than suspected, resolving numerous previously undetected intermediate states consistent with the formation of individual β-strands in the native structure. We identify a stable core of the protein that unfolds last and refolds first, and directly observe several distinct misfolded states that branch off from the native folding pathways at specific points after the formation of the stable core. Partially folded intermediates thus play a crucial role mediating between native and non-native folding. These results suggest an explanation for SOD1’s propensity for prion-like misfolding and point to possible targets for therapeutic intervention.

All-Atom Four-Body Knowledge-Based Statistical Potentials to Distinguish Native Protein Structures from Nonnative Folds.

Recent advances in understanding protein folding have benefitted from coarse-grained representations of protein structures. Empirical energy functions derived from these techniques occasionally succeed in distinguishing native structures from their corresponding ensembles of nonnative folds or decoys which display varying degrees of structural dissimilarity to the native proteins. Here we utilized atomic coordinates of single protein chains, comprising a large diverse training set, to develop and evaluate twelve all-atom four-body statistical potentials obtained by exploring alternative values for a pair of inherent parameters. Delaunay tessellation was performed on the atomic coordinates of each protein to objectively identify all quadruplets of interacting atoms, and atomic potentials were generated via statistical analysis of the data and implementation of the inverted Boltzmann principle. Our potentials were evaluated using benchmarking datasets from Decoys-‘R'-Us, and comparisons were made with twelve other physics- and knowledge-based potentials. Ranking 3rd, our best potential tied CHARMM19 and surpassed AMBER force field potentials. We illustrate how a generalized version of our potential can be used to empirically calculate binding energies for target-ligand complexes, using HIV-1 protease-inhibitor complexes for a practical application. The combined results suggest an accurate and efficient atomic four-body statistical potential for protein structure prediction and assessment.

Markov State Models and tICA Reveal a Nonnative Folding Nucleus in Simulations of NuG2

After reanalyzing simulations of NuG2—a designed mutant of protein G—generated by Lindorff-Larsen et al. with time structure-based independent components analysis and Markov state models as well as performing 1.5 ms of additional sampling on Folding@home, we found an intermediate with a register-shift in one of the β-sheets that was visited along a minor folding pathway. The minor folding pathway was initiated by the register-shifted sheet, which is composed of solely nonnative contacts, suggesting that for some peptides, nonnative contacts can lead to productive folding events. To confirm this experimentally, we suggest a mutational strategy for stabilizing the register shift, as well as an infrared experiment that could observe the nonnative folding nucleus.

Secondary structure encodes a cooperative tertiary folding funnel in the Azoarcus ribozyme.

A requirement for specific RNA folding is that the free-energy landscape discriminate against non-native folds. While tertiary interactions are critical for stabilizing the native fold, they are relatively non-specific, suggesting additional mechanisms contribute to tertiary folding specificity. In this study, we use coarse-grained molecular dynamics simulations to explore how secondary structure shapes the tertiary free-energy landscape of the Azoarcus ribozyme. We show that steric and connectivity constraints posed by secondary structure strongly limit the accessible conformational space of the ribozyme, and that these so-called topological constraints in turn pose strong free-energy penalties on forming different tertiary contacts. Notably, native A-minor and base-triple interactions form with low conformational free energy, while non-native tetraloop/tetraloop–receptor interactions are penalized by high conformational free energies. Topological constraints also give rise to strong cooperativity between distal tertiary interactions, quantitatively matching prior experimental measurements. The specificity of the folding landscape is further enhanced as tertiary contacts place additional constraints on the conformational space, progressively funneling the molecule to the native state. These results indicate that secondary structure assists the ribozyme in navigating the otherwise rugged tertiary folding landscape, and further emphasize topological constraints as a key force in RNA folding.

Optimization of rotamers prior to template minimization improves stability predictions made by computational protein design.

Computational protein design (CPD) predictions are highly dependent on the structure of the input template used. However, it is unclear how small differences in template geometry translate to large differences in stability prediction accuracy. Herein, we explored how structural changes to the input template affect the outcome of stability predictions by CPD. To do this, we prepared alternate templates by Rotamer Optimization followed by energy Minimization (ROM) and used them to recapitulate the stability of 84 protein G domain β1 mutant sequences. In the ROM process, side-chain rotamers for wild-type (WT) or mutant sequences are optimized on crystal or nuclear magnetic resonance (NMR) structures prior to template minimization, resulting in alternate structures termed ROM templates. We show that use of ROM templates prepared from sequences known to be stable results predominantly in improved prediction accuracy compared to using the minimized crystal or NMR structures. Conversely, ROM templates prepared from sequences that are less stable than the WT reduce prediction accuracy by increasing the number of false positives. These observed changes in prediction outcomes are attributed to differences in side-chain contacts made by rotamers in ROM templates. Finally, we show that ROM templates prepared from sequences that are unfolded or that adopt a nonnative fold result in the selective enrichment of sequences that are also unfolded or that adopt a nonnative fold, respectively. Our results demonstrate the existence of a rotamer bias caused by the input template that can be harnessed to skew predictions toward sequences displaying desired characteristics.

Effect of functionalized magnetic MnFe2O4 nanoparticles on fibrillation of human serum albumin.

Pathogenesis of amyloid-related diseases is related to nonnative folding of proteins with the formation of insoluble deposits in the extracellular space of various tissues. Having the unique properties of small size, large surface area, biodegradability, and relative nontoxicity, magnetic nanoparticles have drawn a lot of attention in biomedical applications. Herein, we demonstrate the effect of bare and differently functionalized magnetic MnFe2O4 nanoparticles on fibrillation of human serum albumin in vitro. The process has been monitored using Thioflavin T fluorescence, Congo red binding assay, circular dichroism, fluorescence microscopy, and transmission electron microscopy. From our experimental results, amine functionalized MnFe2O4 nanoparticles are found to inhibit formation of fibrils more effectively than bare ones, while carboxylated nanoparticles do not have a significant effect on fibrillation. This study has explored the prospects of using specific magnetic nanoparticles with appropriate modification to control self-assembly of proteins and may act as a precursor in therapeutic applications.

Strain-Independent Increases of Crystallin Proteins in the Retina of Type 1 Diabetic Rats

Diabetic retinopathy is the leading cause of vision loss in working-age individuals in the United States and is expected to continue growing with the increased prevalence of diabetes. Streptozotocin-induced hyperglycemia in rats is the most commonly used model for diabetic retinopathy. Previous studies have shown that this model can lead to different inflammatory changes in the retina depending on the strain of rat. Our previous work has shown that crystallin proteins, including members of the alpha- and beta/gamma-crystallin subfamilies, are upregulated in the STZ rat retina. Crystallin proteins have been implicated in a number of cellular processes, such as neuroprotection, non-native protein folding and vascular remodeling. In this current study, we have demonstrated that unlike other strain-dependent changes, such as inflammatory cytokines and growth factor levels, in the STZ rat, the protein upregulation of crystallins is consistent across the Brown Norway, Long-Evans and Sprague-Dawley rat strains in the context of diabetes. Taken together, these data illustrate the potential critical role played by crystallins, and especially alpha-crystallins, in the retina in the context of diabetes.