FlexAppeal molecular dynamics, prepared and analysed

FlexAppeal

← All runs

ABL1 kinase, active start (DFG-in)

Simulated
50ns
Frames
500
Final RMSD
1.71Å

mean of the last 10% of frames

Radius of gyration
19.34Å

trajectory mean

Hydrogen bonds
399

occupancy > 10%

Helix content
36%

trajectory mean

Throughput
224ns/day

OpenCL

Wall time
4.9h
Equilibrated
1ns

restrained, before production

Settles after
0.0ns

stationary from the first frame, by pymbar.detect_equilibration

Force field
amber14-all.xml
Platform
OpenCL
Structure
2GQG.pdb

Structure

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

0.0 ns

Loading the structure…

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.

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 over 2.0 ns. 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
A-loop 26, 381 to 409 6.23 9.73 +0.219 +0.204 to +0.205 active active 2HZ4 / 2HZ0
alphaC 11, 282 to 292 1.21 1.73 +0.177 +0.219 to +0.157 intermediate 2GQG / 3QRI
c.l 8, 361 to 368 0.87 1.12 +0.125 -0.067 to +0.018 intermediate 2F4J / 2HYY
g.l 6, 249 to 254 4.93 4.35 -0.062 +0.014 to -0.087 intermediate 4TWP / 2HIW
hinge 3, 316 to 318 0.61 0.55 -0.055 +0.019 to +0.039 intermediate 2GQG / 3QRJ
xDFG 4, 380 to 383 1.09 4.98 +0.642 +0.715 to +0.740 active 2F4J / 3QRJ
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 49.9 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.

What the published criteria measured

Measurement At the start At the end Change State at the end
Lys271 to Glu286 2.83Å 2.73Å -0.10 in
alphaC axis rotation 8.10° 4.90° -3.20
alphaC centroid shift 0.95Å 0.73Å -0.22
P-loop tip to hinge 14.86Å 19.98Å +5.12
R-spine assembled assembled
What this shows

The measurements every state call on this page is made from, at the first frame and the last. Each residue was resolved at run time from a KLIFS pocket position or a GPCRdb generic number, so the labels name the actual residues rather than a motif in the abstract. This is where a reader checks that the score agrees with the geometry.

Every residue here was resolved at run time from KLIFS pocket positions or GPCRdb generic numbers, which is why the labels name actual residues. A measurement with no reading at either end is left out rather than shown empty.

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 residues 43
Core defined by KLIFS pocket scaffold, excluding GK, a.l, alphaC, b.l, c.l, g.l, hinge, xDFG
Alignment target 2HYY
References aligned 19
References dropped 0
Mean core RMSD 0.42 Å
Worst core RMSD 0.84 Å
Rigidity reference 2E2B
Mean rigidity RMSD 0.46 Å
Worst rigidity RMSD 0.856 Å (2V7A)
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 (Å) KLIFS state Our state KLIFS DFG Our DFG Agreement In the ensemble
2E2B B 2.20 out out out out agrees not used
2F4J A 1.91 in in in in agrees not used
2G1T B 1.80 in in in in agrees not used
2GQG A 2.40 in in in in agrees not used
2HIW A 2.20 out out out out agrees not used
2HYY C 2.40 out out out out agrees not used
2HZ0 A 2.10 out out out out agrees not used
2HZ4 B 2.80 in in in in agrees not used
2V7A B 2.50 in in no opinion not used
3CS9 C 2.21 out out out out agrees not used
3PYY A 1.85 out out out out agrees not used
3QRI A 2.10 out out no opinion not used
3QRJ B 1.82 out out out out agrees not used
3QRK A 2.30 out out out out agrees not used
4TWP B 2.40 in in in in agrees not used
5HU9 A 1.53 out out no opinion not used
5MO4 A 2.17 out out out out agrees not used
6NPE A 2.15 out out out out agrees not used
6NPU A 2.33 out out out out agrees not used
6NPV A 1.86 out out out out agrees 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.

20 structures, resolved at run time from the catalogue rather than listed in the code. P1 and P2: our spatial DFG and alphaC labels agree with KLIFS on 17 of 17 structures that carry one. Both labels are measurements here and annotations there. 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.