Molecular modeling, docking and ADMET studies applied to the design of a novel hybrid for treatment of Alzheimer's disease

Metabolites profiling and cheminformatics bioprospection of selected medicinal plants against the main protease and RNA-dependent RNA polymerase of SARS-CoV-2

Despite the existence of some vaccines, SARS-CoV-2 (S-2) infections persist for various reasons relating to vaccine reluctance, rapid mutation rate, and an absence of specific treatments targeted to the infection. Due to their availability, low cost and low toxicity, research into potentially repurposing phytometabolites as therapeutic alternatives has gained attention. Therefore, this study explored the antiviral potential of metabolites of some medicinal plants [Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum (Sesame plant)] identified using liquid chromatography-mass spectrometry (LCMS) as possible inhibitory agents against the S-2 main protease (S-2 MP) and RNA-dependent RNA polymerase (RP) using computational approaches. Molecular docking was used to identify the compounds with the best affinities for the selected therapeutics targets. Afterwards, compounds with poor physicochemical characteristics, pharmacokinetics, and drug-likeness were screened out. The top-ranked compounds were further subjected to a 120-ns molecular dynamics (MD) simulation. Only quercetin 3-O-rhamnoside (-48.77 kcal/mol) had higher binding free energy than the reference standard (zafirlukast) (-44.99 kcal/mol) against S-2 MP. Conversely, all the top-ranked compounds (ellagic acid hexoside, spiraeoside, apigenin-4'-glucoside and chrysoeriol 7-glucuronide) except gnetin L (-24.24 kcal/mol) had higher binding free energy (-55.19 kcal/mol, -52.75 kcal/mol, -47.22 kcal/mol and -43.35 kcal/mol) respectively, against S-2 RP relative to the reference standard (-34.79 kcal/mol). The MD simulations study further revealed that the investigated inhibitors are thermodynamically stable and form structurally compatible complexes that impede the regular operation of the respective S-2 therapeutic targets. Although, these S-2 therapeutic candidates are promising, further in vitro and in vivo evaluation is required and highly recommended.Communicated by Ramaswamy H. Sarma.

Repurposing drugs and identification of inhibitors of integral proteins (spike protein and main protease) of SARS-CoV-2

The outbreak of Coronavirus infection (COVID-19) has prompted the World Health Organisation (WHO) to declare the outbreak, a Public Health Emergency of International concern. As part of the efforts to discover lead compounds for clinical use, 53 molecules were screened using molecular docking and dynamic simulations (MDS) techniques to identify potential inhibitors of SARS-CoV-2 spike protein (COVID-19 Sgp) and main protease (COVID-19 Mpro) or both. Lopinavir (LPV), nelfinavir (NEF), hydroxychloroquine (HCQ), remdesivir (RDV) and an irreversible inhibitor of SARS-CoV (N3) were used as standard drugs for COVID-19 Mpro, while zafirlukast (ZFK) and cefoperazone (CSP)) as standard drugs for COVID-19 Sgp. After 100 ns of MDS, with reference to standard drugs (N3, -52.463 Kcal/mol, NEF, -51.618 Kcal/mol, RDV, -48.780 Kcal/mol, LPV, -46.788 Kcal/mol, DRV, -33.655 Kcal/mol and HCQ, -21.065 Kcal/mol), five molecules, HCR, GRN, C3G, EGCG, and K7G were predicted to be promising inhibitors of COVID-19 Mpro with binding energies of -53.877 kcal/mol, -50.653 Kcal/mol, -48.600 kcal/mol, -47.798 kcal/mol and -46.902 kcal/mol, respectively. These lead molecules were then docked at receptor-binding domain (RBD) of COVID-19 Sgp to examine their inhibitory effects. C3G, GRN and K7G exhibited higher binding energies of -42.310 kcal/mol, -32.210 kcal/mol, -26.922 kcal/mol than the recorded values for the reference drugs (CSP, -35.509 kcal/mol, ZFK, -24.242 kcal/mol), respectively. The results of the binding energy and structural analyses from this study revealed that C3G, GRN and K7G could serve as potential dual inhibitors of COVID-19 Sgp and COVID-19 Mpro, while HCR and EGCG would be inhibitors of COVID-19 Mpro.Communicated by Ramaswamy H. Sarma.

Exploring the binding mode of donepezil with calf thymus DNA using spectroscopic and molecular docking methods

Donepezil (DNP) is one of approved drugs to treat Alzheimer's disease (AD). However, the potential effect of DNP on DNA is still unclear. Therefore, the interaction of DNP with calf thymus DNA (DNA) was studied in vitro using spectroscopic and molecular docking methods. Steady-state and transient fluorescence experiments showed that there was a clear binding interaction between DNP and DNA, resulting from DNP fluorescence being quenched using DNA. DNP and DNA have one binding site between them, and the binding constant (K_b ) was 0.78 × 10⁴ L·mol^-1 at 298 K. In this binding process, hydrophobic force was the main interaction force, because enthalpy change (ΔH) and entropy change (ΔS) of DNP-DNA were 67.92 kJ·mol^-1 and 302.96 J·mol^-1 ·K^-1 , respectively. DNP bound to DNA in a groove-binding mode, which was verified using a competition displacement study and other typical spectroscopic methods. Fourier transform infrared (FTIR) spectrum results showed that DNP interacted with guanine (G) and cytosine (C) bases of DNA. The molecular docking results further supported the results of spectroscopic experiments, and suggested that both Pi-Sigma force and Pi-Alkyl force were the major hydrophobic force functioning between DNP and DNA.

Target-Based Design of Promysalin Analogues Identifies a New Putative Binding Cleft in Succinate Dehydrogenase

Promysalin is a small-molecule natural product that specifically inhibits growth of the Gram-negative pathogen Pseudomonas aeruginosa (PA). This activity holds promise in the treatment of multidrug resistant infections found in immunocompromised patients with chronic illnesses, such as cystic fibrosis. In 2015, our lab completed the first total synthesis; subsequent analogue design and SAR investigation enabled identification of succinate dehydrogenase (Sdh) as the biological target in PA. Herein, we report the target-guided design of new promysalin analogues with varying alkyl chains, one of which is on par with our most potent analogue to date. Computational docking revealed that some analogues have a different orientation in the Sdh binding pocket, placing the terminal carbon proximal to a tryptophan residue. This inspired the design of an extended side chain analogue bearing a terminal phenyl moiety, providing a basis for the design of future analogues.

Facile synthesis of 1,5-disubstituted tetrazoles by reacting a ruthenium acetylide complex with trimethylsilyl azide

Treatment of [Ru]-C[triple bond, length as m-dash]CPh (1, [Ru] = (η⁵-C₅H₅)(dppe)Ru, dppe = Ph₂PCH₂CH₂PPh₂) with trimethylsilyl azide afforded the cationic nitrile complex {[Ru]NCCH₂Ph}[N₃] (2) and the further cycloaddition of 2 with trimethylsilyl azide at 60 °C afforded the N(2)-bound tetrazolato complex [Ru]N₄CCH₂Ph (3). The regiospecific alkylation of 3 gave a series of cationic N(2)-bound N(4)-alkylated-5-benzyl tetrazolato complexes {[Ru]N₄(CH₂R)CCH₂Ph}[Br] (4a, R = C₆F₅; 4b, R = Ph; 4c, R = 4-CN-C₆H₄; 4d, R = 2,6-F₂-C₆H₃; 4e, R = 6-CH₂Br-C₅NH₃) and the subsequent cleavage of the Ru-N bond of 4a-4e gave N(1)-alkylated-5-benzyl tetrazoles N₄(CH₂R)CCH₂Ph (5a-5e) in good to excellent yields and [Ru]-Br, which, on reacting with phenylacetylene, resulted in the formation of 1 thus forming a reaction cycle. The structures of 2, 3, 4a, 4c and 5a were confirmed by single-crystal X-ray diffraction analysis.

Novel tetrahydroacridine derivatives with iodobenzoic moieties induce G0/G1 cell cycle arrest and apoptosis in A549 non-small lung cancer and HT-29 colorectal cancer cells

A series of nine tetrahydroacridine derivatives with iodobenzoic moiety were synthesized and evaluated for their cytotoxic activity against cancer cell lines—A549 (human lung adenocarcinoma), HT-29 (human colorectal adenocarcinoma) and somatic cell line—EA.hy926 (human umbilical vein cell line). All compounds displayed high cytotoxicity activity against A549 (IC₅₀ 59.12–14.87 µM) and HT-29 (IC₅₀ 17.32–5.90 µM) cell lines, higher than control agents—etoposide and 5-fluorouracil. Structure–activity relationship showed that the position of iodine in the substituent in the para position and longer linker most strongly enhanced the cytotoxic effect. Among derivatives, 1i turned out to be the most cytotoxic and displayed IC₅₀ values of 14.87 µM against A549 and 5.90 µM against HT-29 cell lines. In hyaluronidase inhibition assay, all compounds presented anti-inflammatory activity, however, slightly lower than reference compound. ADMET prediction showed that almost all compounds had good pharmacokinetic profiles. 1b, 1c and 1f compounds turned out to act against chemoresistance in cisplatin-resistant 253J B-V cells. Compounds intercalated into DNA and inhibited cell cycle in G0/G1 phase—the strongest inhibition was observed for 1i in A549 and 1c in HT-29. Among compounds, the highest apoptotic effect in both cell lines was observed after treatment with 1i. Compounds caused DNA damage and H2AX phosphorylation, which was detected in A549 and HT-29 cells. All research confirmed anticancer properties of novel tetrahydroacridine derivatives and explained a few pathways of their mechanism of cytotoxic action.

Methodologies Related to Computational Models in View of Developing Anti-Alzheimer Drugs: An Overview

BACKGROUND

Since the last two decades, there has been more focus on the development strategies related to Anti-Alzheimer's drug research. This may be attributed to the fact that most of the Alzheimer's cases are still mostly unknown except for a few cases, where genetic differences have been identified. With the progress of the disease, the symptoms involve intellectual deterioration, memory impairment, abnormal personality and behavioural patterns, confusion, aggression, mood swings, irritability Current therapies available for this disease give only symptomatic relief and do not focus on manipulations of biololecular processes.

METHODS

Nearly all the therapies to treat Alzheimer's disease, target to change the amyloid cascade which is considered to be important in AD pathogenesis. New drug regimens are not able to keep pace with the ever-increasing understanding about dementia at the molecular level. Looking into these aggravated problems, we thought to put forth molecular modeling as a drug discovery approach for developing novel drugs to treat Alzheimer disease. The disease is incurable and it gets worst as it advances and finally causes death. Due to this, the design of drugs to treat this disease has become an utmost priority for research. One of the most important emerging technologies applied for this has been Computer-assisted drug design (CADD). It is a research tool that employs large-scale computing strategies in an attempt to develop a model receptor site which can be used for designing of an anti-Alzheimer drug.

RESULTS

Various models of amyloid-based calcium channels have been computationally optimized. Docking and De novo evolution are used to design the compounds. They are further subjected to absorption, distribution, metabolism, excretion and toxicity (ADMET) studies to finally bring about active compounds that are able to cross BBB. Many novel compounds have been designed which might be promising ones for the treatment of AD.

CONCLUSION

The present review describes the research carried out on various heterocyclic scaffolds that can serve as lead compounds to design Anti-Alzheimer's drugs in the future. The molecular modeling methods can thus become a better alternative for the discovery of newer Anti- Alzheimer agents. This methodology is extremely useful to design drugs in minimum time with enhanced activity keeping balanced ethical considerations. Thus, the researchers are opting for this improved process over the conventional methods hoping to achieve a sure shot way out for the sufferings of people affected by Alzheimer besides other diseases.

From dual binding site acetylcholinesterase inhibitors to allosteric modulators: A new avenue for disease-modifying drugs in Alzheimer's disease

The lack of an effective treatment for Alzheimer' disease (AD), an increasing prevalence and severe neurodegenerative pathology boost medicinal chemists to look for new drugs. Currently, only acethylcholinesterase (AChE) inhibitors and glutamate antagonist have been approved to the palliative treatment of AD. Although they have a short-term symptomatic benefits, their clinical use have revealed important non-cholinergic functions for AChE such its chaperone role in beta-amyloid toxicity. We propose here the design, synthesis and evaluation of non-toxic dual binding site AChEIs by hybridization of indanone and quinoline heterocyclic scaffolds. Unexpectely, we have found a potent allosteric modulator of AChE able to target cholinergic and non-cholinergic functions by fixing a specific AChE conformation, confirmed by STD-NMR and molecular modeling studies. Furthermore the promising biological data obtained on human neuroblastoma SH-SY5Y cell assays for the new allosteric hybrid 14, led us to propose it as a valuable pharmacological tool for the study of non-cholinergic functions of AChE, and as a new important lead for novel disease modifying agents against AD.

QSAR Models towards Cholinesterase Inhibitors for the Treatment of Alzheimer's Disease

Alzheimer's Disease (AD) is a multifactorial neurological syndrome with the combination of aging, genetic, and environmental factors triggering the pathological decline. Interestingly, the importance of the Acetylcholinesterase (AChE) enzyme has increased due to its involvement in the ß-amyloid peptide fibril formation during AD pathogenesis. In silico technique, QSAR has proven its usefulness in pharmaceutical research for the design/optimization of new chemical entities. Further, QSAR method advanced the scope of rational drug design and the search for the mechanism of drug action. It is a well-established fact that the chemical and pharmaceutical effects of a compound are closely related to its physico-chemical properties, which can be calculated by various methods from the compound structure. This chapter focuses on different Quantitative Structure-Activity Relationship (QSAR) studies carried out for a variety of cholinesterase inhibitors for the treatment of AD. These predictive models will be potentially used for further designing better and safer drugs against AD.

A new guanylhydrazone derivative as a potential acetylcholinesterase inhibitor for Alzheimer's disease: synthesis, molecular docking, biological evaluation and kinetic studies by nuclear magnetic resonance

Molecular docking, in silico studies and NMR show that the new guanylhydrazone is a promising compound for the treatment of Alzheimer's disease.

Active Compounds Against Anopheles minimus Carboxypeptidase B for Malaria Transmission-Blocking Strategy

Malaria transmission-blocking compounds have been studied to block the transmission of malaria parasites, especially the drug-resistant Plasmodium. Carboxypeptidase B (CPB) in the midgut of Anopheline mosquitoes has been demonstrated to be essential for the sexual development of Plasmodium in the mosquito. Thus, the CPB is a potential target for blocking compounds. The aim of this research was to screen compounds from the National Cancer Institute (NCI) diversity dataset and U.S. Food and Drug Administration (FDA)-approved drugs that could reduce the Anopheles CPB activity. The cDNA fragment of cpb gene from An. minimus (cpbAmi) was amplified and sequenced. The three-dimensional structure of CPB was predicted from the deduced amino acid sequence. The virtual screening of the compounds from NCI diversity set IV and FDA-approved drugs was performed against CPBAmi. The inhibition activity against CPBAmi of the top-scoring molecules was characterized in vitro. Three compounds-NSC-1014, NSC-332670, and aminopterin with IC50 at 0.99 mM, 1.55 mM, and 0.062 mM, respectively-were found to significantly reduce the CPBAmi activity.

Synthesis and DPPH scavenging assay of reserpine analogues, computational studies and in silico docking studies in AChE and BChE responsible for Alzheimer's disease

Alzheimer's disease (AD) is a fast growing neurodegenerative disorder of the central nervous system and anti-oxidants can be used to help suppress the oxidative stress caused by the free radicals that are responsible for AD. A series of selected synthetic indole derivatives were biologically evaluated to identify potent new antioxidants. Most of the evaluated compounds showed significant to modest antioxidant properties (IC50 value 399.07 140.0±50 µM). Density Functional Theory (DFT) studies were carried out on the compounds and their corresponding free radicals. Differences in the energy of the parent compounds and their corresponding free radicals provided a good justification for the trend found in their IC50 values. In silico, docking of compounds into the proteins acetylcholinesterase (AChE) and butyrylcholinesterase (BChE), which are well known for contributing in AD disease, was also performed to predict anti-AD potential.

Conformational analysis and parallel QM/MM X-ray refinement of protein bound anti-Alzheimer drug donepezil

The recognition and association of donepezil with acetylcholinesterase (AChE) has been extensively studied in the past several decades because of the former's use as a palliative treatment for mild Alzheimer disease. Herein we examine the conformational properties of donepezil and we re-examine the donepezil-AChE crystal structure using combined quantum mechanical/molecular mechanical (QM/MM) X-ray refinement tools. Donepezil's conformational energy surface was explored using the M06 suite of density functionals and with the MP2/complete basis set (CBS) method using the aug-cc-pVXZ (X = D and T) basis sets. The donepezil-AChE complex (PDB 1EVE) was also re-refined through a parallel QM/MM X-ray refinement approach based on an in-house ab initio code QUICK, which uses the message passing interface (MPI) in a distributed SCF algorithm to accelerate the calculation via parallelization. In the QM/MM re-refined donepezil structure, coordinate errors that previously existed in the PDB deposited geometry were improved leading to an improvement of the modeling of the interaction between donepezil and the aromatic side chains located in the AChE active site gorge. As a result of the re-refinement there was a 93% reduction in the donepezil conformational strain energy versus the original PDB structure. The results of the present effort offer further detailed structural and biochemical inhibitor-AChE information for the continued development of more effective and palliative treatments of Alzheimer disease.

βTrCP interacts with the ubiquitin-dependent endocytosis motif of the GH receptor in an unconventional manner

GH (growth hormone) binding to the GHR (GH receptor) triggers essential signalling pathways that promote growth and metabolic regulation. The sensitivity of the cells to GH is mainly controlled by the endocytosis of the receptor via βTrCP (β-transducin repeat-containing protein). In the present study, we show that βTrCP interacts directly via its WD40 domain with the UbE (ubiquitin-dependent endocytosis) motif in GHR, promoting GHR ubiquitination in vitro. NMR experiments demonstrated that the UbE motif is essentially unstructured, and, together with functional mapping of the UbE and βTrCP WD40 residues necessary for binding, led to a unique interaction model of βTrCP with GHR-UbE. This interaction is different from the conventional βTrCP-substrate interactions described to date. This interaction therefore represents a promising specific target to develop drugs that inhibit GHR endocytosis and increase GH sensitivity in cachexia patients.

Molecular interactions of cholinesterases inhibitors using in silico methods: current status and future prospects

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by a low amount of acetylcholine (ACh) in hippocampus and cortex. Acetylcholinesterase (AChE) is one of the most important enzymes in many living organisms including human being and other vertebrates, insects like mosquitoes, among others. Several reports have been published where it has been clearly shown that the genesis of amyloid protein plaques associated with AD is connected to modifications of both AChE and butyrylcholinesterase (BChE), since the plaque is significantly decreased in AD patients using cholinesterase inhibitors (ChEIs). This review gives some examples of these inhibitors discovered during past couple of years that have shown very prominent interactions at the active site triad of the proteins as well as different other parts of the active site like, peripheral anionic site (PAS), oxyanionic hole, anionic subsite or acyl binding pocket (ABP). Most of the inhibition and their interactions have been visualized by X-ray crystallography, but some of the other inhibitors have been studied either by molecular docking or molecular dynamic (MD) simulations or by both the in silico methods. Some of these prominent studies have been crucially observed and reported here.

Acetylcholinesterase-inhibiting activity of pyrrole derivatives from a novel marine gliding bacterium, Rapidithrix thailandica

Acetylcholinesterase-inhibiting activity of marinoquinoline A (1), a new pyrroloquinoline from a novel species of a marine gliding bacterium Rapidithrix thailandica, was assessed (IC₅₀ 4.9 μM). Two related pyrrole derivatives, 3-(2′-aminophenyl)-pyrrole (3) and 2,2-dimethyl-pyrrolo-1,2-dihydroquinoline (4), were also isolated from two other strains of R. thailandica. The isolation of 3 from a natural source is reported here for the first time. Compound 4 was proposed to be an isolation artifact derived from 3. The two isolated compounds were virtually inactive in the acetylcholinesterase-inhibitory assay (enzyme inhibition < 30% at 0.1 g L⁻¹).

Tacrine derivatives-acetylcholinesterase interaction: 1H NMR relaxation study

Two acetylcholinesterase (AChE) inhibitors structurally related to Tacrine, 6-methoxytacrine (1a) and 9-heptylamino-6-methoxytacrine (1b), and their interaction with Electrophorus Electricus AChE were investigated. The complete assignment of the 1H and 13C NMR spectra of 1a and 1b was performed by mono-dimensional and homo- and hetero-correlated two-dimensional NMR experiments. This study was undertaken to elucidate the interaction modes between AChE and 1a and 1b in solution, using NMR. The interaction between the two inhibitors and AChE was studied by the analysis of the motional parameters non-selective and selective spin-lattice relaxation times, thereby allowing the motional state of 1a and 1b, both free and bound with AChE, to be defined. The relaxation data pointed out the ligands molecular moiety most involved in the binding with AChE. The relevant ligand/enzyme interaction constants were also evaluated for both compounds and resulted to be 859 and 5412M(-1) for 1a and1b, respectively.