In silico mechanistic profiling to probe small molecule binding to sulfotransferases

Artificial intelligence to predict inhibitors of drug-metabolizing enzymes and transporters for safer drug design

INTRODUCTION

Drug-metabolizing enzymes (DMEs) and transporters (DTs) play integral roles in drug metabolism and drug-drug interactions (DDIs) which directly impact drug efficacy and safety. It is well-established that inhibition of DMEs and DTs often leads to adverse drug reactions (ADRs) and therapeutic failure. As such, early prediction of such inhibitors is vital in drug development. In this context, the limitations of the traditional in vitro assays and QSAR models methods have been addressed by harnessing artificial intelligence (AI) techniques.

AREAS COVERED

This narrative review presents the insights gained from the application of AI for predicting DME and DT inhibitors over the past decade. Several case studies demonstrate successful AI applications in enzyme-transporter interaction prediction, and the authors discuss workflows for integrating these predictions into drug design and regulatory frameworks.

EXPERT OPINION

The application of AI in predicting DME and DT inhibitors has demonstrated significant potential toward enhancing drug safety and effectiveness. However, critical challenges involve the data quality, biases, and model transparency. The availability of diverse, high-quality datasets alongside the integration of pharmacokinetic and genomic data are essential. Lastly, the collaboration among computational scientists, pharmacologists, and regulatory bodies is pyramidal in tailoring AI tools for personalized medicine and safer drug development.

Polymorphisms in drug-metabolizing genes and urinary bladder cancer susceptibility and prognosis: Possible impacts and future management

Epidemiological studies have shown the association of genetic variants with risks of occupational and environmentally induced cancers, including bladder (BC). The current review summarizes the effects of variants in genes encoding phase I and II enzymes in well-designed studies to highlight their contribution to BC susceptibility and prognosis. Polymorphisms in genes codifying drug-metabolizing proteins are of particular interest because of their involvement in the metabolism of exogenous genotoxic compounds, such as tobacco and agrochemicals. The prognosis between muscle-invasive and non-muscle-invasive diseases is very different, and it is difficult to predict which will progress worse. Web of Science, PubMed, and Medline were searched to identify studies published between January 1, 2010, and February 2023. We included 73 eligible studies, more than 300 polymorphisms, and 46 genes/loci. The most studied candidate genes/loci of phase I metabolism were CYP1B1, CYP1A1, CYP1A2, CYP3A4, CYP2D6, CYP2A6, CYP3E1, and ALDH2, and those in phase II were GSTM1, GSTT1, NAT2, GSTP1, GSTA1, GSTO1, and UGT1A1. We used the 46 genes to construct a network of proteins and to evaluate their biological functions based on the Reactome and KEGG databases. Lastly, we assessed their expression in different tissues, including normal bladder and BC samples. The drug-metabolizing pathway plays a relevant role in BC, and our review discusses a list of genes that could provide clues for further exploration of susceptibility and prognostic biomarkers.

Computational and artificial intelligence-based approaches for drug metabolism and transport prediction

Drug metabolism and transport, orchestrated by drug-metabolizing enzymes (DMEs) and drug transporters (DTs), are implicated in drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Reliable and precise predictions of DDIs and ADRs are critical in the early stages of drug development to reduce the rate of drug candidate failure. A variety of experimental and computational technologies have been developed to predict DDIs and ADRs. Recent artificial intelligence (AI) approaches offer new opportunities for better predicting and understanding the complex processes related to drug metabolism and transport. We summarize the role of major DMEs and DTs, and provide an overview of current progress in computational approaches for the prediction of drug metabolism, transport, and DDIs, with an emphasis on AI including machine learning (ML) and deep learning (DL) modeling.

Role of Conformational Dynamics of Sulfotransferases SULT1A1 and SULT1A3 in Substrate Specificity

Sulfotransferases (SULTs) are phase II metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3'-Phosphoadenosine 5'-Phosphosulfate (PAPS) to a wide variety of endogenous compounds, drugs and natural products. Although SULT1A1 and SULT1A3 share 93% identity, SULT1A1, the most abundant SULT isoform in humans, exhibits a broad substrate range with specificity for small phenolic compounds, while SULT1A3 displays a high affinity toward monoamine neurotransmitters like dopamine. To elucidate the factors determining the substrate specificity of the SULT1 isoenzymes, we studied the dynamic behavior and structural specificities of SULT1A1 and SULT1A3 by using molecular dynamics (MD) simulations and ensemble docking of common and specific substrates of the two isoforms. Our results demonstrated that while SULT1A1 exhibits a relatively rigid structure by showing lower conformational flexibility except for the lip (loop L1), the loop L2 and the cap (L3) of SULT1A3 are extremely flexible. We identified protein residues strongly involved in the recognition of different substrates for the two isoforms. Our analyses indicated that being more specific and highly flexible, the structure of SULT1A3 has particularities in the binding site, which are crucial for its substrate selectivity.

Predicting Regioselectivity of Cytosolic Sulfotransferase Metabolism for Drugs

Cytosolic sulfotransferases (SULTs) are a family of enzymes responsible for the sulfation of small endogenous and exogenous compounds. SULTs contribute to the conjugation phase of metabolism and share substrates with the uridine 5'-diphospho-glucuronosyltransferase (UGT) family of enzymes. UGTs are considered to be the most important enzymes in the conjugation phase, and SULTs are an auxiliary enzyme system to them. Understanding how the regioselectivity of SULTs differs from that of UGTs is essential from the perspective of developing novel drug candidates. We present a general ligand-based SULT model trained and tested using high-quality experimental regioselectivity data. The current study suggests that, unlike other metabolic enzymes in the modification and conjugation phases, the SULT regioselectivity is not strongly influenced by the activation energy of the rate-limiting step of the catalysis. Instead, the prominent role is played by the substrate binding site of SULT. Thus, the model is trained only on steric and orientation descriptors, which mimic the binding pocket of SULT. The resulting classification model, which predicts whether a site is metabolized, achieved a Cohen's kappa of 0.71.

Machine learning and structure-based modeling for the prediction of UDP-glucuronosyltransferase inhibition

UDP-glucuronosyltransferases (UGTs) are responsible for 35% of the phase II drug metabolism. In this study, we focused on UGT1A1, which is a key UGT isoform. Strong inhibition of UGT1A1 may trigger adverse drug/herb-drug interactions, or result in disorders of endobiotic metabolism. Most of the current machine learning methods predicting the inhibition of drug metabolizing enzymes neglect protein structure and dynamics, both being essential for the recognition of various substrates and inhibitors. We performed molecular dynamics simulations on a homology model of the human UGT1A1 structure containing both the cofactor- (UDP-glucuronic acid) and substrate-binding domains to explore UGT conformational changes. Then, we created models for the prediction of UGT1A1 inhibitors by integrating information on UGT1A1 structure and dynamics, interactions with diverse ligands, and machine learning. These models can be helpful for further prediction of drug-drug interactions of drug candidates and safety treatments.

•

UGTs are responsible for 35% of the phase II drug metabolism reactions

•

We created machine learning models for prediction of UGT1A1 inhibitors

•

Our simulations suggested key residues of UGT1A1 involved in the substrate binding

Pharmacogenetics of human sulfotransferases and impact of amino acid exchange on Phase II drug metabolism

Sulfotransferases (SULTs) are Phase II drug-metabolizing enzymes (DMEs) catalyzing the sulfation of a variety of endogenous compounds, natural products, and drugs. Various drugs, such as nonsteroidal anti-inflammatory drugs (NSAIDS) can inhibit SULTs, affecting drug-drug interactions. Several polymorphisms have been identified for SULTs that might be crucial for interindividual variability in drug response and toxicity or for increased disease risk. Here, we review current knowledge on non-synonymous single nucleotide polymorphisms (nsSNPs) of human SULTs, focusing on the coded SULT allozymes and molecular mechanisms explaining their variable activity, which is essential for personalized medicine. We discuss the structural and dynamic bases of key amino acid (AA) variants implicated in the impacts on drug metabolism in the case of SULT1A1, as revealed by molecular modeling approaches.

From Steroid and Drug Metabolism to Glycobiology, Using Sulfotransferase Structures to Understand and Tailor Function

Sulfotransferases are ubiquitous enzymes that transfer a sulfo group from the universal cofactor donor 3'-phosphoadenosine 5'-phosphosulfate to a broad range of acceptor substrates. In humans, the cytosolic sulfotransferases are involved in the sulfation of endogenous compounds such as steroids, neurotransmitters, hormones, and bile acids as well as xenobiotics including drugs, toxins, and environmental chemicals. The Golgi associated membrane-bound sulfotransferases are involved in post-translational modification of macromolecules from glycosaminoglycans to proteins. The sulfation of small molecules can have profound biologic effects on the functionality of the acceptor, including activation, deactivation, or enhanced metabolism and elimination. Sulfation of macromolecules has been shown to regulate a number of physiologic and pathophysiological pathways by enhancing binding affinity to regulatory proteins or binding partners. Over the last 25 years, crystal structures of these enzymes have provided a wealth of information on the mechanisms of this process and the specificity of these enzymes. This review will focus on the general commonalities of the sulfotransferases, from enzyme structure to catalytic mechanism as well as providing examples into how structural information is being used to either design drugs that inhibit sulfotransferases or to modify the enzymes to improve drug synthesis. SIGNIFICANCE STATEMENT: This manuscript honors Dr. Masahiko Negishi's contribution to the understanding of sulfotransferase mechanism, specificity, and roles in biology by analyzing the crystal structures that have been solved over the last 25 years.

Computational Analysis of Chemical Space of Natural Compounds Interacting with Sulfotransferases

The aim of this study was to investigate the chemical space and interactions of natural compounds with sulfotransferases (SULTs) using ligand- and structure-based in silico methods. An in-house library of natural ligands (hormones, neurotransmitters, plant-derived compounds and their metabolites) reported to interact with SULTs was created. Their chemical structures and properties were compared to those of compounds of non-natural (synthetic) origin, known to interact with SULTs. The natural ligands interacting with SULTs were further compared to other natural products for which interactions with SULTs were not known. Various descriptors of the molecular structures were calculated and analyzed. Statistical methods (ANOVA, PCA, and clustering) were used to explore the chemical space of the studied compounds. Similarity search between the compounds in the different groups was performed with the ROCS software. The interactions with SULTs were additionally analyzed by docking into different experimental and modeled conformations of SULT1A1. Natural products with potentially strong interactions with SULTs were outlined. Our results contribute to a better understanding of chemical space and interactions of natural compounds with SULT enzymes and help to outline new potential ligands of these enzymes.

Insights into the substrate binding mechanism of SULT1A1 through molecular dynamics with excited normal modes simulations

Sulfotransferases (SULTs) are phase II drug-metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3′-phosphoadenosine 5′-phosphosulfate (PAPS) to a substrate. It has been previously suggested that a considerable shift of SULT structure caused by PAPS binding could control the capability of SULT to bind large substrates. We employed molecular dynamics (MD) simulations and the recently developed approach of MD with excited normal modes (MDeNM) to elucidate molecular mechanisms guiding the recognition of diverse substrates and inhibitors by SULT1A1. MDeNM allowed exploring an extended conformational space of PAPS-bound SULT1A1, which has not been achieved up to now by using classical MD. The generated ensembles combined with docking of 132 SULT1A1 ligands shed new light on substrate and inhibitor binding mechanisms. Unexpectedly, our simulations and analyses on binding of the substrates estradiol and fulvestrant demonstrated that large conformational changes of the PAPS-bound SULT1A1 could occur independently of the co-factor movements that could be sufficient to accommodate large substrates as fulvestrant. Such structural displacements detected by the MDeNM simulations in the presence of the co-factor suggest that a wider range of drugs could be recognized by PAPS-bound SULT1A1 and highlight the utility of including MDeNM in protein–ligand interactions studies where major rearrangements are expected.

An activity-based fluorescent sensor for the detection of the phenol sulfotransferase SULT1A1 in living cells

Human phenol sulfotransferases mediate the transfer of a sulfuryl moiety from the activated sulfate donor PAPS to hydroxy-containing substrates, altering substrate solubility and charge to affect phase II metabolism and cell signaling. Here, we present the development, computational modeling, in vitro enzymology, and biological application of STS-3, an activity-based fluorescent sensor for the SULT1A1 isoform.

Discovery of small molecule inhibitors of chikungunya virus proteins (nsP2 and E1) using in silico approaches

Chikungunya virus (CHIKV) has emerged as a major viral threat, affecting over a million people worldwide per year. It is a vector borne disease transmitted to the human by Ades mosquitoes and primarily affect people by causing viral fever, severe joint pain and other symptoms, like rash, joint swelling, muscle pain and in rare cases can be fatal. CHIKV is a deadly virus, with its mutation rate found to be significantly higher as compared to other viruses. To date, there has been no reported FDA approved drug against this virus. Thus, keeping in mind the urgent need to scrutinize potential therapies against CHIKV, the present study identified twenty plant bioactive compounds that are available at low price and do not have associated adverse effect. For identification of active potentials molecules the pharmacoinformatics-based perspective was applied against CHIKV structural (E1) and non-structural (nsP2) proteins using molecular docking and scoring. The selected compounds were further studied for pharmacokinetics (PK) and pharmacodynamics (PD) associated parameters such as initial absorption, then distribution and later on metabolism excretion and toxicity (ADMET) profiles based on in silico study. The results reveal five potential lead compounds having high binding energy that can help in the development of commercial drugs with favorable ADMET characteristic.Communicated by Ramaswamy H. Sarma.

Mechanistic Insights into the Effect of Ligands on Structural Stability and Selectivity of Sulfotransferase 2A1 (SULT2A1)

Cytosolic sulfotransferases (SULTs) acting as phase II metabolic enzymes can be used in the sulfonation of small molecules by transferring a sulfonate group from the unique co-factor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to the substrates. In the present study, molecular dynamics (MD) simulations and ensemble docking study were employed to theoretically characterize the mechanism for the effect of co-factor (PAP) and ligands (LCA, raloxifene, α-hydroxytamoxifen, ouabain, and 3'-phosphoadenylyl sulfate) on structural stability and selectivity of SULT2A1 from the perspective of the dynamic behavior of SULT2A1 structures. Structural stability and network analyses indicated that the cooperation between PAP and LCA may enhance the thermal stability and compact communication in enzymes. During the MD simulations, the obviously rigid region and inward displacement were detected in the active-site cap (loop16) of the conformation containing PAP, which may be responsible for the significant changes in substrate accessibility and catalytic activity. The smaller substrates such as LCA could bind stably to the active pocket in the presence of PAP. However, the substrates or inhibitors with a large spatial structure needed to bind to the open conformation (without PAP) prior to PAPS binding.

The Human Skeletal Muscle Proteome Project: a reappraisal of the current literature

Skeletal muscle is a large organ that accounts for up to half the total mass of the human body. A progressive decline in muscle mass and strength occurs with ageing and in some individuals configures the syndrome of 'sarcopenia', a condition that impairs mobility, challenges autonomy, and is a risk factor for mortality. The mechanisms leading to sarcopenia as well as myopathies are still little understood. The Human Skeletal Muscle Proteome Project was initiated with the aim to characterize muscle proteins and how they change with ageing and disease. We conducted an extensive review of the literature and analysed publically available protein databases. A systematic search of peer-reviewed studies was performed using PubMed. Search terms included 'human', 'skeletal muscle', 'proteome', 'proteomic(s)', and 'mass spectrometry', 'liquid chromatography-mass spectrometry (LC-MS/MS)'. A catalogue of 5431 non-redundant muscle proteins identified by mass spectrometry-based proteomics from 38 peer-reviewed scientific publications from 2002 to November 2015 was created. We also developed a nosology system for the classification of muscle proteins based on localization and function. Such inventory of proteins should serve as a useful background reference for future research on changes in muscle proteome assessed by quantitative mass spectrometry-based proteomic approaches that occur with ageing and diseases. This classification and compilation of the human skeletal muscle proteome can be used for the identification and quantification of proteins in skeletal muscle to discover new mechanisms for sarcopenia and specific muscle diseases that can be targeted for the prevention and treatment.

In Silico Prediction of Human Sulfotransferase 1E1 Activity Guided by Pharmacophores from Molecular Dynamics Simulations

Acting during phase II metabolism, sulfotransferases (SULTs) serve detoxification by transforming a broad spectrum of compounds from pharmaceutical, nutritional, or environmental sources into more easily excretable metabolites. However, SULT activity has also been shown to promote formation of reactive metabolites that may have genotoxic effects. SULT subtype 1E1 (SULT1E1) was identified as a key player in estrogen homeostasis, which is involved in many physiological processes and the pathogenesis of breast and endometrial cancer. The development of an in silico prediction model for SULT1E1 ligands would therefore support the development of metabolically inert drugs and help to assess health risks related to hormonal imbalances. Here, we report on a novel approach to develop a model that enables prediction of substrates and inhibitors of SULT1E1. Molecular dynamics simulations were performed to investigate enzyme flexibility and sample protein conformations. Pharmacophores were developed that served as a cornerstone of the model, and machine learning techniques were applied for prediction refinement. The prediction model was used to screen the DrugBank (a database of experimental and approved drugs): 28% of the predicted hits were reported in literature as ligands of SULT1E1. From the remaining hits, a selection of nine molecules was subjected to biochemical assay validation and experimental results were in accordance with the in silico prediction of SULT1E1 inhibitors and substrates, thus affirming our prediction hypotheses.

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Aim: Molecular dynamics simulations and normal mode analysis are well-established approaches to generate receptor conformational ensembles (RCEs) for ligand docking and virtual screening. Here, we report new fast molecular dynamics-based and normal mode analysis-based protocols combined with conformational pocket classifications to efficiently generate RCEs. Materials & Methods: We assessed our protocols on two well-characterized protein targets showing local active site flexibility, dihydrofolate reductase and large collective movements, CDK2. The performance of the RCEs was validated by distinguishing known ligands of dihydrofolate reductase and CDK2 among a dataset of diverse chemical decoys. Results & discussion: Our results show that different simulation protocols can be efficient for generation of RCEs depending on different kind of protein flexibility.[Formula: see text]

The impact of molecular dynamics on drug design: applications for the characterization of ligand-macromolecule complexes

Among all tools available to design new drugs, molecular dynamics (MD) simulations have become an essential technique. Initially developed to investigate molecular models with a limited number of atoms, computers now enable investigations of large macromolecular systems with a simulation time reaching the microsecond range. The reviewed articles cover four years of research to give an overview on the actual impact of MD on the current medicinal chemistry landscape with a particular emphasis on studies of ligand-protein interactions. With a special focus on studies combining computational approaches with data gained from other techniques, this review shows how deeply embedded MD simulations are in drug design strategies and articulates what the future of this technique could be.

Rational design of small-molecule stabilizers of spermine synthase dimer by virtual screening and free energy-based approach

Snyder-Robinson Syndrome (SRS) is a rare mental retardation disorder which is caused by the malfunctioning of an enzyme, the spermine synthase (SMS), which functions as a homo-dimer. The malfunctioning of SMS in SRS patients is associated with several identified missense mutations that occur away from the active site. This investigation deals with a particular SRS-causing mutation, the G56S mutation, which was shown computationally and experimentally to destabilize the SMS homo-dimer and thus to abolish SMS enzymatic activity. As a proof-of-concept, we explore the possibility to restore the enzymatic activity of the malfunctioning SMS mutant G56S by stabilizing the dimer through small molecule binding at the mutant homo-dimer interface. For this purpose, we designed an in silico protocol that couples virtual screening and a free binding energy-based approach to identify potential small-molecule binders on the destabilized G56S dimer, with the goal to stabilize it and thus to increase SMS G56S mutant activity. The protocol resulted in extensive list of plausible stabilizers, among which we selected and tested 51 compounds experimentally for their capability to increase SMS G56S mutant enzymatic activity. In silico analysis of the experimentally identified stabilizers suggested five distinctive chemical scaffolds. This investigation suggests that druggable pockets exist in the vicinity of the mutation sites at protein-protein interfaces which can be used to alter the disease-causing effects by small molecule binding. The identified chemical scaffolds are drug-like and can serve as original starting points for development of lead molecules to further rescue the disease-causing effects of the Snyder-Robinson syndrome for which no efficient treatment exists up to now.