biojs-vis-pdbviewer
Version:
A BioJS 2.0 component to view protein structures
134 lines (128 loc) • 4.71 kB
JavaScript
Clazz.declarePackage ("J.adapter.readers.more");
Clazz.load (["J.adapter.readers.more.ForceFieldReader"], "J.adapter.readers.more.MdTopReader", ["java.lang.Boolean", "JU.Lst", "J.adapter.smarter.Atom", "J.api.JmolAdapter", "JU.Logger"], function () {
c$ = Clazz.decorateAsClass (function () {
this.nAtoms = 0;
this.ac = 0;
this.$atomTypes = null;
this.group3s = null;
Clazz.instantialize (this, arguments);
}, J.adapter.readers.more, "MdTopReader", J.adapter.readers.more.ForceFieldReader);
Clazz.overrideMethod (c$, "initializeReader",
function () {
this.setIsPDB ();
this.setUserAtomTypes ();
});
Clazz.overrideMethod (c$, "checkLine",
function () {
if (this.line.indexOf ("%FLAG ") != 0) return true;
this.line = this.line.substring (6).trim ();
if (this.line.equals ("POINTERS")) this.getPointers ();
else if (this.line.equals ("ATOM_NAME")) this.getAtomNames ();
else if (this.line.equals ("CHARGE")) this.getCharges ();
else if (this.line.equals ("RESIDUE_LABEL")) this.getResidueLabels ();
else if (this.line.equals ("RESIDUE_POINTER")) this.getResiduePointers ();
else if (this.line.equals ("AMBER_ATOM_TYPE")) this.getAtomTypes ();
else if (this.line.equals ("MASS")) this.getMasses ();
return false;
});
Clazz.overrideMethod (c$, "finalizeReader",
function () {
this.finalizeReaderASCR ();
var atoms = this.asc.atoms;
var atom;
for (var i = 0; i < this.ac; i++) {
atom = atoms[i];
atom.isHetero = J.api.JmolAdapter.isHetero (atom.group3);
var atomType = this.$atomTypes[i];
if (!this.getElementSymbol (atom, atomType)) atom.elementSymbol = J.adapter.readers.more.ForceFieldReader.deducePdbElementSymbol (atom.isHetero, atom.atomName, atom.group3);
}
var atoms2 = null;
if (this.filter == null) {
this.nAtoms = this.ac;
} else {
atoms2 = new Array (atoms.length);
this.nAtoms = 0;
for (var i = 0; i < this.ac; i++) if (this.filterAtom (atoms[i], i)) atoms2[this.nAtoms++] = atoms[i];
}for (var i = 0, j = 0, k = 0; i < this.ac; i++) {
if (this.filter == null || this.bsFilter.get (i)) {
if (k % 100 == 0) j++;
this.setAtomCoordXYZ (atoms[i], (i % 100) * 2, j * 2, 0);
}}
if (atoms2 != null) {
this.asc.discardPreviousAtoms ();
for (var i = 0; i < this.nAtoms; i++) this.asc.addAtom (atoms2[i]);
}JU.Logger.info ("Total number of atoms used=" + this.nAtoms);
this.setModelPDB (true);
this.htParams.put ("defaultType", "mdcrd");
});
Clazz.defineMethod (c$, "getDataBlock",
function () {
var vdata = new JU.Lst ();
this.discardLinesUntilContains ("FORMAT");
var n = J.adapter.smarter.AtomSetCollectionReader.getFortranFormatLengths (this.line.substring (this.line.indexOf ("("))).get (0).intValue ();
var i = 0;
var len = 0;
while (true) {
if (i >= len) {
if (this.rd () == null) break;
i = 0;
len = this.line.length;
if (len == 0 || this.line.indexOf ("FLAG") >= 0) break;
}vdata.addLast (this.line.substring (i, i + n).trim ());
i += n;
}
return vdata.toArray ( new Array (vdata.size ()));
});
Clazz.defineMethod (c$, "getPointers",
function () {
var tokens = this.getDataBlock ();
this.ac = this.parseIntStr (tokens[0]);
var isPeriodic = (tokens[27].charAt (0) != '0');
if (isPeriodic) {
JU.Logger.info ("Periodic type: " + tokens[27]);
this.htParams.put ("isPeriodic", Boolean.TRUE);
}JU.Logger.info ("Total number of atoms read=" + this.ac);
this.htParams.put ("templateAtomCount", Integer.$valueOf (this.ac));
for (var i = 0; i < this.ac; i++) this.asc.addAtom ( new J.adapter.smarter.Atom ());
});
Clazz.defineMethod (c$, "getAtomTypes",
function () {
this.$atomTypes = this.getDataBlock ();
});
Clazz.defineMethod (c$, "getCharges",
function () {
var data = this.getDataBlock ();
if (data.length != this.ac) return;
var atoms = this.asc.atoms;
for (var i = this.ac; --i >= 0; ) atoms[i].partialCharge = this.parseFloatStr (data[i]);
});
Clazz.defineMethod (c$, "getResiduePointers",
function () {
var resPtrs = this.getDataBlock ();
JU.Logger.info ("Total number of residues=" + resPtrs.length);
var pt1 = this.ac;
var pt2;
var atoms = this.asc.atoms;
for (var i = resPtrs.length; --i >= 0; ) {
var ptr = pt2 = this.parseIntStr (resPtrs[i]) - 1;
while (ptr < pt1) {
if (this.group3s != null) atoms[ptr].group3 = this.group3s[i];
atoms[ptr++].sequenceNumber = i + 1;
}
pt1 = pt2;
}
});
Clazz.defineMethod (c$, "getResidueLabels",
function () {
this.group3s = this.getDataBlock ();
});
Clazz.defineMethod (c$, "getAtomNames",
function () {
var names = this.getDataBlock ();
var atoms = this.asc.atoms;
for (var i = 0; i < this.ac; i++) atoms[i].atomName = names[i];
});
Clazz.defineMethod (c$, "getMasses",
function () {
});
});