In silico discovery of a compound with nanomolar affinity to antithrombin causing partial activation and increased heparin affinity

Computational analyses reveal fundamental properties of the AT structure related to thrombosis

Summary

Blood coagulation is a vital process for humans and other species. Following an injury to a blood vessel, a cascade of molecular signals is transmitted, inhibiting and activating more than a dozen coagulation factors and resulting in the formation of a fibrin clot that ceases the bleeding. In this process, antithrombin (AT), encoded by the SERPINC1 gene is a key player regulating the clotting activity and ensuring that it stops at the right time. In this sense, mutations to this factor often result in thrombosis-the excessive coagulation that leads to the potentially fatal formation of blood clots that obstruct veins. Although this process is well known, it is still unclear why even single residue substitutions to AT lead to drastically different phenotypes. In this study, to understand the effect of mutations throughout the AT structure, we created a detailed network map of this protein, where each node is an amino acid, and two amino acids are connected if they are in close proximity in the three-dimensional structure. With this simple and intuitive representation and a machine-learning framework trained using genetic information from more than 130 patients, we found that different types of thrombosis have emerging patterns that are readily identifiable. Together, these results demonstrate how clinical features, genetic data and in silico analysis are converging to enhance the diagnosis and treatment of coagulation disorders.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years

Structure-based virtual screening (SBVS), also known as molecular docking, has been increasingly applied to discover small-molecule ligands based on the protein structures in the early stage of drug discovery. In this review, we comprehensively surveyed the prospective applications of molecular docking judged by solid experimental validations in the literature over the past fifteen years. Herein, we systematically analyzed the novelty of the targets and the docking hits, practical protocols of docking screening, and the following experimental validations. Among the 419 case studies we reviewed, most virtual screenings were carried out on widely studied targets, and only 22% were on less-explored new targets. Regarding docking software, GLIDE is the most popular one used in molecular docking, while the DOCK 3 series showed a strong capacity for large-scale virtual screening. Besides, the majority of identified hits are promising in structural novelty and one-quarter of the hits showed better potency than 1 μM, indicating that the primary advantage of SBVS is to discover new chemotypes rather than highly potent compounds. Furthermore, in most studies, only in vitro bioassays were carried out to validate the docking hits, which might limit the further characterization and development of the identified active compounds. Finally, several successful stories of SBVS with extensive experimental validations have been highlighted, which provide unique insights into future SBVS drug discovery campaigns.

Structure based virtual screening of natural product molecules as glycosidase inhibitors

Objective of the present investigation comprised of the application of in silico methods to discover novel natural product (NP) based potential inhibitors for carbohydrate mediated diseases. Structure based drug design studies (molecular docking and structure based pharmacophore analysis)  were carried out on a series of natural product compounds to identify significant bioactive molecules to inhibit α-mannosidase (I and II) and β-galactosidase enzymes. Furthermore, protein ligand interaction fingerprint analysis, molecular dynamics simulations and molecular access system (MACCS) fingerprint analysis were performed to understand the binding behaviors of the studied molecules. The results derived from these analyses showed that the identified compounds exhibit significant binding interactions with the active site residues. The compounds, NP-51, NP-81 and NP-165 have shown significant docking score against the studied enzymes (α-mannosidases-I, α-mannosidases-II and β-galactosidases). The fingerprint studies showed that the presence of rings (aromatic or aliphatic) with sulfur atoms, nitrogen atoms, methyl groups, etc. have favorable effects on the α-mannosidase II inhibitory activity. However, the presence of halogen atoms substituted in the molecules have reduced inhibitory ability against α-mannosidase II. The compound, NP-165 has significant activity against both enzymes (α-mannosidases and β-galactosidases). These studies accomplished that the compounds identified through in silico methodologies can be used to develop semisynthetic derivatives of the glycosidase inhibitors and can be screened for the treatment of different carbohydrate mediated diseases.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-021-00115-9.

Small-molecule modulators of serine protease inhibitor proteins (serpins)

Serine protease inhibitors (serpins) are a large family of proteins that regulate and control crucial physiological processes, such as inflammation, coagulation, thrombosis and thrombolysis, and immune responses. The extraordinary impact that these proteins have on numerous crucial pathways makes them an attractive target for drug discovery. In this review, we discuss recent advances in research on small-molecule modulators of serpins, examine their mode of action, analyse the structural data from crystallised protein-ligand complexes, and highlight the potential obstacles and possible therapeutic perspectives. The application of in silico methods for rational drug discovery is also summarised. In addition, we stress the need for continued research in this field.

Potential inhibitors of the enzyme acetylcholinesterase and juvenile hormone with insecticidal activity: study of the binding mode via docking and molecular dynamics simulations

Models validation in QSAR, pharmacophore, docking and others can ensure the accuracy and reliability of future predictions in design and selection of molecules with biological activity. In this study, pyriproxyfen was used as a pivot/template to search the database of the Maybridge Database for potential inhibitors of the enzymes acetylcholinesterase and juvenile hormone as well. The initial virtual screening based on the 3D shape resulted in 2000 molecules with Tanimoto index ranging from 0.58 to 0.88. A new reclassification was performed on the overlapping of positive and negative charges, which resulted in 100 molecules with Tanimoto's electrostatic score ranging from 0.627 to 0.87. Using parameters related to absorption, distribution, metabolism and excretion and the pivot molecule, the molecules selected in the previous stage were evaluated regarding these criteria, and 21 were then selected. The pharmacokinetic and toxicological properties were considered and for 12 molecules, the DEREK software not fired any alert of toxicity, which were thus considered satisfactory for prediction of biological activity using the Web server PASS. In the molecular docking with insect acetylcholinesterase, the Maybridge3_002654 molecule had binding affinity of -11.1 kcal/mol, whereas in human acetylcholinesterase, the Maybridge4_001571molecule show in silico affinity of -10.2 kcal/mol, and in the juvenile hormone, the molecule MCULE-8839595892 show in silico affinity value of -11.6 kcal/mol. Subsequent long-trajectory molecular dynamics studies indicated considerable stability of the novel molecules compared to the controls.AbbreviationsQSARquantitative structure-activity relationshipsPASSprediction of activity spectra for substancesCommunicated by Ramaswamy H. Sarma.

Size and Flexibility Define the Inhibition of the H3N2 Influenza Endonuclease Enzyme by Calix[n]arenes

Inhibition of H3N2 influenza PA endonuclease activity by a panel of anionic calix[n]arenes and β-cyclodextrin sulfate has been studied. The joint experimental and theoretical results reveal that the larger, more flexible and highly water-soluble sulfonato-calix[n]arenes have high inhibitory activity, with para-sulfonato-calix[8]arene, SC8, having an IC₅₀ value of 6.4 μM. Molecular docking calculations show the SC8 can interact at both the polyanion binding site and also the catalytic site of H3N2 influenza PA endonuclease.

Anxiolytic effects of Formononetin in an inflammatory pain mouse model

Chronic pain is commonly accompanied with anxiety disorder, which complicates treatment. In this study, we investigated the analgesic and anxiolytic effects of Formononetin (FMNT), an active component of traditional Chinese medicine red clover (Trifolium pratense L.) that is capable of protecting neurons from N-methyl-D-aspartate (NMDA)-evoked excitotoxic injury, on mice suffering from complete Freund’s adjuvant (CFA)-induced chronic inflammatory pain. The results show that FMNT administration significantly reduces anxiety-like behavior but does not affect the nociceptive threshold in CFA-injected mice. The treatment reverses the upregulation of NMDA, GluA1, and GABA_A receptors, as well as PSD95 and CREB in the basolateral amygdala (BLA). The effects of FMNT on NMDA receptors and CREB binding protein (CBP) were further confirmed by the potential structure combination between these compounds, which was analyzed by in silico docking technology. FMNT also inhibits the activation of the NF-κB signaling pathway and microglia in the BLA of mice suffering from chronic inflammatory pain. Therefore, the anxiolytic effects of FMNT are partially due to the attenuation of inflammation and neuronal hyperexcitability through the inhibition of NMDA receptor and CBP in the BLA.

Identification of chemically diverse GABAA agonists as potential anti-epileptic agents using structure-guided virtual screening, ADMET, quantum mechanics and clinical validation through off-target analysis

Development of resistance against existing anti-epileptic drugs has alarmed the scientific innovators to find novel potential chemical starting points for the treatment of epilepsy and GABA_A inhibition is a promising drug target strategy against epilepsy. The crystal structure of a subtype-selective β3-homopentameric ligand-gated ion channel of GABA_A receptor has been used for the first time for screening the Asinex library for discovery of GABA_A agonists as potential anti-epileptic agents. Co-crystallized ligand established the involvement of part of the β7-β8 loop (Glu155 and Tyr157) and β9-β10 loop (Phe200 and Tyr205) residues as the crucial amino acids in effective binding, an essential feature, being hydrogen bond or ionic interaction with Glu155 residue. Top ranked hits were further subjected to binding energy estimation, ADMET analysis and ligand efficiency matric calculations as consecutive filters. About 19 compounds qualifying all parameters possessed interaction of one positively charged group with Glu155 with good CNS drug-like properties. Simulation studies were performed on the apo protein, its complex with co-crystallized ligand and the best hit qualifying all screening parameters. The best hit was also analyzed using Quantum mechanical studies, off-target analysis and hit modification. The off-target analysis emphasized that these agents did not have any other predicted side-effects.

The current limits in virtual screening and property prediction

Beyond finding inhibitors that show high binding affinity to the respective target, there is the challenge of optimizing their properties with respect to metabolic and toxicological issues, as well as further off-target effects. To reduce the experimental effort of synthesizing and testing actual substances in corresponding assays, virtual screening has become an indispensable toolbox in preclinical development. The scope of application covers the prediction of molecular properties including solubility, metabolic liability and binding to antitargets, such as the hERG channel. Furthermore, prediction of binding sites and drugable targets are emerging aspects of virtual screening. Issues involved with the currently applied computational models including machine learning algorithms are outlined, such as limitations to the accuracy of prediction and overfitting.

In Silico Designed Axl Receptor Blocking Drug Candidates Against Zika Virus Infection

After a large outbreak in Brazil, novel drugs against Zika virus became extremely necessary. Evaluation of virus-based pharmacological strategies concerning essential host factors brought us to the idea that targeting the Axl receptor by blocking its dimerization function could be critical for virus entry. Starting from experimentally validated compounds, such as RU-301, RU-302, warfarin, and R428, we identified a novel compound 2' (R428 derivative) to be the most potent for this task amongst a number of alternative compounds and leads. The improved affinity of compound 2' was confirmed by molecular docking as well as molecular dynamics simulation techniques using implicit solvation models. The current study summarizes a new possibility for inhibition of the Axl function as a potential target for future antiviral therapies.

Computational investigations of physicochemical, pharmacokinetic, toxicological properties and molecular docking of betulinic acid, a constituent of Corypha taliera (Roxb.) with Phospholipase A2 (PLA2)

BACKGROUND

Betulinic acid (BA) is a natural triterpenoid compound and exhibits a wide range of biological and medicinal properties including anti-inflammatory activity. Therefore, this theoretical investigation is performed to evaluate (a) physicochemical properties such as acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP), aqueous solubility (logS), solvation free energy, dipole moment, polarizability, hyperpolarizability and different reactivity descriptors, (b) pharmacokinetic properties like human intestinal absorption (HIA), cellular permeability, skin permeability (P_Skin), plasma protein binding (PPB), penetration of the blood brain barrier (BBB), (c) toxicological properties including mutagenicity, carcinogenicity, risk of inhibition of hERG gene and (d) molecular mechanism of anti-inflammatory action which will aid the development of analytical method and the synthesis of BA derivatives.

METHODS

The physicochemical properties were calculated using MarvinSketch 15.6.29 and Gaussian 09 software package. The pharmacokinetic and toxicological properties were calculated on online server PreADMET. Further, the molecular docking study was conducted on AutoDock vina in PyRx 0.8.

RESULTS

The aqueous solubility increased with increasing pH due to the ionization of BA leading to decrease in distribution coefficient. The solvation energies in water, dimethyl sulfoxide (DMSO), acetonitrile, n-octanol, chloroform and carbon tetrachloride were - 41.74 kJ/mol, - 53.80 kJ/mol, - 66.27 kJ/mol, - 69.64 kJ/mol, - 65.96 kJ/mol and - 60.13 kJ/mol, respectively. From the results of polarizability and softness, it was clear that BA is less stable and hence, kinetically more reactive in water. BA demonstrated good human intestinal absorption (HIA) and moderate cellular permeability. Further, BA also exhibited positive CNS activity due to high permeability through BBB. The toxicological study revealed that BA was a mutagenic compound but noncarcinogenic in mice model. Moreover, molecular docking study of BA with PLA2 revealed that BA interacts with GLY22 & GLY29 through hydrogen bond formation and LEU2, PHE5, HIS6, ALA17, ALA18, HIS47 and TYR51 through different types of hydrophobic interactions. The binding affinity of BA was - 41.00 kJ/mol which is comparable to the binding affinity of potent inhibitor 6-Phenyl-4(R)-(7-Phenyl-heptanoylamino)-hexanoic acid (BR4) (- 33.89 kJ/mol).

CONCLUSIONS

Our computed properties may assist the development of analytical method to assay BA or to develop BA derivatives with better pharmacokinetic and toxicological profile.

New lipid interaction partners stimulate the inhibition of activated protein C by cell-penetrating protein C inhibitor

Protein C inhibitor (PCI, SerpinA5) is a heparin-binding serpin which can penetrate through cellular membranes. Selected negatively charged phospholipids like unsaturated phosphatidylserine and oxidised phosphatidylethanolamine bind to PCI and stimulate its inhibitory activity towards different proteases. The interaction of phospholipids with PCI might also alter the lipid distribution pattern of blood cells and influence the remodelling of cellular membranes. Here we showed that PCI is an additional binding partner of phosphatidic acid (PA), cardiolipin (CL), and phosphoinositides (PIPs). Protein lipid overlay assays exhibited a unique binding pattern of PCI towards different lipid species. In addition PA, CL, and unsaturated, monophosphorylated PIPs stimulated the inhibitory property of PCI towards activated protein C in a heparin like manner. As shown for kallistatin (SerpinA4) and vaspin (SerpinA12), the incubation of cells with PCI led to the activation of protein kinase B (AKT), which could be achieved through direct interaction of PCI with PIPs. This model is supported by the fact that PCI stimulated the PIP-dependent 5-phosphatase SHIP2 in vitro, which would result in AKT activation. Hence the interaction of PCI with different lipids might not only stimulate the inhibition of potential target protease by PCI, but could also alter intracellular lipid signalling.

A Comprehensive Docking and MM/GBSA Rescoring Study of Ligand Recognition upon Binding Antithrombin

Background:

A high-throughput virtual screening pipeline has been extended from single energetically minimized structure Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) rescoring to ensemble-average MM/GBSA rescoring. The correlation coefficient (R2) of calculated and experimental binding free energies for a series of antithrombin ligands has been improved from 0.36 to 0.69 when switching from the single-structure MM/GBSA rescoring to ensemble-average one. The electrostatic interactions in both solute and solvent are identified to play an important role in determining the binding free energy after the decomposition of the calculated binding free energy. The increasing negative charge of the compounds provides a more favorable electrostatic energy change but creates a higher penalty for the solvation free energy. Such a penalty is compensated by the electrostatic energy change, which results in a better binding affinity. A highly hydrophobic ligand is determined by the docking program to be a non-specific binder.

Results:

Our results have demonstrated that it is very important to keep a few top poses for rescoring, if the binding is non-specific or the binding mode is not well determined by the docking calculation.

Environment sensitive fluorescent analogue of biologically active oxazoles differentially recognizes human serum albumin and bovine serum albumin: Photophysical and molecular modeling studies

An environment sensitive fluorophore, 4-(5-(4-(dimethylamino)phenyl)oxazol-2-yl)benzoic acid (DMOBA), that closely mimics biologically active 2,5-disubstituited oxazoles has been designed to probe two homologous serum proteins, human serum albumin (HSA) and bovine serum albumin (BSA) by means of photophysical and molecular modeling studies. This fluorescent analogue exhibits solvent polarity sensitive fluorescence due to an intramolecular charge transfer in the excited state. In comparison to water, the steady state emission spectra of DMOBA in BSA is characterized by a greater blue shift (~10nm) and smaller Stokes' shift (~5980cm^-1) in BSA than HSA (Stokes'shift~6600cm^-1), indicating less polar and more hydrophobic environment of the dye in the former than the latter. The dye-protein binding interactions are remarkably stronger for BSA than HSA which is evident from higher value of the association constant for the DMOBA-BSA complex (K_a~5.2×10⁶M^-1) than the DMOBA-HSA complex (K_a~1.0×10⁶M^-1). Fӧrster resonance energy transfer studies revealed remarkably less efficient energy transfer (8%) between the donor tryptophans in BSA and the acceptor DMOBA dye than that (30%) between the single tryptophan moiety in HSA and the dye, which is consistent with a much larger distance between the donor (tryptophan)-acceptor (dye) pair in BSA (34.5Å) than HSA (25.4Å). Site specific competitive binding assays have confirmed on the location of the dye in Sudlow's site II of BSA and in Sudlow's site I of HSA, respectively. Molecular modeling studies have shown that the fluorescent analogue is tightly packed in the binding site of BSA due to strong steric complementarity, where, binding of DMOBA to BSA is primarily dictated by the van der Waals and hydrogen bonding interactions. In contrast, in HSA the steric complementarity is less significant and binding is primarily guided by polar interactions and van der Waals interactions appear to be less significant in the formation of the HSA-DMOBA complex. Electrostatic interactions contribute significantly in the binding of DMOBA to HSA (-2.09kcal/mol) compared to BSA (-0.47kcal/mol). Electrostatic surface potential calculation reveals that the DMOBA binding site within HSA is highly charged compared to BSA.

In vitro modulatory effects of functionalized pyrimidines and piperidine derivatives on Aryl hydrocarbon receptor (AhR) and glucocorticoid receptor (GR) activities

The development of biologically active molecules based on molecular recognition is an attractive and challenging task in medicinal chemistry and the molecules that can activate/deactivate certain receptors are of great medical interest. In this contribution, selected pyrimidine/piperidine derivatives were synthesized and tested for the ability to activate/deactivate Aryl hydrocarbon receptor (AhR) and Glucocorticoid receptor (GR). Tested compounds are shown to activate the receptors but to much lesser extent than positive controls, dioxin and dexamethasone for Ahr and GR, respectively. However, some of them antagonized the positive controls action. Although further in vivo studies are needed to fully characterize the bioactivities of these compounds, the reported in vitro evidences demonstrate that they might be used as the modulators of AhR and GR activities.

Chitosan as stabilizing agent for negatively charged nanoparticles

Chitosan is a biocompatible polysaccharide with positive Z potential which can stabilize negative charged nanoparticles. Silk fibroin nanoparticles and citrate gold nanoparticles, both with negative Z potential, but they form aggregates at physiological ionic strength. In this work, we study the behavior of chitosan in solution when the ionic strength of the medium is increased and how the concentration of chitosan and the proportion of the two components (chitosan and AuNP or SFN) significantly affect the stability and size of the nanocomposites formed. In addition to experimental measurements, molecular modeling were used to gain insight into how chitosan interacts with silk fibroin monomers, and to identify the main energetic interactions involved in the process. The optimum values for obtaining the smallest and most homogeneous stable nanocomposites were obtained and two different ways of organization through which chitosan may exert its stabilizing effect were suggested.

Structural grounds for the 2-aminopurine mutagenicity: a novel insight into the old problem of the replication errors

Mutagenic pressure of the 2AP molecule on DNA during its replication is realized via the more intensive generation of the T* mutagenic tautomers through the reaction 2AP·T(WC) → 2AP·T*(w).

Labelling Herceptin with a novel oxaliplatin derivative: a computational approach towards the selective drug delivery

The clinical use of platinum(II)-based drugs has serious side effects due to the non-specific reactions with both malignant and normal cells. To circumvent such major drawback, novel metallodrugs might be combined with suitable carrier molecules, as antibodies, to ensure selective attacks on tumours while sparing healthy tissues. In this contribution, we investigate the stability of a novel oxaliplatin derivate drug embedded in Herceptin (trastuzumab), an antibody which is able to recognise breast cancer cells, by using a wide panel of theoretical tools: docking, molecular dynamics and quantum calculations. Our calculations reveal the binding mechanism: the drug initially interacts non-covalently with the Pro40A and Asp167A residues, and the nitrogen of His171B subsequently replaces one of the water molecules coordinated to the platinum center, where the latter step reversibly fixes the drug into the antibody. These data might be used to further rationalise the synthesis of improved drugs beyond classical platinum(II) derivatives by improving the ligand-protein coupling mode.

Heart health risk assessment system: a nonintrusive proposal using ontologies and expert rules

According to the World Health Organization, the world's leading cause of death is heart disease, with nearly two million deaths per year. Although some factors are not possible to change, there are some keys that help to prevent heart diseases. One of the most important keys is to keep an active daily life, with moderate exercise. However, deciding what a moderate exercise is or when a slightly abnormal heart rate value is a risk depends on the person and the activity. In this paper we propose a context-aware system that is able to determine the activity the person is performing in an unobtrusive way. Then, we have defined ontology to represent the available knowledge about the person (biometric data, fitness status, medical information, etc.) and her current activity (level of intensity, heart rate recommended for that activity, etc.). With such knowledge, a set of expert rules based on this ontology are involved in a reasoning process to infer levels of alerts or suggestions for the users when the intensity of the activity is detected as dangerous for her health. We show how this approach can be accomplished by using only everyday devices such as a smartphone and a smartwatch.

A performance/cost evaluation for a GPU-based drug discovery application on volunteer computing

Bioinformatics is an interdisciplinary research field that develops tools for the analysis of large biological databases, and, thus, the use of high performance computing (HPC) platforms is mandatory for the generation of useful biological knowledge. The latest generation of graphics processing units (GPUs) has democratized the use of HPC as they push desktop computers to cluster-level performance. Many applications within this field have been developed to leverage these powerful and low-cost architectures. However, these applications still need to scale to larger GPU-based systems to enable remarkable advances in the fields of healthcare, drug discovery, genome research, etc. The inclusion of GPUs in HPC systems exacerbates power and temperature issues, increasing the total cost of ownership (TCO). This paper explores the benefits of volunteer computing to scale bioinformatics applications as an alternative to own large GPU-based local infrastructures. We use as a benchmark a GPU-based drug discovery application called BINDSURF that their computational requirements go beyond a single desktop machine. Volunteer computing is presented as a cheap and valid HPC system for those bioinformatics applications that need to process huge amounts of data and where the response time is not a critical factor.

The inhibitory effects of bioactive compounds of tomato juice binding to hepatic HMGCR: in vivo study and molecular modelling

The hypocholesterolemic effect of tomato juice has been investigated in an intervention study with rats, along with the possible inhibition effect of bioactive tomato compounds binding to the HMGCR enzyme. Two experimental groups (n = 8 Sprague-Dawley rats) were fed ad libitum for five weeks, with water or tomato juice provided to the control and intervention groups, respectively. Total, LDL and HDL cholesterol, and total triglycerides were analysed in plasma, and the lycopene content and the expression and activity of the enzyme HMGCR were determined in liver samples. A computational molecular modelling was carried out to determine the interactions between HMGCR and lycopene, chlorogenic acid and naringenin. Total, LDL and HDL cholesterol were significantly lower in the intervention group after the intake of tomato juice. In addition, a significant reduction in HMGCR activity was observed, although this was not accompanied by changes in gene expression. The molecular modelling showed that components of tomato can bind to the active site of the enzyme and compete with the ligand HMGCoA. Lycopene, from tomato juice, accumulates in the liver and can inhibit the activity of the rate-limiting enzyme of cholesterol biosynthesis, HMGCR.