Received 4 April 2016; accepted 26 June 2016; published 29 June 2016
Alzheimer’s Disease (AD), a progressive neurodegenerative disease, is the most common cause of dementia in elderly people   . Clinically, it is characterized by loss of memory, inability to learn new things, loss of language function, a disturbed perception of space, inability to do calculations, and a host of other manifestations  . In AD, there is an overall shrinkage of brain tissue, whereas, it is microscopically characterized by extracellular (neuritic plaques) and/or intracellular (neurofibrillary tangles) deposition of insoluble amyloid beta peptides (Aβs)  . Aβ is a peptide with 42 or 43 amino acids that forms a part of the transmembrane domain of the larger Amyloid Precursor Protein (APP), a transmembrane protein expressed in neurons and other brain cells and is derived from the cleavage of APP by the β- and γ-secretase enzymes  -  . Defective clearance of Aβ from aberrant cleavage of APP and other mechanisms results in its accumulation   . Aβ monomers initially polymerize into soluble oligomers and subsequently into larger insoluble fragments which precipitate as amyloid fibrils  . The beta amyloid hypothesis suggests that AD is caused by the deposition of Aβ in plaques in brain tissue and is the basis of a novel prevention and treatment method for AD  . Indeed, AD animal models can be produced by brain ventricular infusion of Aβ1-42   . Moreover, other fragments of Aβs, such as Aβ1-40  -  and Aβ25-35  , have also been shown to assemble into amyloid fibers in vitro   and have been infused into rat brain ventricles to create AD animal models. Because AD has been shown to be associated with decreased Docosahexaenoic Acid (DHA) levels   , numerous human trials are currently ongoing to determine whether DHA is effective in the treatment and/or delaying the symptoms of AD   . We also reported that DHA decreases amyloid burden in the brains of Aβ-infused model rats   and inhibits in vitro amyloid fibrillation   . These investigations highlight the potential for DHA as an excellent therapeutic agent against Aβ-induced neurodegenerative diseases, including AD.
Previously, we showed that DHA inhibits Aβ fibril formation   -   . In this study, we investigated the computational binding of DHA to Aβ peptides to further support the anti-Aβ fibrillation effect of DHA. Because X-ray crystallographic structures of Aβs are not available, NMR structures were utilized in the present investigation. We performed molecular docking studies, which are useful for identifying agents capable of inhibiting proteins responsible for disease pathology and exploring their possible binding modes. We used thioflavin-T (ThT) as a positive control because it is currently used as a “standard” for selectively binding, staining, and identifying amyloid fibrils both in vivo and in vitro  . We docked DHA and ThT (as ligands) onto Aβ1-42 (as receptors) and compared their binding modes. We aimed to confirm the results of our previous studies, showing the ameliorative effects of DHA and its potential as an AD therapy, in this computational analysis   -   .
2. Methods and Materials
2.1. Preparation of Thioflavin T and Docosahexaenoic Acid Models
DHA (CID: 445580) and ThT (CID: 16953) were downloaded from the PubChem database in sdf format. The sdf files were submitted to the PRODRG server  to give energy-minimized structures of the ligands. Energy minimization is performed to help the docking program identify the bioactive conformer from the local minima. The two-dimensional structures of ligands are illustrated as ball and stick in Figure 1.
2.2. Analyses of Drug-Like Properties of the Ligands
Drug design essentially focuses on optimizing the binding interactions of ligands with their targets. However, a compound with the best binding interactions does not necessarily make the best drug, as other factors are involved. For example, a clinically useful drug must travel through the body to reach its target. Perfecting a compound with good drug-target interactions is irrelevant if it has no chance of reaching its target. The factors that
Figure 1. The ligand molecules thioflavin T (ThT) and docosahexaenoic acid (DHA, C22:6, n-3).
determine whether a drug will reach its target in the body are termed the pharmacokinetic or ADME (Absorp- tion, Distribution, Metabolism, and Excretion) properties. We used the free ADME/tox filtering tool (FAF-Drugs 3), which can predict whether the physicochemical properties of compounds are acceptable by applying several filtering rules, including the well-known Lipinski’s rule of five (i.e., molecular weight > 500, logP or octanol/ water partition coefficient < 5, number of H-bond donors < 5, and number of H-bond acceptors < 10) [Lipinski et al., 2001]  . The main properties computed by FAF-Drugs 3 are the number of rigid and flexible bonds, topological polar surface area (TPSA) value according to Ertl et al. (2000)  , number and maximum size of system rings, and presence of unwanted chemicals or chemical substructures (using SMARTS searches)  .
2.3. Analysis of the Amyloidogenic Regions in the Aβ1-42 and Cross-Beta Aggregation
Sequence-based computational methods, namely FoldAmyloid  , AGGRESCAN  , and ProA  , were used to predict the amyloidogenic regions in Aβ1-42. Furthermore, analysis of the beta-aggregation propensity of Aβ1-42 was achieved by subjecting the primary amino acid sequence of Aβ1-42 to the TANGO algorithm  -  , which is designed to predict cross-beta aggregation in peptides.
2.4. Preparation of Aβ1-42 Receptor Models
The three-dimensional solution structures (3D NMR) of Aβ1-42 (PDB IDs: 2BEG and 1Z0Q) were downloaded from the RCSB Protein Data Bank (PDB) (http://www.rcsb.org/) as receptors for ThT and DHA docking. 2BEG is a homopentamer, namely composed of the A, B, C, D, and E monomers of Aβ1-42. Each monomer (A, B, C, D, and E) of 2BEG comprises 10 coordinate models  . The coordinates of model 1 of the A monomer (A1), model 1 of the B monomer (B1), and model 1 of the C monomer (C1) were loaded in the Molegro virtual docker (MVD) and split from the composite 2BEG PDB file. Accordingly, 1Z0Q is a homo 30-mer composed of 30- monomers (A), each comprising an individual coordinate model. The A1 monomer was split from the composite 1Z0Q PDB file. At least two amyloid molecules are required to achieve the repeating structure of a protofilament fibril. Therefore, the dimer (Figure 2) or the trimer (Figure 6) models were generated by RossettaDock  , as described previously  . Subsequently, the generated dimers (A1B1) or trimers (A1B1C1) were utilized for ThT and DHA docking (both flexible and rigid docking).
2.5. Computational Analysis of Binding Sites in Aβ1-42 Dimers or Trimers
The binding sites or pockets of the dimer (A1B1) and trimer (A1B1C1) were determined by GHECOM, which detects grid-based pockets/binding sites on the surface/interior of the protein  . The program produces a graph of residue-based “pocketness,” and the presence of binding sites was visualized using the 3D molecular viewer Jmol  .
2.6. Inter-Surface Interaction Site Analysis of the Aβ1-42 Dimer (A1B1)
Protein-protein (the monomer-monomer A1B1 dimer) inter-surface interaction sites were analyzed by feeding the A1B1 dimer to the cons-PPISP server  , which predicts the residues that are likely to form the binding site for neighboring proteins. The energy distribution along the protein-protein interface is not homogenous;
Figure 2. The monomeric structure files for Aβ1-42 were prepared from their parent composite PDB files (2BEG and 1Z0Q). Subsequently, dimers of Aβ1-42 were formed in order to dock ligands (ThT and DHA). (a), (b) Monomer; (c), (d) Dimer.
certain residues contribute more to the binding free energy, called “hot spots.” Hot-point servers predict hot spots in protein interfaces using an empirical model that incorporates a few simple rules consisting of occlusion from solvent and total knowledge-based pair potentials of residues  . To further validate the presence of inter-surface hotpoints on the dimer interface binding contacts, the dimer was fed to hot-point prediction servers, including Knowledge-based FADE and Contacts (KFC2)  .
2.7. Molecular Docking Simulations of ThT and DHA onto the Aβ1-42 Dimer
The molecular docking simulations were performed using MVD  and PatchDock  .
2.7.1. Molegro Virtual Docker (MVD)
MVD version 4.3.0 (free academic version), along with its graphical user interface (MVD tools), was utilized to generate a grid, calculate the dock score, and evaluate conformers. The nonpolar hydrogen atoms were removed from the receptor file and their partial charges were added to the corresponding carbon atoms. Docking was performed by following the steps in the MVD user manual  . MolDock score of ligands were calculated during docking. MVD performs flexible ligand docking with optimization of the ligand geometry during docking. Therefore, bond angles, bond lengths, and torsional angles of the ligand are modified during the stages of receptor-ligand complex generation.
PatchDock is a shape complementarity/geometry-based molecular docking algorithm. It is aimed at finding docking transformations that yield good molecular shape complementarity. Such transformations induce both wide interface areas and small amounts of steric clashes. A wide interface is ensured to include several matched local features of the docked molecules that have complementary characteristics. The output of PatchDock is a list of candidate complexes between the receptor and ligand molecules. The list is sorted according to geometric shape complementarity score, approximate interface area of the receptor-ligand complex, atomic contact energy (ACE) between ligand and receptor, and 3D transformation. Finally, the server provides an option to download the ligand-receptor complexes in the PDB format. The PatchDocked ligand-receptor complex (of Rank 1) was visualized using PyMOL  .
The predicted physical (ADME) properties of DHA and ThT, as evaluated by the FAF-Drug 3, are shown in Table 1. The major difference between the properties of DHA and ThT was in the LogP (octanol-water partition coefficient) value, which was >5 for DHA, but <5 for ThT. Thus, DHA violated Lipinski’s rule of five and this may relate to the more hydrophobic acyl chains in DHA compared with those in ThT. Partition coefficients indicate the distribution of drugs within the body. Hydrophobic drugs with high octanol/water partition coefficients are mainly distributed to hydrophobic areas such as lipid bilayers of cells. The number of rotatable flexible bonds helps a molecule to adopt 3D structure with a greater flexibility, while that of the rigid bonds performs the opposite trend. Molecular polar surface area (PSA) is a very useful parameter for prediction of drug transport properties. Polar surface area is defined as a sum of surfaces of polar atoms (usually oxygens, nitrogens and attached hydrogens) in a molecule. This parameter is used to correlate very well with the human intestinal absorption, monolayers permeability, and blood-brain barrier penetration. Usually, the molecules with a polar surface area of greater than 120 - 140 Å2 tend to be poor in permeating cell membranes. The oral bioavailability scores based on solubility in water (Sw) were 2.29 × 10−5 and 2.14 × 10−5 for DHA and ThT, respectively. The log of these two values, i.e. Sw, was −4.64 and −4.67, respectively. This indicates that the water solubility of these two compounds were comparable. Because of the Lipinski’s rules (LR) are based on physicochemical criteria for physiological effect of drugs (ADME mainly) but does not rule out a possible direct interaction between a molecule and a protein outside the cellular context. Moreover, LRs are violated by numerous drugs or physiologically relevant molecules. Notably, the LR is a group of physicochemical properties used to evaluate the probability of a substance to become an effective drug. However, the conclusions that can be drawn from a docking exercise are limited. If a docking is good and a ligand gets a good score, this just suggests that the compound might have activity against a given target. Docking says nothing about physicochemical properties.
Aggregation-propensity analyses of the Aβ1-42 peptides are shown in Table 2. The hot-spot analyses by FoldAmyloid  , AGGRESCAN  , and ProA  revealed that the amino acid residues from Val24 to Lys28/ Gly29 and those from Val39 to Ala42 were the most susceptible to amyloidogenesis.
Beta-aggregation propensity analysis also revealed that the amino acid residues from Lys16 to Asp23, and those from Lys28 to Ala42, displayed the highest beta-aggregation tendency (Figure 3). These two regions (Lys16-Asp23 and Lys28-Ala42) may adopt β-sheet secondary structures, and hence, have a higher beta-aggre- gation propensity during fibrillation.
When a protein molecule binds to another biological polymer (protein) to form a complex, the subset of residues in the interface that account for most of the protein binding free energy are called binding hot spots. Therefore,
Table 1. In silico-predicted ADME properties of DHA and ThT using FAF-Drugs 3 properties.
FAF-Drugs 3 evaluated drug-like properties of ThT and DHA. FAF-Drugs 3 is a program for filtering large compound libraries prior to in silico modeling studies. The tool can perform computational prediction of some ADME-tox properties (Adsorption, Distribution, Metabolism, Excretion, and Toxicity) and offers free online services for calculation of important molecular properties (LogP, polar surface area, number of hydrogen bond donors, and number of hydrogen bond acceptors), as well as prediction of the bioactivity score for the most important drug targets.
Table 2. Analyses of aggregation-prone amino acid residues in the Aβ1-42 peptide sequence.
1FoldAmyloid  , 2AGGRESCAN  , and 3ProA  servers employ algorithms to predict amyloidogenic regions from protein amino acid sequence. The analyses are based on the amino acid physicochemical properties critical to protein aggregation. Color blocks indicate “common” amino acid residues those are prone to aggregation, as analyzed by 1,2,3Different servers.
Figure 3. Beta-aggregation propensity analysis of Aβ1-42. The red line indicates the region of the amino acid sequence (which is shown as single-alphabet strip) that induces beta aggregation.
we analyzed the protein-protein interface interaction sites for the A1B1 dimer. Because amyloid fibrillation starts with the interactions between at least two β strands, identifying these interaction sites is essential for understanding the nature of these sites, and whether DHA affects them, and the subsequent fibrillation. The Phe19 and Phe20 of the A1 monomer interacted with the Phe19 and Phe20 of the B1 monomer. The “Hot-point” server principally considers the solvent accessibility and total contact potential of the interface residues. The output file tabulates the interface residues with the highlighted hot spots and their features. The “hot-spot” with the highest potential was located in Leu34 on chain B1 of 2BEG. The KFC2 showed that the interaction sites were at Phe19 and Phe19 of A1 and B1. Other amino acids from either monomeric chain involved in the inter-surface interactions sites are highlighted in Table 3.
Pocketness is an important parameter that helps to determine the binding sites for ligands to their receptors. The 2BEG dimer or trimer showed ~5 or 6 pockets, whereas the 1Z0Q dimer displayed 4 pockets (Figure 4).
Docking performed by MVD provided the five best poses, and we obtained the corresponding MolDock score values and other thermodynamically calculated values. The 3D structures of the best scoring dockings are shown in Figures 5(a)-(c). Steric interactions of ThT and DHA with Aβ1-42 (the receptor) were also evaluated with MVD. Maps of the interactions of ThT with amino acid residues of the 2BEG dimer (A1B1) are shown in Figure 5(a1). The docking of DHA to dimeric Aβ1-42 (A1B1 of 2BEG and A1A2 of 1Z0Q) and their interaction maps are shown in Figure 5(b), Figure 5(b1), and Figure 5(c), Figure 5(c1), respectively. Using MVD, DHA had a MolDock score of −122.40 Kcal/mol and −140.40 Kcal/mol when docked onto the 2BEG and 1Z0Q dimers, respectively. The ThT reference molecule had a MolDock score of −130.37 Kcal/mol when docked onto the 2BEG dimer. Therefore, the binding of DHA to the 2BEG dimer was comparable with that of ThT, as demonstrated by similar binding energies (Figure 5(b)).
Amino acids involved in the docking of DHA to the Aβ1-42 dimer (2BEG, A1B1) and Aβ1-42 trimer (1Z0Q, A1A2) were visualized by contact maps of the interactions between atoms of DHA and the amyloid amino acid
Table 3. Analyses of momomer (A1)-monomer (B1) inter-surface interaction sites of Aβ1-42 dimer (A1-B1).
1Cons-PPISP is a consensus neural network method for predicting protein-protein interaction sites. Given the structure of a protein, cons-PPISP will predict the residues that will likely form the binding site for another protein. The inputs to the neural network include position-specific sequence profiles and solvent accessibilities of each residue and its spatial neighbors. The neural network is trained on known structures of protein- protein complexes  . 2Hotpoint: Hotpoint, which predicts hot spots in protein interfaces using an empirical model. The empirical model incorporates a few rules consisting of occlusion from solvent and total knowledge-based pair potentials of residues  . 3KFC2 (Knowledge-based FADE and Contacts) server predicts binding “hot spots” within protein-protein interfaces by recognizing structural features indicative of important binding contacts. The server analyzes several chemical and physical features surrounding an interface residue and predicts the classification of the residue using a model trained on prior experimental data. Color blocks indicate “common” amino acid residues those acting at the inter- surface interaction sites, as the Aβ1-42 was analyzed by different servers.
Figure 4. Binding sites (pockets) analysis of Aβ1-42. One of the simplest ways to predict ligand binding sites is to identify pocket-shaped regions on the protein surface. A binding sites or a pocket is defined as a space into which a small probe can enter, but a large probe cannot. A pocket intrinsically has two arbitrary properties, size and depth. The GHECOM analyzes these pockets using probe spheres. The radii of the probe spheres are assumed to correspond to the size and depth of the pockets. These values can be adjusted to individual putative ligand molecule. The binding sites or pockets of the dimer (A1B1) and trimer (A1B1C1) were determined by GHECOM, which detects grid-based pockets/binding sites on the surface/interior of the protein  . The presence of binding sites was visualized using the 3D molecular viewer Jmol  . The program also produces a graph of residue-based “pocketness” (right panel). (a1), (b1), (c1) number and position of amino acids (represented as lines on the primary sequences) involved in each of the pocket of (a), (b) and (c) of the left panel. For details please see Reference  . The individual color indicates the presence of distinct pockets.
Figure 5. The best-scored docking poses of ThT and DHA with the receptor (Aβ1-42) by the MVD. (a) ThT docked onto the 2BEG A1B1 dimer. (b) DHA docked onto the 2BEG A1B1dimer. (c) DHA docked on the 1Z0Q A1A2 dimer. Interaction (binding) energies (∆G) are also shown. Contact (binding) maps of the interactions between atoms of ThT/DHA and amino acid residues of the 2BEG/1Z0Q dimer are shown in the right panel (a1) (b1) and (c1). Contact maps were visualized by the “Ligand map” module of the MVD.
residues using MVD (Figure 5(b1) & Figure 5(c1), and Table 4). The amino acid residues of the A1B1 dimer found to sterically interact with ThT were Ile32, Gly33, and Leu34 of A1 as well as Ile31, Ile32, Gly33, and Leu34 of B1.
We also performed molecular docking of ThT and DHA using PatchDock to confirm the ability of DHA to bind the 2BEG A1B1 dimer and the 2BEG A1B1C1 trimer (Figure 6). The algorithm used in PatchDock performs rigid docking, with surface variability/flexibility implicitly addressed through liberal intermolecular penetration. Geometric shape complementarity scores were similar for DHA and ThT when docked with the dimer (Score: ThT = 3742; DHA = 3732). However, there was higher shape complementarity when DHA was docked onto the trimer (Score: ThT = 3632; DHA = 4326). Approximate complex interface areas (receptor-ligand) were also higher for DHA. ACEs for ThT and DHA were almost similar during docking onto the dimer; however, it was slightly higher for ThT when docked onto the trimer (Table 5). These data suggest that DHA and ThT have similar binding affinities.
Table 4. Patch docking of ThT and DHA with the 2BEG amyloid dimer (A1B1) and trimer (A1B1C1).
Abbreviations: Score, geometric shape complementarity score; Area, approximate interface area of the complex (receptor-ligand); ACE, atomic contact energy; A1B1 dimer, generated by feeding the A1 and B1 model 1 of 2BEG to the Rosetta Server; A1B1C1 trimer, generated by feeding the A1, B1, and C1 coordinate model 1 of 2BEG to the Rosetta Server; ThT, thioflavin T; DHA, docosahexaenoic acid.
Table 5. Summary of dockink of DHA and ThT on to NMR structures of Aβ1-42 dimer structures (A1B1 of 2BEG and A1A2 of 1Z0Q).
When DHA and ThT were docked onto 2BEG dimer, they shared a common binding site at Ile31of A1 peptide. Both the 2BEG and 1Z0Q pdb had common binding sites at Ala21 of A1 and Glu21 of B1 peptide. Color blocks at a given position (of the sequence) indicate common binding sites for the ligand.
Figure 6. The best scoring results of the Patch Dock for ThT (a)-(c) and DHA (b)-(d) with the A1B1 dimer (a) (b) and the A1B1C1 trimer (c) (d). Ligands (ThT and DHA) are represented by green sticks, and the receptors (dimer and trimer) are shown as β-sheet ribbons with the amino acid side chains presented as sticks.
The ligand-docked receptor obtained from PatchDock was visualized using PyMOL as shown in Figure 6. ThT bound with the interior regions of both the dimer and trimer (Figure 6(a) & Figure 6(c)), whereas DHA was found to wrap around the β-sheet of the dimer (Figure 6(b) & Figure 6(d)).
Alzheimer’s disease (AD) is one of the most common chronic dementia-related diseases worldwide, and the number of AD patients is increasing every year. Currently, only two medications (cholinesterase inhibitors and an NMDA receptor antagonist) are approved by the United States Food and Drug Administration  to treat the cognitive symptoms (memory loss and confusion) of AD. However, these medications cannot cure AD or arrest its progression. Vitamins or other supplements, such as omega-3 fatty acids  , and various herbal mixtures are also widely promoted as preparations that may support cognitive health or prevent or delay AD. The NIH panel concluded that there is somewhat stronger data―but not definitive evidence―that omega-3 fatty acids in fish oil may help prevent cognitive decline. The basis of such conclusion might relate to the fact that the concentration of omega-3 polyunsaturated docosahexaenoic acid (DHA, C22:6, n-3) is significantly reduced in the brains of AD patients   . We also reported on the beneficial effects of dietary omega-3 polyunsaturated fatty acid supplementation on age-related cognitive decline in elderly Japanese patients with very mild dementia  . Therefore, DHA supplementation may decelerate the cognitive symptoms of AD.
We previously showed that chronic administration of DHA significantly decreased the amyloid burden in the brains of AD model rats  . DHA also inhibited the in vitro amyloid fibrillations of Aβ1-40  , Aβ1-42  , and Aβ25-35  , indicating that neurotoxicity convened by these amyloid species would be inhibited by the DHA. Moreover, we found that Aβ1-42-induced toxicity of SH-S5Y5 cells was reduced by DHA  , andour previous transmission electron microscopy, laser-scanning fluorescence microscopy, and thioflavin T fluorospectroscopy studies clearly indicated that DHA resists fibrillogenesis by some unknown mechanism    . Despite strong evidences that soluble forms of Aβ are toxic at synapses, the exact forms of the toxic species remain to be determined. There is evidence for toxicity of both low molecular weight dimers and trimers of Aβ  -  . Using Western blot analysis and/or polyacrylamide gel electrophoresis, we have also shown that DHA inhibits the fibrillation of dimeric/trimeric and oligomeric Aβs  ; thus, it ultimately decreases the elongation process of amyloid fibrillation leading to matured fibers. Finally, the anti-Aβ fibrillation effect of DHA was further supported by the reductions of Aβ-induced apoptosis and neurodegeneration  in DHA-treated cell culture. These in vivo/in vitro experimental data strongly suggest that DHA prevents amyloid peptides from forming toxic species that ultimately deteriorate memory and other behavioral aspects in AD patients.
These findings also demonstrate that Aβ peptides are important drug targets, and designing drugs that disrupt the formation of Aβ-Aβ inter-surface interaction sites is particularly important for preventing Aβ fiber formation. Therefore, we performed the present in silico study to support the experimentally proven ameliorative effects of DHA on the Aβ-fibrillation. Accordingly, we modeled the 3D structure of DHA and used docking simulation to determine whether DHA binds to Aβ peptides and to predict the affinity between the ligand (DHA) and the receptor (Aβ). During docking simulations with both MVD and PatchDock, we used ThT as a positive control to verify whether docking would yield ligand geometry similar to that of DHA when binding to Aβ peptides.
Using MVD, DHA had a MolDock score of −122.4 Kcal/mol with 2BEG, and −140.40 Kcal/mol with 1Z0Q, whereas ThT had a MolDock score of −130.37 Kcal/mol. These results suggest considerable binding of DHA with the dimeric Aβ1-42 molecules. This might explain the reduced amyloid burden observed in the brains of DHA-administered AD model rats  and concurrent inhibition of in vitro fibril formation  . Monomer- monomer inter-surface interaction site analyses revealed that Phe19 (A1) and Phe19 (B1) were the common amino acid residues, which acted as monomer-monomer inter-surface interaction sites (Table 3). Other amino acid residues, including Phe20 (A1), Phe20 (B1), Ala21 (A1), Asp23 (A1), Gly25 (A1), Lys28 (B1), Ile32 (B1), Gly33 (A1), Leu34 (B1), Met (A1), Val36 (A1), Val36 (B1), Gly37 (A1), Gly37 (B1), Gly38 (A1), Val39 (B1), Val 40 (A1), and Ile41 (B1) also acted as inter-surface interaction sites (Table 3). These inter-surface interaction sites might have been targeted by the DHA, as indicated by the interactions of DHA with these amino acid regions, particularly, those of Phe19, Ala21, Asp23, and Leu34. Notably, DHA also interacted with Ala21, Glu22, Asp23, and Ser26-Ile31 of 1Z0Q, which is in an a-helix random coil conformation rather than the b-sheet conformation that is the prerequisite for the propagation of fibrillation. As shown in Table 2, the amino acid residues 24 - 29 and 39 - 42 are the most aggregation-sensitive regions of Aβ1-42. All these results suggest that DHA interacts with these aggregation-prone monomer-monomer inter-surface interaction sites. Moreover, ThT docked with 2BEG in a similar manner to DHA and bound to Ile31, Ile32, Gly33, and Leu34.
In the PatchDock system, the docking of DHA and ThT were similar, as indicated by the almost similar docking score. However, when they were docked by PatchDock onto the Aβ1-42 trimer, DHA displayed a greater binding affinity (score) than ThT. When the PatchDock receptor-ligand was visualized with PyMOL (Figure 6), an interesting feature was obtained: DHA spirally wrapped the dimer, whereas ThT bound to the interior of the dimer. Notably, DHA and ThT both bound to the predicted binding pockets (Figure 4). The binding interaction between the ligand molecule (DHA) and receptor (Aβ1-42) has a significant role in determining the anti-amyloid activity. Therefore, the binding potential of DHA in this computational study supports the results of our previous in vivo and in vitro biological studies showing that DHA decreases in vitro amyloid fibrillation and amyloid brain burden in a rat AD model, with concurrent amelioration of Aβ-induced memory loss     . AD drug discovery has been a very lengthy and costly process, and the prevention of AD remains a burning issue worldwide. Polyunsaturated fatty acids, including DHA, have been continually investigated over the last several years as potential drug candidates with minimal side effects in experimental animal models of AD and/or AD human patients. Using a targeted and more specific approach for future high-throughput screens, such as those based on the results of similar in silico experiments performed here, we may be able to more rapidly identify novel AD drugs.
Computational drug design is a promising method for discovering many more drugs in the coming years. Here we demonstrated that molecular docking of DHA is consistent with its anti-amyloid effects, in particular, on in vitro fibrillation and in vivo Aβ-induced memory impairments of AD model rats. Our study further reaffirms that DHA is important for the structure and biological function of Aβs. We also showed that scoring the docking values gives the best prediction for ligand interactions with proteins. Finally, our in silico results provide compelling evidence for the utility of DHA as a preventive medication for neurodegenerative diseases such as Aβ1-42-induced AD. Further studies on DHA-like molecules hold promise for the development of anti-amyloid drugs.
This work was supported in part by a Grant-in-Aid for Scientific Research (C) from the Ministry of Education Culture, Sports, Science and Technology, Japan (26500008 to M.H.). We also gratefully acknowledge the MVD authority for providing us with the free academic version.
Conflict of Interest
The authors declare no-conflict of interest.
AD: Alzheimer’s Disease;
Ab: Amyloid Beta Peptide;
RCSB: Research Collaboratory for Structural Bioinformatics;
PBD: Protein Data Bank;
MVD: Molegro Virtual Docker;
ThT: Thioflavin T;
DHA: Docosahexaenoic Acid
Submit your manuscript at: http://papersubmission.scirp.org/