ADMET Profiler

From ACD Percepta
Revision as of 14:50, 20 December 2013 by Kirilas (talk | contribs) (Added technical section with categorization function info)
Jump to navigation Jump to search

Overview

Drug Profiler is an integrated workspace that aggregates the output of various PhysChem, ADME, and Drug safety-related predictors included in ACD/Percepta and presents the data as the compound's overall ADME/Tox profile with individual predictions categorized according to user-adjustable thresholds and visualized using color-coded indications. Drug Profiler enables the researcher to identify the most suitable structures that meet specific criteria for a wide range of properties:

  • PhysChem
    • LogP
    • Solubility
    • H-bonding
    • Molecular size and flexibility
    • Drug- and Lead-likeness rules
  • ADME
    • Human Intestinal Absorption (HIA) and Caco-2 permeability
    • Plasma protein binding (PPB)
    • Blood-brain pernetration and accessibility to CNS
    • Metabolic stability
  • Drug Safety
    • P-gp substrate specificity
    • Cytochrome P450 inhibitor specificity
    • hERG inhibitor specificity
    • Mutagenicity in Ames test

Interface

drug profiler.png


  1. Navigation Bar
    The Navigation Bar appears in all workspaces except the Spreadsheet View. It displays presented structure and two previous as well as two next stuctures in the corresponding dataset. If no dataset is selected, the Navigation Bar displays the history of performed calculations.
    When the compound outside of the dataset is used as an input structure for calculation, the user can add it to the current dataset:

    Navigation bar.png

    1. Click the History confirm icon.png pictogram to replace the current structure with the new one. All subsequent calculations will be performed automatically for the added compound.
    2. Click the History add compound icon.png pictogram to add new compound to the current dataset. It will appear as the last record in the dataset and all subsequent calculations will be performed automatically.
    3. The structure browse bar displays the active structure number and the total number of records in the current dataset. Scroll through the database using the browse buttons: First/Last record, Previous/Next record or pressing the up/down arrow keys on the keyboard to go to the next/previous record. Type a record number to view it directly.

  2. Structure Pane
    The Structure Pane displays a 2-D representation of the chemical structure used for calculation. A compound name is displayed only for those structures which are found in the Dictionary.
    Paste.png (Ctrl+V): Pastes an input structure from the clipboard to the Structure Pane. Also, if you are pasting a SMILES string from the clipboard, you can save a step by clicking directly on this button to display the 2-D structure. The calculation will be performed automatically for the input structure.
    Copy pictogram.png (Ctrl+C): Copies the active structure in the Structure Pane to the clipboard.
    Opendictionary.png (Ctrl+D): Opens the Dictionary with >117,000 available names to find the corresponding structure and display it. When a compound is selected from the list and the selection confirmed by clicking OK, calculations for the chosen compound will be performed automatically.
    SMILES pictogram.png (Ctrl+S): Type or paste SMILES notation string in this field. The string will automatically be converted to 2-D structure. Calculation will be performed automatically for the input structure.
    Open pictogram.png (Ctrl+O): Opens a file containing a molecular structure. ACD/Percepta supports all versions of MDL *.mol files, ACD/ChemSketch *.sk2, as well as ISIS Draw *.skc, and ChemDraw *.cdx formats.
    Edit pictogram.png (Ctrl+E): Allows to edit the active molecule.

  3. Create pdf pictogram.png Creates a PDF report with the displayed results.
    Copy results pictogram.png Copies the displayed results to the clipboard.

  4. The dropdown list of available threshold templates. Use default thresholds or select any of the previously saved templates from the dropdown, and the property categories will be automatically re-assigned to match the newly selected set of thresholds. If any of the thresholds are altered compared to last saved version, this is indicated by an asterisk displayed next to the name of the template. The changes can be saved under the same template name, by pressing Save button, or alternatively, transfered to a new template by choosing Save as... command.

  5. A list of properties included in the profile, their calculated values, and category indication with a short textual description. The property list is interactive - click the property of interest to display its value scale and view or edit the cut-offs separating the categories in the lower left part of the window.

  6. Displays the value scale of corresponding property and the list of cut-offs defining different categories. A black vertical line indicates the position of the current compound on the property scale. Move the sliders to alter the numerical cut-offs separating the categories.

  7. Hover over the property value to view a screen-tip with more details regarding the particular prediction. For some properties, the screen-tips display additional quantitative parameters of relevance (e.g., CNS), predicted probabilities of positive assay outcome in case of probabilistic properties, and/or the numerical conditions that the compound had fulfilled in obtain to be classified in a particular category.

Technical information

All categorization functions in Drug Profiler can be divided into two large groups:

  • Based on continuous properties
  • Based on probabilistic predictors

Continuous properties

For continuous properties the categories are defined by setting-up certain numerical thresholds on the scale of the corresponding property values. The following is a complete list of default thresholds along with additional comments where appropriate.

PhysChem Profiling

The majority of the PhysChem related categorization functions employ upper property value thresholds introduced by Christopher A. Lipinski in his "Rule of 5". In some cases additional cut-offs on the lower end of the scale have been introduced as well as certain intermediate ranges resulting in the Average outcome. The latter are mostly associated with the concept of lead-likeness which is described later in this text.

Lipophilicity

Value range Category
logP <= -2 Very hydrophilic
-2 < logP <= -1 Hydrophilic
-1 < logP <= 4.2 Optimal
4.2 < logP <= 5 Hydrophobic
logP > 5 Very hydrophobic

Molecular weight

Value range Category
MW <= 460 Good
460 < MW <= 500 Moderate
MW > 500 Bad

Hydrogen bonding acidity

Value range Category
No. of H-Bond donors <= 5 Good
No. of H-Bond donors > 5 Bad

Hydrogen bonding basicity

Value range Category
No. of H-Bond acceptors <= 10 Good
No. of H-Bond acceptors > 10 Bad

"Rule of 5"
The Lipinski categorization function acts as a sort of a summary of the above separate categories, by testing a compound against all rules at once (-2 =< logP < 5; MW < 500; No. of H-Bond donors <= 5; No. of H-Bond acceptors <= 10). Category assignments are the following:

No. of rule violations Category
0 (All 4 rules apply) Good
1..2 (Any 2 or 3 rules apply) Moderate
>2 (Less than 2 rules apply) Bad

Lead-likeness
The Lead-like category assignment functions in exactly the same way, however a little bit lower property value thresholds are usually employed. This concept is based on the assumption that the lead compound undergoes various modifications during the development stage and establishment of certain trends in which way those modifications actually affect the physicochemical properties of the lead compound. In other words, it is an attempt to establish the rules for filtering out the compounds that are unlikely to violate Lipinski’s "Rule of 5" after the development stage (-2 =< logP < 4.2; MW < 460; No. of H-Bond donors <= 5; No. of H-Bond acceptors <= 9).

There are also a few additional categorization functions in the PhysChem Profiling group. Flexibility of the molecule

Value range Category
No. of rotatable bonds <= 10 Good
No. of rotatable bonds > 10 Bad

Number of Rings

Value range Category
No. of Rings <= 4 Good
No. of Rings > 4 Bad

Solubility
The classification scheme in case of the compound solubility is a slightly modified version of the Traffic Lights system published in Lobell, M., et al. Chem. Med. Chem., 2006, 1, 1229-1236. It takes into account the solubility of the compound in buffer at pH of 6.5.

Solubility at pH=6.5 Category
S > 100 mg/L Soluble
10 mg/L < S <= 100 mg/L Insoluble
S <= 10 mg/L Highly insoluble

ADME Profiling

Plasma protein binding
Categorization function for the extent of plasma protein binding takes into account both predicted %PPB values and Reliability Index (RI) of predictions:

Value range Category
Any value with RI < 0.3 Undefined
%PPB <= 10% Not bound
10% < %PPB <= 40% Weakly bound
40% < %PPB <= 80% Moderately bound
80% < %PPB <= 90% Strongly bound
%PPB > 90% Extensively bound

Caco-2 Permeability
Category assignment is based on predicted permeability across Caco-2 monolayers (Pe) at pH = 7.4 and 500 rpm stirring rate:

Value range Category
Pe <= 1 * 10-6 cm/s Poorly permeable
1 * 10-6 cm/s < Pe <= 7 * 10-6 cm/s Moderately permeable
Pe > 7 * 10-6 cm/s Highly permeable

Human Intestinal Absorption

Value range Category
%HIA <= 30% Poorly permeable
30% < %HIA <= 70% Moderately permeable
%HIA > 70% Highly permeable

CNS access
Evaluation of compounds’ brain penetration potential and their classification by CNS access is based on predicted quantitative characteristics of transport across BBB, namely brain/plasma equilibration rate, Log (PS x fu,brain) and steady-state brain/plasma distribution ratio, LogBB. The categories are assigned according to the values of CNS Access Score which is calculated as a simple combination of those parameters: Score = Log(PS x fu,brain) + LogBB.

Score range Category Comment
Score <= -3.5 Non-penetrant Brain delivery too low for a compound to be active in CNS
-3.5 < Score <= -3 Weak penetrant Borderline brain delivery potential
Score > -3 Penetrant Compound sufficiently accessible to CNS to exhibit pharmacological effect in brain

The above mentioned parameters reflect compounds' transport across BBB governed by passive diffusion only. Therefore, CNS Access Score may also be shifted up or down if the respective compound is predicted to undergo carrier-mediated transport at the BBB. The Score is increased by 0.5 units for probable substrates of influx transporters, and decreased by 1 unit for compounds potentially susceptible to P-gp efflux. These are the final Score values used in category assignment.

Probabilistic categories

In probabilistic categorization functions the numerical predictions are first converted to the so called classification Scores on the basis of the obtained probability (p) and Reliability Index (RI) values. The conversion is performed according to the following equation: Score = (p - 0.5) * RI + 0.5

The main idea lying behind this transformation is that predicted probability close to 0.5 and low RI are both indicators of an inconclusive prediction. Therefore, the amount by which the initially calculated probability differs from 0.5 is multiplied by RI, so that predictions with low RI values ultimately yield scores in the intermediate range even if the original p-value is quite high or low. The final category assignment takes place on the basis of the defined classification score ranges.

This class of probabilistic categorization functions covers the Metabolic Stability category in the ADME Profiling group and the entire Drug Safety Profiling group:

  • Ames categories
  • hERG categories
  • Pgp Substrates categories
  • CYP450 Inhibitor categories (3A4, 2D6, 2C9, 2C19, 1A2)

In case of AMES, hERG and Metabolic Stability categories everything applies as explained. Since the classification criteria for all of these properties are such that a positive outcome is considered unfavorable for a drug-candidate (i.e., Ames positive, hERG inhibitor with IC50 < 10 uM, and t1/2 in Human Liver Microsomes < 30 min., respectively), high probabilities result in a classification as Bad, i.e.:

Score range Category label
Ames hERG HLM
Score <= 0.33 Non-mutagenic Non-inhibitor Stable in HLM
0.33 < Score <= 0.67 Undefined Undefined Undefined
Score > 0.67 Mutagenic Inhibitor Unstable in HLM

P-gp Substrate and Cytochrome P450 inhibition categories constitute special cases as they are assigned on the basis of two predictions each. In case of P-gp:

  • Score1 is derived from Good P-gp substrate model (compounds are only considered positive if they undergo significant P-gp efflux in vivo)
  • Score2 is derived from General P-gp substrate model

In case of all considered cytochrome P450 isoforms, the two models correspond to different IC50 thresholds for classification:

  • Score1 - Efficient inhibition model (IC50 < 10 uM)
  • Score2 - General inhibition model (IC50 < 50 uM)

The categories are assigned as follows:

Score1 Score2 Category label
P-gp CYP450
Score1 > 0.67 Score2 - not relevant Good substrate Efficient inhibitor
Score1 <= 0.67 Score2 > 0.67 Substrate Inhibitor
Score1 <= 0.67 0.33 < Score2 <= 0.67 Undefined Undefined
Score1 <= 0.67 Score2 <= 0.33 Non-substrate Non-inhibitor