Skip to content
Snippets Groups Projects
Commit 4111d22e authored by Marco Biasini's avatar Marco Biasini
Browse files

properly mark ligands

All residues of chains that don't have a SEQRES entry and
are not type water are marked as ligand
parent d80effd5
No related branches found
No related tags found
No related merge requests found
......@@ -1080,6 +1080,15 @@ void MMCifParser::OnEndData()
if (edm_it->second.seqres.length() > 0) {
seqres_.AddSequence(seq::CreateSequence(css->first.GetName(),
edm_it->second.seqres));
} else if (edm_it->second.type!=mol::CHAINTYPE_WATER) {
// mark everything that doesn't have SEQRES as ligand and isn't of type
// water as ligand
mol::ChainHandle chain=css->first;
mol::ResidueHandleList residues=chain.GetResidueList();
for (mol::ResidueHandleList::iterator
i=residues.begin(), e=residues.end(); i!=e; ++i) {
(*i).SetIsLigand(true);
}
}
} else {
LOG_WARNING("No entity description found for atom_site.label_entity_id '"
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment