DrugCentral 2023 extends human clinical data and integrates veterinary drugs.

DrugCentral monitors new drug approvals and standardizes drug information. The current update contains 285 drugs (131 for human use). New additions include: (i) the integration of veterinary drugs (154 for animal use only), (ii) the addition of 66 documented off-label uses and iii) the identification of adverse drug events from pharmacovigilance data for pediatric and geriatric patients. Additional enhancements include chemical substructure searching using SMILES and 'Target Cards' based on UniProt accession codes. Statistics of interests include the following: (i) 60% of the covered drugs are on-market drugs with expired patent and exclusivity coverage, 17% are off-market, and 23% are on-market drugs with active patents and exclusivity coverage; (ii) 59% of the drugs are oral, 33% are parenteral and 18% topical, at the level of the active ingredients; (iii) only 3% of all drugs are for animal use only; however, 61% of the veterinary drugs are also approved for human use; (iv) dogs, cats and horses are by far the most represented target species for veterinary drugs; (v) the physicochemical property profile of animal drugs is very similar to that of human drugs. Use cases include azaperone, the only sedative approved for swine, and ruxolitinib, a Janus kinase inhibitor.

Exploring DrugCentral: from molecular structures to clinical effects.

DrugCentral, accessible at https://drugcentral.org , is an open-access online drug information repository. It covers over 4950 drugs, incorporating structural, physicochemical, and pharmacological details to support drug discovery, development, and repositioning. With around 20,000 bioactivity data points, manual curation enhances information from several major digital sources. Approximately 724 mechanism-of-action (MoA) targets offer updated drug target insights. The platform captures clinical data: over 14,300 on- and off-label uses, 27,000 contraindications, and around 340,000 adverse drug events from pharmacovigilance reports. DrugCentral encompasses information from molecular structures to marketed formulations, providing a comprehensive pharmaceutical reference. Users can easily navigate basic drug information and key features, making DrugCentral a versatile, unique resource. Furthermore, we present a use-case example where we utilize experimentally determined data from DrugCentral to support drug repurposing. A minimum activity threshold t should be considered against novel targets to repurpose a drug. Analyzing 1156 bioactivities for human MoA targets suggests a general threshold of 1 µM: t = 6 when expressed as - log[Activity(M)]). This applies to 87% of the drugs. Moreover, t can be refined empirically based on water solubility (S): t = 3 - logS, for logS < - 3. Alongside the drug repurposing classification scheme, which considers intellectual property rights, market exclusivity protections, and market accessibility, DrugCentral provides valuable data to prioritize candidates for drug repurposing programs efficiently.

DrugCentral 2021 supports drug discovery and repositioning.

DrugCentral is a public resource (http://drugcentral.org) that serves the scientific community by providing up-to-date drug information, as described in previous papers. The current release includes 109 newly approved (October 2018 through March 2020) active pharmaceutical ingredients in the US, Europe, Japan and other countries; and two molecular entities (e.g. mefuparib) of interest for COVID19. New additions include a set of pharmacokinetic properties for ∼1000 drugs, and a sex-based separation of side effects, processed from FAERS (FDA Adverse Event Reporting System); as well as a drug repositioning prioritization scheme based on the market availability and intellectual property rights forFDA approved drugs. In the context of the COVID19 pandemic, we also incorporated REDIAL-2020, a machine learning platform that estimates anti-SARS-CoV-2 activities, as well as the 'drugs in news' feature offers a brief enumeration of the most interesting drugs at the present moment. The full database dump and data files are available for download from the DrugCentral web portal.

Off-Patent Drug Repositioning.

Drug repositioning aims to reuse "old" drugs to treat diseases outside their approved indication(s). Composition-of-matter patents and FDA exclusivities can hinder the immediate availability of some drugs to be repositioned (repurposed). Here, we analyze data from the FDA Orange Book and use current on-market patent validity and exclusivities to classify drugs into on-patent (ONP), off-patent (OFP), and off-market (OFM) sets. In the absence of an unanimously accepted definition for small molecules, these sets include organic molecules and peptides with molecular weight between 100 and 1250, which resulted in 237 ONP drugs, 320 OFM, and 996 OFP drugs, respectively. We discuss the differences between the three categories in terms of primary molecular properties, chemical diversity, mechanism-of-action target classes, and therapeutic areas and comment on the enrichment of OFP drugs in the near future. Given the intellectual property landscape, and in the absence of specific property rights, we suggest that drugs should be prioritized as follows, to improve the repositioning strategy: (i) OFP, (ii) OFM, and (iii) ONP, respectively.

Enhancing Molecular Promiscuity Evaluation Through Assay Profiles.

PURPOSE: The growing amount of heterogeneous bioactivity data requires effective strategies to assess the promiscuity/selectivity of small-molecules and aid drug discovery. In the current study, we aim to evaluate the potential of assay profiles (APs, i.e., unique combinations of assay-related features describing how activity determinations were performed and reported) in molecular promiscuity analysis. METHODS: Using PubChem bioactivity data, we computed for all Molecular Libraries Small Molecule Repository (MLSMR library) compounds the frequency of hits score (FoH, i.e., the ratio between the number of times the compound was found active and the number of times it was tested), which were subsequently fit into 32 theoretical APs. The promiscuity of drugs and non-drugs was compared at different levels of test results. RESULTS: We found 8 dominant APs, indicating that compounds tested in more than ten assays (or against ten targets) and found active at least once tend to reach near to maximum hit rates in scientific literature and confirmatory assays (e.g., 95% of the drugs show FoH scores >0.93). Primary and high-throughput screening testing results in very low hit rates (e.g., 95% of the compounds show FoH scores <0.11), promoting a different perspective of promiscuity. In general, drugs exert higher promiscuity compared to non-drugs. Targets and classes of drugs are also discussed within the main APs. CONCLUSION: APs contain relevant features and are suited for big data promiscuity analysis. The activity data of the main APs are freely available on www.chembioinf.ro .

Modeling Kinase Inhibition Using Highly Confident Data Sets.

Protein kinases form a consistent class of promising drug targets, and several efforts have been made to predict the activities of small molecules against a representative part of the kinome. This study continues our previous work ( Bora , A. ; Avram , S. ; Ciucanu , I. ; Raica , M. ; Avram , S. Predictive Models for Fast and Effective Profiling of Kinase Inhibitors . J. Chem. Inf. MODEL: 2016 , 56 , 895 - 905 ; www.chembioinf.ro ) aiming to build and measure the performance of ligand-based kinase inhibitor prediction models. Here we analyzed kinase-inhibitor pairs with multiple activity points extracted from the ChEMBL database and identified the main sources of inconsistency. Our results indicate that lower IC50 values are usually less affected by errors and reflect more accurately the structure-activity relationship of the molecules against the target, ideally for quantitative structure-activity relationship studies. Further, we modeled the activities of 104 kinases using unbiased target-specific activity points. The performance of predictors built on extended connectivity fingerprints (ECFP4) and two-dimensional pharmacophore fingerprints (PFPs) are compared by means of tolerance intervals (TIs) (95%/95%) in virtual screening (VS) and classification tasks using external random ( RandSets) and diversity-based ( DivSets) test sets. We found that the two encodings perform superior to each other on different kinases in VS and that PFP models perform consistently better in classifying actives (higher sensitivity). Next, we combined the two encodings into a single one (PFPECFP) and demonstrated that especially in VS (as indicated by the exponential receiver operating curve enrichment metric (eROCE)), for the vast majority of kinases the model performance increased compared with the individual fingerprint models. These findings are highlighted in the more challenging DivSets compared with RandSets. The current paper explores the boundaries of inhibitor predictors for individual kinases to enhance VS and ultimately aid the discovery of novel compounds with desirable polypharmacology.

Pharmacophore-based screening and drug repurposing exemplified on glycogen synthase kinase-3 inhibitors.

The current study was conducted to elaborate a novel pharmacophore model to accurately map selective glycogen synthase kinase-3 (GSK-3) inhibitors, and perform virtual screening and drug repurposing. Pharmacophore modeling was developed using PHASE on a data set of 203 maleimides. Two benchmarking validation data sets with focus on selectivity were assembled using ChEMBL and PubChem GSK-3 confirmatory assays. A drug repurposing experiment linking pharmacophore matching with drug information originating from multiple data sources was performed. A five-point pharmacophore model was built consisting of a hydrogen bond acceptor (A), hydrogen bond donor (D), hydrophobic (H), and two rings (RR). An atom-based 3D quantitative structure-activity relationship (QSAR) model showed good correlative and satisfactory predictive abilities (training set [Formula: see text]; test set: [Formula: see text]; whole data set: stability [Formula: see text]). Virtual screening experiments revealed that selective GSK-3 inhibitors are ranked preferentially by Hypo-1, but fail to retrieve nonselective compounds. The pharmacophore and 3D QSAR models can provide assistance to design novel, potential GSK-3 inhibitors with high potency and selectivity pattern, with potential application for the treatment of GSK-3-driven diseases. A class of purine nucleoside antileukemic drugs was identified as potential inhibitor of GSK-3, suggesting the reassessment of the target range of these drugs.

Predictive Models for Fast and Effective Profiling of Kinase Inhibitors.

In this study we developed two-dimensional pharmacophore-based random forest models for the effective profiling of kinase inhibitors. One hundred seven prediction models were developed to address distinct kinases spanning over all kinase groups. Rigorous external validation demonstrates excellent virtual screening and classification potential of the predictors and, more importantly, the capacity to prioritize novel chemical scaffolds in large chemical libraries. The models built upon more diverse and more potent compounds tend to exert the highest predictive power. The analysis of ColBioS-FlavRC (Collection of Bioselective Flavonoids and Related Compounds) highlighted several potentially promiscuous derivatives with undesirable selectivity against kinases. The prediction models can be downloaded from www.chembioinf.ro .

ColBioS-FlavRC: a collection of bioselective flavonoids and related compounds filtered from high-throughput screening outcomes.

Flavonoids, the vastest class of natural polyphenols, are extensively investigated for their multiple benefits on human health. Due to their physicochemical or biological properties, many representatives are considered to exhibit low selectivity among various protein targets or to plague high-throughput screening (HTS) outcomes. The aim of this study is to highlight reliable, bioselective compounds sharing flavonoidic scaffolds in HTS experiments. A filtering scheme was applied to remove undesired flavonoids (and related compounds) from confirmatory PubChem bioassays. A number of 433 compounds addressing various protein targets form the core of the collection of bioselective flavonoids and related compounds (ColBioS-FlavRC). With an additional set of 2908 inactive related compounds, ColBioS-FlavRC offers the grounds for method optimization and validation. We exemplified the use of ColBioS-FlavRC by pharmacophore modeling, subsequently (externally) validated for virtual screening purposes. The early enrichment capabilities of the pharmacophore hypotheses were measured by means of the median exponential retriever operating curve enrichment (MeROCE), a suited metric in comparative evaluations of virtual screening methods. ColBioS-FlavRC is available in the Supporting Information and is freely accessible for further studies.

Exploring the biological promiscuity of high-throughput screening hits through DFT calculations.

The goal of this study is the understanding of biologically promiscuous compounds (frequent hitters) in HTS outcomes through their chemical behavior estimated via reactivity descriptors. Chemical reactivity is often an undesirable property due to the lack in biological selectivity of compounds comprised in HTS libraries. In this study the reactivity indexes have been computed within the DFT formalism, at different levels of theory, for two classes of representative compounds compiled from PubChem database, one comprising frequent hitters and the second one comprising rare hitters (biologically more selective compounds). We found that frequent hitters exert increased reactivity, mainly due to their electrophilic character, compared to the more selective class of compounds.

PLS and shape-based similarity analysis of maleimides--GSK-3 inhibitors.

CONTEXT: Glycogen synthase kinase-3 (GSK-3) overactivity was correlated with several pathologies including type 2 diabetes mellitus, Alzheimer's disease, cancer, inflammation, obesity, etc. OBJECTIVE: The aim of the current investigation was to model the inhibitory activity of maleimide derivatives--inhibitors of GSK-3, to evaluate the impact of alignment on statistical performances of the Quantitative Structure-Activity Relationship (QSAR) and the effect of the template on shape-similarity--binding affinity relationship. MATERIALS AND METHODS: Dragon descriptors were used to generate Projection to Latent Structures (PLS) models in order to identify the structural prerequisites of maleimides to inhibit GSK-3. Additionally, shape/volume structural analysis of binding site interactions was evaluated. RESULTS: Reliable statistics R(2)(Y(CUM)) = 0.938/0.920, Q((2)(Y)(CUM)) = 0.866/0.838 for aligned and alignment free QSAR models and significant (Pearson, Kendall and Spearman) correlations between shape/volume similarity and affinities were obtained. DISCUSSION AND CONCLUSIONS: The crucial structural features modulating the activity of maleimides include topology, charge, geometry, 2D autocorrelations, 3D-MoRSE as well as shape/volume and molecular flexibility.

Overview of the Knowledge Management Center for Illuminating the Druggable Genome.

The Knowledge Management Center (KMC) for the Illuminating the Druggable Genome (IDG) project aims to aggregate, update, and articulate protein-centric data knowledge for the entire human proteome, with emphasis on the understudied proteins from the three IDG protein families. KMC collates and analyzes data from over 70 resources to compile the Target Central Resource Database (TCRD), which is the web-based informatics platform (Pharos). These data include experimental, computational, and text-mined information on protein structures, compound interactions, and disease and phenotype associations. Based on this knowledge, proteins are classified into different Target Development Levels (TDLs) for identification of understudied targets. Additional work by the KMC focuses on enriching target knowledge and producing DrugCentral and other data visualization tools for expanding investigation of understudied targets.

Retrospective group fusion similarity search based on eROCE evaluation metric.

In this study, a simple evaluation metric, denoted as eROCE was proposed to measure the early enrichment of predictive methods. We demonstrated the superior robustness of eROCE compared to other known metrics throughout several active to inactive ratios ranging from 1:10 to 1:1000. Group fusion similarity search was investigated by varying 16 similarity coefficients, five molecular representations (binary and non-binary) and two group fusion rules using two reference structure set sizes. We used a dataset of 3478 actives and 43,938 inactive molecules and the enrichment was analyzed by means of eROCE. This retrospective study provides optimal similarity search parameters in the case of ALDH1A1 inhibitors.

GIS utilization in determining the limits between low morphological units.

The limit between low morphological units is quite difficult to be established due to the fact that the lithological, elevation and/or structure differences are hardly visible. Thus, the contact between piedmont and plain units is represented, in many cases, by areas where limits are hard to be drawn. In order to establish the way geographic information system (GIS) may contribute to the determination of the aforementioned limits, we used the situation of the limit between the Getic Piedmont and the Romanian plain and compared the obtained results to examples from other areas with different formation and evolution conditions. Presently, the Getic Piedmont undergoes a modelling process as valleys widen, torrential enlargements laterally advance, and slope denudation develop, which make difficult to establish a clear southern limit. The fan texture of the valleys, their age, the flowing direction, as well as the differences of altitude, lithology, and structure disposal represent the main elements when establishing limits even in the case of GIS analysis. The drawing of thematic layers rendering elements such as hydrographic system, hypsometry, DEM on a traditional cartographic support, orthophotoplans, satellite images, for a better interpretation of the geomorphologic aspects allow a more objective attempt of reconstructing the morphological limits.

PLS-DA - Docking Optimized Combined Energetic Terms (PLSDA-DOCET) protocol: a brief evaluation.

Docking studies have become popular approaches in drug design, where the binding energy of the ligand in the active site of the protein is estimated by a scoring function. Many promising techniques were developed to enhance the performance of scoring functions including the fusion of multiple scoring functions outcomes into a so-called consensus scoring function. Hereby, we evaluated the target oriented consensus technique using the energetic terms of several scoring functions. The approach was denoted PLSDA-DOCET. Optimization strategies for consensus energetic terms and scoring functions based on ROC metric were compared to classical rigid docking and to ligand-based similarity search methods comprising 2D fingerprints and ROCS. The ROCS results indicate large performance variations depending on the biological target. The AUC-based strategy of PLSDA-DOCET outperformed the other docking approaches regarding simple retrieval and scaffold-hopping. The superior performance of PLSDA-DOCET protocol relative to single and combined scoring functions was validated on an external test set. We found a relative low mean correlation of the ranks of the chemotypes retrieved by the PLSDA-DOCET protocol and all the other methods employed here.

Comprehensive investigation of selectivity landscape of glycogen synthase kinase-3 inhibitors.

Interaction signatures of drug candidates are characteristic to off-target (neutral) and antitarget (negative) effects, inferring reduced efficiency, side-effects and high attrition rate. Today's retroactive scaled-down virtual screening (VS) experiments relying on benchmarking datasets are extensively involved to assess ligand enrichment in the real-world problem. In recent years, unbiased benchmarking sets turned into a tremendous need to assist virtual screening methodologies for emerging drug targets. To date, the benchmarking datasets are quite limited, whereas glycogen synthase kinase-3 (GSK-3) is not included into directories of benchmarking datasets such as DUD-e, MUV, etc. Herein we introduced our in-house algorithm to build an unbiased benchmarking dataset, including highly selective, moderately selective and nonselective inhibitors for a significant therapeutic target - GSK-3, suitable for both ligand-based and structure-based VS approaches. These datasets are unbiased in terms of physico-chemical properties and topological descriptors, as resulted from mean(ROC-AUC) leave-one-out cross-validation (LOO CV). and additional 2 D similarity search. Moreover, we investigated the gradual selectivity dataset by application of multiple 2 D similarity coefficients and distances, 3 D similarity and docking. Besides the resulted links between the enrichment of selective GSK-3 inhibitors and their chemical structures, a database of compounds and their 3 D similarity signatures including cut-off thresholds for enhanced selectivity was generated. 2 D similarity space analysis revealed that selectivity problem cannot be evaluated appropriately with 2 D similarity searching alone. The current analysis provided useful, comprehensive insights, which may facilitate the knowledge-based identification of novel selective GSK-3 inhibitors.Communicated by Ramaswamy H. Sarma.

Applying a Complex Integrated Method for Mapping and Assessment of the Degraded Ecosystem Hotspots from Romania.

To meet the global challenges of climate change and human activity pressure on biodiversity conservation, it has become vital to map such pressure hotspots. Large areas, such as nation-wide regions, are difficult to map from the point of view of the resources needed for such mapping (human resources, hard and soft resources). European biodiversity policies have focused on restoring degraded ecosystems by at least 10% by 2020, and new policies aim to restore up to 30% of degraded ecosystems by 2030. In this study, methods developed and applied for the assessment of the degradation state of the ecosystems in a semi-automatic manner for the entire Romanian territory (238,391 km2) are presented. The following ecosystems were analyzed: forestry, grassland, rivers, lakes, caves and coastal areas. The information and data covering all the ecoregions of the Romania (~110,000 km2) were analyzed and processed, based on GIS and remote sensing techniques. The largest degraded areas were identified within the coastal area (49.80%), grassland ecosystems (38.59%) and the cave ecosystems (2.66%), while 27.64% of rivers ecosystems were degraded, followed by 8.52% of forest ecosystems, and 14.05% of lakes ecosystems. This analysis can contribute to better definition of the locations of the most affected areas, which will yield a useful spatial representation for future ecological reconstruction strategy.

Evolution of Romania's Economic Structure and Environment Degradation - An Assessment Through LMDI Decomposition Approach.

PURPOSE: This paper studies the relationships between air pollutants (PM10, PM2.5, N2O) and different diseases (tumors, skin and respiratory) and the factors influencing air pollutant emissions in Romania. METHODS: The methods are Toda-Yamamoto procedure of non-causality Granger test, grey relational analysis and logarithmic mean Divisia index method (LMDI). RESULTS: Air pollutants intensities dropped significantly over 2008-2017 period due to structural changes. The only economic activity that showed an increase both in volume and intensity of air pollutants, despite a downward trend of farming activities output is agriculture. Technology improvements play a significant role in mitigation of PM2.5 emissions and a moderate role in mitigation of PM10 emissions. For N2O emissions technology used contributed to an increase of N2O intensities. CONCLUSION: Health policy makers should address the issue of technology improvements and mitigation of agriculture emissions to improve health of individuals and air quality.

Evaluation of the quality of lentic ecosystems in Romania by a GIS based WRASTIC model.

Globally, ecosystems are constantly degrading as a result of pressures derived from human activities and climate change. For working towards the restoration of the natural balance, it is necessary to evaluate the deviations induced in the ecosystems, to identify where the changes took place, to know what is their amplitude and to decide where it is possible to get involved. Many aquatic ecosystems are depreciated and their restoration is often difficult. Development of appropriate assessment methodologies will improve the decision-making process in public policies for environmental protection and conservation of biodiversity. This study presents an assessment of the degradation level of lentic ecosystems in Romania, performed through a multi-criteria analysis. An extension of the WRASTIC index (Wastewater-Recreational-Agricultural-Size-Transportations-Indutrial-Cover) was generated, namely WRASTIC-HI. The new index was obtained by including values derived from the Potential Pollutant Load index. The analysis showed that 13% of the evaluated lakes are natural, 56.5% are semi-degraded and 30.5% are degraded. The proposed methodology allows to determine the spatial distribution of the degradation sources and to calculate the corresponding indicators. The results obtained provide a useful tool for diagnostic step that can be used as a cornerstone to further identification of environmental conflicts and proposals for improvement of the ecological status of the lentic ecosystems.