The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein

The amyloid formation of the folded segment of a variant of Exenatide (a marketed drug for type-2 diabetes mellitus) was studied by electronic circular dichroism (ECD) and NMR spectroscopy. We found that the optimum temperature for E5 protein amyloidosis coincides with body temperature and requires well below physiological salt concentration. Decomposition of the ECD spectra and its barycentric representation on the folded-unfolded-amyloid potential energy surface allowed us to monitor the full range of molecular transformation of amyloidogenesis. We identified points of no return (e.g.; T=37 °C, pH 4.1, c_E5=250 μm , c_NaCl=50 mm , t>4–6 h) that will inevitably gravitate into the amyloid state. The strong B-type far ultraviolet (FUV)-ECD spectra and an unexpectedly strong near ultraviolet (NUV)-ECD signal (Θ_{≈275–285 nm}) indicate that the amyloid phase of E5 is built from monomers of quasi-elongated backbone structure (φ≈−145°, ψ≈+145°) with strong interstrand Tyr↔Trp interaction. Misfolded intermediates and the buildup of “toxic” early-stage oligomers leading to self-association were identified and monitored as a function of time. Results indicate that the amyloid transition is triggered by subtle misfolding of the α-helix, exposing aromatic and hydrophobic side chains that may provide the first centers for an intermolecular reorganization. These initial clusters provide the spatial closeness and sufficient time for a transition to the β-structured amyloid nucleus, thus the process follows a nucleated growth mechanism.     

An Off‐Pathway Folding Intermediate of an Acyl Carrier Protein Domain Coexists with the Folded and Unfolded States under Native Conditions

A protein can exist in multiple states under native conditions and those states with low populations are often critical to biological function and self‐assembly. To investigate the role of the minor states of an acyl carrier protein, NMR techniques were applied to determine the number of minor states and characterize their structures and kinetics. The acyl carrier protein from Micromonospora echinospora was found to exist in one major folded state (95.2 %), one unfolded state (4.1 %), and one intermediate state (0.7 %) under native conditions. The three states are in dynamic equilibrium and the intermediate state very likely adopts a native‐like structure and is an off‐pathway folding product. The intermediate state may mediate the formation of oligomers in vitro and play an important role in the recognition of partner enzymes in vivo.     

Reduction of the Search Space for the Folding of Proteins at the Level of Formation and Assembly of Secondary Structures: A New View on the Solution of Levinthal′s Paradox

The complete volume of the protein conformation space is, by many orders of magnitude, smaller at the level of secondary structure elements than that at the level of amino acid residues; the latter, according to Levinthal′s estimate, scales approximately as 10^{2 L}, with L being the number of residues in the chain, whereas the former, as demonstrated in this paper, scales no faster than ～L^N, with N being the number of the secondary structure elements, which is approximately equal to L/15. This drastic decrease in the exponent (L/15 instead of 2 L) explains why sampling of the conformation space does not contradict the ability of the protein chain to find its most stable fold.     

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.     

Common folding processes of mini‐proteins: Partial formations of secondary structures initiate the immediate protein folding

The folding processes of mini‐proteins (FSD‐EY/FBPWW28 domain) were computationally investigated by an enhanced conformational sampling method. Through the analyses of trajectories, these mini‐proteins had multiple folding pathways different from a simple two‐state folding, and the multiple folding processes were initiated by partial formations of secondary structures. These findings can be used to understand the folding of large proteins, that is, which secondary structures are partially folded in the early process, and how the remaining parts are sequentially folded. It is found that FSD‐EY (α/β topology) folds by a simple diffusion‐collision mechanism, while the folding process of the FBPWW28 domain (all‐β topology) requires a modification of the diffusion‐collision theory to adequately treat the coil‐sheet transition for the β sheet formation. © 2017 Wiley Periodicals, Inc.     

Dynamic folding pathway models of the villin headpiece subdomain (HP‐36) structure

We have investigated the folding pathway of the 36‐residue villin headpiece subdomain (HP‐36) by action‐derived molecular dynamics simulations. The folding is initiated by hydrophobic collapse, after which the concurrent formation of full tertiary structure and α‐helical secondary structure is observed. The collapse is observed to be associated with a couple of specific native contacts contrary to the conventional nonspecific hydrophobic collapse model. Stable secondary structure formation after the collapse suggests that the folding of HP‐36 follows neither the framework model nor the diffusion‐collision model. The C‐terminal helix forms first, followed by the N‐terminal helix positioned in its native orientation. The short middle helix is shown to form last. Signs for multiple folding pathways are also observed. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Solvent viscosity dependence of the folding rate of a small protein: distributed computing study

By using distributed computing techniques and a supercluster of more than 20,000 processors we simulated folding of a 20-residue Trp Cage miniprotein in atomistic detail with implicit GB/SA solvent at a variety of solvent viscosities (gamma). This allowed us to analyze the dependence of folding rates on viscosity. In particular, we focused on the low-viscosity regime (values below the viscosity of water). In accordance with Kramers' theory, we observe approximately linear dependence of the folding rate on 1/gamma for values from 1-10(-1)x that of water viscosity. However, for the regime between 10(-4)-10(-1)x that of water viscosity we observe power-law dependence of the form k approximately gamma(-1/5). These results suggest that estimating folding rates from molecular simulations run at low viscosity under the assumption of linear dependence of rate on inverse viscosity may lead to erroneous results.     

Relationships between unfolded configurations of proteins and dynamics of folding to the native state

We compare folding trajectories of chymotrypsin inhibitor (CI2) using a dynamic Monte Carlo scheme with Go-type potentials. The model considers the four backbone atoms of each residue and a sphere centered around C^β the diameter of which is chosen according to the type of the side group. Bond lengths and bond angles are kept fixed. Folding trajectories are obtained by giving random increments to the φ and ψ torsion angles with some bias toward the native state. Excluded volume effects are considered. Two sets of 20 trajectories are obtained, with different initial configurations. The first set is generated from random initial configurations. The initial configurations of the second set are generated according to knowledge-based neighbor dependent torsion probabilities derived from triplets in the Protein Data Bank. Compared to chains with randomly generated initial configurations, those generated with neighbor-dependent probabilities (i) fold faster, (ii) have better defined secondary structure elements, and (iii) have less number of non-native contacts during folding. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3667–3678, 2006     

Dual Binding Mode of Metallacarborane Produces a Robust Shield on Proteins

An inorganic sandwich molecule, Na[Co(C₂B₉H₁₁)₂], able to produce vesicles through self-assembly and known to produce strong dihydrogen-bond interactions with amine groups is capable of interacting with proteins. This dual non-bonding ability of Na[Co(C₂B₉H₁₁)₂] is what makes this molecule unique: it can be firmly anchored to a protein surface and is capable of extending over it. To prove this, the widely available bovine serum albumin (BSA), which has many pendant amino groups in its structure, has been taken as the model protein. It has been found that around 100 molecules of Na[Co(C₂B₉H₁₁)₂] preserve the native structure of BSA, while endorsing it with a significantly increased stability with respect to chemical- and thermal-induced denaturation due to efficient encapsulation. The advantages of this encapsulation technique are two-fold; the first is its simplicity as it relies on the anchoring capacity of Na[Co(C₂B₉H₁₁)₂] to the surface of the protein through the amine-containing residues and the second is its self-assembling capacity allowing it to spread across the surface. The dense shield of protection offered by Na[Co(C₂B₉H₁₁)₂] has been demonstrated by the inhibition of BSA pseudo-esterase activity, which indicates that the inorganic corset around BSA protects its reactive surface residues, thereby preventing their acetylation.     

Flexible Folding: Disulfide-Containing Peptides and Proteins Choose the Pathway Depending on the Environments

In the last few decades, development of novel experimental techniques, such as new types of disulfide (SS)-forming reagents and genetic and chemical technologies for synthesizing designed artificial proteins, is opening a new realm of the oxidative folding study where peptides and proteins can be folded under physiologically more relevant conditions. In this review, after a brief overview of the historical and physicochemical background of oxidative protein folding study, recently revealed folding pathways of several representative peptides and proteins are summarized, including those having two, three, or four SS bonds in the native state, as well as those with odd Cys residues or consisting of two peptide chains. Comparison of the updated pathways with those reported in the early years has revealed the flexible nature of the protein folding pathways. The significantly different pathways characterized for hen-egg white lysozyme and bovine milk α-lactalbumin, which belong to the same protein superfamily, suggest that the information of protein folding pathways, not only the native folded structure, is encoded in the amino acid sequence. The application of the flexible pathways of peptides and proteins to the engineering of folded three-dimensional structures is an interesting and important issue in the new realm of the current oxidative protein folding study.     

Frustration Sculpts the Early Stages of Protein Folding

The funneled energy landscape theory implies that protein structures are minimally frustrated. Yet, because of the divergent demands between folding and function, regions of frustrated patterns are present at the active site of proteins. To understand the effects of such local frustration in dictating the energy landscape of proteins, here we compare the folding mechanisms of the two alternative spliced forms of a PDZ domain (PDZ2 and PDZ2as) that share a nearly identical sequence and structure, while displaying different frustration patterns. The analysis, based on the kinetic characterization of a large number of site‐directed mutants, reveals that although the late stages for folding are very robust and biased by native topology, the early stages are more malleable and dominated by local frustration. The results are briefly discussed in the context of the energy‐landscape theory.     

Comparison of stochastic optimization methods for all-atom folding of the Trp-Cage protein.

The performances of three different stochastic optimization methods for all-atom protein structure prediction are investigated and compared. We use the recently developed all-atom free-energy force field (PFF01), which was demonstrated to correctly predict the native conformation of several proteins as the global optimum of the free energy surface. The trp-cage protein (PDB-code 1L2Y) is folded with the stochastic tunneling method, a modified parallel tempering method, and the basin-hopping technique. All the methods correctly identify the native conformation, and their relative efficiency is discussed.     

Randomizing the Unfolded State of Peptides (and Proteins) by Nearest Neighbor Interactions between Unlike Residues

To explore the influence of nearest neighbors on conformational biases in unfolded peptides, we combined vibrational and 2D NMR spectroscopy to obtain the conformational distributions of selected “GxyG” host–guest peptides in aqueous solution: GDyG, GSyG, GxLG, GxVG, where x/y=A, K, L, V. Large changes of conformational propensities were observed due to nearest‐neighbor interactions, at variance with the isolated pair hypothesis. We found that protonated aspartic acid and serine lose their above‐the‐average preference for turn‐like structures in favor of polyproline II (pPII) populations in the presence of neighbors with bulky side chains. Such residues also decrease the above‐the‐average pPII preference of alanine. These observations suggest that the underlying mechanism involves a disruption of the hydration shell. Thermodynamic analysis of ³J(H^N,H^α) (T) data for each x,y residue reveals that modest changes in the conformational ensemble masks larger changes of enthalpy and entropy governing the pPII?β equilibrium indicating a significant residue dependent temperature dependence of the peptides’ conformational ensembles. These results suggest that nearest‐neighbor interactions between unlike residues act as conformational randomizers close to the enthalpy–entropy compensation temperature, eliminating intrinsic biases in favor of largely balanced pPII/β dominated ensembles at physiological temperatures.     

Amyloid Polymorphism in the Protein Folding and Aggregation Energy Landscape

Protein folding involves a large number of steps and conformations in which the folding protein samples different thermodynamic states characterized by local minima. Kinetically trapped on‐ or off‐pathway intermediates are metastable folding intermediates towards the lowest absolute energy minima, which have been postulated to be the natively folded state where intramolecular interactions dominate, and the amyloid state where intermolecular interactions dominate. However, this view largely neglects the rich polymorphism found within amyloid species. We review the protein folding energy landscape in view of recent findings identifying specific transition routes among different amyloid polymorphs. Observed transitions such as twisted ribbon→crystal or helical ribbon→nanotube, and forbidden transitions such helical ribbon?crystal, are discussed and positioned within the protein folding and aggregation energy landscape. Finally, amyloid crystals are identified as the ground state of the protein folding and aggregation energy landscape.     

Effects of Ionic Liquids on Metalloproteins

In the past decade, innovative protein therapies and bio-similar industries have grown rapidly. Additionally, ionic liquids (ILs) have been an area of great interest and rapid development in industrial processes over a similar timeline. Therefore, there is a pressing need to understand the structure and function of proteins in novel environments with ILs. Understanding the short-term and long-term stability of protein molecules in IL formulations will be key to using ILs for protein technologies. Similarly, ILs have been investigated as part of therapeutic delivery systems and implicated in numerous studies in which ILs impact the activity and/or stability of protein molecules. Notably, many of the proteins used in industrial applications are involved in redox chemistry, and thus often contain metal ions or metal-associated cofactors. In this review article, we focus on the current understanding of protein structure-function relationship in the presence of ILs, specifically focusing on the effect of ILs on metal containing proteins.