Emergence of Alternative Structures in Amyloid Beta 1-42 Monomeric Landscape by N-terminal Hexapeptide Amyloid Inhibitors

Amyloidogenesis: What Do We Know So Far?

The study of protein aggregation, and amyloidosis in particular, has gained considerable interest in recent times. Several neurodegenerative diseases, such as Alzheimer's (AD) and Parkinson's (PD) show a characteristic buildup of proteinaceous aggregates in several organs, especially the brain. Despite the enormous upsurge in research articles in this arena, it would not be incorrect to say that we still lack a crystal-clear idea surrounding these notorious aggregates. In this review, we attempt to present a holistic picture on protein aggregation and amyloids in particular. Using a chronological order of discoveries, we present the case of amyloids right from the onset of their discovery, various biophysical techniques, including analysis of the structure, the mechanisms and kinetics of the formation of amyloids. We have discussed important questions on whether aggregation and amyloidosis are restricted to a subset of specific proteins or more broadly influenced by the biophysiochemical and cellular environment. The therapeutic strategies and the significant failure rate of drugs in clinical trials pertaining to these neurodegenerative diseases have been also discussed at length. At a time when the COVID-19 pandemic has hit the globe hard, the review also discusses the plausibility of the far-reaching consequences posed by the virus, such as triggering early onset of amyloidosis. Finally, the application(s) of amyloids as useful biomaterials has also been discussed briefly in this review.

Study on the Kinetics and Mechanism of Ferrocene-Tripeptide Inhibiting Insulin Aggregation

Peptides are gaining popularity as neurodegenerative disease-targeted drugs due to their medicinal value and the simplicity in the biomedicine and pharmaceutical industry field. In this study, based on previously studied...

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.

Sequence-dependent aggregation-prone conformations of islet amyloid polypeptide

Amyloid proteins, which aggregate to form highly ordered structures, play a crucial role in various disease pathologies. Despite many previous studies on amyloid fibrils, which are an end product of protein aggregation, the structural characteristics of amyloid proteins in the early stage of aggregation and their related aggregation mechanism still remain elusive. The role of the amino acid sequence in the aggregation-prone structures of amyloid proteins at such a stage is not understood. Here, we have studied the sequence-dependent structural characteristics of islet amyloid polypeptide based on atomistic simulations and spectroscopic experiments. We show that the amino acid sequence determines non-bonded interactions that play a leading role in the formation of aggregation-prone conformations. Specifically, a single point mutation critically changes the population of aggregation-prone conformations, resulting in a change of the aggregation mechanism. Our simulation results were supported by experimental results suggesting that mutation affects the kinetics of aggregation and the structural characteristics of amyloid aggregates. Our study provides an insight into the role of sequence-dependent aggregation-prone conformations in the underlying mechanisms of amyloid aggregation.

Molecular insight into the early stage of amyloid-β(1-42) Homodimers aggregation influenced by histidine tautomerism

Aggregated amyloid β-peptide (Aβ) in small oligomeric forms inside the brain causes synaptic function disruption and the development of Alzheimer's disease (AD). Histidine is an important amino acid that may lead to structural changes. Aβ42 monomer chain includes 3 histidine residues that considering two ε and δ tautomers 8 isomers, including (εεε) and (εδδ) could be formed. Molecular dynamics simulation on homodimerization of (εεε) (the most common type of tautomers) and (εδδ) tautomers with different initial configurations using monomer chains from our previous work were performed to uncover the tautomeric behavior of histidine on Aβ42 aggregation in a physiological pH which is still largely unknown and impossible to observe experimentally. We found a higher propensity of forming β-sheet in (εδδ) homodimers and specifically in a greater amount from Aβ42 than from Aβ40. A smaller amount of β-sheet formation was observed for (εεε) homodimers compared with (εδδ). Additionally, interactions in (εδδ) homodimers may indicate the importance of the hydrophobic core and C-/N-terminals during oligomerization. Our findings indicate the important role of the tautomeric effect of histidine and (εδδ) homodimers at the early stage of Aβ aggregation.

Structure-based virtual screening for novel potential selective inhibitors of class IIa histone deacetylases for cancer treatment

The fundamental cause of human cancer is strongly influenced by down- or up-regulations of epigenetic factors. Upregulated histone deacetylases (HDAC) have been shown to be effectively neutralized by the action of HDACs inhibitors (HDACi). However, cytotoxicity has been reported in normal cells because of non-specificity of several available HDACis that are in clinical use or at different phases of clinical trials. Because of the high amino acid sequence and structural similarity among HDAC enzymes, it is believed to be a challenging task to obtain isoform-selectivity. The essential aim of the present research work was to identify isoform-selective inhibitors against class IIa HDACs via structure-based drug design. Based on the highest binding affinity and isoform-selectivity, the top-ranked inhibitors were in silico tested for their absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties, which were classified as drug-like compounds. Later, molecular dynamics simulation (MD) was carried out for all compound-protein complexes to evaluate the structural stability and the biding mode of the inhibitors, which showed high stability throughout the 100 ns simulation. Free binding energy predictions by MM-PBSA method showed the high binding affinity of the identified compounds toward their respective targets. Hence, these inhibitors could be used as drug candidates or as lead compounds for more in silico or in vitro optimization to design safe isoform-selective HDACs inhibitors.

Effects of Ions and Small Compounds on the Structure of Aβ42 Monomers

Aggregation of amyloid-β (Aβ) proteins in the brain is a hallmark of Alzheimer's disease. This phenomenon can be promoted or inhibited by adding small molecules to the solution where Aβ is embedded. These molecules affect the ensemble of conformations sampled by Aβ monomers even before aggregation starts. Here, we perform extensive all-atom replica exchange molecular dynamics (REMD) simulations to provide a comparative study of the ensemble of conformations sampled by Aβ₄₂ monomers in solutions that promote (i.e., aqueous solution containing NaCl) and inhibit (i.e., aqueous solutions containing scyllo-inositol or 4-aminophenol) aggregation. Simulations performed in pure water are used as our reference. We find that secondary-structure content is only affected in an antagonistic manner by promoters and inhibitors at the C-terminus and the central hydrophilic core. Moreover, the end of the C-terminus binds more favorably to the central hydrophobic core region of Aβ₄₂ in NaCl adopting a type of strand-loop-strand structure that is disfavored by inhibitors. Nonpolar residues that form the dry core of larger aggregates of Aβ₄₂ (e.g., PDB ID 2BEG) are found at close proximity in these strand-loop-strand structures, suggesting that their formation could play an important role in initiating nucleation. In the presence of inhibitors, the C-terminus binds the central hydrophilic core with a higher probability than in our reference simulation. This sensitivity of the C-terminus, which is affected in an antagonistic manner by inhibitors and promoters, provides evidence for its critical role in accounting for aggregation.

Neurodegenerative Proteinopathies in the Proteoform Spectrum-Tools and Challenges

Proteinopathy refers to a group of disorders defined by depositions of amyloids within living tissue. Neurodegenerative proteinopathies, including Alzheimer's disease, Parkinson's disease, Creutzfeldt-Jakob disease, and others, constitute a large fraction of these disorders. Amyloids are highly insoluble, ordered, stable, beta-sheet rich proteins. The emerging theory about the pathophysiology of neurodegenerative proteinopathies suggests that the primary amyloid-forming proteins, also known as the prion-like proteins, may exist as multiple proteoforms that contribute differentially towards the disease prognosis. It is therefore necessary to resolve these disorders on the level of proteoforms rather than the proteome. The transient and hydrophobic nature of amyloid-forming proteins and the minor post-translational alterations that lead to the formation of proteoforms require the use of highly sensitive and specialized techniques. Several conventional techniques, like gel electrophoresis and conventional mass spectrometry, have been modified to accommodate the proteoform theory and prion-like proteins. Several new ones, like imaging mass spectrometry, have also emerged. This review aims to discuss the proteoform theory of neurodegenerative disorders along with the utility of these proteomic techniques for the study of highly insoluble proteins and their associated proteoforms.

Time-Resolved Observation of Evolution of Amyloid-β Oligomer with Temporary Salt Crystals

The aggregation behavior of amyloid-β (Aβ) peptides remains unclarified despite the fact that it is closely related to the pathogenic mechanism of Alzheimer's disease. Aβ peptides form diverse oligomers with various diameters before nucleation, making clarification of the mechanism involved a complex problem with conventional macroscopic analysis methods. Time-resolved single-molecule level analysis in bulk solution is thus required to fully understand their early stage aggregation behavior. Here, we perform time-resolved observation of the aggregation dynamics of Aβ oligomers in bulk solution using liquid-state transmission electron microscopy. Our observations reveal previously unknown behaviors. The most important discovery is that a salt crystal can precipitate even with a concentration much lower than its solubility, and it then dissolves in a short time, during which the aggregation reaction of Aβ peptides is significantly accelerated. These findings will provide new insights in the evolution of the Aβ oligomer.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

Exploring the Multi-Target Performance of Mitochondriotropic Antioxidants against the Pivotal Alzheimer's Disease Pathophysiological Hallmarks

Alzheimer disease (AD) is the most common neurodegenerative disease featuring progressive and degenerative neurological impairments resulting in memory loss and cognitive decline. The specific mechanisms underlying AD are still poorly understood, but it is suggested that a deficiency in the brain neurotransmitter acetylcholine, the deposition of insoluble aggregates of fibrillar β-amyloid 1–42 (Aβ₄₂), and iron and glutamate accumulation play an important role in the disease progress. Despite the existence of approved cholinergic drugs, none of them demonstrated effectiveness in modifying disease progression. Accordingly, the development of new chemical entities acting on more than one target is attracting progressively more attention as they can tackle intricate network targets and modulate their effects. Within this endeavor, a series of mitochondriotropic antioxidants inspired on hydroxycinnamic (HCA’s) scaffold were synthesized, screened toward cholinesterases and evaluated as neuroprotectors in a differentiated human SH-SY5Y cell line. From the series, compounds 7 and 11 with a 10-carbon chain can be viewed as multi-target leads for the treatment of AD, as they act as dual and bifunctional cholinesterase inhibitors and prevent the neuronal damage caused by diverse aggressors related to protein misfolding and aggregation, iron accumulation and excitotoxicity.

Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson–Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ_RMSE = 1.25 kcal mol⁻¹. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses.

Modulation of Amyloid-β42 Conformation by Small Molecules Through Nonspecific Binding

Aggregation of amyloid-β (Aβ) peptides is a crucial step in the progression of Alzheimer's disease (AD). Identifying aggregation inhibitors against AD has been a great challenge. We report an atomistic simulation study of the inhibition mechanism of two small molecules, homotaurine and scyllo-inositol, which are AD drug candidates currently under investigation. We show that both small molecules promote a conformational change of the Aβ42 monomer toward a more collapsed phase through a nonspecific binding mechanism. This finding provides atomistic-level insights into designing potential drug candidates for future AD treatments.

Molecular Simulations Reveal Terminal Group Mediated Stabilization of Helical Conformers in Both Amyloid-β42 and α-Synuclein

The presence of partially structured helices in natively unfolded amyloid-β42 (Aβ42) and α-synuclein (αS) has been shown to accelerate fibrillation in the onset of Alzheimer's and Parkinson's disease, respectively. At the other extreme, folded stable helical conformers have also been reported to resist amyloid formation. Recent studies indicate that amyloidogenic aggregation can be impeded using small molecules that stabilize the α-helical monomers and switch off the neurotoxic pathway. We predict a common intrapeptide route to stabilization based on the plasticity of helical conformations of Aβ42 and αS as assessed through extensive atomistic molecular dynamics (MD) computer simulations (∼36 μs) across ten distinct protein force field and water model combinations. Computed free energies and interaction maps (not obtainable from experiments alone) show that flexible terminal groups (N-terminus of Aβ42 and C-terminus of αS) show a tendency to stabilize folded helical conformations in both peptides via primary hydrophobic interactions with central hydrophobic domains, and secondary salt bridges with other domains. These interactions confer aggregation resistance by decreasing the population of partially structured helices and are absent in control simulations of complete unfolding. Computed helical stability is also significantly reduced in terminal-deleted variants. The models suggest new strategies to tackle neurodegeneration by rationally re-engineering terminal groups to optimize their predicted ability to deactivate helical monomers.

Adequate prediction for inhibitor affinity of Aβ40 protofibril using the linear interaction energy method

The search for efficient inhibitors targeting Aβ oligomers and fibrils is an important issue in Alzheimer's disease treatment. As a consequence, an accurate and computationally cheap approach to estimate the binding affinity for many ligands interacting with Aβ peptides is very important. Here, the calculated binding free energies of 30 ligands interacting with 12Aβ11-40 peptides using the linear interaction energy (LIE) approach are found to be in good correlation with experimental data (R = 0.79). The binding affinities of these complexes are also calculated by using free energy perturbation (FEP) and molecular mechanic/Poisson-Boltzmann surface area (MM/PBSA) methods. The time-consuming FEP method provides results with similar correlation (R = 0.72), whereas MM/PBSA calculations show very low correlation with experimental data (R = 0.27). In all complexes, van der Waals interactions contribute much more than electrostatic interactions. The LIE model, which is much less time-consuming than both the FEP and MM/PBSA methods, opens the door to accurate and rapid affinity prediction of ligands with Aβ peptides and the design of new ligands.

Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.

Looking at the Disordered Proteins through the Computational Microscope

Intrinsically disordered proteins (IDPs) have attracted wide interest over the past decade due to their surprising prevalence in the proteome and versatile roles in cell physiology and pathology. A large selection of IDPs has been identified as potential targets for therapeutic intervention. Characterizing the structure-function relationship of disordered proteins is therefore an essential but daunting task, as these proteins can adapt transient structure, necessitating a new paradigm for connecting structural disorder to function. Molecular simulation has emerged as a natural complement to experiments for atomic-level characterizations and mechanistic investigations of this intriguing class of proteins. The diverse range of length and time scales involved in IDP function requires performing simulations at multiple levels of resolution. In this Outlook, we focus on summarizing available simulation methods, along with a few interesting example applications. We also provide an outlook on how these simulation methods can be further improved in order to provide a more accurate description of IDP structure, binding, and assembly.