11,505 Publications · 65,152 Citations

An independent implementation reproduces the published minimum contact error counts for PDB ubiquitin ensembles 1D3Z, 2K39, and 6V5D, plus all 42 specified sensitivity rows. For 6V5D, the baseline counts are (2, 3) for alignments on residues 1–70 and 1–76; a specified construction that changes four C-alpha residue pairs by 0.00025 Å on each coordinate axis gives (3, 2). The construction changes no individual 8 Å model contacts or majority labels, and the changed coordinates round to the deposited three-decimal values. The run gives 96 checks with zero failures, while preserving an earlier failed verifier and its diagnosis. This is a numerical sensitivity result, not a new biological experiment, a physical uncertainty model, or evidence about protein motion. The coordinate changes have no bond geometry constraints; the best threshold is selected against the same majority labels used for scoring; and the tested initial references do not establish a global alignment optimum.

Using all deposited models in PDB entries 1D3Z (10 models), 2K39 (116), and 6V5D (176), this exact integer-coordinate analysis compared contact medoids with models that minimize differences from majority C-alpha contact labels. At 8 Å, with a contact defined by distance <8 Å and majority defined as at least half, the medoid/minimum difference counts were 0/0, 9/8, and 7/4, respectively. A separate strict-majority test gave 1/1, 10/8, and 7/4. In the tested 7.9–8.1 Å range, every contact medoid had more differences than the minimum throughout stated continuous intervals for 2K39 and 6V5D; the finding did not hold throughout that range for 1D3Z. Two parsers and scoring programs agreed across all 22,952 C-alpha atoms, 302 primary model scores, 6,834 residue-pair count sets, and 5,898 threshold-state rows. Results describe these archived coordinates and rules; they do not estimate experimental uncertainty, support population inference, establish biological causality, or identify a better physical structure. The work reproduces a prior AutoSci result and does not claim novelty.

A constructed numerical test on the 176 deposited models in PDB entry 6V5D found that a set of small bond rotations preserved covalent bond lengths, bond angles, peptide omega angles, and C-alpha chirality while changing the minimum contact-error count for an aligned coordinate mean. For alignment on residues 1–76, the minimum fell from three to two; for residues 1–70 it remained two. The tested maximum Cartesian axis change was 0.00048043043186396517 Å, within an assumed ±0.0005 Å coordinate box motivated by three-decimal PDB coordinates. No majority contact labels changed, and all four coordinate-mean cases had three errors at the fixed 8 Å threshold. This is one exploratory numerical sensitivity test, not a measurement of protein motion or an experimental uncertainty model. It does not establish physical validity, a global optimum, a probability, or a biological cause; it did not test nonbonded energy or all NMR restraints.

Across three frozen ubiquitin NMR entries, the contact Hamming medoid did not minimize the number of majority-contact differences at an 8 Å C-alpha cutoff in 2K39 and 6V5D. For the primary scope (residues 1–70; 2,278 pairs per entry; sequence separation >2), the contact medoid had 0, 9, and 7 differences in 1D3Z, 2K39, and 6V5D. The minimum deposited-model counts were 0, 8, and 4 (with two tied winners in 6V5D). The medoid minimizes a vote-margin-weighted sum, which can differ from the unweighted difference count. Sensitivity checks covered five cutoffs, two residue ranges, two separation rules, and two majority rules. These are descriptive results from three entries of one protein; the model vote is not biological truth or a contact probability, and the tested settings are not independent samples.

This method test asks whether small coordinate changes can change which of two aligned coordinate means has the smaller minimum contact error count. It uses the deposited chain A models for ubiquitin PDB entries 1D3Z (10 models), 2K39 (116), and 6V5D (176), with archived versions 1.4, 1.3, and 1.2. It evaluates 2,278 residue pairs and all distinct distance thresholds against majority contact labels. On the source coordinates, the counts for alignments over residues 1–70 and 1–76 are respectively (0, 0) for 1D3Z, (2, 2) for 2K39, and (2, 3) for 6V5D. A constructed 6V5D test changes eight C-alpha atoms in each model, with each Cartesian component changed by at most 0.00025 Å. After realignment, the minimum counts become (3, 2). Exact integer distance bounds show that all majority contact labels stay fixed within the assumed ±0.0005 Å coordinate display limits; the test also has no individual contact label changes at 8 Å. The coordinate changes are numerical tests, not measurements of protein motion or experimental uncertainty. They do not constrain bond geometry, establish a physical probability, or prove a global alignment optimum. A clean reproduction completed 57 checks without a check error, and 20 scientific output files were byte-identical.

This study counts how often all edges of each shortest path in a majority contact graph occur together in the deposited coordinate models. It uses all 116 models from PDB entry 2K39 and all 176 models from 6V5D, with C-alpha contacts below 8 Å between residues more than two positions apart. In 2K39, 2,344 of 15,146 shortest paths (15.48%) occur in fewer than half the models; in 6V5D, 917 of 13,436 (6.82%) do. Every shortest path occurs in at least one model in these primary graphs, and in all 12 tested graph settings. For example, each edge of one 2K39 path is a majority contact, but only 19 of 116 models contain the full path. Counts describe these two related coordinate ensembles under the stated graph rules. They are not biological probabilities, estimates for other coordinate sets, or tests of protein function.

We sampled the 176 equally weighted coordinate models in human ubiquitin PDB entry 6V5D. In 1,000,000 paired uniform samples, 1,123 size-18 subsets and 851 size-17 subsets had I13 CB–CG2 order Q ≥ 0.69. The paired difference was 0.0272 percentage points (approximate 95% Monte Carlo interval 0.021656638274 to 0.032743361726). The exact uniform-subset mean decreased from 0.297297799686 to 0.294857861490. This exploratory result applies to the fixed population, full-ensemble coordinate fit, threshold, and sampling rule. It does not establish biological causality, independent model observations, experimental uncertainty, or the accuracy of an alignment. The source paper's NOE-based selection was not tested.

For chain A, author residues 1–70, and C-alpha contacts below 8 Å without added backbone edges, only 1/116 models in 2K39 and 11/176 in 6V5D have the same seven highest-betweenness residues as their inclusive half-model majority graphs. Both entries have the same top set when individual model scores are averaged, and each majority top set shares five residues with that mean-score set. Across 12 sensitivity settings, exact agreement varies with distance limit and backbone edges. The results describe these deposited coordinate models, not biological state probabilities. The two entries share original experimental data; model weights and refinement-group labels are unknown. The analysis does not establish biological causes or validate functional residue labels. All scoped C-alpha counts and majority scores from the referenced analysis reproduce; no general novelty claim is made.

This analysis compares chain A C-alpha distances in the representative calcium-free calmodulin structure 1CFD with the 25 coordinate models in its related 1CFC ensemble. For residues 5–75, 594 of 2,346 within-domain distances in 1CFD were more than 0.01 Å below the ensemble minimum, and 9 were more than 0.01 Å above the ensemble maximum. The corresponding counts were 50 and 1 of 2,016 for residues 82–147, and 0 and 0 of 4,686 between domains. The largest difference from an ensemble limit was 0.685028 Å for N-domain residues 38 and 63. At 0.10 Å tolerance, the combined counts were 379, 7, and 0 for the N domain, C domain, and between domains. These results apply to the recorded coordinate files and selections. The ranges from 25 archive models are not physical limits or confidence intervals; models are not independent physical samples or time points. The comparison does not test experimental restraint satisfaction or establish a cause for shorter distances. Independent arithmetic, synthetic boundary and rigid-motion checks, and offline reproducibility checks passed in the original work.

Across two related ubiquitin coordinate entries, each primary contact rule (C-alpha, C-beta, and minimum heavy-atom distance) shares six of seven top betweenness residues between entries. Within an entry, C-alpha and heavy-atom rules share three of seven in 2K39 and four of seven in 6V5D. This rank agreement differs from edge agreement and depends on graph construction: adding consecutive-residue backbone edges makes the C-alpha top sets identical across entries. The analysis uses all deposited models in chain A, author residues 1–70, nonlocal residue pairs separated by more than two positions, majority contacts with half-count ties, and fixed distance cutoffs. The entries share source experimental information, model counts are not independent samples or physical probabilities, and the graph results do not establish biological causes or residue function.

We sampled 100,000 uniform subsets of each size from 13 to 18 from the 176 human ubiquitin models in PDB 6V5D. For the CB–CG2 axis of isoleucine 13, 572 to 84 subsets per size reached order 0.69 in the primary fixed frame; counts in the second fitting frame were similar. Across 1,000 groups per size, every mean of the first 31 subset values and every mean of the largest 31 values was below 0.69 in both frames. The largest primary-frame means were 0.3913 and 0.4943, respectively. Zero of 1,000 groups per size and frame crossed the threshold for either mean statistic; each 95% Clopper–Pearson interval was [0, 0.003682]. These intervals describe Monte Carlo uncertainty conditional on the deposited coordinates and sampler. They do not include measurement error or multiple-comparison correction. Direct ranking by I13 order does not reproduce the article’s NOE-based selection, and the single-subset tail does not give a probability for its selected mean. The calculation does not estimate biological population weights, motion rates, or causes.

This analysis compares distance standard deviation (SD) between and inside the two domains of calcium-free calmodulin in PDB entries 1CFC and 1DMO. For each entry, it used all 25 or 30 deposited models, respectively, and matched between-domain residue pairs to inside-domain pairs with the nearest mean distance, using a maximum gap of 0.25 Å. Among the matched subsets, the between-domain median SD was 8.6674 times the reference median in 1CFC and 4.7671 times the reference median in 1DMO. The matching retained 24.82% and 17.97% of candidate pairs, respectively. The result describes spread across these deposited coordinate models; it does not establish molecular motion, a biological cause, or an inferential effect. The pairs and models are not independent, and the study reports no p-value or confidence interval. Sensitivity, reproduction checks, source hashes, methods, and limitations are in the main text and original package.

For chain A coordinate models from ubiquitin entries 1D3Z (10 models), 2K39 (116), and 6V5D (176), this analysis compares 8 Å majority-contact labels with contacts from equally weighted mean coordinates after rigid alignment. For 2,278 residue pairs, it tests every graph attainable by a positive distance limit on each mean structure. The minimum error counts for fits on residues 1–70 and 1–76 were 0 and 0 for 1D3Z, 2 and 2 for 2K39, and 2 and 3 for 6V5D. Thus, for these specified coordinates and alignments, no single distance limit recovers every target label for 2K39 or 6V5D; both 1D3Z fits recover all labels. The result is limited to three selected ubiquitin entries, equal model weights, the stated fitting procedure, and calculated float64 distances. Model counts are not independent biological observations or physical state probabilities, and the analysis does not establish experimental contact truth or a coordinate-rounding bound after refitting.

This calculation tests how bounded changes in model weights affect heavy-atom (HA) and side-chain (SC) contact labels for fixed ubiquitin coordinates in PDB entries 2K39 and 6V5D. At largest-to-smallest weight ratio R=1.10, two of five original HA differences and one of two SC differences retain positive-width distance intervals in [3.5, 5.0] Å for which labels are opposite for every permitted pair of entry-specific weight vectors. Across the six tested R values (1, 1.05, 1.10, 1.25, 1.50, and 2), the original eight pairs have 5, 5, 2, 1, 0, and 0 HA intervals and 2, 2, 1, 0, 0, and 0 SC intervals. For all eligible pairs at R=1.10, the counts are 105 HA and 83 SC. These are geometric labels conditional on deposited coordinate numbers and a mathematical weight constraint. The permitted weights have no experimental-fit constraint; the result does not establish physical state probabilities, model independence, biological causes, or experimental replication. The deposits share experimental information. Coordinate and model-selection uncertainty are not included.

For ubiquitin coordinate ensembles 2K39 and 6V5D, generalized Procrustes fits used C-alpha residues 1–70 or 1–76, and N-H coordinate order was compared with the two experimental S² lists in BMRB 6470. Both lists contained 61 shared residues for the primary comparison. RMSE and Pearson correlation selected different alignments for both ensembles; removing any one residue did not change these directions. Fixed scale sensitivity checks for 2K39 retained the RMSE direction. However, BMRB 6470 lacks list-to-model identities and experimental error values, and its recorded conditions (300 K, pH 4.7) differ from the PDB headers (308 K, pH 6.5). The analysis has no statistical confidence intervals and does not establish physical alignment accuracy, motion rates, populations, or biological causes. It reports numerical agreement for these files and proposes a future check using independently obtained, condition-matched N-H data.

In three selected calmodulin coordinate entries, all 1296 planned one-residue endpoint and axis-rule comparisons kept the calcium-free helix-angle ranges above the 1X02 ranges. The smallest gap was 0.6629351897280458 degrees for C/D in 1DMO versus 1X02. A later C/D check allowed each endpoint to change by up to two residues. Ranges overlapped in 39 of 1250 rules for 1CFC and 129 of 1250 for 1DMO. The rule counts describe analysis choices, not biological probabilities. Undirected axes changed the direction of range separation for G/H. These are results from selected archive coordinate models and specified rules. They do not establish calcium causality, protein dynamics, or a correct physical helix axis. The downloaded files are identified by their recorded URLs and SHA-256 checksums; the archive URLs serve current files rather than immutable historical snapshots.

Across 302 deposited ubiquitin NMR models from PDB entries 1D3Z, 2K39, and 6V5D, exhaustive single-model deletion changed the 70-residue correlation matrix more than the core-versus-full alignment mask in 1D3Z at a 0.03 Å variance floor: all 10 removals exceeded the full-data tail-mask change, by at least 4.74705 times. The corresponding counts were 13/116 for 2K39 and 5/176 for 6V5D. The seven 1D3Z removals meeting the 0.2 sign-change size cutoff affected 92 unique residue pairs; no 2K39 or 6V5D removal met it. Refitting alignment changed the matrices further, but the counts above the tail effect were the same with fixed full-data alignment. These are observations on selected deposited coordinates, not independent biological observations, confidence intervals, or evidence of physical interactions. Shared experimental data and refinement dependence limit comparisons between entries, and the analysis does not establish which correlation matrix is physically accurate.