|Year : 2020 | Volume
| Issue : 4 | Page : 157-162
The coronavirus disease 2019 main protease inhibitor from Andrographis paniculata (Burm.f) Ness
Sukardiman1, Martha Ervina2, Mohammad Rizki Fadhil Pratama3, Hadi Poerwono4, Siswandono Siswodihardjo4
1 Department of Pharmacognosy and Phytochemistry, Faculty of Pharmacy, Universitas Airlangga, Surabaya, Indonesia
2 Doctoral Program of Pharmaceutical Sciences, Faculty of Pharmacy, Universitas Airlangga; Department of Pharmaceutical Biology, Faculty of Pharmacy, Widya Mandala Catholic University, Surabaya, Indonesia
3 Doctoral Program of Pharmaceutical Sciences, Faculty of Pharmacy, Universitas Airlangga, Surabaya; Department of Pharmacy, Faculty of Health Sciences, Universitas Muhammadiyah Palangkaraya, Palangkaraya, Indonesia
4 Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Universitas Airlangga, Surabaya, Indonesia
|Date of Submission||18-Jun-2020|
|Date of Decision||23-Jul-2020|
|Date of Acceptance||10-Aug-2020|
|Date of Web Publication||10-Oct-2020|
Prof. DR. Sukardiman
Faculty of Pharmacy, Universitas Airlangga, Campus C Unair Jl. DR. Ir. H. Soekarno Mulyorejo, Surabaya 60115
Source of Support: None, Conflict of Interest: None
The coronavirus disease 2019 (COVID-19) pandemic has attracted worldwide attention. Andrographis paniculata (Burm. f) Ness (AP) is naturally used to treat various diseases, including infectious diseases. Its Andrographolide has been clinically observed for anti-HIV and has also in silico tested for COVID-19 main protease inhibitors. Meanwhile, the AP phytochemicals content also provides insight into the molecular structures diversity for the bioactive discovery. This study aims to find COVID-19 main protease inhibitor from AP by the molecular docking method and determine the toxicity profile of the compounds. The results obtained two compounds consisting of flavonoid glycosides 5,4'-dihydroxy-7-O-β -D-pyran-glycuronate butyl ester and andrographolide glycoside 3-O-β-D-glucopyranosyl-andrographolide have lower free binding energy and highest similarity in types of interaction with amino acid residues compared to its co-crystal ligands (6LU7) and Indinavir or Remdesivir. The toxicity prediction of the compounds also reveals their safety. These results confirm the probability of using AP phytochemical compounds as COVID-19 main protease inhibitors, although further research must be carried out.
Keywords: 6 LU7, Andrographis paniculata, coronavirus disease 2019, main protease inhibitor, molecular docking
|How to cite this article:|
Sukardiman, Ervina M, Fadhil Pratama MR, Poerwono H, Siswodihardjo S. The coronavirus disease 2019 main protease inhibitor from Andrographis paniculata (Burm.f) Ness. J Adv Pharm Technol Res 2020;11:157-62
|How to cite this URL:|
Sukardiman, Ervina M, Fadhil Pratama MR, Poerwono H, Siswodihardjo S. The coronavirus disease 2019 main protease inhibitor from Andrographis paniculata (Burm.f) Ness. J Adv Pharm Technol Res [serial online] 2020 [cited 2022 May 16];11:157-62. Available from: https://www.japtr.org/text.asp?2020/11/4/157/297701
| Introduction|| |
The coronavirus (CoV) has attracted attention, especially since its virulence features changed to human infection. Overall, six CoV varieties four of them induce common colds symptoms, while two varieties of SARS and middle-east respiratory syndrome produce more severe to deadly infections. The SARS-CoV-2 caused severe acute respiratory diseases (coronavirus disease 2019 [COVID-19]); which is responsible for the global pandemic and caused 5.67% deaths in confirmed cases. To date, the exponential rates in confirmed cases were observed, in addition, the occurrence of second wave erase.
The development of COVID-19 therapies, show promising results, but is still far from convincing. Vaccine development began, and some claimed have good results. Furthermore, repurposing drugs such as Remdesivir (REM) or Chloroquine is considered to have a good influence on the COVID-19 patients' recovery. Moreover, the researchers also looked at the medicinal plant metabolites. Natural products provide extraordinary in structure and complexity for drug discovery by providing a drift for synthetic or semi-synthetic antiviral compounds. As summarized some phytochemicals worked on the virulence proteins key. In particular, in silico study of several promising Chinese and Indian medicinal plants for COVID-19 treatment were reported., These studies showed that plant metabolites are promising in the discovery of COVID-19 compounds.
Among receptor targets, previous research has developed therapies for COVID-19, the 3C-like protease/3CLpro (SARS-CoV-2 Mpro). is an interesting and quite developed one for the COVID-19 therapeutic target. Not only for drug repurposing but it is also a potential point for tracing inhibitors derived from phytochemicals. Screening of medicinal plants library showed potential phytochemical compounds such as Isoflavone and Myricitrin compared to antiviral nelfinavir.
Among the various phytochemicals that have been studied, Andrographis paniculata (AP) is one very interesting to be further investigated is. It is commonly used in several countries for common cold treatment and many more diseases,, This Acanthaceae family plant has been reported as SARS-CoV-2 Mpro inhibitors. A separate study, explained that AP suppressed the increase in NOD-like receptor protein 3, caspase-1, and interleukin-1β, which are extensively involved in the pathogenesis of SARS-CoV and possibly SARS-CoV-2 as well. Enmozhi et al. even specifically examined the Andrographolide, as SARS-CoV-2 Mpro inhibitor with molecular docking method. However, the study only tested andrographolide without considering other contents, and the docking protocol is also not fully listed, making it difficult to reproduce the results.
Moreover, the AP consists of diterpene lactone and flavonoids as well as their aglycone or glycoside derivatives form. Some were including andrograpanin, neoandrographolide, 14-deoxyandrographolide, 3-O-β-D-glucopyranosyl-14, 19-dideoxyandrographolide, 5,4-dihydroxy-7-O-β-D-pyran-glycuronate-butyl-ester, apigenin, 7-O-methylwogonin, and onysilin., These compounds have become an object for synthetic or semisynthetic modification for various diseases that have experienced resistance, for example, as anti-HIV drugs.
From the background and previous research, this study aims to reveal the potency of AP phytoconstituents as SARS-CoV-2 Mpro inhibitors. The method used is molecular docking and in silico toxicity to compare safety profiles of AP phytoconstituents. The results of this study are expected not only to contribute to the development of new drugs or potential derivatives from AP phytoconstituents but also to provide a scientific basis for the traditional use of AP in COVID-19 therapy.
| Materials and Methods|| |
Preparation of ligands and receptor
The 45 AP's secondary metabolites ligands are listed in [Table 1] base on reference information.,,, The ligands were prepared with hardware, and software as describe Pratama and Gusdinar. The programs were included HyperChem 7.5., OpenBabel 2.4.1 and AutoDockTools 1.5.6. The molecular information in the SMILES format is used in the prediction of toxicity properties.
|Table 1: Energy binding of Andrographis paniculata's and reference ligands|
Click here to view
The structure of 3-chymotrypsin-like protease (3CL-protease, protein data bank [PDB] ID: 6LU7) COVID-19 retrieved from PDB (https://www.rcsb.org/) in.pdb format.
Validation of the docking process
The validation prior docking process was done by extracted the co-crystal ligands from its receptors, and re-docking as well as their point of binding on its protein active sites. The compound co-crystal ligand was prepared. The RMSD redocking obtained with PyMOL 2.3.1.
Autodock Vina 1.1.2 and Chimera 1.13.1 programs were used for the docking process, while Discovery Studio Visualizer 19.1.0 was used for interaction obtaining. The fixed parameters were the size and position of a grid box as resulted in the validation process.,
The toxicity prediction of the most potent substances
The toxicity study of the three most potential substances was processed. The results obtained then compared and analyzed further.
| Results|| |
The receptor 6LU7 contains a monomer protein (chain A)-N3 (CL) with a resolution of 2.16 Š, while the 6LU7-AP's complex ligands structure resulted in 1.9 Š, r-value free and work of 0.262 and 0.229, respectively. These values were good parameter for molecular docking experimental studies, which serves the overlapping/accuracy of the redocking atom's positions to its original ligand.
[Table 1] is listed the energy binding (ΔG) of 45 AP's compounds and 7 reference ligands docking with the 6 LU7. The ΔG presents relative binding ligand-receptor affinities. The AP's ligands of 5,4-dihydroxy-7-O-β-D-pyran-glycuronate-butyl-ester (DGE) and 3-O-β-D-glucopyranosyl-andrographolide (GAD) have the lowest ΔG compare to Indinavir and REM. These ligands and 7,8-dimethoxy-2-hydroxy-5-O-β-D-glucopyranosyloxyflavone (DGF) have a better affinity compare to its original N3. Furthermore, the Hydroxychloroquine, Chloroquine, and Favipiravir were on the bottom list compare to AP's and others. The three of AP's most potent and two of reference ligands were proceeded for further examination. [Figure 1]a shows the redocking co-crystal ligand position that almost overlaps with the crystallographic ligand position, while [Figure 1]b, [Figure 1]c, [Figure 1]d, [Figure 1]e, [Figure 1]f presents their visualization, respectively. [Figure 2] presents the three most active AP's chemical structures.
|Figure 1: Two-dimension interactions of (a) co-crystal ligand (N3). (b) 5,4'-Dihydroxy-7-O-β-D-pyran-glycuronate-butyl-ester. (c) 3-O-β-D-glucopyranosyl-andrographolide. (d) 7,8-Dimethoxy-2'-hydroxy-5-O-β-D-glucopyranosyloxyflavone. (e) Indinavir. (f) Remdesivir with binding site of 6LU7 receptor (24-Thr) Threonine amino acid number 24, (Ala) Alanine, (Arg) Arginine, (Asn) Asparagine, (Asp) Aspartic acid, (Cys) Cysteine, (Gln) Glutamine, (Glu) glutamic acid, (Gly) Glycine, (His) Histidine, (Leu) Leucine, (Met) Methionine, (Phe) Phenylalanine, (Pro) Proline, (Thr) Threonine, (Tyr) Tyrosine, (Ser) Serine|
Click here to view
|Figure 2: The three most AP's 6LU7 Mpro. (a) 5,4-Dihydroxy-7-O-β-D-pyran-glycuronate butyl ester. (b) 3-O-β-D-glucopyranosyl-andrographolide. (c) 7,8-Dimethoxy-2'-hydroxy-5-O-β-D-glucopyranosyloxyflavone|
Click here to view
The toxicity parameter of DGE, GAD, DGF, two reference compounds INDI and REM were obtained with ProTox II describe in [Figure 3].
|Figure 3: Toxicity prediction probability of AP's ligands and reference, 5,4'-Dihydroxy-7Oβ-D-pyran-glycuronate-butyl-ester, 3Oβ-D-glucopyranosyl-andrographolide, 7,8-Dimethoxy-2'-hydroxy-5-O-β-D-glucopyranosyloxyflavone, Indinavir, Remdesivir|
Click here to view
| Discussion|| |
The molecular docking of 45 AP's ligands and 7 reference synthetic drugs, which were consist of 4 anti-HIV Food and Drug Administration approved and 3 most suggested drugs were presented in this study to reveal their potency as COVID-19 main protease inhibitor. The use of reference synthetic drugs in this study was as internal validation process intended, so the results can be compared to others finding. The binding site to the 6LU7 Mpro receptor was performed using the X-ray crystal structure of CL as a co-crystal ligand, as can be seen in [Figure 1]a. This shows that the docking process is validated by producing ligand binding that approaches crystallographic ligands positions. The 6LU7 docking results in the RdRp binding pocket suggests that REM as the anti-HIV drug isin vitro highly effective in fighting COVID-19 infection. The REM was used for Ebola and Marburg virus outbreaks (proposed for COVID-19 clinical trials by the Gilead company).
The ΔG of the reference ligand is −8.10 ± 0.10 kcal/mol, with the RMSD value is 1.981 Š. Based on the data presented in [Table 1], the two AP ligands (DGE, and GAD) have lower ΔG compare to anti-HIV Indinavir and REM. All of them and DGF, have lower ΔG than its original CL. Interestingly observed similarities among the active ligands, that all of them were glycoside form, while DGE and DGF were share flavone structure with the difference in glycosyl side chain [Figure 2]. The AP's flavone glycoside was reported as HL-60 antiproliferative, while the Andrographolide was known as a broad-spectrum antiviral. Anti COVID-19 of Andrographolide was presented in silico study,, while this research counted its ΔG was −6.80 kcal/mol, which was on least order among the AP's compounds.
Research obtained ΔG phytochemicals to 6LU7 compare to Lopinavir; which were in range of-11.82 to 13.51 compare to −11.62 kcal/mol. While in this research, the Lopinavir (−7.27 ± 0.06 kcal/mol) was listed on the 26th rank among AP's. It showed as good reference, by yielding in range ratio to DGE and Hesperidin (1.151 to 1.620) 1.151 to 1.62. The difference of ΔG values was supposed of varieties on the analysis parameters, which is common in the molecular docking study. Although, these results highlighted the potency of AP's ligands to bind the 6LU7, so as their potency to inhibit COVID-19 replication.
Another finding of this study also obtained the GAD (Andrographolide glycoside form) possessed lower ΔG than its aglycone. The physicochemical of the glycoside influenced the solubility of the ligands. Kren wrote the glycoside influence the polarity of the compounds. This result provides insight of glycoside may induce active site of the 6LU7 receptors. Although Andrographolide has an advantage as main phytoconstituents and easier to absorb due to its small molecule compare to its glycoside form. Furthermore, this research finding was the three most suggested drugs of hydrochloroquine, chloroquine, and favipiravir have the least ΔG compare to most AP's ligands and anti-HIV drugs.
There was 25 total interaction between the N3-6LU7 receptor, with 36% consisted of hydrogen and Van der Walls bonds [Figure 1]a. [Figure 1]b, presents the interaction of DGE which has 19 bonding. DGE has the highest similarity (76% of amino acid and 36% of bonding type) to the N3 interaction among all ligands.
[Figure 1]c describes 17% of hydrogen bonds of GAD. This glycoside ligand represented the least interaction compare to others. [Figure 1]d illustrates DGF interactions. Although it only possessed 2 hydrogen bonds on 143-Gly, 163-His (18%), it has a similarity higher (72% and 32%) than GAD.
The reference antiviral ligands of INDI and REM are defined in [Figure 1]e and [Figure 1]f, respectively. These ligands have 20–21 of bonding interaction out of 25 compares to N3, and have 80% similarities of amino acids residues. REM has a similar type of interaction higher than INDI (48% to 20%). These docking studies suggest that DGE and REM have the most overlapping binding pose with N3.
LD50 is an important quantitative parameters for toxicity. All ligands have LD50 not <5000 mg/kg/day, except for GAD (590 mg/kg) and REM (1000 mg/kg). Therefore, the toxicity class was as may be harmful and slightly toxic if orally taken. The further potential toxicity was including the immunotoxicity, carcinogen-and mutagenicity. The probability of the potential toxicity is relatively low, except for INDI hepatotoxicity. Furthermore, in the 15model receptors toxicity resulted, most ligands have potential toxicity to PG/H1, except for DGE and REM. Interestingly the GAD was active to many receptors, including AhR, AND-LBD, ARO, and MMP, while DGF was active only to AhR. The GAD and INDI have the toxicity of AOA, while DGF was on AdR-β2 [Figure 3]. These results suggested AP's glycoside flavone of DGE was relatively safer than others. Moreover, these results determine the side effects which may arise. These indicate the association between these docking results and the effectiveness of treating COVID-19 still needs further examination.
| Conclusions|| |
The A. paniculata (Burm. f) Ness. have potential as 6LU7 COVID-19 main protease inhibitors. The DGE, GAD, DGF; possess low energy and overlap binding interactions on docking sites to the 6LU7 pockets. The DGE was the safest AP's compound compare to others. This finding needs further and extended examination of its use for COVID-19 treatment.
To Faculty of Pharmacy, Universitas Airlangga.
Financial support and sponsorship
RnD Project (PPTI Risbang RistekDikti 2019).
Conflicts of interest
There are no conflicts of interest.
| References|| |
Liu DX, Liang JQ, Fung TS. Human Coronavirus-229E, -OC43, -NL63, and -HKU1. Reference Module in Life Sciences 2020:B978-0-12-809633-8.21501-X.
Zhang N, Li C, Hu Y, Li K, Liang J, Wang L, et al
. Current development of COVID-19 diagnostics, vaccines and therapeutics. Microbes Infect 2020;22:231-5.
Wang M, Cao R, Zhang L, Yang X, Liu J, Xu M, et al
. Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro
. Cell Res 2020;30:269-71.
Zhang DH, Wu KL, Zhang X, Deng SQ, Peng B. In silico
screening of Chinese herbal medicines with the potential to directly inhibit 2019 novel coronavirus. J Integr Med 2020;18:152-8.
Vellingiri B, Jayaramayya K, Iyer M, Narayanasamy A, Govindasamy V, Giridharan B, et al
. COVID-19: A promising cure for the global panic. Sci Total Environ 2020;725:138277.
Jin Z, Du X, Xu Y, Deng Y, Liu M, Zhao Y, et al
. Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors. Nature 2020;582:289-93.
Ul Qamar MT, Alqahtani SM, Alamri MA, Chen LL. Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants. J Pharm Anal 2020;10:313-19.
Jayakumar T, Hsieh CY, Lee JJ, Sheu JR. Experimental and clinical pharmacology of Andrographis paniculata
and its major bioactive phytoconstituent andrographolide. Evid Based Complement Alternat Med 2013;2013:846740.
Sukardiman H, Widyawaruyanti A, Sismindari, Zaini NC. Apoptosis inducing effect of andrographolide on TD-47 human breast cancer cell line. Afr J Tradit Complement Altern Med 2007;4:345-51.
Hossain MS, Urbi Z, Sule A, Hafizur Rahman KM. Andrographis paniculata
(Burm. f.) Wall. ex Nees: A review of ethnobotany, phytochemistry, and pharmacology. Sci W J 2014;274905.
Liu Z, Xiao X, Wei X, Li J, Yang J, Tan H, et al
. Composition and divergence of coronavirus spike proteins and host ACE2 receptors predict potential intermediate hosts of SARS-CoV-2. J Med Virol 2020;92:595-601.
Enmozhi SK, Raja K, Sebastine I, Joseph J. Andrographolide as a potential inhibitor of SARS-CoV-2 main protease: An in silico
approach. J Biomol Struc Dynamics 2020;1-7.
Chen LX, He H, Xia GY, Zhou KL, Qiu F. A new flavonoid from the aerial parts of Andrographis paniculata
. Nat Prod Res 2014;28:138-43.
Uttekar MM, Das T, Pawar RS, Bhandari B, Menon V, et al
. Anti-HIV activity of semisynthetic derivatives of andrographolide and computational study of HIV-1 gp120 protein binding. Eur J Med Chem 2012;56:368-74.
U.S National Library Medicine, National Center for Biotechnology Information. PubChem Substance and Compound Databases; 2020. Available from: https://pubchem.ncbi.nlm.nih.gov/
. [Last accessed on 2020 Mar 24].
Pratama MR, Gusdinar T. Docking study of secondary metabolites from Glycyrrhiza glabra
as PPAR-γ agonist. Biointerface Res Appl Chem 2019;9:4006-10.
Forli S, Huey R, Pique ME, Sanner MF, Goodsell DS, Olson AJ. Computational protein-ligand docking and virtual drug screening with the AutoDock suite. Nat Protoc 2016;11:905-19.
Jin Z, Du X, Xu Y, Deng Y, Liu M, Zhao Y, et al.
Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors Nature 2020; 582:289-93.
Calligari P, Bobone S, Ricci G, Bocedi A. Molecular investigation of SARS-CoV-2 proteins and their interactions with antiviral drugs. Viruses 2020;12:1-13.
Křen V. Glycoside vs. Aglycon: The role of glycosidic residue in biological activity. In: Fraser-Reid BO, Tatsuta K, Thiem J, editors. Glycoscience. Springer-Verlag Berlin, Heidelberg; 2008.
Banerjee P, Eckert AO, Schrey AK, Preissner R. ProTox-II: A webserver for the prediction of toxicity of chemicals. Nucleic Acids Res 2018;46:W257-W263.
Wahl J, Smieško M. Endocrine disruption at the androgen receptor: Employing molecular dynamics and docking for improved virtual screening and toxicity prediction. Int J Mol Sci 2018;19:1784.
[Figure 1], [Figure 2], [Figure 3]
|This article has been cited by|
||Efficacy of add-on Ayurveda and Yoga intervention in health care workers of tertiary care hospital during COVID-19: Randomized controlled trial
| ||Vandna Sharma, Bharat Krushna Khuntia, Manish Soneja, Vitthal G. Huddar, S. Ramakrishnan, Payal Sharma, Shubhangi Rathore, Varun V, Mohit Wadhawan, Varun Chhabra, Aman Agrawal, Man Singh, Arvind Kumar, Tanuja Manoj Nesari, Gautam Sharma |
| ||Complementary Therapies in Clinical Practice. 2022; : 101601 |
|[Pubmed] | [DOI]|
||Molecular docking unveils the potential of andrographolide derivatives against COVID-19: an in silico approach
| ||Ravichandran Veerasamy, Rohini Karunakaran |
| ||Journal of Genetic Engineering and Biotechnology. 2022; 20(1) |
|[Pubmed] | [DOI]|
||The Effect of Plant Metabolites on Coronaviruses: A Comprehensive
Review Focusing on their IC50 Values and Molecular Docking Scores
| ||Parastou Farshi, Eda Ceren Kaya, Fataneh Hashempour-Baltork, Kianoush Khosravi-Drani |
| ||Mini-Reviews in Medicinal Chemistry. 2022; 22(3): 457 |
|[Pubmed] | [DOI]|
||Pharmacogenetics and Precision Medicine Approaches for the Improvement of COVID-19 Therapies
| ||Mohitosh Biswas, Nares Sawajan, Thanyada Rungrotmongkol, Kamonpan Sanachai, Maliheh Ershadian, Chonlaphat Sukasem |
| ||Frontiers in Pharmacology. 2022; 13 |
|[Pubmed] | [DOI]|
||Andrographis paniculata (Burm. f.) Wall. ex Nees: An Updated Review of Phytochemistry, Antimicrobial Pharmacology, and Clinical Safety and Efficacy
| ||Sanower Hossain,Zannat Urbi,Hidayah Karuniawati,Ramisa Binti Mohiuddin,Ahmed Moh Qrimida,Akrm Mohamed Masaud Allzrag,Long Chiau Ming,Ester Pagano,Raffaele Capasso |
| ||Life. 2021; 11(4): 348 |
|[Pubmed] | [DOI]|
||Andrographolide: A Herbal-Chemosynthetic Approach for Enhancing Immunity, Combating Viral Infections, and Its Implication on Human Health
| ||Archana Mishra, Haq Abdul Shaik, Rakesh Kumar Sinha, Bakht Ramin Shah |
| ||Molecules. 2021; 26(22): 7036 |
|[Pubmed] | [DOI]|
||Artemether, Artesunate, Arteannuin B, Echinatin, Licochalcone B and Andrographolide Effectively Inhibit SARS-CoV-2 and Related Viruses In Vitro
| ||Yunjia Hu,Meiqin Liu,Hongbo Qin,Haofeng Lin,Xiaoping An,Zhengli Shi,Lihua Song,Xinglou Yang,Huahao Fan,Yigang Tong |
| ||Frontiers in Cellular and Infection Microbiology. 2021; 11 |
|[Pubmed] | [DOI]|
||Andrographis paniculata (Burm. F.) Wall. Ex Nees, Andrographolide, and Andrographolide Analogues as SARS-CoV-2 Antivirals? A Rapid Review
| ||Xin Yi Lim, Janice Sue Wen Chan, Terence Yew Chin Tan, Bee Ping Teh, Mohd Ridzuan Mohd Abd Razak, Saharuddin Mohamad, Ami Fazlin Syed Mohamed |
| ||Natural Product Communications. 2021; 16(5): 1934578X21 |
|[Pubmed] | [DOI]|
||Activity Prediction of Bioactive Compounds Contained in Etlingera elatior Against the SARS-CoV-2 Main Protease: An In Silico Approach
| ||Dwi Syah Fitra Ramadhan,Taufik Muhammad Fakih,Arfan Arfan |
| ||Borneo Journal of Pharmacy. 2020; 3(4): 235 |
|[Pubmed] | [DOI]|