PDB IO

PDB IO

SASMOL: Copyright (C) 2011 Joseph E. Curtis, Ph.D.

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see <http://www.gnu.org/licenses/>.

class sasmol.pdb_io.PDB

Bases: object

check_error(error)

check_for_all_zero_columns(coor, frame=0)

Make sure there are no two all zero columns in coordinates This is needed for the alignment module

create_conect_pdb_lines()

Output stored conect information in PDB record format. Indices are converted from those read in - original_index - to those used in output. Format of CONECT line: Record name CONECT followed by list of atom indexes in 6 character columns. First is the base atom, following atoms are those connected to it.

element_filter()

This function filters the PDB file to sraighten out various things (add to this as you go along).

Filter element list…

get_element(name, resname)

Get elment from charmm27 atom name list plus extras from periodic table and resolve name conflicts

get_resnames()

This method holds names of residues to use to set moltype. Based on Charmm 27 naming.

initialize_children()

This initializes the following “children” with their masks already defined to the “parent” molecule

names() : names_mask() resnames() : resnames_mask() resids() : resids_mask() chains() : chains_mask() segnames() : segnames_mask() occupancies() : occupancies_mask() betas() : betas_mask() elements() : elements_mask()

The objects on the left contain the unique values and the objects on the right contain the masks that have the indices to extract the information for each unique value from the parent molecule.

print_error(name, my_message)

read_pdb(filename, **kwargs)

This method reads a PDB file.

write_pdb(filename, frame, flag, **kwargs)

This method writes the PDB file