Hydrophobic cooperativity as a mechanism for amyloid nucleation

Marine-Derived Bioactive Peptides Self-Assembled Multifunctional Materials: Antioxidant and Wound Healing

Peptide self-assembling materials have received significant attention from researchers in recent years, emerging as a popular field in biological, environmental, medical, and other new materials studies. In this study, we utilized controllable enzymatic hydrolysis technology (animal proteases) to obtain supramolecular peptide self-assembling materials (CAPs) from the Pacific oyster (Crassostrea gigas). We conducted physicochemical analyses to explore the pro-healing mechanisms of CAPs on skin wounds in both in vitro and in vivo experiments through a topical application. The results demonstrated that CAPs exhibit a pH-responsive behavior for self-assembly and consist of peptides ranging from 550 to 2300 Da in molecular weight, with peptide chain lengths of mainly 11-16 amino acids. In vitro experiments indicated that CAPs display a procoagulant effect, free radical scavenging activity, and promote the proliferation of HaCaTs (112.74% and 127.61%). Moreover, our in vivo experiments demonstrated that CAPs possess the ability to mitigate inflammation, boost fibroblast proliferation, and promote revascularization, which accelerates the epithelialization process. Consequently, a balanced collagen I/III ratio in the repaired tissue and the promotion of hair follicle regeneration were observed. With these remarkable findings, CAPs can be regarded as a natural and secure treatment option with high efficacy for skin wound healing. The potential of CAPs to be further developed for traceless skin wound healing is an exciting area for future research and development.

Molecular Modification of Kex2 P1' Site Enhances Expression and Druggability of Fungal Defensin

Pichia pastoris is the widely used expression system for producing recombinant secretory proteins. It is known that Kex2 protease plays a vital role in the process of protein secretion, in which the P1' site affects its cleavage efficiency. To enhance the expression level of fungal defensin-derived peptide NZ2114, this work attempts to optimize the P1' site of Kex2 by replacing it with 20 amino acids in turn. The results showed that when the amino acid of the P1' site was changed to Phe (F), the yield of target peptide significantly increased from 2.39 g/L to 4.81 g/L. Additionally, the novel peptide F-NZ2114 (short for FNZ) showed strong antimicrobial activity against Gram-positive (G⁺) bacteria, especially for Staphylococcus aureus and Streptococcus agalactiae (MIC: 4-8 μg/mL). The FNZ was very stable and retained high activity in various conditions; in addition, a low cytotoxicity and no hemolysis were observed even at a high concentration of 128 μg/mL, and a longer postantibiotic effect was reached. The above results indicate that this engineering strategy provided a feasible optimization scheme for enhancing the expression level and druggability of this antimicrobial peptide from fungal defensin and other similar targets by this updated recombinant yeast.

Microscopic Understanding of the Conformational Stability of the Aggregated Nonamyloid β Components of α-Synuclein

Self-association of α-synuclein peptides into oligomeric species and ordered amyloid fibrils is associated with Parkinson's disease, a progressive neurodegenerative disorder. In particular, the peptide domain formed between the residues Glu-61 (or E61) and Val-95 (or V95) of α-synuclein, typically termed the "nonamyloid β component" (NAC), is known to play critical roles in forming aggregated structures. In this work, we have employed molecular dynamics simulations to explore the conformational properties and relative stabilities of aggregated protofilaments of different orders, namely, tetramer (P(4)), hexamer (P(6)), octamer (P(8)), decamer (P(10)), dodecamer (P(12)), and tetradecamer (P(14)), formed by the NAC domains of α-synuclein. Besides, center-of-mass pulling and umbrella sampling simulation methods have also been employed to characterize the mechanistic pathway of peptide association/dissociation and the corresponding free energy profiles. Structural analysis showed that the disordered C-terminal loop and the central core regions of the peptide units lead to more flexible and distorted structures of the lower order protofilaments (P(4) and P(6)) as compared to the higher order ones. Interestingly, our calculation shows the presence of multiple distinctly populated conformational states for the lower order protofilament P(4), which may drive the oligomerization process along multiple pathways to form different polymorphic α-synuclein fibrillar structures. It is further observed that the nonpolar interaction between the peptides and the corresponding nonpolar solvation free energy play a dominant role in stabilizing the aggregated protofilaments. Importantly, our result showed that reduced cooperativity during the binding of a peptide unit beyond a critical size of the protofilament (P(12)) leads to less favorable binding free energy of a peptide.

Modulating Nanodroplet Formation En Route to Fibrillization of Amyloid Peptides with Designed Flanking Sequences

Soluble oligomers populating early amyloid aggregation can be regarded as nanodroplets of liquid-liquid phase separation (LLPS). Amyloid peptides typically contain hydrophobic aggregation-prone regions connected by hydrophilic linkers and flanking sequences, and such a sequence hydropathy pattern drives the formation of supramolecular structures in the nanodroplets and modulates subsequent fibrillization. Here, we studied LLPS and fibrillization of coarse-grained amyloid peptides with increasing flanking sequences. Nanodroplets assumed lamellar, cylindrical micellar, and spherical micellar structures with increasing peptide hydrophilic/hydrophobic ratios, and such morphologies governed subsequent fibrillization processes. Adding glycine-serine repeats as flanking sequences to Aβ16-22, the amyloidogenic core of amyloid-β, our computational predictions of morphological transitions were corroborated experimentally. The uncovered inter-relationships between the peptide sequence pattern, oligomer/nanodroplet morphology, and fibrillization pathway, kinetics, and structure may contribute to our understanding of pathogenic amyloidosis in aging, facilitate future efforts ameliorating amyloidosis through peptide engineering, and aid in the design of novel amyloid-based functional nanobiomaterials and nanocomposites.

Scan-Find-Scan-Model: Discrete Site-Targeted Suppressor Design Strategy for Amyloid-β

Alzheimer's disease is undoubtedly the most well-studied neurodegenerative disease. Consequently, the amyloid-β (Aβ) protein ranks at the top in terms of getting attention from the scientific community for structural property-based characterization. Even after decades of extensive research, there is existing volatility in terms of understanding and hence the effective tackling procedures against the disease that arises due to the lack of knowledge of both specific target- and site-specific drugs. Here, we develop a multidimensional approach based on the characterization of the common static-dynamic-thermodynamic trait of the monomeric protein, which efficiently identifies a small target sequence that contains an inherent tendency to misfold and consequently aggregate. The robustness of the identification of the target sequence comes with an abundance of a priori knowledge about the length and sequence of the target and hence guides toward effective designing of the target-specific drug with a very low probability of bottleneck and failure. Based on the target sequence information, we further identified a specific mutant that showed the maximum potential to act as a destabilizer of the monomeric protein as well as enormous success as an aggregation suppressor. We eventually tested the drug efficacy by estimating the extent of modulation of binding affinity existing within the fibrillar form of the Aβ protein due to a single-point mutation and hence provided a proof of concept of the entire protocol.

A coarse-grained mechanical model for folding and unfolding of tropoelastin with possible mutations

We propose a simple general framework to predict folding, native states, energy barriers, protein unfolding, as well as mutation induced diseases and other protein structural analyses. The model should not be considered as an alternative to classical approaches (Molecular Dynamics or Monte Carlo) because it neglects low scale details and rather focuses on global features of proteins and structural information. We aim at the description of phenomena that are out of the range of classical molecular modeling approaches due to the large computational cost: multimolecular interactions, cyclic behavior under variable external interactions, and similar. To demonstrate the effectiveness of the approach in a real case, we focus on the folding and unfolding behavior of tropoelastin and its mutations. Specifically, we derive a discrete mechanical model whose structure is deduced based on a coarse graining approach that allows us to group the amino acids sequence in a smaller number of `equivalent' masses. Nearest neighbor energy terms are then introduced to reproduce the interaction of such amino acid groups. Nearest and non-nearest neighbor energy terms, inter and intra functional blocks are phenomenologically added in the form of Morse potentials. As we show, the resulting system reproduces important properties of the folding-unfolding mechanical response, including the monotonic and cyclic force-elongation behavior, representing a physiologically important information for elastin. The comparison with the experimental behavior of mutated tropoelastin confirms the predictivity of the model. STATEMENT OF SIGNIFICANCE: Classical approaches to the study of phenomena at the molecular scale such as Molecular Dynamics (MD) represent an incredible tool to unveil mechanical and conformational properties of macromolecules, in particular for biological and medical applications. On the other hand, due to the computational cost, the time and spatial scales are limited. Focusing of the real case of tropoelastin, we propose a new approach based on a careful coarse graining of the system, able to describe the overall properties of the macromolecule and amenable of extension to larger scale effects (protein bundles, protein-protein interactions, cyclic loading). The comparison with tropoelastin behavior, also for mutations, is very promising.

A theoretical study of polymorphism in VQIVYK fibrils

The VQIVYK fragment from the Tau protein, also known as PHF6, is essential for aggregation of Tau into neurofibrillary lesions associated with neurodegenerative diseases. VQIVYK itself forms amyloid fibrils composed of paired β-sheets. Therefore, the full Tau protein and VQIVYK fibrils have been intensively investigated. A central issue in these studies is polymorphism, the ability of a protein to fold into more than one structure. Using all-atom molecular simulations, we generate five stable polymorphs of VQIVYK fibrils, establish their relative free energy with umbrella sampling methods, and identify the side chain interactions that provide stability. The two most stable polymorphs, which have nearly equal free energy, are formed by interdigitation of the mostly hydrophobic VIY "face" sides of the β-sheets. Another stable polymorph is formed by interdigitation of the QVK "back" sides. When we turn to examine structures from cryo-electron microscopy experiments on Tau filaments taken from diseased patients or generated in vitro, we find that the pattern of side chain interactions found in the two most stable face-to-face as well as the back-to-back polymorphs are recapitulated in amyloid structures of the full protein. Thus, our studies suggest that the interactions stabilizing PHF6 fibrils explain the amyloidogenicity of the VQIVYK motif within the full Tau protein and provide justification for the use of VQIVYK fibrils as a test bed for the design of molecules that identify or inhibit amyloid structures.

Dual effect of 2D WS2 nanoparticles on the lysozyme conformation

In the present work we studied the effect of 2D WS₂ nanoparticles on the conformational changes in lysozyme protein at different pH values (2.0-11.5). The contributions of various structural conformations (α-helix, β-sheets parallel and antiparallel, unordered structure and side groups) were determined by decomposition of Amid I absorbance bands. The 2D WS₂ were shown to have different impact on secondary structure depending on pH of the solution and protein concentration. The amyloid fibril presence was confirmed with confocal microscopy enhanced by gold support, and fluorescent spectroscopy with amyloid-sensitive dye Thioflavin T. Our data show that WS₂ can both inhibit and stimulate amyloid formation. Additionally, we have also reported an unusual spectroscopic behavior displayed by lysozyme, indicated by narrowing of Amide I and Amide II bands at pH 2.5 and 3.5 when incubated with 2D WS2 nanoparticles.

Protein Supramolecular Structures: From Self-Assembly to Nanovaccine Design

Life-inspired protein supramolecular assemblies have recently attracted considerable attention for the development of next-generation vaccines to fight against infectious diseases, as well as autoimmune diseases and cancer. Protein self-assembly enables atomic scale precision over the final architecture, with a remarkable diversity of structures and functionalities. Self-assembling protein nanovaccines are associated with numerous advantages, including biocompatibility, stability, molecular specificity and multivalency. Owing to their nanoscale size, proteinaceous nature, symmetrical organization and repetitive antigen display, protein assemblies closely mimic most invading pathogens, serving as danger signals for the immune system. Elucidating how the structural and physicochemical properties of the assemblies modulate the potency and the polarization of the immune responses is critical for bottom-up design of vaccines. In this context, this review briefly covers the fundamentals of supramolecular interactions involved in protein self-assembly and presents the strategies to design and functionalize these assemblies. Examples of advanced nanovaccines are presented, and properties of protein supramolecular structures enabling modulation of the immune responses are discussed. Combining the understanding of the self-assembly process at the molecular level with knowledge regarding the activation of the innate and adaptive immune responses will support the design of safe and effective nanovaccines.

Critical nucleus of Greek-key-like core of α-synuclein protofibril and its disruption by dopamine and norepinephrine

Protofibrillar trimer is the critical nucleus for the αS fibril formation, and the tetramer is the minimal stable nucleus. The interactions of DA/NE with αS trimer/tetramer disrupt the backbone H-bonds and destabilize the αS protofibrils.

A computational approach to get insights into multiple faces of additives in modulation of protein aggregation pathways

An enormous population worldwide is presently confronted with debilitating neurodegenerative diseases. The etiology of the disease is connected to protein aggregation and the events involved therein. Thus, a complete understanding of an inhibitor at different stages in the process is imperative for the formulation of a drug molecule. This review presents a detailed summary of the current status of different cosolvents. It further develops how the complex aggregation pathway can be simplified into three steps common to all proteins and the way computer simulations can be exploited to gain insights into the ways by which known inhibitors can affect all these stages. Computation of theoretical parameters in this regard and their correlation with experimental techniques is accentuated. In addition to providing an outline of the scope of different additives, this review showcases the way by which the problem of analyzing an effect of an additive can be addressed effectively via MD simulations.

Adsorption of Amyloidogenic Peptides to Functionalized Surfaces Is Biased by Charge and Hydrophilicity

Surfaces are abundant in living systems, such as in the form of cellular membranes, and govern many biological processes. In this study, the adsorption of the amyloidogenic model peptides GNNQQNY, NNFGAIL, and VQIVYK as well as the amyloid-forming antimicrobial peptide uperin 3.5 (U3.5) were studied at low concentrations (100 μM) to different surfaces. The technique of a quartz crystal microbalance with dissipation monitoring (QCM-D) was applied as it enables the monitoring of mass binding to sensors at nanogram sensitivity. Gold-coated quartz sensors were used as unmodified gold surfaces or functionalized with self-assembled monolayers (SAMs) of alkanethiols (terminated as methyl, amino, carboxyl, and hydroxyl) resulting in different adsorption affinities of the peptides. Our objective was to evaluate the underlying role of the nature and feature of interfaces in biological systems which could concentrate peptides and impact or trigger peptide aggregation processes. In overall, the largely hydrophobic peptides adsorbed with preference to hydrophobic or countercharged surfaces. Further, the glycoprotein lubricin (LUB) was tested as an antiadhesive coating. Despite its hydrophilicity, the adsorption of peptides to LUB coated sensors was similar to the adsorption to unmodified gold surfaces, which indicates that some peptides diffused through the LUB layer to reach the underlying gold sensor surface. The LUB protein-antiadhesive is thus more effective as a biomaterial coating against larger biomolecules than small peptides under the conditions used here. This study provides directions toward a better understanding of amyloid peptide adsorption to biologically relevant interfaces, such as cellular membranes.

The initial stage of structural transformation of Aβ42 peptides from the human and mole rat in the presence of Fe2+ and Fe3+: Related to Alzheimer's disease

The early stage of secondary structural conversion of amyloid beta (Aβ) to misfolded aggregations is a key feature of Alzheimer's disease (AD). Under normal physiological conditions, Aβ peptides can protect neurons from the toxicity of highly concentrated metals. However, they become toxic under certain conditions. Under conditions of excess iron, amyloid precursor proteins (APP) become overexpressed. This subsequently increases Aβ production. Experimental studies suggest that Aβ fibrillation (main-pathway) and amorphous (off-pathway) aggregate formations are two competitive pathways driven by factors such as metal binding, pH and temperature. In this study, we performed molecular dynamic (MD) simulations to examine the initial stage of conformational transformations of human Aβ (hAβ) and rat Aβ (rAβ) peptides in the presence of Fe²⁺ and Fe³⁺ ions. Our results demonstrated that Fe²⁺ and Fe³⁺ play key roles in Aβs folding and aggregation. Fe³⁺ had a greater effect than Fe²⁺on Aβs' folding during intermolecular interactions and subsequently, had a greater effect in decreasing structural diversity. Fe²⁺ was observed to be more likely than Fe³⁺ to interact with nitrogen atoms from the residues of imidazole rings of His. rAβ peptides are more energetically favorable than hAβ for intermolecular interactions and amorphous aggregations. We concluded that most hAβ structures were energetically unfavorable. However, hAβs with intermolecular β-sheet formations in the C-terminal were energetically favorable. It is notable that Fe²⁺ can change the surface charge of hAβ. Furthermore, Fe³⁺ can promote C-terminal folding by binding to Glu22 and Ala42, and by forming stable β-sheet formations on the C-terminal. Fe³⁺ can also pause the main-pathway by inducing random aggregations.

Molecular Recognition of Structures Is Key in the Polymerization of Patterned Barnacle Adhesive Sequences

The permanent adhesive produced by adult barnacles is held together by tightly folded proteins that form amyloid-like materials distinct among marine foulants. In this work, we link stretches of alternating charged and noncharged linear sequences from a family of adhesive proteins to their role in forming fibrillar nanomaterials. Using recombinant proteins and short barnacle cement derived peptides (BCPs), we find a central sequence with charged motifs of the pattern [Gly/Ser/Val/Thr/Ala-X], where X are charged amino acids, to exert specific control over timing, structure, and morphology of fibril formation. While most BCPs remain dormant, the core segment demonstrates rapid polymerization as well as an ability to template other peptides with no propensity for self-assembly. Patterned charge domains assemble dormant peptides through a specific antiparallel β-sheet structure as measured by FTIR. While charged domains favor an antiparallel structure, BCPs without charged domains switch fibril assembly to favor simpler parallel β-sheet aggregates. In addition to activation, charged domains direct nanofibers to grow into discrete microns long fibrils similar to the natural adhesive, while segments without such domains only form short branched aggregates. The assembly of adhesive sequences through recognition of structured templates outlines a strategy used by barnacles to control physical mechanisms of underwater adhesive delivery, activation, and curing based on molecular recognition between proteins.

Influence of crowding and surfaces on protein amyloidogenesis: A thermo-kinetic perspective

The last few decades have irreversibly implicated protein self-assembly and aggregation leading to amyloid fibril formation in proteopathies that include several neurodegenerative diseases. Emerging studies recognize the importance of eliciting the pathways leading to protein aggregation in the context of the crowded intracellular environment rather than in conventional in vitro conditions. It is found that crowded environments can have acceleratory as well as inhibitory effects on protein aggregation, depending on the interplay of underlying factors on the crucial rate limiting steps. The aggregation mechanism and transient species formed along the pathway are further altered when they interface with natural and artificial surfaces in the cellular milieu. An increasing number of studies probe the autocatalytic nature of amyloid surfaces as well as membrane bilayer effects on amyloidogenesis. Moreover, exposure to modern nanosurfaces via nanomedicines and other sources potentially invokes beneficial or deleterious biological response that needs rigorous investigation. Mounting evidences indicate that nanoparticles can either promote or impede amyloid aggregation, spurring efforts to tune their interactions for developing effective anti-amyloid strategies. Mechanistic insights into nanoparticle mediated aggregation pathways are therefore crucial for engineering anti-amyloid nanoparticle strategies that are biocompatible and sustainable. This review is a compilation of studies that contribute to the current understanding of the altering effects of molecular crowding as well as natural and artificial surfaces on protein amyloidogenesis.

Biomolecular Assemblies: Moving from Observation to Predictive Design

Biomolecular assembly is a key driving force in nearly all life processes, providing structure, information storage, and communication within cells and at the whole organism level. These assembly processes rely on precise interactions between functional groups on nucleic acids, proteins, carbohydrates, and small molecules, and can be fine-tuned to span a range of time, length, and complexity scales. Recognizing the power of these motifs, researchers have sought to emulate and engineer biomolecular assemblies in the laboratory, with goals ranging from modulating cellular function to the creation of new polymeric materials. In most cases, engineering efforts are inspired or informed by understanding the structure and properties of naturally occurring assemblies, which has in turn fueled the development of predictive models that enable computational design of novel assemblies. This Review will focus on selected examples of protein assemblies, highlighting the story arc from initial discovery of an assembly, through initial engineering attempts, toward the ultimate goal of predictive design. The aim of this Review is to highlight areas where significant progress has been made, as well as to outline remaining challenges, as solving these challenges will be the key that unlocks the full power of biomolecules for advances in technology and medicine.

Insights into Structure-Activity Relationships of 3-Arylhydrazonoindolin-2-One Derivatives for Their Multitarget Activity on β-Amyloid Aggregation and Neurotoxicity

Despite the controversial outcomes of clinical trials executed so far, the prevention of β-amyloid (Aβ) deposition and neurotoxicity by small molecule inhibitors of Aβ aggregation remains a target intensively pursued in the search of effective drugs for treating Alzheimer’s disease (AD) and related neurodegeneration syndromes. As a continuation of previous studies, a series of new 3-(2-arylhydrazono)indolin-2-one derivatives was synthesized and assayed, investigating the effects of substitutions on both the indole core and arylhydrazone moiety. Compared with the reference compound 1, we disclosed equipotent derivatives bearing alkyl substituents at the indole nitrogen, and fairly tolerated bioisosteric replacements at the arylhydrazone moiety. For most of the investigated compounds, the inhibition of Aβ₄₀ aggregation (expressed as pIC₅₀) was found to be correlated with lipophilicity, as assessed by a reversed-phase HPLC method, through a bilinear relationship. The N¹-cyclopropyl derivative 28 was tested in cell-based assays of Aβ₄₂ oligomer toxicity and oxidative stress induced by hydrogen peroxide, showing significant cytoprotective effects. This study confirmed the versatility of isatin in preparing multitarget small molecules affecting different biochemical pathways involved in AD.

Sequence dependent aggregation of peptides and fibril formation

Deciphering the links between amino acid sequence and amyloid fibril formation is key for understanding protein misfolding diseases. Here we use Monte Carlo simulations to study the aggregation of short peptides in a coarse-grained model with hydrophobic-polar (HP) amino acid sequences and correlated side chain orientations for hydrophobic contacts. A significant heterogeneity is observed in the aggregate structures and in the thermodynamics of aggregation for systems of different HP sequences and different numbers of peptides. Fibril-like ordered aggregates are found for several sequences that contain the common HPH pattern, while other sequences may form helix bundles or disordered aggregates. A wide variation of the aggregation transition temperatures among sequences, even among those of the same hydrophobic fraction, indicates that not all sequences undergo aggregation at a presumable physiological temperature. The transition is found to be the most cooperative for sequences forming fibril-like structures. For a fibril-prone sequence, it is shown that fibril formation follows the nucleation and growth mechanism. Interestingly, a binary mixture of peptides of an aggregation-prone and a non-aggregation-prone sequence shows the association and conversion of the latter to the fibrillar structure. Our study highlights the role of a sequence in selecting fibril-like aggregates and also the impact of a structural template on fibril formation by peptides of unrelated sequences.

Mechanisms and rates of nucleation of amyloid fibrils

The classical nucleation theory finds the rate of nucleation proportional to the monomer concentration raised to the power, which is the "critical nucleus size," n_c. The implicit assumption, that amyloids nucleate in the same way, has been recently challenged by an alternative two-step mechanism, when the soluble monomers first form a metastable aggregate (micelle) and then undergo conversion into the conformation rich in β-strands that are able to form a stable growing nucleus for the protofilament. Here we put together the elements of extensive knowledge about aggregation and nucleation kinetics, using a specific case of Aβ_1-42 amyloidogenic peptide for illustration, to find theoretical expressions for the effective rate of amyloid nucleation. We find that at low monomer concentrations in solution and also at low interaction energy between two peptide conformations in the micelle, the nucleation occurs via the classical route. At higher monomer concentrations, and a range of other interaction parameters between peptides, the two-step "aggregation-conversion" mechanism of nucleation takes over. In this regime, the effective rate of the process can be interpreted as a power of monomer concentration in a certain range of parameters; however, the exponent is determined by a complicated interplay of interaction parameters and is not related to the minimum size of the growing nucleus (which we find to be ∼7-8 for Aβ_1-42).

Peptide strand length controls the energetics of self-assembly and morphology of β-sheet fibrils

Self-assembling peptides can be used as versatile, natural, and multifunctional building blocks to produce a variety of well-defined nanostructures, materials and devices for applications in medicine and nanotechnology. Here, we concentrate on the 1D self-assembly of de novo designed Px-2 peptide β-strands into anti-parallel β-sheet tapes and higher order aggregates. We study six members of the Px-2 family, ranging from 3 amino acids (aa) to 13 aa in length, using a range of complementary experimental techniques, computer simulation and theoretical statistical mechanics. The critical concentration for self-assembly (c*) is found to increase systematically with decreasing peptide length. The shortest peptide found to self-assemble into soluble β-tapes in water is a 5 amino acid residue peptide. These investigations help decipher the role of the peptide length in controlling self-assembly, aggregate morphology, and material properties. By extracting free energies from these data using a statistical mechanical analysis and combining the results with computer simulations at the atomistic level, we can extract the entropy of association for individual β-strands.

Mannich base approach to 5-methoxyisatin 3-(4-isopropylphenyl)hydrazone: A water-soluble prodrug for a multitarget inhibition of cholinesterases, beta-amyloid fibrillization and oligomer-induced cytotoxicity

Targeting protein aggregation for the therapy of neurodegenerative diseases remains elusive for medicinal chemists, despite a number of small molecules known to interfere in amyloidogenesis, particularly of amyloid beta (Aβ) protein. Starting from previous findings in the antiaggregating activity of a class of indolin-2-ones inhibiting Aβ fibrillization, 5-methoxyisatin 3-(4-isopropylphenyl)hydrazone 1 was identified as a multitarget inhibitor of Aβ aggregation and cholinesterases with IC₅₀s in the low μM range. With the aim of increasing aqueous solubility, a Mannich-base functionalization led to the synthesis of N-methylpiperazine derivative 2. At acidic pH, an outstanding solubility increase of 2 over the parent compound 1 was proved through a turbidimetric method. HPLC analysis revealed an improved stability of the Mannich base 2 at pH2 along with a rapid release of 1 in human serum as well as an outstanding hydrolytic stability of the parent hydrazone. Coincubation of Aβ_1-42 with 2 resulted in the accumulation of low MW oligomers, as detected with PICUP assay. Cell assays on SH-SY5Y cells revealed that 2 exerts strong cytoprotective effects in both cell viability and radical quenching assays, mainly related to its active metabolite 1. These findings show that 2 drives the formation of non-toxic, off-pathway Aβ oligomers unable to trigger the amyloid cascade and toxicity.

Lattice model for amyloid peptides: OPEP force field parametrization and applications to the nucleus size of Alzheimer's peptides

Coarse-grained protein lattice models approximate atomistic details and keep the essential interactions. They are, therefore, suitable for capturing generic features of protein folding and amyloid formation at low computational cost. As our aim is to study the critical nucleus sizes of two experimentally well-characterized peptide fragments Aβ16-22 and Aβ37-42 of the full length Aβ1-42 Alzheimer's peptide, it is important that simulations with the lattice model reproduce all-atom simulations. In this study, we present a comprehensive force field parameterization based on the OPEP (Optimized Potential for Efficient protein structure Prediction) force field for an on-lattice protein model, which incorporates explicitly the formation of hydrogen bonds and directions of side-chains. Our bottom-up approach starts with the determination of the best lattice force parameters for the Aβ16-22 dimer by fitting its equilibrium parallel and anti-parallel β-sheet populations to all-atom simulation results. Surprisingly, the calibrated force field is transferable to the trimer of Aβ16-22 and the dimer and trimer of Aβ37-42. Encouraged by this finding, we characterized the free energy landscapes of the two decamers. The dominant structure of the Aβ16-22 decamer matches the microcrystal structure. Pushing the simulations for aggregates between 4-mer and 12-mer suggests a nucleus size for fibril formation of 10 chains. In contrast, the Aβ37-42 decamer is largely disordered with mixed by parallel and antiparallel chains, suggesting that the nucleus size is >10 peptides. Our refined force field coupled to this on-lattice model should provide useful insights into the critical nucleation number associated with neurodegenerative diseases.

The Levinthal Problem in Amyloid Aggregation: Sampling of a Flat Reaction Space

The formation of amyloid fibrils has been associated with many neurodegenerative disorders, yet the mechanism of aggregation remains elusive, partly because aggregation time scales are too long to probe with atomistic simulations. A microscopic theory of fibril elongation was recently developed that could recapitulate experimental results with respect to the effects of temperature, denaturants, and protein concentration on fibril growth kinetics (Schmit, J. D. J. Chem. Phys. 2013, 138 (18), 185102). The theory identifies the conformational search over H-bonding states as the slowest step in the aggregation process and suggests that this search can be efficiently modeled as a random walk on a rugged one-dimensional energy landscape. This insight motivated the multiscale computational algorithm for simulating fibril growth presented in this paper. Briefly, a large number of short atomistic simulations are performed to compute the system diffusion tensor in the reaction coordinate space predicted by the analytic theory. Ensemble aggregation pathways and growth kinetics are then computed from Markov state model (MSM) trajectories. The algorithm is deployed here to understand the fibril growth mechanism and kinetics of Aβ_16-22 and three of its mutants. The order of growth rates of the wild-type and two single mutation peptides (CHA19 and CHA20) predicted by the MSM trajectories is consistent with experimental results. The simulation also correctly predicts that the double mutation (CHA19/CHA20) would reduce the fibril growth rate, even though the degree of rate reduction with respect to either single mutation is overestimated. This artifact may be attributed to the simplistic implicit solvent model. These trends in the growth rate are not apparent from inspection of the rate constants of individual bonds or the lifetimes of the mis-registered states that are the primary kinetic traps but only emerge in the ensemble of trajectories generated by the MSM.

Silicone oil promotes amyloid-like aggregation of αB-crystallin

Silicone oil induced changes of secondary structures of αB-crystallin (CRYAB) and promoted its aggregation into amyloid-like nanostructures.

An Atomistic View of Amyloidogenic Self-assembly: Structure and Dynamics of Heterogeneous Conformational States in the Pre-nucleation Phase

The formation of well-defined filamentous amyloid structures involves a polydisperse collection of oligomeric states for which relatively little is known in terms of structural organization. Here we use extensive, unbiased explicit solvent molecular dynamics (MD) simulations to investigate the structural and dynamical features of oligomeric aggregates formed by a number of highly amyloidogenic peptides at atomistic resolution on the μs time scale. A consensus approach has been adopted to analyse the simulations in multiple force fields, yielding an in-depth characterization of pre-fibrillar oligomers and their global and local structure properties. A collision cross section analysis revealed structurally heterogeneous aggregate ensembles for the individual oligomeric states that lack a single defined quaternary structure during the pre-nucleation phase. To gain insight into the conformational space sampled in early aggregates, we probed their substructure and found emerging β-sheet subunit layers and a multitude of ordered intermolecular β-structure motifs with growing aggregate size. Among those, anti-parallel out-of-register β-strands compatible with toxic β-barrel oligomers were particularly prevalent already in smaller aggregates and formed prior to ordered fibrillar structure elements. Notably, also distinct fibril-like conformations emerged in the oligomeric state and underscore the notion that pre-nucleated oligomers serve as a critical intermediate step on-pathway to fibrils.

Critical Nucleus Structure and Aggregation Mechanism of the C-terminal Fragment of Copper-Zinc Superoxide Dismutase Protein

The aggregation of the copper-zinc superoxide dismutase (SOD1) protein is linked to familial amyotrophic lateral sclerosis, a progressive neurodegenerative disease. A recent experimental study has shown that the (147)GVIGIAQ(153) SOD1 C-terminal segment not only forms amyloid fibrils in isolation but also accelerates the aggregation of full-length SOD1, while substitution of isoleucine at site 149 by proline blocks its fibril formation. Amyloid formation is a nucleation-polymerization process. In this study, we investigated the oligomerization and the nucleus structure of this heptapeptide. By performing extensive replica-exchange molecular dynamics (REMD) simulations and conventional MD simulations, we found that the GVIGIAQ hexamers can adopt highly ordered bilayer β-sheets and β-barrels. In contrast, substitution of I149 by proline significantly reduces the β-sheet probability and results in the disappearance of bilayer β-sheet structures and the increase of disordered hexamers. We identified mixed parallel-antiparallel bilayer β-sheets in both REMD and conventional MD simulations and provided the conformational transition from the experimentally observed parallel bilayer sheets to the mixed parallel-antiparallel bilayer β-sheets. Our simulations suggest that the critical nucleus consists of six peptide chains and two additional peptide chains strongly stabilize this critical nucleus. The stabilized octamer is able to recruit additional random peptides into the β-sheet. Therefore, our simulations provide insights into the critical nucleus formation and the smallest stable nucleus of the (147)GVIGIAQ(153) peptide.

Identification of fibrillogenic regions in human triosephosphate isomerase

Background. Amyloid secondary structure relies on the intermolecular assembly of polypeptide chains through main-chain interaction. According to this, all proteins have the potential to form amyloid structure, nevertheless, in nature only few proteins aggregate into toxic or functional amyloids. Structural characteristics differ greatly among amyloid proteins reported, so it has been difficult to link the fibrillogenic propensity with structural topology. However, there are ubiquitous topologies not represented in the amyloidome that could be considered as amyloid-resistant attributable to structural features, such is the case of TIM barrel topology.

Methods. This work was aimed to study the fibrillogenic propensity of human triosephosphate isomerase (HsTPI) as a model of TIM barrels. In order to do so, aggregation of HsTPI was evaluated under native-like and destabilizing conditions. Fibrillogenic regions were identified by bioinformatics approaches, protein fragmentation and peptide aggregation.

Results. We identified four fibrillogenic regions in the HsTPI corresponding to the β3, β6, β7 y α8 of the TIM barrel. From these, the β3-strand region (residues 59–66) was highly fibrillogenic. In aggregation assays, HsTPI under native-like conditions led to amorphous assemblies while under partially denaturing conditions (urea 3.2 M) formed more structured aggregates. This slightly structured aggregates exhibited residual cross-β structure, as demonstrated by the recognition of the WO1 antibody and ATR-FTIR analysis.

Discussion. Despite the fibrillogenic regions present in HsTPI, the enzyme maintained under native-favoring conditions displayed low fibrillogenic propensity. This amyloid-resistance can be attributed to the three-dimensional arrangement of the protein, where β-strands, susceptible to aggregation, are protected in the core of the molecule. Destabilization of the protein structure may expose inner regions promoting β-aggregation, as well as the formation of hydrophobic disordered aggregates. Being this last pathway kinetically favored over the thermodynamically more stable fibril aggregation pathway.

The Rho Termination Factor of Clostridium botulinum Contains a Prion-Like Domain with a Highly Amyloidogenic Core

Prion-like proteins can switch between a soluble intrinsically disordered conformation and a highly ordered amyloid assembly. This conformational promiscuity is encoded in specific sequence regions, known as prion domains (PrDs). Prions are best known as the causative factors of neurological diseases in mammals. However, bioinformatics analyses reveal that proteins bearing PrDs are present in all kingdoms of life, including bacteria, thus supporting the idea that they serve conserved beneficial cellular functions. Despite the proportion of predicted prion-like proteins in bacterial proteomes is generally low, pathogenic species seem to have a higher prionic load, suggesting that these malleable proteins may favor pathogenic traits. In the present work, we performed a stringent computational analysis of the Clostridium botulinum pathogen proteome in the search for prion-like proteins. A total of 54 candidates were predicted for this anaerobic bacterium, including the transcription termination Rho factor. This RNA-binding protein has been shown to play a crucial role in bacterial adaptation to changing environments. We show here that the predicted disordered PrD domain of this RNA-binding protein contains an inner, highly polar, asparagine-rich short sequence able to spontaneously self-assemble into amyloid-like structures, bearing thus the potential to induce a Rho factor conformational switch that might rewire gene expression in response to environmental conditions.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

Recent advances in the understanding of two-step nucleation of protein crystals

The two-step mechanism of nucleation of crystals in solutions posits that the formation of crystal nuclei occurs within structures of extended lifetimes, in which the nucleating solute is at high concentration. The validity of this mechanism has been demonstrated for proteins, small-molecule organic and inorganic materials, colloids, and polymers. Due to large molecule sizes, proteins are an ideal system to study the details of this nucleation pathway, in particular the formation mechanisms of the nucleation precursors and the associated physico-chemical rules. The precursors of protein crystal nuclei are protein-rich clusters of sizes ∼100 nm that contain 10 000–100 000 molecules and occupy less than 10⁻³of the total solution volume. Here we demonstrate, using oblique illumination microscopy, the liquid nature of the clusters of the protein lysozyme and reveal their inhomogeneous structure. We test a hypothesis put forth by theory that clusters primarily consist of transient protein oligomers. For this, we explore how varying the strength of the Coulomb interaction affects the cluster characteristics. We find that the cluster’s size is insensitive to variations of pH and ionic strength. In contrast, the addition of urea, a chaotropic agent that leads to protein unfolding, strongly decreases the cluster size. Shear stress, a known protein denaturant, induced by bubbling of the solutions with an inert gas, elicits a similar response. These observations support partial protein unfolding, followed by dimerization, as the mechanism of cluster formation. The amide hydrogen–deuterium exchange, monitored by nuclear magnetic resonance, highlights that lysozyme conformational flexibility is a condition for the formation of the protein-rich clusters and facilitates the nucleation of protein crystals.

A carboxylated Zn-phthalocyanine inhibits fibril formation of Alzheimer's amyloid β peptide

Amyloid β (Aβ), a 39-42 amino acid peptide derived from amyloid precursor protein, is deposited as fibrils in Alzheimer's disease brains, and is considered to play a major role in the pathogenesis of the disease. We have investigated the effects of a water-soluble Zn-phthalocyanine, ZnPc(COONa)₈, a macrocyclic compound with near-infrared optical properties, on Aβ fibril formation in vitro. A thioflavin T fluorescence assay showed that ZnPc(COONa)₈ significantly inhibited Aβ fibril formation, increasing the lag time and dose-dependently decreasing the plateau level of fibril formation. Moreover, it destabilized pre-formed Aβ fibrils, resulting in an increase in low-molecular-weight species. After fibril formation in the presence of ZnPc(COONa)₈, immunoprecipitation of Aβ₁₋₄₂ using Aβ-specific antibody followed by near-infrared scanning demonstrated binding of ZnPc(COONa)₈ to Aβ₁₋₄₂. A study using the hydrophobic fluorescent probe 8-anilino-1-naphthalenesulfonic acid showed that ZnPc(COONa)8 decreased the hydrophobicity during Aβ₁₋₄₂ fibril formation. CD spectroscopy showed an increase in the α helix structure and a decrease in the β sheet structure of Aβ₁₋₄₀ in fibril-forming buffer containing ZnPc(COONa)₈. SDS/PAGE and a dot-blot immunoassay showed that ZnPc(COONa)₈ delayed the disappearance of low-molecular-weight species and the appearance of higher-molecular-weight oligomeric species of Aβ₁₋₄₂. A cell viability assay showed that ZnPc(COONa)₈ was not toxic to a neuronal cell line (A1), but instead protected A1 cells against Aβ₁₋₄₂-induced toxicity. Overall, our results indicate that ZnPc(COONa)₈ binds to Aβ and decreases the hydrophobicity, and this change is unfavorable for Aβ oligomerization and fibril formation.

Cysteine inhibits the fibrillisation and cytotoxicity of amyloid-β 40 and 42: implications for the contribution of the thiophilic interaction

Inhibitors of amyloid fibril formation have been at the centre of intense research efforts for the prevention of amyloidosis. Here, we hypothesise that a specific non-covalent interaction, the thiophilic interaction between the side chain of an aromatic residue in a polypeptide and a sulphur atom of the compound, effectively inhibits amyloid fibril formation. Fluorescence spectroscopy and transmission electron microscopy revealed that sulphur compounds, particularly Cys, inhibit the fibrillisation of amyloid-β 1-40 (Aβ40) and 1-42 (Aβ42). Interestingly, aggregates of Aβ40 and Aβ42 induced by Cys were less cytotoxic than those induced by catechin, which is the most typical inhibitor of amyloid fibril formation. Because the essential amino acid, Cys, is an abundant molecule in the blood and cytosol, our data provide a new basis for the prevention of amyloid-related diseases and the elucidation of the mechanism of these diseases.

Understanding amyloid fibril nucleation and aβ oligomer/drug interactions from computer simulations

Evolution has fine-tuned proteins to accomplish a variety of tasks. Yet, with aging, some proteins assemble into harmful amyloid aggregates associated with neurodegenerative diseases, such as Alzheimer's disease (AD), which presents a complex and costly challenge to our society. Thus, far, drug after drug has failed to slow the progression of AD, characterized by the self-assembly of the 39-43 amino acid β-amyloid (Aβ) protein into extracellular senile plaques that form a cross-β structure. While there is experimental evidence that the Aβ small oligomers are the primary toxic species, standard tools of biology have failed to provide structures of these transient, inhomogeneous assemblies. Despite extensive experimental studies, researchers have not successfully characterized the nucleus ensemble, the starting point for rapid fibril formation. Similarly scientists do not have atomic data to show how the compounds that reduce both fibril formation and toxicity in cells bind to Aβ42 oligomers. In this context, computer simulations are important tools for gaining insights into the self-assembly of amyloid peptides and the molecular mechanism of inhibitors. This Account reviews what analytical models and simulations at different time and length scales tell us about the dynamics, kinetics, and thermodynamics of amyloid fibril formation and, notably, the nucleation process. Though coarse-grained and mesoscopic protein models approximate atomistic details by averaging out unimportant degrees of freedom, they provide generic features of amyloid formation and insights into mechanistic details of the self-assembly process. The thermodynamics and kinetics vary from linear peptides adopting straight β-strands in fibrils to longer peptides adopting in parallel U shaped conformations in fibrils. In addition, these properties change with the balance between electrostatic and hydrophobic interactions and the intrinsic disorder of the system. However, simulations suggest that the critical nucleus size might be on the order of 20 chains under physiological conditions. The transition state might be characterized by a simultaneous change from mixed antiparallel/parallel β-strands with random side-chain packing to the final antiparallel or parallel states with the steric zipper packing of the side chains. Second, we review our current computer-based knowledge of the 3D structures of inhibitors with Aβ42 monomer and oligomers, a prerequisite for developing new drugs against AD. Recent extensive all-atom simulations of Aβ42 dimers with known inhibitors such as the green tea compound epigallocatechin-3-gallate and 1,4-naphthoquinon-2-yl-l-tryptophan provide a spectrum of initial Aβ42/inhibitor structures useful for screening and drug design. We conclude by discussing future directions that may offer opportunities to fully understand nucleation and further AD drug development.

Combination of theoretical and experimental approaches for the design and study of fibril-forming peptides

Self-assembling peptides that can form supramolecular structures such as fibrils, ribbons, and nanotubes are of particular interest to modern bionanotechnology and materials science. Their ability to form biocompatible nanostructures under mild conditions through non-covalent interactions offers a big biofabrication advantage. Structural motifs extracted from natural proteins are an important source of inspiration for the rational design of such peptides. Examples include designer self-assembling peptides that correspond to natural coiled-coil motifs, amyloid-forming proteins, and natural fibrous proteins. In this chapter, we focus on the exploitation of structural information from beta-structured natural fibers. We review a case study of short peptides that correspond to sequences from the adenovirus fiber shaft. We describe both theoretical methods for the study of their self-assembly potential and basic experimental protocols for the assessment of fibril-forming assembly.

Computational assembly of polymorphic amyloid fibrils reveals stable aggregates

Amyloid proteins aggregate into polymorphic fibrils that damage tissues of the brain, nerves, and heart. Experimental and computational studies have examined the structural basis and the nucleation of short fibrils, but the ability to predict and precisely quantify the stability of larger aggregates has remained elusive. We established a complete classification of fibril shapes and developed a tool called CreateFibril to build such complex, polymorphic, modular structures automatically. We applied stability landscapes, a technique we developed to reveal reliable fibril structural parameters, to assess fibril stability. CreateFibril constructed HET-s, Aβ, and amylin fibrils up to 17 nm in length, and utilized a novel dipolar solvent model that captured the effect of dipole-dipole interactions between water and very large molecular systems to assess their aqueous stability. Our results validate experimental data for HET-s and Aβ, and suggest novel (to our knowledge) findings for amylin. In particular, we predicted the correct structural parameters (rotation angles, packing distances, hydrogen bond lengths, and helical pitches) for the one and three predominant HET-s protofilaments. We reveal and structurally characterize all known Aβ polymorphic fibrils, including structures recently classified as wrapped fibrils. Finally, we elucidate the predominant amylin fibrils and assert that native amylin is more stable than its amyloid form. CreateFibril and a database of all stable polymorphic fibril models we tested, along with their structural energy landscapes, are available at http://amyloid.cs.mcgill.ca.

Modelling and analysis of early aggregation events of BMHP1-derived self-assembling peptides

Despite the increasing use and development of peptide-based scaffolds in different fields including that of regenerative medicine, the understanding of the factors governing the self-assembly process and the relationship between sequence and properties have not yet been fully understood. BMHP1-derived self-assembling peptides (SAPs) have been developed and characterized showing that biotinylation at the N-terminal cap corresponds to better performing assembly and scaffold biomechanics. In this study, the effects of biotinylation on the self-assembly dynamics of seven BMHP1-derived SAPs have been investigated by molecular dynamics simulations. We confirmed that these SAPs self-assemble into β-structures and that proline acts as a β-breaker of the assembled aggregates. In biotinylated peptides, the formation of ordered β-structured aggregates is triggered by both the establishment of a dense and dynamic H-bonds network and the formation of a 'hydrophobic wall' available to interact with other peptides. Such conditions result from the peculiar chemical composition of the biotinyl-cap, given by the synergic cooperation of the uracil function of the ureido ring with the high hydrophobic portion consisting of the thiophenyl ring and valeryl chain. The inbuilt propensity of biotinylated peptides towards the formation of ordered small aggregates makes them ideal precursors of higher hierarchically organized self-assembled nanostructures as experimentally observed.

Conformational ensemble and polymorphism of the all-atom Alzheimer's Aβ(37-42) amyloid peptide oligomers

Although the Aβ(37-42) peptide has two opposite terminal charges, counterintuitively its current fibril amyloid structure reveals in register parallel β-strands, as formed by the full length Aβ peptide. In this study, we carried out a replica exchange molecular dynamics simulation of 16 all-atom Aβ(37-42) peptides in explicit water starting from randomized and dispersed chains. The extensive conformational sampling (48 replicas, 460 ns/replica) with a total simulation time of 23 μs allows us to obtain a full picture on the equilibrium conformational distribution of oligomers and β-sheet sizes and gain some insights into the oligomerization process at 300 K. At the peptide concentration of 12 mM, self-assembly is described by the condensation-polymerization mechanism with conversion from micelle-like to high β-sheet structures. At equilibrium, the oligomer distribution consists of large aggregates and free monomers, representing 70% and 25% of all species, respectively. Though the formation/dissociation of β-strand is high, the population of 4-5 fully parallel β-strands, consistent with the arrangement in the current fibril, is marginal and that of 4-5 fully antiparallel β-strands, consistent with amyloid polymorphism, is non-negligible. However, the system adopts essentially mixed parallel/antiparallel β-strands. This indicates that a system of 16 Aβ(37-42) chains in explicit solvent still does not form more stable species that will irreversibly grow to a fibril, independently of polymorphism. Our results also suggest that the Aβ(37-42) fibril may display packing polymorphism with antiparallel β-strands, in addition to the experimentally observed in register parallel β-strands.

β-sheet propensity controls the kinetic pathways and morphologies of seeded peptide aggregation

The effect of seeds in templating the morphology of peptide aggregates is examined using molecular dynamics simulations and a coarse-grained peptide representation. Varying the nature of the aggregate seed between β-sheet, amorphous, and β-barrel seeds leads to different aggregation pathways and to morphologically different aggregates. Similar effects are seen by varying the β-sheet propensity of the free peptides. For a fibrillar seed and free peptides of high β-sheet propensity, fibrillar growth occurred by means of direct attachment (without structural rearrangement) of free individual peptides and small ordered oligomers onto the seed. For a fibrillar seed and free peptides of low β-sheet propensity, fibrillar growth occurred through a dock-lock mechanism, in which the free peptides first docked onto the seed, and then locked on, extending and aligning to join the fibril. Amorphous seeds absorbed free peptides into themselves indiscriminately, with any fibrillar rearrangement subsequent to this absorption by means of a condensation-ordering transition. Although the mechanisms observed by varying peptide β-sheet propensity are diverse, the initial pathways can always be broken down into the following steps: (i) the free peptides diffuse in the bulk and attach individually to the seed; (ii) the free peptides diffuse and aggregate among themselves; (iii) the free peptide oligomers collide with the seed; and (iv) the free oligomers merge with the seed and rearrange in a manner dependent on the backbone flexibility of both the free and seed peptides. Our simulations indicate that it is possible to sequester peptides from amorphous aggregates into fibrils, and also that aggregate morphology (and thus cytoxicity) can be controlled by introducing seeds of aggregate-compatible peptides with differing β-sheet propensities into the system.

New method for determining size of critical nucleus of fibril formation of polypeptide chains

A new method for determining the size of critical nucleus of fibril formation of polypeptide chains is proposed. Based on the hypothesis that the fibril grows by addition of a nascent peptide to the preformed template, the nucleus size N(c) is defined as the number of forming template peptides above which the time to add a new monomer becomes independent of the template size. Using lattice models one can show that our method and the standard method which is based on calculation of the free energy, provide the same result for N(c).

Interplay between drying and stability of a TIM barrel protein: a combined simulation-experimental study

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.