key: cord-0924195-g8gkr6jw authors: Goyzueta-Mamani, Luis Daniel; Barazorda-Ccahuana, Haruna Luz; Chávez-Fumagalli, Miguel Angel; F. Alvarez, Karla Lucia; Aguilar-Pineda, Jorge Alberto; Vera-Lopez, Karin Jannet; Lino Cardenas, Christian Lacks title: In Silico Analysis of Metabolites from Peruvian Native Plants as Potential Therapeutics against Alzheimer’s Disease date: 2022-01-28 journal: Molecules DOI: 10.3390/molecules27030918 sha: b046ab67bdc25a76451ec50e104028a39e23e186 doc_id: 924195 cord_uid: g8gkr6jw Background: Despite research on the molecular bases of Alzheimer’s disease (AD), effective therapies against its progression are still needed. Recent studies have shown direct links between AD progression and neurovascular dysfunction, highlighting it as a potential target for new therapeutics development. In this work, we screened and evaluated the inhibitory effect of natural compounds from native Peruvian plants against tau protein, amyloid beta, and angiotensin II type 1 receptor (AT1R) pathologic AD markers. Methods: We applied in silico analysis, such as virtual screening, molecular docking, molecular dynamics simulation (MD), and MM/GBSA estimation, to identify metabolites from Peruvian plants with inhibitory properties, and compared them to nicotinamide, telmisartan, and grapeseed extract drugs in clinical trials. Results: Our results demonstrated the increased bioactivity of three plants’ metabolites against tau protein, amyloid beta, and AT1R. The MD simulations indicated the stability of the AT1R:floribundic acid, amyloid beta:rutin, and tau:brassicasterol systems. A polypharmaceutical potential was observed for rutin due to its high affinity to AT1R, amyloid beta, and tau. The metabolite floribundic acid showed bioactivity against the AT1R and tau, and the metabolite brassicasterol showed bioactivity against the amyloid beta and tau. Conclusions: This study has identified molecules from native Peruvian plants that have the potential to bind three pathologic markers of AD. Alzheimer's disease (AD) is a chronic neurodegenerative disorder characterized by the progressive manifestations of disturbances on rational thinking, memory loss, cognitive decline, and mood changes [1] [2] [3] . These clinical manifestations are consequences of the formation of neurofibrillary tangles, senile plaques, glial cell activation, and cerebrovascular dysregulation [4] . The importance of understanding the mechanisms of cerebrovascular alterations and their relation to AD has gained more attention, since cerebrovascular dysfunction can cause the degenerative processes of smooth muscle cells, astrocytes, pericytes, and endothelial cells [5] [6] [7] [8] . Furthermore, several studies have demonstrated the reduced resting cerebral blood flow, low vasoreactivity, and neurovascular coupling dysregulation in AD patients [9] [10] [11] , and that more than 50% of AD-diagnosed patients also have a cerebrovascular lesion. Besides the tau protein and amyloid beta, the constant search for novel molecular targets has implicated the renin-angiotensin system (RAS), which is involved in cere-native plants: Smallanthus sonchifolius (Yacon) of the family Asteraceae, showing that the antihypertensive and anti-inflammatory properties of the leaf extract on LPS-stimulated mouse microglial cells in vitro [42] prevented the deposition of amyloid plaques and neurotoxicity because of its antioxidant effect on the hippocampus [43] . Lepidium meyenii (Maca) is a hypocotyl of the family Brassicaceae with neuroprotective effects due to the inhibition of acetylcholinesterase (AChE) activity, improving scopolamine-induced memory deficits [38, 44] . Secondary metabolites, such as macamides have shown a reduction of Mninduced mitochondrial toxicity in glioblastoma U/87 MG cells, regulating inflammation, and glucose homeostasis, hence beneficial for AD treatment [45, 46] . Croton lechleri (Sangre de Drago), of the family Euphorbiaceae, a sap rich in phenolic compounds with the ability to reduce ROS formation, LDL oxidations, and albumin glycation, thereby reducing the risk of vascular diseases [47, 48] . Uncaria tomentosa (Uña de Gato or Cat's claw) is a vine of the Rubiaceae family with antimutagenic and neuroprotective properties. The polyphenols, sterols, and alkaloids present in Cat's Claw have shown a positive effect in cholinergic dysfunction, inhibiting AChE activity and DNA repairment, avoiding oxidative stress and inflammation, and thus, cardiovascular and autoimmune diseases [49] [50] [51] . Physalis peruviana (Aguaymanto) is a berry of the Solanaceae family with anti-inflammatory and cytoprotective properties in astrocytic cells (T98G) exposed to neurotoxic stimuli, preserving mitochondrial functions and nuclei damage [52, 53] . Minthostachys mollis (Muña) of the Lamiaceae family is rich in flavonoids that inhibit COX1 and COX2 enzymatic production, hence avoiding the proinflammatory process. According to Benites et al., the essential oil of this species has a potential cytotoxic effect against cancer cells (T24, DU-145, and MCF-7) [54, 55] . Virtual screening was performed to screening potential molecules against tau, amyloid beta, and AT1R ( Figure 1 and Tables S7-S9) . Among these molecules, we selected the top three with the lowest binding affinity energy for our study: rutin (effective against the three receptors), brassicasterol (effective against tau and amyloid beta), and floribundic acid (effective against tau and AT1R). Figure 1 also shows the drug-likeness score calculation of the complete set of molecules against the three targets using OSIRIS DataWarrior, where more potential molecules are represented for possibly further analysis. In parallel, studied compounds from ongoing clinical and preclinical studies (phase II and III) were also analyzed and used as controls (Table 1) . We chose one compound per receptor as control: nicotinamide for tau, grapeseed extract for amyloid beta, and telmisartan for AT1R. In parallel, studied compounds from ongoing clinical and preclinical studies (phase II and III) were also analyzed and used as controls (Table 1) . We chose one compound per receptor as control: nicotinamide for tau, grapeseed extract for amyloid beta, and telmisartan for AT1R. The ADME/TOX study was performed to investigate and evaluate the therapeutic potential of metabolites derived from Peruvian plants to predict their metabolism in the human body. In Table 2 , the toxicity assessment of the chosen compounds was also analyzed to evaluate their hepatotoxicity, carcinogenicity, mutagenicity, and cytotoxicity. Most compounds showed a favorable effect in this screening, except floribundic acid and kaurenoic acid, which showed a cytotoxic and hepatotoxic effect, respectively. The ADME/TOX study was performed to investigate and evaluate the therapeutic potential of metabolites derived from Peruvian plants to predict their metabolism in the human body. In Table 2 , the toxicity assessment of the chosen compounds was also analyzed to evaluate their hepatotoxicity, carcinogenicity, mutagenicity, and cytotoxicity. Most compounds showed a favorable effect in this screening, except floribundic acid and kaurenoic acid, which showed a cytotoxic and hepatotoxic effect, respectively. The active sites or regions of the receptors against the ligands were determined using the PATCHDOCK server (http://bioinfo3d.cs.tau.ac.il/PatchDock/, accessed on 10 August 2021). The snapshot obtained from 100 ns of NPT simulation was used for the thermodynamics parameter calculations. The root mean squared deviation (RMSD) and root mean squared fluctuations (RMSF) per residue of the backbone were analyzed. The structure variation was calculated by RMSD values of protein:ligand from 0 to 100 ns, reaching stability at 70 ns ( Figure 2 ). We can see in Figure 2A The tendency of the RMSD plot for the AT1R-rutin system shows that the protein does not constantly deviate from its original conformation even though this compound generates the largest fluctuations in AT1R (RMSF analysis); (B) the RMSD plot for amyloid beta in its native state (black line) shows a high structural instability due to the size and lack of other secondary conformations. However, the coupled systems reduce this deviation; (C) tau with ligands shows a tendency of high conformational stability. When we looked at the RMSD plot of amyloid beta ( Figure 2B ), we noticed a constantly deviate from its original conformation even though this compound generates the largest fluctuations in AT1R (RMSF analysis); (B) the RMSD plot for amyloid beta in its native state (black line) shows a high structural instability due to the size and lack of other secondary conformations. However, the coupled systems reduce this deviation; (C) tau with ligands shows a tendency of high conformational stability. When we looked at the RMSD plot of amyloid beta ( Figure 2B ), we noticed a significant deviation in the RMSD trend in all systems analyzed. Indeed, the RMSD average values for the native receptor (1.33 nm), grapeseed extract (0.84 nm), rutin (0.88 nm), kaurenoic acid (1.03 nm), and brassicasterol (1.14 nm) were obtained. This was caused by the lack of secondary regions in the peptides. However, when it is docked, the fluctuations are reduced (RMSF plot from Figure 2B ). Figure 2C shows the RMSD graph concerning the initial structure of tau, where the comparison of the average RMSDs between the tau receptor is less than 0.36 nm, and an average RMSD for tau-rutin and tau-nicotinamide is equal to 0.26 nm. At the same time, the RMSF per residue of the tau indicates the reduction of fluctuations in coupled systems (RMSF plot from Figure 2C ). The molecular mechanics generalized born surface area (MM/GBSA) binding free energy estimation was carried out considering the last 500 snapshots from MD simulations of receptor-ligand systems. Likewise, the generalized born (GB) model and solvent accessibility methods are used for the MM/GBSA calculation. The results obtained of this estimation are shown in Table 3 , in which among the ligand-receptor systems, the best option for amyloid beta, tau, and AT1R was rutin (48.05 ± 7.0 kcal·mol −1 ), brassicasterol (−38.05 ± 3.5 kcal·mol −1 ), and floribundic acid (58.98 ± 4.4 kcal·mol −1 ), respectively. The second-best option, and shared among the three receptors, were rutin against amyloid beta, tau, and AT1R, in which the binding energy was found to be equal to 48.05 ± 7.0 kcal·mol −1 , −27.25 ± 6.1 kcal·mol −1 , and −39.56 ± 3.3 kcal·mol −1 , respectively. There is significant energy contribution by the Van der Waals energies in the systems regarding electrostatic and generalized born energies. Substances under clinical trials were run in order to be used as controls, acting as references point to evaluate the potential anti-Alzheimer effect of the Peruvian natural plants. Thereof, based on these calculations, the ligands rutin for amyloid beta, brassicasterol for tau, and floribundic acid for AT1R were deduced with considerable binding affinity by comparison to their control reference, which was grapeseed extract, nicotinamide, and telmisartan, respectively. For the last three decades, over 900 billion USD has been invested in drug development and clinical trials with no proven effective benefit for treating patients with AD [56] . Multiple studies have recently demonstrated a strong relationship between the pathogenesis of AD and some vascular-associated diseases, including atherosclerosis and hypertension [57, 58] . Moreover, Aguilar et al., 2021, have shown that vascular smooth muscle cells contribute to the neuroinflammation and tau hyperphosphorylation in the early and late stages of the disease [6] . These findings identify the vascular tissue as a novel therapeutic target for the development of drugs to reduce the impact of vascular dysfunction on the initiation and progression of AD. In this context, the neuroprotective effects of AT1R blockers (Losartan and Resveratrol) indicate that RAS components can be plausible targets to develop pharmacological innovations for AD [59, 60] . For these reasons, the research for new compounds with neurovascular or antineurodegenerative bioactivity is gaining considerable attention in the field [61] . In this regard, a total of 94 options between drugs and antibodies are in clinical trials (phases II and III) [62] . Among those compounds, some natural products, such as grapeseed extract, nicotinamide, nicotine, saffron, colostrinin, and others, are considered as promising therapeutic candidates [63] [64] [65] . Furthermore, in vivo experiments using animal models of dementia and AD have been used to evaluate the potential therapeutic potential of other plant-based molecules, such as curcumin, quercetin, resveratrol, piperine, and epigallocatechin gallate. Peru is known worldwide for its incredible phytogenic biodiversity, in which~10% of the world's flora grows [37] . Polypharmacology is emerging as a new paradigm to treat complex diseases by regulating multiple targets to achieve desired physiological responses [66] . Additionally, a multi-target drug can decrease the risk of drug−drug interactions and diminish the number of pharmacokinetic and safety profile tests [67] . The ethnopharmacological knowledge of the Peruvian population helped to identify and categorize more than 5000 Peruvian plants with multiple pharmacological effects. In this sense, we aimed to characterize the polypharmacological properties of metabolites from Peruvian plants known by their neuroprotective and vasoactive properties and thus identify potential therapeutic candidates for treating patients with AD. Based on our in silico and computer-aided drug screening, we found remarkable bioactivity against three pathologic markers of AD, including tau protein, amyloid beta and AT1R. Our results demonstrated the highest binding energy for the systems as follows: AT1R:floribundic acid (a compound present in S. sonchilofolius), amyloid beta:brassicasterol (a compound present in L. meyenii), and tau-rutin (a compound present in C. lechleri). Interestingly, the rutin molecule has shown strong polypharmacological activity against all three pathogenic markers, brassicasterol against tau and amyloid beta, and floribundic acid against tau and AT1R. It is worth mentioning that the metabolite rutin showed a higher affinity for tau protein, amyloid beta, and AT1R compared to the controls, drugs in clinical trials, such as nicotinamide, grapeseed extract, and telmisartan, respectively. In conclusion, we identified three metabolites with the potential to bind three pathologic markers of AD and inform the development of future therapies against AD and other neurodegenerative diseases. In addition, the beneficial effects of these metabolites have been reported by different authors. Rutin has shown an inhibitory effect on α-glucosidase, effective for the treatment of diabetes [68] , an antiviral effect on SARS-CoV-2 [69] , and anti-inflammatory [70] and anticarcinogenic properties [71] . The sterol brassicasterol has shown a dual anti-infective property against herpes simplex virus type 1 and Mycobacterium tuberculosis [72] , an antiviral effect against human adenovirus [73] , and an anti-tumoral [74] and potential application as an AD biomarker [75] . For floribundic acid, besides its antioxidant effect, there is no other research based on the effects or bioactivity of this clerodane diterpenoid [76] . The binding site of the proteins is essential during the protein-ligand interaction and crucial for drug design. Our study revealed the catalytic dyad for AT1R composed of hydrogen bonds with floribundic acid (ARG167, TYR113), glucobrassicin (TYR184), rutin (TYR113) telmisartan (LYS199), and hydrophobic bonds (TYR87, TYR175, ILE288, TRP84, THR260, and PHE182 as common interacting residues (Figure 3) ). Similar binding sites were observed with Olmesartan, where LYS 199, ARG167, TYR113, and TRP84 were noticed [77, 78] . For amyloid beta, the hydrogen bonds were observed with rutin (LYS28, ASP23, and PHE20) and hydrophobic bonds (HIS14, VAL39, ALA21, ILE32, and ILE41) as common For amyloid beta, the hydrogen bonds were observed with rutin (LYS28, ASP23, and PHE20) and hydrophobic bonds (HIS14, VAL39, ALA21, ILE32, and ILE41) as common interacting residues (Figure 4) . Our results agree with other in silico evidence of the interaction between statins and amyloid beta, where authors evaluated the protective role of these drugs against amyloidogenesis and neurodegeneration, such as atorvastatin, rosuvastatin, pravastatin, and lovastatin, where ASP23, VAL39, and ALA21 residues were identified, respectively [79] . According to Castro et al., 2020, fucosterol from the algae Sargassum horridum showed an inhibitory aggregation effect of amyloid beta fibrils and stabilized and interacted with HIS14, PHE20, ALA21, ILE32, VAL39, ILE41, and LYS28 residues [80] . The residue HIS14 was also observed in the binding site of memantine and amyloid beta [81] . interacting residues (Figure 4) . Our results agree with other in silico evidence of the interaction between statins and amyloid beta, where authors evaluated the protective role of these drugs against amyloidogenesis and neurodegeneration, such as atorvastatin, rosuvastatin, pravastatin, and lovastatin, where ASP23, VAL39, and ALA21 residues were identified, respectively [79] . According to Castro et al., 2020, fucosterol from the algae Sargassum horridum showed an inhibitory aggregation effect of amyloid beta fibrils and stabilized and interacted with HIS14, PHE20, ALA21, ILE32, VAL39, ILE41, and LYS28 residues [80] . The residue HIS14 was also observed in the binding site of memantine and amyloid beta [81] . For tau, hydrogen bonds were observed with rutin (SER352), brassicasterol (GLN351), ILE360, and THR361 as common interacting residues ( Figure 5 ). In fact, the For tau, hydrogen bonds were observed with rutin (SER352), brassicasterol (GLN351), ILE360, and THR361 as common interacting residues ( Figure 5 ). In fact, the residue GLN351 was also an active binding site for the interaction with titanium dioxide nanoparticles, where the alteration of secondary and tertiary structures of tau was aimed at to avoid further aggregation [82] . In addition, some second-generation positron emission tomography tau tracers were developed and engineered to target the active site found in our work, composed of GLN351and ILE360 [83, 84] . residue GLN351 was also an active binding site for the interaction with titanium dioxide nanoparticles, where the alteration of secondary and tertiary structures of tau was aimed at to avoid further aggregation [82] . In addition, some second-generation positron emission tomography tau tracers were developed and engineered to target the active site found in our work, composed of GLN351and ILE360 [83, 84] . The MD simulation analyzes the stability and dynamic states of the protein-ligand systems. The protein backbone RMSD of tau, amyloid beta, and AT1R, bound with the The MD simulation analyzes the stability and dynamic states of the protein-ligand systems. The protein backbone RMSD of tau, amyloid beta, and AT1R, bound with the best option of the ligands are represented in Figure 6 . We observed that tau, amyloid beta, and AT1R started to stabilize and equilibrate at 60 ns, and all backbones remained stable compared to their respective controls (proteins without ligand binding). The stability of the compound related to AT1R was: rutin > floribundic acid >> telmisartan > glucobrassicin; the compound related to amyloid beta was: grapeseed extract > rutin > kaurenoic acid > brassicasterol; and the compound related to tau was: rutin > nicotinamide > brassicasterol > floribundic acid. The slight variation may be attributed to the complex conformations. > brassicasterol > floribundic acid. The slight variation may be attributed to the complex conformations. In the same figure, the RMSF analysis showed a fluctuation in the residue regions of 200-250 for AT1R, 10-35 for amyloid beta, and 340-360 for tau. The RMSF values of the systems relative to the amino acids of the complexes were evaluated to compare their flexibilities. The more fluctuated amino acids in the active sites, related to their ligands, were for AT1R, ILE228, and GLU227 attributed to the EIQKN sequence, an intracellular loop of the AT1R and a cytoplasmatic domain essential for receptor activation and Gprotein selection and coupling [85] . For amyloid beta, the residues GLN15, VAL24, and SER26 were attributed to the α-helices. The ALA30 and GLY29 regions were less structured and formed an α-helix (unfolded structures). This α-helix stabilization may be crucial to prevent the formation of the insoluble β-sheet form of amyloid beta [86, 87] . For tau, the amino acids the ASP348 and ARG349 residues are responsible for saline bonds formation and β-helical folding, arrangements found in the AD-tau fibrils. Surprisingly, the molecules considered as controls, such as grapeseed extract, nicotinamide, and telmisartan for the receptors amyloid beta, tau protein, and AT1R, respectively, were not the best options with the highest binding energy [88, 89] . In the same figure, the RMSF analysis showed a fluctuation in the residue regions of 200-250 for AT1R, 10-35 for amyloid beta, and 340-360 for tau. The RMSF values of the systems relative to the amino acids of the complexes were evaluated to compare their flexibilities. The more fluctuated amino acids in the active sites, related to their ligands, were for AT1R, ILE228, and GLU227 attributed to the EIQKN sequence, an intracellular loop of the AT1R and a cytoplasmatic domain essential for receptor activation and G-protein selection and coupling [85] . For amyloid beta, the residues GLN15, VAL24, and SER26 were attributed to the α-helices. The ALA30 and GLY29 regions were less structured and formed an α-helix (unfolded structures). This α-helix stabilization may be crucial to prevent the formation of the insoluble β-sheet form of amyloid beta [86, 87] . For tau, the amino acids the ASP348 and ARG349 residues are responsible for saline bonds formation and β-helical folding, arrangements found in the AD-tau fibrils. Surprisingly, the molecules considered as controls, such as grapeseed extract, nicotinamide, and telmisartan for the receptors amyloid beta, tau protein, and AT1R, respectively, were not the best options with the highest binding energy [88, 89] . The literature search strategy for collecting the Peruvian natural products and ongoingclinical drugs was found on the National Center for Biotechnology Information (NCBI) databases, PubMed (accessed on 15 July 2021), and ChEMBL. Duplicated compounds were removed, and only registered, catalogued, and well-characterized compounds were considered for analysis. The Medical Subject Headings (MeSH terms) were employed in the string query to improve the search accuracy. The dataset was retrieved from PubMed on 16 July 2021, based on the following search string: "Peruvian" [MeSH Terms] AND "Natural products" [MeSH Terms] AND "Alzheimer disease". After the first step of screening of the compounds from natural sources (Physalis peruvianus [90] [91] [92] , Minthostachys mollis [54, 93] , Uncaria tomentosa [94] , Croton lechleri [76, 95, 96] , Smallanthus sonchifolius [97] [98] [99] [100] , and Lepidium meyenii [101, 102] ), the simplified molecularinput line-entry system (SMILE) was searched and retrieved from DrugBank (https://go. drugbank.com, accessed on 17 July 2021) or PubChem (https://pubchem.ncbi.nlm.nih.gov, accessed on 18 July 2021) servers. The physicochemical properties such as total molecular weight (MW), octanol/water partition coefficient (cLogP), number of H-bond acceptors (HBAs), number of H-bond donors (HBDs), and the molecular polar surface area (PSA) for each compound were calculated within the Osiris DataWarrior v05.02.01 software [103] . The FASTA sequences of the compound targets the angiotensin II type 1 receptor (AT1R) (PDB ID: 4ZUD), tau (region from the 306 to 378 residues) (PDB ID: 5O3T), and amyloid beta (PDB ID: 1IYT) were retrieved from Protein Data Bank (https://www.rcsb.org, accessed on 10 September 2021). Furthermore, the compounds were imported into Open-Babel within the Python Prescription Virtual Screening Tool (PyRx) [104] and subjected to energy minimization. PyRx performs structure-based virtual screening, applying docking simulations using the AutoDock Vina tool [105] , whereas the drug targets were uploaded as macromolecules. For the analysis, the search space encompassed the whole of the modelled 3D models; the docking simulation was then run at an exhaustiveness of 8 and set to only output the lowest energy pose. The Osiris DataWarrior software was employed to calculate the drug-likeness score of each compound; the calculation is based on a library of~5300 substructure fragments and their associated drug-likeness scores. This library was prepared by fragmenting 3300 commercial drugs and 15,000 commercial non-drug-like Fluka compounds. Furthermore, the potential tumorigenic, mutagenic, and irritant actions of each compound were predicted by comparison to a precompiled fragment library derived from the RTECS (Registry of Toxic Effects of Chemical Substances) database [103] . The Lipinski rule of five (RO5) was used to evaluate compounds that violated more than one of the rules: absorption, distribution, metabolism, excretion, toxicity properties, and predicted medial lethal dose (LD50) of the compounds were predicted using Swiss-ADME (http://www.swissadme.ch, accessed on 11 October 2021) [106] and PROTOX-II web server (https://tox-new.charite.de/protox_II/, accessed on 14 October 2021) [107] . Ligands were downloaded from PubChem server. Each input file was generated on the ACPYPE server [108] . For this calculation, the semi-empirical charge method AM1-BCC and parameters to reproduce HF/6-31G * RESP charges were selected and implemented in ANTECHAMBER [109] . AMBER to define atom type and zero total net charges were considered. The molecules retrieved are shown in Tables S1-S6, and natural compounds in clinical phase II and III studies were considered as controls, as shown in Table S7 , Supplementary Materials [62] . The protein (AT1R) and peptides (tau and amyloid beta) structures were edited in a Text Editor, removing not-essential molecules for calculation. Moreover, the best systems obtained by virtual screening were considered to analyze the coupled systems. Consequently, molecular dynamics (MD) simulations were carried out in GROMACS v. 2020 [110] . We considered an AMBER99SB-ILDN force field, explicit solvent (TIP4P water model), and ions to neutralize the system. The systems were minimized with the steepest descent algorithm for 50,000 steps. Then, the equilibrium MD calculations in the NVT (number of molecules, volume, and temperature constant) canonical ensemble with the V-rescale thermostat at 309.65 K with a calculation of 10 ns were performed. The final analysis was the MD production without restraint condition in the isobaricisothermal ensemble with V-rescale thermostat at 309.65 K and Parrinello-Rahman barostat at 1 bar of reference pressure with a trajectory of 100 ns. The thermodynamic parameters of each system were analyzed using the Gromacs tools and the plotting graph realized with Gnuplot v5 package [111] . The molecular mechanics generalized born surface area (MM/GBSA) was used to determine the binding free energy of receptor-ligand using gmx MMPBSA v1.4.1 [112] tool based on MMPBSA.py [113] , from AmberTools20 suite. In general, this estimation calculates the free energy difference between the bound and unbound form of the receptor as the following equation: where ∆G bind is the binding free energy and G complex (free energy of complex), G protein (free energy of protein), and G ligand (free energy ligand). The molecular visualizations of the interaction of the complex receptor-ligand structure and output files were carried out using the VMD v1.9.4 (Visual Molecular Dynamics) software [114] , and the 2D diagrams of receptor-ligand were generated and visualized using the Ligplot + software [115] . This study identified three natural substances with potential effects against tau (rutin), amyloid beta (brassicasterol), and AT1R receptors (floribundic acid) using in silico analysis from the scarce available data and research of Peruvian natural products. The ADME/TOX analysis predicted no toxic effects for these compounds, showing them as possible alternatives for drug formulations. The RMSD simulations and predictions have demonstrated higher docking stabilities for the systems with their optimal ligands compared to their controls. We observed the potential uses of Peruvian native plants, such as S. sonchilofolius, L. meyenii, and C. lechleri, against Alzheimer's disease or dementia-related neurovascular issues. However, further in vitro and in vivo studies of these compounds are required to confirm these results. In addition, several vegetal sources might be further studied or characterized in order to analyze their potential against more neurovascular or neurodegenerative diseases. Supplementary Materials: The following are available online. Table S1 : Natural compounds description of Smallanthus sonchifolius; Table S2 : Natural compounds description of Lepidium meyenii; Table S3 : Natural compounds description of Croton lechleri; Table S4 : Natural compounds description of Uncaria tomentosa; Table S5 : Natural compounds description of Minthostachys mollis; Table S6 : Natural compounds description of Physalis peruvianus; Alzheimer's: A Progressive Brain Disease: Causes, Symptoms, and Prevention Update on Alzheimer's disease therapy and prevention strategies Late-onset Alzheimer disease Neuropathological alterations in Alzheimer disease Inflammatory pathology markers (activated microglia and reactive astrocytes) in early and late onset Alzheimer disease: A post mortem study Vascular smooth muscle cell dysfunction contribute to neuroinflammation and Tau hyperphosphorylation in Alzheimer disease Brain microvascular pericytes in vascular cognitive impairment and dementia. Front Blood-brain barrier breakdown in Alzheimer disease and other neurodegenerative disorders Neurovascular unit dysregulation, white matter disease, and executive dysfunction: The shared triad of vascular cognitive impairment and Alzheimer disease Early role of vascular dysregulation on late-onset Alzheimer's disease based on multifactorial data-driven analysis The pathobiology of vascular dementia Angiotensin-receptor blockers (ARBs) and risk of Alzheimer's disease: A systematic review and metaanalysis Brain Renin-Angiotensin System at the Intersect of Physical and Cognitive Frailty Angiotensin receptor (AT2R) agonist C21 prevents cognitive decline after permanent stroke in aged animals-A randomized double-blind pre-clinical study The brain renin-angiotensin system: Location and physiological roles The renin angiotensin system and Alzheimer's disease Brain renin-angiotensin-a new look at an old system An Overview of the Potential Use of Ethno-Medicinal Plants Targeting the Renin-Angiotensin System in the Treatment of Hypertension Resurgence of herbal antihypertensives in management of hypertension In vitro inhibition of [3H]-angiotensin II binding on the human AT1 receptor by proanthocyanidins from Guazuma ulmifolia bark Cardioprotective effects of Qishenyiqi mediated by angiotensin II type 1 receptor blockade and enhancing angiotensin-converting enzyme 2 Pharmacological effects of Astragaloside IV: A literature review Effects of gastrodin injection on blood pressure and vasoactive substances in treatment of old patients with refractory hypertension: A randomized controlled trial Gastrodin reduces blood pressure by intervening with RAAS and PPARγ in SHRs Silibinin inhibits the production of pro-inflammatory cytokines through inhibition of NF-κB signaling pathway in HMC-1 human mast cells Chemical Composition and Hypotensive Effect of Campomanesia xanthocarpa. Evid. Based Complement Chlorogenic acids in cardiovascular disease: A review of dietary consumption, pharmacology, and pharmacokinetics et Chilensis: Sive, Descriptiones, et Icones Plantarum Peruvianarum, et Chilensium, Secundum Systema Linnaeanum Digestae, cum Characteribus Plurium Generum Evulgatorum Reformatis; Typis Gabrielis de Sancha Ethnopharmacological approaches for dementia therapy and significance of natural products and herbal drugs. Front Effect of three different cultivars of Lepidium meyenii (Maca) on learning and depression in ovariectomized mice Imidazole alkaloids from lepidium m eyenii Constituents of Lepidium meyenii 'maca' Isorhynchophylline: A plant alkaloid with therapeutic potential for cardiovascular and central nervous system diseases Uncaria rhynchophylla and its major constituents on central nervous system: A review on their pharmacological actions Vasodilatory effect of kaempferol-7-O-α-L-rhamnopyranoside via NO-cGMP-PKG signaling Bioactive compounds from plants used in Peruvian traditional medicine Aqueous and hydroalcoholic extracts of Black Maca (Lepidium meyenii) improve scopolamine-induced memory impairment in mice Multi-target approach for natural products in inflammation Computational polypharmacology comes of age An up-to-date overview of computational polypharmacology in modern drug discovery Smallanthus sonchifolius leaf attenuates neuroinflammation Yacon leaf extract supplementation demonstrates neuroprotective effect against memory deficit related to β-amyloid-induced neurotoxicity Lepidium meyenii Supplemented Diet Modulates Neurobehavioral and Biochemical Parameters in Mice Fed High-Fat High-Sugar Diet Natural compounds for the prevention and treatment of cardiovascular and neurodegenerative diseases Neuroprotective activity of macamides on manganese-induced mitochondrial disruption in U-87 MG glioblastoma cells Bioactive compounds extraction of Croton lechleri barks from Amazon forest using chemometrics tools Effects of Croton lechleri sap (Sangre de Drago) on AGEs formation, LDL oxidation and oxidative stress related to vascular diseases Comparison of the chemical constituents and anti-Alzheimer's disease effects of Uncaria rhynchophylla and Uncaria tomentosa Neuroprotective effects of aqueous extracts of Uncaria tomentosa: Insights from 6-OHDA induced cell damage and transgenic Caenorhabditis elegans model Uncaria tomentosa improves cognition, memory and learning in middle-aged rats Actividad antinflamatoria de extractos y fracciones obtenidas de cálices de Physalis peruviana L Extracts of Physalis peruviana Protect Astrocytic Cells Under Oxidative Stress With Rotenone Chemical composition, in vitro cytotoxic and antioxidant activities of the essential oil of Peruvian Minthostachys mollis Griseb. Boletín Latinoam Efecto Antiinflamatorio de un Gel Elaborado a Base de Extracto Seco de Hojas de Minthostachys Mollis (Muña) en Rattus Rattus FDA position statement "Early Alzheimer's disease: Developing drugs for treatment, Guidance for Industry Alzheimer's and Aducanumab: Unjust Profits and False Hopes Common brain structural alterations associated with cardiovascular disease risk factors and Alzheimer's dementia: Future directions and implications Effect of resveratrol on reactive oxygen species-induced cognitive impairment in rats with angiotensin II-induced early Alzheimer's disease angiotensin II type 1 receptor blocker is able to relieve age-related cognitive impairment Brain angiotensin II and angiotensin IV receptors as potential Alzheimer's disease therapeutic targets Alzheimer's disease drug development pipeline: 2020. Alzheimer's Dement Effects of acute subcutaneous nicotine on attention, information processing and short-term memory in Alzheimer's disease Six-month randomized, placebo-controlled, double-blind, pilot clinical trial of curcumin in patients with Alzheimer disease New insights into clinical trial for colostrinin TM in Alzheimer's disease A perspective on multi-target drug discovery and design for complex diseases Computational multitarget drug design LC-QTOF-MS/MS Based Molecular Networking Approach for the Isolation of α-Glucosidase Inhibitors and Virucidal Agents from Coccinia grandis Antiviral Activity of Metabolites from Peruvian Plants against SARS-CoV-2: An In Silico Approach Djulis (Chenopodium formosanum) and its bioactive compounds affect vasodilation, angiotensin converting enzyme activity, and hypertension Rutin: A Flavonoid as an Effective Sensitizer for Anticancer Therapy; Insights into Multifaceted Mechanisms and Applicability for Combination Therapy Brassicasterol with Dual Anti-Infective Properties against HSV-1 and Mycobacterium tuberculosis, and Cardiovascular Protective Effect: Nonclinical In Vitro and In Silico Assessments The novel antiviral properties of brassicasterol against human adenovirus Mechanistic analysis of bladder cancer suppression by brassicasterol in a rat model of N-butyl-N-(4-hydroxybutyl) nitrosamine-induced bladder cancer The plant sterol brassicasterol as additional CSF biomarker in Alzheimer's disease Identification of minor secondary metabolites from the latex of Croton lechleri (Muell-Arg) and evaluation of their antioxidant activity The dynamic properties of angiotensin II type 1 receptor inverse agonists in solution and in the receptor site Structural Basis for Ligand Recognition and Functional Selectivity at Angiotensin Receptor In silico evidence of direct interaction between statins and β-amyloid Fucosterol from Sargassum horridum as an amyloid-beta (Aβ1-42) aggregation inhibitor: In vitro and in silico studies In silico approach to discover the role of metals for the treatment of Alzheimer's disease amyloid-beta (Aβ) peptide Amorphous aggregation of tau in the presence of titanium dioxide nanoparticles: Biophysical, computational, and cellular studies Different positron emission tomography tau tracers bind to multiple binding sites on the tau fibril: Insight from computational modeling Computational Insight into the Binding Profile of the Second-Generation PET Tracer PI2620 with Tau Fibrils Structure of two fragments of the third cytoplasmic loop of the rat angiotensin II AT1A receptor: Implications with respect to receptor activation and G-protein selection and coupling Interaction between amyloid beta peptide and an aggregation blocker peptide mimicking islet amyloid polypeptide Free energy landscape for alpha-helix to beta-sheet interconversion in small amyloid forming peptide under nanoconfinement Analysis of the relationship between metalloprotease-9 and tau protein in Alzheimer's disease. J. Alzheimer's Dis Elucidating Tau function and dysfunction in the era of cryo-EM Screening of phytometabolites in oil from cape gooseberry (Physalis peruviana L.) fruit pomace Oil goldenberry (Physalis peruviana L.) Bioactive phytochemicals of cape gooseberry (Physalis peruviana L.). In Bioactive Compounds in Underutilized Fruits and Nuts Essential oil of Minthostachys mollis: Extraction and chemical composition of fresh and stored samples. Arab A review from patents inspired by two plant genera: Uncaria and Hamelia Polyphenolic compounds from Croton lechleri SP-303, an antiviral oligomeric proanthocyanidin from the latex of Croton lechleri (Sangre de Drago) Smallanthus sonchifolius) as a Food Supplement: Health-Promoting Benefits of Fructooligosaccharides Topical antiinflammatory activity of yacon leaf extracts Smallanthus sonchifolius Poepp. & Endl.) as a novel source of health promoting compounds: Antioxidant activity, phytochemicals and sugar content in flesh, peel, and whole tubers of seven cultivars Hypoglycaemic Property of Yacon (Smallanthus sonchifolius (Poepp. and Hendl.) H. Robinson): A Review Lepidium meyenii Walp), una revisión sobre sus propiedades biológicas Medicinal effects of Peruvian maca (Lepidium meyenii): A review DataWarrior: An open-source program for chemistry aware data visualization and analysis Small-molecule library screening by docking with PyRx AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules ProTox-II: A webserver for the prediction of toxicity of chemicals ACPYPE-Antechamber python parser interface Automatic atom type and bond type perception in molecular mechanical calculations GROMACS: Fast, flexible, and free 0: A portable interactive plotting utility gmx_MMPBSA: A new tool to perform end-state free energy calculations with GROMACS py: An efficient program for end-state free energy calculations VMD: Visual molecular dynamics Multiple ligand-protein interaction diagrams for drug discovery The authors declare no conflict of interest.