Identifying the Template for Oligomer to Fibril Conversion for Amyloid-β (1-42) Oligomers using Hamiltonian Replica Exchange Molecular Dynamics

The toxicity of amyloid-β (Aβ) oligomers has been known to be higher compared to mature fibrils. Yet the presence of plaques in Alzheimer's disease patients indicates the significance of oligomer to fibril conversion for Aβ aggregates. In this study, we investigate Aβ_13-42 oligomers having two to five peptide chains using extensive all-atom molecular dynamics simulations to identify the on- or off-pathway intermediates in fibril formation pathway. Hamiltonian replica exchange method through solute tempering (REST2) has been employed to explore the different structures attained by these aggregates. Using intra-chain and inter-chain contacts as reaction coordinates, we obtain the free energy surface for the Aβ_13-42 oligomers. Consequently, their stable conformations and structural features have been identified. The found conformations belonging to most probable structures possess both parallel and anti-parallel β-sheets, characteristic of on- and off-pathway intermediates, respectively. Further, we have measured the tendency to form fibril like interactions among the β-sheet forming residues. Our analysis finds that residues 30-36 possess higher tendency to form fibril like contacts among all the residues. While we find stronger interaction among residues 30-36, these amino acids are also found to be more shielded from water compared to others. With previous experimental studies finding these residues to be more crucial for the stability of Aβ42 oligomers, we propose that interactions within this patch could trigger seed formation that leads to conversion of on-pathway oligomers into disease relevant fibrils.

Interactions of the Aβ(1-42) Peptide with Boron Nitride Nanoparticles of Varying Curvature in an Aqueous Medium: Different Pathways to Inhibit β-Sheet Formation

The aggregation of amyloid β (Aβ) peptide triggered by its conformational changes leads to the commonly known neurodegenerative disease of Alzheimer's. It is believed that the formation of β sheets of the peptide plays a key role in its aggregation and subsequent fibrillization. In the current study, we have investigated the interactions of the Aβ(1-42) peptide with boron nitride nanoparticles and the effects of the latter on conformational transitions of the peptide through a series of molecular dynamics simulations. In particular, the effects of curvature of the nanoparticle surface are studied by considering boron nitride nanotubes (BNNTs) of varying diameter and also a planar boron nitride nanosheet (BNNS). Altogether, the current study involves the generation and analysis of 9.5 μs of dynamical trajectories of peptide-BNNT/BNNS pairs in an aqueous medium. It is found that BN nanoparticles of different curvatures that are studied in the present work inhibit the conformational transition of the peptide to its β-sheet form. However, such an inhibition effect follows different pathways for BN nanoparticles of different curvatures. For the BNNT with the highest surface curvature, i.e., (3,3) BNNT, the nanoparticle is found to inhibit β-sheet formation by stabilizing the helical structure of the peptide, whereas for planar BNNS, the β-sheet formation is prevented by making more favorable pathways available for transitions of the peptide to conformations of random coils and turns. The BNNTs with intermediate curvatures are found to exhibit diverse pathways of their interactions with the peptide, but in all cases, essentially no formation of the β sheet is found whereas substantial β-sheet formation is observed for Aβ(1-42) in water in the absence of any nanoparticle. The current study shows that BN nanoparticles have the potential to act as effective tools to prevent amyloid formation from Aβ peptides.

The Pathogenesis Mechanism, Structure Properties, Potential Drugs and Therapeutic Nanoparticles against the Small Oligomers of Amyloid-β

Alzheimer's Disease (AD) is a devastating neurodegenerative disease that affects millions of people in the world. The abnormal aggregation of amyloid β protein (Aβ) is regarded as the key event in AD onset. Meanwhile, the Aβ oligomers are believed to be the most toxic species of Aβ. Recent studies show that the Aβ dimers, which are the smallest form of Aβ oligomers, also have the neurotoxicity in the absence of other oligomers in physiological conditions. In this review, we focus on the pathogenesis, structure and potential therapeutic molecules against small Aβ oligomers, as well as the nanoparticles (NPs) in the treatment of AD. In this review, we firstly focus on the pathogenic mechanism of Aβ oligomers, especially the Aβ dimers. The toxicity of Aβ dimer or oligomers, which attributes to the interactions with various receptors and the disruption of membrane or intracellular environments, were introduced. Then the structure properties of Aβ dimers and oligomers are summarized. Although some structural information such as the secondary structure content is characterized by experimental technologies, detailed structures are still absent. Following that, the small molecules targeting Aβ dimers or oligomers are collected; nevertheless, all of these ligands have failed to come into the market due to the rising controversy of the Aβ-related "amyloid cascade hypothesis". At last, the recent progress about the nanoparticles as the potential drugs or the drug delivery for the Aβ oligomers are present.

Benchmark of Popular Free Energy Approaches Revealing the Inhibitors Binding to SARS-CoV-2 Mpro

The COVID-19 pandemic has killed millions of people worldwide since its outbreak in December 2019. The pandemic is caused by the SARS-CoV-2 virus whose main protease (Mpro) is a promising drug target since it plays a key role in viral proliferation and replication. Currently, developing an effective therapy is an urgent task, which requires accurately estimating the ligand-binding free energy to SARS-CoV-2 Mpro. However, it should be noted that the accuracy of a free energy method probably depends on the protein target. A highly accurate approach for some targets may fail to produce a reasonable correlation with the experiment when a novel enzyme is considered as a drug target. Therefore, in this context, the ligand-binding affinity to SARS-CoV-2 Mpro was calculated via various approaches. The molecular docking approach was manipulated using Autodock Vina (Vina) and Autodock4 (AD4) protocols to preliminarily investigate the ligand-binding affinity and pose to SARS-CoV-2 Mpro. The binding free energy was then refined using the fast pulling of ligand (FPL), linear interaction energy (LIE), molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA), and free energy perturbation (FEP) methods. The benchmark results indicated that for docking calculations, Vina is more accurate than AD4, and for free energy methods, FEP is the most accurate method, followed by LIE, FPL, and MM-PBSA (FEP > LIE ≈ FPL > MM-PBSA). Moreover, atomistic simulations revealed that the van der Waals interaction is the dominant factor. The residues Thr26, His41, Ser46, Asn142, Gly143, Cys145, His164, Glu166, and Gln189 are essential elements affecting the binding process. Our benchmark provides guidelines for further investigations using computational approaches.

The F19W mutation reduces the binding affinity of the transmembrane Aβ11-40 trimer to the membrane bilayer

Alzheimer's disease is linked to the aggregation of the amyloid-β protein (Aβ) of 40 or 42 amino acids. Lipid membranes are known to modulate the rate and mechanisms of the Aβ aggregation. Point mutations in Aβ can alter these rates and mechanisms. In particular, experiments show that F19 mutations influence the aggregation rate, but maintain the fibril structures. Here, we used molecular dynamics simulations to examine the effect of the F19W mutation in the 3Aβ11-40 trimer immersed in DPPC lipid bilayers submerged in aqueous solution. Substituting Phe by its closest (non-polar) aromatic amino acid Trp has a dramatic reduction in binding affinity to the phospholipid membrane (measured with respect to the solvated protein) compared to the wild type: the binding free energy of the protein-DPPC lipid bilayer increases by 40-50 kcal mol-1 over the wild-type. This is accompanied by conformational changes and loss of salt bridges, as well as a more complex free energy surface, all indicative of a more flexible and less stable mutated trimer. These results suggest that the impact of mutations can be assessed, at least partially, by evaluating the interaction of the mutated peptides with the lipid membranes.

Thermodynamics and kinetics in antibody resistance of the 501Y.V2 SARS-CoV-2 variant

Understanding the thermodynamics and kinetics of the binding process of an antibody to the SARS-CoV-2 receptor-binding domain (RBD) of the spike protein is very important for the development of COVID-19 vaccines. In particular, it is essential to understand how the binding mechanism may change under the effects of RBD mutations. In this context, we have demonstrated that the South African variant (B1.351 or 501Y.V2) can resist the neutralizing antibody (NAb). Three substitutions in the RBD including K417N, E484K, and N501Y alter the free energy landscape, binding pose, binding free energy, binding kinetics, hydrogen bonding, nonbonded contacts, and unbinding pathway of RBD + NAb complexes. The low binding affinity of NAb to 501Y.V2 RBD confirms the antibody resistance of the South African variant. Moreover, the fragment of NAb + RBD can be used as an affordable model to investigate changes in the binding process between the mutated RBD and antibodies.

Increasing FEL minima of 501Y.V2 RBD + antibody in comparison with the WT RBD systems imply that the complex 501Y.V2 RBD + antibody is more unstable than the WT one.

Emergence of Barrel Motif in Amyloid-β Trimer: A Computational Study

Amyloid-β (Aβ) peptides form assemblies that are pathological hallmarks of Alzheimer's disease. Aβ oligomers are soluble, mobile, and toxic forms of the peptide that act in the extracellular space before assembling into protofibrils and fibrils. Therefore, oligomers play an important role in the mechanism of Alzheimer's disease. Since it is difficult to determine by experiment the atomic structures of oligomers, which accumulate fast and are polymorphic, computer simulation is a useful tool to investigate elusive oligomers' structures. In this work, we report extended all-atom molecular dynamics simulations, both canonical and replica exchange, of Aβ(1-42) trimer starting from two different initial conformations: (i) the pose produced by the best docking of a monomer aside of a dimer (simulation 1), representing oligomers freshly formed by assembling monomers, and (ii) a configuration extracted from an experimental mature fibril structure (simulation 2), representing settled oligomers in equilibrium with extended fibrils. We showed that in simulation 1, regions with small β-barrels are populated, indicating the chance of spontaneous formation of domains resembling channel-like structures. These structural domains are alternative to those more representative of mature fibrils (simulation 2), the latter showing a stable bundle of C-termini that is not sampled in simulation 1. Moreover, trimer of Aβ(1-42) can form internal pores that are large enough to be accessed by water molecules and Ca2+ ions.

On the Conformational Dynamics of β-Amyloid Forming Peptides: A Computational Perspective

Understanding the conformational dynamics of proteins and peptides involved in important functions is still a difficult task in computational structural biology. Because such conformational transitions in β-amyloid (Aβ) forming peptides play a crucial role in many neurological disorders, researchers from different scientific fields have been trying to address issues related to the folding of Aβ forming peptides together. Many theoretical models have been proposed in the recent years for studying Aβ peptides using mathematical, physicochemical, and molecular dynamics simulation, and machine learning approaches. In this article, we have comprehensively reviewed the developmental advances in the theoretical models for Aβ peptide folding and interactions, particularly in the context of neurological disorders. Furthermore, we have extensively reviewed the advances in molecular dynamics simulation as a tool used for studying the conversions between polymorphic amyloid forms and applications of using machine learning approaches in predicting Aβ peptides and aggregation-prone regions in proteins. We have also provided details on the theoretical advances in the study of Aβ peptides, which would enhance our understanding of these peptides at the molecular level and eventually lead to the development of targeted therapies for certain acute neurological disorders such as Alzheimer's disease in the future.

Impact of the Astaxanthin, Betanin, and EGCG Compounds on Small Oligomers of Amyloid Aβ40 Peptide

There is experimental evidence that the astaxanthin, betanin, and epigallocatechin-3-gallate (EGCG) compounds slow down the aggregation kinetics and the toxicity of the amyloid-β (Aβ) peptide. How these inhibitors affect the self-assembly at the atomic level remains elusive. To address this issue, we have performed for each ligand atomistic replica exchange molecular dynamic (REMD) simulations in an explicit solvent of the Aβ_11-40 trimer from the U-shape conformation and MD simulations starting from Aβ_1-40 dimer and tetramer structures characterized by different intra- and interpeptide conformations. We find that the three ligands have similar binding free energies on small Aβ₄₀ oligomers but very distinct transient binding sites that will affect the aggregation of larger assemblies and fibril elongation of the Aβ₄₀ peptide.

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.

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.

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.

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.

In vitroandin silicodetermination of glutaminyl cyclase inhibitors

Alzheimer's disease (AD) is the most common form of neurodegenerative disease currently. It is widely accepted that AD is characterized by the self-assembly of amyloid beta (Aβ) peptides. The human glutaminyl cyclase (hQC) enzyme is characterized by association with Aβ peptide generation. The development of hQC inhibitors could prevent the self-aggregation of Aβ peptides, resulting in impeding AD. Utilizing structural knowledge of the hQC substrates and known hQC inhibitors, new heterocyclic and peptidomimetic derivatives were synthesized and were able to inhibit the hQC enzyme. The inhibiting abilities of these compounds were evaluated using a fluorometric assay. The binding mechanism at the atomic level was estimated using molecular docking, free energy perturbation, and quantum chemical calculation methods. The predicted log(BBB) and human intestinal absorption values indicated that these compounds are able to permeate the blood–brain barrier and be well-absorbed through the gastrointestinal tract. Overall, 5,6-dimethoxy-N-(3-(5-methyl-1H-imidazol-1-yl)propyl)-1H-benzo[d]imidazol-2-amine (1_2) was indicated as a potential drug for AD treatment.

Rational design of new hQC inhibitors.

Interaction of carbohydrate binding module 20 with starch substrates

CBM20s are starch-binding domains found in many amylolytic enzymes, including glucoamylase, alpha-amylase, beta-amylases, and a new family of starch-active polysaccharide monooxygenases (AA13 PMOs). Previous studies of CBM20–substrate interaction only concerned relatively small or soluble amylose molecules, while amylolytic enzymes often work on extended chains of insoluble starch molecules. In this study, we utilized molecular simulation techniques to gain further insights into the interaction of CBM20 with substrates of various sizes via its two separate binding sites, termed as BdS1 and BdS2. Results show that substrate binding at BdS1 involving two conserved tryptophan residues is about 2–4 kcal mol⁻¹ stronger than that at BdS2. CBM20 exhibits about two-fold higher affinity for helical substrates than for the amylose random coils. The affinity for amylose individual double helices does not depend on the helices' length. At least three parallel double helices are required for optimal binding. The binding affinity for a substrate containing 3 or more double helices is ∼−15 kcal mol⁻¹, which is 2–3 kcal mol⁻¹ larger than that for individual double helices. 100 ns molecular dynamics simulations were carried out for the binding of CBM20 to an extended substrate containing 3 layers of 9 60-unit double helices (A3L). A stable conformation of CBM20–A3L was found at BdS1. However, when CBM20 binds A3L viaBdS2, it moves across the surface of the substrate and does not form a stable complex. MD simulations show that small amylose helices are quickly disrupted upon binding to CBM20. Our results provide some important molecular insights into the interactions of CBM20 with starch substrates, which will serve as the basis for further studies of CBM20-containing enzymes, including AA13 PMOs.

CBM20 quickly disrupts small helical amylose substrates and exhibits optimal binding affinity when the substrate has three or more parallel double helices.

Etersalate prevents the formations of 6Aβ16-22 oligomer: An in silico study

Oligomerization of amyloid beta (Aβ) peptides has been considered as the crucially causative agent in the development of Alzheimer's disease. Etersalate, a nonsteroidal anti-inflammatory oral drug (United State Food and Drug Administration—Unique Ingredient Identifier: 653GN04T2G) was previously suggested to bind well to proto-fibrils of Aβ peptides in silico. Here, the effect of etersalate on the oligomerization of soluble Aβ_16–22 hexamer (6Aβ_16–22) were extensively investigated using temperature replica exchange molecular dynamics (REMD) simulations over ~16.8 μs in total for 48 replicas (350 ns per replica). The results reveal that etersalate can enter the inner space or bind on the surface of 6Aβ_16–22 conformations, which destabilizes the hexamer. Etersalate was predicted to able to cross the blood brain barrier using prediction of absorption, distribution, metabolism, and excretion—toxicity (preADMET) tools. Overall, although the investigation was performed with the low concentration of trial inhibitor, the obtained results indicate that etersalate is a potential drug candidate for AD through inhibiting formation of Aβ oligomers with the average binding free energy of -11.7 kcal/mol.

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.

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

Amyloid beta-peptide (Aβ) is the key to developing Alzheimer's disease. Experiments have confirmed that different acidity influences directly not only the structural morphology and population of Aβ oligomers, but also the toxicity. The atomic-level association between the pH, charged residues, and Aβ properties remains obscure. Herein, conformational changes of Aβ₄₂ monomer, fibril-like trimer, and pentamer in the medium pH range of 4.0-7.5 are studied. The results reveal that, as the pH changes from 7.5 to the isoelectric pH, His6, His13, and His14 are protonated in turn, successively approach the center of mass of folded Aβ monomer, trigger ionic interactions and changes of neighboring turns (Asp7-Ser8, His14-Lys16) and even a distant one (Leu34-Met35), as well as concomitant changes of secondary structure, and promote the conformation transition from unfolded to folded. This observation discloses that protonation can convert these charged residues from originally hydrophilic to "hydrophobic-like". For fibril-like oligomers, the pH susceptibility essentially stems from the pK values of charged residues in the context of the Aβ fibril, and in turn one can predict the dynamic behavior of these residues in the processes of dissociation or stabilization of a fibril by comparing the pK values of residues involved in salt bridges in the normal state with those in the current context. This idea is justified by two fibril models and appears to be applicable to other peptides and their fibril systems.

Conformational Dynamics and Stability of U-Shaped and S-Shaped Amyloid β Assemblies

Alzheimer’s disease is the most fatal neurodegenerative disorder characterized by the aggregation and deposition of Amyloid β (Aβ) oligomers in the brain of patients. Two principal variants of Aβ exist in humans: Aβ_1–40 and Aβ_1–42. The former is the most abundant in the plaques, while the latter is the most toxic species and forms fibrils more rapidly. Interestingly, fibrils of Aβ_1–40 peptides can only assume U-shaped conformations while Aβ_1–42 can also arrange as S-shaped three-stranded chains, as recently discovered. As alterations in protein conformational arrangement correlate with cell toxicity and speed of disease progression, it is important to characterize, at molecular level, the conformational dynamics of amyloid fibrils. In this work, Replica Exchange Molecular Dynamics simulations were carried out to compare the conformational dynamics of U-shaped and S-shaped Aβ_17–42 small fibrils. Our computational results provide support for the stability of the recently proposed S-shaped model due to the maximized interactions involving the C-terminal residues. On the other hand, the U-shaped motif is characterized by significant distortions resulting in a more disordered assembly. Outcomes of our work suggest that the molecular architecture of the protein aggregates might play a pivotal role in formation and conformational stability of the resulting fibrils.

Atomistic investigation of an Iowa Amyloid-β trimer in aqueous solution

The self-assembly of Amyloid beta (Aβ) peptides are widely accepted to associate with Alzheimer's disease (AD) via several proposed mechanisms. Because Aβ oligomers exist in a complicated environment consisting of various forms of Aβ, including oligomers, protofibrils, and fibrils, their structure has not been well understood. The negatively charged residue D23 is one of the critical residues of the Aβ peptide as it is located in the central hydrophobic domain of the Aβ N-terminal and forms a salt-bridge D23-K28, which helps stabilize the loop domain. In the familial Iowa (D23N) mutant, the total net charge of Aβ oligomers decreases, resulting in the decrease of electrostatic repulsion between D23N Aβ monomers and thus the increase in their self-aggregation rate. In this work, the impact of the D23N mutation on 3Aβ_11–40 trimer was characterized utilizing temperature replica exchange molecular dynamics (REMD) simulations. Our simulation reveals that D23N mutation significantly enhances the affinity between the constituting chains in the trimer, increases the β-content (especially in the sequence 21–23), and shifts the β-strand hydrophobic core from crossing arrangement to parallel arrangement, which is consistent with the increase in self-aggregation rate. Molecular docking indicates that the Aβ fibril-binding ligands bind to the D23N and WT forms at different poses. These compounds prefer to bind to the N-terminal β-strand of the D23N mutant trimer, while they mostly bind to the N-terminal loop region of the WT. It is important to take into account the difference in the binding of ligands to mutant and wild type Aβ peptides in designing efficient inhibitors for various types of AD.

Amyloid beta peptide oligomers are believed to play key roles in Alzheimer's disease pathogenesis. D23N mutation significantly changes their structure and how they bind potential inhibitors.

Evaluation of the absolute affinity of neuraminidase inhibitor using steered molecular dynamics simulations

The absolute free energy difference of binding (ΔG) between neuraminidase and its inhibitor was evaluated using fast pulling of ligand (FPL) method over steered molecular dynamics (SMD) simulations. The metric was computed through linear interaction approximation. Binding nature was described by free energy differences of electrostatic and van der Waals (vdW) interactions. The finding indicates that vdW metric is dominant over electrostatics in binding process. The computed values are in good agreement with experimental data with a correlation coefficient of R=0.82 and error of σΔG_exp=2.2kcal/mol. The results were observed using Amber99SB-ILDN force field in comparison with CHARMM27 and GROMOS96 43a1 force fields. Obtained results may stimulate the search for an Influenza therapy.

EGCG inhibits the oligomerization of amyloid beta (16-22) hexamer: Theoretical studies

An extensive replica exchange molecular dynamics (REMD) simulation was performed to investigate the progress patterns of the inhibition of (-)-epigallocatechin-3-gallate (EGCG) on the Aβ_16-22 hexamer. Structural variations of the oligomers without and with EGCG were monitored and analyzed in detail. It has been found that EGCG prevents the formation of Aβ oligomer through two different ways by either accelerating the Aβ oligomerization or reducing the β-content of the hexamer. It also decreases the potential "highly toxic" conformations of Aβ oligomer, which is related to the conformations having high order β-sheet sizes. Both electrostatic and van der Waals interaction energies are found to be involved to the binding process. Computed results using quantum chemical methods show that the π-π stacking is a critical factor of the interaction between EGCG and the peptides. As a result, the binding free energy of the EGCG to the Aβ peptides is slightly larger than that of the curcumin.

Replica exchange molecular dynamics study of the amyloid beta (11–40) trimer penetrating a membrane

The transmembrane Aβ_11–40 trimer is investigated for the first time using REMD and FEP.

In silico studies of solvated F19W amyloid β (11–40) trimer

REMD studies shows that F19W mutation does not change in the overall structure of Aβ_11–40 trimer significantly but increases it flexibility, consistent with the observed formation of the same fibril structures at slower rates.