This command writes the current MODEL to a file in the selected format.
file can be a filename or a writeable file handle (see modfile.File()).
If you want to only write out a subset of the atoms, see selection.write().
'PDB' writes out files in the Protein Data Bank (PDB) format. Note that the isotropic temperature factor (Biso) field can be set by selection.energy() or model.write_data().
'MMCIF' writes out files in the new PDBx or Macromolecular Crystallographic Information File (mmCIF) format.
Note that if the model contains 100,000 or more atoms, it cannot be represented in traditional PDB format (as the atom serial numnber is a fixed 5 digits). MODELLER will work around this by using hybrid-36 notation, which is understood by many PDB readers. In general, however, mmCIF output is recommended for large systems.
The 'GRASP' format is the same as the 'PDB' format, except that it includes two special lines at the top of the file and the atomic radii and charges in the columns following the Cartesian coordinates of atoms. This format is useful for input to program GRASP, written by Anthony Nicholls in the group of Barry Honig at Columbia University [Nicholls et al., 1991]. For GRASP output, the atomic radii are needed. This usually means using the complete_pdb() script rather than model.read() to read in any original PDB file.
If the model contains any disulfide bridges, CONECT and SSBOND records for each bridge are included when writing out PDB files (note that CONECT records for other contacts, such as between HETATM residues, are not included). For mmCIF files, equivalent _struct_conn records are written.
A number of MODELLER-specific REMARKs are added when writing out PDB files (and are parsed when reading such files back in with model.read()). For mmCIF files, the same information is stored in the '_modeller' or '_modeller_blk' categories (shown in parentheses below):
The contents of model.remark are also added to PDB files verbatim.
If extra_data is given, its contents are also added to PDB or mmCIF files verbatim (so it must be formatted correctly for each format).
Example: See model.read() command.