Message Passing Neural Networks Improve Prediction of Metabolite Authenticity.

Cytochrome P450 enzymes aid in the elimination of a preponderance of small molecule drugs, but can generate reactive metabolites that may adversely react with protein and DNA and prompt drug candidate attrition or market withdrawal. Previously developed models help understand how these enzymes modify molecule structure by predicting sites of metabolism or characterizing formation of metabolite-biomolecule adducts. However, the majority of reactive metabolites are formed by multiple metabolic steps, and understanding the progenitor molecule's network-level behavior necessitates an integrative approach that blends multiple site of metabolism and structure inference models. Our previously developed tool, XenoNet 1.0, generates metabolic networks, where nodes are molecules and weighted edges are metabolic transformations. We extend XenoNet with a bidirectional message passing neural network that integrates edge feature information and local network structure using edge-conditioned graph convolutions and jumping knowledge to predict the authenticity of inferred Phase I metabolite structures. Our model significantly outperformed prior work and algorithmic baselines on a data set of 311 networks and 6606 intermediates annotated using a chemically diverse set of 20 736 individual in vitro and in vivo reaction records accounting for 92.3% of all human Phase I metabolism in the Accelrys Metabolite Database. Cross-validated predictions resulted in area under the receiver operating characteristic curves of 88.5% and 87.6% for separating experimentally observed and unobserved metabolites at global and network levels, respectively. Further analysis verified robustness to networks of varying depth and breadth, accurate detection of metabolites, such as d,l-methamphetamine, that are experimentally observed or unobserved in different network contexts, extraction of important metabolic subnetworks, and identification of known bioactivation pathways, such as for nimesulide and terbinafine. By exploiting network structures, our approach accurately suggests unreported metabolites for experimental study and may rationalize modifications for avoiding deleterious pathways antecedent to reactive metabolite formation.

Membrane Adaptations and Cellular Responses of Sulfolobus acidocaldarius to the Allylamine Terbinafine.

Cellular membranes are essential for compartmentalization, maintenance of permeability, and fluidity in all three domains of life. Archaea belong to the third domain of life and have a distinct phospholipid composition. Membrane lipids of archaea are ether-linked molecules, specifically bilayer-forming dialkyl glycerol diethers (DGDs) and monolayer-forming glycerol dialkyl glycerol tetraethers (GDGTs). The antifungal allylamine terbinafine has been proposed as an inhibitor of GDGT biosynthesis in archaea based on radiolabel incorporation studies. The exact target(s) and mechanism of action of terbinafine in archaea remain elusive. Sulfolobus acidocaldarius is a strictly aerobic crenarchaeon thriving in a thermoacidophilic environment, and its membrane is dominated by GDGTs. Here, we comprehensively analyzed the lipidome and transcriptome of S. acidocaldarius in the presence of terbinafine. Depletion of GDGTs and the accompanying accumulation of DGDs upon treatment with terbinafine were growth phase-dependent. Additionally, a major shift in the saturation of caldariellaquinones was observed, which resulted in the accumulation of unsaturated molecules. Transcriptomic data indicated that terbinafine has a multitude of effects, including significant differential expression of genes in the respiratory complex, motility, cell envelope, fatty acid metabolism, and GDGT cyclization. Combined, these findings suggest that the response of S. acidocaldarius to terbinafine inhibition involves respiratory stress and the differential expression of genes involved in isoprenoid biosynthesis and saturation.

Antibiofilm and staphyloxanthin inhibitory potential of terbinafine against Staphylococcus aureus: in vitro and in vivo studies.

BACKGROUND: Antimicrobial resistance is growing substantially, which necessitates the search for novel therapeutic options. Terbinafine, an allylamine antifungal agent that exhibits a broad spectrum of activity and is used in the treatment of dermatophytosis, could be a possible option to disarm S. aureus virulence. METHODS: Terbinafine inhibitory effect on staphyloxanthin was characterized by quantitative measurement of staphyloxanthin intermediates and molecular docking. The effect of terbinafine on S. aureus stress survival was characterized by viable counting. The anti-biofilm activity of terbinafine on S. aureus was assessed by the crystal violet assay and microscopy. Changes in S. aureus membrane following treatment with terbinafine were determined using Fourier transform infrared (FTIR) analysis. The synergistic action of terbinafine in combination with conventional antibiotics was characterized using the checkerboard assay. qRT-PCR was used to evaluate the impact of terbinafine on S. aureus gene expression. The influence of terbinafine on S. aureus pathogenesis was investigated in mice infection model. RESULTS: Terbinafine inhibits staphyloxanthin biosynthesis through targeting dehydrosqualene desaturase (CrtN). Docking analysis of terbinafine against the predicted active site of CrtN reveals a binding energy of - 9.579 kcal/mol exemplified by the formation of H-bonds, H-arene bonds, and hydrophobic/hydrophilic interactions with the conserved amino acids of the receptor pocket. Terbinafine treated S. aureus was more susceptible to both oxidative and acid stress as well as human blood killing as compared to untreated cells. Targeting staphyloxanthin by terbinafine rendered S. aureus more sensitive to membrane acting antibiotics. Terbinafine interfered with S. aureus biofilm formation through targeting cell autoaggregation, hydrophobicity, and exopolysaccharide production. Moreover, terbinafine demonstrated a synergistic interaction against S. aureus when combined with conventional antibiotics. Importantly, terbinafine attenuated S. aureus pathogenesis using mice infection model. qRT-PCR revealed that terbinafine repressed expression of the transcriptional regulators sigB, sarA, and msaB, as well as icaA in S. aureus. CONCLUSIONS: Present findings strongly suggest that terbinafine could be used safely and efficiently as an anti-virulent agent to combat S. aureus infections.

Miconazole and terbinafine induced reactive oxygen species accumulation and topical toxicity in human keratinocytes.

There are an estimated 1 billion cases of superficial fungal infection globally. Fungal pathogens form biofilms within wounds and delay the wound healing process. Miconazole and terbinafine are commonly used to treat fungal infections. They induce the accumulation of reactive oxygen species (ROS) in fungi, resulting in the death of fungal cells. ROS are highly reactive molecules, such as oxygen (O2), superoxide anion (O2•-), hydrogen peroxide (H2O2) and hydroxyl radicals (•OH). Although ROS generation is useful for killing pathogenic fungi, it is cytotoxic to human keratinocytes. To the best of our knowledge, the effect of miconazole and terbinafine on HaCaT cells has not been studied with respect to intracellular ROS stimulation. We hypothesized that miconazole and terbinafine have anti-wound healing effects on skin cells when used in antifungal treatment because they generate ROS in fungal cells. We used sulforhodamine B protein staining to investigate cytotoxicity and 2',7'-dichlorofluorescein diacetate to determine ROS accumulation at the 50% inhibitory concentrations of miconazole and terbinafine in HaCaT cells. Our preliminary results showed that topical treatment with miconazole and terbinafine induced cytotoxic responses, with miconazole showing higher cytotoxicity than terbinafine. Both the treatments stimulated ROS in keratinocytes, which may induce oxidative stress and cell death. This suggests a negative correlation between intracellular ROS accumulation in keratinocytes treated with miconazole or terbinafine and the healing of fungi-infected skin wounds.

Engineering of Saccharomyces cerevisiae for the production of (+)-ambrein.

The triterpenoid (+)-ambrein is the major component of ambergris, a coprolite of the sperm whale that can only be rarely found on shores. Upon oxidative degradation of (+)-ambrein, several fragrance molecules are formed, amongst them (-)-ambrox, one of the highest valued compounds in the perfume industry. In order to generate a Saccharomyces cerevisiae whole-cell biocatalyst for the production of (+)-ambrein, intracellular supply of the squalene was enhanced by overexpression of two central enzymes in the mevalonate and sterol biosynthesis pathway, namely the N-terminally truncated 3-hydroxy-3-methylglutaryl-CoA reductase 1 (tHMG) and the squalene synthase (ERG9). In addition, another key enzyme in sterol biosynthesis, squalene epoxidase (ERG1) was inhibited by an experimentally defined amount of the inhibitor terbinafine in order to reduce flux of squalene towards ergosterol biosynthesis while retaining sufficient activity to maintain cell viability and growth. Heterologous expression of a promiscuous variant of Bacillus megaterium tetraprenyl-ß-curcumene cyclase (BmeTC-D373C), which has been shown to be able to catalyse the conversion of squalene to 3-deoxyachillol and then further to (+)-ambrein resulted in production of these triterpenoids in S. cerevisiae for the first time. Triterpenoid yields are comparable with the best microbial production chassis described in literature so far, the methylotrophic yeast Pichia pastoris. Consequently, we discuss similarities and differences of these two yeast species when applied for whole-cell (+)-ambrein production.

Trichophyton rubrum Azole Resistance Mediated by a New ABC Transporter, TruMDR3.

The mechanisms of terbinafine resistance in a set of clinical isolates of Trichophyton rubrum have been studied recently. Of these isolates, TIMM20092 also showed reduced sensitivity to azoles. The azole resistance of TIMM20092 could be inhibited by milbemycin oxime, prompting us to examine the potential of T. rubrum to develop resistance through multidrug efflux transporters. The introduction of a T. rubrum cDNA library into Saccharomyces cerevisiae allowed the isolation of one transporter of the major facilitator superfamily (MFS) conferring resistance to azoles (TruMFS1). To identify more azole efflux pumps among 39 ABC and 170 MFS transporters present within the T. rubrum genome, we performed a BLASTp analysis of Aspergillus fumigatus, Candida albicans, and Candida glabrata on transporters that were previously shown to confer azole resistance. The identified candidates were further tested by heterologous gene expression in S. cerevisiae Four ABC transporters (TruMDR1, TruMDR2, TruMDR3, and TruMDR5) and a second MFS transporter (TruMFS2) proved to be able to operate as azole efflux pumps. Milbemycin oxime inhibited only TruMDR3. Expression analysis showed that both TruMDR3 and TruMDR2 were significantly upregulated in TIMM20092. TruMDR3 transports voriconazole (VRC) and itraconazole (ITC), while TruMDR2 transports only ITC. Disruption of TruMDR3 in TIMM20092 abolished its resistance to VRC and reduced its resistance to ITC. Our study highlights TruMDR3, a newly identified transporter of the ABC family in T. rubrum, which can confer azole resistance if overexpressed. Finally, inhibition of TruMDR3 by milbemycin suggests that milbemycin analogs could be interesting compounds to treat dermatophyte infections in cases of azole resistance.

CYP2C19 and 3A4 Dominate Metabolic Clearance and Bioactivation of Terbinafine Based on Computational and Experimental Approaches.

Lamisil (terbinafine) is an effective, widely prescribed antifungal drug that causes rare idiosyncratic hepatotoxicity. The proposed toxic mechanism involves a reactive metabolite, 6,6-dimethyl-2-hepten-4-ynal (TBF-A), formed through three N-dealkylation pathways. We were the first to characterize them using in vitro studies with human liver microsomes and modeling approaches, yet knowledge of the individual enzymes catalyzing reactions remained unknown. Herein, we employed experimental and computational tools to assess terbinafine metabolism by specific cytochrome P450 isozymes. In vitro inhibitor phenotyping studies revealed six isozymes were involved in one or more N-dealkylation pathways. CYP2C19 and 3A4 contributed to all pathways, and so, we targeted them for steady-state analyses with recombinant isozymes. N-Dealkylation yielding TBF-A directly was catalyzed by CYP2C19 and 3A4 similarly. Nevertheless, CYP2C19 was more efficient than CYP3A4 at N-demethylation and other steps leading to TBF-A. Unlike microsomal reactions, N-denaphthylation was surprisingly efficient for CYP2C19 and 3A4, which was validated by controls. CYP2C19 was the most efficient among all reactions. Nonetheless, CYP3A4 was more selective at steps leading to TBF-A, making it more effective in terbinafine bioactivation based on metabolic split ratios for competing pathways. Model predictions did not extrapolate to quantitative kinetic constants, yet some results for CYP3A4 and CYP2C19 agreed qualitatively with preferred reaction steps and pathways. Clinical data on drug interactions support the CYP3A4 role in terbinafine metabolism, while CYP2C19 remains understudied. Taken together, knowledge of P450s responsible for terbinafine metabolism and TBF-A formation provides a foundation for investigating and mitigating the impact of P450 variations in toxic risks posed to patients.

Identification of novel glutathione conjugates of terbinafine in liver microsomes and hepatocytes across species.

1. Terbinafine (TBF), a common antifungal agent, has been associated with rare incidences of hepatotoxicity. It is hypothesized that bioactivation of TBF to reactive intermediates and subsequent binding to critical cellular proteins may contribute to this toxicity. In the present study, we have characterized the bioactivation pathways of TBF extensively in human, mouse, monkey, dog and rat liver microsomes and hepatocytes. 2. A total of twenty glutathione conjugates of TBF were identified in hepatocytes; thirteen of these conjugates were also detected in liver microsomes. To the best of our knowledge, only two of these conjugates have been reported previously. The conjugates were categorized into three groups based on their mechanism of formation: (a) alkene/alkyne oxidation followed by glutathione conjugation, with or without N-demethylation, (b) arene oxidation followed by glutathione conjugation, with or without N-demethylation, and (c) N-dealkylation followed by glutathione conjugation of the allylic aldehyde, alcohol and acid intermediates. 3. Differences were observed across species in the contributions of these pathways toward overall metabolic turnover. We conclude that, in addition to the glutathione conjugates known to form by Michael addition to the allylic aldehyde, there are other pathways involving the formation of arene oxides and alkene/alkyne epoxides that may be relevant to the discussion of TBF-mediated idiosyncratic drug reactions.

Quantitative Investigation of Terbinafine Hydrochloride Absorption into a Living Skin Equivalent Model by MALDI-MSI.

The combination of microspotting of analytical and internal standards, matrix sublimation, and recently developed software for quantitative mass spectrometry imaging has been used to develop a high-resolution method for the determination of terbinafine hydrochloride in the epidermal region of a full thickness living skin equivalent model. A quantitative assessment of the effect of the addition of the penetration enhancer (dimethyl isosorbide (DMI)) to the delivery vehicle has also been performed, and data have been compared to those obtained from LC-MS/MS measurements of homogenates of isolated epidermal tissue. At 10% DMI, the levels of signal detected for the drug in the epidermis were 0.20 ± 0.072 mg/g tissue for QMSI and 0.28 ± 0.040 mg/g tissue for LC-MS/MS at 50% DMI 0.69 ± 0.23 mg/g tissue for QMSI and 0.66 ± 0.057 mg/g tissue for LC-MS/MS. Comparison of means and standard deviations indicates no significant difference between the values obtained by the two methods.

Modulation of yeast Erg1 expression and terbinafine susceptibility by iron bioavailability.

Ergosterol is a specific sterol component of yeast and fungal membranes. Its biosynthesis is one of the most effective targets for antifungal treatments. However, the emergent resistance to multiple sterol-based antifungal drugs emphasizes the need for new therapeutic approaches. The allylamine terbinafine, which selectively inhibits squalene epoxidase Erg1 within the ergosterol biosynthetic pathway, is mainly used to treat dermatomycoses, whereas its effectiveness in other fungal infections is limited. Given that ergosterol biosynthesis depends on iron as an essential cofactor, in this report, we used the yeast Saccharomyces cerevisiae to investigate how iron bioavailability influences Erg1 expression and terbinafine susceptibility. We observed that both chemical and genetic depletion of iron decrease ERG1 expression, leading to an increase in terbinafine susceptibility. Deletion of either ROX1 transcriptional repressor or CTH1 and CTH2 post-transcriptional repressors of ERG1 expression led to an increase in Erg1 protein levels and terbinafine resistance. On the contrary, overexpression of CTH2 led to the opposite effect, lowering Erg1 levels and increasing terbinafine susceptibility. Although strain-specific particularities exist, opportunistic pathogenic strains of S. cerevisiae displayed a response similar to the laboratory strain. These data indicate that iron bioavailability and particular regulatory factors could be used to modulate susceptibility to terbinafine.

Co-delivery of terbinafine hydrochloride and urea with an in situ film-forming system for nail targeting treatment.

Onychomycosis is a chronic nail disorder consisting of a fungal infection that causes physical and psychosocial discomfort to patients. However, its treatment remains challenging owing to the barrier of the highly keratinized nail plate and the short time that conventional formulations reside on nails. In this work, we developed an in situ film-forming system(IFFS) based on Eudragit® RLPO to co-deliver terbinafine hydrochloride (TBH) and urea, i.e., TBH-urea-RLPO IFFS, with the aim of overcoming the nail barrier, prolonging the residence time, and efficiently treating onychomycosis. The IFFS formulation formed a thin film with good appearance and adhesion upon application in situ. The physical states of TBH and urea in the film were evaluated with polarization microscopy and powder X-ray diffraction. TBH and urea were both amorphousmiscible components within the RLPO film. TBH release from TBH-urea-RLPO IFFS fitted to the Korsmeyer-Pappas model, and the cumulative release at 72 h was significantly higher than that from commercial preparations (Lamisil Pedisan® once). In vitro permeation of TBH from TBH-urea-RLPO IFFS through bovine hoof membranes was evaluated in comparison with the film containing TBH alone (TBH-RLPO) and commercial preparations. The retention and cumulative permeated amount of TBH were significantly enhanced for the TBH-urea-RLPO IFFS (170.80 ± 44.63 µg/cm2vs 75.49 ± 21.50 µg/cm2vs 60.25 ± 27.38 µg/cm2; 61.81 ± 16.09 µg/cm2vs 21.80 ± 11.56 µg/cm2vs 7.91 ± 1.03 µg/cm2, respectively), and the membranes treated with different formulations were observed with SEM and FTIR to identify the denaturing effect of urea on bovine hoof keratin. In vitro antifungal tests against Trichophyton rubrum,Microsporum canis, Fusarium, and Aspergillus fumigatus were cultured on Muller-Hinton agar; the findings indicated that TBH-urea-RLPO IFFS enhanced TBH antifungal activity. Overall, the results support that TBH-urea-RLPO IFFS is an efficient and promising approach for onychomycosis targeting treatment.

Comprehensive kinetic and modeling analyses revealed CYP2C9 and 3A4 determine terbinafine metabolic clearance and bioactivation.

Terbinafine N-dealkylation pathways result in formation of 6,6-dimethyl-2-hepten-4-ynal (TBF-A), a reactive allylic aldehyde, that may initiate idiosyncratic drug-induced liver toxicity. Previously, we reported on the importance of CYP2C19 and 3A4 as major contributors to TBF-A formation. In this study, we expanded on those efforts to assess individual contributions of CYP1A2, 2B6, 2C8, 2C9, and 2D6 in terbinafine metabolism. The combined knowledge gained from these studies allowed us to scale the relative roles of the P450 isozymes in hepatic clearance of terbinafine including pathways leading to TBF-A, and hence, provide a foundation for assessing their significance in terbinafine-induced hepatotoxicity. We used in vitro terbinafine reactions with recombinant P450s to measure kinetics for multiple metabolic pathways and calculated contributions of all individual P450 isozymes to in vivo hepatic clearance for the average human adult. The findings confirmed that CYP3A4 was a major contributor (at least 30% total metabolism) to all three of the possible N-dealkylation pathways; however, CYP2C9, and not CYP2C19, played a critical role in terbinafine metabolism and even exceeded CYP3A4 contributions for terbinafine N-demethylation. A combination of their metabolic capacities accounted for at least 80% of the conversion of terbinafine to TBF-A, while CYP1A2, 2B6, 2C8, and 2D6 made minor contributions. Computational approaches provide a more rapid, less resource-intensive strategy for assessing metabolism, and thus, we additionally predicted terbinafine metabolism using deep neural network models for individual P450 isozymes. Cytochrome P450 isozyme models accurately predicted the likelihood for terbinafine N-demethylation, but overestimated the likelihood for a minor N-denaphthylation pathway. Moreover, the models were not able to differentiate the varying roles of the individual P450 isozymes for specific reactions with this particular drug. Taken together, the significance of CYP2C9 and 3A4 and to a lesser extent, CYP2C19, in terbinafine metabolism is consistent with reported drug interactions. This finding suggests that variations in individual P450 contributions due to other factors like polymorphisms may similarly contribute to terbinafine-related adverse health outcomes. Nevertheless, the impact of their metabolic capacities on formation of reactive TBF-A and consequent idiosyncratic hepatotoxicity will be mitigated by competing detoxification pathways, TBF-A decay, and TBF-A adduction to glutathione that remain understudied.

Lamisil (terbinafine) toxicity: Determining pathways to bioactivation through computational and experimental approaches.

Lamisil (terbinafine) may cause idiosyncratic liver toxicity through a proposed toxicological mechanism involving the reactive metabolite 6,6-dimethyl-2-hepten-4-ynal (TBF-A). TBF-A toxicological relevance remains unclear due to a lack of identification of pathways leading to and competing with TBF-A formation. We resolved this knowledge gap by combining computational modeling and experimental kinetics of in vitro hepatic N-dealkylation of terbinafine. A deep learning model of N-dealkylation predicted a high probability for N-demethylation to yield desmethyl-terbinafine followed by N-dealkylation to TBF-A and marginal contributions from other possible pathways. We carried out steady-state kinetic experiments with pooled human liver microsomes that relied on development of labeling methods to expand metabolite characterization. Those efforts revealed high levels of TBF-A formation and first order decay during metabolic reactions; actual TBF-A levels would then reflect the balance between those processes as well as reflect the impact of stabilizing adduction with glutathione and other biological molecules. Modeling predictions and experimental studies agreed on the significance of N-demethylation and insignificance of N-denaphthylation in terbinafine metabolism, yet differed on importance of direct TBF-A formation. Under steady-state conditions, the direct pathway was the most important source of the reactive metabolite with a Vmax/Km of 4.0 pmol/min/mg protein/µM in contrast to model predictions. Nevertheless, previous studies show that therapeutic dosing leads to accumulation of desmethyl-terbinafine in plasma, which means that likely sources for TBF-A would draw from metabolism of both the major metabolite and parent drug based on our modeling and experimental studies. Through this combination of novel modeling and experimental approaches, we are the first to identify pathways leading to generation of TBF-A for assessing its role in idiosyncratic adverse drug interactions.

Terbinafine hydrochloride nail lacquer for the management of onychomycosis: formulation, characterization and in vitro evaluation.

AIM: The present investigation's intention was to develop an optimized nail lacquer (NL) for the management of onychomycosis. MATERIALS & METHODS: The NL was optimized statistically adopting 32 full factorial design having different polymer ratios and solvent ratios. The formulations were assessed for drug permeation drying time and peak adhesive strength of the film. Characterization was done using techniques including attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR), scanning electron microscopy (SEM), differential scanning calorimetry (DSC), x-ray diffraction (XRD), etc. RESULTS & CONCLUSION: The formulation that had 1:1 polymer ratio and 80:20 solvent ratio was chosen as the optimized formulation. In vitro permeation studies showed better penetration (∼3.25-fold) as well as retention (∼11-fold) of the optimized NL formulation in the animal hoof as compared with the commercial formulation. The findings of in vitro and ex vivo studies elucidated the potential of the optimized formulation. [Formula: see text].

Two strategies to enhance ungual drug permeation from UV-cured films: Incomplete polymerisation to increase drug release and incorporation of chemical enhancers.

UV-curable gels, which polymerise into long-lasting films upon exposure to UVA, have been identified as potential topical drug carriers for the treatment of nail diseases. Limitations of such films include incomplete drug release and low ungual drug permeation. The aim of the work herein was therefore to investigate two strategies, namely: (1) increasing drug release from the film, and (2) increasing nailplate permeability, with the ultimate goal of enhancing ungual drug permeation. To increase drug release via Strategy 1, a UV-LED lamp (whose emitted light was suboptimal for gel polymerisation) was used, and it was hypothesised that such a lamp would result in films that are less polymerised/cross-linked and where the drugs are less 'trapped'. Indeed, the suboptimal lamp influenced polymerisation, such that the films were thinner, had lower glass transition temperatures and enabled a slightly greater (by 15%) drug release of one of the two drugs tested. However, the greater drug release had only a modest impact on ungual drug permeation. To evaluate Strategy 2, i.e. increase nailplate permeability, chemical ungual enhancers, 2-mercaptoethanol (ME), 2-methyl pyrrolidone (NMP), PEG 200 and water were incorporated within the UV-cured films. These chemicals caused increased ungual drug permeation, with ME showing the greatest (by 140%), and water showing the least (by 20%) increase in the amount of drug permeated by day 30. Surprisingly, these chemicals also caused increased drug release from the films, with ME once again having the greatest effect (by 51%) and water the least effect (by 12%). It seems that these chemicals were increasing ungual drug permeation via their influence on drug release (i.e. via their impact on the film) as well as via their influence on the nail itself. We conclude that, of the two strategies tested, the second strategy proved to be more successful at enhancing ungual drug permeation.