mol-script.ts 3.8 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394959697
  1. import { MolScriptBuilder } from 'mol-script/language/builder';
  2. import { compile, QuerySymbolRuntime, DefaultQueryRuntimeTable } from 'mol-script/runtime/query/compiler';
  3. import { QueryContext, Structure, StructureQuery, CustomPropertyDescriptor } from 'mol-model/structure';
  4. import { readCifFile, getModelsAndStructure } from '../apps/structure-info/model';
  5. import { CustomPropSymbol } from 'mol-script/language/symbol';
  6. import Type from 'mol-script/language/type';
  7. import { parseMolScript } from 'mol-script/language/parser';
  8. import * as util from 'util'
  9. import { transpileMolScript } from 'mol-script/script/mol-script/symbols';
  10. import { formatMolScript } from 'mol-script/language/expression-formatter';
  11. import { StructureQualityReport } from 'mol-model-props/pdbe/structure-quality-report';
  12. import fetch from 'node-fetch';
  13. // import Examples from 'mol-script/script/mol-script/examples'
  14. // import { parseMolScript } from 'mol-script/script/mol-script/parser'
  15. // //import { compileAST } from 'mol-script/script/mol-script/compile';
  16. // for (const e of Examples) {
  17. // const expr = parseMolScript(e.value)[0];
  18. // console.log(e.name, util.inspect(expr, true, 10, true));
  19. // }
  20. // const exprs = parseMolScript(`(sel.atom.atom-groups
  21. // :residue-test (= atom.auth_comp_id ALA)
  22. // ;; ho ho ho
  23. // :atom-test (set.has { _C _N } atom.el)) ; comm
  24. // ;; this is a comment
  25. // ((hi) (ho))`);
  26. // ;; :residue-test (= atom.label_comp_id REA)
  27. const exprs = parseMolScript(`(sel.atom.atom-groups
  28. :residue-test (> pdbe.structure-quality.issue-count 0)
  29. :atom-test (= atom.el _C))`);
  30. const tsp = transpileMolScript(exprs[0]);
  31. // console.log(util.inspect(exprs, true, 10, true));
  32. console.log(util.inspect(tsp, true, 10, true));
  33. console.log(formatMolScript);
  34. console.log(formatMolScript(tsp));
  35. // //console.log(expr);
  36. const expr = MolScriptBuilder.core.math.add([1, 2, 3]);
  37. const compiled = compile<number>(expr);
  38. const result = compiled(new QueryContext(Structure.Empty));
  39. console.log(result);
  40. const CustomProp = CustomPropertyDescriptor({
  41. name: 'test_prop',
  42. isStatic: true,
  43. cifExport: { prefix: '', categories: [ ]},
  44. symbols: {
  45. residueIndex: QuerySymbolRuntime.Dynamic(CustomPropSymbol('custom.test-prop', 'residue-index', Type.Num), ctx => {
  46. const e = ctx.element;
  47. // console.log(e.element, e.unit.model.atomicHierarchy.residueAtomSegments.index[e.element])
  48. return e.unit.model.atomicHierarchy.residueAtomSegments.index[e.element];
  49. })
  50. }
  51. });
  52. DefaultQueryRuntimeTable.addCustomProp(CustomProp);
  53. DefaultQueryRuntimeTable.addCustomProp(StructureQualityReport.Descriptor);
  54. export async function testQ() {
  55. const frame = await readCifFile('e:/test/quick/1cbs_updated.cif');
  56. const { structure } = await getModelsAndStructure(frame);
  57. await StructureQualityReport.attachFromCifOrApi(structure.models[0], {
  58. PDBe_apiSourceJson: async model => {
  59. const rawData = await fetch(`https://www.ebi.ac.uk/pdbe/api/validation/residuewise_outlier_summary/entry/${model.label.toLowerCase()}`, { timeout: 1500 });
  60. return await rawData.json();
  61. }
  62. })
  63. let expr = MolScriptBuilder.struct.generator.atomGroups({
  64. 'atom-test': MolScriptBuilder.core.rel.eq([
  65. MolScriptBuilder.struct.atomProperty.core.elementSymbol(),
  66. MolScriptBuilder.es('C')
  67. ]),
  68. // 'residue-test': MolScriptBuilder.core.rel.eq([
  69. // MolScriptBuilder.struct.atomProperty.macromolecular.label_comp_id(),
  70. // 'REA'
  71. // ])
  72. 'residue-test': MolScriptBuilder.core.rel.inRange([CustomProp.symbols.residueIndex.symbol(), 1, 5])
  73. });
  74. expr = tsp;
  75. const compiled = compile<StructureQuery>(expr);
  76. const result = compiled(new QueryContext(structure));
  77. console.log(result);
  78. }
  79. testQ();