Creating Maps on SmallWorld

From DISI
Revision as of 13:24, 12 July 2019 by Jyoung (Talk | contribs)

Jump to: navigation, search

Written by Jennifer Young on April 10, 2019

Install: The developer need only do this once. For the user, skip to the Run section

Step 1: Create the /sw directory

SmallWorld requires a /sw directory which contains /anon and /maps The /sw directory is located on n-9-22 in

   /srv/db4

This directory contains:

The /anon directory contains 12 TB of pre-computed subgraphs

The /maps directory contains text files computed by SmallWorld that map your molecules of choice to the corresponding nodes in the anonymous graph index

Step 2: Find a location for the sw.jar file

Copy the sw.jar file to a reasonable location such as

   /opt/nextmove/sw.jar

Run

Important! Currently, we only run SmallWorld map generation on n-9-22

Step 1: Set SWDIR environment variable

Determine where the /sw directory with the /anon and /maps directory and if using csh, set the SWDIR variable using

   setenv SWDIR /srv/db4

If using bash, set the SWDIR variable using

   export SWDIR=/srv/db4

Step 2: Generate the map files

Run sw.jar to convert your smiles file into a map file.

(Optional: use the time command to see how long this takes)

(Optional: tee your results into a log file to check progress)

Here is an example where wait-ok.smi is your .smi file to map:

   (time java -jar /opt/nextmove/sw.jar map /mnt/nfs/ex9/work/smallworld/wait-ok.smi 
   -T /mnt/nfs/ex9/work/smallworld/scratch_sw 
   -o /mnt/nfs/ex9/work/smallworld/wait-ok.anon.map 
   -u /mnt/nfs/ex9/work/smallworld/wait-ok.anon.unmapped)
    |& tee /mnt/nfs/ex9/work/smallworld/wait-ok_log

This is a long command, so we can break down the different parts.

map The SMILES file you want to map follows the word "map"

-T This is the temp directory you want to write temp files to, rather than writing to /tmp

-o The final name of your output map file

-u Argument to keep track of the unmapped compounds.

Important: The unmapped compounds file may be used to generate extensions and is the same as the .anon.ext.map file mentioned in the SmallWorld documentation


For now, generate the map files in the

   /mnt/nfs/ex9/work/smallworld 

directory so they can be easily moved to the maps directory when the server is restarted

Processing large SMILES files in SmallWorld

If your SMILES file contains more than 1 billion molecules, I think it is better to split the input in to chunks of 1 billion to avoid waiting too long to get intermediate results in case the job fails for some reason.

Step 1: Split the input file into chunks of 1 billion smiles

   split -l 1000000000 <your-smiles-file>

Step 2: Process the chunks in SmallWorld using the command above

Write a helper script to process these using a loop, usually one or two at a time to be mindful of disk overload. An example is shown below. Just modify the path and pattern to match your files.

   #!/bin/csh
   foreach j(/mnt/nfs/ex9/work/smallworld/REAL_SPACE_400_to_500/x??)
       echo $j
       echo '(time java -jar /opt/nextmove/sw.jar map' $j '-T /mnt/nfs/ex9/work/smallworld/scratch_sw -o '${j}'.anon.map -u '${j}'.anon.unmapped) |& tee '${j}'_log'
       (time java -jar /opt/nextmove/sw.jar map $j -T /mnt/nfs/ex9/work/smallworld/scratch_sw -o ${j}.anon.map -u ${j}.anon.unmapped) |& tee ${j}_log
   end

Important: Make sure that the file you want to process only contains two columns. The first column contains the SMILES and the second column contains the ID. If this is not the case, you need to use

   awk '{print $1"\t"$2}'  > your_file_smi_ID

to get only these columns before processing with SmallWorld.

Step 3: Copy to SSD and Concatenate

It is much faster and avoids overloading the disk if you concatenate the results on an SSD rather than a regular hard drive. I have been using /nfs/db5/jyoung as a place to concatenate the files. You can use another helper script to perform that concatenation one file at a time.

   #!/bin/csh
   foreach j(*.anon.map)
       ls $j
       cat $j >> enamine_private_400_to_500_combined
   end

Step 2.5 Incremental Updates

There are several ways to perform incremental updates for SmallWorld maps which are equivalent.

If you just have a new SMILES file input.smi you wish to add to an already existing map named <already_existing>.anon.map, then you can use the sw map command with the --append option provided in the SmallWorld Java Command Line Interface (CLI)

First, if you have not already done so, set an alias "sw" to the SmallWorld Java file as shown below: Make sure the SWDIR environment variable is also set

   setenv SWDIR /srv/db4
   alias sw 'java -jar /opt/nextmove/sw.jar'

Now you can run the sw map command using --append

   sw map input.smi --append --o <already_existing>.anon.map

Another method for performing incremental updates is discussed below. A worked example is located in the SmallWorld version 4 documentation on pages 17-18 in the Incremental Updates Section

What if you already created a .anon.map file from your new smiles, and simply want to combine the two separate .anon.map files into a single .anon.map file?

Make a copy and modify the bash script below with the name of your new map file

   /nfs/ex9/work/smallworld/bash_merge_maps

This script simply combines two map files that have already been created

Replace the OLD_MAPFILE environment variable with the old map

Replace the ADD_MAPFILE environment variable with the map you wish to incrementally add

Replace the NEW_MAPFILE environment variable with the name you want for the new combined map

Updating a new map file in this way does not automatically update the .anon.map.cfg file with the new total number of molecules mapped and other statistics. The create_cfg.py program will automatically generate a new .anon.map.cfg file for the combined map

The create_cfg.py program is located in

   /mnt/nfs/ex9/work/smallworld/create_cfg.py

The create_cfg.py script will automatically generate a cfg file for a map that is composed of two existing SmallWorld Maps. Run as follows:

   python create_cfg.py first second combined

Note: first is first.anon.map second is second.anon.map combined is combined.anon.map

The program will look for the .anon.map files and the .anon.map.cfg for the first and second before proceeding

The new_cfg.py program is located in

   python /mnt/nfs/ex9/work/smallworld/new_cfg.py <combined cfg file name>

This Python program will find all the cfg folders in a directory and combine them into a single combined cfg file

Step 3: How to Make the New Database Available. Stop the server first!

First, stop the server. You do not need to restart Tomcat every time you want to restart SmallWorld. The TomCat GUI is available at

   http://10.20.9.22:8080/

You need the tunnel to access the private network or this link will not work! Also, the admin username and password is required.

In the Manager App, in the Applications section, use the Stop and Start buttons for the sw-ws app to restart the SmallWorld server

Step 4: Move the Map File into the /maps directory

If you followed the example above, your .anon.map file should be located in

   /mnt/nfs/ex9/work/smallworld

Now, simply move the map file into the "/sw/maps" directory so it can be recognized by SmallWorld.

   mv <your-file-name>.anon.map /mnt/nfs/ex9/work/smallworld/maps

Or, simply

   mv <your-file-name>.anon.map maps


Step 5: Restart the Server

Again, go to the TomCat manager under applications and click Start once you are finished moving around the map files.