|This code has been put under version control in the project Pymol-script-repo|
LiSiCA is a software for 2D and 3D ligand based virtual screening. It uses a fast maximum clique algorithm to find two- and three- dimensional similarities between reference compound and a database of target compounds in Mol2 format. The similarities are expressed using Tanimoto coefficients and the target compounds are ranked accordingly. LiSiCA is developed and maintained at National Institute of Chemistry, Slovenia.
LiSiCA Plugin is free for Academic (NON-COMMERCIAL) use. However, the users have to acquire an activation key. Please see Licensing section below for further details.
LiSiCA Plugin Features
- Graphical User Interface which facilitate the choice of parameters for LiSiCA.
- Results displayed according to the ranking in similarities as measured by Tanimoto coefficients.
- Structural similarity between molecules visualized using PyMOL viewer.
The plugin should work on Windows and all variety of Unix like systems. Most of the tests were performed on Windows 7, Ubuntu 14.04 and Linux Mint. To work properly, the plugin requires Python with Tkinter and PyMOL. The plugin was mainly tested with PyMOL 1.7.x.
This plugin is ready "out-of-box" for both Linux and Windows users through the project Pymol-script-repo.
The lisica.py script initializes the installation of the plugin. It performs the following functionalities:
- Downloads and installs all the required files. On successful installation, a directory named .lisicagui is downloaded and saved in the home directory. This folder contains the executables, log files, icon files, python modules etc. Also, a text file named ".insilab-license.txt" will be created in the home directory.
- Prompts user for Activation on installation. (See section: Licensing).
Once the files are properly installed and the license activated, LiSiCA plugin version will be ready to use.
Important note: For complete installation of the plugin, a stable INTERNET connection is required at the user end. Note for Linux users : It is advisable to start PyMOL from the terminal using sudo pymol command.
The plugin window has four tabs.
Both the reference and the target compound files need to be in the Tripos mol2 format. The reference file should contain only one (reference) compound; the target file may contain many compounds. If the reference file contains more than one compound, the first molecule is used as a reference. The molecules in the target file having the same name are considered as different conformers of the same molecule. By default only the best-scoring (by Tanimoto coefficient) conformer will be shown in the final output for the 3D screening option.
In the mol2 input files, the name of the molecule (ZINC ID) must be specified under the line for @<TRIPOS>MOLECULE tag. For example,
@<TRIPOS>MOLECULE ZINC73655097 46 48 0 0 0 SMALL USER_CHARGES @<TRIPOS>ATOM ...
LiSiCA checks similarities based on the mol2 atom types. Hence SYBYL atom types have to be specified in the mol2 files. For example, under the @<TRIPOS>ATOM tag, each atom specification should include the SYBYL atom type, as in:
1 C1 -0.0647 1.4496 -0.0592 C.3 1 <0> -0.167
In the above line from a mol2 file, in bold, in column 5, the SYBYL atom type is specified.
In 2D screening, the option Maximum Allowed Shortest path Difference corresponds to the maximum allowed difference in shortest-path length between atoms of the two compared product graph vertices. Lesser values correspond to a more rigorous screening. By default this value is unit bond.
In 3D screening, the option Maximum Allowed Spatial Distance Difference corresponds to the maximum allowed difference in distances between atoms of the two compared product graph vertices. Lesser values correspond to a more rigorous screening. By default this value is 1 Å. The Number of Conformations option corresponds to the maximum number of outputted files of one molecule in different conformations and is to be used only for 3D screening.
According to the value of Number of highest ranked molecules to be written to the output (say W), LiSiCA will create mol2 files of that many (W) highest scoring target molecules with a comment section at the end of the file where the matching atom pairs are displayed. This value is by default 100. The resulting mol2 files are written into a time-stamp directory in the folder specified in Save results in:. By default this folder is the user's home directory. Also, a text file named lisica_results.txt with the target molecules with (in the descending order of) Tanimoto coefficients is written to the time-stamp folder.
The Number of CPU cores to be used allows selection of CPU threads used for LiSiCA. By default, it tries to detect the number of CPUs available. The Consider Hydrogen options lets the user to choose if the hydrogen atoms are to be considered for the calculation of the similarity using the maximum clique algorithm. By default, hydrogen atoms are not considered in finding the largest substructure common to the reference and target molecules, so as to obtain faster results.
Load Project Tab
The plugin also has a feature to load saved results. On the Load Project tab, the user can choose the directory with the saved results (mol2 files of each target and the reference) and the lisica_results.txt file. When the load button is clicked, the results will be loaded onto the output tab and the PyMOL Viewer window.
In the output tab, there are two listboxes:
- One contains ZINC ID and Tanimoto Coefficients of target molecules in the decreasing order of the Tanimoto coefficient values.
Any single target molecule can be selected on this listbox using a mouse click or using up/down arrow keys. The selected target molecule is displayed with the reference molecule on the PyMOL viewer window.
- Depending on the target molecule chosen on the first listbox, the corresponding atoms from reference and the target molecules are displayed on the other listbox.
Any single pair of corresponding atoms can be selected on this listbox using a mouse click or using up/down arrow keys. The selected pair is highlighted on the PyMOL viewer window.
For 2D Screening, the two molecules (the reference and the selected target) are visualized side by side on the PyMOL viewer screen.
For 3D Screening, the 3D structures of the two molecules (the reference and the selected target) are superimposed on one another to visualize the similarity on the PyMOL viewer screen.
The about tab gives License, Version and Contact information. To deactivate LiSiCA, click on the deactivate button on the About tab. The users can get information on new updates if available on the About tab.
If you are using LiSiCA in your work, please cite:
S. Lesnik, T. Stular, B. Brus, D. Knez, S. Gobec, D. Janezic, J. Konc, LiSiCA: A Software for Ligand-Based Virtual Screening and Its Application for the Discovery of Butyrylcholinesterase Inhibitors, J. Chem. Inf. Model., 2015, 55, 1521–1528.
If you receive the error <class '_tkinter.TclError'>: can't find package treectrl during the LiSiCA installation process, you should fix this by installing tktreectrl: sudo apt-get install tktreectrl.
LiSiCA software is copyrighted by :
- Janez Konc
- National Institute of Chemistry
- Laboratory for Molecular Modeling
- Hajdrihova 19
- SI-1000 Ljubljana
Important note: For using LiSiCA plugin, an Activation Key is required. To obtain the key for NON-COMMERCIAL usage of LiSiCA, email to email@example.com providing your title, full name, institutuion, laboratory, supervisor name. Also, kindly state your intended usage of the software.
The terms stated in the following agreement apply to all files associated with the software unless explicitly disclaimed in individual files.
LiSiCA SOFTWARE LICENSE AGREEMENT 1. Grant Of Limited License; Software Use Restrictions. The programs received by you will be used only for NON COMMERCIAL purposes. This license is issued to you as an individual. For COMMERCIAL use of the software, please contact Janez Konc for details about commercial usage license agreements. For any question regarding license agreements, please contact: Janez Konc National Institute of Chemistry Laboratory for Molecular Modeling Hajdrihova 19 SI-1000 Ljubljana Slovenia. 2. COMMERCIAL USAGE is defined as revenues generating activities. These include using this software for consulting activities and selling applications built on top of, or using this software. Scientific research in an academic environment and teaching are considered NON COMMERCIAL. 3. Copying Restrictions. You will not sell or otherwise distribute commercially these programs or derivatives to any other party, whether with or without consideration. 4. Ownership of Software. You will not obtain, and will not attempt to obtain copyright coverage thereon without the express purpose written consent of Janez Konc. 5. IN NO EVENT SHALL THE AUTHORS OR DISTRIBUTORS BE LIABLE TO ANY PARTY FOR DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE USE OF THIS SOFTWARE, ITS DOCUMENTATION, OR ANY DERIVATIVES THEREOF, EVEN IF THE AUTHORS HAVE BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 6. THE AUTHORS AND DISTRIBUTORS SPECIFICALLY DISCLAIM ANY WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, AND NON-INFRINGEMENT. THIS SOFTWARE IS PROVIDED ON AN "AS IS" BASIS, AND THE AUTHORS AND DISTRIBUTORS HAVE NO OBLIGATION TO PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.