EVALUATION OF NOVEL LIGAND FOR THE MAINTENANCE OF NATIVE Aβ(1-42) PEPTIDE CONFORMATION: RELEVANCE TO ALZHEIMER’S DISEASE
HTML Full TextEVALUATION OF NOVEL LIGAND FOR THE MAINTENANCE OF NATIVE Aβ(1-42) PEPTIDE CONFORMATION: RELEVANCE TO ALZHEIMER’S DISEASE
Dilawar Ahmad Mir and R. Boopathy Rathanam *
Department of Biotechnology, Bharathiar University, Coimbatore - 641046, Tamil Nadu, India.
ABSTRACT: Alzheimer's disease (AD) is a disorder of the central nervous system with progressive neurodegeneration, cognition, and memory loss. A major molecular hallmark of the disorder includes extracellular deposition of the Amyloid beta-peptide (Aβ) in senile plaques, the appearance of intracellular neurofibrillary tangles. The Aβ peptide is produced by sequential cleavage of the Amyloid precursor protein (APP) by α, β, and γ-secretase. The secretase β and γ generate a number of isoforms of peptides containing 36-43 amino acid residues in length. The mainly universal isoforms are Aβ(1-40) and Aβ(1-42). Aβ(1-42) C-terminal domain (29-42 amino acid residues) adopts β conformation, and the N-terminal domain (10-24 amino acid residues) facilitates a dynamic equilibrium between α- helix and β-strand. Thus Aβ polymerization requires unfolding of the native α-helical structure of Aβ. Stabilization of the Aβ central α- helix is an efficient step to foil the Aβ polymerization. Here we report test compounds which to bind and stabilize the 11-30 amino acid regions of Aβ(1-42) in α-helical conformation. On docking, the ligands with Aβ peptide followed by Molecular dynamics simulation for 20ns favored the identification of test compound NIAID. That postulated to bind and stabilize the Aβ central α-helix. So stabilization of Aβ secondary structure with a ligand that maintains the α-helical conformation may provide clues in developing drugs to control the Aβ deposition observed in AD patients.
Keywords: |
Alzheimer’s disease, NIAID, Molecular Docking and Molecular dynamics simulation, Aβ(1-42) α-helix stabilization
INTRODUCTION: AD is caused by fibril and plaque formation of the Amyloid beta-peptide (Aβ) in the central nervous system 1, 2. Aβ assemblies include intermediate to monomeric and fibrillary forms which are considered as the source of cytotoxicity 3. Such Aβ assemblies include low-number oligomers and larger assemblies known as protofibrils, globulomers 4, 5.
The some of the Aβ peptides are an integral membrane protein, called the amyloid precursor protein (AβPP), cleaved predominantly into elongated 40-residue peptide (Aβ 1-40).
Also, a C-terminally elongated 42-residue version which also can be excised (Aβ1-42) and this longer variant is the main constituent of parenchymal amyloid deposits 6, 7. The link between Aβ aggregation and AD implies that any inhibitors of their aggregation should be able to slow down disease progression. Several low molecular mass Aβ aggregation inhibitors have been identified by screening library of the compound as well as rational design strategies. They include chemically diverse compounds such as Rifampicin, Curcumine, Inositol, Melatonin, D737, and Dobutamine 8, 9, 10. These inhibitors are predicted to bind to Aβ in an elongated β-strand-like conformation and prevent its polymerization. A potential problem with this strategy is that blocking the later stages of fibril formation will not favor the formation of prefibrillar oligomeric forms that are cytotoxic 11, 12, 13. No Aβ oligomers have been structurally determined, and thus these oligomers cannot yet be targeted by rational design. A further problem with some of the Aβ polymerization inhibitors is that they can also act as aggregators 14. Chemical aggregators can physically sequester proteins in a promiscuous and nonspecific manner which is typically made up of conjugated aromatics, hydrophobic and dye-like features 14, 15 and also a library of peptide inhibitor candidates have formed which systematically are mutating the RGTFEGKF amino acid residues 16. In light of these circumstances, an alternative strategy to reduce Aβ aggregation and toxicity is warranted. One possible approach, explored herein, is to trap Aβ in a state similar to its perceived structure in membrane-embedded AβPP, by targeting the discordant α-helix region.
Nuclear magnetic resonance (NMR) data showed that Aβ (1-40) adopts a folded structure including two α-helical regions (residues 15-24 and 29-35) in water or sodium dodecyl sulphate (SDS) micelles which provide a water membrane interface mimicking environment and that Aβ(1-42) adopts an unfolded structure including two β-strands (residues 17-21 and 31-36) in aqueous solution 17, 18. Using NMR, it has also been shown that an Aβ(1-42) fibril is a β-sheet composed of two β-strands (18-26 and 31-42) 19. This structure of Aβ is the unfolded elongated β-strand that departs from the membrane structure, (α-helical) region like forms and that the β-strand of Aβ enable the formation of β-sheets of fibrils and prefibrillar aggregates 20, 21. Recent experimental studies demonstrated that trapping Aβ in a state similar to its native structure by stabilizing the Aβ central helix (residues 15-24) is an effective strategy to reduce Aβ polymerization and Aβ toxicity 22, 23. The 3D structure of a disease-relevant Aβ (1-42) fibril polymorph has also been deciphered 24. The 3D structure is composed of two molecules per fibril layer, forming a double-horseshoe–like cross–a β-sheet entity with maximally buried hydrophobic side chains.
Residues 1–14 are partially ordered and in a β-strand conformation 24. The actual secondary structure of Aβ(1-42) have been determined it is about, 45% in two α-helical forms (α1: residues from 11–24, α2: residues from 27–33) and are joined by a flexible hinge.
In the present study, the effect of the test compound NIAID screened from PubChem library on the unfolding process of the Aβ central helix (11-24) was investigated by Docking and Molecular Dynamic (MD) simulations. The difference between native-type Aβ with and without ligand has been studied during MD simulation, and here we demonstrated that the given ligand is effective in stabilizing the Aβ central helix. The NIAID compound was finally identified as potential lead molecule based on stabilizing the Aβ central helix.
This compounds can be easily synthesized and have structural novelty, which allows for further examination of their abilities to inhibit Aβ confirmation through in-vitro and in-vivo biological tests and will act as potential agents to treat Alzheimer’s disease.
MATERIALS AND METHODS: The α-helical form of Aβ peptide (14–23 residues) provides a suitable target for stabilizing ligands. The Peptide Aβ(1-42) (1ZOQ) has retrieved from PDB, which has Resolution of 2.37 Å, R-Value of 0.214 (obs.), R-Free energy of 0.226. For possible interaction of human Aβ, we concentrated on 2 partial surfaces of this Aβ helix. One surface of the helix is largely hydrophobic and adjacent to Glu-22 and Asp-23. A compound conjugated with hydrocarbon chain could potentially interact with this surface. Also in Aβ helix adjacent to Glu-22 and Asp-23 is a surface containing Phe-20, which connects these residues with Lys-16 and His-14. To investigate the concept, we have evaluated possible chemical interactions of Aβ peptide by 350 test compounds that could be initial potential ligands for these surfaces, all compounds were retrieved from PubChem. The potential test compounds and the Aβ peptide were submitted to AutoDocking 25.
Finally, docked complex of each test compounds and Aβ peptide was used to do further analysis and molecular dynamics simulation. In AutoDock analysis, we have selected 30 test compounds, on the bases of their lowest binding energy score and interaction with Aβ middle amino acids residues (E11VHHQKLVFFAEDVGS26). The selected test compounds are given in Table 1.
TABLE 1: INTERACTION ANALYSIS OF TEST COMPOUNDS WITH Aβ1-42. THE COMPOUNDS WERE INDEPENDENTLY DOCKED ONTO Aβ1-42 USING AUTODOCK. THE 2D STRUCTURE OF THE TEST COMPOUNDS ALONG WITH THEIR BINDING ENERGY AND RESIDUES OF Aβ1-42 PARTICIPATING IN THE H-BONDING AND HYDROPHOBIC INTERACTION ARE SHOWN IN TABLE. THE COMPOUNDS WERE RETRIEVED FROM PUBCHEM. RESIDUES OF Aβ1-42 INVOLVED IN THE BINDING OF THESE TEST COMPOUNDS ARE GIVEN ACCORDING TO THE ASCENDING ORDER OF BINDING ENERGY
Name of Compounds | Residues of Aβ(1-42) involved in the interaction with test compounds | Binding
Energy (Kcal/mol) |
|
Hydrogen Bonds | Hydrophobic interactions | ||
24721112 | K16, F19, D23 | V12, Q15, F20, V24, N27, A30, I32, G33, V36 | -14.4 |
11313 | H14, F20 | V12, H13, Q15, K16, F19, D23, V24, N27, A30, I32 | -10 |
688035 | K16, G33 | Q15, F19, F20, V24 | -8.7 |
3708931 | D23, I32 | Q15, F19, E22, G29, A30, G33, L34 | -8.1 |
8786203 | H14, F19 | D23, N27, G29, A30, I31, G33, | -7.5 |
552906 | K16, F20 | Q15, F19, D23, V24, N27, G29, A30 | -7.2 |
62770 | F19 | G15, K16, D23, N27, G29, A30I32, G33, L34 | -7.0 |
4534 | V12, D23 | H13, Q15, K16, F19, F20 | -7.1 |
91895552 | F19, D23 | V12, Q15, K16, V24, N27, A30, I32, G33 | -7.0 |
2585 | G15 | K16, F19, F20, G22, D23, A30, I32, G33 | -6.9 |
5757 | D23 | L17, F19, V24, A30, G33 | -6.7 |
4534 | V12, H13, D23 | G15, K16, F19, F20 | -6.1 |
9990 | V12, K16 | Q15, L17, F19, F20, D23, A30 | -6.0 |
15818601 | D23, G29 | Q15, K16, F19, F20 | -6.7 |
896 | Q15, K16, F19, F20, D23 | H14, L17, A30, I32, G33 | -6.1 |
181139 | H14, D23 | G29, A30, I32, G33 | -5.8 |
128419 | K16, F20, D23, I32 | F19, D23, V24, A30 | -5.5 |
3018172 | K16, F19 | H13, G15, A30 | -5.3 |
8786203 | H13, H14, L17 | Q15, F19, F20, V24, G29, A30, G33 | -5.2 |
3853025 | H13, H14, D23, I32 | Y10, K16 | -5.0 |
681 | V12, H13, L16, D23 | Q15, L17, F19, F20, A30, I32, G38 | -5.0 |
23043846 | H14, F19 | Q15, V24, A30, G33 | -5.0 |
5381226 | H13, D23 | H14, L17, A30, I32, G33 | -5.0 |
638015 | G15, D23 | V12, K16, F20, L17, F19, V24 | -5.0 |
5281672 | Q15, K16, F20, D23 | F19, F20, V24, N27, A30, I32 | -5.0 |
55290627 | Q15, D23, I32 | K16, F19, F20 | -5.0 |
54267694 | V12, D23 | F19, F20 | -5.0 |
46856272 | L16, D23 | F19, V24, A30 | -5.0 |
18347908 | V12, K16 | V24, DN27, A30 | -4.9 |
Molecular Dynamic Simulation: The entire MD simulation was performed using the Macro model package from the Schrodinger Program. The force-field parameters for the ligand and protein were picked from the AMBER force field. For Aβ(1-42) stochastic dynamic method is used. The SHAKE algorithm was applied to fix all covalent bonds containing a hydrogen atom allowing a 1.5fs time step to be used in the integration of Newton’s equation. To observe the structural changes of Aβ quantitatively after MD simulation, the root means– square deviation (RMSD) and deviation of Torsion angle were calculated. The RMSD is calculated for backbone heavy atoms against the time coordinates. Therefore, the simulation for each system was performed at 300-330K to accelerate the dynamics of Aβ. RMSD analyses were carried out by using the trajectories obtained by MD simulations at 320K for native Aβ(1-42) and simulations analyses for Aβ (1-42) with the test compounds was carried out for the trajectories obtained up to 20ns at 320K. The data of every 2ns of the trajectories after the heating time of the MD simulations were used for the revelation of the structural changes of Aβ with and without test compounds during simulations was carried out by using the visual molecular dynamics (VMD) Software. The test compound NIAID (PubChem SID = 24721112) was the only ligand which retains the native α-helical conformation of Aβ(1-42) during MD simulation upto 20ns at 320K while as other test compounds have failed to retain the Aβ in alpha helical conformation. To discriminate the type or pattern in the Aβ structure, the number of α-helical backbone hydrogen bonds (αHB) in the middle region (11-24) was calculated, using the criterion acceptor hydrogen ≤ 2.4Aº to define the existence of hydrogen.
Result: Targeting the Aβ central α- helix (15-24 residues) is an effective strategy to reduce Aβ polymerization and toxicity. In AutoDock analysis test compound with lowest binding energy and interaction with Aβ middle amino acids residues (E11VHHQKLVFFAEDVGS26) was criteria to predict the favored compounds for dynamic molecular studies. The Docking simulation interactions of Aβ(1-42) with test compounds are given in Table 1.
The docking studies of various test compounds to the fragments of Aβ reveal that the stretch of residues from K16, to F20 broadly interact with the compounds and interestingly they are the key coordinating residues that are proposed to be responsible for conversion of Aβ alpha-helix into beta-sheet events. Among the entire test compounds, NIAID exhibits a favored binding to Aβ(1-42) segment and have shown the highest docking-score ≥ -14 Kcal/mol. The docking interaction analysis showed NIAID is interacting with the middle amino acid residues of Aβ(1-42) (K16L17V18F19 F20D23). The average numbers of hydrophobic interaction with Aβ might have increased the binding affinity of NIAID test compound with Aβ interfaces.
FIG. 1: STRUCTURAL DISTORTIONS IN THE Aβ(1-42) PROTEIN AT DIFFERENT TIME INTERVALS DURING THE MD SIMULATION: A) NATIVE AΒ(1-42), B) Aβ(1-42) WITH NIAID. FIGURE A SHOWS THE UNFOLDING OF HELICAL Aβ WHEREAS SUCH EVENT IS ABSENT IN FIGURE B
The further signify the docking interaction of diverse test compounds; these could be exploited as competitive inhibitors to facilitate the prevention of AD events. The AutoDock built complexes between Aβ (1-42), and the test compounds were submitted to MD simulations. The MD simulation of Aβ (1-42) with and without test compounds up to 20ns was analyzed to understand the structural distortion of Aβ (1-42) during the simulation.
From a thousand trajectories we have mentioned here five trajectories of different nano-seconds (0, 5, 10, 15, and 20). The MD simulations of Aβ(1-42) with and without test compound were performed to analyze the structural alteration of native Aβ(1-42) are shown in Fig. 1A. The structure of these trajectories showed the native Aβ central helix unfolded during the simulation at different nano-seconds.
Simulations structure at different periods showed that the central α-helix of Aβ(1-42) is completely unfolded without the ligand up to 20ns at 300K Fig. 1A. The amino acid K28G29A30 and E22D23V24G25S26 in the middle region of Aβ(1-42) are relatively disordered and converted into β-sheet at 5ns and 10ns respectively. MD simulation of Aβ(1-42) peptide with NIAID compounds was performed up to 20ns at 320K showed Aβ peptide retains its α-helix form The NIAID was only tested compound which retained the Aβ native confirmation Fig. 1B. During MD simulations, the secondary structures adopted or retained by Aβ with and without ligand are shown in Fig. 2.
FIG. 2: SECONDARY STRUCTURES ADOPTED BY Aβ(1-42) PROTEIN INITIAL (A) AND FINAL (B-C) DURING A 20ns MD SIMULATION WITH AND WITHOUT THE LIGANDS A: NATIVE Aβ(1-42) at 0ns, B: NATIVE Aβ(1-42) at 20ns IN ABSENCE OF LIGAND AND C: Aβ(1-42) at 20ns IN PRESENCE OF NIAID
It is proportionately illustrated that the stability of Aβ central α-helix is retained by NIAID. Starting with Aβ(1-42) in α-helical conformation, Aβ alone starts to unfold and lose its helical structure during the simulation, were in the presence of NIAID the α-helical conformation of Aβ is retained or refolded after partial unfold at the end of the simulations. This is postulated to roll and stabilize the Aβ central α- helix. MD simulation result analysis favored the identification of NIAID compounds as the best inhibitors for the conversion of α-helix of the Amyloid peptide into β-sheet conformation. The NIAID 2D chemical structure is shown in Fig. 3.
FIG. 3: THE 2D STRUCTURE OF THE TEST COMPOUND NIAID
To examine the structural change of Aβ quantitatively at different nano-seconds during simulation, the RMSD was calculated. The RMSD was calculated for heavy backbone atoms against the initial energy-minimized coordinates calculated for all atoms along the MD simulation time. During MD simulation of native Aβ (without ligand) we have analyzed that the central helix of Aβ is completely unfolded up to 20ns at 300K. The Mean RMSD value of native Aβ at 300K increased from 2ns to 10 ns and showed there is large instability in the structure of native Aβ α-helix after 5ns and changed completely into β-sheet up to 10ns. The Aβ without ligand showed central α-helix becomes fully extended and more flexible at the end of simulations. RMSD of native Aβ is shown in Fig. 4A. It has been analyzed that the average backbone RMSD of the middle region and the average number of alpha-helix hydrogen bonds (αHBs) showed fluctuation in the backbone. In the presence of NIAID backbone, RMSD of the Aβ up to 20ns was relatively stable during simulations, RMSD was small in every trajectory, and Aβ α-helix was relatively stable during simulation Fig. 4B.
FIG. 4: RMSD PLOT SHOWING THE STRUCTURAL CHANGES IN Aβ DURING MD SIMULATION. A) NATIVE Aβ42 IS HIGHLY UNSTABLE; B) Aβ42 WITH NIAID BECOMES STABLE AFTER 5ns
To examine whether the Aβ central helix eventually unfolded by the end of the simulations, an average number of the αHBs of the Aβ middle region (residues 15-30) was calculated for the last 2ns of the 20ns simulation. The analysis showed the fluctuation of the Aβ backbone RMSD is relatively small in every trajectory. The trajectories were classified into three groups shown in Table 2 group A (RMSD < 2.0Aº, 2 ≤ αHB ≤ 6), group B (2.0 Aº ≤ RMSD < 4.0Aº, 1 ≤ αHB ≤ 4) group C (RMSD ≥2.0 Aº, αHB = 0). On diagram inspection, in group A trajectories it was ascertained that the Aβ central α-helix maintained its helical conformation during the whole simulation or refolded after partial unfolding by the end of the simulations, in group B trajectories partially unfolded at the end of the simulations, and in group C trajectories completely unfolded at the end of the simulations. The helical Aβ (group A) is observed in only two trajectories in the absence of a ligand up to 20ns, whereas it is observed in six trajectories in the presence of NIAID. In contrast, the completely unfolded Aβ (group C) is observed in three trajectories in the absence of a ligand, whereas it is not observed in any trajectory in the present compound. By analyzing the backbone RMSD of the whole simulation of each trajectory, it was found that Aβ helix was relatively stable during simulation in all trajectories in the presence of NIAID. This result indicates that the addition of NIAID compound is effective in stabilizing the Aβ central helix.
TABLE 2: AVERAGE RMSD (Aº) AND AVERAGE NUMBER OF INTRA-MOLECULAR HYDROGEN BONDS DURING LAST 2ns OF 20ns MD SIMULATIONS WAS CALCULATED FOR THE AΒ MIDDLE REGION (RESIDUES 15-30) IN THE ABSENCE AND PRESENCE OF THE NIAID COMPOUND. THE TRAJECTORIES ARE CLASSIFIED INTO THREE GROUPS: GROUP A (RMSD < 2.0 Aº, 2 ≤ αHB ≤ 6), GROUP B (2.0 Aº ≤ RMSD < 4.0 Aº 1 ≤ αHB ≤ 4) GROUP C (RMSD ≥ 2.0 Aº, αHB = 0)
Time (ns) | Average RMSD ( Aº) | Average number of hydrogen bonds in the α-helix of Aβ | Group | |||
Aβ | Aβ + NIAID | Aβ | Aβ + NIAID | Aβ | Aβ + NIAID | |
18.2 | 2.30 | 0.95 | 2.4 | 4.2 | B | A |
18.4 | 2.69 | 1.53 | 2.8 | 3.9 | B | A |
18.6 | 2.75 | 0.81 | 0.6 | 3.5 | C | A |
18.8 | 1.40 | 1.85 | 1.1 | 2.7 | A | B |
19.0 | 1.80 | 1.46 | 1.2 | 4.2 | A | A |
19.2 | 5.30 | 2.21 | 2.2 | 4.8 | B | B |
19.4 | 5.10 | 1.36 | 0.5 | 3.7 | C | A |
19.6 | 3.13 | 2.12 | 3.75 | 4.1 | B | B |
19.8 | 4.82 | 3.15 | 1.46 | 3.4 | B | B |
20.0 | 5.30 | 1.49 | 0.5 | 3.8 | C | A |
Mean Value | 3.23 | 1.7 | 1.65 | 3.83 | - | - |
S.D | 1.25 | 0.68 | 1.09 | 0.56 | - | - |
To inspect whether the ligand was in contact with Aβ(1-42) during MD simulation, the contact maps at different time scale during simulation were analyzed by visual assessment of the trajectories from 2ns to 20ns. All contact maps up to 20ns show ≥ 1 hydrogen bond interaction between Aβ (1-42) peptide and NIAID compound. We found that Aβ (1-42) and NIAID are quite flexible during simulation and forms hydrogen-bond interaction shown in Fig. 5. We have mentioned only three trajectories (10, 15, and 20 ns) among a thousand trajectories. The contact maps showed bonds between two basic functional groups (N1and N2) of NIAID compound with Aβ acidic amino acid residues (E22 and D23). We have analyzed the contacts between the acidic functional groups (O1) of given compounds and the basic residues (H14, K16, and K34) of Aβ middle region and conjugate aromatic hydrocarbon region of NIAID compound contacts with Aβ middle non-polar part localized from L17 to V24. The yellow doted lines in Fig. 5 are hydrogen-bonds confirm that ligand binds with Aβ during the simulation. The amino acid residues of Aβ(1-42) which takes part in hydrogen bonding are shown in Table 3. It is predicted that hydrogen bonds between hydroxyl group of amino acid K16, F19 and D23 bring the test compound and peptide closer. The hydrophobic properties of Aβ amino acids (residues K16L17V18F19F20A21 and V24) are responsible for its unfolding and aggregation. Analysis of simulation trajectories reveals that NIAID compounds forms the hydrogen bond with the respective amino acids (residue16-24) and prevents them to change the α-helix into β-sheet.
FIG. 5: THE RESIDUES INVOLVED IN THE INTRA-MOLECULAR HYDROGEN BOND FORMATION IN THE MIDDLE REGION OF Aβ42 (10-30 RESIDUES) DURING A 20ns MD SIMULATIONS. HYDROGEN BONDING PATTERNS IN Aβ42 WITH NIAID
Torsion angle value was also calculated to observe the alteration of Aβ structure quantitatively during the simulation. The Psi ψ) and Phi (φ) torsion angles as degree freedom provide an efficient way to study the conformations of simulated proteins. The online software ProCheck was used to plot the Ramachandran-graph of native Aβ(1-42) at 0ns Fig. 6A which showed maximum amino acids residues of Aβ are in α-helical conformation only a few of the amino acids residues are in β-sheet. After simulation of native Aβ up to 20ns, amino acids in β-sheet conformation has increased whereas amino acid in α-helical conformation has decreased.
FIG. 6: RAMACHANDRAN PLOT OF THE ψ AND φ DISTRIBUTION PRODUCED BY PROCHECK AFTER MD SIMULATION [A, B, L] MOST FAVOURED REGION [a, b, l, p] ADDITIONAL ALLOWED REGION, [~a, b, ~l, ~ p] GENEROUSLY ALLOWED REGIONS; WHITE ARE DISALLOWED REGIONS. A) NATIVE Aβ42 AT 0ns AND 20ns, B) Aβ42 COMPLEXES WITH NIAID AT 0ns AND 20ns
The Ramachandran-plot of Aβ(1-42) in the presence of NIAID compound Fig. 6B showed after MD simulation only 5 amino acid of Aβ has changed into β-sheet while remaining amino acid of Aβ are in α-helical conformation. In native Aβ, there is a fundamental difference between initial and final torsion angles values after simulations. The difference in torsion angle values is conformation variation in the core structure of Aβ during simulations. During dynamics, NIAID compound has prevented much deviation in Aβ angle values and not allowed the Aβ middle region to unfold and change into β –sheet.
TABLE 3: TABLE SHOWING THE RESIDUES OF Aβ42 INVOLVED IN HYDROGEN BONDING WITH A NIAID TEST COMPOUND. THE HYDROGEN BONDING ANALYSIS WAS PERFORMED FOR 10 STRUCTURES OF Aβ42-LIGAND COMPLEX OBTAINED AT DIFFERENT TIME INTERVALS OF A 20ns MD SIMULATION. CONTRIBUTING SIGNIFICANTLY TO STABILIZE THE Aβ α-HELIX. MOSTLY THE AMINO ACID RESIDUES (K16, F19, F20, AND D23) RESPONSIBLE FOR TO ALLOW Aβ BACKBONE TO UNWIND BINDS WITH TEST COMPOUND WITH HYDROGEN BONDS USING THE CRITERION DISTANCE (≤ 3.0 Aº)
Time (ns) | During dynamics residues in Aβ42 involved in hydrogen bonding with NIAID |
0 | K16, D23 |
2 | D23, L34 |
4 | V12, K16, L34 |
6 | D23, G25 |
8 | H14, D23, L34 |
10 | D23, G37 |
12 | F20, D23, G37 |
14 | L17, F19, L34 |
16 | D23, I32, G37 |
18 | D23, I32, G37 |
20 | D23, G37 |
DISCUSSION AND CONCLUSION: The amino acid sequence of Aβ(11-24) has a strong preference for β-strand structure. The α-helical form of this region of the Aβ peptide provides a suitable target for stabilizing by ligands. The residue of human Aβ E11VHHQKLVFFAEDV24, including the discordant regions is in α-helical conformation 26. The evaluated test ligand NIAID shows hydrogen and hydrophobic chemical interactions with Aβ (1-42) amino acid residues which are responsible for the α-helical conversion into β-sheet. The α-helical propensity of the five amino acid residues (K16L17V18F19 and F20) is the dominating factor for the stability of the Aβ central helix during the molecular dynamics simulation.
The effect of the ligand on the stability of the Aβ(1-42) central helix (residue 11-24) was investigated by using the MD simulations. Detailed information on the structural changes upon loss of helicity in the presence and absence of the ligand was examined.
The Molecular dynamics results signify that the addition of given ligand is efficient in stabilizing the Aβ central helix. The investigation also showed that NIAID compounds bind to the highly elongated form of Aβ and was found to be proficient at forming parallel conformations with Aβ. Thus binding of the test compound to Aβ(1-42) was found to agitate the extension of β-sheet. The MD simulation in presence and absence of the test compounds up to 20ns showed how amino acid of native Aβ loses the helicity at different nano-seconds and also shows test compound NIAID retains the helicity of these amino acids during dynamics.
As indicated mainly by the Aβ backbone RMSD vs. initial structure and by the existence of αHB's of the Aβ. RMSD graph shows fluctuations of the native Aβ backbone, RMSD is relatively small in every trajectory, which infers protein is losing its conformation slowly. The Aβ central helix completely unfolded at the end of the simulation without test compound NIAID, whereas remains folded in the presence of given NIAID by the end of the simulation.
Compared to Aβ alone, the probability of the α -helical state for Aβ during simulations is higher in presence given test compound, and also stability of the Aβ central α-helix was retained. RMSD, αHBs and HB interaction between Aβ and test compound data thus indicate the ability of the compound to stabilize the Aβ central helix. Stabilization of the Aβ central α- helix is an effective step in preventing Aβ polymerization. It has been suggested that complete unfolding occurs via three steps mechanism: Sufficient loss of α-helical backbone hydrogen bonds (αHBs), strong interaction between non-polar side chains and strong interaction between polar side chains 27, 28, 29. Thus, we suggest that test compound NIAID prevent the unfolding of Aβ helix by preventing the breakage of the αHBs and but also disturbing the interactions between polar side chains and between nonpolar side chains.
To understand the polar interactions, the existence of hydrogen bonding between NIAID and Aβ was analyzed. In overall ligand formed as a minimum as one HB with Aβ and the main reason for this is NIAID ligand has several acidic and basic functional groups which can interact to the acidic and basic residues of Aβ respectively. The intermolecular interactions between the Aβ polar residues and the NIAID polar functional groups are significant in stabilizing the Aβ central helix because they can prevent intra-molecular interactions between the Aβ polar residues that stimulate absolute unfolding of the Aβ central helix. Also, benzene-ring of NIAID straddles the Aβ middle non-polar part (residue 17-21). The inter-molecular interactions between the Aβ middle non-polar part and the ligand non-polar part are important in stabilizing the Aβ central helix since the Aβ middle non-polar part includes three non-polar residues (VFF) that have low α-helical propensity and high β-strand propensity 30, 31. Therefore test compounds NIAID by establishing hydrophobic interaction with Aβ stabilize its native formation. Thus could act as a lead molecule to be developed as a drug against AD in preventing fibril formation.
ACKNOWLEDGEMENT: We thank Department of Biotechnology Bharathiar University for providing the computational and Bioinformatics facility (Funded by UGC Government of India).
CONFLICTS OF INTEREST: This article does not contain any studies about human participants or any other animals. The authors R. Boopathy and Dilawar Mir declare that they have no conflict of interest concerning Human and Animal Rights and Informed Consent.
REFERENCES:
- Favaro E: Ageing and Neurologic Disease. Imagining Ageing: Representations of Age and Ageing in Anglophone Literatures 2018; 18: 183.
- Ricciarelli R and Fedele E: The amyloid cascade hypothesis in Alzheimer's disease: it's time to change our mind. Current Neuropharmacology 2017; 15(6): 926-35.
- Kametani F and Hasegawa M: Reconsideration of amyloid hypothesis and tau hypothesis in Alzheimer's disease. Frontiers in Neuroscience 2018; 12: 25.
- Mucke L and Selkoe DJ: Neurotoxicity of amyloid β-protein: synaptic and network dysfunction. Cold Spring Harbor Perspectives in Medicine 2012; 2(7): a006338.
- Lambert MP, Barlow AK, Chromy BA, Edwards C, Freed R, Liosatos M, Morgan TE, Rozovsky I, Trommer B, Viola KL and Wals P: Diffusible, nonfibrillar ligands derived from Aβ1–42 are potent central nervous system neurotoxins. Proceedings of the National Academy of Sciences 1998; 95(11): 6448-53.
- Borlikova GG, Trejo M, Mably AJ, Mc Donald JM, Frigerio CS, Regan CM, Murphy KJ, Masliah E and Walsh DM: Alzheimer brain-derived amyloid β-protein impairs synaptic remodeling and memory consolidation. Neurobiology of Aging 2013; 34(5): 1315-27.
- Michno W, Nyström S, Wehrli P, Lashley T, Brinkmalm G, Guerard L, Syvänen S, Sehlin D, Kaya I, Brinet D and Nilsson KP: Pyroglutamation of amyloid-βx-42 (Aβx-42) followed by Aβ1-40 deposition underlies plaque polymorphism in progressing Alzheimer’s disease pathology. Journal of Biological Chemistry 2019: jbc-RA118.
- Han X and He G: Toward a rational design to regulate β-amyloid fibrillation for Alzheimer's disease treatment. ACS Chemical Neuroscience 2018; 9(2): 198-10.
- Levine III H: Small molecule inhibitors of A beta assembly. Amyloid 2007; 14(3): 185-97.
- Salahuddin P, Siddiqi MK, Khan S, Abdelhameed AS and Khan RH: Mechanisms of protein misfolding: Novel therapeutic approaches to protein-misfolding diseases. Journal of Molecular Structure 2016; 1123: 311-26.
- Re F, Airoldi C, Zona C, Masserini M, Ferla BL, Quattrocchi N and Nicotra F: Beta-amyloid aggregation inhibitors: small molecules as candidate drugs for therapy of Alzheimer's disease. Current Medicinal Chemistry 2010; 17(27): 2990-06.
- Feng BY, Toyama BH, Wille H, Colby DW, Collins SR, May BC, Prusiner SB, Weissman J and Shoichet BK: Small-molecule aggregates inhibit amyloid polymerization. Nature Chemical Biology 2008; 3: 197.
- McGovern SL, Caselli E, Grigorieff N and Shoichet BK: A common mechanism underlying promiscuous inhibitors from virtual and high-throughput screening. Journal of Medicinal Chemistry 2002; 45(8): 1712-22.
- McGovern SL, Helfand BT, Feng B, Shoichet BK. A specific mechanism of nonspecific inhibition. Journal of Medicinal Chemistry 2003; 46(20): 4265-72.
- McKoy AF, Chen J, Schupbach T and Hecht MH: A novel inhibitor of Amyloid β (Aβ) peptide aggregation from high throughput screening to efficacy in an animal model of Alzheimer disease. Journal of Biological Chemistry 2012; 287(46): 38992-00.
- Zhou ZL, Ho Y, Liu HL, Elumalai P and Chen WH: Computer-aided discovery of novel non-peptide inhibitors against amyloid-beta (Aβ) peptide aggregation for treating Alzheimer's disease. Molecular Simulation 2015; 41(8): 622-32.
- Scheidt HA, Morgado I, Rothemund S and Huster D: Dynamics of amyloid β fibrils revealed by solid-state NMR. Journal of Biological Chemistry 2012; 287(3): 2017-21.
- Yan Y and Wang C: Aβ42 is more rigid than Aβ40 at the C terminus: implications for Aβ aggregation and toxicity. Journal of Molecular Biology 2006; 364(5): 853-62.
- Hou L, Lee HG, Han F, Tedesco JM, Perry G, Smith MA and Zagorski MG: Modification of amyloid-β 1-42 fibril structure by methionine-35 oxidation. Journal of Alzheimer's Disease 2013; 37(1): 9-18.
- Economou NJ, Giammona MJ, Do TD, Zheng X, Teplow DB, Buratto SK and Bowers MT: Amyloid β-protein assembly and Alzheimer’s disease: dodecamers of Aβ42, but not of Aβ40, seed fibril formation. Journal of the American Chemical Society 2016; 138(6): 1772-5.
- Jarvet J, Danielsson J, Damberg P, Oleszczuk M and Gräslund A: Positioning of the Alzheimer Aβ (1–40) peptide in SDS micelles using NMR and paramagnetic probes. Journal of biomolecular NMR 2007; 39(1): 63-72.
- Hoyer W, Grönwall C, Jonsson A, Ståhl S and Härd T: Stabilization of a β-hairpin in monomeric Alzheimer's amyloid-β peptide inhibits amyloid formation. Proceedings of the National Academy of Sciences 2008; 105(13): 5099-104.
- Nerelius C, Sandegren A, Sargsyan H, Raunak R, Leijonmarck H, Chatterjee U, Fisahn A, Imarisio S, Lomas DA, Crowther DC and Strömberg R: α-Helix targeting reduces amyloid-β peptide toxicity. Proceedings of the National Academy of Sciences 2009; 106(23): 9191-6.
- Wälti MA, Ravotti F, Arai H, Glabe CG, Wall JS, Böckmann A, Güntert P, Meier BH and Riek R: Atomic-resolution structure of a disease-relevant Aβ (1–42) amyloid fibril. Proceedings of the National Academy of Sciences 2016; 113(34): E4976-84.
- Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK, Goodsell DS and Olson AJ: AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. Journal of Computational Chemistry 2009(16): 2785-91.
- Ito M, Johansson J, Strömberg R and Nilsson L: Unfolding of the amyloid β-peptide central helix: mechanistic insights from molecular dynamics simulations. PloS One 2011; 6(3): e17587.
- Goldschmidt L, Teng PK, Riek R and Eisenberg D: Identifying the amylome, proteins capable of forming amyloid-like fibrils. Proceedings of the National Academy of Sciences 2010; 107(8): 3487-92.
- Berryman JT, Radford SE and Harris SA: Systematic examination of polymorphism in amyloid fibrils by molecular-dynamics simulation. Biophysical Journal 2011; 100(9): 2234-42.
- Gnanakaran S, Nussinov R and García AE: Atomic-level description of amyloid β-dimer formation. Journal of the American Chemical Society 2006; 128(7): 2158-9.
- Goldschmidt L, Teng PK, Riek R and Eisenberg D: Identifying the amylome, proteins capable of forming amyloid-like fibrils. Proceedings of the National Academy of Sciences 2010; 107(8): 3487-92.
- Johansson J, Nerelius C, Willander H and Presto J: Conformational preferences of non-polar amino acid residues: an additional factor in amyloid formation. Biochemical and Biophysical Research Communications 2010; 402(3): 515-8.
How to cite this article:
Mir DA and Rathanam RB: Evaluation of novel ligand for the maintenance of native Aβ(1-42) peptide conformation: relevance to Alzheimer’s disease. Int J Pharm Sci & Res 2019; 10(12): 5489-98. doi: 10.13040/IJPSR.0975-8232.10(12).5489-98.
All © 2013 are reserved by International Journal of Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
32
5489-5498
1018
621
English
IJPSR
D. A. Mir and R. B. Rathanam *
Department of Biotechnology, Bharathiar University, Coimbatore, Tamil Nadu, India.
boopathybiotech@gmail.com
20 March 2019
22 June 2019
17 July 2019
10.13040/IJPSR.0975-8232.10(12).5489-98
01 December 2019