Computer-based de novo designs of tripeptides as novel neuraminidase inhibitors

The latest influenza A (H1N1) pandemic attracted worldwide attention and called for the urgent development of novel antiviral drugs. Here, seven tripeptides are designed and explored as neuraminidase (NA) inhibitors on the structural basis of known inhibitors. Their interactions with NA are studied and compared with each other, using flexible docking and molecular dynamics simulations. The various composed tripeptides have respective binding specificities and their interaction energies with NA decrease in the order of FRI > FRV > FRT > FHV > FRS > FRG > YRV (letters corresponding to amino acid code). The Arg and Phe portions of the tripeptides play important roles during the binding process: Arg has strong electrostatic interactions with the key residues Asp151, Glu119, Glu227 and Glu277, whereas Phe fits well in the hydrophobic cave within the NA active site. Owing to the introduction of hydrophobic property, the interaction energies of FRV and FRI are larger; in particular, FRI demonstrates the best binding quality and shows potential as a lead compound. In addition, the influence of the chemical states of the terminal amino acids are clarified: it is revealed that the charged states of the N-terminus (NH(3) (+)) and C-terminus (COO(-)) are crucial for the tripeptide inhibitory activities and longer peptides may not be appropriate. In addition, the medium inhibiting activity by acetylation of the N-terminus indicates the possible chemical modifications of FRI. Experimental efforts are expected in order to actualize the tripeptides as potent NA inhibitors in the near future.     

Dynophore-Based Approach in Virtual Screening: A Case of Human DNA Topoisomerase IIα

In this study, we utilized human DNA topoisomerase IIα as a model target to outline a dynophore-based approach to catalytic inhibitor design. Based on MD simulations of a known catalytic inhibitor and the native ATP ligand analog, AMP-PNP, we derived a joint dynophore model that supplements the static structure-based-pharmacophore information with a dynamic component. Subsequently, derived pharmacophore models were employed in a virtual screening campaign of a library of natural compounds. Experimental evaluation identified flavonoid compounds with promising topoisomerase IIα catalytic inhibition and binding studies confirmed interaction with the ATPase domain. We constructed a binding model through docking and extensively investigated it with molecular dynamics MD simulations, essential dynamics, and MM-GBSA free energy calculations, thus reconnecting the new results to the initial dynophore-based screening model. We not only demonstrate a new design strategy that incorporates a dynamic component of molecular recognition, but also highlight new derivates in the established flavonoid class of topoisomerase II inhibitors.     

Molecular Docking Studies and Biological Evaluation of Berberine–Benzothiazole Derivatives as an Anti-Influenza Agent via Blocking of Neuraminidase

In this study, we have introduced newly synthesized substituted benzothiazole based berberine derivatives that have been analyzed for their in vitro and in silico biological properties. The activity towards various kinds of influenza virus strains by employing the cytopathic effect (CPE) and sulforhodamine B (SRB) assay. Several berberine–benzothiazole derivatives (BBDs), such as BBD1, BBD3, BBD4, BBD5, BBD7, and BBD11, demonstrated interesting anti-influenza virus activity on influenza A viruses (A/PR/8/34, A/Vic/3/75) and influenza B viral (B/Lee/40, and B/Maryland/1/59) strain, respectively. Furthermore, by testing neuraminidase activity (NA) with the neuraminidase assay kit, it was identified that BBD7 has potent neuraminidase activity. The molecular docking analysis further suggests that the BBD1–BBD14 compounds’ antiviral activity may be because of interaction with residues of NA, and the same as in oseltamivir.     

Molecular Mechanism of Food-Derived Polyphenols on PD-L1 Dimerization: A Molecular Dynamics Simulation Study

In cancer immunotherapy, an emerging approach is to block the interactions of programmed cell death-1 (PD-1) and programmed cell death-ligand 1 (PD-L1) using small-molecule inhibitors. The food-derived polyphenols curcumin (CC), resveratrol (RSV) and epigallocatechin gallate (EGCG) have anticancer immunologic functions, which, recently, have been proposed to act via the downregulation of PD-L1 expression. However, it remains unclear whether they can directly target PD-L1 dimerization and, thus, interrupt the PD-1/PD-L1 pathway. To elucidate the molecular mechanism of such compounds on PD-L1 dimerization, molecular docking and nanosecond molecular dynamics simulations were performed. Binding free energy calculations show that the affinities of CC, RSV and EGCG to the PD-L1 dimer follow a trend of CC > RSV > EGCG. Hence, CC is the most effective inhibitor of the PD-1/PD-L1 pathway. Analysis on contact numbers, nonbonded interactions and residue energy decomposition indicate that such compounds mainly interact with the C-, F- and G-sheet fragments of the PD-L1 dimer, which are involved in interactions with PD-1. More importantly, nonpolar interactions between these compounds and the key residues Ile54, Tyr56, Met115, Ala121 and Tyr123 play a dominant role in binding. Free energy landscape and secondary structure analyses further demonstrate that such compounds can stably interact with the binding domain of the PD-L1 dimer. The results provide evidence that CC, RSV and EGCG can inhibit PD-1/PD-L1 interactions by directly targeting PD-L1 dimerization. This provides a novel approach to discovering food-derived small-molecule inhibitors of the PD-1/PD-L1 pathway with potential applications in cancer immunotherapy.     

Comparative analysis of putative agonist-binding modes in the human A1 adenosine receptor

A recent study reported a model of the human A(1) adenosine receptor and its agonist binding site, proposing two putative binding modes in the same binding site for the natural agonist, adenosine. The present work investigates the flexibility of this binding site by exhaustive exploration with the natural agonist and with three other adenosine derivatives: N6-cyclopentyladenosine (CPA), 2-chloro-N6-cyclopentyladenosine (CCPA), and 5'-N-ethylcarboxamidoadenosine (NECA). Our aim was to find a common binding mode for agonists that would explain the role in the binding process of the different substitutions allowed at the 2, N6, and 5' positions of adenosine. This problem was addressed through docking simulations, molecular dynamics studies, and estimations of the ligand-binding free energy with both the AUTODOCK scoring function and the linear interaction energy (LIE) approach. The results point to a single receptor-binding position that explains the effects of the different chemical modifications on the adenosine derivatives considered here.     

The Fight against the Influenza A Virus H1N1: Synthesis,Molecular Modeling,and Biological Evaluation of Benzofurazan Derivatives as Viral RNA Polymerase Inhibitors

The influenza RNA polymerase complex, which consists of the three subunits PA, PB1, and PB2, is a promising target for the development of new antiviral drugs. A large library of benzofurazan compounds was synthesized and assayed against influenza virus A/WSN/33 (H1N1). Most of the new derivatives were found to act by inhibiting the viral RNA polymerase complex through disruption of the complex formed between subunits PA and PB1. Docking studies were also performed to elucidate the binding mode of benzofurazans within the PB1 binding site in PA and to identify amino acids involved in their mechanism of action. The predicted binding pose is fully consistent with the biological data and lays the foundation for the rational development of more effective PA–PB1 inhibitors.     

Valproic Acid as a Potential Inhibitor of Plasmodium falciparum Histone Deacetylase 1 (PfHDAC1): An in Silico Approach

A new Plasmodium falciparum histone deacetylase1 (PfHDAC1) homology model was built based on the highest sequence identity available template human histone deacetylase 2 structure. The generated model was carefully evaluated for stereochemical accuracy, folding correctness and overall structure quality. All evaluations were acceptable and consistent. Docking a group of hydroxamic acid histone deacetylase inhibitors and valproic acid has shown binding poses that agree well with inhibitor-bound histone deacetylase-solved structural interactions. Docking affinity dG scores were in agreement with available experimental binding affinities. Further, enzyme-ligand complex stability and reliability were investigated by running 5-nanosecond molecular dynamics simulations. Thorough analysis of the simulation trajectories has shown that enzyme-ligand complexes were stable during the simulation period. Interestingly, the calculated theoretical binding energies of the docked hydroxamic acid inhibitors have shown that the model can discriminate between strong and weaker inhibitors and agrees well with the experimental affinities reported in the literature. The model and the docking methodology can be used in screening virtual libraries for PfHDAC1 inhibitors, since the docking scores have ranked ligands in accordance with experimental binding affinities. Valproic acid calculated theoretical binding energy suggests that it may inhibit PfHDAC1.     

Molecular Dynamics Simulation of Tryptophan Hydroxylase-1: Binding Modes and Free Energy Analysis to Phenylalanine Derivative Inhibitors

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.     

Discovery of Tyrosine Kinase Inhibitors by Docking into an Inactive Kinase Conformation Generated by Molecular Dynamics

Several small molecules that bind to the inactive DFG‐out conformation of tyrosine kinases (called type II inhibitors) have shown a good selectivity profile over other kinase targets. To obtain a set of DFG‐out structures, we performed an explicit solvent molecular dynamics (MD) simulation of the complex of the catalytic domain of a tyrosine kinase receptor, ephrin type‐A receptor 3 (EphA3), and a manually docked type II inhibitor. Automatic docking of four previously reported type II inhibitors was used to select a single snapshot from the MD trajectory for virtual screening. High‐throughput docking of a pharmacophore‐tailored library of 175 000 molecules resulted in about 4 million poses, which were further filtered by van der Waals efficiency and ranked according to a force‐field‐based energy function. Notably, around 20 % of the compounds with predicted binding energy smaller than ?10 kcal mol^?1 are known type II inhibitors. Moreover, a series of 5‐(piperazine‐1‐yl)isoquinoline derivatives was identified as a novel class of low‐micromolar inhibitors of EphA3 and unphosphorylated Abelson tyrosine kinase (Abl1). The in silico predicted binding mode of the new inhibitors suggested a similar affinity to the gatekeeper mutant T315I of Abl1, which was verified in vitro by using a competition binding assay. Additional evidence for the type II binding mode was obtained by two 300 ns MD simulations of the complex between N‐(3‐chloro‐4‐(difluoromethoxy)phenyl)‐2‐(4‐(8‐nitroisoquinolin‐5‐yl)piperazin‐1‐yl)acetamide and EphA3.     

Exploring the Molecular Interactions of 7,8-Dihydroxyflavone and Its Derivatives with TrkB and VEGFR2 Proteins

7,8-Dihydroxyflavone (7,8-DHF) is a TrkB receptor agonist, and treatment with this flavonoid derivative brings about an enhanced TrkB phosphorylation and promotes downstream cellular signalling. Flavonoids are also known to exert an inhibitory effect on the vascular endothelial growth factor receptor (VEGFR) family of tyrosine kinase receptors. VEGFR2 is one of the important receptors involved in the regulation of vasculogenesis and angiogenesis and has also been implicated to exhibit various neuroprotective roles. Its upregulation and uncontrolled activity is associated with a range of pathological conditions such as age-related macular degeneration and various proliferative disorders. In this study, we investigated molecular interactions of 7,8-DHF and its derivatives with both the TrkB receptor as well as VEGFR2. Using a combination of molecular docking and computational mapping tools involving molecular dynamics approaches we have elucidated additional residues and binding energies involved in 7,8-DHF interactions with the TrkB Ig2 domain and VEGFR2. Our investigations have revealed for the first time that 7,8-DHF has dual biochemical action and its treatment may have divergent effects on the TrkB via its extracellular Ig2 domain and on the VEGFR2 receptor through the intracellular kinase domain. Contrary to its agonistic effects on the TrkB receptor, 7,8-DHF was found to downregulate VEGFR2 phosphorylation both in 661W photoreceptor cells and in retinal tissue.     

Insights into the Interaction of LVV-Hemorphin-7 with Angiotensin II Type 1 Receptor

Hemorphins are known for their role in the control of blood pressure. Recently, we revealed the positive modulation of the angiotensin II (AngII) type 1 receptor (AT1R) by LVV-hemorphin-7 (LVV-H7) in human embryonic kidney (HEK293) cells. Here, we examined the molecular binding behavior of LVV-H7 on AT1R and its effect on AngII binding using a nanoluciferase-based bioluminescence resonance energy transfer (NanoBRET) assay in HEK293FT cells, as well as molecular docking and molecular dynamics (MD) studies. Saturation and real-time kinetics supported the positive effect of LVV-H7 on the binding of AngII. While the competitive antagonist olmesartan competed with AngII binding, LVV-H7 slightly, but significantly, decreased AngII’s k_D by 2.6 fold with no effect on its B_max. Molecular docking and MD simulations indicated that the binding of LVV-H7 in the intracellular region of AT1R allosterically potentiates AngII binding. LVV-H7 targets residues on intracellular loops 2 and 3 of AT1R, which are known binding sites of allosteric modulators in other GPCRs. Our data demonstrate the allosteric effect of LVV-H7 on AngII binding, which is consistent with the positive modulation of AT1R activity and signaling previously reported. This further supports the pharmacological targeting of AT1R by hemorphins, with implications in vascular and renal physiology.     

Pentacycloundecane-diol-based HIV-1 protease inhibitors: biological screening, 2D NMR, and molecular simulation studies

Novel compounds incorporating a pentacycloundecane (PCU) diol moiety were designed, synthesized, and evaluated as inhibitors of the wild-type C-South African (C-SA) HIV-1 protease. Seven compounds are reported herein, three of which displayed IC(50) values in the 0.5-0.6 μM range. The cytotoxicity of PCU cage peptides toward human MT-4 cells appears to be several orders of magnitude less toxic than the current antiviral medications ritonavir and lopinavir. NMR studies based on the observed through-space (1)H,(1)H distances/contacts in the EASY-ROESY spectra of three of the considered PCU peptide inhibitors enabled us to describe their secondary solution structure. Conserved hydrogen bonding interactions were observed between the hydroxy group of the PCU diol inhibitors and the catalytic triad (Asp25, Ile26, Gly27) of HIV protease in docking and molecular dynamics simulations. The biological significance and possible mode of inhibition by PCU-based HIV protease inhibitors discussed herein facilitates a deeper understanding of this family of inhibitors and their potential application to a vast number of alternative diseases related to proteases.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.     

Analysis of Binding Determinants for Different Classes of Competitive and Noncompetitive Inhibitors of Glycine Transporters

Glycine transporters are interesting therapeutic targets as they play significant roles in glycinergic and glutamatergic systems. The search for new selective inhibitors of particular types of glycine transporters (GlyT-1 and GlyT-2) with beneficial kinetics is hampered by limited knowledge about the spatial structure of these proteins. In this study, a pool of homology models of GlyT-1 and GlyT-2 in different conformational states was constructed using the crystal structures of related transporters from the SLC6 family and the recently revealed structure of GlyT-1 in the inward-open state, in order to investigate their binding sites. The binding mode of the known GlyT-1 and GlyT-2 inhibitors was determined using molecular docking studies, molecular dynamics simulations, and MM-GBSA free energy calculations. The results of this study indicate that two amino acids, Gly373 and Leu476 in GlyT-1 and the corresponding Ser479 and Thr582 in GlyT-2, are mainly responsible for the selective binding of ligands within the S1 site. Apart from these, one pocket of the S2 site, which lies between TM3 and TM10, may also be important. Moreover, selective binding of noncompetitive GlyT-1 inhibitors in the intracellular release pathway is affected by hydrophobic interactions with Ile399, Met382, and Leu158. These results can be useful in the rational design of new glycine transporter inhibitors with desired selectivity and properties in the future.     

Influenza A Virus NS1 Protein Structural Flexibility Analysis According to Its Structural Polymorphism Using Computational Approaches

Influenza A viruses are highly contagious RNA viruses that cause respiratory tract infections in humans and animals. Their non-structural protein NS1, a homodimer of two 230-residue chains, is the main viral factor in counteracting the antiviral defenses of the host cell. Its RNA-binding domain is an obligate dimer that is connected to each of the two effector domains by a highly flexible unstructured linker region of ten amino acids. The flexibility of NS1 is a key property that allows its effector domains and its RNA binding domain to interact with several protein partners or RNAs. The three-dimensional structures of full-length NS1 dimers revealed that the effector domains could adopt three distinct conformations as regards their mutual interactions and their orientation relative to the RNA binding domain (closed, semi-open and open). The origin of this structural polymorphism is currently being investigated and several hypotheses are proposed, among which one posits that it is a strain-specific property. In the present study, we explored through computational molecular modeling the dynamic and flexibility properties of NS1 from three important influenza virus A strains belonging to three distinct subtypes (H1N1, H6N6, H5N1), for which at least one conformation is available in the Protein Data Bank. In order to verify whether NS1 is stable in three forms for the three strains, we constructed homology models if the corresponding forms were not available in the Protein Data Bank. Molecular dynamics simulations were performed in order to predict the stability over time of the three distinct sequence variants of NS1, in each of their three distinct conformations. Our results favor the co-existence of three stable structural forms, regardless of the strain, but also suggest that the length of the linker, along with the presence of specific amino acids, modulate the dynamic properties and the flexibility of NS1.     

Molecular Docking and Molecular Dynamics Simulation Studies of Triterpenes from Vernonia patula with the Cannabinoid Type 1 Receptor

A molecular docking approach was employed to evaluate the binding affinity of six triterpenes, namely epifriedelanol, friedelin, α-amyrin, α-amyrin acetate, β-amyrin acetate, and bauerenyl acetate, towards the cannabinoid type 1 receptor (CB1). Molecular docking studies showed that friedelin, α-amyrin, and epifriedelanol had the strongest binding affinity towards CB1. Molecular dynamics simulation studies revealed that friedelin and α-amyrin engaged in stable non-bonding interactions by binding to a pocket close to the active site on the surface of the CB1 target protein. The studied triterpenes showed a good capacity to penetrate the blood–brain barrier. These results help to provide some evidence to justify, at least in part, the previously reported antinociceptive and sedative properties of Vernonia patula.     

Computational Investigation Identified Potential Chemical Scaffolds for Heparanase as Anticancer Therapeutics

Heparanase (Hpse) is an endo-β-D-glucuronidase capable of cleaving heparan sulfate side chains. Its upregulated expression is implicated in tumor growth, metastasis and angiogenesis, thus making it an attractive target in cancer therapeutics. Currently, a few small molecule inhibitors have been reported to inhibit Hpse, with promising oral administration and pharmacokinetic (PK) properties. In the present study, a ligand-based pharmacophore model was generated from a dataset of well-known active small molecule Hpse inhibitors which were observed to display favorable PK properties. The compounds from the InterBioScreen database of natural (69,034) and synthetic (195,469) molecules were first filtered for their drug-likeness and the pharmacophore model was used to screen the drug-like database. The compounds acquired from screening were subjected to molecular docking with Heparanase, where two molecules used in pharmacophore generation were used as reference. From the docking analysis, 33 compounds displayed higher docking scores than the reference and favorable interactions with the catalytic residues. Complex interactions were further evaluated by molecular dynamics simulations to assess their stability over a period of 50 ns. Furthermore, the binding free energies of the 33 compounds revealed 2 natural and 2 synthetic compounds, with better binding affinities than reference molecules, and were, therefore, deemed as hits. The hit compounds presented from this in silico investigation could act as potent Heparanase inhibitors and further serve as lead scaffolds to develop compounds targeting Heparanase upregulation in cancer.