Ligand-based virtual screening, molecular dynamics, and biological evaluation of repurposed drugs as inhibitors of Trypanosoma cruzi proteasome.

Chagas disease is a well-known Neglected Tropical Disease, mostly endemic in continental Latin America, but that has spread to North America and Europe. Unfortunately, current treatments against such disease are ineffective and produce known and undesirable side effects. To find novel effective drug candidates to treat Chagas disease, we uniquely explore the Trypanosoma cruzi proteasome as a recent biological target and, also, apply drug repurposing through different computational methodologies. For this, we initially applied protein homology modeling to build a robust model of proteasome ß4/ß5 subunits, since there is no crystallographic structure of this target. Then, we used it on a drug repurposing via a virtual screening campaign starting with more than 8,000 drugs and including the methodologies: ligand-based similarity, toxicity predictions, and molecular docking. Three drugs were selected concerning their favorable interactions at the protein binding site and subsequently submitted to molecular dynamics simulations, which allowed us to elucidate their behavior and compare such theoretical results with experimental ones, obtained in biological assays also described in this paper.Communicated by Ramaswamy H. Sarma.

Chromatographic profile, in silico and in vivo study of the pharmacokinetic and toxicological properties of major constituent present in kefir, the kefiran.

Kefiran is a polysaccharide present in kefir grains that have been widely explored due to its potential health benefits. The objective of this work was to characterize and quantify the components present in the ethanolic extract of milk kefir grains; to study its pharmacokinetic and toxicological properties in silico and evaluate the acute toxicity of the kefiran in zebrafish. The prediction of pharmacokinetic properties was performed by QikProp software. In silico toxicity assessment was performed using the DEREK (deductive estimate of risk from existing knowledge) software. In the chromatographic, kefiran was identified as the major component. Results showed that the kefiran had low human oral absorption and intestinal absorption its due poor solubility profile; low logP value, indicating its lipophilicity and the low MDCK and Caco-2 cells permability, and unable to cross the blood-brain barrier. Kefiran did not present any structural warning for in silico toxicity. In zebrafish, the dose of 2,000 mg/kg of kefiran produced nonsignificant alterations in the analyzed organs. It can be said then that kefiran has an acceptable degree of safety for use in the development of drugs or functional foods. Further research such as in vivo testing to confirm its pharmacological potential is currently underway.

Synthesis, antitumor activity and in silico analyses of amino acid derivatives of artepillin C, drupanin and baccharin from green propolis.

Breast cancer has the highest incidence and mortality in females, while prostate cancer has the second-highest incidence in males. Studies have shown that compounds from Brazilian green propolis have antitumor activities and can selectively inhibit the AKR1C3 enzyme, overexpressed in hormone-dependent prostate and breast tumors. Thus, in an attempt to develop new cytotoxic inhibitors against these cancers, three prenylated compounds, artepillin C, drupanin and baccharin, were isolated from green propolis to synthesize new derivatives via coupling reactions with different amino acids. All obtained derivatives were submitted to antiproliferative assays against four cancer cells (MCF-7, MDA MB-231, PC-3, and DU145) and two normal cell lines (MCF-10A and PNT-2) to evaluate their cytotoxicity. In general, the best activity was observed for compound6e, derived from drupanin, which exhibited half-maximal inhibitory concentration (IC50) of 9.6 ± 3 µM and selectivity index (SI) of 5.5 against MCF-7 cells.In silicostudies demonstrated that these derivatives present coherent docking interactions and binding modes against AKR1C3, which might represent a possible mechanism of inhibition in MCF-7 cells.

Potential beneficial effects of kefir and its postbiotic, kefiran, on child food allergy.

Food allergies are known as the public health problem, affecting people of all age groups, but more commonly in babies and children, with consequences for nutritional status and quality of life. The increase in the consumption of healthy foods has consequently led to an increased demand for functional foods with specific health benefits. Thus, the pharmaceutical industry's interest in natural products has grown every time and is therefore considered as an alternative to synthetic drugs. Kefir has been outstanding for several years as promising in the manufacture of various pharmaceutical products, due to its nutritional and therapeutic properties for the treatment of many diseases. Currently, a wide variety of new functional foods are appearing on the market, representing an important segment. Postbiotics, for example, has stood out for being a product with action similar to probiotics, without offering side effects. The kefiran is the postbiotic from kefir that promotes potential beneficial effects on food allergy from the intestinal microbiome to the immune system. In this context, it is necessary to know the main promoting component of this functional effect. This review compiles the benefits that kefir, and especially its postbiotic, kefiran, can bring to food allergy. In addition, it serve as a subsidy for studies on the development of innovative nutraceutical products, including the use of kefiran as an alternative therapy in food allergy.

Virtual screening, ADME/Tox predictions and the drug repurposing concept for future use of old drugs against the COVID-19.

The new Coronavirus (SARS-CoV-2) is the cause of a serious infection in the respiratory tract called COVID-19. Structures of the main protease of SARS-CoV-2 (Mpro), responsible for the replication of the virus, have been solved and quickly made available, thus allowing the design of compounds that could interact with this protease and thus to prevent the progression of the disease by avoiding the viral peptide to be cleaved, so that smaller viral proteins can be released into the host's plasma. These structural data are extremely important for in silico design and development of compounds as well, being possible to quick and effectively identify potential inhibitors addressed to such enzyme's structure. Therefore, in order to identify potential inhibitors for Mpro, we used virtual screening approaches based with the structure of the enzyme and two compounds libraries, targeted to SARS-CoV-2, containing compounds with predicted activity against Mpro. In this way, we selected, through docking studies, the 100 top-ranked compounds, which followed to subsequent studies of pharmacokinetic and toxicity predictions. After all the simulations and predictions here performed, we obtained 10 top-ranked compounds that were again in silico analyzed inside the Mpro catalytic site, together some drugs that are being currently investigated for treatment of COVID-19. After proposing and analyzing the interaction modes of these compounds, we submitted one molecule then selected as template to a 2D similarity study in a database containing drugs approved by FDA and we have found and indicated Apixaban as a potential drug for future treatment of COVID-19.

The Aspergillus fumigatus Phosphoproteome Reveals Roles of High-Osmolarity Glycerol Mitogen-Activated Protein Kinases in Promoting Cell Wall Damage and Caspofungin Tolerance.

The filamentous fungus Aspergillus fumigatus can cause a distinct set of clinical disorders in humans. Invasive aspergillosis (IA) is the most common life-threatening fungal disease of immunocompromised humans. The mitogen-activated protein kinase (MAPK) signaling pathways are essential to the adaptation to the human host. Fungal cell survival is highly dependent on the organization, composition, and function of the cell wall. Here, an evaluation of the global A. fumigatus phosphoproteome under cell wall stress caused by the cell wall-damaging agent Congo red (CR) revealed 485 proteins potentially involved in the cell wall damage response. Comparative phosphoproteome analyses with the ΔsakA, ΔmpkC, and ΔsakA ΔmpkC mutant strains from the osmotic stress MAPK cascades identify their additional roles during the cell wall stress response. Our phosphoproteomics allowed the identification of novel kinases and transcription factors (TFs) involved in osmotic stress and in the cell wall integrity (CWI) pathway. Our global phosphoproteome network analysis showed an enrichment for protein kinases, RNA recognition motif domains, and the MAPK signaling pathway. In contrast to the wild-type strain, there is an overall decrease of differentially phosphorylated kinases and phosphatases in ΔsakA, ΔmpkC, and ΔsakA ΔmpkC mutants. We constructed phosphomutants for the phosphorylation sites of several proteins differentially phosphorylated in the wild-type and mutant strains. For all the phosphomutants, there is an increase in the sensitivity to cell wall-damaging agents and a reduction in the MpkA phosphorylation upon CR stress, suggesting these phosphosites could be important for the MpkA modulation and CWI pathway regulation.IMPORTANCEAspergillus fumigatus is an opportunistic human pathogen causing allergic reactions or systemic infections, such as invasive pulmonary aspergillosis in immunocompromised patients. The mitogen-activated protein kinase (MAPK) signaling pathways are essential for fungal adaptation to the human host. Fungal cell survival, fungicide tolerance, and virulence are highly dependent on the organization, composition, and function of the cell wall. Upon cell wall stress, MAPKs phosphorylate multiple target proteins involved in the remodeling of the cell wall. Here, we investigate the global phosphoproteome of the ΔsakA and ΔmpkCA. fumigatus and high-osmolarity glycerol (HOG) pathway MAPK mutants upon cell wall damage. This showed the involvement of the HOG pathway and identified novel protein kinases and transcription factors, which were confirmed by fungal genetics to be involved in promoting tolerance of cell wall damage. Our results provide understanding of how fungal signal transduction networks modulate the cell wall. This may also lead to the discovery of new fungicide drug targets to impact fungal cell wall function, fungicide tolerance, and virulence.

Allosteric Modulators of Potential Targets Related to Alzheimer's Disease: a Review.

Among neurodegenerative disorders, Alzheimer's disease (AD) is the most common type of dementia, and there is an urgent need to discover new and efficacious forms of treatment for it. Pathological patterns of AD include cholinergic dysfunction, increased ß-amyloid (Aß) peptide concentration, the appearance of neurofibrillary tangles, among others, all of which are strongly associated with specific biological targets. Interactions observed between these targets and potential drug candidates in AD most often occur by competitive mechanisms driven by orthosteric ligands that sometimes result in the production of side effects. In this context, the allosteric mechanism represents a key strategy; this can be regarded as the selective modulation of such targets by allosteric modulators in an advantageous manner, as this may decrease the likelihood of side effects. The purpose of this review is to present an overview of compounds that act as allosteric modulators of the main biological targets related to AD.

New Horizons on Molecular Pharmacology Applied to Drug Discovery: When Resonance Overcomes Radioligand Binding.

One of the cornerstones of rational drug development is the measurement of molecular parameters derived from ligand-receptor interaction, which guides therapeutic windows definition. Over the last decades, radioligand binding has provided valuable contributions in this field as key method for such purposes. However, its limitations spurred the development of more exquisite techniques for determining such parameters. For instance, safety risks related to radioactivity waste, expensive and controlled disposal of radioisotopes, radiotracer separation-dependence for affinity analysis, and one-site mathematical models-based fitting of data make radioligand binding a suboptimal approach in providing measures of actual affinity conformations from ligands and G proteincoupled receptors (GPCR). Current advances on high-throughput screening (HTS) assays have markedly extended the options of sparing sensitive ways for monitoring ligand affinity. The advent of the novel bioluminescent donor NanoLuc luciferase (Nluc), engineered from Oplophorus gracilirostris luciferase, allowed fitting bioluminescence resonance energy transfer (BRET) for monitoring ligand binding. Such novel approach named Nluc-based BRET (NanoBRET) binding assay consists of a real-time homogeneous proximity assay that overcomes radioligand binding limitations but ensures the quality in affinity measurements. Here, we cover the main advantages of NanoBRET protocol and the undesirable drawbacks of radioligand binding as molecular methods that span pharmacological toolbox applied to Drug Discovery. Also, we provide a novel perspective for the application of NanoBRET technology in affinity assays for multiple-state binding mechanisms involving oligomerization and/or functional biased selectivity. This new angle was proposed based on specific biophysical criteria required for the real-time homogeneity assigned to the proximity NanoBRET protocol.

3D descriptors calculation and conformational search to investigate potential bioactive conformations, with application in 3D-QSAR and virtual screening in drug design.

The knowledge of the bioactive conformation for an active hit is relevant because of the easier interpretation and the general quality of the recognition models of protein and ligand. With the aim of investigating potential bioactive conformations without previous structural knowledge of the molecular target, we present herewith a 'protocol' that could be used which includes generation of low-energy conformations, calculations of tridimensional descriptors and investigation of structural similarity via principal component analysis. The protocol was used in the search for potential bioactive conformations. An initial selection of targets was made from a set of protein-ligand complexes with structure deposited in the Protein Data Bank, which was systematically filtered by lead-like rules, resulting in 45 ligands of 8 important therapeutic targets. After extensive optimization of the protocol and parameters of both OMEGA and Pentacle softwares, the best results were obtained for series of compounds such as the beta-trypsin and urokinase inhibitors, which are more structurally related among each other, inside the respective therapeutic class. Future improvements of the protocol, including a suitable choice and combination of robust 3D descriptors, could yield more reliable and less restrictive results, with general and diverse applications in drug design, in particular for improving the 3D-QSAR methodologies as well as virtual screening experiments for a more reliable selection of new lead compounds for different molecular targets.

In silico design and search for acetylcholinesterase inhibitors in Alzheimer's disease with a suitable pharmacokinetic profile and low toxicity.

Alzheimer's disease is a complex neurodegenerative disorder of the central nervous system, characterized by amyloid-ß deposits, τ-protein aggregation, oxidative stress and reduced levels of acetylcholine in the brain. One pharmacological approach is to restore acetylcholine level by inhibiting acetylcholinesterase (AChE) with reversible inhibitors, such as galanthamine, thus helping to improve the cognitive symptoms of the disease. In order to design new galanthamine derivatives and search for novel, potential inhibitors with improved interactions, as well as a suitable pharmacokinetic profile and low toxicity, several molecular modeling techniques were applied. These techniques included the investigation of AChE-drug complexes (1QT1 and 1ACJ Protein Data Bank codes), ligand-binding sites calculation within the active site of the enzyme, pharmacophore perception of galanthamine derivatives, virtual screening, toxicophorical analysis and estimation of pharmacokinetics properties. A total of four galanthamine derivatives having a N-alkyl-phenyl chain were designed, since the tertiary amine substituents could reach the peripheral anionic site that is not occupied by galanthamine. In addition, 12 drug-like compounds from the Ilibdiverse database were selected by virtual screening as novel, hypothetical AChE inhibitors. The toxicophorical analysis revealed that only four proposed inhibitors have chemical groups able to develop mutagenicity and chromosome damage. The remaining compounds showed only mild or none toxicophorical alerts. At least three screened compounds presented theoric parameters consistent with good oral bioavailability. The designed molecules have the potential to become new lead compounds that might guide the design of drugs with optimized pharmacodynamic and pharmacokinetic properties in order to improve the treatment of Alzheimer's disease by creating new pharmacotherapeutic options.

Toxicophoric and metabolic in silico evaluation of benzimidazole and phenylbenzamide derivatives with potential application as anticancer agents.

Poor pharmacokinetics and toxicity are responsible for most drug candidate failures. In order to attempt to some degree of ADMET (Absorption, Distribution, Metabolism, Excrection and Toxicity) information, in silico predictions arise currently as an interesting alternative to evaluate prototypes during early stages of the drug design processes, especially for anticancer candidates that constitute a class of therapeutic agents that exhibit substantial toxicity. A benzimidazole and a phenylbenzamide derivatives, previously identified as novel anticancer lead compounds able to prevent DNA binding to hnRNP K protein, were evaluated in silico regarding their metabolic profile and toxicity potential in order to give insights to the design of drug candidates with an adequate pharmaceutical profile. Considering the structure of proposed metabolites for both molecules, the phenylbenzamide derivative seems to be a molecule with better pharmaceutic profile, since its possible metabolites present a milder degree of chemical structure toxic alerts than the benzimidazole derivative that can cause chromosome damage induced by the benzimidazole group. It would be desirable during optimization of the phenylbenzamide derivative to maintain these characteristics during generation of analogues with substituents that are not known as potent toxicophoric groups. For the benzimidazole derivative, if the toxic events are really severe as it seems, one possible strategy would be replace the benzimidazole ring system by bioisosteres with lower toxic potential, hoping to maintain or enhance biological activity.

Computer-aided drug design and ADMET predictions for identification and evaluation of novel potential farnesyltransferase inhibitors in cancer therapy.

We have used various computational methodologies including molecular dynamics, density functional theory, virtual screening, ADMET predictions and molecular interaction field studies to design and analyze four novel potential inhibitors of farnesyltransferase (FTase). Evaluation of two proposals regarding their drug potential as well as lead compounds have indicated them as novel promising FTase inhibitors, with theoretically interesting pharmacotherapeutic profiles, when compared to the very active and most cited FTase inhibitors that have activity data reported, which are launched drugs or compounds in clinical tests. One of our two proposals appears to be a more promising drug candidate and FTase inhibitor, but both derivative molecules indicate potentially very good pharmacotherapeutic profiles in comparison with Tipifarnib and Lonafarnib, two reference pharmaceuticals. Two other proposals have been selected with virtual screening approaches and investigated by us, which suggest novel and alternatives scaffolds to design future potential FTase inhibitors. Such compounds can be explored as promising molecules to initiate a research protocol in order to discover novel anticancer drug candidates targeting farnesyltransferase, in the fight against cancer.

Psoralen and bergapten: in silico metabolism and toxicophoric analysis of drugs used to treat vitiligo.

PURPOSE: To discuss the contribution of psoralen and bergapten metabolites on psoralens toxicity. METHODS: Computational chemistry prediction of metabolic reactions and toxicophoric groups based on the expert systems Derek and Meteor. RESULTS: a total of 15 metabolites were suggested for both psoralen and bergapten based on phase 1 and 2 biotransformations until the 3rd generation. Five toxicophoric substructures were shared among psoralen, bergapten and their corresponding metabolites; one toxicophoric marker (resorcinol) was only identified in bergapten and its biotransformation products. CONCLUSION: Although the toxic effects of psoralens are well known and documented, there is little information concerning the role of their metabolites in this process. We believe this work add to the knowledge of which molecular substructures are relevant to the process of metabolism and toxicity induction, thus guiding the search and development of more effective and less toxic drugs to treat vitiligo.

Molecular dynamics, density functional, ADMET predictions, virtual screening, and molecular interaction field studies for identification and evaluation of novel potential CDK2 inhibitors in cancer therapy.

In this work, we have used molecular dynamics, density functional theory, virtual screening, ADMET predictions, and molecular interaction field studies to design and propose eight novel potential inhibitors of CDK2. The eight molecules proposed showed interesting structural characteristics that are required for inhibiting the CDK2 activity and show potential as drug candidates for the treatment of cancer. The parameters related to the Rule of Five were calculated, and only one of the molecules violated more than one parameter. One of the proposals and one of the drug-like compounds selected by virtual screening indicated to be promising candidates for CDK2-based cancer therapy.

Current topics in computer-aided drug design.

The addition of computer-aided drug design (CADD) technologies to the research and drug discovery approaches could lead to a reduction of up to 50% in the cost of drug design. Designing a drug is the process of finding or creating a molecule which has a specific activity on a biological organism. Development and drug discovery is a time-consuming, expensive, and interdisciplinary process whereas scientific advancements during the past two decades have altered the way pharmaceutical research produces new bioactive molecules. Advances in computational techniques and hardware solutions have enabled in silico methods to speed up lead optimization and identification. We will review current topics in computer-aided molecular design underscoring some of the most recent approaches and interdisciplinary processes. We will discuss some of the most efficient pathways and design.

Transdermal delivery of ketoprofen: the influence of drug-dioleylphosphatidylcholine interactions.

PURPOSE: Considering that most inflammatory diseases occur locally and near the body surface, transdermal delivery of non-steroidal anti-inflammatory drugs (NSAIDs) may be an interesting strategy for delivering these drugs directly to the diseased site. To optimize ketoprofen (KP) transdermal delivery we investigated the influence of dioleylphosphatidylcholine (DOPC) on skin permeation. MATERIALS AND METHODS: The formulations studied were: i) a physical mixture of KP and DOPC and ii) DOPC and KP complex, in a molar ratio of 1:3, obtained by dissolution of the components in chloroform followed by drying under a N2 atmosphere. Both systems were dispersed in mineral oil and the in vitro percutaneous was assayed by absorption using a flow through diffusion cell. Differential Scanning Calorimetry (DSC) and 1H NMR studies were carried out to characterize KP and DOPC interactions. Geometry optimizations using Density Functional Theory and semiempirical methods, as well as a flexible docking procedure were carried out to obtain a binding model for KP with DOPC. KP solubility and partition studies in the formulations, as well as skin irritation and hypersensitivity assays were also carried out. RESULTS: DSC determinations in the complex showed enthalpy and temperature depressions, indicating KP and DOPC interaction. In addition, dipole-dipole interactions between the KP carboxylic acid and OH groups in phospholipids were shown by 1H NMR studies. Based on the NMR studies, a KP-DOPC binding model is proposed, in which KP is involved by the two long aliphatic chains of the phospholipid. Solubility studies indicated that DOPC improved drug solubility. KP permeation was enhanced by both formulations tested, but the complex also increased its skin uptake. Such behavior could be attributed to the solubilizing, melting and enhancing effects of DOPC. Skin irritation and hypersensitivity were not significantly changed compared to control, suggesting that the formulation may be therapeutically explored for KP transdermal delivery.

Fluorescence properties of tryptophan residues in the monomeric d-chain of Glossoscolex paulistus hemoglobin: an interpretation based on a comparative molecular model.

The primary structure of the 142 residue Glossoscolex paulistus d-chain hemoglobin has been determined from Edman degradation data of 11 endo-Glu-C peptides and 11 endo-Lys-C peptides, plus the results of Edman degradation of the intact globin. Tryptophan occupies positions 15, 33 and 129. Homology modeling allowed us to assign the positions of these Trp residues relative to the heme and its environment. The reference coordinates of the indole rings (average coordinates of the C(varepsilon2) and C(delta2) atoms) for W15 and W129 were 16.8 and 18.5 A, respectively, from the geometric center of the heme, and W33 was located in close proximity to the heme group at a distance which was approximately half of that for W15 and W129. It was possible to identify three rotamers of W33 on the basis of electrostatic and Van der Waals energy criteria. The calculated distances from the center of the heme were 8.3, 8.4 and 9.1 A for Rot1, Rot2 and Rot3, respectively. Radiationless energy transfer from the excited indole to the heme was calculated on the basis of Förster theory. For W33, the distance was more important than the orientation factor, kappa(2), due to its proximity to the heme. However, based on kappa(2), Rot2 (kappa(2)=0.945) was more favorable for the energy transfer than Rot1 (kappa(2)=0.433) or Rot3 (kappa(2)=0.125). In contrast, despite its greater distance from the heme, the kappa(2) of W129 (2.903) established it as a candidate to be more efficiently quenched by the heme than W15 (kappa(2)=0.191). Although the Förster approach is powerful for the evaluation of the relative efficiency of quenching, it can only explain pico- and sub-nanosecond lifetimes. With the average lifetime, =3 ns, measured for the apomonomer as the reference, the lifetimes calculated for each emitter were: W33-1 (1 ps), W33-2 (2 ps), W33-3 (18 ps), W129 (100 ps), and W15 (600 ps). Experimentally, there are four components for oxymonomers at pH 7: two long ones of 4.6 and 2.1 ns, which contribute approximately 90% of the total fluorescence, one of 300 ps (4%), and the last one of 33 ps (7.4%). It is clear that the equilibrium structure resulting from homology modeling explains the sub-nanosecond fluorescence lifetimes, while the nanosecond range lifetimes require more information about the protein in solution, since there is a significant contribution of lifetimes that resemble the apo molecule.