A Computational Study of Green Tea Extracts and their Derivatives as Potential Inhibitors for Squalene Monooxygenase.

BACKGROUND: According to the World Health Organisation, cardiovascular complications have been recognized as the leading course of death between 2000 and 2019. Cardiovascular complications are caused by excess LDL cholesterol in the body or arteries that can build up to form a plaque. There are drugs currently in clinical use called statins that target HMGCoA reductase. However, these drugs result in several side effects. This work investigated using computational approaches to lower cholesterol by investigating green tea extracts as an inhibitors for squalene monooxygenase (the second-rate-controlling step in cholesterol synthesis). METHODS: Pharmacophore modeling was done to identify possible pharmacophoric sites based on the pIC50 values. The best hypothesis generated by pharmacophore modeling was further validated by atom-based 3D QSAR, where 70% of the data set was treated as the training set. Prior molecular docking ADMET studies were done to investigate the physiochemical properties of these molecules. Glide docking was performed, followed by molecular dynamics to evaluate the protein conformational changes. RESULTS: Pharmacophore results suggest that the best molecules to interact with the biological target should have at least one hydrogen acceptor (A5), two hydrogen donors (D9 and D10), and two benzene rings (R14 and R15) for green tea polyphenols and theasinensin A. ADMET result shows that all molecules in this class have low oral adsorption. Molecular docking results showed that some green tea polyphenols have good binding affinities, with most of these structures having a docking score of less than -10 kcal/mol. Molecular dynamics further illustrated that the best-docked ligands perfectly stay within the active site over a 100 ns simulation. CONCLUSION: The results obtained from this study suggest that green tea polyphenols have the potential for inhibition of squalene monooxygenase, except for theasinensin A.

Inhibitory effects of selected cannabinoids against dipeptidyl peptidase IV, an enzyme linked to type 2 diabetes.

Ethnopharmacological relevance: In recent times the decriminalisation of cannabis globally has increased its use as an alternative medication. Where it has been used in modern medicinal practises since the 1800s, there is limited scientific investigation to understand the biological activities of this plant. Aim of the study: Dipeptidyl peptidase IV (DPP-IV) plays a key role in regulating glucose homeostasis, and inhibition of this enzyme has been used as a therapeutic approach to treat type 2 diabetes. However, some of the synthetic inhibitors for this enzyme available on the market may cause undesirable side effects. Therefore, it is important to identify new inhibitors of DPP-IV and to understand their interaction with this enzyme. Methods: In this study, four cannabinoids (cannabidiol, cannabigerol, cannabinol and Δ9-tetrahydrocannabinol) were evaluated for their inhibitory effects against recombinant human DPP-IV and their potential inhibition mechanism was explored using both in vitro and in silico approaches. Results: All four cannabinoids resulted in a dose-dependent response with IC50 values of between 4.0 and 6.9 µg/mL. Kinetic analysis revealed a mixed mode of inhibition. CD spectra indicated that binding of cannabinoids results in structural and conformational changes in the secondary structure of the enzyme. These findings were supported by molecular docking studies which revealed best docking scores at both active and allosteric sites for all tested inhibitors. Furthermore, molecular dynamics simulations showed that cannabinoids formed a stable complex with DPP-IV protein via hydrogen bonds at an allosteric site, suggesting that cannabinoids act by either inducing conformational changes or blocking the active site of the enzyme. Conclusion: These results demonstrated that cannabinoids may modulate DPP-IV activity and thereby potentially assist in improving glycaemic regulation in type 2 diabetes.

Rationalising the retro-Diels-Alder fragmentation pattern of viscutins using electrospray interface-tandem mass spectrometry coupled to theoretical modelling.

RATIONALE: Although mass spectrometry (MS) is a powerful tool in structural elucidation of unknown flavonoids based on their unique fragmentation patterns, proposing the correct fragmentation mechanism is still a challenge from tandem mass spectrometry data only. In recent years, computational tools such as molecular networking and MS2LDA have played a major role in the identification of structurally related compounds through an in-depth survey of their fragmentation patterns. METHODS: Therefore, in this study, three viscutin molecules in Viscum combreticola Engl. crude extracts were characterised using ultra-high-performance liquid chromatography-quadrupole-time-of-flight tandem mass spectrometry and MS2LDA, a computational tool. Ion-trap mass spectrometry and density functional theoretical modelling were used as confirmatory tools to rationalise the unique fragmentation patterns observed for these molecules. RESULTS: Here, MS2LDA revealed the presence of a unique Mass2Motif in all the three viscutin molecules at m/z 137, which was confirmed to be a 1,3 A- RDA (retro-Diels-Alder) fragmentation product using liquid chromatography-ion-trap mass spectrometry and density functional theoretical modelling. Moreover, MS2LDA proved to be useful in differentiating this spectral feature that was specific to viscutin molecules in the presence of other isobaric ions at m/z 137 occurring in compounds in other molecular families. CONCLUSION: Therefore, the results of the current study showed that computational tools such as MS2LDA are essential in uncovering some gas-phase fragmentation reactions of molecules in MS and that theoretical modelling is a powerful tool in rationalising these reactions in metabolite identification.

In Silico Substrate-Binding Profiling for SARS-CoV-2 Main Protease (Mpro) Using Hexapeptide Substrates.

The SARS-CoV-2 main protease (Mpro) is essential for the life cycle of the COVID-19 virus. It cleaves the two polyproteins at 11 positions to generate mature proteins for virion formation. The cleavage site on these polyproteins is known to be Leu-Gln↓(Ser/Ala/Gly). A range of hexapeptides that follow the known sequence for recognition and cleavage was constructed using RDKit libraries and complexed with the crystal structure of Mpro (PDB ID 6XHM) through extensive molecular docking calculations. A subset of 131 of these complexes underwent 20 ns molecular dynamics simulations. The analyses of the trajectories from molecular dynamics included principal component analysis (PCA), and a method to compare PCA plots from separate trajectories was developed in terms of encoding PCA progression during the simulations. The hexapeptides formed stable complexes as expected, with reproducible molecular docking of the substrates given the extensiveness of the procedure. Only Lys-Leu-Gln*** (KLQ***) sequence complexes were studied for molecular dynamics. In this subset of complexes, the PCA analysis identified four classifications of protein motions across these sequences. KLQ*** complexes illustrated the effect of changes in substrate on the active site, with implications for understanding the substrate recognition of Mpro and informing the development of small molecule inhibitors.

Insights into the Dynamics and Binding of Two Polyprotein Substrate Cleavage Points in the Context of the SARS-CoV-2 Main and Papain-like Proteases.

It is well known that vital enzymes in the replication process of the coronavirus are the SARS-CoV-2 PLpro and SARS-CoV-2 3CLpro, both of which are important targets in the search for anti-coronavirus agents. These two enzymes are responsible for cleavage at various polyprotein sites in the SARS-CoV-2 lifecycle. Herein, the dynamics of the polyprotein cleavage sequences for the boundary between non-structural proteins Nsp1 and Nsp2 (CS1) and between Nsp2 and Nsp3 (CS2) in complex with both the papain-like protein PLpro and the main protease 3CLpro were explored using computational methods. The post dynamics analysis reveals that CS1 and CS2 both have greater stability when complexed with PLpro. Of these two, greater stability is observed for the CS1-PLpro complex, while destabilization resulting in loss of CS2 from the PLpro active site is observed for CS2-PLpro, suggesting the rate of exchange by the papain-like protease is faster for CS2 compared to CS1. On the other hand, the 3CLpro main protease also reveals stability for CS1 suggesting that the main protease could also play a potential role in the cleavage at point CS1. However, destabilization occurs early in the simulation for the complex CLpro-CS2 suggesting a poor interaction and non-plausible protease cleavage of the polyprotein at CS2 by the main protease. These findings could be used as a guide in the development and design of potent COVID-19 antiviral inhibitors that mimic the CS1 cleavage site.

A multiscale ONIOM study of the buckminsterfullerene (C60) Diels-Alder reaction: from model design to reaction path analysis.

The hybrid ONIOM (Our own N-layered Integrated molecular Orbital and molecular Mechanics) formalism is employed to investigate the Diels-Alder reaction of the buckminsterfullerene C60. Our computations suggest that the ONIOM2(M06-2X/6-31G(d): SVWN/STO-3G) model, enclosing both the diene and the pyracyclene fragment of C60 in the higher-layer, provides a reasonable trade-off between accuracy and computational cost as it comes to predicting reaction energetics. Moreover, the frontier molecular orbital (FMO) theory and activation strain model (ASM) are jointly relied on to rationalize the effect of -OH and -CN substituents on the activation barrier of this reaction. Finally, reaction paths are scrutinized to get insight into the various forces underpinning the process of cycloadduct formation.

Solvent promoted tautomerism in thione-containing tetraazatricyclics: evidence from 1H NMR spectroscopy and transition state studies.

Tautomerism in the nitro substituted thione-containing traazatricyclics has been investigated. Evidence from 1H NMR indicating the existence of the tautomers has been augmented with computational studies providing evidence of the stability or otherwise of these tautomers. The role of water and DMSO in the formation of the tautomers has been explained. The role of the nitro group in assisting in the formation of the tautomers has been discussed.

AMADAR: a python-based package for large scale prediction of Diels-Alder transition state geometries and IRC path analysis.

Predicting transition state geometries is one of the most challenging tasks in computational chemistry, which often requires expert-based knowledge and permanent human intervention. This short communication reports technical details and preliminary results of a python-based tool (AMADAR) designed to generate any Diels-Alder (DA) transition state geometry (TS) and analyze determined IRC paths in a (quasi-)automated fashion, given the product SMILES. Two modules of the package are devoted to performing, from IRC paths, reaction force analyses (RFA) and atomic (fragment) decompositions of the reaction force F and reaction force constant [Formula: see text]. The performance of the protocol has been assessed using a dataset of 2000 DA cycloadducts retrieved from the ZINC database. The sequential location of the corresponding TSs was achieved with a success rate of 95%. RFA plots confirmed the reaction force constant [Formula: see text] to be a good indicator of the (non)synchronicity of the associated DA reactions. Moreover, the atomic decomposition of [Formula: see text] allows for the rationalization of the (a)synchronicity of each DA reaction in terms of contributions stemming from pairs of interacting atoms. The source code of the AMADAR tool is available on GitHub [ CMCDD/AMADAR(github.com) ] and can be used directly with minor customizations, mostly regarding the local working environment of the user.

The Molecular Basis of the Effect of Temperature on the Structure and Function of SARS-CoV-2 Spike Protein.

The remarkable rise of the current COVID-19 pandemic to every part of the globe has raised key concerns for the current public healthcare system. The spike (S) protein of SARS-CoV-2 shows an important part in the cell membrane fusion and receptor recognition. It is a key target for vaccine production. Several researchers studied the nature of this protein under various environmental conditions. In this work, we applied molecular modeling and extensive molecular dynamics simulation approaches at 0°C (273.15 K), 20°C (293.15 K), 40°C (313.15 K), and 60°C (333.15 K) to study the detailed conformational alterations in the SARS-CoV-2 S protein. Our aim is to understand the influence of temperatures on the structure, function, and dynamics of the S protein of SARS-CoV-2. The structural deviations, and atomic and residual fluctuations were least at low (0°C) and high (60°C) temperature. Even the internal residues of the SARS-CoV-2 S protein are not accessible to solvent at high temperature. Furthermore, there was no unfolding of SARS-CoV-2 spike S reported at higher temperature. The most stable conformations of the SARS-CoV-2 S protein were reported at 20°C, but the free energy minimum region of the SARS-CoV-2 S protein was sharper at 40°C than other temperatures. Our findings revealed that higher temperatures have little or no influence on the stability and folding of the SARS-CoV-2 S protein.

New Insights into the (A)Synchronicity of Diels-Alder Reactions: A Theoretical Study Based on the Reaction Force Analysis and Atomic Resolution of Energy Derivatives.

In the present manuscript, we report new insights into the concept of (a)synchronicity in Diels-Alder (DA) reactions in the framework of the reaction force analysis in conjunction with natural population calculations and the atomic resolution of energy derivatives along the intrinsic reaction coordinate (IRC) path. Our findings suggest that the DA reaction transitions from a preferentially concerted mechanism to a stepwise one in a 0.10 Å window of synchronicity indices ranging from 0.90 to 1.00 Å. We have also shown that the relative position of the global minimum of the reaction force constant with respect to the TS is an alternative and quantifiable indicator of the (a)synchronicity in DA reactions. Moreover, the atomic resolution of energy derivatives reveals that the mechanism of the DA reaction involves two inner elementary processes associated with the formation of each of the two C-C bonds. This resolution goes on to indicate that, in asynchronous reactions, the driving and retarding components of the reaction force are mostly due to the fast and slow-forming C-C bonds (elementary processes) respectively, while in synchronous reactions, both elementary processes retard and drive the process concomitantly and equivalently.

Drug Resistance in the HIV-1 Subtype C Protease Enzyme: A High Throughput Virtual Screening Approach in Search of New Ligands with Activity.

BACKGROUND: HIV-1 subtype C protease is a strategic target for antiretroviral treatment. However, resistance to protease inhibitors appears after months of treatment. Chromones and 2- biscoumarin derivatives show potential for inhibition of the HIV- subtype C protease. OBJECTIVE: Different heterocyclic structures from the ZINC database were docked against Human Immunodeficiency Virus-1 (HIV) subtype C protease crystal structure 2R5Q and 2R5P. The 5 best molecules were selected to be docked against 62 homology models based on HIV-protease sequences from infants failing antiretroviral protease treatment. This experimentation was performed with two molecular docking programs: Autodock and Autodock Vina. These molecules were modified by substituting protons with different moieties, and the derivatives were docked against the same targets. Ligand-protein interactions, physical/chemical proprieties of the molecules, and dynamics simulations were analyzed. METHODS: Docking of all of the molecules was performed to find out the binding sites of HIV-1 subtype C proteases. An in-house script was made to substitute protons of molecules with different moieties. According to the Lipinski rule of five, physical and chemical properties were determined. Complexes of certain ligands-protease were compared to the protein alone in molecular dynamics simulations. RESULTS: From the first docking results, the 5 best (lowest energy) ligands (dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[a]acridine) were selected. The binding energy of the modified ligands increased, including the poorest-performing molecules. A correlation between nature, the position, and the resulting binding energy was observed. According to the Lipinski rules, the physico-chemical characteristics of the five best-modified ligands are ideal for oral bioavailability. Molecular dynamics simulations show that some lead-protease complexes were stable. CONCLUSION: Dibenz[a,h]acridine, dibenz[a, i]acridine, NSC114903, dibenz[c,h]acridine, benzo[ a]acridine and their derivatives might be considered as promising HIV-1 subtype C protease inhibitors. This could be confirmed through synthesis and subsequent in vitro assays.

Force Field Parameters for Fe2+4S2-4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies.

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe2+4S2-4 clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe2+4S2-4 clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe2+ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe2+ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe2+4S2-4 cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

Interaction of silver nanoparticles with catechol O-methyltransferase: Spectroscopic and simulation analyses.

Catechol O-methyltransferase, an enzyme involved in the metabolism of catechol containing compounds, catalyzes the transfer of a methyl group between S-adenosylmethionine and the hydroxyl groups of the catechol. Furthermore it is considered a potential drug target for Parkinson's disease as it metabolizes the drug levodopa. Consequently inhibitors of the enzyme would increase levels of levodopa. In this study, absorption, fluorescence and infrared spectroscopy as well as computational simulation studies investigated human soluble catechol O-methyltransferase interaction with silver nanoparticles. The nanoparticles form a corona with the enzyme and quenches the fluorescence of Trp143. This amino acid maintains the correct structural orientation for the catechol ring during catalysis through a static mechanism supported by a non-fluorescent fluorophore-nanoparticle complex. The enzyme has one binding site for AgNPs in a thermodynamically spontaneous binding driven by electrostatic interactions as confirmed by negative ΔG and ΔH and positive ΔS values. Fourier transform infrared spectroscopy within the amide I region of the enzyme indicated that the interaction causes relaxation of its ß-structures, while simulation studies indicated the involvement of six polar amino acids. These findings suggest AgNPs influence the catalytic activity of catechol O-methyltransferase, and therefore have potential in controlling the activity of the enzyme.

SANCDB: an update on South African natural compounds and their readily available analogs.

BACKGROUND: South African Natural Compounds Database (SANCDB; https://sancdb.rubi.ru.ac.za/ ) is the sole and a fully referenced database of natural chemical compounds of South African biodiversity. It is freely available, and since its inception in 2015, the database has become an important resource to several studies. Its content has been: used as training data for machine learning models; incorporated to larger databases; and utilized in drug discovery studies for hit identifications. DESCRIPTION: Here, we report the updated version of SANCDB. The new version includes 412 additional compounds that have been reported since 2015, giving a total of 1012 compounds in the database. Further, although natural products (NPs) are an important source of unique scaffolds, they have a major drawback due to their complex structure resulting in low synthetic feasibility in the laboratory. With this in mind, SANCDB is, now, updated to provide direct links to commercially available analogs from two major chemical databases namely Mcule and MolPort. To our knowledge, this feature is not available in other NP databases. Additionally, for easier access to information by users, the database and website interface were updated. The compounds are now downloadable in many different chemical formats. CONCLUSIONS: The drug discovery process relies heavily on NPs due to their unique chemical organization. This has inspired the establishment of numerous NP chemical databases. With the emergence of newer chemoinformatic technologies, existing chemical databases require constant updates to facilitate information accessibility and integration by users. Besides increasing the NPs compound content, the updated SANCDB allows users to access the individual compounds (if available) or their analogs from commercial databases seamlessly.

Potential repurposing of four FDA approved compounds with antiplasmodial activity identified through proteome scale computational drug discovery and in vitro assay.

Malaria elimination can benefit from time and cost-efficient approaches for antimalarials such as drug repurposing. In this work, 796 DrugBank compounds were screened against 36 Plasmodium falciparum targets using QuickVina-W. Hits were selected after rescoring using GRaph Interaction Matching (GRIM) and ligand efficiency metrics: surface efficiency index (SEI), binding efficiency index (BEI) and lipophilic efficiency (LipE). They were further evaluated in Molecular dynamics (MD). Twenty-five protein-ligand complexes were finally retained from the 28,656 (36 × 796) dockings. Hit GRIM scores (0.58 to 0.78) showed their molecular interaction similarity to co-crystallized ligands. Minimum LipE (3), SEI (23) and BEI (7) were in at least acceptable thresholds for hits. Binding energies ranged from -6 to -11 kcal/mol. Ligands showed stability in MD simulation with good hydrogen bonding and favorable protein-ligand interactions energy (the poorest being -140.12 kcal/mol). In vitro testing showed 4 active compounds with two having IC50 values in the single-digit µM range.

Synthesis and anti-parasitic activity of N-benzylated phosphoramidate Mg2+-chelating ligands.

A series of N-benzylated phosphoramidate esters, containing a 3,4-dihydroxyphenyl Mg2+-chelating group, has been synthesised in five steps as analogues of fosmidomycin, a Plasmodium falciparum 1-deoxy-1-d-xylulose-5-phosphate reductoisomerase (PfDXR) inhibitor. The 3,4-dihydroxyphenyl group effectively replaces the Mg2+-chelating hydroxamic acid group in fosmidomycin. The compounds showed very encouraging anti-parasitic activity with IC50 values of 5.6-16.4 µM against Plasmodium falciparum parasites and IC50 values of 5.2 - 10.2 µM against Trypanosoma brucei brucei (T.b.brucei). Data obtained from in silico docking of the ligands in the PfDXR receptor cavity (3AU9)5 support their potential as PfDXR inhibitors.

Synthesis and anti-parasitic activity of achiral N-benzylated phosphoramidic acid derivatives.

Synthetic pathways have been developed to access a series of N-benzylated phosphoramidic acid derivatives as novel, achiral analogues of the established Plasmodium falciparum 1-deoxy-d-xylulose-5-phosphate reductase (PfDXR) enzyme inhibitor, FR900098. Bioassays of the targeted compounds and their synthetic precursors have revealed minimal antimalarial activity but encouraging anti-trypanosomal activity - in one case with an IC50 value of 5.4 µM against Trypanosoma brucei, the parasite responsible for Nagana (African cattle sleeping sickness). The results of relevant in silico modelling and docking studies undertaken in the design and evaluation of these compounds are discussed.

Synthesis, Characterization and Biological Activity of Some Dithiourea Derivatives.

Novel dithiourea derivatives have been designed as HIV-1 protease inhibitors using Autodock 4.2, synthesized and characterized by spectroscopic methods and microanalysis. 1-(3-Bromobenzoyl)-3-[2-([(3-bromophenyl)formami-do]methanethioylamino)phenyl]thiourea (10) and 3-benzoyl-1[(phenylformamido)methanethioyl]aminothiourea (12) gave a percentage viability of 17.9 ± 5.6% and 11.2 ± 0.9% against Trypanosoma brucei. Single crystal X-ray dif-fraction analysis of 1-benzoyl-3-(5-methyl-2-[(phenylformamido)methanethioyl]aminophenyl)thiourea (1), 3-ben-zoyl-1-(2-[(phenylformamido)methanethioyl]aminoethyl)thiourea (11), 3-benzoyl-1-[(phenylformamido)methan-ethioyl]aminothiourea (12) and 3-benzoyl-1-(4-[(phenylformamido)methanethioyl]aminobutyl)thiourea (14) have been presented. 1-(3-Bromobenzoyl)-3-[2-([(3-bromophenyl)formamido]methanethioylamino)phenyl]thiourea (10) gave a percentage inhibition of 97.03 ± 0.37% against HIV-1 protease enzyme at a concentration of 100 ?M.

Establishing Computational Approaches Towards Identifying Malarial Allosteric Modulators: A Case Study of Plasmodium falciparum Hsp70s.

Combating malaria is almost a never-ending battle, as Plasmodium parasites develop resistance to the drugs used against them, as observed recently in artemisinin-based combination therapies. The main concern now is if the resistant parasite strains spread from Southeast Asia to Africa, the continent hosting most malaria cases. To prevent catastrophic results, we need to find non-conventional approaches. Allosteric drug targeting sites and modulators might be a new hope for malarial treatments. Heat shock proteins (HSPs) are potential malarial drug targets and have complex allosteric control mechanisms. Yet, studies on designing allosteric modulators against them are limited. Here, we identified allosteric modulators (SANC190 and SANC651) against P. falciparum Hsp70-1 and Hsp70-x, affecting the conformational dynamics of the proteins, delicately balanced by the endogenous ligands. Previously, we established a pipeline to identify allosteric sites and modulators. This study also further investigated alternative approaches to speed up the process by comparing all atom molecular dynamics simulations and dynamic residue network analysis with the coarse-grained (CG) versions of the calculations. Betweenness centrality (BC) profiles for PfHsp70-1 and PfHsp70-x derived from CG simulations not only revealed similar trends but also pointed to the same functional regions and specific residues corresponding to BC profile peaks.

Identification and evaluation of bioactive natural products as potential inhibitors of human microtubule affinity-regulating kinase 4 (MARK4).

Microtubule affinity-regulating kinase 4 (MARK4) has recently been identified as a potential drug target for several complex diseases including cancer, diabetes and neurodegenerative disorders. Inhibition of MARK4 activity is an appealing therapeutic option to treat such diseases. Here, we have performed structure-based virtual high-throughput screening of 100,000 naturally occurring compounds from ZINC database against MARK4 to find its potential inhibitors. The resulted hits were selected, based on the binding affinities, docking scores and selectivity. Further, binding energy calculation, Lipinski filtration and ADMET prediction were carried out to find safe and better hits against MARK4. Best 10 compounds bearing high specificity and binding efficiency were selected, and their binding pattern to MARK4 was analyzed in detail. Finally, 100 ns molecular dynamics simulation was performed to evaluate; the dynamics stability of MARK4-compound complex. In conclusion, these selected natural compounds from ZINC database might be potential leads against MARK4, and can further be exploited in drug design and development for associated diseases.