Changes between Version 1 and Version 2 of CreateMicelle


Ignore:
Timestamp:
Feb 9, 2011, 5:49:53 PM (14 years ago)
Author:
FrederikHeber
Comment:

added parsing of .potentials file and typo fixes

Legend:

Unmodified
Added
Removed
Modified
  • CreateMicelle

    v1 v2  
    55
    66=== First Micelle ===
    7 First, we will specify an empty '''TREMOLO''' data file (''-i test.data''), we give '''PDB''' and '''TREMOLO''' as output formats (''-o pdb tremolo''). Then we load the tenside molecule (''-l tensid.data''). We select it by its name, which is automatically set to from the filename (''--select-molecule-by-name tensid''). And then we call the !CreateMicelleAction with the following parameters (''--create-micelle 200 --position "0,0,0" --radius 10''):
     7First, we will specify an empty '''TREMOLO''' data file (''-i test.data''), we give '''PDB''' and '''TREMOLO''' as output formats (''-o pdb tremolo''). Then we first parse the potentials file (''--parse-tremolo-potentials tensid.potentials''), then load the tenside molecule (''-l tensid.data''). ^(Note: The potentials file gives us information about which type represents which element, e.g. HLA2 actually means hydrogen or H. So, this part of the command is only necessary if your Type column in tensid.data does not contain valid element symbols.)^ We select the tensid molecule by its name, which is automatically set to from the filename (''--select-molecule-by-name tensid''). And then we call the !CreateMicelleAction with the following parameters (''--create-micelle 200 --position "0,0,0" --radius 10''):
    88 * 200 points on the surface of the sphere equals 200 resulting tenside molecules (MoleCuilder will state 400 because the ionic compound is regarded as a distinct molecule)
    99 * The sphere is centered at (0,0,0)
     
    1111
    1212{{{
    13 ~/workspace_C/molecuilder/debug64/src/molecuilder -v 3 -i test.data -o pdb tremolo -l tensid.data --select-molecule-by-name tensid --create-micelle 200 --position "0,0,0" --radius 10
     13~/workspace_C/molecuilder/debug64/src/molecuilder -v 3 -i test.data -o pdb tremolo --parse-tremolo-potentials tensid.potentials -l tensid.data --select-molecule-by-name tensid --create-micelle 200 --position "0,0,0" --radius 10
    1414}}}
    1515
     
    1717The second call is as follows:
    1818{{{
    19 ~/workspace_C/molecuilder/debug64/src/molecuilder -v 3 -i test.data -o pdb tremolo -l tensid.data --select-molecule-by-name tensid --create-micelle 200 --position "100,0,0" --radius 10
     19~/workspace_C/molecuilder/debug64/src/molecuilder -v 3 -i test.data -o pdb tremolo --parse-tremolo-potentials -l tensid.data --select-molecule-by-name tensid --create-micelle 200 --position "100,0,0" --radius 10
    2020}}}
    2121
     
    2323
    2424=== Filling empty space with water molecules ===
    25 Now, we run into a basic problem of the '''PDB''' format. The atomic ids only go up to 99999, i.e. 5 digits. Filling the void space with water however will result in roughly 400.000 atoms. Hence, we set the output format to '''XYZ''' (''-o xyz'') instead. First, we have to specify a domain that we fill. Therefore, we call only we now call !AddEmptyBoundaryAction (''--boundary "20,20,20"'') which causes every atom to be last 20 Angstroem away from the boundary. Then, all is set for calling !FillVoidWithMoleculeAction (''--fill-void water.xyz --!MaxDistance -1 --distances "3.1, 3.1, 3.1"  --lengths "2.1, 0., 0." --!DoRotate 0'') with the following parameters:
     25Now, we run into a basic problem of the '''PDB''' format. The atomic ids only go up to 99999, i.e. 5 digits. Filling the void space with water however will result in roughly 400.000 atoms. Hence, we set the output format to '''XYZ''' (''-o xyz'') instead. First, we have to specify a domain that we fill. Therefore, we now call !AddEmptyBoundaryAction (''--boundary "20,20,20"'') which causes every atom to be last 20 Angstroem away from the boundary. Then, all is set for calling !FillVoidWithMoleculeAction (''--fill-void water.xyz --!MaxDistance -1 --distances "3.1, 3.1, 3.1"  --lengths "2.1, 0., 0." --!DoRotate 0'') with the following parameters:
    2626 * The water molecule (the ''filler'') is obtained from parsing ''water.xyz''
    2727 * There is no maximum distance given (e.g. to some surface)
    28  * The distances in between the water molecules is 3.1 Angstroem, relating to a density of 1g/cm^3 (try to calculate with '''units''' yourself).
     28 * The distances in between the water molecules is 3.1 Angstroem, relating to a density of 1g/cm^3^ (try to calculate with '''units''' yourself).
    2929 * The lengths here are: 2.1 Angstroem distance to every other already present atom found. The two remaining numbers are for randomizing atomic and molecule's central position.
    3030 * Finally, we do not rotate randomly the water molecule.
    3131
    3232{{{
    33 ~/workspace_C/molecuilder/debug64/src/molecuilder -v 1 -i test.data -o xyz --boundary "20,20,20" --fill-void water.xyz --MaxDistance -1 --distances "3.1, 3.1, 3.1"  --lengths "2.1, 0., 0." --DoRotate 0
     33~/workspace_C/molecuilder/debug64/src/molecuilder -v 1 --parse-tremolo-potentials -i test.data -o xyz --boundary "20,20,20" --fill-void water.xyz --MaxDistance -1 --distances "3.1, 3.1, 3.1"  --lengths "2.1, 0., 0." --DoRotate 0
    3434}}}
    3535