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

In Silico View of Crowding: Biomolecular Processes to Nanomaterial Design

It is widely accepted that deciphering biomolecular structure and function requires going beyond the single-molecule or single-complex paradigm. The densely packed macromolecules, cosolutes, and metabolites in the living cell impose crowding effects on the biomolecular structure and dynamics that need to be accounted for. Molecular simulations have proven to be a powerful tool to advance the current molecular-level understanding of such a highly concentrated, complex milieu. This Mini-Review focuses on summarizing the understanding achieved so far for the effects of crowding on biomolecular processes using computational methods, along with highlighting a new direction in employing crowding as a tool for tunable nanomaterial design. The two schools of thought that form the pillars of the current understanding of crowding effects are discussed. The investigation of crowded solutions using physics-based models that encompass different time and length scales to mimic the intracellular environment are described. The limitations and challenges faced by the current models and simulation methods are addressed, highlighting the gaps to be filled for better agreement with experiments. Crowding can also act as an effective tool to modulate the structure-property-function relationships of nanomaterials, leading to the development of novel functional materials. A few recent studies, mostly experimental, have been summarized in this direction. The Mini-Review concludes with an outlook for future developments in this field in order to enable accurate mimicking of the intracellular environment using simulations and to bridge the gap between biological processes and nanomaterial design.

Ethylene glycol energetically disfavours oligomerization of pseudoisocyanine dyestuffs at crowded concentrations

The intriguing role of the intracellular crowded environment in regulating protein aggregation remains elusive. The convolution of several factors such as the protein sequence-dependence, crowder's shape and size and diverse intermolecular interactions makes it complex to identify systematic trends. One of the ways to simplify the problem is to study a synthetic model for self-assembling proteins. In this study, we examine the aggregation behaviour of the cationic pseudoisocyanine chloride (PIC) dyestuff which is known to self-assemble and form fibril-like J-aggregates in aqueous solutions, similar to those formed by amyloid-forming proteins. Prior experimental studies have shown that polyethylene glycol impedes and Ficoll-400 promotes the self-assembly of PIC dyes. To achieve molecular insights, we examine the effect of crowding by ethylene glycol on the solvation thermodynamics of oligomerization of dyes into H-type and J-type oligomers using extensive molecular dynamics simulations. The binding free energy calculations show that the formation of J-oligomers is more favourable than that of H-oligomers in water. The stability of H- and J- tetramers and pentamers decreases in crowded solutions. The formation of oligomers is supported by the favourable change in dye-solvent interaction energy in both pure water and aqueous ethylene glycol solution although it is opposed by the reduced dye-solvent entropy. Ethylene glycol, as a molecular crowder, disfavours the H- as well as J-oligomerization via preferential binding to the dye oligomers. An unfavourable change in dye-crowder and dye-dye interaction energy on dye association makes the H-oligomer formation less favourable in crowded solution than in pure water solution. In the case of J-oligomers, however, the unfavourable change in dye-crowder interaction energy primarily contributes to making total dye-solvent energy unfavourable. The results are supported by isothermal titration calorimetry measurements where the binding of ethylene glycol to PIC molecules is found to be endothermic. The results provide an emerging view that a crowded environment can disfavour self-assembly of PIC dyes by interactions with the oligomeric states. The findings have implications in understanding the role of a crowded environment in shaping the free energy landscapes of proteins.

The conformational phase diagram of neutral polymers in the presence of attractive crowders

Extensive coarse-grained molecular dynamics simulations are performed to investigate the conformational phase diagram of a neutral polymer in the presence of attractive crowders. We show that, for low crowder densities, the polymer predominantly shows three phases as a function of both intra-polymer and polymer-crowder interactions: (1) weak intra-polymer and weak polymer-crowder attractive interactions induce extended or coil polymer conformations (phase E), (2) strong intra-polymer and relatively weak polymer-crowder attractive interactions induce collapsed or globular conformations (phase CI), and (3) strong polymer-crowder attractive interactions, regardless of intra-polymer interactions, induce a second collapsed or globular conformation that encloses bridging crowders (phase CB). The detailed phase diagram is obtained by determining the phase boundaries delineating the different phases based on an analysis of the radius of gyration as well as bridging crowders. The dependence of the phase diagram on strength of crowder-crowder attractive interactions and crowder density is clarified. We also show that when the crowder density is increased, a third collapsed phase of the polymer emerges for weak intra-polymer attractive interactions. This crowder density-induced compaction is shown to be enhanced by stronger crowder-crowder attraction and is different from the depletion-induced collapse mechanism, which is primarily driven by repulsive interactions. We also provide a unified explanation of the observed re-entrant swollen/extended conformations of the earlier simulations of weak and strongly self-interacting polymers in terms of crowder-crowder attractive interactions.

Influencing factors and characterization methods of nanoparticles regulating amyloid aggregation

Human disorders associated with amyloid aggregation, such as Alzheimer's disease and Parkinson's disease, afflict the lives of millions worldwide. When peptides and proteins in the body are converted to amyloids, which have a tendency to aggregate, the toxic oligomers produced during the aggregation process can trigger a range of diseases. Nanoparticles (NPs) have been found to possess surface effects that can modulate the amyloid aggregation process and they have potential application value in the treatment of diseases related to amyloid aggregation and fibrillary tangles. In this review, we discuss recent progress relating to studies of nanoparticles that regulate amyloid aggregation. The review focuses on the factors influencing this regulation, which are important as guidelines for the future design of NPs for the treatment of amyloid aggregation. We describe the characterization methods that have been utilized so far in such studies. This review provides research information and characterization methods for the rational design of NPs, which should result in therapeutic strategies for amyloid diseases.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Chaperone potential of erythroid spectrin: Effects of hemoglobin interaction, macromolecular crowders, phosphorylation and glycation

Spectrin, the major protein component of the erythrocyte membrane skeleton has chaperone like activity and is known to bind membrane phospholipids and hemoglobin. We have probed the chaperone activity of spectrin in presence of hemoglobin and phospholipid SUVs of different compositions to elucidate the effect of phospholipid/hemoglobin binding on chaperone function. It is seen that spectrin displays a preference for hemoglobin over other substrates leading to a decrease in chaperone activity in presence of hemoglobin. A competition is seen to exist between phospholipid binding and chaperone function of spectrin, in a dose dependent manner with the greatest extent of decrease being seen in case of phospholipid vesicles containing aminophospholipids e.g. PS and PE which may have implications in diseases like hereditary spherocytosis where mutation in spectrin is implicated in its detachment from cell membrane. To gain a clearer understanding of the chaperone like activity of spectrin under in-vivo like conditions we have investigated the effect of macromolecular crowders as well as phosphorylation and glycation states on chaperone activity. It is seen that the presence of non-specific, protein and non-protein macromolecular crowders do not appreciably affect chaperone function. Phosphorylation also does not affect the chaperone function unlike glycation which progressively diminishes chaperone activity. We propose a model where chaperone clients adsorb onto spectrin's surface and processes that bind to and occlude these surfaces decrease chaperone activity.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

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.

Structural Analysis of a Trimer of β2-Microgloblin Fragment by Molecular Dynamics Simulations

A peptide β₂-m_21-31, which is a fragment from residue 21 to residue 31 of β₂-microgloblin, is experimentally known to self-assemble and form amyloid fibrils. In order to understand the mechanism of amyloid fibril formations, we applied the replica-exchange molecular dynamics method to the system consisting of three fragments of β₂-m_21-31. From the analyses on the temperature dependence, we found that there is a clear phase transition temperature in which the peptides aggregate with each other. Moreover, we found by the free energy analyses that there are two major stable states: One of them is like amyloid fibrils and the other is amorphous aggregates.

Seeding and cross-seeding fibrillation of N-terminal prion protein peptides PrP(120-144)

Prion diseases are infectious neurodegenerative diseases that are capable of cross-species transmission, thus arousing public health concerns. Seed-templating propagation of prion protein is believed to underlie prion cross-species transmission pathology. Understanding the molecular fundamentals of prion propagation is key to unravelling the pathology of prion diseases. In this study, we use coarse-grained molecular dynamics to investigate the seeding and cross-seeding aggregation of three prion protein fragments PrP(120-144) originating from human (Hu), bank vole (BV), and Syrian hamster (SHa). We find that the seed accelerates the aggregation of the monomer peptides by eliminating the lag phase. The monomer aggregation kinetics are mainly determined by the structure of the seed. The stronger the hydrophobic residues on the seed associate with each other, the higher the probability that the seed recruits monomer peptides to its surface/interface. For cross-seeding aggregation, we show that Hu has a strong tendency to adopt the conformation of the BV seed and vice versa; the Hu and BV monomers have a weak tendency to adopt the conformation of the SHa seed. These two findings are consistent with Apostol et al.'s experimental findings on PrP(138-143) and partially consistent with Jones et al.'s finding on PrP(23-144). We also identify several conformational mismatches when SHa cross-seeds BV and Hu peptides, indicating the existence of a cross-seeding barrier between SHa and the other two sequences. This study sheds light on the molecular mechanism of seed-templating aggregation of prion protein fragments underlying the sequence-dependent transmission barrier in prion diseases.

Assemblies of amyloid-β30-36 hexamer and its G33V/L34T mutants by replica-exchange molecular dynamics simulation

The aggregation of amyloid-β peptides is associated with the pathogenesis of Alzheimer’s disease, in which the 30–36 fragments play an important part as a fiber-forming hydrophobic region. The fibrillar structure of Aβ_30–36 has been detected by means of X-ray diffraction, but its oligomeric structural determination, biophysical characterization, and pathological mechanism remain elusive. In this study, we have investigated the structures of Aβ_30–36 hexamer as well as its G33V and L34T mutants in explicit water environment using replica-exchange molecular dynamics (REMD) simulations. Our results show that the wild-type (WT) Aβ_30–36 hexamer has a preference to form β-barrel and bilayer β-sheet conformations, while the G33V or L34T mutation disrupts the β-barrel structures: the G33V mutant is homogenized to adopt β-sheet-rich bilayers, and the structures of L34T mutant on the contrary get more diverse. The hydrophobic interaction plays a critical role in the formation and stability of oligomeric assemblies among all the three systems. In addition, the substitution of G33 by V reduces the β-sheet content in the most populated conformations of Aβ_30–36 oligomers through a steric effect. The L34T mutation disturbs the interpeptide hydrogen bonding network, and results in the increased coil content and morphological diversity. Our REMD runs provide structural details of WT and G33V/L34T mutant Aβ_30–36 oligomers, and molecular insight into the aggregation mechanism, which will be helpful for designing novel inhibitors or amyloid-based materials.

N-terminal Prion Protein Peptides (PrP(120-144)) Form Parallel In-register β-Sheets via Multiple Nucleation-dependent Pathways

The prion diseases are a family of fatal neurodegenerative diseases associated with the misfolding and accumulation of normal prion protein (PrP^C) into its pathogenic scrapie form (PrP^Sc). Understanding the fundamentals of prion protein aggregation and the molecular architecture of PrP^Sc is key to unraveling the pathology of prion diseases. Our work investigates the early-stage aggregation of three prion protein peptides, corresponding to residues 120-144 of human (Hu), bank vole (BV), and Syrian hamster (SHa) prion protein, from disordered monomers to β-sheet-rich fibrillar structures. Using 12 μs discontinuous molecular dynamics simulations combined with the PRIME20 force field, we find that the Hu-, BV-, and SHaPrP(120-144) aggregate via multiple nucleation-dependent pathways to form U-shaped, S-shaped, and Ω-shaped protofilaments. The S-shaped HuPrP(120-144) protofilament is similar to the amyloid core structure of HuPrP(112-141) predicted by Zweckstetter. HuPrP(120-144) has a shorter aggregation lag phase than BVPrP(120-144) followed by SHaPrP(120-144), consistent with experimental findings. Two amino acid substitutions I138M and I139M retard the formation of parallel in-register β-sheet dimers during the nucleation stage by increasing side chain-side chain association and reducing side chain interaction specificity. On average, HuPrP(120-144) aggregates contain more parallel β-sheet content than those formed by BV- and SHaPrP(120-144). Deletion of the C-terminal residues 138-144 prevents formation of fibrillar structures in agreement with the experiment. This work sheds light on the amyloid core structures underlying prion strains and how I138M, I139M, and S143N affect prion protein aggregation kinetics.

Peptide folding in the presence of interacting protein crowders

Using Monte Carlo methods, we explore and compare the effects of two protein crowders, BPTI and GB1, on the folding thermodynamics of two peptides, the compact helical trp-cage and the β-hairpin-forming GB1m3. The thermally highly stable crowder proteins are modeled using a fixed backbone and rotatable side-chains, whereas the peptides are free to fold and unfold. In the simulations, the crowder proteins tend to distort the trp-cage fold, while having a stabilizing effect on GB1m3. The extent of the effects on a given peptide depends on the crowder type. Due to a sticky patch on its surface, BPTI causes larger changes than GB1 in the melting properties of the peptides. The observed effects on the peptides stem largely from attractive and specific interactions with the crowder surfaces, and differ from those seen in reference simulations with purely steric crowder particles.

Graphene oxide inhibits hIAPP amyloid fibrillation and toxicity in insulin-producing NIT-1 cells

Human islet amyloid polypeptide (hIAPP or amylin) aggregation is directly associated with pancreatic β-cell death and subsequent insulin deficiency in type 2 diabetes (T2D). Since no cure is currently available for T2D, it is of great benefit to devise new anti-aggregation molecules, which protect β-cells against hIAPP aggregation-induced toxicity. Engineered nanoparticles have been recently exploited as anti-aggregation nanomedicines. In this work, we studied graphene oxide (GO) nanosheets for their potential for hIAPP aggregation inhibition by combining computational modeling, biophysical characterization and cell toxicity measurements. Using discrete molecular dynamics (DMD) simulations and in vitro studies, we showed that GO exhibited an inhibitory effect on hIAPP aggregation. DMD simulations indicated that the strong binding of hIAPP to GO nanosheets was driven by hydrogen bonding and aromatic stacking and that the strong peptide-GO binding efficiently inhibited hIAPP self-association and aggregation on the nanosheet surface. Secondary structural changes of hIAPP upon GO binding derived from DMD simulations were consistent with circular dichroism (CD) spectroscopy measurements. Transmission electron microscopy (TEM) images confirmed the reduction of hIAPP aggregation in the presence of GO. Furthermore, we carried out a cell toxicity assay and found that these nanosheets protected insulin-secreting NIT-1 pancreatic β-cells against hIAPP-induced toxicity. Our multidisciplinary study suggests that GO nanosheets have the potential to be utilized as an anti-aggregation nanomedicine itself in addition to a biosensor or delivery vehicle for the mitigation of T2D progression.