FlexAppeal molecular dynamics, prepared and analysed

FlexAppeal

← All runs

beta2 adrenergic receptor, active start

Simulated
50ns
Frames
499
Final RMSD
4.61Å

mean of the last 10% of frames

Radius of gyration
24.12Å

trajectory mean

Hydrogen bonds
508

occupancy > 10%

Helix content
81%

trajectory mean

Throughput
60ns/day

OpenCL

Wall time
13.3h
Equilibrated
10ns

restrained, before production

Settles after
3.2ns

6% of production, by pymbar.detect_equilibration

Force field
amber14-all.xml
Platform
OpenCL
Structure
4LDE.pdb

Structure

The decimated trajectory, 499 frames of 5244 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.

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 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
H8 11, 329 to 339 0.68 1.21 +0.279 +0.199 to +0.214 active 8GE3 / 8W1V
NPxxY 5, 322 to 326 0.34 3.09 +0.804 +0.718 to +0.755 active 7DHI / 8W1V
TM5 35, 196 to 230 1.15 1.11 -0.019 +0.000 to -0.087 intermediate 7DHI / 3PDS
TM6 31, 269 to 299 1.54 3.43 +0.380 +0.432 to +0.353 active 3P0G / 9RKI
TM7 25, 304 to 328 0.70 2.15 +0.511 +0.411 to +0.427 active active 7DHI / 6PS4
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.8 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
R3.50 (R131) - L6.34 (L272) 12.41Å 12.41Å +0.00
R3.50 (R131) - E6.30 (E268) 13.07Å 6.52Å -6.55 broken
Y5.58–Y7.53 4.90Å 4.64Å -0.26
TM6 axis rotation 7.20° 23.10° +15.90
TM6 centroid shift 1.32Å 2.80Å +1.48
PIF connector RMSD 1.09Å 0.50Å -0.59
Helix 8 to the membrane plane 18.90° 14.00° -4.90
D2.50 to S3.39 2.73Å 2.66Å -0.07
D2.50 to N7.49 2.71Å 2.73Å +0.02
Nearest sodium to D2.50 26.57Å 26.56Å -0.01
W6.48 (W286) χ₁ -65.20° -61.50° +3.70
W6.48 (W286) χ₂ 113.40° 99.40° -14.00
Y7.53 (Y326) χ₁ -62.30° -73.60° -11.30
Y7.53 (Y326) χ₂ -76.50° -68.90° +7.60
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.

Bilayer quality, chunk by chunk

What this shows

The bilayer QC, chunk by chunk, against both of section 5.8's windows: the expected one it should sit in and the wider acceptable one outside which the membrane invalidates the run. A chunk between them is reported rather than failed.

Superposition core

Core residues 82
Core defined by generic 38 to 58 of TM1, TM2, TM3, TM4
Alignment target 8UNW
References aligned 100
References dropped 0
Mean core RMSD 0.52 Å
Worst core RMSD 1.08 Å
Rigidity reference 2RH1
Mean rigidity RMSD 0.84 Å
Worst rigidity RMSD 1.368 Å (3SN6)
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
2R4R A 3.40 inactive inactive 8.65 agrees not used
2R4S A 3.40 inactive inactive 8.65 agrees not used
2RH1 A 2.40 inactive inactive 8.40 agrees not used
3D4S A 2.80 inactive inactive 8.15 agrees not used
3KJ6 A 3.40 inactive inactive 8.53 agrees not used
3NY8 A 2.84 inactive inactive 8.30 agrees not used
3NY9 A 2.84 inactive inactive 8.29 agrees not used
3NYA A 3.16 inactive inactive 8.23 agrees not used
3P0G A 3.50 active active 14.13 agrees not used
3PDS A 3.50 inactive inactive 8.27 agrees not used
3SN6 R 3.20 active active 15.41 agrees not used
4LDE A 2.79 active no opinion not used
4LDL A 3.10 active no opinion not used
4LDO A 3.20 active no opinion not used
4QKX A 3.30 active no opinion not used
5D5A A 2.48 inactive inactive 8.35 agrees not used
5D6L A 3.20 inactive inactive 8.43 agrees not used
5JQH A 3.20 inactive no opinion not used
5X7D A 2.70 inactive inactive 8.34 agrees not used
6MXT A 2.96 active no opinion not used
6N48 A 3.20 active no opinion not used
6OBA A 3.10 inactive inactive 8.06 agrees not used
6PRZ A 2.80 inactive inactive 8.04 agrees not used
6PS0 A 3.40 inactive inactive 8.18 agrees not used
6PS1 A 3.20 inactive inactive 8.22 agrees not used
6PS2 A 2.40 inactive inactive 8.01 agrees not used
6PS3 A 2.50 inactive inactive 8.22 agrees not used
6PS4 A 2.60 inactive inactive 8.12 agrees not used
6PS5 A 2.90 inactive inactive 8.27 agrees not used
6PS6 A 2.70 inactive inactive 8.14 agrees not used
7DHI R 3.26 active active 15.52 agrees not used
7XK9 A 3.40 active no opinion not used
7XKA A 3.10 active no opinion not used
8GDZ R 3.50 active active 15.72 agrees not used
8GE1 R 3.40 active active 15.72 agrees not used
8GE2 R 3.30 active active 15.74 agrees not used
8GE3 R 3.30 active active 15.77 agrees not used
8GE4 R 3.30 active active 15.65 agrees not used
8GE5 R 3.20 active active 15.55 agrees not used
8GE6 R 3.40 active active 15.69 agrees not used
8GE7 R 3.40 active active 15.70 agrees not used
8GE8 R 3.40 active active 15.51 agrees not used
8GE9 R 3.50 active active 15.62 agrees not used
8GEA R 3.50 active active 15.36 agrees not used
8GEB R 3.60 active active 15.36 agrees not used
8GEC R 3.50 active active 15.50 agrees not used
8GED R 3.50 active active 15.64 agrees not used
8GFV R 3.10 active active 15.85 agrees not used
8GFW R 3.00 active active 15.83 agrees not used
8GFX R 3.00 active active 15.89 agrees not used
8GFY R 3.00 active active 15.74 agrees not used
8GFZ R 3.00 active active 15.74 agrees not used
8GG0 R 2.90 active active 15.71 agrees not used
8GG1 R 3.00 active active 15.74 agrees not used
8GG2 R 3.00 active active 15.84 agrees not used
8GG3 R 3.10 active active 15.79 agrees not used
8GG4 R 3.20 active active 15.84 agrees not used
8GG5 R 3.20 active active 15.89 agrees not used
8GG6 R 3.30 active active 15.89 agrees not used
8GG7 R 3.30 active active 15.92 agrees not used
8GG8 R 3.30 active active 16.30 agrees not used
8GG9 R 3.40 active active 15.67 agrees not used
8GGA R 3.50 active active 16.04 agrees not used
8GGB R 3.50 active active 15.92 agrees not used
8GGC R 3.40 active active 16.09 agrees not used
8GGE R 3.50 active active 16.13 agrees not used
8GGF R 3.60 active active 16.52 agrees not used
8GGI R 3.50 active active 16.31 agrees not used
8GGJ R 3.50 active active 16.13 agrees not used
8GGK R 3.30 active active 16.02 agrees not used
8GGL R 3.20 active active 16.01 agrees not used
8GGM R 3.20 active active 15.95 agrees not used
8GGN R 3.30 active active 15.78 agrees not used
8GGO R 3.20 active active 15.82 agrees not used
8GGP R 3.20 active active 16.13 agrees not used
8GGQ R 3.40 active active 16.01 agrees not used
8GGR R 3.40 active active 16.06 agrees not used
8GGS R 3.50 active active 15.97 agrees not used
8GGT R 3.50 active active 16.15 agrees not used
8GGU R 3.40 active active 15.90 agrees not used
8GGV R 3.50 active active 16.16 agrees not used
8GGW R 3.60 active active 15.83 agrees not used
8GGX R 3.60 active active 15.66 agrees not used
8JJ8 F 3.20 active no opinion not used
8JJL A 3.20 inactive no opinion not used
8JJO F 3.40 inactive no opinion not used
8UNL R 3.30 active active 16.08 agrees not used
8UNM R 3.40 active active 15.80 agrees not used
8UNN R 3.40 active active 15.94 agrees not used
8UNO R 3.30 active active 15.84 agrees not used
8UNP R 3.30 active active 15.87 agrees not used
8UNQ R 3.30 active active 15.89 agrees not used
8UNR R 3.40 active active 16.11 agrees not used
8UNS R 3.40 active active 15.83 agrees not used
8UNT R 3.40 active active 15.96 agrees not used
8UNU R 3.50 active active 16.19 agrees not used
8UNV R 3.30 active active 15.85 agrees not used
8UNW R 3.40 active active 15.62 agrees not used
8UNX R 3.60 active active 16.37 agrees not used
8UNY R 3.60 active active 15.97 agrees not used
8UO0 R 3.50 active active 16.26 agrees not used
8UO1 R 3.60 active active 16.24 agrees not used
8UO2 R 3.60 active active 15.89 agrees not used
8UO3 R 3.50 active active 16.74 agrees not used
8W1V A 3.00 inactive inactive 8.12 agrees not used
9BUY R 2.90 active active 16.28 agrees not used
9CHU A 3.49 inactive no opinion not used
9CHX A 3.50 inactive no opinion not used
9L8L R 3.22 active active 15.97 agrees not used
9LW5 R 2.87 active active 15.38 agrees not used
9RKF A 2.45 inactive no opinion not used
9RKG A 2.45 inactive inactive 8.39 agrees not used
9RKH A 2.50 inactive no opinion not used
9RKI A 2.60 inactive inactive 8.48 agrees not used
9U9V R 2.85 active active 15.82 agrees not used
9W3F A 3.04 inactive inactive 8.36 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.

116 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 100 of 100 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.