Advancing Drug Safety in Drug Development: Bridging Computational Predictions for Enhanced Toxicity Prediction.

The attrition rate of drugs in clinical trials is generally quite high, with estimates suggesting that approximately 90% of drugs fail to make it through the process. The identification of unexpected toxicity issues during preclinical stages is a significant factor contributing to this high rate of failure. These issues can have a major impact on the success of a drug and must be carefully considered throughout the development process. These late-stage rejections or withdrawals of drug candidates significantly increase the costs associated with drug development, particularly when toxicity is detected during clinical trials or after market release. Understanding drug-biological target interactions is essential for evaluating compound toxicity and safety, as well as predicting therapeutic effects and potential off-target effects that could lead to toxicity. This will enable scientists to predict and assess the safety profiles of drug candidates more accurately. Evaluation of toxicity and safety is a critical aspect of drug development, and biomolecules, particularly proteins, play vital roles in complex biological networks and often serve as targets for various chemicals. Therefore, a better understanding of these interactions is crucial for the advancement of drug development. The development of computational methods for evaluating protein-ligand interactions and predicting toxicity is emerging as a promising approach that adheres to the 3Rs principles (replace, reduce, and refine) and has garnered significant attention in recent years. In this review, we present a thorough examination of the latest breakthroughs in drug toxicity prediction, highlighting the significance of drug-target binding affinity in anticipating and mitigating possible adverse effects. In doing so, we aim to contribute to the development of more effective and secure drugs.

POSEIDON: Peptidic Objects SEquence-based Interaction with cellular DOmaiNs: a new database and predictor.

Cell-penetrating peptides (CPPs) are short chains of amino acids that have shown remarkable potential to cross the cell membrane and deliver coupled therapeutic cargoes into cells. Designing and testing different CPPs to target specific cells or tissues is crucial to ensure high delivery efficiency and reduced toxicity. However, in vivo/in vitro testing of various CPPs can be both time-consuming and costly, which has led to interest in computational methodologies, such as Machine Learning (ML) approaches, as faster and cheaper methods for CPP design and uptake prediction. However, most ML models developed to date focus on classification rather than regression techniques, because of the lack of informative quantitative uptake values. To address these challenges, we developed POSEIDON, an open-access and up-to-date curated database that provides experimental quantitative uptake values for over 2,300 entries and physicochemical properties of 1,315 peptides. POSEIDON also offers physicochemical properties, such as cell line, cargo, and sequence, among others. By leveraging this database along with cell line genomic features, we processed a dataset of over 1,200 entries to develop an ML regression CPP uptake predictor. Our results demonstrated that POSEIDON accurately predicted peptide cell line uptake, achieving a Pearson correlation of 0.87, Spearman correlation of 0.88, and r2 score of 0.76, on an independent test set. With its comprehensive and novel dataset, along with its potent predictive capabilities, the POSEIDON database and its associated ML predictor signify a significant leap forward in CPP research and development. The POSEIDON database and ML Predictor are available for free and with a user-friendly interface at https://moreiralab.com/resources/poseidon/ , making them valuable resources for advancing research on CPP-related topics. Scientific Contribution Statement: Our research addresses the critical need for more efficient and cost-effective methodologies in Cell-Penetrating Peptide (CPP) research. We introduced POSEIDON, a comprehensive and freely accessible database that delivers quantitative uptake values for over 2,300 entries, along with detailed physicochemical profiles for 1,315 peptides. Recognizing the limitations of current Machine Learning (ML) models for CPP design, our work leveraged the rich dataset provided by POSEIDON to develop a highly accurate ML regression model for predicting CPP uptake.

DELFOS-drug efficacy leveraging forked and specialized networks-benchmarking scRNA-seq data in multi-omics-based prediction of cancer sensitivity.

MOTIVATION: Cancer is currently one of the most notorious diseases, with over 1 million deaths in the European Union alone in 2022. As each tumor can be composed of diverse cell types with distinct genotypes, cancer cells can acquire resistance to different compounds. Moreover, anticancer drugs can display severe side effects, compromising patient well-being. Therefore, novel strategies for identifying the optimal set of compounds to treat each tumor have become an important research topic in recent decades. RESULTS: To address this challenge, we developed a novel drug response prediction algorithm called Drug Efficacy Leveraging Forked and Specialized networks (DELFOS). Our model learns from multi-omics data from over 65 cancer cell lines, as well as structural data from over 200 compounds, for the prediction of drug sensitivity. We also evaluated the benefits of incorporating single-cell expression data to predict drug response. DELFOS was validated using datasets with unseen cell lines or drugs and compared with other state-of-the-art algorithms, achieving a high prediction performance on several correlation and error metrics. Overall, DELFOS can effectively leverage multi-omics data for the prediction of drug responses in thousands of drug-cell line pairs. AVAILABILITY AND IMPLEMENTATION: The DELFOS pipeline and associated data are available at github.com/MoreiraLAB/delfos.

MUG: A mutation overview of GPCR subfamily A17 receptors.

G protein-coupled receptors (GPCRs) mediate several signaling pathways through a general mechanism that involves their activation, upholding a chain of events that lead to the release of molecules responsible for cytoplasmic action and further regulation. These physiological functions can be severely altered by mutations in GPCR genes. GPCRs subfamily A17 (dopamine, serotonin, adrenergic and trace amine receptors) are directly related with neurodegenerative diseases, and as such it is crucial to explore known mutations on these systems and their impact in structure and function. A comprehensive and detailed computational framework - MUG (Mutations Understanding GPCRs) - was constructed, illustrating key reported mutations and their effect on receptors of the subfamily A17 of GPCRs. We explored the type of mutations occurring overall and in the different families of subfamily A17, as well their localization within the receptor and potential effects on receptor functionality. The mutated residues were further analyzed considering their pathogenicity. The results reveal a high diversity of mutations in the GPCR subfamily A17 structures, drawing attention to the considerable number of mutations in conserved residues and domains. Mutated residues were typically hydrophobic residues enriched at the ligand binding pocket and known activating microdomains, which may lead to disruption of receptor function. MUG as an interactive web application is available for the management and visualization of this dataset. We expect that this interactive database helps the exploration of GPCR mutations, their influence, and their familywise and receptor-specific effects, constituting the first step in elucidating their structures and molecules at the atomic level.

SYNPRED: prediction of drug combination effects in cancer using different synergy metrics and ensemble learning.

BACKGROUND: In cancer research, high-throughput screening technologies produce large amounts of multiomics data from different populations and cell types. However, analysis of such data encounters difficulties due to disease heterogeneity, further exacerbated by human biological complexity and genomic variability. The specific profile of cancer as a disease (or, more realistically, a set of diseases) urges the development of approaches that maximize the effect while minimizing the dosage of drugs. Now is the time to redefine the approach to drug discovery, bringing an artificial intelligence (AI)-powered informational view that integrates the relevant scientific fields and explores new territories. RESULTS: Here, we show SYNPRED, an interdisciplinary approach that leverages specifically designed ensembles of AI algorithms, as well as links omics and biophysical traits to predict anticancer drug synergy. It uses 5 reference models (Bliss, Highest Single Agent, Loewe, Zero Interaction Potency, and Combination Sensitivity Score), which, coupled with AI algorithms, allowed us to attain the ones with the best predictive performance and pinpoint the most appropriate reference model for synergy prediction, often overlooked in similar studies. By using an independent test set, SYNPRED exhibits state-of-the-art performance metrics either in the classification (accuracy, 0.85; precision, 0.91; recall, 0.90; area under the receiver operating characteristic, 0.80; and F1-score, 0.91) or in the regression models, mainly when using the Combination Sensitivity Score synergy reference model (root mean square error, 11.07; mean squared error, 122.61; Pearson, 0.86; mean absolute error, 7.43; Spearman, 0.87). Moreover, data interpretability was achieved by deploying the most current and robust feature importance approaches. A simple web-based application was constructed, allowing easy access by nonexpert researchers. CONCLUSIONS: The performance of SYNPRED rivals that of the existing methods that tackle the same problem, yielding unbiased results trained with one of the most comprehensive datasets available (NCI ALMANAC). The leveraging of different reference models allowed deeper insights into which of them can be more appropriately used for synergy prediction. The Combination Sensitivity Score clearly stood out with improved performance among the full scope of surveyed approaches and synergy reference models. Furthermore, SYNPRED takes a particular focus on data interpretability, which has been in the spotlight lately when using the most advanced AI techniques.

Decoding Partner Specificity of Opioid Receptor Family.

This paper describes an exciting big data analysis compiled in a freely available database, which can be applied to characterize the coupling of different G-Protein coupled receptors (GPCRs) families with their intracellular partners. Opioid receptor (OR) family was used as case study in order to gain further insights into the physiological properties of these important drug targets, known to be associated with the opioid crisis, a huge socio-economic issue directly related to drug abuse. An extensive characterization of all members of the ORs family (µ (MOR), δ (DOR), κ (KOR), nociceptin (NOP)) and their corresponding binding partners (ARRs: Arr2, Arr3; G-protein: Gi1, Gi2, Gi3, Go, Gob, Gz, Gq, G11, G14, G15, G12, Gssh, Gslo) was carried out. A multi-step approach including models' construction (multiple sequence alignment, homology modeling), complex assembling (protein complex refinement with HADDOCK and complex equilibration), and protein-protein interface characterization (including both structural and dynamics analysis) were performed. Our database can be easily applied to several GPCR sub-families, to determine the key structural and dynamical determinants involved in GPCR coupling selectivity.

MENSAdb: a thorough structural analysis of membrane protein dimers.

Membrane proteins (MPs) are key players in a variety of different cellular processes and constitute the target of around 60% of all Food and Drug Administration-approved drugs. Despite their importance, there is still a massive lack of relevant structural, biochemical and mechanistic information mainly due to their localization within the lipid bilayer. To help fulfil this gap, we developed the MEmbrane protein dimer Novel Structure Analyser database (MENSAdb). This interactive web application summarizes the evolutionary and physicochemical properties of dimeric MPs to expand the available knowledge on the fundamental principles underlying their formation. Currently, MENSAdb contains features of 167 unique MPs (63% homo- and 37% heterodimers) and brings insights into the conservation of residues, accessible solvent area descriptors, average B-factors, intermolecular contacts at 2.5 Å and 4.0 Å distance cut-offs, hydrophobic contacts, hydrogen bonds, salt bridges, π-π stacking, T-stacking and cation-π interactions. The regular update and organization of all these data into a unique platform will allow a broad community of researchers to collect and analyse a large number of features efficiently, thus facilitating their use in the development of prediction models associated with MPs. Database URL: http://www.moreiralab.com/resources/mensadb.

Predicting Hot Spots Using a Deep Neural Network Approach.

Targeting protein-protein interactions is a challenge and crucial task of the drug discovery process. A good starting point for rational drug design is the identification of hot spots (HS) at protein-protein interfaces, typically conserved residues that contribute most significantly to the binding. In this chapter, we depict point-by-point an in-house pipeline used for HS prediction using only sequence-based features from the well-known SpotOn dataset of soluble proteins (Moreira et al., Sci Rep 7:8007, 2017), through the implementation of a deep neural network. The presented pipeline is divided into three steps: (1) feature extraction, (2) deep learning classification, and (3) model evaluation. We present all the available resources, including code snippets, the main dataset, and the free and open-source modules/packages necessary for full replication of the protocol. The users should be able to develop an HS prediction model with accuracy, precision, recall, and AUROC of 0.96, 0.93, 0.91, and 0.86, respectively.

The Central Role of Non-Structural Protein 1 (NS1) in Influenza Biology and Infection.

Influenza (flu) is a contagious viral disease, which targets the human respiratory tract and spreads throughout the world each year. Every year, influenza infects around 10% of the world population and between 290,000 and 650,000 people die from it according to the World Health Organization (WHO). Influenza viruses belong to the Orthomyxoviridae family and have a negative sense eight-segment single-stranded RNA genome that encodes 11 different proteins. The only control over influenza seasonal epidemic outbreaks around the world are vaccines, annually updated according to viral strains in circulation, but, because of high rates of mutation and recurrent genetic assortment, new viral strains of influenza are constantly emerging, increasing the likelihood of pandemics. Vaccination effectiveness is limited, calling for new preventive and therapeutic approaches and a better understanding of the virus-host interactions. In particular, grasping the role of influenza non-structural protein 1 (NS1) and related known interactions in the host cell is pivotal to better understand the mechanisms of virus infection and replication, and thus propose more effective antiviral approaches. In this review, we assess the structure of NS1, its dynamics, and multiple functions and interactions, to highlight the central role of this protein in viral biology and its potential use as an effective therapeutic target to tackle seasonal and pandemic influenza.

An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population.

Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.

A Complete Assessment of Dopamine Receptor- Ligand Interactions through Computational Methods.

Background: Selectively targeting dopamine receptors (DRs) has been a persistent challenge in the last years for the development of new treatments to combat the large variety of diseases involving these receptors. Although, several drugs have been successfully brought to market, the subtype-specific binding mode on a molecular basis has not been fully elucidated. Methods: Homology modeling and molecular dynamics were applied to construct robust conformational models of all dopamine receptor subtypes (D1-like and D2-like). Fifteen structurally diverse ligands were docked. Contacts at the binding pocket were fully described in order to reveal new structural findings responsible for selective binding to DR subtypes. Results: Residues of the aromatic microdomain were shown to be responsible for the majority of ligand interactions established to all DRs. Hydrophobic contacts involved a huge network of conserved and non-conserved residues between three transmembrane domains (TMs), TM2-TM3-TM7. Hydrogen bonds were mostly mediated by the serine microdomain. TM1 and TM2 residues were main contributors for the coupling of large ligands. Some amino acid groups form electrostatic interactions of particular importance for D1R-like selective ligands binding. Conclusions: This in silico approach was successful in showing known receptor-ligand interactions as well as in determining unique combinations of interactions, which will support mutagenesis studies to improve the design of subtype-specific ligands.

Structural Characterization of Membrane Protein Dimers.

Membrane proteins are essential vessels for cell communication both with other cells and noncellular structures. They modulate environment responses and mediate a myriad of biological processes. Dimerization and multimerization processes have been shown to further increase the already high specificity of these processes. Due to their central role in various cell and organism functions, these multimers are often associated with health conditions, such as Alzheimer's disease (AD), Parkinson's disease (PD), and diabetes, among others.Understanding the membrane protein dimers' interface takes advantage of the specificity of the structure, for which we must pinpoint the most relevant interfacial residues, since they are extremely likely to be crucial for complex formation. Here, we describe step by step our own in silico protocol to characterize these residues, making use of known experimental structures. We detail the computational pipeline from data acquisition and pre-processing to feature extraction. A molecular dynamics simulation protocol to further study membrane dimer proteins and their interfaces is also illustrated.

Computational Approaches in Antibody-drug Conjugate Optimization for Targeted Cancer Therapy.

Cancer has become one of the main leading causes of morbidity and mortality worldwide. One of the critical drawbacks of current cancer therapeutics has been the lack of the target-selectivity, as these drugs should have an effect exclusively on cancer cells while not perturbing healthy ones. In addition, their mechanism of action should be sufficiently fast to avoid the invasion of neighbouring healthy tissues by cancer cells. The use of conventional chemotherapeutic agents and other traditional therapies, such as surgery and radiotherapy, leads to off-target interactions with serious side effects. In this respect, recently developed target-selective Antibody-Drug Conjugates (ADCs) are more effective than traditional therapies, presumably due to their modular structures that combine many chemical properties simultaneously. In particular, ADCs are made up of three different units: a highly selective Monoclonal antibody (Mab) which is developed against a tumour-associated antigen, the payload (cytotoxic agent), and the linker. The latter should be stable in circulation while allowing the release of the cytotoxic agent in target cells. The modular nature of these drugs provides a platform to manipulate and improve selectivity and the toxicity of these molecules independently from each other. This in turn leads to generation of second- and third-generation ADCs, which have been more effective than the previous ones in terms of either selectivity or toxicity or both. Development of ADCs with improved efficacy requires knowledge at the atomic level regarding the structure and dynamics of the molecule. As such, we reviewed all the most recent computational methods used to attain all-atom description of the structure, energetics and dynamics of these systems. In particular, this includes homology modelling, molecular docking and refinement, atomistic and coarse-grained molecular dynamics simulations, principal component and cross-correlation analysis. The full characterization of the structure-activity relationship devoted to ADCs is critical for antibody-drug conjugate research and development.

SpotOn: High Accuracy Identification of Protein-Protein Interface Hot-Spots.

We present SpotOn, a web server to identify and classify interfacial residues as Hot-Spots (HS) and Null-Spots (NS). SpotON implements a robust algorithm with a demonstrated accuracy of 0.95 and sensitivity of 0.98 on an independent test set. The predictor was developed using an ensemble machine learning approach with up-sampling of the minor class. It was trained on 53 complexes using various features, based on both protein 3D structure and sequence. The SpotOn web interface is freely available at: http://milou.science.uu.nl/services/SPOTON/ .

Membrane proteins structures: A review on computational modeling tools.

BACKGROUND: Membrane proteins (MPs) play diverse and important functions in living organisms. They constitute 20% to 30% of the known bacterial, archaean and eukaryotic organisms' genomes. In humans, their importance is emphasized as they represent 50% of all known drug targets. Nevertheless, experimental determination of their three-dimensional (3D) structure has proven to be both time consuming and rather expensive, which has led to the development of computational algorithms to complement the available experimental methods and provide valuable insights. SCOPE OF REVIEW: This review highlights the importance of membrane proteins and how computational methods are capable of overcoming challenges associated with their experimental characterization. It covers various MP structural aspects, such as lipid interactions, allostery, and structure prediction, based on methods such as Molecular Dynamics (MD) and Machine-Learning (ML). MAJOR CONCLUSIONS: Recent developments in algorithms, tools and hybrid approaches, together with the increase in both computational resources and the amount of available data have resulted in increasingly powerful and trustworthy approaches to model MPs. GENERAL SIGNIFICANCE: Even though MPs are elementary and important in nature, the determination of their 3D structure has proven to be a challenging endeavor. Computational methods provide a reliable alternative to experimental methods. In this review, we focus on computational techniques to determine the 3D structure of MP and characterize their binding interfaces. We also summarize the most relevant databases and software programs available for the study of MPs.