FlexAppeal molecular dynamics, prepared and analysed

FlexAppeal

← All runs

CB1 cannabinoid receptor with taranabant, in POPC (200 ps verification)

Simulated
0.200ns
Frames
20
Final RMSD
0.99Å

mean of the last 10% of frames

Radius of gyration
29.83Å

trajectory mean

Hydrogen bonds
680

occupancy > 10%

Helix content
64%

trajectory mean

Throughput
64ns/day

OpenCL

Wall time
0.1h
Equilibrated
0.02ns

restrained, before production

Settles after
0.0ns

11% of production, by pymbar.detect_equilibration

Force field
amber14-all.xml
Platform
OpenCL
Structure
5u09_opm.pdb

Structure

The decimated trajectory, 20 frames of 7775 atoms. Press play, or scrub the frame slider: the plots below highlight the frame you are on.

Two-dimensional structure of the ligand, each atom shaded by how far it moved over the scored window
mobility Each atom is shaded by how far it moved over the scored window, in six bands of half an angstrom: deep indigo for an atom that held still, through to amber for 3.0 angstrom and over. The bands are the same on every diagram, so a still ligand reads as still rather than being stretched to fill the scale.
0.0 ns

Loading the structure…

The binding site

A 96.5/100
Pose 100.0 33%
Occupancy 100.0 22%
Interactions 99.8 22%
Energy -60.2 kcal/mol ± 3.2 100.0 17%
Convergence not measured here — 10%
Conformation 37.4 6%

Ligand pose

What this shows

Two fits of the same ligand. Against the pocket, the curve is the ligand moving in a receptor held still; against itself, it is the ligand changing shape. The shaded head is the equilibration that was measured and not scored.

Centroid drift

What this shows

Displacement of the ligand centroid from where it started. The dotted line at 5 Å is the drift beyond which DeepFlex caps a grade.

Pocket occupancy

What this shows

Receptor residues in contact with the ligand. The dotted line is a quarter of the contacts present at the start, which is the threshold for counting a frame as still bound.

Interaction occupancy

What this shows

Which contacts carried the pose, and for how much of the scored window. Typed from heavy-atom geometry, so these are polar, ionic and hydrophobic contacts rather than a full interaction fingerprint.

Contact barcode

What this shows

When each contact was present, busiest first. Gaps are the contact breaking and reforming.

Pose clusters

What this shows

How the scored window divides into distinct poses. One tall bar near zero is a ligand that held its pose.

Ligand strain

What this shows

MMFF94 strain in excess of the starting pose, never the absolute energy. Below zero means the ligand relaxed out of the pose it started in, which is a good outcome rather than a neutral one.

Ligand torsions

What this shows

Each rotatable torsion over the run. Flat traces are a rigid ligand; a trace jumping between bands is a rotamer flip.

Ligand atom fluctuation

What this shows

Which part of the ligand moves, about its own mean position.

Receptor flexibility

What this shows

Per-residue backbone fluctuation, with the binding-site residues marked. A mobile pocket makes every pose measurement noisier.

Frame-to-frame RMSD

What this shows

Every frame against every other. Blocks along the diagonal are distinct states; a uniform field is a run that stayed in one.

Residue cross-correlation

What this shows

Which parts of the receptor move together. Blue is correlated, red is anticorrelated.

MM/GBSA, term by term

What this shows

Where the binding energy comes from. The bars are the difference block, complex minus receptor minus ligand, which is the only part of an MM/GBSA calculation that is a binding quantity. Error bars are the standard deviation over the sampled frames, not the standard error. igb=2, single trajectory, no entropy term: good for ranking ligands against each other, not an affinity.

The whole system

RMSD

What this shows

How far the structure has moved from its reference. A plateau means the trajectory has settled; a steady climb means it has not.

Radius of gyration

What this shows

Global compactness. A rise suggests unfolding or expansion.

Per-residue fluctuation

What this shows

Which parts of the molecule move. Peaks are usually loops and termini.

Solvent exposure

What this shows

Total solvent-accessible surface area over time.

Secondary structure content

What this shows

Helix and sheet content over the trajectory.

Secondary structure timeline

What this shows

Per-residue assignment at every frame. Vertical stripes are transient melting; horizontal bands are stable elements.

Contact map

What this shows

How often each residue pair is in contact across the whole trajectory.

Essential dynamics

What this shows

The two dominant collective motions, with each frame coloured by time. Distinct clouds mean distinct conformational states.

Conformational clusters

What this shows

Representative states found by RMSD clustering.

Ligand pose stability

What this shows

How far the ligand moves within the binding site, with the protein held fixed. A flat trace means the pose held; a jump means it moved or left.

Protein-ligand contacts

What this shows

Which residues the ligand touches, and for what fraction of the run.

Area per lipid

What this shows

The standard bilayer equilibration diagnostic: it should fall and then plateau. POPC sits near 0.63 nm² at 310 K. Two caveats on the absolute value: the protein's own cross-section is counted as membrane area, which inflates it for a large complex, and a freshly built bilayer starts loose and needs nanoseconds to condense.

Lipid chain order

What this shows

Order parameters down the acyl chain. A plateau near the headgroup falling away toward the tail is what a fluid bilayer looks like.

Bilayer thickness

What this shows

Phosphate-to-phosphate distance across the bilayer.

Convergence

What this shows

Each measure on its own axis. Energies, temperature and density differ by orders of magnitude, so sharing a plot between any two of them would make their alignment arbitrary and any apparent correlation fictional.

Reference states

The same trajectory measured against the experimental structures of this same protein in both of its known states, motif by motif, rather than against itself. Every threshold below is either published and cited or refitted and reported as refitted, and every reference structure was resolved from KLIFS or GPCRdb at run time.

State score over time

What this shows

The whole protein as one number, against the experimental active and inactive structures of this same protein: +1 is the active references, -1 the inactive ones. The thin line is per frame and the bold one is smoothed. The shaded bands start at s = ±0.2, which is where a label is decided; inside them the protein reads as intermediate.

Motif timelines

What this shows

Each structural motif scored separately, in the colours section 6.1 fixes for them. The whole-protein score above is these weighted by residue count, which is why a small motif moving a long way can leave it almost unchanged. Motifs of unequal discriminating power carry equal residue weight, so read this panel rather than the headline number.

Motifs at the end of the run

Motif Residues To active (Å) To inactive (Å) s First half to second Label Secondary structure leans Nearest active / inactive
H8 8, 401 to 408 1.86 0.85 -0.375 -0.418 to -0.380 inactive 8WRZ / 5U09
NPxxY 5, 393 to 397 2.03 0.33 -0.718 -0.608 to -0.639 inactive 8GHV / 9B9Z
TM5 30, 272 to 301 1.99 0.55 -0.568 -0.455 to -0.493 inactive 8IKG / 9B9Y
TM6 33, 337 to 369 4.54 0.93 -0.659 -0.638 to -0.636 inactive 5XRA / 5U09
TM7 29, 372 to 400 1.99 0.99 -0.336 -0.340 to -0.303 inactive inactive 6KPG / 5U09
What this shows

The last frame, motif by motif, with the distances the score is computed from and the nearest reference structure in each ensemble.

At 0.2 ns, the last frame analysed. A label is decided beyond s = ±0.2; inside that the motif reads as intermediate, which is a measurement and not a missing answer. The two half means say whether a motif was still moving: one number at the end cannot distinguish a motif that settled from one that was on its way somewhere. The last column is left blank where too few residues distinguish the reference sets to read a state from their secondary structure.

Secondary structure, motifs and the reference consensus

What this shows

Secondary structure per residue per frame, with each motif marked down the left edge and the reference ensembles' own consensus in the two columns on the right. A motif whose secondary structure matches the inactive references and not the active ones is the clearest evidence this analysis produces, and it is visible along a row.

Reference state map

What this shows

The reference structures on the plane their published criterion is drawn on, coloured by what the catalogue calls them, with this run's path across it in dark grey. This is the panel where a threshold that does not separate the clusters is obvious rather than argued about.

Superposition core

Core defined by generic 38 to 58 of TM1, TM2, TM3, TM4
Alignment target 23IV
References aligned 23
References dropped 0
Mean core RMSD 1.00 Å
Worst core RMSD 1.87 Å
Rigidity reference 23IV
Mean rigidity RMSD 1.05 Å
Worst rigidity RMSD 1.869 Å (5TGZ)
What this shows

What every distance on this page was measured after superposing on, and how well the reference structures agree about it.

The core is the part of the fold that does not move between states, so superposing on it leaves the motifs free to.

Reference structures

PDB Chain Resolution (Å) GPCRdb state Our state TM3-TM6 (Å) Agreement In the ensemble
23IV D 2.42 active 13.88 no opinion not used
23IW R 2.42 active 13.83 no opinion not used
5TGZ A 2.80 inactive 8.21 no opinion not used
5U09 A 2.60 inactive 6.30 no opinion not used
5XRA A 2.80 active 13.05 no opinion not used
6KPG R 3.00 active 13.80 no opinion not used
6KQI A 3.25 inactive 6.32 no opinion not used
6N4B R 3.00 active 14.14 no opinion not used
7FEE A 2.70 inactive 8.91 no opinion not used
7WV9 R 3.36 active 14.21 no opinion not used
8GHV D 2.80 active 14.30 no opinion not used
8IKG R 3.40 active 14.02 no opinion not used
8IKH R 3.30 active 13.86 no opinion not used
8K8J R 2.88 active 14.24 no opinion not used
8WRZ R 3.60 active 15.46 no opinion not used
8WU1 R 3.20 active 12.94 no opinion not used
9B65 R 3.03 active 13.99 no opinion not used
9B9Y R 3.50 inactive 6.35 no opinion not used
9B9Z R 3.30 inactive 6.56 no opinion not used
9BA0 R 3.13 inactive 6.29 no opinion not used
9ERX R 2.90 active 13.60 no opinion not used
What this shows

The ensemble this run was scored against, and whether our spatial label agrees with the annotation each structure came with. Counted, that agreement is the positive control for the whole method: if the labels did not reproduce the catalogue on structures whose state is known, nothing this page says about a trajectory would mean anything.

21 structures, resolved at run time from the catalogue rather than listed in the code. P3: the threshold on the TM3-TM6 distance, fitted on this set under leave-one-out, reproduces GPCRdb's own state call on 0 of 0 structures. No structure helped decide the line it is judged by. A structure reads "not used" when it has no usable state label or the alignment dropped it; the core panel says how many that was.