martimunicoy / peleffy-benchmarks

A compilation of benchmarks to validate the parameterization of peleffy (PELE Force Field Yielder).
MIT License
0 stars 0 forks source link

Improve functions to call PELE's minimization #9

Closed martimunicoy closed 4 years ago

martimunicoy commented 4 years ago

There are several improvements and bug-fixes that we need to solve on the functions and classes that call PELE's minimization.

martimunicoy commented 4 years ago

Improvements list:

martimunicoy commented 4 years ago

Regarding the use of offpele's API, below is an example of how to use it:

from offpele.topology import Molecule
from offpele.template import Impact
from offpele.solvent import OBC2

# Create representation of a particular molecule
molecule = Molecule(path_to_pdb_file)

# Generate its rotamer library
rotamer_library = offpele.topology.RotamerLibrary(molecule)
rotamer_library.to_file(rotamer_library_output_path)

# Generate its parameters and template file
molecule.parameterize(forcefield, charges_method=charges_method)
impact = Impact(molecule)
impact.write(impact_output_path)

# Generate its solvent parameters
solvent = OBC2(molecule)
solvent.to_json_file(solvent_output_path)

Besides, the DataLocal's output hierarchy can be constructed through the function handle_output_paths() from main.py

laumalo commented 4 years ago

The following PDB files: benzene.pdb and toluene.pdb present the following error when computing the PELE's minimization.

what(): /!\ Error in MacroBuilder::createMacro: /!\ Error in ImpactTemplateReader::processRESXLine: Incorrect RESX line (wrong format) in file benz!!!! And the same error for tolz.