Drug-likeness

Check a structure against the Lipinski, Veber, Ghose, Egan and lead-likeness rules, with QED, a bioavailability radar and every descriptor behind the verdict.

Compound

Examples:

Structure

NHONNNHNNN
ImatinibC29H31N7O · 493.62 g/mol

raises logPlowers logP

Rule sets

Passes 3 of 5
  • Lipinski rule of fivePasses

    All 4 criteria met.

    Lipinski et al., Adv. Drug Deliv. Rev. 1997, 23, 3

  • VeberPasses

    All 2 criteria met.

    Veber et al., J. Med. Chem. 2002, 45, 2615

  • GhoseFails
    • Molar refractivity 146.89 cm³/mol, above 130
    • Molar mass 493.62 g/mol, above 480

    Ghose et al., J. Comb. Chem. 1999, 1, 55

  • EganPasses

    All 2 criteria met.

    Egan et al., J. Med. Chem. 2000, 43, 3867

  • Lead-likenessFails
    • Molar mass 493.62 g/mol, above 350
    • logP 4.59, above 3.5

    Teague et al., Angew. Chem. Int. Ed. 1999, 38, 3743

  • QED0.39

    Quantitative estimate of drug-likeness from 0 to 1: the weighted geometric mean of eight property desirabilities.

    Bickerton et al., Nat. Chem. 2012, 4, 90

Bioavailability radar

Lipophilicity: logP 4.59Size: Molar mass 493.62Polarity: TPSA 86.28Insolubility: logS (ESOL) −5.81Unsaturation: Fraction Csp³ 0.24Flexibility: Rotatable bonds 7 LipophilicitySizePolarityInsolubilityUnsaturationFlexibility

Shaded: the range suited to oral bioavailability. Red points fall outside it.

AxisValueOptimal
LipophilicitylogP4.59−0.7 to 5
SizeMolar mass493.62 g/mol150 to 500
PolarityTPSA86.28 Ų20 to 130
InsolubilitylogS (ESOL)−5.81−6 to 0
UnsaturationFraction Csp³0.240.25 to 1
FlexibilityRotatable bonds70 to 9

Descriptors

PropertyValueLimits
Molar mass493.62 g/molLipinski ≤ 500Ghose 160 to 480Lead-like 250 to 350
logPWildman–Crippen4.59Lipinski ≤ 5Ghose −0.4 to 5.6Egan ≤ 5.88Lead-like ≤ 3.5
TPSAtopological polar surface area, N and O86.28 ŲVeber ≤ 140Egan ≤ 131.6
NH and OHLipinski's donor count2Lipinski ≤ 5
N and OLipinski's acceptor count8Lipinski ≤ 10
H-bond donors2—
H-bond acceptors7—
Rotatable bonds7Veber ≤ 10Lead-like ≤ 7
Molar refractivityWildman–Crippen146.89 cm³/molGhose 40 to 130
Atomshydrogens included68Ghose 20 to 70
Heavy atoms37—
Fraction Csp³0.24—
Aromatic rings4—
Ringssmallest set of smallest rings5—
Formal charge0—
logSESOL estimate, log(mol/L)−5.81—

Details

logP and molar refractivity by atom

Each atom and its hydrogens get one of the Wildman–Crippen atom types; logP and molar refractivity are sums over the types. Point at a row to find the atom in the structure.

AtomTypeEnvironmentlogPMR
C1+3 HC8H1methyl on an aromatic carbon+0.4545.635
C2C21aromatic carbon bearing an alkyl group+0.1363.509
C3+1 HC18H1aromatic CH+0.2814.407
C4+1 HC18H1aromatic CH+0.2814.407
C5C22aromatic carbon bearing nitrogen+0.4624.067
N6+1 HN4H3secondary aromatic amine−0.3053.963
C7C5carbon double-bonded to a heteroatom−0.2785.007
O8O10aryl carbonyl oxygen+0.1130.222
C9C21aromatic carbon bearing an alkyl group+0.1363.509
C10+1 HC18H1aromatic CH+0.2814.407
C11+1 HC18H1aromatic CH+0.2814.407
C12C21aromatic carbon bearing an alkyl group+0.1363.509
C13+2 HC10H1benzylic CH2+0.1944.602
N14N7tertiary amine−0.3191.839
C15+2 HC3H1aliphatic CH3 or CH2 bonded to a heteroatom+0.0434.867
C16+2 HC3H1aliphatic CH3 or CH2 bonded to a heteroatom+0.0434.867
N17N7tertiary amine−0.3191.839
C18+3 HC3H1aliphatic CH3 or CH2 bonded to a heteroatom+0.1665.924
C19+2 HC3H1aliphatic CH3 or CH2 bonded to a heteroatom+0.0434.867
C20+2 HC3H1aliphatic CH3 or CH2 bonded to a heteroatom+0.0434.867
C21+1 HC18H1aromatic CH+0.2814.407
C22+1 HC18H1aromatic CH+0.2814.407
C23+1 HC18H1aromatic CH+0.2814.407
C24C22aromatic carbon bearing nitrogen+0.4624.067
N25+1 HN4H3secondary aromatic amine−0.3053.963
C26C22aromatic carbon bearing nitrogen+0.4624.067
N27N11aromatic nitrogen−0.3242.202
C28+1 HC18H1aromatic CH+0.2814.407
C29+1 HC18H1aromatic CH+0.2814.407
C30C20aromatic carbon of a biaryl bond+0.2713.904
C31C20aromatic carbon of a biaryl bond+0.2713.904
C32+1 HC18H1aromatic CH+0.2814.407
C33+1 HC18H1aromatic CH+0.2814.407
C34+1 HC18H1aromatic CH+0.2814.407
N35N11aromatic nitrogen−0.3242.202
C36+1 HC18H1aromatic CH+0.2814.407
N37N11aromatic nitrogen−0.3242.202
Sum4.590146.894
Other definitions and the QED breakdown
H-bond acceptors, RDKit NumHAcceptors
7
H-bond acceptors, RDKit Lipinski module SMARTS
7
N and O, Lipinski
8
H-bond donors, RDKit NumHDonors
2
NH and OH, Lipinski
2
TPSA with S and P
86.28 Ų
Exact mass
493.2590 Da
Components
1

QED, property by property

PropertyValueDesirability
Molar mass493.620.188
logP4.590.652
H-bond acceptors (QED definition)70.332
H-bond donors20.792
TPSA86.280.821
Rotatable bonds70.533
Aromatic rings (QED definition)40.035
Structural alerts00.842
How the descriptors are defined

logP and molar refractivity follow Wildman and Crippen (J. Chem. Inf. Comput. Sci. 1999, 39, 868): every atom, hydrogens included, is assigned one of 72 atom types, defined by 110 SMARTS patterns, and the property is the sum of the type contributions. The type table and its order are RDKit's, so the numbers match RDKit's MolLogP and MolMR.

TPSA is Ertl's topological polar surface area (J. Med. Chem. 2000, 43, 3714) from nitrogen and oxygen only, as in Ertl's own implementation. The value with sulfur and phosphorus is listed under Details.

H-bond donors and acceptors are RDKit's definitions: donors are NH, OH, SH and aromatic NH; acceptors are N, O and S atoms with a free lone pair, which leaves out carboxylic acid OH groups, amide-type nitrogens and aromatic nitrogens that carry a hydrogen or a substituent, as in pyrrole or indole. Lipinski counted every NH and OH as a donor and every N and O as an acceptor, and the rule of five below uses those counts.

Rotatable bonds are acyclic single bonds between non-terminal heavy atoms, excluding bonds to triple-bonded atoms, amide, thioamide and amidine C–N bonds, and bonds to CF3, CCl3, CBr3, tert-butyl and methyl groups (RDKit's strict definition).

Across the 31,358 compounds of our catalogue, logP, molar refractivity, TPSA, donors, acceptors, rotatable bonds and fraction Csp³ agree with RDKit for more than 99.9 % of structures; the exceptions are aromaticity edge cases such as fused N-oxides and large polycyclic aromatics.

What the rule sets say

Lipinski (Adv. Drug Deliv. Rev. 1997, 23, 3): poor absorption or permeation is more likely when more than one of these holds: molar mass above 500, logP above 5, more than 5 NH and OH, more than 10 N and O. Lipinski used ClogP; this page uses Wildman–Crippen logP against the same limit.

Veber (J. Med. Chem. 2002, 45, 2615): compounds with 10 or fewer rotatable bonds and a polar surface area of 140 Ų or less were likely to show good oral bioavailability in the rat.

Ghose (J. Comb. Chem. 1999, 1, 55): the qualifying range covering 80 % of known drugs, logP −0.4 to 5.6, molar refractivity 40 to 130, molar mass 160 to 480 and 20 to 70 atoms.

Egan (J. Med. Chem. 2000, 43, 3867): the edges of the absorption egg, logP 5.88 and TPSA 131.6 Ų.

Lead-likeness (Teague, Davis, Leeson and Oprea, Angew. Chem. Int. Ed. 1999, 38, 3743): leads should leave room for optimisation; molar mass 250 to 350, logP up to 3.5 and at most 7 rotatable bonds, the limits SwissADME applies.

QED (Bickerton et al., Nat. Chem. 2012, 4, 90) maps molar mass, logP, acceptors, donors, polar surface area, rotatable bonds, aromatic rings and structural alerts onto desirability functions and combines them with the published mean weights.

How the radar is drawn

The radar follows SwissADME (Daina, Michielin and Zoete, Sci. Rep. 2017, 7, 42717). The shaded band on each axis is the range considered suitable for oral bioavailability: logP −0.7 to 5.0, molar mass 150 to 500 g/mol, TPSA 20 to 130 Ų, logS −6 to 0, fraction Csp³ 0.25 or more and at most 9 rotatable bonds. A molecule drawn entirely inside the band is predicted to be orally bioavailable.

SwissADME uses XLOGP3; this page uses Wildman–Crippen logP. Solubility is estimated with ESOL (Delaney, J. Chem. Inf. Comput. Sci. 2004, 44, 1000), logS = 0.16 − 0.63 logP − 0.0062 MW + 0.066 RB − 0.74 AP, again with Wildman–Crippen logP in place of ClogP, so treat it as a rough estimate.