SIMONAPreprocessor
POEM++ - First Steps

How to create your first molecule simulation.

This is just an example to get to know how this section works and what it does roughly. For more detailed information about all functions of POEM++ see here.

  • Get a molecule
    In order to simulate a molecule you need an xml-datasheet. This sheet should be created by the SIMONA Preprocessor. The preprocessor needs an input file with the structure of a molecule. For example you can get various types of known molecules from PDB. Make sure you chose the *.pdb file for download. Often in these files is more then one molecule. Select one for example by using PyMol.
  • Startwindow
    Now start the SIMONA Preprocessor by typing
    1 python simonagui.py
    into the shell (you have to be in the directory where simonagui.py is located). In opening window you can choose of which type your simulation should be. Select POEM++ for a protein simulation.

    1s_Mainpage.png

  • Input files
    To read in the data of the molecule, select Configuration at the menu on the left hand side. Type the path of the input-file in the box beside Input configuration. You may also browse the path by clicking the ... button beside the box. Also type the output-path in the box beside Output xml. To receive a file readable by SIMONA the file extension has to be *.xml.

    1s_Configuration_Input.png

  • Environment
    In this section the temperature has to be set. Temperature defines how free to move the molecule is. Either select Constant temperature or Geometric annealing. If constant temperature is selected only one temperature is to be entered in the box on the right. The temperature scale is in Kelvin. For a geometric annealed temperature a Start temperature and an End temperature are required. A standard simulation starts with a rather high, e.g. 500, while End temperature should be rather low, e.g. 5. But feel free to experiment with other settings.

    1s_Configuration_Environment.png

  • Simulation Settings
    For the outcome of a simulation the most important thing is the number of Simulation steps. The more you take, the better the result will be. But the needed calculation time will grow as well. Don't set more than 100.000 steps for your first try. Random seed should be selected.

    1s_Configuration_Simulation.png

  • Forcefields
    Now we proceed to the next section of the left menu. Use one of the predefined forcefields for your first simulation, for example PFF02 [1] .

    1s_Forcefields.png

  • Moves
    In the third section of the side menu you have to define which moves your molecule should be able to do. For a single molecule you only need Dihedral rotations . Uncheck Rigid body translations and Rigid body rotations . Check Dihedral rotations and set Delta [°] to 10.00. In the selection column select Mainchain dihedrals.

    1s_Moves.png

  • Output
    The fourth section of the side menu provides options for output settings. Select Std-,Trajectory-, Energy- and RMSD-output. Stepmod should be 1000. The Output filename of the trajectory has to have the ending *.pdb or *.vcf. A *.vcf file contains coordinates only, so is smaller, but can't be read by visualization tools. A *.pdb file contains additional information, so can be easily visualized.

    1s_Output.png

  • Generate XML
    Last step on the way to your *.xml is to let the Preprocessor do its work. On the top menu select Generate XML and wait until the output file has been created. The output file will be saved at the position you selected earlier. It is shown in the side menu section Preprocessor output.
  • Start SIMONA
    So far you have created an input file for SIMONA. To run SIMONA with your molecule and your settings type
    1 poempp <outputfile.xml>
    in the shell. <outpufile.xml> is here the outputfile you created with SIMONA Preprocessor. Make sure you are in the right directory. SIMONA will start its work and after some time creates the requested files.
  • Result
    To see the outcome use a visualization tool. In the trajectory file you can see the evolution of the simulation. The last step of the calculation is saved in *_snap.pdb.