Effect of pH on Aβ42 Monomer and Fibril-like Oligomers-Decoding in Silico of the Roles of pK Values of Charged Residues

Aβ40 Aggregation under Changeable Conditions

Homeostasis is crucial for cell function, and disturbances in homeostasis can lead to health disorders. Under normal conditions, intracellular pH is maintained between 7.35 and 7.45. Altered endosomal and lysosomal pH together with a general drop in brain pH are associated with the aggregation of amyloid-β-peptide (Aβ) and the development of Alzheimer's disease. Under acidic conditions, close to the Aβ isoelectric point, the absence of charges favors the formation of intermolecular contacts and promotes aggregation. Here, we analyzed how pH levels affect the aggregation of Aβ40 considering the variations in brain pH and the coexistence of different aggregated conformations. Our results suggest that different macromolecular conformations can interact with each other and influence the aggregation process. In addition, we showed that neutral pH and physiological salt concentrations favor a slow aggregation, resulting in ordered, stable fibrils, with low cytotoxic effects. Overall, we highlight the complexity of the aggregation processes occurring in different physiological and pathological environments.

New Insights into the Structural and Binding Properties on Aβ Mature Fibrils Due to Histidine Protonation Behaviors

Histidine tautomeric behaviors have been considered origin factors for controlling the structure and aggregation properties of misfolding peptides. Except for tautomeric behaviors, histidine protonation behaviors definitely have the same capacities due to the net charge changes and the various N/N-H orientations on imidazole rings. However, such phenomena are still unknown. In the current study, Aβ mature fibrils substituted with various protonation states were performed by molecular dynamics simulations to investigate the structure and binding properties. Our results show that all kinds of protonation states can increase the ΔG1 stability and decrease ΔG2 and ΔG3 stabilities. A significantly higher averaged β-sheet content was detected in (εεp), (εpp), and (ppp) fibrils in one, two, and three protonation stages, respectively. Impressively, we found that the substituted fibril with specific protonated states can control the N-terminus structural properties. Further analysis confirmed that H6 and H13 are more important than H14 since the H-bond donor and receptor cooperate among C1/C3/C8_H6, C1/C3/C8_H13, and C1/C3/C8_E11. Furthermore, the mechanism of protonation behaviors was discussed. The current study is helpful for understanding the histidine protonation behaviors on one, two, and three protonation stages, which provides new horizons for exploring the origin of protein folding and misfolding.

Interactive patches over amyloid-β oligomers mediate fractal self-assembly

The monomeric units of intrinsically disordered proteins self-assemble into oligomers, protofilaments, and eventually fibrils which may turn into amyloid. The aggregation of these proteins is primarily studied in bulk with no restriction on their degrees of freedom. Herein we experimentally demonstrate that amyloid-β (Aβ) aggregation under diffusion-limited conditions leads to its fractal self-assembly. Confocal microscopy and scanning electron microscopy with energy dispersion x-ray analysis were used to confirm that the fractal self-assemblies were formed from Aβ rather than the salt present in the two supporting media: deionized water and phosphate buffered saline. The results from the molecular docking experiments implicated that electrostatic and hydrophobic patches on the solvent-accessible surface area of the Aβ oligomers mediate the fractal self-assembly. These implications were tested with laser light scattering experiments on the oligomers formed by breaking mature fibrils of Aβ through sonication, which were observed to self-assemble into fractals when sonicated solutions were drop casted. The electrostatic interactions modulate the fractal morphologies with pH of the solution, which leads to a morphological phase transition observed through the variation in their fractal dimension. These transitions provide experimental evidence for the existing theoretical framework in terms of different kinetic models. The higher surface-to-volume ratio of these fractal self-assemblies may have applications in drug delivery, biosensing, and other biomedical applications.

Arginine and Arginine-Rich Peptides as Modulators of Protein Aggregation and Cytotoxicity Associated With Alzheimer's Disease

A substantial body of evidence indicates cationic, arginine-rich peptides (CARPs) are effective therapeutic compounds for a range of neurodegenerative pathologies, with beneficial effects including the reduction of excitotoxic cell death and mitochondrial dysfunction. CARPs, therefore, represent an emergent class of promising neurotherapeutics with multimodal mechanisms of action. Arginine itself is a known chaotrope, able to prevent misfolding and aggregation of proteins. The putative role of proteopathies in chronic neurodegenerative diseases such as Alzheimer's disease (AD) warrants investigation into whether CARPs could also prevent the aggregation and cytotoxicity of amyloidogenic proteins, particularly amyloid-beta and tau. While monomeric arginine is well-established as an inhibitor of protein aggregation in solution, no studies have comprehensively discussed the anti-aggregatory properties of arginine and CARPs on proteins associated with neurodegenerative disease. Here, we review the structural, physicochemical, and self-associative properties of arginine and the guanidinium moiety, to explore the mechanisms underlying the modulation of protein aggregation by monomeric and multimeric arginine molecules. Arginine-rich peptide-based inhibitors of amyloid-beta and tau aggregation are discussed, as well as further modulatory roles which could reduce proteopathic cytotoxicity, in the context of therapeutic development for AD.

Destabilization of the Alzheimer's amyloid-β protofibrils by THC: A molecular dynamics simulation study

Alzheimer's disease is a leading cause of dementia in the elderly population for which there is no cure at present. Deposits of neurotoxic plaques are found in the brains of patients which are composed of fibrils of the amyloid-β peptide. Molecules which can disrupt these fibrils have gained attention as potential therapeutic agents. Δ-tetrahydrocannabidiol (THC) is a cannabinoid, which can bind to the receptors in the brain, and has shown promise in reducing the fibril content in many experimental studies. In our present study, by employing all atom molecular dynamics simulations, we have investigated the mechanism of the interaction of the THC molecules with the amyloid-β protofibrils. Our results show that the THC molecules disrupt the protofibril structure by binding strongly to them. The driving force for the binding was the hydrophobic interactions with the hydrophobic residues in the fibrils. As a result of these interactions, the tight packing of the hydrophobic core of the protofibrils was made loose, and salt bridges, which were important for stability were disrupted. Hydrogen bonds between the chains of the protofibrils which are important for stability were disrupted, as a result of which the β-sheet content was reduced. The destabilization of the protofibrils by the THC molecules leads to the conclusion that THC molecules may be considered for the therapy in treating Alzheimer's disease.

Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease

The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b, the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.

Dual promoter strategy enhances co-expression of α-L-rhamnosidase and enhanced fluorescent protein for whole-cell catalysis and bioresource valorization

Developing circular economy is the only way to improve the efficiency of resource utilization. Whole-cell catalysis is an effective method to recycle enzymes, improve catalytic efficiency, and reduce production costs. The enzyme, α-L-rhamnosidase has considerable application prospects in the field of biocatalysis as it can hydrolyze a variety of α-L rhamnoses. In the present study, the genes for α-L-rhamnosidase (rhaB1) and enhanced fluorescent protein (EGFP) were co-expressed using a bi-promoter expression vector pRSFDuet1 and their enzymatic properties were evaluated. To our knowledge, this study has established an effective rhamnosidase-fluorescent indicator and whole-cell catalytic system for the first time. Moreover, we analyzed the change in the activity of the crude rhaB1-EGFP as well as its whole-cell during the biocatalysis process using fluorescence intensity. Recombinant rhaB1-EGFP as a product which contains rhaB1 and EGFP showed higher thermal stability, pH stability, and conversion efficiency than rhaB1, and its optimum temperature for rutin catalysis was ideal for industrial applications. Moreover, under the optimal conditions of a rutin concentration of 0.05 g/L, pH of 6.0, temperature of 40 °C, a yield of 92.5% was obtained. Furthermore, we demonstrated the relationship between the fluorescence intensity and enzyme activity. This study established a highly efficient whole-cell catalytic system whose activity can be evaluated by fluorescence intensity, providing a reference for enzyme recycling.

Novel Disassembly Mechanisms of Sigmoid Aβ42 Protofibrils by Introduced Neutral and Charged Drug Molecules

Alzheimer's disease (AD) is characterized by fibrillar deposits of amyloid-β (Aβ) peptides and neurofibrillary tangles of Tau proteins. Aβ peptides are composed of 37-49 residues, among which the Aβ₄₂ isoform is particularly toxic and aggregation-prone and is enriched in the plaques of AD brains and thus considered central to the development of AD. Therefore, disaggregation and disruption provide potential therapeutic approaches to reduce, inhibit, and even reverse Aβ aggregation. Here we capture the atomic-level details of the interactions between sigmoid Aβ₄₂ fibril 2MXU or 5KK3 and either natural tanshinone compounds TS1 or TS0 or negatively charged ER, proposing two unprecedented disassembly mechanisms. Natural TS1 or TS0 prefers to insert into the cavity together with part at the surface of the 2MXU to open up the mouth and twist the conformation, destroying the ordered growth of subsequent monomers along the fibril axis. For the more compact two-fold 5KK3 , attachment of TS1 or TS0 at the surface including some inserted in cavity results in the separation of the two folds. In the two sigmoid fibril systems, it is no longer applicable for the routine criteria to assess Aβ₄₂ fibril disassembly by introduction of these drugs, such as either reduced H-bond number, decreased β-sheet contents, or both. ER, like-charged to Aβ₄₂ fibril, is especially exceptional, and departs utterly from the neutral ones to disassemble Aβ₄₂ fibril. Besides the inapplicable routine criteria, positive binding energy between ER and Aβ₄₂ fibril also deviates from the hypotheses of "ligands exhibiting greater affinity for the β-amyloid peptide are effective at altering its aggregation and inhibiting cell toxicity" ( Cairo et al. , Biochemistry 2002 , 41 , 8620 - 8629 ) but results in stronger disassembly effect on the two kinds of sigmoid Aβ₄₂ fibrils than neutral TS0 or TS1. The disassembly power of charged ER molecules derives from its stronger deformation ability to the conformation of Aβ₄₂ fibril than the neutral ones, twisting the one-fold 2MXU into tapered-shape and separating two-fold 5KK3 in two parts further, which is in great agreement with experimental observations ( Irwin et al. Biomacromolecules 2013 , 14 ( 1 ), 264 - 274 ). The unusual disassembly mechanisms fill the gaps and offer an alternative direction in engineering new inhibitors to treat AD.

Folding dynamics of Aβ42 monomer at pH 4.0–7.5 with and without physiological salt conditions – does the β1 or β2 region fold first?

The seeding region of Aβ42 monomer is jointly affected by the solution acidity, ionic distribution of the salt, and charged residues.

Amyloid beta oligomers: how pH influences over trimer and pentamer structures?

The aggregation of proteins in the brain is one of the main features of neurodegenerative diseases. In Alzheimer's disease, the abnormal aggregation of Aβ-42 is due to intrinsic and extrinsic factors. The latter is due to variations in the environment, such as temperature, salt concentration, and pH. We evaluated the effect of protonation/deprotonation of residues that are part of trimeric and pentameric oligomers at pH 5, pH 6, and pH 7. Molecular dynamics simulation at 200 ns in the canonical ensemble was implemented. The results have revealed that histidine, glutamic acid, and aspartic acid residues showed a protonation/deprotonation effect in oligomers. The root mean square deviation analysis was used to analyze the structural stability at different pHs. We found an increase in hydrophobicity in the side chains of the trimer, while in the pentamer, the structural instability of a compact structure at pH 5 caused the hydrophobic core to open, revealing the hydrophobic region to the environment. At this point, we believe that conformational changes mediated by pH are essential in the aggregation of Aβ-42 oligomers.

C-Terminal Plays as the Possible Nucleation of the Self-Aggregation of the S-Shape Aβ11-42 Tetramer in Solution: Intensive MD Study

Amyloid beta (Aβ) peptides are characterized as the major factors associated with neuron death in Alzheimer's disease, which is listed as the most common form of neurodegeneration. Disordered Aβ peptides are released from proteolysis of the amyloid precursor protein. The Aβ self-assembly process roughly takes place via five steps: disordered forms → oligomers → photofibrils → mature fibrils → plaques. Although Aβ fibrils are often observed in patient brains, oligomers were recently indicated to be major neurotoxic elements. In this work, the neurotoxic compound S-shape Aβ11-42 tetramer (S4Aβ11-42) was investigated over 10 μs of unbiased MD simulations. In particular, the S4Aβ11-42 oligomer adopted a high dynamics structure, resulting in unsuccessful determination of their structures in experiments. The C-terminal was suggested as the possible nucleation of the Aβ42 aggregation. The sequences 27-35 and 39-40 formed rich β-content, whereas other residues mostly adopted coil structures. The mean value of the β-content over the equilibrium interval is ∼42 ± 3%. Furthermore, the dissociation free energy of the S4Aβ11-42 peptide was predicted using a biased sampling method. The obtained free energy is ΔG US = -58.44 kcal/mol which is roughly the same level as the corresponding value of the U-shape Aβ17-42 peptide. We anticipate that the obtained S4Aβ11-42 structures could be used as targets for AD inhibitor screening over the in silico study.

The influences of E22Q mutant on solvated 3Aβ11-40 peptide: A REMD study

The residue E22 plays a critical role in the aggregation process of Amyloid beta (Aβ) peptides. The effect of E22Q mutant on the shapes of the solvated Aβ_11-40 trimer is clarified using a replica exchange molecular dynamics (REMD) simulation employing ∼20.6 μs of MD simulations with 48 disparate replicas. The increase of intramolecular polar contacts and salt bridge between the residue D23 to residues (24-29) was observed. The residual secondary structure of the mutated trimer is shifted in a similar way to the picture observed in previous investigations of F19W mutant. The free energy surface (FES) of the mutated E22Q system has a fewer number of minima in comparison with the wild-type trimer. The optimized shapes of the mutated E22Q form a significant increase in beta structure (47%) and serious decrease in coil content (46%) compared with the wild-type (of 36 and 56%, respectively). The binding affinity of constituting chains to the rest is of -43.7 ± 6.5 kcal/mol, implying that the representative structure of E22Q is more stable than the wild-type one. Furthermore, the E22Q mutant increases the size of stable structures due to larger collision cross section (CCS) and solvent accessible area (SASA). The observed results may enhance the Aβ inhibition throughout the contribution to the knowledge of the Aβ oligomerization/aggregation.