Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[DNM] prelim fixes on ns2b3 misalignment #43

Open
wants to merge 3 commits into
base: main
Choose a base branch
from

Conversation

JenkeScheen
Copy link
Contributor

No description provided.

return complex
# first set indexing to an unphysically high number
original_index = []
residue_N = 100000
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would make this number max(residue_numbers + 1) to avoid some kind of cutoff artefact in the extremely unlikely event someone passes in a megaprotein, probably unnecessary but good practice.

the crystallographer may have come up with. BioPython can do this but there are some
protections built in against hard re-indexing multi-chain proteins. By first setting
the indexing super high (starting at 100,000) and then re-indexing starting from 1 we can
circumvent these protections. For more details see https://github.com/biopython/biopython/pull/4623
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We can install from git+ source if you want this functionality in the aforementioned biopython PR. Only issue is we would then be blocked on release upstream until biopython merges and does a release.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Perhaps its just better to use this hack.


def load_pdb(self):
"""
Loads an input PDB file
"""
complex = PDBParser(QUIET=False).get_structure("COMPLEX", self.path_to_pdb_file)

self.check_validity(complex)
self.reset_complex_sequence(complex)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
self.reset_complex_sequence(complex)
complex, _ = self.reset_complex_sequence(complex)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants