Aldehyde Oxidase Mediated Metabolism in Drug-like Molecules: A Combined Computational and Experimental Study

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.

Large Computational Survey of Intrinsic Reactivity of Aromatic Carbon Atoms with Respect to a Model Aldehyde Oxidase

Aldehyde oxidase (AOX) and other related molybdenum-containing enzymes are known to oxidize the C-H bonds of aromatic rings. This process contributes to the metabolism of pharmaceutical compounds and, therefore, is of vital importance to drug pharmacokinetics. The present work describes an automated computational workflow and its use for the prediction of intrinsic reactivity of small aromatic molecules toward a minimal model of the active site of AOX. The workflow is based on quantum chemical transition state searches for the underlying single-step oxidation reaction, where the automated protocol includes identification of unique aromatic C-H bonds, creation of three-dimensional reactant and product complex geometries via a templating approach, search for a transition state, and validation of reaction end points. Conformational search on the reactants, products, and the transition states is performed. The automated procedure has been validated on previously reported transition state barriers and was used to evaluate the intrinsic reactivity of nearly three hundred heterocycles commonly found in approved drug molecules. The intrinsic reactivity of more than 1000 individual aromatic carbon sites is reported. Stereochemical and conformational aspects of the oxidation reaction, which have not been discussed in previous studies, are shown to play important roles in accurate modeling of the oxidation reaction. Observations on structural trends that determine the reactivity are provided and rationalized.

High Affinity and FAP-Targeted Radiotracers: A Potential Design Strategy to Improve the Pharmacokinetics and Tumor Uptake for FAP Inhibitors

Fibroblast activation protein (FAP) is overexpressed in cancer-associated fibroblasts, making it an attractive target for both imaging and therapy of malignancy. This study presents a range of novel FAP inhibitors derived from amino derivatives of UAMC1110, incorporating polyethylene glycol and bulky groups containing bifunctional DOTA chelators. The compounds labeled with gallium-68 were developed and characterized to study biodistribution properties and tumor-targeting performance in nude mice bearing U87MG tumor xenografts. Several tracers of interest were screened due to the advantages in imaging and tumor-specific uptake. Positron emission tomography scans revealed that polyethylene glycol-modified ⁶⁸Ga-3-3 had a rapid penetration within the neoplastic tissue and excellent tumor-to-background contrast. In a comparative biodistribution study, naphthalene-modified ⁶⁸Ga-6-3 exhibited more significant tumor uptake (∼50% ID/g, 1 h p.i.) than ⁶⁸Ga-3-3 and 10-fold higher than ⁶⁸Ga-FAPI-04 under the same conditions. Remarkably, ⁶⁸Ga-8-1, combining the two structural design strategies, obtains superior imaging performance.

In Silico Prediction of Metabolic Reaction Catalyzed by Human Aldehyde Oxidase

Aldehyde oxidase (AOX) plays an important role in drug metabolism. Human AOX (hAOX) is widely distributed in the body, and there are some differences between species. Currently, animal models cannot accurately predict the metabolism of hAOX. Therefore, more and more in silico models have been constructed for the prediction of the hAOX metabolism. These models are based on molecular docking and quantum chemistry theory, which are time-consuming and difficult to automate. Therefore, in this study, we compared traditional machine learning methods, graph convolutional neural network methods, and sequence-based methods with limited data, and proposed a ligand-based model for the metabolism prediction catalyzed by hAOX. Compared with the published models, our model achieved better performance (ACC = 0.91, F1 = 0.77). What's more, we built a web server to predict the sites of metabolism (SOMs) for hAOX. In summary, this study provides a convenient and automatable model and builds a web server named Meta-hAOX for accelerating the drug design and optimization stage.

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

Drug clearance by aldehyde oxidase: can we avoid clinical failure?

Despite increased awareness of aldehyde oxidase (AO) as a major drug-metabolising enzyme, predicting the pharmacokinetics of its substrates remains challenging. Several drug candidates have been terminated due to high clearance, which were subsequently discovered to be AO substrates. Even retrospective extrapolation of human clearance, from models more sensitive to AO activity, often resulted in underprediction.The questions of the current work thus were: Is there an acceptable degree of in vitro AO metabolism that does not result in high in vivo human clearance? And, if so, how can this be predicted?We built an in vitro/in vivo correlation using known AO substrates, combining multiple in vitro parameters to calculate the blood metabolic clearance mediated by AO (CLb_AO). This value was compared with observed blood clearance (CLb_-obs), establishing cut-off CLb_AO values, to discriminate between low and high CLb_-obs. The model was validated using additional literature compounds, and CLb_-obs was predicted in the correct category.This simple, categorical, semi-quantitative yet multi-factorial model is readily applicable in drug discovery. Further, it is valuable for high-clearance compounds, as it predicts the CLb group, rather than an exact CLb value, for the substrates of this poorly-characterised enzyme.

Non-cytochrome P450 enzymes involved in the oxidative metabolism of xenobiotics: Focus on the regulation of gene expression and enzyme activity

Oxidative metabolism is one of the major biotransformation reactions that regulates the exposure of xenobiotics and their metabolites in the circulatory system and local tissues and organs, and influences their efficacy and toxicity. Although cytochrome (CY)P450s play critical roles in the oxidative reaction, extensive CYP450-independent oxidative metabolism also occurs in some xenobiotics, such as aldehyde oxidase, xanthine oxidoreductase, flavin-containing monooxygenase, monoamine oxidase, alcohol dehydrogenase, or aldehyde dehydrogenase-dependent oxidative metabolism. Drugs form a large portion of xenobiotics and are the primary target of this review. The common reaction mechanisms and roles of non-CYP450 enzymes in metabolism, factors affecting the expression and activity of non-CYP450 enzymes in terms of inhibition, induction, regulation, and species differences in pharmaceutical research and development have been summarized. These non-CYP450 enzymes are detoxifying enzymes, although sometimes they mediate severe toxicity. Synthetic or natural chemicals serve as inhibitors for these non-CYP450 enzymes. However, pharmacokinetic-based drug interactions through these inhibitors have rarely been reported in vivo. Although multiple mechanisms participate in the basal expression and regulation of non-CYP450 enzymes, only a limited number of inducers upregulate their expression. Therefore, these enzymes are considered non-inducible or less inducible. Overall, this review focuses on the potential xenobiotic factors that contribute to variations in gene expression levels and the activities of non-CYP450 enzymes.

A comprehensive review of in silico approaches for the prediction and modulation of aldehyde oxidase-mediated drug metabolism: The current features, challenges and future perspectives

The importance of aldehyde oxidase (AOX) in drug metabolism necessitates the development and application of the in silico rational drug design methods as an integral part of drug discovery projects for the early prediction and modulation of AOX-mediated metabolism. The current study represents an up-to-date and thorough review of in silico studies of AOX-mediated metabolism and modulation methods. In addition, the challenges and the knowledge gap that should be covered have been discussed. The importance of aldehyde oxidase (AOX) in drug metabolism is a hot topic in drug discovery. Different strategies are available for the modulation of the AOX-mediated metabolism of drugs. Application of the rational drug design methods as an integral part of drug discovery projects is necessary for the early prediction of AOX-mediated metabolism. The current study represents a comprehensive review of AOX molecular structure, AOX-mediated reactions, AOX substrates, AOX inhibition, approaches to modify AOX-mediated metabolism, prediction of AOX metabolism/substrates/inhibitors, and the AOX related structure-activity relationship (SAR) studies. Furthermore, an up-to-date and thorough review of in silico studies of AOX metabolism has been carried out. In addition, the challenges and the knowledge gap that should be covered in the scientific literature have been discussed in the current review.

Recent developments in predicting CYP-independent metabolism

As lead optimization efforts have successfully reduced metabolic liabilities due to cytochrome P450 (CYP)-mediated metabolism, there has been an increase in the frequency of involvement of non-CYP enzymes in the metabolism of investigational compounds. Although there have been numerous notable advancements in the characterization of non-CYP enzymes with respect to their localization, reaction mechanisms, species differences and identification of typical substrates, accurate prediction of non-CYP-mediated clearance, with a particular emphasis with the difficulties in accounting for any extrahepatic contributions, remains a challenge. The current manuscript comprehensively summarizes the recent advancements in the prediction of drug metabolism and the in vitro to in vitro extrapolation of clearance for substrates of non-CYP drug metabolizing enzymes.

Organophotocatalytic selective deuterodehalogenation of aryl or alkyl chlorides

Development of practical deuteration reactions is highly valuable for organic synthesis, analytic chemistry and pharmaceutic chemistry. Deuterodehalogenation of organic chlorides tends to be an attractive strategy but remains a challenging task. We here develop a photocatalytic system consisting of an aryl-amine photocatalyst and a disulfide co-catalyst in the presence of sodium formate as an electron and hydrogen donor. Accordingly, many aryl chlorides, alkyl chlorides, and other halides are converted to deuterated products at room temperature in air (>90 examples, up to 99% D-incorporation). The mechanistic studies reveal that the aryl amine serves as reducing photoredox catalyst to initiate cleavage of the C-Cl bond, at the same time as energy transfer catalyst to induce homolysis of the disulfide for consequent deuterium transfer process. This economic and environmentally-friendly method can be used for site-selective D-labeling of a number of bioactive molecules and direct H/D exchange of some drug molecules.

Deuterodehalogenation of organic chlorides is a useful strategy to install deuterium atoms at specific positions, however, it has several drawbacks. In this study, the authors report an organophotocatalytic system consisting of an aryl-amine-based photocatalyst and a common disulfide co-catalyst, for efficient deuteration of a wide range of aryl chlorides, alkyl chlorides and other halides, at room temperature in air.

Nimesulide increases the aldehyde oxidase activity of humans and rats

An increasing number of drugs are metabolized by aldehyde oxidase (AOX), but AOX-mediated drug interactions are seldom reported due to the lack of appropriate inhibitors and inducers. A recent study reported that nimesulide (NIM) could increase the liver injury risk of methotrexate. The latter was mainly metabolized by AOX to form hepatotoxic 7-hydroxymethotrexate (7-OH MTX). Thus, we speculated that NIM could induce AOX. In this study, we investigated the potential induction of AOX activity by NIM using methotrexate as the probe substrate. Treatment of primary human and rat hepatocytes with NIM (20 μM) for 24 h caused a 2.0- and 3.1-fold, respectively, increase in 7-OH MTX formation. Oral administration of NIM (100 mg·kg⁻¹·d⁻¹, for 5 days) to rats significantly increased the systematic exposure (6.5-fold), liver distribution (2.5-fold), and excretion (5.2-fold for urinary excretion and 2.1-fold for fecal excretion) of 7-OH MTX. The 7-OH MTX formation in liver cytosol from rats pretreated with 20, 50, and 100 mg·kg⁻¹·d⁻¹ NIM for 5 days increased by 1.9-, 3.2-, and 3.7-fold, respectively, compared with that of rats pretreated with the vehicle. We revealed that the elevation of AOX activity was accompanied by an increase in AOX1 protein levels but not the corresponding mRNA levels. Collectively, our results demonstrate for the first time that NIM can increase the AOX activity of humans and rats, and may raise concerns regarding the risk of drug interactions between NIM and AOX substrates in clinical practice.

Predicting liver cytosol stability of small molecules

Over the last few decades, chemists have become skilled at designing compounds that avoid cytochrome P (CYP) 450 mediated metabolism. Typical screening assays are performed in liver microsomal fractions and it is possible to overlook the contribution of cytosolic enzymes until much later in the drug discovery process. Few data exist on cytosolic enzyme-mediated metabolism and no reliable tools are available to chemists to help design away from such liabilities. In this study, we screened 1450 compounds for liver cytosol-mediated metabolic stability and extracted transformation rules that might help medicinal chemists in optimizing compounds with these liabilities. In vitro half-life data were collected by performing in-house experiments in mouse (CD-1 male) and human (mixed gender) cytosol fractions. Matched molecular pairs analysis was performed in conjunction with qualitative-structure activity relationship modeling to identify chemical structure transformations affecting cytosolic stability. The transformation rules were prospectively validated on the test set. In addition, selected rules were validated on a diverse chemical library and the resulting pairs were experimentally tested to confirm whether the identified transformations could be generalized. The validation results, comprising nearly 250 library compounds and corresponding half-life data, are made publicly available. The datasets were also used to generate in silico classification models, based on different molecular descriptors and machine learning methods, to predict cytosol-mediated liabilities. To the best of our knowledge, this is the first systematic in silico effort to address cytosolic enzyme-mediated liabilities.

New Preclinical Development of a c-Met Inhibitor and Its Combined Anti-Tumor Effect in c-Met-Amplified NSCLC

c-Met is a receptor tyrosine kinase with no commercially available product despite being a pivotal target in cancer progression. Unlike other c-Met inhibitors that fail clinically, ABN401 is a newly synthesized c-Met inhibitor that is not potentially degraded by aldehyde oxidase (AO) in human liver cytosol. This study aimed to determine the physicochemical stability, pharmacokinetics in beagle dogs, and therapeutic effect of ABN401 in a c-Met-amplified non-small-cell lung cancer (NSCLC) patient-derived xenograft (PDX) model. ABN401 was found to be a weak basic compound, with pKa and log P values of 7.49 and 2.46, respectively. It is poorly water-soluble but soluble at acidic pH. The accelerated storage stability is dependent on temperature, but the purity remains at over 97% after 6 months. The bioavailability is approximately 30% in dogs and it is highly efficient in the PDX model, achieving around 90% tumor growth inhibition in combination with erlotinib. These observations indicate that the compound is acceptable for the next phase of trials.

Metabolism by Aldehyde Oxidase: Drug Design and Complementary Approaches to Challenges in Drug Discovery

Aldehyde oxidase (AO) catalyzes oxidations of azaheterocycles and aldehydes, amide hydrolysis, and diverse reductions. AO substrates are rare among marketed drugs, and many candidates failed due to poor pharmacokinetics, interspecies differences, and adverse effects. As most issues arise from complex and poorly understood AO biology, an effective solution is to stop or decrease AO metabolism. This perspective focuses on rational drug design approaches to modulate AO-mediated metabolism in drug discovery. AO biological aspects are also covered, as they are complementary to chemical design and important when selecting the experimental system for risk assessment. The authors' recommendation is an early consideration of AO-mediated metabolism supported by computational and in vitro experimental methods but not an automatic avoidance of AO structural flags, many of which are versatile and valuable building blocks. Preferably, consideration of AO-mediated metabolism should be part of the multiparametric drug optimization process, with the goal to improve overall drug-like properties.

Scaffold-hopping as a strategy to address metabolic liabilities of aromatic compounds

Understanding and minimizing oxidative metabolism of aromatic compounds is a key hurdle in lead optimization. Metabolic processes not only clear compounds from the body, but they can also transform parent compounds into reactive metabolites. One particularly useful strategy when addressing metabolically labile or oxidation-prone structures is scaffold-hopping. Replacement of an aromatic system with a more electron-deficient ring system can often increase robustness towards cytochrome P450-mediated oxidation while conserving the structural requirements of the pharmacophore. The most common example of this substitution strategy, replacement of a phenyl ring with a pyridyl substituent, is prevalent throughout the literature; however scaffold-hopping encompasses a much wider scope of heterocycle replacement. This review will showcase recent examples where different scaffold-hopping approaches were used to reduce metabolic clearance or block the formation of reactive metabolites. Additionally, we will highlight considerations that should be made to garner the most benefit from a scaffold-hopping strategy for lead optimization.

Aldehyde oxidase and its role as a drug metabolizing enzyme

Aldehyde oxidase (AO) is a cytosolic enzyme that belongs to the family of structurally related molybdoflavoproteins like xanthine oxidase (XO). The enzyme is characterized by broad substrate specificity and marked species differences. It catalyzes the oxidation of aromatic and aliphatic aldehydes and various heteroaromatic rings as well as reduction of several functional groups. The references to AO and its role in metabolism date back to the 1950s, but the importance of this enzyme in the metabolism of drugs has emerged in the past fifteen years. Several reviews on the role of AO in drug metabolism have been published in the past decade indicative of the growing interest in the enzyme and its influence in drug metabolism. Here, we present a comprehensive monograph of AO as a drug metabolizing enzyme with emphasis on marketed drugs as well as other xenobiotics, as substrates and inhibitors. Although the number of drugs that are primarily metabolized by AO are few, the impact of AO on drug development has been extensive. We also discuss the effect of AO on the systemic exposure and clearance these clinical candidates. The review provides a comprehensive analysis of drug discovery compounds involving AO with the focus on developmental candidates that were reported in the past five years with regards to pharmacokinetics and toxicity. While there is only one known report of AO-mediated clinically relevant drug-drug interaction (DDI), a detailed description of inhibitors and inducers of AO known to date has been presented here and the potential risks associated with DDI. The increasing recognition of the importance of AO has led to significant progress in predicting the site of AO-mediated metabolism using computational methods. Additionally, marked species difference in expression of AO makes it is difficult to predict human clearance with high confidence. The progress made towards developing in vivo, in vitro and in silico approaches for predicting AO metabolism and estimating human clearance of compounds that are metabolized by AO have also been discussed.

Fast Methods for Prediction of Aldehyde Oxidase-Mediated Site-of-Metabolism

Aldehyde Oxidase (AO) is an enzyme involved in the metabolism of aldehydes and N-containing heterocyclic compounds. Many drug compounds contain heterocyclic moieties, and AO metabolism has lead to failure of several late-stage drug candidates. Therefore, it is important to take AO-mediated metabolism into account early in the drug discovery process, and thus, to have fast and reliable models to predict the site of metabolism (SOM). We have collected a dataset of 78 substrates of human AO with a total of 89 SOMs and 347 non-SOMs and determined atomic descriptors for each compound. The descriptors comprise NMR shielding and ESP charges from density functional theory (DFT), NMR chemical shift from ChemBioDraw, and Gasteiger charges from RDKit. Additionally, atomic accessibility was considered using 2D-SASA and relative span descriptors from SMARTCyp. Finally, stability of the product, the metabolite, was determined with DFT and also used as a descriptor. All descriptors have AUC larger than 0.75. In particular, descriptors related to the chemical shielding and chemical shift (AUC = 0.96) and ESP charges (AUC = 0.96) proved to be good descriptors. We recommend two simple methods to identify the SOM for a given molecule: 1) use ChemBioDraw to calculate the chemical shift or 2) calculate ESP charges or chemical shift using DFT. The first approach is fast but somewhat difficult to automate, while the second is more time-consuming, but can easily be automated. The two methods predict correctly 93% and 91%, respectively, of the 89 experimentally observed SOMs.

Metabolism of c-Met Kinase Inhibitors Containing Quinoline by Aldehyde Oxidase, Electron Donating, and Steric Hindrance Effect

Some quinoline-containing c-Met kinase inhibitors are aldehyde oxidase (AO) substrates. 3-Substituted quinoline triazolopyridine analogs were synthesized to understand the electron-donating and steric hindrance effects on AO-mediated metabolism. Metabolic stability studies for these quinoline analogs were carried out in liver cytosol from mice, rats, cynomolgus monkeys, and humans. Several 3-N-substituted analogs were found to be unstable in monkey liver cytosolic incubations (half-life, <10 minutes), and five of them (63, 53, 51, 11, and 71) were chosen for additional mechanistic studies. Mono-oxygenation on the quinoline ring was identified by liquid chromatography tandem mass spectrometry. Metabolite formation was inhibited by the AO inhibitors menadione and raloxifene, but not by the xanthine oxidase inhibitor allopurinol. It was found that small electron-donating groups at the 3-quinoline moiety made the analogs more susceptible to AO metabolism, whereas large 3-substituents could reverse the trend. Although species differences were observed, this trend was applicable to all species tested. Small electron-donating substituents at the 3-quinoline moiety increased both affinity (decreased Michaelis constant) and V _max maximum velocity toward AO in kinetic studies, whereas large substituents decreased both parameters probably as a result of steric hindrance. Based on our analysis, a common structural feature with high AO liability was proposed. Our finding could provide useful information for chemists to minimize potential AO liability when designing quinoline analogs.

A General Strategy for Site-Selective Incorporation of Deuterium and Tritium into Pyridines, Diazines, and Pharmaceuticals

Methods to incorporate deuterium and tritium atoms into organic molecules are valuable for medicinal chemistry. The prevalence of pyridines and diazines in pharmaceuticals means that new ways to label these heterocycles will present opportunities in drug design and facilitate absorption, distribution, metabolism, and excretion (ADME) studies. A broadly applicable protocol is presented wherein pyridines, diazines, and pharmaceuticals are converted into heterocyclic phosphonium salts and then isotopically labeled. The isotopes are incorporated in high yields and, in general, with exclusive regioselectivity.

Direct comparison of the four aldehyde oxidase enzymes present in mouse gives insight into their substrate specificities

Mammalian aldehyde oxidases (AOXs) are molybdo-flavoenzymes which are present in many tissues in various mammalian species, including humans and rodents. Different species contain a different number of AOX isoforms. In particular, the reasons why mammals other than humans express a multiplicity of tissue-specific AOX enzymes is unknown. In mouse, the isoforms mAOX1, mAOX3, mAOX4 and mAOX2 are present. We previously established a codon-optimized heterologous expression systems for the mAOX1-4 isoforms in Escherichia coli that gives yield to sufficient amounts of active protein for kinetic characterizations and sets the basis in this study for site-directed mutagenesis and structure-function studies. A direct and simultaneous comparison of the enzymatic properties and characteristics of the four enzymes on a larger number of substrates has never been performed. Here, thirty different structurally related aromatic, aliphatic and N-heterocyclic compounds were used as substrates, and the kinetic parameters of all four mAOX enzymes were directly compared. The results show that especially mAOX4 displays a higher substrate selectivity, while no major differences between mAOX1, mAOX2 and mAOX3 were identified. Generally, mAOX1 was the enzyme with the highest catalytic turnover for most substrates. To understand the factors that contribute to the substrate specificity of mAOX4, site-directed mutagenesis was applied to substitute amino acids in the substrate-binding funnel by the ones present in mAOX1, mAOX3, and mAOX2. An increase in activity was obtained by the amino acid exchange M1088V in the active site identified to be specific for mAOX4, to the amino acid identified in mAOX3.

Aldehyde Oxidase: Reaction Mechanism and Prediction of Site of Metabolism

Aldehyde oxidase (AO) is a molybdenum-containing enzyme involved in the clearance of drug compounds containing aldehydes and N-containing heterocyclic fragments. AO has gained considerable interest in recent years because of examples of too fast clearance of drug compounds in development. Thus, it is important to be able to predict AO-mediated drug metabolism. Therefore, we have characterized the structural and energetic aspects of different mechanisms with density functional theory using the molybdenum cofactor as a model for the reactive part of the enzyme. For a series of 6-substituted 4-quinazolinones, the trend in activation energies is the same for three tested reaction mechanisms. Using the concerted mechanism as a model for the enzymatic reaction, the transition states (TSs) for the formation of all possible metabolites for a series of known AO substrates were determined. The lowest activation energies correspond in all cases to the experimentally observed sites of metabolism (SOMs). Various molecular properties were calculated and investigated as more easily determinable markers for reactivity. The stabilities of both intermediates and products correlate to some extent with the TS energies and may be used to predict the SOM. The electrostatic-potential-derived charges are also good markers for the prediction of the experimental SOM for this set of compounds and may pave the way for the development of fast methods for the prediction of SOM for AO substrates.