Skip to content
Plate XIMethod

Science

Every quantity this archive shows is either something a laboratory measured or something derived here by arithmetic. This page defines both kinds: the notation, the mathematics, and the limits of what each number supports. Nothing on the site uses a formula that is not written out below.

Part One

The compound

What a molecule is, as this archive records it: how tightly it binds, and what it is built from.

Section 1

Affinity, and the pX scale

A binding constant is a concentration: the amount of compound needed to occupy a receptor. Lower means stronger. Because these concentrations span ten orders of magnitude — from picomolar to millimolar — they are compared on a logarithmic scale, which also turns the relationship the right way round, so that a larger number means tighter binding.

pX = −log₁₀(CM) = 9 − log₁₀(CnM)
CM
the measured concentration in molar
CnM
the same concentration in nanomolar, the unit the archive stores
pX
the negative base-10 logarithm of that concentration expressed in molar
9
the conversion from nanomolar to molar: 1 nM = 10−9 M
So a Ki of 1 nM is pX 9; 10 nM is pX 8; 1 µM is pX 6. Each whole unit is a tenfold difference in concentration. Every pX on this site is computed with this expression from the normalized measurement, never read from a stored column — the database holds a partially populated pX field, and an axis that is sometimes one quantity and sometimes another is not an axis.
Fig. 1One ruler, two scales
121 pM111091 nM8761 µM5431 mMpXconc.1 nM = 10⁻⁹ M = pX 9tighter binding
Computed · every point evaluated from the expression above. Binding constants span ten orders of magnitude, so they are read as logarithms. Each tick is a tenfold change in concentration and exactly one pX unit. The scales run in opposite directions because the logarithm is negated: the smaller the concentration needed, the tighter the binding, and the larger the pX.

The measurements, and why they are not interchangeable

Ki — inhibition constant
The concentration at which a compound occupies half the receptors, corrected for the radioligand competing with it. Derived from IC₅₀ through the Cheng–Prusoff relation below. An equilibrium property of the compound and the receptor, largely independent of assay conditions — which is what makes two Ki values comparable across laboratories.
Kd — dissociation constant
The same equilibrium quantity measured by direct saturation binding rather than competition. Kd = koff / kon: the ratio of the rates at which the complex falls apart and forms.
IC₅₀ — half-maximal inhibitory concentration
The concentration producing half the maximum inhibition in that assay. It depends on how much radioligand was used and on that ligand’s own affinity, so it is not portable between experiments.
EC₅₀ — half-maximal effective concentration
A functional readout: the concentration producing half the maximum response, not half the maximum occupancy. Receptor reserve means a compound can produce a full response while occupying a small fraction of receptors, so EC₅₀ can sit far below Ki at the same target.
AC₅₀ — half-maximal activity concentration
The concentration producing half the maximum change in an assay’s readout, whichever direction that change runs — activation or inhibition. It is the quantity reported by high-throughput concentration–response screens, where a curve is fitted to every compound tested and the assay does not say in advance which way a compound will move it. Like IC₅₀ and EC₅₀ it belongs to its assay, not to the compound and the receptor, and it is never read as a binding constant.
Potency
ChEMBL’s term for concentrations from PubChem’s screening assays, where the value is the concentration at which the compound was active. It is a screening result, not a binding experiment designed for that target.
Ki = IC₅₀ ⁄ (1 + [L] ⁄ Kd,L)
IC₅₀
the half-maximal inhibitory concentration measured in the competition assay
[L]
the concentration of radioligand used
Kd,L
the dissociation constant of that radioligand at the same receptor
The Cheng–Prusoff equation, which is why an IC₅₀ and a Ki are different numbers for the same interaction. This archive does not apply it: the radioligand identity and concentration are not recorded for most rows, so the correction cannot be made without inventing its inputs. Values are shown as the type they were published as, and the two are never pooled into one chart, ranking or average.
Fig. 2The same compound, two apparent potencies
050100IC₅₀ at low [L*]IC₅₀ at high [L*]an artifact of the assay, not a property of the compoundlabeled ligand still bound (%)competitor concentration →
Diagram · drawn to explain. No marks encode data. A competition assay reports the concentration that displaces half the labeled ligand. That concentration depends on how much labeled ligand it had to compete against, so the raw IC₅₀ belongs to the experiment, not to the compound. Cheng–Prusoff divides the assay conditions back out, and both curves yield one Ki. This is why the archive stores Ki and IC₅₀ as separate quantities and never averages across them.

Section 2

Structure: descriptors and scaffolds

Physicochemical descriptors are computed here with RDKit from the structure PubChem returns, never from a stored local copy. They are marked R throughout, because they are calculated rather than measured.

Molecular weight
Sum of atomic weights, in g/mol. Computed from the molecular formula.
cLogP
Calculated partition coefficient between octanol and water, by the Crippen method: a sum of contributions from each atom in its bonding context. Positive means lipophilic. It is an estimate from structure, not a measurement — an experimental logP is a different number.
TPSA
Topological polar surface area, in Ų. The summed surface contribution of nitrogen and oxygen atoms and their attached hydrogens, from a fragment table. Used as a proxy for membrane permeability; values under roughly 90 Ų are associated with CNS penetration, though that is a correlation, not a rule.
H-bond donors / acceptors
Counts of N–H and O–H groups, and of nitrogen and oxygen atoms respectively, by Lipinski's definitions.
Rotatable bonds
Single, acyclic, non-terminal bonds between heavy atoms — a measure of conformational flexibility.
Fraction sp³
The proportion of carbon atoms that are sp³-hybridised: the count of sp³ carbons over all carbons. A measure of three-dimensionality; flat aromatic molecules score near zero.

Bemis–Murcko scaffolds

A compound’s scaffold is what remains after every side chain is stripped away: the ring systems, and the linkers joining them. Formally, delete all atoms that are not in a ring and not on a path between two rings, and repeat until nothing more can be removed. Compounds reducing to the same scaffold share a framework.

The Families page groups by this rather than by curated class names, because a label like “tryptamine” is an assertion somebody has to stand behind, while a scaffold is a derivation anyone can repeat from the structure. The scaffold SMILES is printed with each family so the grouping can be checked rather than trusted.

Bemis GW, Murcko MA. The properties of known drugs. 1. Molecular frameworks. J Med Chem. 1996;39(15):2887–93.

Fig. 3What a scaffold keeps
ONreducestructure as recordedscaffold: rings + linkers
Diagram · drawn to explain. No marks encode data. The Bemis–Murcko reduction strips terminal substituents and keeps the ring systems with the bonds that join them. Compounds sharing a scaffold share a skeleton, which is a fact about the structures rather than an opinion about them — which is why the families page groups by scaffold instead of by curated class names. The skeleton drawn here is generic; it is not any compound in the archive.
Part Two

The receptor

What it binds to. Occupancy is not effect, transcript is not protein, and neither gap is closed by arithmetic.

Section 3

Receptors, occupancy and efficacy

Binding and effect are separate questions. A compound’s affinity says how readily it occupies a receptor; its efficacy says what the receptor does once occupied. A high-affinity antagonist and a high-affinity agonist can have identical Ki values and opposite consequences, which is why the mechanism records on each dossier carry an action type and the binding data does not stand in for one.

occupancy = [C] ⁄ ([C] + Kd)
[C]
free concentration of the compound at the receptor
Kd
the dissociation constant at that receptor
The Hill–Langmuir relation at equilibrium. It follows that at [C] = Kd, half the receptors are occupied — which is what makes a dissociation constant interpretable as a concentration. This archive plots affinities, not occupancies: computing occupancy needs a free brain concentration, which is not measured for these compounds and would have to be invented.
Fig. 4Occupancy is gradual
01050901000.01×0.1×Kd​10×100×10 % → 90 % occupancy spans 81× in concentrationreceptors occupied (%)free concentration, relative to Kd →
Computed · every point evaluated from the expression above. Every point is the expression above, evaluated. Half the receptors are occupied when the free concentration equals Kd — which is what makes a dissociation constant readable as a concentration. Going from a tenth of the receptors occupied to nine tenths takes an eighty-one-fold change in concentration, so there is no dose at which a receptor population switches on. The archive plots affinities and not occupancies, because occupancy needs a free brain concentration that nobody has measured for these compounds.
Fig. 5Identical affinity, opposite consequence
same Kd, same occupancy, same concentrationagonistreceptor signalsantagonistreceptor occupied, silentoutin
Diagram · drawn to explain. No marks encode data. Both ligands occupy the same fraction of receptors at the same concentration — the curve in Fig. 4 is the same curve for both. What differs is what the receptor does once occupied, and no binding constant records that. A Ki table therefore cannot tell you whether a compound activates or blocks a target, and this archive never infers a direction of effect from an affinity.

Target classification

Assay target names arrive as free text from ChEMBL and mix three different kinds of thing: named molecular targets (5-HT2A), whole organisms (Plasmodium falciparum), and cell lines (HepG2). They are sorted by an explicit pattern list, and each classified target carries the pattern that classified it, visible in the interface. A target matching no rule is shown as unclassified rather than assigned a class by guesswork.

Only molecular targets carry an affinity interpretation. A cytotoxicity value against a cell line is a real measurement of something else, and is never placed on a pX axis or drawn in the affinity rosette.

Section 4

Expression: nTPM and what it is not

The anatomy plates shade a region by how strongly the gene for a compound’s target is transcribed there, from the Human Protein Atlas consensus dataset.

nTPM — normalized transcripts per million
TPM
transcripts per million: read counts per transcript, corrected for transcript length, scaled so each sample sums to one million
n
trimmed-mean normalization across samples, so values are comparable between tissues
Shading runs on √nTPM rather than nTPM. Expression spans three orders of magnitude, and on a linear ramp every region except the strongest would render as blank — reading as “not expressed” rather than “expressed less”. The square root compresses that range without reordering anything.

Transcript is not protein, and protein is not activation. mRNA abundance correlates imperfectly with receptor density, and a receptor present in a region is a receptor available to be bound — not evidence that a compound reaches it, occupies it, or changes what that region does. The plates therefore claim only that: this compound binds this target, and this target’s gene is transcribed here. Shading is never affinity, and never affinity multiplied by expression, because that product is not a quantity anyone measures.

Part Three

Time

Binding, acting, and the hours that follow — three processes on clocks that differ by six orders of magnitude.

Section 5

Kinetics of drug action

Binding, acting, and the behavior that follows — three processes on three timescales, linked but not interchangeable.

Every quantity defined so far is an equilibrium: a state the system settles into given enough time. Nothing anyone experiences is at equilibrium. A dose arrives, receptors fill, something downstream happens, and hours later it stops — and the three stages run on clocks that differ by six orders of magnitude. This section defines each stage, and says which of them the archive can put a number on.

1 · Binding — what the affinity column leaves out

A compound and a receptor associate at a rate proportional to how much of each is present, and the complex falls apart at a rate that depends on nothing but itself. Affinity is the ratio of those two rates, and a ratio is a quotient of two numbers neither of which it reports.

R + L ⇌ RL,   Kd = koff ⁄ kon
kon
association rate constant, in M−1s−1. Second order: the rate of binding rises with concentration
koff
dissociation rate constant, in s−1. First order: an occupied receptor releases its ligand at a rate independent of what is around it
Kd
the equilibrium the two rates balance at — the concentration at which half the receptors are occupied
Fig. 6Two rates, one constant
outinR + Lkon · [R][L]second order · rises with concentrationkoff · [RL]first order · independent of concentrationRL
Diagram · drawn to explain. No marks encode data. Association depends on how much compound is present; dissociation does not — once a complex exists, it comes apart at a rate set by the complex alone. Affinity is the ratio of the two, which means it is silent about either. Two compounds can share a Kd and differ a hundredfold in how long a single binding event lasts.

The quantity that governs how long an occupied receptor stays occupied is the off-rate alone. It has its own half-life, and it is not the drug’s.

τ = 1 ⁄ koff,   t½,occupancy = ln(2) ⁄ koff
τ
residence time: the mean lifetime of a single receptor–ligand complex
t½,occupancy
the time for half the occupied receptors to release their ligand, once free compound is gone
This archive holds no rate constants. ChEMBL and the PDSP database report equilibrium quantities — Ki, Kd, IC₅₀ — and no arithmetic recovers kon and kofffrom their ratio. Residence time is defined here because it explains what follows, not because any page displays one. Where you see a number on this site, it is an equilibrium constant.
Fig. 7Same affinity, different persistence
050100010203040fast off-rateslow off-ratehalf the receptors still occupiedreceptors still occupied (%)minutes after the free compound is removed →
Computed · every point evaluated from the expression above. Occupancy after the free compound is removed, evaluated as e−kofft for two dissociation rates. Give the faster one a proportionally faster association rate and the two have identical Kd values — they are indistinguishable in every binding table on this site, and they hold the receptor for very different lengths of time. Residence time is not recorded in the archive: the measurements it holds are equilibrium constants, and no amount of arithmetic recovers a rate from a ratio.

2 · Acting — occupancy is not response

An occupied receptor changes shape, couples a transducer, and sets a second messenger moving. Each step takes its own time — milliseconds for the conformational change, seconds to minutes for the messenger — and, more consequentially, each step multiplies. One receptor activates many transducers; one enzyme makes many molecules of product. By the time a response is measured, the signal has been amplified enough that a small occupied fraction can produce a nearly maximal effect.

E ⁄ Emax = τA·occupancy ⁄ (τA·occupancy + 1)
τA
the operational efficacy: how much of the system's amplification a given agonist can recruit. Large for a full agonist in an amplifying tissue, near zero for an antagonist
occupancy
the fraction bound, from the Hill–Langmuir expression in section 3
E
the measured response
The Black–Leff operational model. It makes the separation explicit: affinity sets occupancy, τA sets what the tissue does with it. Two consequences follow for how this site handles data. An EC₅₀ measured in an amplifying system sits far below the Kd at the same receptor, so functional and equilibrium values are never pooled. And τA is a property of the agonist and the assay — the same compound reports different efficacies in different systems, which is why no efficacy is displayed as though it were a property of the compound.
Fig. 8Most of the response, few of the receptors
050100Kd · half-maximal occupancyEC₅₀ · half-maximal responsegain contributed by the signaling cascadeper cent of maximumconcentration →
Computed · every point evaluated from the expression above. Occupancy and response against the same concentration axis, evaluated. Each step between binding and readout multiplies, so the response curve sits to the left of the binding curve: half-maximal effect arrives well below half-maximal occupancy. An EC₅₀ is therefore a property of the assay’s amplification as much as of the compound, which is why this archive keeps functional readouts in a separate table from equilibrium constants and never averages the two into one number.

A further complication the archive records but does not resolve: a receptor can couple more than one transducer, and an agonist can favor one over another. A compound biased toward G-protein coupling and one biased toward β-arrestin recruitment can share a Ki and drive different downstream courses. The mechanism records on each dossier carry an action type where a source states one; where none does, the dossier says so rather than inferring a direction from a binding constant.

3 · Behavioral consequences — why the clocks diverge

The compound is in plasma; the receptors are not. Between the two sits distribution across the blood–brain barrier, and the compartment where the effect is generated fills and empties more slowly than the blood does. The standard treatment adds a hypothetical compartment whose concentration lags plasma by a single rate constant.

dCe ⁄ dt = ke0 · (Cp − Ce)
Cp
plasma concentration — what a blood sample measures
Ce
concentration in the effect compartment — what the receptors experience
ke0
the rate at which the effect compartment equilibrates with plasma. Small values mean a long lag
The Sheiner link model. Plot effect against plasma concentration through a whole exposure and the result is a loop rather than a curve, because the same blood level corresponds to two different effects depending on whether it is rising or falling.
Fig. 9Effect lags concentration
counterclockwiseeffect lags concentrationeffectplasma concentration →plasma concentrationeffecttime →
Computed · every point evaluated from the expression above. Left: effect against plasma concentration, traced through one exposure. It is a loop, not a curve, and it runs counterclockwise — at any given concentration the effect is smaller on the way up than on the way down. Right: the same two quantities against time, where the lag is visible directly. Peak effect arrives after peak concentration and outlasts it, so a blood level neither predicts nor times an experience. The constants are chosen to show the shape; they belong to no compound.

Three mechanisms stack in the same direction, and together they are why subjective effects outlast the plasma curve so consistently: the effect compartment empties after plasma does; receptors stay occupied after free compound is cleared, at a rate set by koffrather than by elimination; and amplification means a falling occupancy still drives a near-maximal response until it drops a long way. In some compounds a further term applies — an active metabolite with its own kinetics, which is why psilocybin and psilocin are filed as two records rather than one.

None of this is used to produce a number. The mechanisms above explain a direction, not a magnitude: no ke0, no koff and no free brain concentration is recorded for these compounds, so a duration cannot be calculated here even in principle. Every duration on this site is quoted from a trial that measured it in people, and where no trial has, the dossier shows an absence.

Section 6

Half-life and duration

These are different quantities, they diverge sharply, and conflating them is the single easiest way to misread a dossier.

t½ = ln(2) ⁄ ke ≈ 0.693 ⁄ ke
t½
the time for plasma concentration to fall by half
ke
the first-order elimination rate constant
Elimination half-life describes how fast the body clears a drug. Where a compound is cleared in two phases, an α (distribution) and a β (terminal) half-life are reported separately, as for ergonovine.

Duration of subjective effects is not derived from half-life. LSD is cleared with a plasma half-life of about 2.6 h while its subjective effects run 8–12 h — a factor of three or four. Effects persist through receptor occupancy and downstream signaling that outlast the plasma curve. Every duration on this site is quoted from a controlled human trial that measured it, never inferred from pharmacokinetics and never borrowed from a chemical relative. Why the two diverge so reliably is set out in section 5.

Fig. 10Two clocks, and they do not agree
Measured · every value cited on the dossier it links to. Elimination half-life beside duration of acute subjective effects, on one logarithmic axis. Of the compounds shown, 3 carry a cited value for both, and in each the effects outrun the plasma curve severalfold. That gap is the reason no duration on this site is ever calculated from a half-life. Bars span the range a trial reported; a single reported value is drawn as a mark, not widened into a range nobody measured. Follow a compound to the dossier that quotes its source.
Part Four

Analysis

Comparing incomplete profiles, and what each way of handling a gap costs.

Section 7

Similarity, distance and clustering

Two plates, two questions, and one decision — what to do about a gap — that separates them.

The analysis page carries two matrices. The first asks which targets a compound binds — one compound per row, one target per column. The second asks how alike two compounds’ whole profiles are — compounds on both axes. Both are built from the same evidence, either the archive’s measured affinities or DRIFT’s structural scores, and the two sources are never combined into one grid.

Both matrices are mostly empty: 95% of the measured grid and 92% of DRIFT’s. Any comparison between two compounds therefore has to decide what an empty cell means, and there are only two coherent answers. Leave it out, and the comparison rests on whatever the pair happens to share. Fill it with a zero, and the comparison spans everything — at the cost of asserting a value nobody measured. This site does both, in different plates, and says which is which.

Fig. 11One profile pair, two ways to handle a gap
pairwise-complete, in pX · drug–target platecompound A7.9—8.46.2——9.1—compound B8.26.6—5.9—7.48.8—3 of 8 targets contribute · the other 5 are set aside, not filledcomplete vector, as 1/Kᵢ · compound–compound platecompound A0000compound B0008 of 8 contribute · absence is 0 because that is the limit of 1/Kᵢ
Diagram · drawn to explain. No marks encode data. Above, in pX: only the targets both compounds were measured at enter the distance, and the rest are set aside. Nothing is invented, and in exchange the comparison rests on three columns. Below, the same values as 1/Ki: every target enters, and an empty one enters as zero — not as a substitution but because zero is where 1/Ki goes as Ki goes to infinity. Notice what else the reciprocal does. The two weakest measurements, at pX 6.2 and 5.9, are drawn to scale and have all but vanished, so a target both compounds merely happened to be tested at contributes almost nothing on its own.

1 · Leaving the gap out

The drug–target plate orders its axes this way. Only cells where both members carry a value contribute, and nothing is imputed.

d(a,b) = √( Σi ∈ S (va,i − vb,i)² ⁄ |S| )
S
the set of targets at which BOTH members have a value
va,i
compound a's value at target i — a pX, or a DRIFT score
|S|
the number of shared targets — the overlap
Root-mean-square over shared observations only, so a pair is not judged more distant merely for having been tested more often. Treating “not measured” as “does not bind” would cluster compounds by how thoroughly they have been studied. Pairs sharing fewer than three values are treated as incomparable rather than similar, and are placed at the far end of the observed range — so they end up adjacent without that adjacency meaning anything.

The cost is coverage, and the plate reports it rather than hiding it. On DRIFT’s full grid only 12% of target pairs share three values, so most of that axis is ordered by tie-breaking; the page says so in place, next to the ordering it qualifies. The compound axis fares much better at 82%.

2 · Giving the gap a value, in the coordinate where zero means something

The compound–compound plate takes the other option, because a cosine is defined over a whole vector rather than over an intersection. Each compound becomes a vector across every target of its source — all 172 assayed, or all 489 DRIFT ranked — and every cell has to hold a number, including the empty ones. Which number is admissible depends entirely on the coordinate.

cos θ = (a · b) ⁄ (‖a‖ ‖b‖)
a, b
two compounds’ vectors over the complete target space, carrying 1/Ki at each target and zero where absent
a · b
the dot product, Σi aibi over every target
‖a‖
the vector's length, √(a · a)
1 is the same direction, 0 is orthogonal. The diagonal is 1 by construction, not by finding. A compound whose profile is entirely empty has no direction at all, so its row is left blank rather than scored as similar to nothing. Cosine normalizes each vector by its own length, so the units of the component cancel and 1/nM and 1/M give the same answer.

The component is 1/Ki, and that is what makes the zero admissible. A reciprocal binding constant tends to zero as Ki tends to infinity, so an empty cell holds the limit of the quantity itself at “does not bind”. Nothing is substituted for a missing value; the coordinate simply has a value there.

Measuring in pX instead — which is what this plate did first — puts a real and absurd number in the same place, because pX 0 is a binding constant of one molar. The consequence was not only cosmetic. Every component was then of the same order, so a target a pair had merely both been tested at contributed as much as one either of them binds, and the cosine largely tracked co-assay rather than pharmacology. In reciprocal space a target bound at 10 µM contributes about a hundred-thousandth of what one bound at 100 pM does. The effect is easy to see on real pairs: LSD against lisuride, two ergolines, scores 0.43 where an unrelated pair like LSD against atropine scores 0.00, against 0.75 and 0.08 under pX.

DRIFT’s score needs no reciprocal. It is already bounded at zero, and a target absent from a compound’s ranking is one the model found no structural resemblance to, so zero is already what absence means.

What the reciprocal does not fix. An empty cell still means “nobody ran this assay”, which is a different fact from “this compound does not bind here”. Zero is the right value for the second and says nothing about the first, and no method that requires a complete vector can tell them apart. That is the standing reason the drug–target plate above uses the other treatment, and the reason both are on the page.

The two plates can therefore disagree about which compounds belong together, and that disagreement is informative rather than a fault. Where they agree, the grouping survives both treatments of the missing data. Where they do not, the difference is usually telling you about the assay record rather than about the molecules.

3 · Building the tree

Both plates order their axes by average-linkage agglomeration: begin with every item alone, repeatedly join the two closest clusters, and read the leaves off the finished tree. Distances between merged clusters are updated by the Lance–Williams recurrence rather than recomputed from their members.

d(i∪j, k) = ( |i|·d(i,k) + |j|·d(j,k) ) ⁄ ( |i| + |j| )
i, j
the two clusters being joined
k
any other live cluster
|i|
the number of leaves in cluster i — the weighting that makes this an average over all member pairs
For the compound–compound plate the distance is 1 − cos θ. Because zero-filling made every vector complete, every pair there is comparable and no fallback distance enters the tree at all — the exact trade the drug–target plate refuses.

4 · Reading the tree, and the shaded groups

A join is drawn at the distance its two sides merged at, so depth reads as disagreement. Joins at the fallback distance — pairs too sparse to compare — are drawn as hairline dashes, because they are the algorithm running out of alternatives rather than a finding.

Fig. 12Reading the tree, and where the bands come from
tightness0.00.51.0distance at joinABCDEFgroupgroupnot a group
Diagram · drawn to explain. No marks encode data. A join sits at the distance its two sides merged at, so depth reads as disagreement: A with B and C with D are tight, E with F much less so. The plate shades a group when its members joined below the tightness line — the closest quartile of the distances actually observed on that axis — which here keeps two groups and leaves E and F unbanded. Two other cuts were tried first and both misreport this data; the text below says why.

A group is shaded when its members joined below the tightness line: the closest quartile of the distances actually observed on that axis. Runs of one are never banded, so a singleton is not dressed up as a cluster. The threshold is derived from each matrix rather than fixed, so it means the same thing on a pX axis and on a cosine axis.

Two more obvious cuts were tried first and both misreport this data, which is worth recording because both are standard. Cutting at the fallback distance returns the entire axis as one group: average linkage blends real distances with the fallback as clusters grow, so even a tree whose pairs were mostly incomparable joins below it. Cutting into a fixed number of top-level groups fails the same way and for the same reason — with this much missing data the tree chains, absorbing one leaf at a time, so splitting the tallest join repeatedly peels off singletons and leaves one cluster holding 106 of 112 compounds. Cutting on tightness says the same thing whatever the shape of the rest of the tree.

5 · What DRIFT’s score is

FP2 similarity = |A ∩ B| ⁄ |A ∪ B|
A, B
the sets of FP2 fingerprint bits set for two structures
∩ , ∪
bits present in both, and bits present in either
The Tanimoto coefficient over path-based fingerprints, which is what DRIFT’s score reports. It measures structural resemblance to known ligands of a target — not predicted affinity. DRIFT’s own API states that its target-screen predictions do not affect this ranking, and this archive keeps predicted values in a separate grid from measurements for that reason. DRIFT also returns an estimated pK for some rows; none of them are shaded anywhere on this site, and where one appears in a cell readout it is marked as an estimate.

Just over half of what DRIFT returns is not a molecular target at all — cell lines, organisms, tissue preparations, and the placeholders “No relevant target”, “Unchecked” and “ADMET”. All of it is excluded by the same rule that keeps such rows off the measured grid, so the DRIFT axes carry 489 molecular targets out of the 19,715 rows returned.

Part Five

Clinical description

The vocabulary a clinical account uses, and how this archive measures what its literature says under each heading.

Section 8

Domains, and the evidence measure

A glossary of clinical terms, and a count of how much of the recorded literature speaks to each.

A dossier quotes trials and monographs that describe what a compound does to people, and those descriptions use clinical vocabulary in its clinical sense. Two of the terms are ordinary English words with narrower meanings here — affect is not a synonym for mood, and neither is a synonym for feeling — so the domains are defined once, below, rather than left to context.

These are the domains of a mental status examination: the headings a clinician describes a person under. None of them is a property of a molecule, and nothing in this glossary asserts that any compound produces any state.

Affect
The observable, moment-to-moment expression of emotion — described by its range and quality, as flat, blunted, restricted, labile or reactive.
/affective|(?:flat|blunted|restricted|constricted|labile|inappropriate|incongruent|pseudobulbar)\s+affect|affect\s+(?:lability|blunting)|emotional (?:blunting|lability|numbing)|tearful/i
Mood
The sustained, internally experienced emotional tone the person reports, as euthymic, dysphoric, euphoric, anxious or irritable.
/\bmood|euphor(?:ia|ic)|dysphor(?:ia|ic)|anxi(?:ety|ous|olytic|ogenic)|depress(?:ion|ed|ive)\b|irritab|panic|elation|apath/i
Perception
Sensory experience and its disturbances: illusions, in which a real stimulus is misinterpreted, and hallucinations, which occur without one.
/hallucinat|hallucinogen|illusion|perceptu|percept\b|visual distort|synaesthes|synesthes|depersonaliz|derealiz|sensory distort|altered (?:state|perception)/i
Cognition
Attention, concentration, orientation, learning and memory, and executive function — what a formal examination tests directly.
/cognitiv|memory|amnesi|attention|concentrat|confus|disorient|executive function|working memory/i
Thought process
The form and organization of thinking rather than its content: linear and goal-directed, or circumstantial, tangential or loosened.
/thought (?:process|disorder|form)|disorganiz|tangential|incoheren|flight of ideas|loosening of association/i
Thought content
What is being thought about — preoccupations, overvalued ideas, delusions, and thoughts of harm to self or others.
/delusion|paranoi|persecut|suicid|obsess|grandios/i
Psychomotor activity
The rate and quality of movement and speech, from retardation through to agitation.
/psychomotor|agitat|sedat|catatoni|restless|stupor|somnolen|hyperactiv|akathisi/i
Insight and judgment
Whether the person recognizes the experience as drug-induced or as externally real, and the soundness of decisions made while it lasts.
/\binsight|judg(?:e)?ment/i

The measure

Each compound page carries a web with one spoke per domain. The spoke is a count of records, not a rating of severity.

Ed = | { s ∈ S : match(td, text(s)) } |
Ed
the evidence count for domain d — the length of its spoke
S
every record this archive holds as prose for the compound: PubChem's pharmacology passages, its toxicology passages, and the curated trial quotations
td
the term list for that domain, printed above and matched case-insensitively
Distinct records, not word occurrences. A passage using “hallucination” six times is one record saying something about perception, and counting the sixes would measure prose style. Every count resolves to the records behind it, which are listed under the web.

It measures the literature, not the compound. A long spoke means the subject has been written about; it says nothing about how strong the effect is. A spoke at the center means nothing recorded here discusses that domain, not that nothing happens — Thought process is currently zero for every compound in the archive, which is a fact about what these sources chose to report. Comparing two compounds compares how much has been published about each, and comparing spokes is not comparing effects.

A severity profile would be a different object and would need per-domain magnitudes from a validated instrument — 5D-ASC, ARCI, HRS or MEQ30 subscale scores from human trials. This archive holds none, and the previous build’s attempt at one is listed under held back for having had no source column at all. Counting sources is what the evidence supports, so counting sources is what the web does.

Two deliberate exclusions. Paper titles from the bibliography are not scanned: a title states a study’s subject, not its findings, and one well-titled paper would inflate a domain it never reported on. And the term lists are conservative about clinical homonyms — in this corpus the everyday verb “affects” outnumbered the clinical noun eight to one, so the Affect pattern requires a qualifier and bare “affects” matches nothing.

Reference

Notation used throughout

  • PPubChem record — identity, structure, synonyms
  • CChEMBL record — a measured bioactivity
  • Rcomputed here with RDKit from an upstream structure
  • Dderived here by arithmetic on a stored measurement
  • —an explicit absence: no value exists. Never a zero, never a default