Dual effect of crowders on fibrillation kinetics of polypeptide chains revealed by lattice models

Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation

Alzheimer’s disease is understood to be caused by amyloid fibrils and oligomers formed by aggregated amyloid-β (Aβ) peptides. This review article presents molecular dynamics (MD) simulation studies of Aβ peptides and Aβ fragments on their aggregation, aggregation inhibition, amyloid fibril conformations in equilibrium, and disruption of the amyloid fibril by ultrasonic wave and infrared laser irradiation. In the aggregation of Aβ, a β-hairpin structure promotes the formation of intermolecular β-sheet structures. Aβ peptides tend to exist at hydrophilic/hydrophobic interfaces and form more β-hairpin structures than in bulk water. These facts are the reasons why the aggregation is accelerated at the interface. We also explain how polyphenols, which are attracting attention as aggregation inhibitors of Aβ peptides, interact with Aβ. An MD simulation study of the Aβ amyloid fibrils in equilibrium is also presented: the Aβ amyloid fibril has a different structure at one end from that at the other end. The amyloid fibrils can be destroyed by ultrasonic wave and infrared laser irradiation. The molecular mechanisms of these amyloid fibril disruptions are also explained, particularly focusing on the function of water molecules. Finally, we discuss the prospects for developing treatments for Alzheimer’s disease using MD simulations.

Computational Models for the Study of Protein Aggregation

Protein aggregation has been studied by many groups around the world for many years because it can be the cause of a number of neurodegenerative diseases that have no effective treatment. Obtaining the structure of related fibrils and toxic oligomers, as well as describing the pathways and main factors that govern the self-organization process, is of paramount importance, but it is also very difficult. To solve this problem, experimental and computational methods are often combined to get the most out of each method. The effectiveness of the computational approach largely depends on the construction of a reasonable molecular model. Here we discussed different versions of the four most popular all-atom force fields AMBER, CHARMM, GROMOS, and OPLS, which have been developed for folded and intrinsically disordered proteins, or both. Continuous and discrete coarse-grained models, which were mainly used to study the kinetics of aggregation, are also summarized.

Effect of Surface Roughness on Aggregation of Polypeptide Chains: A Monte Carlo Study

The self-assembly of amyloidogenic peptides and proteins into fibrillar structures has been intensively studied for several decades, because it seems to be associated with a number of neurodegenerative diseases, such as Alzheimer’s and Parkinson’s disease. Therefore, understanding the molecular mechanisms of this phenomenon is important for identifying an effective therapy for the corresponding diseases. Protein aggregation in living organisms very often takes place on surfaces like membranes and the impact of a surface on this process depends not only on the surface chemistry but also on its topology. Our goal was to develop a simple lattice model for studying the role of surface roughness in the aggregation kinetics of polypeptide chains and the morphology of aggregates. We showed that, consistent with the experiment, an increase in roughness slows down the fibril formation, and this process becomes inhibited at a very highly level of roughness. We predicted a subtle catalytic effect that a slightly rough surface promotes the self-assembly of polypeptide chains but does not delay it. This effect occurs when the interaction between the surface and polypeptide chains is moderate and can be explained by taking into account the competition between energy and entropy factors.

Heat-induced degradation of fibrils: Exponential vs logistic kinetics

The degradation of fibrils under the influence of thermal fluctuations was studied experimentally by various groups around the world. In the first set of experiments, it was shown that the decay of fibril content, which can be measured by the ThT fluorescence assay, obeys a bi-exponential function. In the second series of experiments, it was demonstrated that when the monomers separated from the aggregate are not recyclable, the time dependence of the number of monomers belonging to the dominant cluster is described by a single-exponential function if the fraction of bound chains becomes less than a certain threshold. Note that the time dependence of the fraction of bound chains can be measured by tryptophan fluorescence. To understand these interesting experimental results, we developed a phenomenological theory and performed molecular simulation. According to our theory and simulations using the lattice and all-atom models, the time dependence of bound chains is described by a logistic function, which slowly decreases at short time scales but becomes a single exponential function at large time scales. The results, obtained by using lattice and all-atom simulations, ascertained that the time dependence of the fibril content can be described by a bi-exponential function that decays faster than the logistic function on short time scales. We have uncovered the molecular mechanism for the distinction between the logistic and bi-exponential behavior. Since the dissociation of the chain from the fibrils requires the breaking of a greater number of inter-chain contacts as compared to the breaking of the beta sheet structure, the decrease in the number of connected chains is slower than the fibril content. Therefore, the time dependence of the aggregate size is logistic, while the two-exponential behavior is preserved for the content of fibrils. Our results are in agreement with the results obtained in both sets of experiments.

Interplay between Amyloid Fibrillation Delay and Degradation by Magnetic Zinc-Doped Ferrite Nanoparticles

Amyloidosis, the aggregation of naturally soluble proteins into fibrils, is the main pathological hallmark of central nervous system (CNS) disorders, and new therapeutic approaches can be introduced through nanotechnology. Herein, magnetic nanoparticles (MNPs) are proposed to combat amyloidosis and act as CNS theranostic (therapy and diagnosis) candidates through magnetomechanical forces that can be induced under a low-frequency magnetic field. In that vein, a modified one-step microwave-assisted polyol process has been employed to synthesize hybrid organic/inorganic zinc ferrite (Zn_xFe_3-xO₄) MNPs with different levels of zinc doping (0.30 < x < 0.6) derived from the utilized polyol. The lowest doped (x = 0.30) MNPs exhibited high magnetization (127 emu/g), high T₂ imaging ability (r₂ = 432 mM^-1 s^-1), and relatively small hydrodynamic size (180 nm), decisive characteristics to further evaluate their CNS theranostic potential. Their effect on the fibrillation/degradation was monitored in two model proteins, insulin and albumin, in the presence/absence of variant external magnetic fields (static, rotating, or alternating) via Thioflavin T (ThT) fluorescence assay and optical fluorescence microscopy. The MNPs were injected either in oligomer solution where significant fibrillation delay was observed, boosted by zinc ionic leaching of MNPs, or in already formed amyloid plaques where up to 86% amyloid degradation was recorded in the presence of magnetic fields, unveiling magnetomechanical antifibrillation properties. The alternating magnetic field (4 Hz) allows the bouncing of the MNPs into the amyloid net driven by the magnetic forces, and thus is featured as the preferred "dancing mode", which strengthens the degrading efficacy of MNPs.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ₄₀ and Aβ₄₂ and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

Physical mechanism for biopolymers to aggregate and maintain in non-equilibrium states

Many human or animal diseases are related to aggregation of proteins. A viable biological organism should maintain in non-equilibrium states. How protein aggregate and why biological organisms can maintain in non-equilibrium states are not well understood. As a first step to understand such complex systems problems, we consider simple model systems containing polymer chains and solvent particles. The strength of the spring to connect two neighboring monomers in a polymer chain is controlled by a parameter s with s → ∞ for rigid-bond. The strengths of bending and torsion angle dependent interactions are controlled by a parameter s A with s A  → -∞ corresponding to no bending and torsion angle dependent interactions. We find that for very small s A , polymer chains tend to aggregate spontaneously and the trend is independent of the strength of spring. For strong springs, the speed distribution of monomers in the parallel (along the direction of the spring to connect two neighboring monomers) and perpendicular directions have different effective temperatures and such systems are in non-equilibrium states.

Finite-size corrections and scaling for the dimer model on the checkerboard lattice

Lattice models are useful for understanding behaviors of interacting complex many-body systems. The lattice dimer model has been proposed to study the adsorption of diatomic molecules on a substrate. Here we analyze the partition function of the dimer model on a 2M×2N checkerboard lattice wrapped on a torus and derive the exact asymptotic expansion of the logarithm of the partition function. We find that the internal energy at the critical point is equal to zero. We also derive the exact finite-size corrections for the free energy, the internal energy, and the specific heat. Using the exact partition function and finite-size corrections for the dimer model on a finite checkerboard lattice, we obtain finite-size scaling functions for the free energy, the internal energy, and the specific heat of the dimer model. We investigate the properties of the specific heat near the critical point and find that the specific-heat pseudocritical point coincides with the critical point of the thermodynamic limit, which means that the specific-heat shift exponent λ is equal to ∞. We have also considered the limit N→∞ for which we obtain the expansion of the free energy for the dimer model on the infinitely long cylinder. From a finite-size analysis we have found that two conformal field theories with the central charges c=1 for the height function description and c=-2 for the construction using a mapping of spanning trees can be used to describe the dimer model on the checkerboard lattice.

Effects of hydrophobic macromolecular crowders on amyloid β (16-22) aggregation

In Alzheimer's disease (AD), the amyloid β (Aβ) peptide aggregates in the brain to form progressively larger oligomers, fibrils, and plaques. The aggregation process is strongly influenced by the presence of other macromolecular species, called crowders, that can exert forces on the proteins. One very common attribute of macromolecular crowders is their hydrophobicity. We examined the effect of hydrophobic crowders on protein aggregation by using discontinuous molecular dynamics (DMD) simulations in combination with an intermediate resolution protein model, PRIME20. The systems considered contained 48 Aβ (16-22) peptides and crowders with diameters of 5 Å, 20 Å, and 40 Å, represented by hard spheres or spheres with square-well/square-shoulder interactions, at a crowder volume fraction of ϕ = 0.10. Results show that low levels of crowder hydrophobicity are capable of increasing the fibrillation lag time and high levels of crowder hydrophobicity can fully prevent the formation of fibrils. The types of structures that remain during the final stages of the simulations are summarized in a global phase diagram that shows fibril, disordered oligomer, or β-sheet phases in the space spanned by crowder size and crowder hydrophobicity. In particular, at high levels of hydrophobicity, simulations with 5 Å crowders result in only disordered oligomers and simulations with 40 Å crowders result in only β-sheets. The presence of hydrophobic crowders reduces the antiparallel β-sheet content of fibrils, whereas hard sphere crowders increase it. Finally, strong hydrophobic crowders alter the secondary structure of the Aβ (16-22) monomers, bending them into a shape that is incapable of forming ordered β-sheets or fibrils. These results qualitatively agree with previous theoretical and experimental work.

Preformed template fluctuations promote fibril formation: insights from lattice and all-atom models

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite the fact that the fibril formation process is very slow and thus poses a significant challenge for theoretical and experimental studies, a number of alternative pictures of molecular mechanisms of amyloid fibril formation have been recently proposed. What seems to be common for the majority of the proposed models is that fibril elongation involves the formation of pre-nucleus seeds prior to the creation of a critical nucleus. Once the size of the pre-nucleus seed reaches the critical nucleus size, its thermal fluctuations are expected to be small and the resulting nucleus provides a template for sequential (one-by-one) accommodation of added monomers. The effect of template fluctuations on fibril formation rates has not been explored either experimentally or theoretically so far. In this paper, we make the first attempt at solving this problem by two sets of simulations. To mimic small template fluctuations, in one set, monomers of the preformed template are kept fixed, while in the other set they are allowed to fluctuate. The kinetics of addition of a new peptide onto the template is explored using all-atom simulations with explicit water and the GROMOS96 43a1 force field and simple lattice models. Our result demonstrates that preformed template fluctuations can modulate protein aggregation rates and pathways. The association of a nascent monomer with the template obeys the kinetics partitioning mechanism where the intermediate state occurs in a fraction of routes to the protofibril. It was shown that template immobility greatly increases the time of incorporating a new peptide into the preformed template compared to the fluctuating template case. This observation has also been confirmed by simulation using lattice models and may be invoked to understand the role of template fluctuations in slowing down fibril elongation in vivo.

Nanopore Diameters Tune Strain in Extruded Fibronectin Fibers

Fibronectin is present in the extracellular matrix and can be assembled into nanofibers in vivo by undergoing conformational changes. Here, we present a novel approach to prepare fibronectin nanofibers under physiological conditions using an extrusion approach through nanoporous aluminum oxide membranes. This one-step process can prepare nanofiber bundles up to a millimeter in length and with uniform fiber diameters in the nanometer range. Most importantly, by using different pore diameters and protein concentrations in the extrusion process, we could induce varying lasting structural changes in the fibers, which were monitored by Förster resonance energy transfer and should impose different physiological functions.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.

Effects of macromolecular crowding on amyloid beta (16-22) aggregation using coarse-grained simulations

To examine the effect of crowding on protein aggregation, discontinuous molecular dynamics (DMD) simulations combined with an intermediate resolution protein model, PRIME20, were applied to a peptide/crowder system. The systems contained 192 Aβ(16–22) peptides and crowders of diameters 5, 20, and 40 Å, represented here by simple hard spheres, at crowder volume fractions of 0.00, 0.10, and 0.20. Results show that both crowder volume fraction and crowder diameter have a large impact on fibril and oligomer formation. The addition of crowders to a system of peptides increases the rate of oligomer formation, shifting from a slow ordered formation of oligomers in the absence of crowders, similar to nucleated polymerization, to a fast collapse of peptides and subsequent rearrangement characteristic of nucleated conformational conversion with a high maximum in the number of peptides in oligomers as the total crowder surface area increases. The rate of conversion from oligomers to fibrils also increases with increasing total crowder surface area, giving rise to an increased rate of fibril growth. In all cases, larger volume fractions and smaller crowders provide the greatest aggregation enhancement effects. We also show that the size of the crowders influences the formation of specific oligomer sizes. In our simulations, the 40 Å crowders enhance the number of dimers relative to the numbers of trimers, hexamers, pentamers, and hexamers, while the 5 Å crowders enhance the number of hexamers relative to the numbers of dimers, trimers, tetramers, and pentamers. These results are in qualitative agreement with previous experimental and theoretical work.

Transformation between α-helix and β-sheet structures of one and two polyglutamine peptides in explicit water molecules by replica-exchange molecular dynamics simulations

Aggregation of polyglutamine peptides with β-sheet structures is related to some important neurodegenerative diseases such as Huntington's disease. However, it is not clear how polyglutamine peptides form the β-sheets and aggregate. To understand this problem, we performed all-atom replica-exchange molecular dynamics simulations of one and two polyglutamine peptides with 10 glutamine residues in explicit water molecules. Our results show that two polyglutamine peptides mainly formed helix or coil structures when they are separated, as in the system with one-polyglutamine peptide. As the interpeptide distance decreases, the intrapeptide β-sheet structure sometimes appear as an intermediate state, and finally the interpeptide β-sheets are formed. We also find that the polyglutamine dimer tends to form the antiparallel β-sheet conformations rather than the parallel β-sheet, which is consistent with previous experiments and a coarse-grained molecular dynamics simulation.

Transition-metal-substituted polyoxometalate derivatives as functional anti-amyloid agents for Alzheimer's disease

Inhibitions of amyloid β (Aβ) aggregation and Aβ-haem peroxidase-like activity have received much attention because these two symptoms can be the primary targets of therapeutic strategies for Alzheimer's disease (AD). Recently, our group found that polyoxometalate (POM) with a Wells-Dawson structure can efficiently inhibit Aβ aggregation. However, the interaction between POMs and Aβ is robust, but still needs to improve Aβ binding affinity. More importantly, it is unclear whether POMs can cross the blood-brain barrier and decrease Aβ-haem peroxidase-like activity. Here we show that our designed series of transition metal-functionalized POM derivatives with a defined histidine-chelated binding site have much better Aβ inhibition and peroxidase-like activity inhibition effects than the parent POM. More intriguingly, we show that these compounds can cross the blood-brain barrier and are metabolized after 48 h. Our work provides insights into the design, synthesis and screening of inorganic metal compounds as multifunctional therapeutic agents against AD.