Skip to content

Commit

Permalink
remove the temporary pdb file written for propka
Browse files Browse the repository at this point in the history
  • Loading branch information
stefdoerr committed Jan 3, 2022
1 parent c48a97a commit 1818a51
Showing 1 changed file with 6 additions and 4 deletions.
10 changes: 6 additions & 4 deletions pdb2pqr/main.py
Original file line number Diff line number Diff line change
Expand Up @@ -509,13 +509,15 @@ def run_propka(args, biomolecule):
lines = io.print_biomolecule_atoms(
atomlist=biomolecule.atoms, chainflag=args.keep_chain, pdbfile=True
)
with NamedTemporaryFile("wt", suffix=".pdb", delete=False) as pdb_file:
with NamedTemporaryFile("wt", suffix=".pdb") as pdb_file:
for line in lines:
pdb_file.write(line)
pdb_path = pdb_file.name
parameters = pk_in.read_parameter_file(args.parameters, Parameters())
molecule = MolecularContainer(parameters, args)
molecule = pk_in.read_molecule_file(pdb_path, molecule)

parameters = pk_in.read_parameter_file(args.parameters, Parameters())
molecule = MolecularContainer(parameters, args)
molecule = pk_in.read_molecule_file(pdb_path, molecule)

molecule.calculate_pka()

# Extract pKa information from PROPKA
Expand Down

0 comments on commit 1818a51

Please sign in to comment.