Molecular Docking and Dynamic Simulation of AZD3293 and Solanezumab Effects Against BACE1 to Treat Alzheimer's Disease

The design of novel inhibitors to target BACE1 with reduced cytotoxicity effects is a promising approach to treat Alzheimer's disease (AD). Multiple clinical drugs and antibodies such as AZD3293 and Solanezumab are being tested to investigate their therapeutical potential against AD. The current study explores the binding pattern of AZD3293 and Solanezumab against their target proteins such as β-secretase (BACE1) and mid-region amyloid-beta (Aβ) (PDBIDs: 2ZHV & 4XXD), respectively using molecular docking and dynamic simulation (MD) approaches. The molecular docking results show that AZD3293 binds within the active region of BACE1 by forming hydrogen bonds against Asp32 and Lys107 with distances 2.95 and 2.68 Å, respectively. However, the heavy chain of Solanezumab interacts with Lys16 and Asp23 of amyloid beta having bond length 2.82, 2.78, and 3.00 Å, respectively. The dynamic cross correlations and normal mode analyses show that BACE1 depicted good residual correlated motions and fluctuations, as compared to Solanezumab. Using MD, the Root Mean Square Deviation and Fluctuation (RMSD/F) graphs show that AZD3293 residual fluctuations and RMSD value (0.2 nm) was much better compared to Solanezumab (0.7 nm). Moreover, the radius of gyration (Rg) results also depicts the significance of AZD3293 docked complex compared to Solanezumab through residual compactness. Our comparative results show that AZD3293 is a better therapeutic agent for treating AD than Solanezumab.

Introduction

The β-site APP cleaving enzyme 1 (BACE1) is among the most significant targets for novel drugs to treat Alzheimer's disease (AD) (Huang et al., 2014). The thinning of lipid bilayer has been observed due to protein amyloid-Aβ (Aβ) accumulation and oxidative stress in AD. One study showed that pathologically thin bilayers may play in Aβ aggregation on neuronal bilayers and pathological lipid oxidation may contribute to Aβ misfolding (Korshavn et al., 2017). The aggregation of Aβ peptides results in cellular toxicities due to the formation of polymorphic oligomers, protofibrils, and fibrils (Kotler et al., 2014; Rajasekhar et al., 2016).

Metals accumulation has also deep effect on brain activity and causes neurodegenerations. For example, copper and zinc are associated with the prevalence of AD. In the brains of AD patients, the copper homeostasis is changed with elevated extracellular and low intracellular copper levels. Animals and cell culture studies reported that increasing intracellular copper can cause AD like symptoms through the accumulation of amyloid plaques and tau phosphorylation. Therefore, by changing the copper homeostasis may results in improved cognitive function in animal models of AD (Filiz et al., 2008).

The first generation of BACE1 inhibitors are peptide-based transition state analogs synthesized on the basis of residual sequence of amyloid precursor protein (APP) cleaved by β-secretase (Sinha et al., 1999; Hong et al., 2002; Vassar, 2014). Multiple drugs and antibodies are under different phases of clinical trials to treat AD such as Astra Zenica BACE1 inhibitor (AZD3293) and Solanezumab. AZD3293 is a potent inhibitor against AD and mostly administered orally. AZD3293 prevents the accumulation of β-amyloid and helps slow or cure AD symptoms (Cebers et al., 2017). Clinical studies showed that in both Phase I and II AZD3293 showed good therapeutical potential against AD. Some reports showed that first-generation BACE1 inhibitors were hampered by blood brain barrier (BBB) penetration (Butini et al., 2013; Oehlrich et al., 2014), whereas AZD3293 remarkably has a good penetration against BBB in humans (Cebers et al., 2017).

AZD3293 targets BACE1 and helps reduce Aβ peptide generation (Vassar, 2014). One study showed that AZD3293 has direct effects on plasma and CSF Aβ levels in healthy young and elderly individuals having age 18–55 and 55–80 years, respectively (Alexander et al., 2014). Phase 1 results showed that AZD3293 is well tolerated with no serious adverse effects observed up to the 750 mg in single ascending dose (SAD) study (Alexander et al., 2014). Moreover, a multiple ascending dose (MAD) study reports showed that AZD3293 reduces the CSF Aβ40 and Aβ42 concentrations up to 50 to 75% at 15 or 50 mg doses (Hoglund et al., 2014). Presently, Phase 1 studies of AZD3293 in healthy subjects and in AD patients have been completed, and combined Phase 2/3 trials in 1,551 mild cognitive impairment (MCI) and mild AD patients are planned at different doses from 20 to 50 mg, for 104 weeks duration. Recently, Astra Zenica and Lilly entered into a partnership to jointly develop AZD3293 to treat AD.

Solanezumab (Eli Lilly) is the leading clinical antibody targeting amyloid peptides. Presently, it has been under Phase III clinical trials for the prevention of AD. Solanezumab, a monoclonal antibody IgG1 binds with amyloid-β peptides that aggregate and form plaques in the brain which are considered as basic pathological feature of AD (Villemagne et al., 2013). Solanezumab binds specifically at the monomeric amyloid-β, motif KLVFFAD (Crespi et al., 2015) with pico-molar affinity (Watt et al., 2014). The biding epitope of amyloid-β is known as the nucleation site for Aβ oligomerization, and it is these oligomers of Aβ that are thought to be toxic to neurons. Solanezumab is thought to act as an amyloid beta sink that is “facilitating flux of amyloid beta from a central to peripheral compartment” (DeMattos et al., 2001). Amyloid β-plaques mostly consist of amyloid β-42. Solanezumab binds to free amyloid beta which causes amyloid β-42 to solubilize to re-establish the equilibrium in the cerebrospinal fluid (Farlow et al., 2012). Phase 1–3 results of Solanezumab showed good results against mild to moderate AD.

In the present study, we use docking and dynamic simulation approaches to study the potential of both drugs to cure AD. Research data showed the significance of our proposed methodology in the prediction of drug designing by target various enzymes (Hassan et al., 2017, 2018a,b). We selected 79 drugs and antibodies collectively from clinical drug database (https://clinicaltrials.gov/) and classified them on the basis of active clinical phases. AZD3293 and Solanezumab were selected and further analyzed using various computational tools to calculate their pharmacokinetic properties. Structure-based analyses were performed to evaluate the Root Mean Square Deviation and Fluctuations (RMSD/F), Radius of gyration (Rg) and Solvent Accessible Surface Area (SASA) through dynamic simulation.

Methodology

Repossession of Target Proteins Structure

The human crystal structure of beta-secretase (BACE1) and mid-region amyloid-beta (Aβ) in complex with Solanezumab (PDBIDs: 2ZHV & 4XXD), respectively were retrieved form the Protein Data Bank (PDB) (http://www.rcsb.org). The energy minimization of target proteins was conducted by using online tool Chiron to resolve the steric clashes from protein structures (Ramachandran et al., 2011). The stereo-chemical properties, Ramachandran graph and values (Lovell et al., 2003) of targeted proteins were assessed by MolProbity server (Chen et al., 2010), whereas the hydrophobicity graphs were generated by Discovery Studio 4.1 Client (Studio Discovery, 2008). The protein architecture and statistical percentage values of helices, β-sheets, coils and turns were accessed by using the online server VADAR 1.8 (Willard et al., 2003).

Candidate Structure

The AZD3293 and Solanezumab were selected as clinical drugs for present study. The AZD3293 drug was sketched in drawing ACD/ChemSketch tool. 79 drugs and antibodies were selected from clinical drug database (https://clinicaltrials.gov/) and classified on the basis of active clinical phases (Table S1). The selected drug molecule further minimized by UCSF Chimera 1.10.1 (Pettersen et al., 2004). Multiple online drug assessment computational tools such as Molinspiration (http://www.molinspiration.com/) and Molsoft (http://www.molsoft.com/) were used to predict the drug-likeness and biological properties of AZD3293 molecule. Lipinski's rule of five was analyzed using Molsoft and Molinspiraion tools. Furthermore, their predicted Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties were evaluated by pkCSM online tool (Pires et al., 2015). Solanezumab is a monoclonal IgG1 antibody which is directed against Aβ peptide. The three dimensional (3D) structure of Solanezumab was also accessed from protein data bank as mentioned above.

Molecular Docking Simulation

Molecular docking of AZD3293 drug against BACE1 was carried out using diverse AutoDock 4.2 tool according to the specified instructions (Morris et al., 2009). In brief, for receptor protein, the polar hydrogen atoms and Kollman charges were assigned. For ligand, Gasteiger partial charges were designated and non-polar hydrogen atoms were merged. All the torsion angles for AZD3293 were set free to rotate through docking experiment. A grid map of 80 × 80 × 80 Å was adjusted on whole protein structure to generate the grid map and to get the best conformational state of docking. The 100 number of runs were adjusted using docking experiments. The Lamarckian genetic algorithm (LGA) and empirical free energy function were applied by taking docking parameters default. All the docked complexes were further evaluated on lowest binding energy (Kcal/mol) values and hydrogen and hydrophobic interactions analysis using Discovery Studio (2.1.0) and UCSF Chimera 1.10.1. The two dimensional graphical depiction of best docked complexes were accessed by LIGPLOT tool (Wallace et al., 1995).

Molecular Dynamics Simulations

Based on docking results, we performed structural dynamic analysis studies on the selected lowest energy valued and best posed docking complex. The MD simulations were carried out by Groningen Machine for Chemicals Simulations (GROMACS) 4.5.4 package (Pronk et al., 2013) with GROMOS 96 force field (Chiu et al., 2009). The PRODRG Server was employed to generate ligands topology files (Schüttelkopf and van Aalten, 2004). Before minimization, the overall system charge was neutralized by adding ions. The energy minimization (nsteps = 50,000) was conducted using the steepest descent approach (1,000 ps) for each protein-ligand complex. The Particle Mesh Ewald (PME) method was employed for energy calculation and for electrostatic and Van der Waals interactions; cut-off distance for the short-range VdW (rvdw) was set to 14 Å, where Coulomb cut-off (r coulomb) and neighbor list (rlist) were fixed at 9 Å (Wang et al., 2010). It permits the use of the Ewald summation at a computational cost comparable with that of a simple truncation method of 10 Å or less, and the linear constraint solver (LINCS) (Amiri et al., 2007) algorithm was used for covalent bond constraints and the time step was set to 0.002 ps. Finally, 20 ns molecular dynamics simulation was carried out for all the complexes with nsteps 1,000,000. The Root Mean Square Deviation and fluctuation (RMSD/F), Soluble Accessible Surface Area (SASA) and Radius of gyration (Rg) analysis were carried out using Xmgrace (http://plasma-gate.weizmann.ac.il/Grace/) and UCSF Chimera 1.10.1 software.

Results and Discussion

Structural Assessment of BACE1, Aβ, and Solanezumab

BACE1 is a class of hydrolase protein comprises having single chain and 411 amino acids. The structural architecture of BACE1 showed that it consists of 9% helices, 49% β sheets, 41% coil, and 21% turns. The X-Ray diffraction study confirmed its resolution 1.85 Å, R-value 0.290 and unit cell crystal dimensions like length and angles of coordinates for a = 102.33, b = 102.33, and c = 170.52 with angles 90, 90, and 120° and for all α, β, and γ dimensions respectively. The Ramachandran plots and values indicated that 98% of residues were in favored regions and >99.8% residues were present in allowed regions. Ramachandran graph and values showed the good accuracy of phi (ϕ) and psi (ψ) angles among the coordinates of receptor molecules and most of residues plummeted in acceptable region. The Aβ is clump of multiple chains (A, B, C, D, and E) which contain 42 amino acids.

The 3D crystal structure of Solanezumab showed that it consists of two Fab light chains (A, D) and two Fab heavy chains (B, E) having 219 and 223 amino acids, respectively. The X-Ray diffraction study indicated its resolution 2.41Å and R-value 0.290. The unit cell crystal dimensions such as length and angles of coordinates for a = 38.80, b = 73.56, and c = 92.12 with angles α = 109.91°, β = 93.64°, and γ = 93.31° for all dimensions respectively. Moreover, the VADAR 1.8 analysis showed that Solanezumab contain 1% helices, 59% β sheets, 38% coil, and 4% turns. Furthermore, Ramachandran analysis depicts the 98% residues are present in favored region while 9.8% amino acids are lie in allowed region. The graphical depiction of Ramachandran plots of BACE1 and Solanezumab are mentioned in detail in Supplementary Data (Figures S2, S3), respectively. The detailed structural analysis of BACE1, Aβ and Solanezumab are discussed in Table S4.

Computational Evaluation of AZD3293

Chemoinformatics Properties and Lipinski's Rule (RO5) Validation of AZD3293

Multiple computational approaches were employed to predict basic chemoinformatics and basic molecular properties of AZD3293 (Figure 1). The predicted properties such as molecular weight (g/mol) polar surface area (PSA, A²), molar volume (cm³), density (cm³), molar refractivity (cm³) and RO5 were evaluated to justify their drug likeness behavior (Prerana et al., 2015). One study revealed that PSA is a significant parameter for drug absorption prediction in drug discovery (Ertl et al., 2000). The molar refractivity and molecular lipophilicity properties of drug molecules are also important for receptor binding, bioavailability and cellular uptake within the body. One report justifies the standard values for molar refractivity (40–130 cm3) and molecular weight (160 to 480 g/mol) and PSA (<89 Å2) (Ghose et al., 2012). Moreover, the qualifying range for total number of atoms in the drug molecule is between 20 and 70 atoms. Table 1 results showed that molar refractivity and PSA predicted values were comparable with standard values. AZD3293 showed higher molar refractivity and logP value (122.60 cm³ and 4.82) respectively, compared to standard values. Comparative results showed that AZD3293 confirm its significant and good candidate molecule.

FIGURE 1

Figure 1. Chemical structure of AZD3293.

TABLE 1

Table 1. Chemoinformatic properties of AZD3293.

Furthermore, AZD3293 was validated by RO5 and contains no more and < 10 hydrogen bond acceptors (HBA) and 5 (hydrogen bond acceptors) HBD, respectively. Moreover, the logP and molecular mass value also be <5 and 500 (g/mol), respectively. Literature study revealed that the exceed values of HBA and HBD results in poor permeation (Kadam and Roy, 2007). The hydrogen bonding ability has been considered a significant parameter for drug permeability. Our results justified that the AZD3293 possess < 10 HBA, < 5 HBD, < 500 (g/mol) molecular weight and <5 logP values which were comparable with standard values. The reported study showed that molecules with poor absorption are more likely to be observed upon Lipinski violation. However, multiple examples are available for RO5 violation amongst the existing drugs (Bakht et al., 2010; Tian et al., 2015). The predicted drug score (0.40) and bioactivity score values are also significant for further analysis. The predicted score values of G-protein couple receptor (GPCR) (0.51), protease and enzymes inhibition score (0.53 and 0.56), respectively showed their good lead like behavior.

Pharmacokinetic Properties of AZD3293

The designing of novel drugs require a high attention rate with good pharmacokinetic properties. The Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties were assessed to confirm the efficacy of candidate molecules. In ADMET evaluation, absorption like water and intestinal solubility (log mol/L & % absorbed) and skin permeability (logKp) predicted values justified the strong therapeutic potential of chemical compounds. One report justified that compounds with good absorption values have potency to cross gut barrier by passive penetration to reach the target molecule (Selick et al., 2002). The water solubility results justified that AZD3293 showed good absorption value (−4.956 log mol/L). Moreover, the intestinal solubility prediction value (96.90) also justified its good efficacy compared to a standard value (>30% abs). Any chemical lead like structure with <30% absorbance value is considered as poorly absorbed compound (Pires et al., 2015). The predicted skin permeability value (−2.902 log Kp) of AZD3293 was also comparable with standard value (−2.5 logKp) which showed their significance as a good lead structures and justified their drug likeness behavior. The p-glycoprotein inhibition behavior was also confirmed for AZD3293. Moreover, in distribution properties, the Blood Brain Barrier (BBB) and Central Nervous System (CNS) permeability values of AZD3293 were also evaluated and compared with the standard values (>0.3 to < −1 log BB and >−2 to < −3 logPS) respectively. It has been observed that compounds with a >0.3 log BB value have potential to cross BBB, while with < −1 value are poor distributed to brain. The predicted results showed that AZD3293 have poor BBB value (−0.164 log BB). However, the CNS permeability value (−1.72 log PS) is quite comparable with standard value. Similarly, the compounds have > −2 logPS value are considered to penetrate the CNS, while with < −3 are difficult to move in the CNS.

Moreover, metabolic behavior of AZD3293 was confirmed by CYP3A4, which is isoform of cytochrome P450. The excretion and toxicity predicted values were also justified the drug likeness behavior of AZD3293 on the basis of total clearance (log ml/min/kg), AMES toxicity, maximum tolerated dose (MTD) and LD₅₀ values. The AMES toxicity prediction for AZD3293 also confirmed there is non-mutagenic and non-toxic behavior. The hepatotoxicity positive effect showed its lethal behavior while skin sensitivity negative behavior presents their non-toxic and less sensitive effects. Disruption of normal liver function is commonly associated with hepatotoxicity (Table 2).

TABLE 2

Table 2. Pharmacokinetic assessment of AZD3293.

Molecular Docking Analysis

AZD3293 Binding Energy Analysis Against BACE1

The AZD3293-BACE1 docked complexes were analyzed on the basis of lowest binding energy values (Kcal/mol) and hydrogen/hydrophobic interaction analyses. The best pose selection from all the docking complexes was also conducted on the basis of lowest binding energy values and bonding interaction pattern within the active region of target protein. Results showed that AZD3293 with pose 8 was the most active conformational position and predicts the best energy values (−5.33 kcal/mol) as compared to others docking complexes (Table 3). Furthermore, the intermolecular energy value (−7.42 Kcal/mol) was also good compared to other docking pose conformations. The Autodock energy calculation was conducted by using equation 1. However, the standard energy error for Autodock is 2.5 Kcal/mol (Morris et al., 2013). In all docking complexes poses the binding energy difference was less 2.5 Kcal/mol.

TABLE 3

Table 3. Docking energy of AZD3293 against BACE1.

ΔGbinding = ΔGgauss + ΔGrepulsion + ΔGhbond + ΔGhydrophobic + ΔGtors…..(i)

Here, ΔG gauss: attractive term for dispersion of two gaussian functions, ΔGrepulsion: square of the distance if closer than a threshold value, ΔGhbond: ramp function-also used for interactions with metal ions, ΔGhydrophobic: ramp function, ΔGtors: proportional to the number of rotatable bonds.

Binding Conformational Analysis

The best docked energy complexes were further deep analyzed on the basis of hydrogen and hydrophobic interactions pattern between ligand and target protein. The active binding region of BACE1 was accessed through one study (Hernández-Rodríguez et al., 2016). Results showed that AZD3293 perfectly binds within the active region of target protein by forming couple of hydrogen bonds. The graphical depiction of best docking complex along with receptor and active binding pocket region is mentioned in Figure 2.

FIGURE 2

Figure 2. (A) The BACE1 3D structure in gray color with interior yellow color having surface format. (B) The close depiction of active region of target protein with embedded AZD3293. The AZD3293 is mentioned in purple color while the functional groups such as amino and oxygen are highlighted with blue and red colors, respectively. (C) Docking interaction between AZD3293 and target protein (BACE1). The target protein is highlighted in line ribbon format with gray color. The active binding amino acids are highlighted in dark maroon color around the ligand. Two hydrogen bonds were observed between AZD3293 and BACE1 amino acids at Asp32 and Lys107 with distances 2.95 and 2.68 Å, respectively. The black dotted lines show the binding distance in angstrom (Å).

The structure activity relationship (SAR) analysis shows that AZD3293 forms two hydrogen bonds at specific residues (Asp32 and Lys107) with target protein. The hydroxyl group of benzene ring of AZD3293 interacts with Asp32 having bonds length 2.95 Å. Similarly, amino group of other AZD3293 forms another hydrogen bond with bond length 2.68 Å. One study also justified that these interacted residues are significant in the downstream signaling pathways (Hernández-Rodríguez et al., 2016). The graphical representations of all other docking poses are mentioned in Supplementary Data (Figures S5–S13).

Binding Analysis of Solanezumab With Aβ

Solanezumab is a monoclonal IgG1 antibody directed against Aβ peptide. Solanezumab exerts its effect by sequestering Aβ, shifting equilibria between different species of Aβ, and removing small soluble species of Aβ that are directly toxic to synaptic function. The binding interaction pattern was observed to residual involvement of both proteins. The results showed that Fab region of both light and heavy chains of Solanezumab collapse with Aβ peptide. However, the observed interacted residue were present at heavy chain. The residues of Aβ peptide such as Lys16 and Asp23 are directly involved in hydrogen bonding with heavy chain of Asp96 and Ser33 respectively. The Lys16 forms hydrogen bond with Asp96 with distance 2.82 Å. While Ser33 forms two hydrogen bonds with Asp23 with bond length 2.78 and 3.00 Å, respectively (Figure 3). The Aβ peptide region is considered as most active part and multiple hydrophobic interaction were also observed against Solanezumab. Literature report also justified that peptide region from Lys16 to Val24 is most significant in the binding with Solanezumab and further downstream signaling pathways (Crespi et al., 2015). The hydrophobic interacted residues of Solanezumab are as follows: Phe27, Phe55, Phe36, Ser94, Gly95, His34, Tyr32, Lys50, Asp28, Thr92, and Val94 which binds with Aβ peptide. The hydrophobic interacted residues of both Solanezumab and Aβ peptide are mentioned in Supplementary Data (Figure S14).

FIGURE 3

Figure 3. The binding interaction of Solanezumab with Aβ. (A) The complex structure of Solanezumab with Aβ. The light chain is represented in purple color while heavy chain is highlighted in green color. Aβ is depicted in red color at Fab region of Solanezumab. The whole structure is represented in surface format. (B) The epitope region residues of Solanezumab of both light and heavy chain are represented in purple and green colors, respectively. The Aβ interacted residues are represented in red color while the black dotted lines are for bonding distance in angstrom (Å).

Dynamical Cross-Correlation Matrix (DCCM) of BACE1 and Solanezumab

Bio3D-web is an online application to analyze the sequence, structure and conformational heterogeneity of protein families (Skjærven et al., 2016). The residual fluctuations of target protein structures were analyzed using Bio3D server. The generated dynamical cross-correlation graphs depicts positive and negative correlation effect of amino acids. The pairwise correlated graphs were generated between the function of residue indices i and j. The predicted map results were analyzed on the basis of colors such as dark cyan, white and pink, respectively. The fully correlated pairs are represented by cyan color, while the anti-correlated are justified as pink color. However, the moderately and un-correlated regions are highlighted by yellow and cyan, respectively. The comparative results showed that BACE1 depicted good residual correlated motions as compared Solanezumab complex (Figure 4).

FIGURE 4

Figure 4. Dynamical cross-correlation analysis of both target structures. (A) The residual cross-correlation analysis of BACE1 (PDBID: 2ZHV). (B) The residual cross-correlation analysis of solanezumab antibody (PDBID: 4XXD).

Normal Mode Analysis (NMA) Fluctuations Analysis BACE1 and Solanezumab

Bio-3D NMA web analysis displays the dynamic behavior of various residues in the protein structures (Yao et al., 2016). The predicted graphs results showed that BACE1 has less residual fluctuations compared to Solanezumab antibody complex. The comparative analysis reveals that BACE1 residual fluctuation range start from 0 to 1.0, while Solanezumab has 0–2.5. The N-terminus region of BACE1 has less amino acids fluctuations compared to C-terminus. Whereas, in Solanezumab whole structure is fluctuated (Figure 5).

FIGURE 5

Figure 5. Ensemble normal mode analysis. (A) The residual fluctuations analysis of BACE1 (PDBID: 2ZHV). (B) The residual fluctuations analysis of solanezumab antibody (PDBID: 4XXD). The fluctuations peaks in both graphs less than <0.05 showed the accuracy of crystal structural behavior.

Molecular Dynamic Simulations Analysis

Root Mean Square Deviation and Fluctuation of Target Structures

To evaluate the flexibility and overall stability of docked complexes, we conducted time dependent MD simulation at 20 ns using Gromacs 4.5.4. The residual deviations and fluctuation in the complexes were determined by using RMSD and RMSF graphs generated by using Xmgrace software. Figures 6, 7 exhibited the residual deviation and fluctuations of AZD3293 and Solanezumab docked complexes, respectively. The increasing trend was observed in both (AZD3293 and Solanezumab) complexes having diverse RMSD values 0–0.25 and 0.5 nm at equilibrium state (starting) from 0 to 2,500 ps in the simulation period. Initially, the AZD3293 complex shows little variations while Solanezumab depicted higher fluctuations. From 2,500 to 5,000 ps, AZD3293 complex indicates little fluctuations while Solanezumab exhibited more deviations with increased value of RMSD value 0.6 nm. From 5,000 to 20,000 ps AZD3293 still remain static depiction with constant RMSD value 0.25 nm while Solanezumab continuously increases with enhanced RMSD value. The overall MD results showed that AZD3293 fluctuations was much better as compared to Solanezumab in the simulation time. The predicted graph suggested that AZD3293 complex is much better and little fluctuated throughout the simulation period. The RMSF results of both AZD3293 and Solanezumab reflect the loops fluctuations throughout the simulation period. The comparative analyses showed that Solanezumab exhibited higher fluctuated peaks as compared to AZD3293 which may depict the significance of AZD3293 over Solanezumab. Insight from MD simulations stable behaviors of AZD3293 docked complex throughout MD trajectories thus increasing the efficacy of docking results.

FIGURE 6

Figure 6. RMSD graphs of AZD3293 and Solanezumab at 20 ns. The graph lines with red and purple represents AZD3293 and Solanezumab complexes, respectively.

FIGURE 7

Figure 7. RMSF graphs of AZD3293 and Solanezumab at 20 ns. The graph lines with red and purple represents AZD3293 and Solanezumab complexes, respectively.

Radius of Gyration Analyses

The compactness of protein is measured by radius of gyration (Rg). The predicted results showed that AZD3293 showed much static depiction and constant Rg value at 2.05 nm throughout the simulation time 0–20,000 ps. Whereas, Solanezumab depicts variations at initiating while later shows stability with Rg value 2.5 nm. The comparative results justified that the residual backbone and proper conformation of AZD3293 docked receptor is much batter compared to Solanezumab complex (Figure 8).

FIGURE 8

Figure 8. Radius of gyration (Rg) graph of AZD3293 and Solanezumab at 20 ns. The graph lines with red and purple represents AZD3293 and Solanezumab complexes, respectively.

The comparative analyses showed that AZD3293 has better therapeutical potential than Solanezumab. Our analyses explore the residual fluctuations after binding the drug with BACE1 receptors. The predicted results depicted the significance of AZD3293 over Solanezumab in treating AD. One limitation with MD simulations is that the protein covalent bonds remain unbroken which may disturb the rmsd and rmsf graph fluctuations. However, there are few things that can done in future work to enhance the accuracy and reliability of MD simulation, such as improvement of simulation time, MD force field, and novel algorithms.

In the future, a molecular dynamics simulation study may be used as good computational approach to study the interactions of Aβ40 and Aβ42 with model neuronal membranes. The higher running time (i.e., 200–300 ns) scale with best force field (i.e., a collection of equations and associated constants designed to reproduce molecular geometry and selected properties of tested structures) simulation can help investigate the effect of Aβ peptides on neuronal membrane. Moreover, some key factors should also be investigated such as surface area per lipid, bilayer thickness and lipid order parameter to obtain the reliable results. The advancement in computer hardware can further mature the simulation method.

Conclusion

Computational interpretation of newly designed compounds and their binding analysis in the active region of target proteins allow pharmaceutical industries and labs to test the efficacy of drugs before starting experimental lab work. In our in-silico approach, we compared two clinical drugs (AZD3293 and Solanezumab) are being testing against AD. Multiple online drug analysis computational tools and server were employed to predict the efficacy and possible lethality of selected compounds. Molsoft and molinspiration results depict the validity of RO5 for AZD3293. Our results show that AZD3293 may be considered as a good therapeutic agent and have good drug like behavior. Moreover, their pharmacokinetic properties such as ADMET properties also justified their good lead like behavior and drug potential. In all parameters such as Absorption, Distribution and Excretion, AZD3293 gave positive results except little lethality in Metabolism and Toxicity. AZD3293 depicts its hepatotoxicity behavior which may cause serious effect in the body. Furthermore, the docking based hydrogen binding and structural analyses show the significance of both selected drugs with proper conformation within the active region of target proteins. The cross correlation and normal mode results depict the structural stability of receptors molecules. The comparative analyses showed that AZD3293-receptor complex has more stable behavior compared to Solanezumab. Finally, a detail simulation study was performed to get the deeper insight of backbone fluctuations and structural stability through RMSD/F results. The generated RMSD/F graphs results showed that the AZD3293-BACE1 docked complex predicts stable behavior compared to Solanezumab regarding the simulation time. The comparative RMSD value (0.2 nm) of AZD3293 was more significant compared to Solanezumab (0.7 nm). The radius of gyration (Rg) results also represent the proper residual conformation and compactness AZD3293 receptor compared to Solanezumab. Based on our computational and comparative results, we may conclude the significance of AZD3293 over Solanezumab on the basis of good binding and structural stability behavior against target proteins of AD.

Author Contributions

MH conceived the study under the guidelines of SYS and AM. MH and SS collected data and performed the experimental work. MH wrote the initial draft of manuscript. HA and NZ make corrections in initial draft. SYS and AM edited the manuscript and compiled it into final format.

Funding

This work was supported by the research grant of the Kongju National University in 2014.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2018.00034/full#supplementary-material

References

Alexander, R., Budd, S., Russell, M., Kugler, A., Cebers, G., Ye, N., et al. (2014). AZD3293 a novel bace1 inhibitor: safety, tolerability and effects on plasma and CSF Ab peptides following single-and multiple-dose administration. Neurobiol. Aging 35:S2. doi: 10.1016/j.neurobiolaging.2014.01.033

CrossRef Full Text

Amiri, S., Sansom, M. S., and Biggin, P. C. (2007). Molecular dynamics studies of AChBP with nicotine and carbamylcholine: the role of water in the binding pocket. Protein Eng. Des. Sel. 20, 353–359. doi: 10.1093/protein/gzm029

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakht, M. A., Yar, M. S., Abdel-Hamid, S. G., Al Qasoumi, S. I., and Samad, A. (2010). Molecular properties prediction, synthesis and antimicrobial activity of some newer oxadiazole derivatives. Eur. J. Med. Chem. 45, 5862–5869. doi: 10.1016/j.ejmech.2010.07.069

PubMed Abstract | CrossRef Full Text | Google Scholar

Butini, S., Brogi, S., Novellino, E., Campiani, G., Ghosh, A. K., Brindisi, M., et al. (2013). The structural evolution of beta-secretase inhibitors: a focus on the development of small-molecule inhibitors. Curr. Top. Med. Chem. 13, 1787–1807. doi: 10.2174/15680266113139990137

PubMed Abstract | CrossRef Full Text | Google Scholar

Cebers, G., Alexander, R. C., Haeberlein, S. B., Han, D., Goldwater, R., Ereshefsky, L., et al. (2017). AZD3293: pharmacokinetic and pharmacodynamic effects in healthy subjects and patients with Alzheimer's disease. J. Alzheimers. Dis. 55, 1039–1053. doi: 10.3233/JAD-160701

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, V. B., Arendall, W. B. III., Headd, J. J., Keedy, D. A., Immormino, R. M., Kapral, G. J., et al. (2010). MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. D Biol. Crystallogr. 66(Pt 1), 12–21. doi: 10.1107/S0907444909042073

PubMed Abstract | CrossRef Full Text | Google Scholar

Chiu, S. W., Pandit, S. A., Scott, H. L., and Jakobsson, E. (2009). An improved united atom force field for simulation of mixed lipid bilayers. J. Phys. Chem. B 113, 2748–2763. doi: 10.1021/jp807056c

PubMed Abstract | CrossRef Full Text | Google Scholar

Crespi, G. A., Hermans, S. J., Parker, M. W., and Miles, L. A. (2015). Molecular basis for mid-region amyloid-β capture by leading Alzheimer's disease immunotherapies. Sci. Rep. 5:9649. doi: 10.1038/srep09649

PubMed Abstract | CrossRef Full Text | Google Scholar

DeMattos, R. B., Bales, K. R., Cummins, D. J., Dodart, J. C., Paul, S. M., and Holtzman, D. M. (2001). Peripheral anti-A beta antibody alters CNS and plasma A beta clearance and decreases brain A beta burden in a mouse model of Alzheimer's disease. Proc. Natl. Acad. Sci. U.S.A. 98, 8850–8855. doi: 10.1073/pnas.151261398

PubMed Abstract | CrossRef Full Text | Google Scholar

Ertl, P., Rohde, B., and Selzer, P. (2000). Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J. Med. Chem. 43, 3714–3717. doi: 10.1021/jm000942e

PubMed Abstract | CrossRef Full Text | Google Scholar

Farlow, M., Arnold, S. E., van Dyck, C. H., Aisen, P. S., Snider, B. J., Porsteinsson, A. P., et al. (2012). Safety and biomarker effects of solanezumab in patients with Alzheimer's disease. Alzheimers Dement. 8, 261–271. doi: 10.1016/j.jalz.2011.09.224

PubMed Abstract | CrossRef Full Text | Google Scholar

Filiz, G., Price, K. A., Caragounis, A., Du, T., Crouch, P. J., and White, A. R. (2008). The role of metals in modulating metalloprotease activity in the AD brain. Eur. Biophys. J. 37, 315–321. doi: 10.1007/s00249-007-0244-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghose, A. K., Herbertz, T., Hudkins, R. L., Dorsey, B. D., and Mallamo, J. P. (2012). Knowledge-Based, Central Nervous System (CNS) lead selection and lead optimization for CNS drug discovery. ACS Chem. Neurosci. 3, 50–68. doi: 10.1021/cn200100h

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, M., Abbas, Q., Ashraf, Z., Moustafa, A. A., and Seo, S. Y. (2017). Pharmacoinformatics exploration of polyphenol oxidases leading to novel inhibitors by virtual screening and molecular dynamic simulation study. Comput. Biol. Chem. 68, 131–142. doi: 10.1016/j.compbiolchem.2017.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, M., Azhar, M., Abbas, Q., Raza, H., Moustafa, A. A., Shahzadi, S., et al. (2018a). Finding novel anti-carcinomas compounds by targeting SFRP4 through molecular modeling, docking and dynamic simulation studies. Curr. Comput. Aided. Drug Des. 14, 160–173. doi: 10.2174/1573409914666180112100122

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, M., Ashraf, Z., Abbas, Q., Raza, H., and Seo, S. Y. (2018b). Exploration of novel human tyrosinase inhibitors by molecular modeling, docking and simulation studies. Interdiscip. Sci. 10, 68–80. doi: 10.1007/s12539-016-0171-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernández-Rodríguez, M., Correa-Basurto, J., Gutiérrez, A., Vitorica, J., and Rosales-Hernández, M. C. (2016). Asp32 and Asp228 determine the selective inhibition of BACE1 as shown by docking and molecular dynamics simulations. Eur. J. Med. Chem. 124, 1142–1154. doi: 10.1016/j.ejmech.2016.08.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoglund, K., Salter, H., Zetterberg, H., Andreason, U., Olsson, T., Alexander, R., et al. (2014). Monitoring the Soluble Amyloid Precursor Protein Alpha (SAPPA) and Beta (SAPPB) fragments in plasma and CSF from healthy individuals treated with bace inhibitor AZD3293 in a multiple ascending dose study: Pharmacokinetic and Pharmacodynamic correlate. Alzheimers. Dement. 10:P447. doi: 10.1016/j.jalz.2014.05.605

CrossRef Full Text | Google Scholar

Hong, L., Turner, R. T. III., Koelsch, G., Shin, D., Ghosh, A. K., and Tang, J. (2002). Crystal structure of memapsin 2 (beta-secretase) in complex with an inhibitor OM00-3. Biochemistry 41, 10963–10967. doi: 10.1021/bi026232n

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H. J., Lee, C. C., and Chen, C. Y. (2014). In silico design of BACE1 inhibitor for Alzheimer's disease by traditional Chinese medicine. Biomed Res. Int. 2014:741703. doi: 10.1155/2014/741703

PubMed Abstract | CrossRef Full Text | Google Scholar

Kadam, R. U., and Roy, N. (2007). Recent trends in drug-likeness prediction: a comprehensive review of in silico methods. Indian J. Pharm. Sci. 69, 609–615 doi: 10.4103/0250-474X.38464

CrossRef Full Text | Google Scholar

Korshavn, K. J., Satriano, C., Lin, Y., Zhang, R., Dulchavsky, M., Bhunia, A., et al. (2017). Reduced lipid bilayer thickness regulates the aggregation and cytotoxicity of amyloid-β. J. Biol. Chem. 292, 4638–4650. doi: 10.1074/jbc.M116.764092

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotler, S. A., Walsh, P., Brender, J. R., and Ramamoorthy, A. (2014). Differences between amyloid-β aggregation in solution and on the membrane: insights towards elucidation of the mechanistic details of Alzheimer's disease. Chem. Soc. Rev. 43, 6692–6700. doi: 10.1039/c3cs60431d

PubMed Abstract | CrossRef Full Text | Google Scholar

Lovell, S. C., Davis, I. W., Arendall, W. B. III., de Bakker, P. I., Word, J. M., Prisant, M. G., et al. (2003). Structure validation by C-alpha geometry: phi, psi and C beta deviation. Proteins 50, 437–450. doi: 10.1002/prot.10286

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, G. M., Green, L. G., Radić, Z., Taylor, P., Sharpless, K. B., Olson, A. J., et al. (2013). Automated docking with protein flexibility in the design of femtomolar “click chemistry” inhibitors of acetylcholinesterase. J. Chem. Inf. Model. 53, 898–906. doi: 10.1021/ci300545a

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and autodocktools4: automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi: 10.1002/jcc.21256

PubMed Abstract | CrossRef Full Text | Google Scholar

Oehlrich, D., Prokopcova, H., and Gijsen, H. J. (2014). The evolution of amidine-based brain penetrant BACE1 inhibitors. Bioorg. Med. Chem. Lett. 24, 2033–2045. doi: 10.1016/j.bmcl.2014.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Pettersen, E. F., Goddard, T. D., Huang, C. C., Couch, G. S., Greenblatt, D. M., Meng, E. C., et al. (2004). UCSF chimera – a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612. doi: 10.1002/jcc.20084

PubMed Abstract | CrossRef Full Text | Google Scholar

Pires, D. E., Blundell, T. L., and Ascher, D. B. (2015). pkCSM: predicting small-molecule pharmacokinetic properties using graph-based signatures. J. Med. Chem. 58, 4066–4072. doi: 10.1021/acs.jmedchem.5b00104

PubMed Abstract | CrossRef Full Text | Google Scholar

Prerana, B. J., Akshay, R. Y., and Megha, G. G. (2015). Concept of drug likeness in pharmaceutical research. Int. J. Pharm. Bio. Sci. 6, 142–154.

Google Scholar

Pronk, S., Páll, S., Schulz, R., Larsson, P., Bjelkmar, P., Apostolov, R., et al. (2013). GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 29, 845–854. doi: 10.1093/bioinformatics/btt055

PubMed Abstract | CrossRef Full Text | Google Scholar

Rajasekhar, K., Narayanaswamy, N., Murugan, N. A., Kuang, G., Ågren, H., and Govindaraju, T. (2016). A high affinity red fluorescence and colorimetric probe for amyloid β aggregates. Sci. Rep. 6:23668. doi: 10.1038/srep23668

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramachandran, S., Kota, P., Ding, F., and Dokholyan, N. V. (2011). Automated minimization of steric clashes in protein structures. Proteins 79, 261–270. doi: 10.1002/prot.22879

PubMed Abstract | CrossRef Full Text | Google Scholar

Schüttelkopf, A. W., and van Aalten, D. M. (2004). PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr. D Biol. Crystallogr. 60(Pt 8), 1355–1363. doi: 10.1107/S0907444904011679

PubMed Abstract | CrossRef Full Text | Google Scholar

Selick, H. E., Beresford, A. P., and Tarbit, M. H. (2002). The emerging importance of predictive ADME simulation in drug discovery. Drug Discov. Today 7, 109–116. doi: 10.1016/S1359-6446(01)02100-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sinha, S., Anderson, J. P., Barbour, R., Basi, G. S., Caccavello, R., Davis, D., et al. (1999). Purification and cloning of amyloid precursor protein beta-secretase from human brain. Nature 402, 537–540. doi: 10.1038/990114

PubMed Abstract | CrossRef Full Text | Google Scholar

Skjærven, L., Jariwala, S., Yao, X. Q., and Grant, B. J. (2016). Online interactive analysis of protein structure ensembles with Bio3D-web. Bioinformatics 32, 3510–3512. doi: 10.1093/bioinformatics/btw482

PubMed Abstract | CrossRef Full Text | Google Scholar

Studio Discovery (2008). “version 2.1.” San Diego, CA: Accelrys.

Tian, S., Wang, J., Li, Y., Li, D., Xu, L., and Hou, T. (2015). The application of in silico drug-likeness predictions in pharmaceutical research. Adv. Drug Deliv. Rev. 86, 2–10. doi: 10.1016/j.addr.2015.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Vassar, R. (2014). BACE1 inhibitor drugs in clinical trials for Alzheimer's disease. Alzheimers. Res. Ther. 6:89. doi: 10.1186/s13195-014-0089-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Villemagne, V. L., Burnham, S., Bourgeat, P., Brown, B., Ellis, K. A., Salvado, O., et al. (2013). Amyloid β deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer's disease: a prospective cohort study. Lancet Neurol. 12, 357–367. doi: 10.1016/S1474-4422(13)70044-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wallace, A. C., Laskowski, R. A., and Thornton, J. M. (1995). LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 8, 127–134. doi: 10.1093/protein/8.2.127

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Dommert, F., and Holm, C. (2010). Optimizing working parameters of the smooth particle mesh Ewald algorithm in terms of accuracy and efficiency. J. Chem. Phys. 133:034117. doi: 10.1063/1.3446812

PubMed Abstract | CrossRef Full Text | Google Scholar

Watt, A. D., Crespi, G. A., Down, R. A., Ascher, D. B., Gunn, A., Perez, K. A., et al. (2014). Do current therapeutic anti-Aβ antibodies for Alzheimer's disease engage the target? Acta Neuropathol. 127, 803–810. doi: 10.1007/s00401-014-1290-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Willard, L., Ranjan, A., Zhang, H., Monzavi, H., Boyko, R. F., Sykes, B. D., et al. (2003). VADAR: a web server for quantitative evaluation of protein structure quality. Nucleic Acids Res. 31, 3316–3319.

PubMed Abstract | Google Scholar

Yao, X. Q., Skjærven, L., and Grant, B. J. (2016). Rapid characterization of allosteric networks with ensemble normal mode analysis. J. Phys. Chem. B 120, 8276–8288. doi: 10.1021/acs.jpcb.6b01991

PubMed Abstract | CrossRef Full Text | Google Scholar