S1 Subsite Research Articles

Identification of novel and potent inhibitors of SARS-CoV-2 main protease from DNA-encoded chemical libraries.

In vitro screening of large compound libraries with automated high-throughput screening is expensive and time-consuming and requires dedicated infrastructures. Conversely, the selection of DNA-encoded chemical libraries (DECLs) can be rapidly performed with routine equipment available in most laboratories. In this study, we identified novel inhibitors of SARS-CoV-2 main protease (Mpro) through the affinity-based selection of the DELopen library (open access for academics), containing 4.2 billion compounds. The identified inhibitors were peptide-like compounds containing an N-terminal electrophilic group able to form a covalent bond with the nucleophilic Cys145 of Mpro, as confirmed by x-ray crystallography. This DECL selection campaign enabled the discovery of the unoptimized compound SLL11 (IC50 = 30 nM), proving that the rapid exploration of large chemical spaces enabled by DECL technology allows for the direct identification of potent inhibitors avoiding several rounds of iterative medicinal chemistry. As demonstrated further by x-ray crystallography, SLL11 was found to adopt a highly unique U-shaped binding conformation, which allows the N-terminal electrophilic group to loop back to the S1' subsite while the C-terminal amino acid sits in the S1 subsite. MP1, a close analog of SLL11, showed antiviral activity against SARS-CoV-2 in the low micromolar range when tested in Caco-2 and Calu-3 (EC50 = 2.3 µM) cell lines. As peptide-like compounds can suffer from low cell permeability and metabolic stability, the cyclization of the compounds will be explored in the future to improve their antiviral activity.

A Structural Investigation of the Interaction between a GC-376-Based Peptidomimetic PROTAC and Its Precursor with the Viral Main Protease of Coxsackievirus B3

The conservation of the main protease in viral genomes, combined with the absence of a homologous protease in humans, makes this enzyme family an ideal target for developing broad-spectrum antiviral drugs with minimized host toxicity. GC-376, a peptidomimetic 3CL protease inhibitor, has shown significant efficacy against coronaviruses. Recently, a GC-376-based PROTAC was developed to target and induce the proteasome-mediated degradation of the dimeric SARS-CoV-2 3CLPro protein. Extending this approach, the current study investigates the application of the GC-376 PROTAC to the 3CPro protease of enteroviruses, specifically characterizing its interaction with CVB3 3CPro through X-ray crystallography, NMR (Nuclear Magnetic Resonance) and biochemical techniques. The crystal structure of CVB3 3CPro bound to the GC-376 PROTAC precursor was obtained at 1.9 Å resolution. The crystallographic data show that there are some changes between the binding of CVB3 3CPro and SARS-CoV-2 3CLPro, but the overall similarity is strong (RMSD on C-alpha 0.3 Å). The most notable variation is the orientation of the benzyloxycarbonyl group of GC-376 with the S4 subsite of the proteases. NMR backbone assignment of CVB3 3CPro bound and unbound to the GC-376 PROTAC precursor (80% and 97%, respectively) was obtained. This information complemented the investigation, by NMR, of the interaction of CVB3 3CPro with the GC-376 PROTAC, and its precursor allows us to define that the GC-376 PROTAC binds to CVB3 3CPro in a mode very similar to that of the precursor. The NMR relaxation data indicate that a quench of dynamics of a large part of the protein backbone involving the substrate-binding site and surrounding regions occurs upon GC-376 PROTAC precursor binding. This suggests that the substrate cavity, by sampling different backbone conformations in the absence of the substrate, is able to select the suitable one necessary to covalently bind the substrate, this being the latter reaction, which is the fundamental step required to functionally activate the enzymatic reaction. The inhibition activity assay showed inhibition potency in the micromolar range for GC-376 PROTAC and its precursor. Overall, we can conclude that the GC-376 PROTAC fits well within the binding sites of both proteases, demonstrating its potential as a broad-spectrum antiviral agent.

Crystal structure of exfoliative toxin D from Staphylococcus aureus

Staphylococcus aureus exfoliative toxins (ETs) are serine proteases responsible for staphylococcal scalded skin syndrome. Four ETs, ETA, ETB, ETD, and ETE, have been identified, all of which cleave desmoglein-1. This study presents the crystal structure of ETD at 1.75 Å resolution. The protein exhibits a structure composed of two β-barrels and two α-helices as described in previous studies of ETs. A predicted model of ETD in complex with Ile380–Glu381–Gly382–Pro383 (IEGP), a segment of human desmoglein-1 (hDsg1), was constructed. Glu381 of hDsg1 was predicted to interact with as many as six amino acid residues in ETD, whereas two amino acid residues in ETD primarily constituted subsite S1′, and a space near subsite S1′ was noted. It is likely that polypeptide chains located near the IEGP segment in the predicted structure of hDsg1 bind to this space. The structure of loop D, which was predicted to participate in subsite S2′, in ETD was markedly different from those in other ETs.

Bridging the Gap in the Structure-Function Paradigm of Enzymatic PET Degradation-Aromatic Residue Driven Balanced Interactions with Catalytic and Anchoring Subsite.

Understanding all parameters contributing to enzyme activity is crucial in enzyme catalysis. For enzymatic PET degradation, this involves examining the formation of the enzyme-PET complex. In IsPETase (WT), a PET-degrading enzyme from Ideonella sakaiensis, mutating two non-catalytic residues (DM) significantly enhances activity. Such mutations, depending on their position in the tertiary structure, fine-tune enzyme function. However, detailed molecular insights into these mutations' structure-function relationship for PET degradation are lacking. This study characterizes IsPETase's catalytic ability compared to WT TfCut2 using molecular dynamics simulations and quantum mechanical methods. We explore the conformational landscape of the enzyme-PET complex and quantify residue-wise interaction energy. Notably, aromatic and hydrophobic residues Tyr, Trp, and Ile in the catalytic subsite S1, and aromatic Phe and polar Asn in the anchoring subsite S3, crucially optimize PET binding. These residues enhance PET specificity over non-aromatic plastics. Our findings suggest that the balance between binding at subsite S1 and subsite S3, which is influenced by cooperative mutations, underlies catalytic activity. This balance shows a positive correlation with experimentally obtained kcat/Km values: WT TfCut2<WT IsPETase≪DM IsPETase.

Robustaflavone as a novel scaffold for inhibitors of native and auto-proteolysed human neutrophil elastase

ABSTRACT Human neutrophil elastase (HNE) plays a key role in initiating inflammation in the cardiopulmonary and systemic contexts. Pathological auto-proteolysed two-chain (tc) HNE exhibits reduced binding affinity with inhibitors. Using AutoDock Vina v1.2.0, 66 flavonoid inhibitors, sivelestat and alvelestat were docked with single-chain (sc) HNE and tcHNE. Schrodinger PHASE v13.4.132 was used to generate a 3D-QSAR model. Molecular dynamics (MD) simulations were conducted with AMBER v18. The 3D-QSAR model for flavonoids with scHNE showed r 2 = 0.95 and q 2 = 0.91. High-activity compounds had hydrophobic A/A2 and C/C2 rings in the S1 subsite, with hydrogen bond donors at C5 and C7 positions of the A/A2 ring, and the C4’ position of the B/B1 ring. All flavonoids except robustaflavone occupied the S1’-S2’ subsites of tcHNE with decreased AutoDock binding affinities. During MD simulations, robustaflavone remained highly stable with both HNE forms. Principal Component Analysis suggested that robustaflavone binding induced structural stability in both HNE forms. Cluster analysis and free energy landscape plots showed that robustaflavone remained within the sc and tcHNE binding site throughout the 100 ns MD simulation. The robustaflavone scaffold likely inhibits both tcHNE and scHNE. It is potentially superior to sivelestat and alvelestat and can aid in developing therapeutics targeting both forms of HNE.

Novel selective proline-based peptidomimetics for human cathepsin K inhibition

Human cathepsin K (CatK) stands out as a promising target for the treatment of osteoporosis, considering its role in degrading the bone matrix. Given the small and shallow S2 subsite of CatK and considering its preference for proline or hydroxyproline, we now propose the rigidification of the leucine fragment found at the P2 position in a dipeptidyl-based inhibitor, generating rigid proline-based analogs. Accordingly, with these new proline-based peptidomimetics inhibitors, we selectively inhibited CatK against other human cathepsins (B, L and S). Among these new ligands, the most active one exhibited a high affinity (pKi = 7.3 – 50.1 nM) for CatK and no inhibition over the other cathepsins. This specific inhibitor harbors two novel substituents never employed in other CatK inhibitors: the trifluoromethylpyrazole and the 4-methylproline at P3 and P2 positions. These results broaden and advance the path toward new potent and selective inhibitors for CatK.

Computational evaluation and benchmark study of 342 crystallographic holo-structures of SARS-CoV-2 Mpro enzyme

The coronavirus disease 19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to a global health crisis with millions of confirmed cases and related deaths. The main protease (Mpro) of SARS-CoV-2 is crucial for viral replication and presents an attractive target for drug development. Despite the approval of some drugs, the search for effective treatments continues. In this study, we systematically evaluated 342 holo-crystal structures of Mpro to identify optimal conformations for structure-based virtual screening (SBVS). Our analysis revealed limited structural flexibility among the structures. Three docking programs, AutoDock Vina, rDock, and Glide were employed to assess the efficiency of virtual screening, revealing diverse performances across selected Mpro structures. We found that the structures 5RHE, 7DDC, and 7DPU (PDB Ids) consistently displayed the lowest EF, AUC, and BEDROCK scores. Furthermore, these structures demonstrated the worst pose prediction results in all docking programs. Two structural differences contribute to variations in docking performance: the absence of the S1 subsite in 7DDC and 7DPU, and the presence of a subpocket in the S2 subsite of 7DDC, 7DPU, and 5RHE. These findings underscore the importance of selecting appropriate Mpro conformations for SBVS, providing valuable insights for advancing drug discovery efforts.

A novel method for synthesizing authentic SARS-CoV-2 main protease

The SARS-CoV-2 main protease (Mpro) plays a crucial role in virus amplification and is an ideal target for antiviral drugs. Currently, authentic Mpro is prepared through two rounds of proteolytic cleavage. In this method, Mpro carries a self-cleavage site at the N-terminus and a protease cleavage site followed by an affinity tag at the C-terminus. This article proposes a novel method for producing authentic Mpro through single digestion. Mpro was constructed by fusing a His tag containing TEV protease cleavage sites at the N-terminus. The expressed recombinant protein was digested by TEV protease, and the generated protein had a decreased molecular weight and significantly increased activity, which was consistent with that of authentic Mpro generated by the previous method. These findings indicated that authentic Mpro was successfully obtained. Moreover, the substrate specificity of Mpro was investigated. Mpro had a strong preference for Phe at position the P2, which suggested that the S2 subsite was an outstanding target for designing inhibitors. This article also provides a reference for the preparation of Mpro for sudden coronavirus infection in the future.

Modulation of proteasome subunit selectivity of syringolins

Development of selective or dual proteasome subunit inhibitors based on syringolin B as a scaffold is described. We focused our efforts on a structure–activity relationship study of inhibitors with various substituents at the 3-position of the macrolactam moiety of syringolin B analogue to evaluate whether this would be sufficient to confer subunit selectivity by using sets of analogues with hydrophobic, basic and acidic substituents, which were designed to target Met45, Glu53 and Arg45 embedded in the S1 subsite, respectively. The structure–activity relationship study using systematic analogues provided insight into the origin of the subunit-selective inhibitory activity. This strategy would be sufficient to confer subunit selectivity regarding β5 and β2 subunits.

Collagen derived Gly-Pro-type DPP-IV inhibitory peptides: Structure-activity relationship, inhibition kinetics and inhibition mechanism

Our previous study has demonstrated that both the amino acid at N3 position and peptide length affected the DPP-IV inhibitory activity of Gly-Pro-type peptides. To further elucidate their molecular mechanism, a combined approach of QSAR modeling, enzymatic kinetics and molecular docking was used. Results showed that the QSAR models of Gly-Pro-type tripeptides and Gly-Pro-type peptides containing 3–12 residues were successfully constructed by 5z-scale descriptor with R2 of 0.830 and 0.797, respectively. The lower values of electrophilicity, polarity, and side-chain bulk of amino acid at N3 position caused higher DPP-IV inhibitory activity of Gly-Pro-type peptides. Moreover, an appropriate increase in the length of Gly-Pro-type peptides did not change their competitive inhibition mode, but decreased their inhibition constants (Ki values) and increased interactions with DPP-IV. More importantly, the interactions between the residues at C-terminal of Gly-Pro-type peptides containing 5 ∼ 6 residues with S2 extensive subsites (Ser209, Phe357, Arg358) of DPP-IV increased the interactions of Gly residue at N1 position with the S2 subsites (Glu205, Glu206, Asn710, Arg125, Tyr662) and decreased the acylation level of DPP-IV-peptide complex, and thereby increasing peptides’ DPP-IV inhibitory activity.

Design of substituted tetrahydrofuran derivatives for HIV-1 protease inhibitors: synthesis, biological evaluation, and X-ray structural studies.

Substituted tetrahydrofuran derivatives were designed and synthesized to serve as the P2 ligand for a series of potent HIV-1 protease inhibitors. Both enantiomers of the tetrahydrofuran derivatives were synthesized stereoselectivity in optically active forms using lipase-PS catalyzed enzymatic resolution as the key step. These tetrahydrofuran derivatives are designed to promote hydrogen bonding and van der Waals interactions with the backbone atoms in the S2 subsite of the HIV-1 protease active site. Several inhibitors displayed very potent HIV-1 protease inhibitory activity. A high-resolution X-ray crystal structure of an inhibitor-bound HIV-1 protease provided important insight into the ligand binding site interactions in the active site.

Allosteric modulation of serotonin and dopamine transporters: New insights from computations and experiments

Human monoamine transporters (MATs) are critical to regulating monoaminergic neurotransmission by translocating their substrates from the synaptic space back into the presynaptic neurons. As such, their primary substrate binding site S1 has been targeted by a wide range of compounds for treating neuropsychiatric and neurodegenerative disorders including depression, ADHD, neuropathic pain, and anxiety disorders. We present here a comparative study of the structural dynamics and ligand-binding properties of two MATs, dopamine transporter (DAT) and serotonin transporter (SERT), with focus on the allosteric modulation of their transport function by drugs or substrates that consistently bind a secondary site S2, proposed to serve as an allosteric site. Our systematic analysis of the conformational space and dynamics of a dataset of 50 structures resolved for DAT and SERT in the presence of one or more ligands/drugs reveals the specific residues playing a consistent role in coordinating the small molecules bound to subsites S2–I and S2-II within S2, such as R476 and Y481 in dDAT and E494, P561, and F556 in hSERT. Further analysis reveals how DAT and SERT differ in their two principal modes of structural changes, PC1 and PC2. Notably, PC1 underlies the transition between outward- and inward-facing states of the transporters as well as their gating; whereas PC2 supports the rearrangements of TM helices near the S2 site. Finally, the examination of cross-correlations between structural elements lining the respective sites S1 and S2 point to the crucial role of coupled motions between TM6a and TM10. In particular, we note the involvement of hSERT residues F335 and G338, and E493-E494-T497 belonging to these two respective helices, in establishing the allosteric communication between S1 and S2. These results help understand the molecular basis of the action of drugs that bind to the S2 site of DAT or SERT. They also provide a basis for designing allosteric modulators that may provide better control of specific interactions and cellular pathways, rather than indiscriminately inhibiting the transporter by targeting its orthosteric site.

Analysis of some biochemical properties of recombinant siberian roe deer (Capreolus pygargus) chymosin obtained in the mammalian cell culture (CHO-K1)

The structure of the Siberian roe deer (Capreolus pygargus) chymosin gene has been established for the first time and its exon/intron organization has been determined. The coding part of the C. pygargus chymosin gene was reconstructed and obtained as a DNA clone using the Golden Gate method. Comparative analysis of the sequences of prochymosins of roe deer, cow and single-humped camel revealed a number of amino acid substitutions in the sites forming the substrate-binding cavity of the enzyme and affecting the specificity subsites S4 and S1′ + S3′. The recombinant plasmid pIP1-Cap was constructed using the integration vector pIP1 for the expression of the roe deer prochymosin gene in CHO-K1 cells. A polyclone of CHO-K1-CYM-Cap cells was obtained, providing synthesis and secretion of recombinant prochymosin into the culture fluid of the producer. As a result of zymogen activation, a recombinant roe deer chymosin preparation with a total milk-clotting activity of 468.4 ± 11.1 IMCU/ml was obtained. The yield of recombinant roe deer chymosin was 500 mg/liter or ≈ 468,000 IMCU/liter, which exceeds the yield of genetically engineered chymosins in most of the expression systems used. The main biochemical properties of the obtained enzyme were compared with commercial preparations of recombinant chymosins of single-humped camel (Camelus dromedarius) and cow (Bos taurus). The specific milk-clotting activity of recombinant C. pygargus chymosin was 938 ± 22 IMCU/mg of protein and was comparable with the indicators of comparison enzymes. The nonspecific proteolytic activity of recombinant roe deer chymosin was 1.4-4.5 times higher than that of cow and camel enzymes. In terms of coagulation specificity, the recombinant C. pygargus chymosin occupied an intermediate position between the genetically engineered analogues of B. taurus and C. dromedarius chymosins. The threshold of thermal stability of recombinant roe deer chymosin was equal to 55°C. At 60°C, the enzyme retained

AlphaFold modeling of nepovirus 3C-like proteinases provides new insights into their diverse substrate specificities

The majority of picornaviral 3C proteinases (3Cpro) cleavage sites possess glutamine at the P1 position. Plant nepovirus 3C-like proteinases (3CLpro) show however much broader specificity, cleaving not only after glutamine, but also after several basic and hydrophobic residues. To investigate this difference, we employed AlphaFold to generate structural models of twelve selected 3CLpro, representing six substrate specificities. Generally, we observed favorable correlations between the architecture and charge of nepovirus proteinase S1 subsites and their ability to accept or restrict larger residues. The models identified a conserved aspartate residue close to the P1 residue in the S1 subsites of all nepovirus proteinases examined, consistent with the observed strong bias against negatively-charged residues at the P1 position of nepovirus cleavage sites. Finally, a cramped S4 subsite along with the presence of two unique histidine and serine residues explains the strict requirement of the grapevine fanleaf virus proteinase for serine at the P4 position.

Structure-Based Optimization and Biological Evaluation of Potent and Selective MMP-7 Inhibitors for Kidney Fibrosis.

Matrix metalloproteinase-7 (MMP-7) has been shown to play important roles in pathophysiological processes involved in the development/progression of diseases such as cancer and fibrosis. We discovered selective MMP-7 inhibitors composed of arylsulfonamide, carboxylate, and short peptides by a molecular hybridization approach. These compounds interacted with MMP-7 via multiple hydrogen bonds in the cocrystal structures. To obtain compounds for in vivo evaluation, we attempted structural optimization, particularly targeting Tyr167 at the S3 subsite through structure-based drug design, and identified compound 15 as showing improved MMP-7 potency and MMP subtype selectivity. A novel π-π stacking interaction with Tyr167 was achieved when 4-pyridylalanine was introduced as the P3 residue. Compound 15 suppressed the progression of kidney fibrosis in a dose-dependent manner in a mouse model of unilateral ureteral obstruction. Thus, we demonstrated, for the first time, that potent and selective MMP-7 inhibitors could prevent the progression of kidney fibrosis.

Molecular Dynamics Simulation Study of the Selective Inhibition of Coagulation Factor IXa over Factor Xa.

Thromboembolic disorders, arising from abnormal coagulation, pose a significant risk to human life in the modern world. The FDA has recently approved several anticoagulant drugs targeting factor Xa (FXa) to manage these disorders. However, these drugs have potential side effects, leading to bleeding complications in patients. To mitigate these risks, coagulation factor IXa (FIXa) has emerged as a promising target due to its selective regulation of the intrinsic pathway. Due to the high structural and functional similarities of these coagulation factors and their inhibitor binding modes, designing a selective inhibitor specifically targeting FIXa remains a challenging task. The dynamic behavior of protein-ligand interactions and their impact on selectivity were analyzed using molecular dynamics simulation, considering the availability of potent and selective compounds for both coagulation factors and the co-crystal structures of protein-ligand complexes. Throughout the simulations, we examined ligand movements in the binding site, as well as the contact frequencies and interaction fingerprints, to gain insights into selectivity. Interaction fingerprint (IFP) analysis clearly highlights the crucial role of strong H-bond formation between the ligand and D189 and A190 in the S1 subsite for FIXa selectivity, consistent with our previous study. This dynamic analysis also reveals additional FIXa-specific interactions. Additionally, the absence of polar interactions contributes to the selectivity for FXa, as observed from the dynamic profile of interactions. A contact frequency analysis of the protein-ligand complexes provides further confirmation of the selectivity criteria for FIXa and FXa, as well as criteria for binding and activity. Moreover, a ligand movement analysis reveals key interaction dynamics that highlight the tighter binding of selective ligands to the proteins compared to non-selective and inactive ligands.

Molecular mechanisms of SARS-CoV-2 resistance to nirmatrelvir.

Nirmatrelvir is a specific antiviral drug that targets the main protease (Mpro) of SARS-CoV-2 and has been approved to treat COVID-191,2. As an RNA virus characterized by high mutation rates, whether SARS-CoV-2 will develop resistance to nirmatrelvir is a question of concern. Our previous studies have shown that several mutational pathways confer resistance to nirmatrelvir, but some result in a loss of viral replicative fitness, which is then compensatedfor by additional alterations3. The molecular mechanisms for this observed resistance areunknown. Here we combined biochemical and structural methods to demonstrate that alterations at the substrate-binding pocket of Mpro can allow SARS-CoV-2 to develop resistance to nirmatrelvir in two distinct ways. Comprehensive studies of thestructures of14 Mpro mutantsin complex with drugs or substrate revealed that alterations at the S1 and S4 subsites substantially decreasedthe level of inhibitor binding, whereas alterations at the S2 and S4' subsites unexpectedly increased protease activity. Both mechanisms contributed to nirmatrelvir resistance, withthe latter compensating for the loss in enzymatic activity of the former, which in turn accounted for the restoration of viral replicative fitness, as observed previously3. Such a profile was also observed for ensitrelvir, another clinically relevant Mpro inhibitor. These results shed light on the mechanisms by which SARS-CoV-2 evolves to develop resistance to the current generation of protease inhibitors and provide the basis for the design of next-generation Mpro inhibitors.

A water-mediated hydrogen-bond network facilitates an exceptionally strong preference for asparagine in the S2 subsite of dipeptidyl peptidase 7 from Stenotrophomonas maltophilia

A water-mediated hydrogen-bond network facilitates an exceptionally strong preference for asparagine in the S2 subsite of dipeptidyl peptidase 7 from <i>Stenotrophomonas maltophilia</i>

Structure of puromycin-sensitive aminopeptidase and polyglutamine binding.

Puromycin-sensitive aminopeptidase (E.C. 3.4.11.14, UniProt P55786), a zinc metallopeptidase belonging to the M1 family, degrades a number of bioactive peptides as well as peptides released from the proteasome, including polyglutamine. We report the crystal structure of PSA at 2.3 Ǻ. Overall, the enzyme adopts a V-shaped architecture with four domains characteristic of the M1 family aminopeptidases, but it is in a less compact conformation compared to most M1 enzymes of known structure. A microtubule binding sequence is present in a C-terminal HEAT repeat domain of the enzyme in a position where it might serve to mediate interaction with tubulin. In the catalytic metallopeptidase domain, an elongated active site groove lined with aromatic and hydrophobic residues and a large S1 subsite may play a role in broad substrate recognition. The structure with bound polyglutamine shows a possible interacting mode of this peptide, which is supported by mutation.

Identification of potential inhibitors of cholinergic and β-secretase enzymes from phytochemicals derived from Gongronema latifolium Benth leaf: an integrated computational analysis.

Neurodegenerative disorders (NDDs) are associated with increased activities of the brain acetylcholinesterase (AChE), butyrylcholinesterase (BChE) and β-secretase enzyme (BACE1). Inhibition of these enzymes affords therapeutic option for managing NDDs such as Alzheimer's disease (AD) and Parkinson's disease (PD). Although, Gongronema latifolium Benth (GL) has been widely documented in ethnopharmacological and scientific reports for the management of NDDs, there is paucity of information on its underlying mechanism and neurotherapeutic constituents. Herein, 152 previously reported Gongronema latifolium derived-phytochemicals (GLDP) were screened against hAChE, hBChE and hBACE-1 using molecular docking, molecular dynamics (MD) simulations, free energy of binding calculations and cluster analysis. The result of the computational analysis identified silymarin, alpha-amyrin and teraxeron with the highest binding energies (-12.3, -11.2, -10.5 Kcal/mol) for hAChE, hBChE and hBACE-1 respectively as compared with those of the reference inhibitors (-12.3, -9.8 and - 9.4 for donepezil, propidium and aminoquinoline compound respectively). These best docked phytochemicals were found to be orientated in the hydrophobic gorge where they interacted with the choline-binding pocket in the A-site and P-site of the cholinesterase and subsites S1, S3, S3' and flip (67-75) residues of the pocket of the BACE-1. The best docked phytochemicals complexed with the target proteins were stable in a 100 ns molecular dynamic simulation. The interactions with the catalytic residues were preserved during the simulation as observed from the MMGBSA decomposition and cluster analyses. The presence of these phytocompounds most notably silymarin, which demonstrated dual high binding tendencies to both cholinesterases, were identified as potential neurotherapeutics subject to further investigation.