CCG SVL Exchange

An automated method for generating pharmacophore models from protein binding pockets

Installation:

1. Extract the package named MOEWeb_QuickDock to a temporay folder.
   Follow the instructions below to install the extracted files.
2. Extract the archives named dock.svl under $MOE/web on the MOEWeb server.
3. Download JSMOL viewer from the following link and extract under $MOE/web on the MOEWeb server. https://sourceforge.net/projects/jsmol/
4. Install the svl script from the 'svl' folder of the extracted files to the following location, both on MOEWeb server, and your MOE desktop client. If folders shown below do not exist, please create them.

   Filename	Location
   dock_ui.svl	$HOME/moefiles/svl/run
   moeweb_quickdock_util.svl	$HOME/moefiles/svl
   soap_web_fileutil.svl	$HOME/moefiles/svl

5. Place the menu file menu_extra_moewebdockutil from the 'menu' folder of the extracted files under $HOME/moefiles/menu
6. Restart the MOEWeb server by explicitely loading the soap_web_fileutil.svl, for example by adding the following line to MOEWeb configuration file ($MOE/web/moewebservice.cfg on Windows or equivalent on Linux):
   -load /home/$USERNAME/moefiles/svl/soap_web_fileutil.svl
   
7. Browse to the main page of the MOEWeb to see the icon for the newly added application named QuickDock
   You may need to force a refresh of the MOEWeb page, the first time after upgrade. This is often done by pressing ctrl+F5 or holding ctrl and pressing the refresh button on the browser.

Usage:

1. Browse to the main page of the MOEWeb server, locate the icon for the newly added application 'QuickDock' and click on it to start the application.
2. Browse to locate a *.sdf file containing the ligand, or use the sketcher to sketch the input molecule. You may also paste the ligand in *.mol format into the skether.
3. Under the Docking Options choose one of the targets prepared by the MOEWeb administrator.
4. Press MOE-Dock to start the docking.
5. Once docking is finished, use the links from the results page to download individual files.

Adding new targets:

1. Start desktop MOE which has the extensions from Step 2 of the installation.
2. Open the target structure from any source and prepare it as desired.
3. Open the Docking application via MOE | Compute | Dock .
4. Configure the docking as desired.
5. Once ready to run, press the 'Batch...' button and save the batch file.
   Ensure to use a descriptive name.
   This file includes the structure of the target, cocrystallized ligand (for comparison), and all the required settings.
6. Copy this file to the MOEWeb server under $MOE/web/dock/scripts.
7. Alternatively, open the panel via MOE | Extra | MOEWeb Dock File Transfer and use it to upload the batch script to the MOEWeb server.
8. Browse to QuickDock application page. Refresh the page if needed by pressing ctrl+F5. The new target name must appear under the list.

Sketch ligands or load from *.PDB/*.SDF files to perform conformational search.

Installation:

1. Extract confsearch.zip under $MOE/web
2. Download JSMol viewer from the following link and extract it under $MOE/web:
   http://jsmol.sourceforge.net/
   
3. Restart MOEWeb. A new icon titled Conformational Search appears under MOEWeb starting page.

This SVL code allows to create a sequence consensus logo.
It implements the approach published in:
T.D.Schneider and R.M.Stephens
(1990) Nucl.Acids Res. 18:6097-100
doi: 10.1093/nar/18.20.6097
Installation:

1. Copy the file to $HOME/moefiles/svl/run, where $HOME is your home directory

Usage:

1. In the sequence editor select the chains to create a consensus logo for.
2. Open an SVL command line and type: run 'consensus_logo.svl'
3. In the Consensus Logo panel press Apply.
4. To generate an image of the logo press the button Export.

Description:
This package will predict how molecules in MOE (or in a MOE Database) might be metabolized by Cytochrome P450s, aka CYPs. Specifically, the following properties are evaluated:

• Selectivity -- Using two different QSAR methods, predict which of the five main human isoforms of CYP (1A2, 2C9, 2C19, 2D6, and 3A4) might interact with a provided ligand
• Reactivity -- Using two different QSAR methods, predict whether a provided ligand might act as a Substrate or an Inhibitor for one of the five main human CYP isoforms
• Sites of Metabolism -- Using a pharmacophoric approach, predict which atom(s) on a provided ligand might be metabolized by one of three reactions -- Aliphatic Oxidation, Aromatic Oxidation, and N-Dealkylation -- catalyzed by the two most prevalent CYP isoforms -- 2D6 and 3A4.

To install, copy cyp_predictor.svl and the CYP_Files directory to $HOME/moefiles/svl. To launch CYP Predictor, double-click on cyp_predictor.svl from the MOE File | Open panel.

Description:
The tool performs principal component analysis of atomic coordinates on trajectory data obtained by molecular dynamics and show the motion of molecules along each pricipal component axis in MOE window. Please see pdf for installation directions and a tutorial.

This tool has been designed to bring together a variety of clustering methods which share a common framework.
Here is an overview of what this tool can do:

• Jarvis-Patrick clustering
• K-Means Map-Reduce clustering
• Parallelized K-Means clustering for large datasets
• Hierarchical clustering (Single-, complete- and average-linkage
• Fingerprint cluster model creation
• Cluster assignment of an independent dataset with a pre-existing fingerprint cluster model
• Clustering of subsets (selected entries only, random, first N, or descriptor-based criteria)
• Creation of reusable similarity matrices for exploratory analysis (ASCII and binary format)

Please see pdf for installation directions, a description of the application and a tutorial

This submission is an implementation of:
J.Dunbar, A.Fuchs, J.Shi, C.M.Deane (2013) PEDS 26(10):611-20
https://doi.org/10.1093/protein/gzt020

Installation:

1. Copy abangle.svl to $HOME/moefiles/svl/, where $HOME is your home directory
2. Restart MOE to load the file

Usage:

1. The SVL code implements 6 descriptors that match the ABangle descriptors: ABangle_dc, ABangle_HL, ABangle_HC1, ABangle_HC2, ABangle_LC1, ABangle_LC2.
   See the Dunbar paper for definition of these descriptors.
2. The SVL code also implements a graphical mode that renders the vectors. This only works if exactly one Fv is present in MOE.
   Call this SVL function:
   ABangle [Chains [], [annotate: 1, render: 1, render_color: 0x00ffff]]

This submission contains an implementation of:

Baell JB, Holloway GA. New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays.
J Med Chem. 2010 Apr 8;53(7):2719-40.
PubMed PMID: 20131845.
http://dx.doi.org/10.1021/jm901137j

Installation:

1. Copy q_pains.svl to $HOME/moefiles/svl, where $HOME is the users home directory
2. Restart MOE to load the new file

Notes:

• The directory must be created if it does not already exist.
• $HOME can be determined at the SVL Commands windows. Type:
  svl> HOME

Usage:

1. In the Database Viewer select the Compute | Descriptors menu. The new descriptors will be added to the list. If you type PAINS in the "Filter" box at the bottom of the panel, then just these descriptors will be shown. Select the descriptor that you want to calculate and click the OK button.
2. To annotate a database of compounds with the rules that match compounds:
   svl> db_PAINS ['in.mdb', 'mol']
3. The PAINS will also be added to the ligand properties in the main window

KNIME:

Copy the included node_pains_splitter.svl and node_pains_annotate.svl in your custom node directory eg. $MOE/custom/knime and restart KNIME. The PAINS Splitter node can be used to filter PAINS out of a data table. The PAINS Annotation node will give you flags for the rule hit by the molecules (firts hitted rule only!).
Remarks:
The rules were defined as SLN in the original publication. For the usage in MOE they have been translated into SMARTS and manually checked. Due to different aromaticity perception slight differences are expected. If you find any problem with this implementation please contact CCG Support.

Peptide and peptide-containing polymer structures are often stored in SMILES strings or other formats which do not retain information about amino-acid monomer identities and their relative positions in the polymeric structure. The "split_into_peptide_residues.svl" application attempts to regenerate the monomeric residues and their polymeric arrangement given a molecular structure where this information is absent. The amino acids are detected and named using the currently loaded rotamer library in MOE, and additional optional libraries of non-amino acid functionalities can also be provided for detection and naming of residue units.

This script will detect and encode the sequence information.

Last Update: February 2020

Calculates CNS BBB Score from "The Blood-Brain Barrier (BBB) Score"
Gupta, M, Lee, H J, Barden, C J, and Weaver, D F.,
J. Med. Chem., DOI: 10.1021/acs.jmedchem.9b01220 (2019)

1. Save this file in moefiles/svl
2. Use DBV | Compute | Descriptors | Calculate

This script is meant to automate the generation of multiple batch docking jobs where a single ligand is being docked into a database of proteins. This script outputs multiple .svl files that can then be docked on another computer. This script requires user input during the generation of the batch files but greatly automates the process.

If you run this SVL in MOE, it will ask you to select a MOE-Grid file to create a molecular surface. You can adjust the surface isolevel with a slider. Then you can use another grid file for the coloring of the surface. A typical application is creating a surface based on an electron density grid and using an electrostatic potential grid for the coloring.

Locating 'hot-spot' for binding to a protein structure is the goal of many computational chemistry applications.
Hot-spot analysis can help validate the importance of known interactions in holo complexes and suggest regions to target when designing ligands for apo structures. In the context of formulating biological-based medicines such as antibodies, hot-spots can be used to determine locations where excipients and other formulation molecules may aggregate on the protein surface and affect bulk properties of the medicine.

Open the SVL file in the MOE text editor then save and load the script from the text editor window.
Enter this command:
ligandRMSD ['reference structure', 'database', 'mol', 'receptor']
where "reference structure" is a .pdb or .moe file containing a crystallized reference structure with ligand bound and "database" is a .mdb file containing the results of a MOE docking run.
Ligand RMSD values will be written to a new column in the docking results database.

Description:
The five Methods in this PROTAC modeling package are designed to enable the in silico screening, comparison, and potentially the design of PROTAC molecules. Specifically, these tools have been designed to tackle two fundamental challenges in PROTAC modeling:

1. The prediction of the conformational ensemble accessible for a constructed ternary complex. The makeup of this conformational ensemble depends not only on the conformational flexibility of the PROTAC (primarily imparted via its linker), but also on direct protein-protein interactions. Conceptually, this challenge is similar to pose generation in traditional small molecule docking.
2. The identification of ternary complex poses generated by the conformational sampling above that are deemed biologically relevant, i.e., that resemble known ternary complexes. Conceptually, this challenge is similar to pose scoring in traditional small molecule docking.

A third fundamental challenge — the fruitful interaction of each binding end of the PROTAC with its respective protein — is readily addressable in MOE using standard Structure-Based Drug Design techniques. This challenge is not addressed by these tools, i.e., it is assumed that the binder moieties are active against their respective targets.
Please note: This download contains only the instruction manual and worked tutorial examples. To download the scripts and other necessary files, please email mdrummond@chemcomp.com

Description:
This script contains fragmentation and visualization GUI for ABINIT-MP Fragment Molecular Orbital method program in MOE.

Installation:

• Extract zip archive into your favorite directory ($FMOE_DIR)
• write the code running loader.svl to $HOME/moefiles/start.svl (local) or $MOE/custom/start.svl (global)

local function main []
 run '$FMOE_DIR/loader.svl';
endfunction

Source code repository:
https://github.com/drugdesign/FMOe

Description:
This script accepts HDX data from HDExaminer and DynamX software and maps it onto protein structures into MOE. Multiple protein chains can be colored with different HDX input files, so that changes in the HDX signal during complexation can be visualized for multimeric proteins. An automated alignment procedure ensures that differences in the numbering schemes between the HDX data and the residue numbers within MOE can be handled. Information for the individual peptides in included, as is the ability to color by difference or in side-by-side mode. Finally, the poses in a protein-protein docking database can be colored and scored by how well they correlate with the experimental HDX results.

Description:
This package contains more than 150 MOE nodes for KNIME written in SVL, the documentation and an installation video.
To use the nodes the MOE Extensions for KNIME have to be installed in KNIME either from the KNIME Partner Extension Update Site or from the MOE Extensions for KNIME (Offline Installation) available for download from the SVL Exchange.

Installation:

• Extract the zip archive from this download into $MOE/custom.
• Make sure the path to $MOE is set in the KNIME Preferences (KNIME|File|Preferences|Chemical Computing Group)
• Restart KNIME

These nodes have been tested with MOE 2019 but are also backward compatible to MOE 2018.

Example Workflows are available from the KNIME Hub.

For any questions contact support@chemcomp.com.

Description:
This package contains more than 150 MOE nodes for KNIME written in SVL, the documentation and an installation video.
To use the nodes the MOE Extensions for KNIME have to be installed in KNIME either from the KNIME Partner Extension Update Site or from the MOE Extensions for KNIME (Offline Installation) available for download from the SVL Exchange.

Installation:

• Extract the zip archive from this download into $MOE/custom.
• Make sure the path to $MOE is set in the KNIME Preferences (KNIME|File|Preferences|Chemical Computing Group)
• Restart KNIME

These nodes have been tested with MOE 2018 but are also backward compatible to MOE 2016.

Example Workflows are available from the KNIME Hub.

For any questions contact support@chemcomp.com.

Description:
This package contains the zipped Update Site for the MOE Extensions for KNIME Version 2.4.2, the documentation and an installation video.
It is meant for installation on computers that cannot access the KNIME Partner Extension Update Site.
Before installation make sure the KNIME Base Chemistry Types&Nodes Extension is installed.

Installation:

• Start KNIME
• Click on Help|Install New Software.
• Click the "Add" button.
• Click on "archive" and browse to the zip file from this download.
• Click OK and follow the instructions on the screen.
• Restart KNIME (MOE category in the Node Repository will be shown)
• Set the path to MOE in the KNIME Perferences (KNIME|File|Preferences|Chemical Computing Group)

More nodes are available as separate download from the SVL Exchange.

Example Workflows are available from the KNIME Hub.

For any questions contact support@chemcomp.com.

Description:
This package contains the documentation for the MOE Extensions for KNIME as pdf file. It also includes three videos showing the different steps for on-, offline installation and activation of SVL based nodes.

The MOE Extensions for KNIME and the SVL based nodes are NOT included and have to be downloaded separately.

The workflow examples from the tutorial are available from the KNIME Hub.

For any questions contact support@chemcomp.com.

Introduction:
The Bio-MOE package contains a series of custom SVL applications that are used for biologics modeling. These applications, accessible from the Sequence Editor and the Database Viewer menus, include tools for antibody humanization and back mutation prioritization, liability detection, patch visualization, and PTM motif annotation, aid importing and exporting collections of sequences, as well as visualizing protein properties and coloring antibody motifs.

This download contains only the manual.

Installation:

• Download and installation instructions are page 3 of the document

QFSS, Quick Federated Structure Search, provides means to search huge database stored across multiple mdb files. Based on our original substructure fingerprint and indexing technology, you can handle huge molecular database containing even hundreds of millions of entries stored separately in multiple mdb files. QFSS does not require any database engine, such as Oracle or SQL Server, nor high spec hardware. Even a laptop PC, you can search tens of millions entries in your mdb files.

QFSS SOAP version is also available. Via SOAP functions, you can run QFSS on remote server from your MOE client or from internet browser.

QFSS package contains SVL code, SSFP definition file and sample databases(ChEMBL and Leadlike Conformer Database) with index files for your convenience. The attached file contains only the manual. To receive the full download instructions for this package please contact ccg@molsis.co.jp

“BOILED-Egg” model by Daina et. al (2016) for predicting gastrointestinal absorption and brain penetration was adapted for MOEsaic. The y-axis (WlogP) of the model was substituted for the SlogP descriptor. The Model was re-optimized according to the dataset's SlogP vs TPSA.

A Monte Carlo (MC) method was used for computing the model, similar to the methodology described by Daina et al. To obtain optimal ellipses, 10⁹ independent MC runs were performed. True and false positives (TP, FP) as well as true and false negatives (FP, FN) were used for the binning statistics of each ellipse and evaluated using Matthews correlation coefficient (MCC ranges from -1, total disagreement between observation and prediction, to 1, perfect agreement).

The model achieve an MCC of 0.682 and 0.791 for the human instestinal absorption (HIA) and blood-brain-barrier (BBB) datasets, respectively. The original model reported similar MCC (0.653 and 0.747 for the HIA and BBB datasets, respectively).

How to include the model in MOEsaic:
This adapted model is available via this SVL function provided in the package. Copy the SVL file to the custom/svl directory of the moe installation. Restart the MOE/web server. The model should now be available in any MOEsaic session.

This script enables SVL highlighting in a katepart texteditor such as KWrite/Kate/KDevelop (part of KDE).
SVL files are recognized using their file extension .svl.

Kate is available in your repos (Linux) or here: https://kate-editor.org/ (Linux, Windows, Mac).


Installation (version numbers might need to be adapted to your KDE/katepart version):

Linux:
Copy the extracted xml file to

/usr/share/kde4/apps/katepart/syntax (all users installation; needs root privileges)

or

$HOME/.kde4/share/apps/katepart/syntax (single user installation)


Windows:
$HOME\AppData\Local\katepart5\syntax (single user installation)


*** Update ***
function list: MOE2019.01

Warning: This version does NOT work on moe-type database fields.

These protein descriptors are calculated on sequence only. I.e. there is no need for atomic data.

INSTRUCTIONS:
- copy this file to $HOME/moefiles/svl/
($HOME is your home folder)
- restart MOE
- the new descriptors will be available from:
DBV | Compute | Descriptors | Calculate...

List of descriptors:
seq_nChains: Number of chains, regardless of type.
seq_nChainsAA: Number of chains that contain at least one amino acid residue. (protein or peptidic ligand)
seq_nRes: Total number of residues.
seq_nResAA: Total number of amino acid residues.
seq_n***: Overall number of amino acid residues of type ***. *** = amino acid 3-letter code.
seq_MW: Sequence based protein mass in kilodalton. I.e. sum of internal amino acids masses + free termini. The masses correspond to the mass of the most abundant protomeric and isomeric species at pH = 7.0. I.e. Glu,Asp are deprotonated, Lys,Arg are protonated, His and all other amino acids are neutral. One pair of free termini is added per amino acid chain.
seq_n+: Number of positively charged amino acids. I.e. seq_n+ = seq_nLys + seq_nArg
seq_n-: Number of negatively charged amino acids. I.e. seq_n- = seq_nAsp + seq_nGlu
seq_PI: Sequence based isoelectric point as calculated in MOE release code ($MOE/svl/biopolymer/proprop.svl).
seq_aliphIdx: Aliphatic index. Ref: A. Ikai (1980) J Biochem 88: 1895-8
seq_GRAVY: Grand Average of Hydrophathy score. Ref: J. Kyte, R.F. Doolittle (1982) J Mol Biol 157: 105-32
seq_polariz: Sum of residue polarizabilities based on M. Swart (2004) J Comp Meth Sci Eng 4: 419-25. Please note that polarizabity varies with conformation!
seq_InstabIdx: Instability Index. Ref: K. Guruprasad, B.V.B. Reddy, M.W. Pandit (1990) Prot Eng 4: 155-61
seq_PI_Bellqvist: Sequence based isoelectric point as published by B. Bjellqvist et al. (1993) Electrophoresis 14: 1023-31
seq_ExtCoeff_280: Extinction Coefficient at 280 nm assuming that half of the cysteins form cystines.
seq_Absorb_280: Absorbance at 280 nm (0.1%).

--------------------------------------------------------------------
Description:
--------------------------------------------------------------------
Code to compute the PBSA estimate of the solvation free energy
of selected atoms in the MOE window (or on a database of molecules)

--------------------------------------------------------------------
Usage:
--------------------------------------------------------------------
MOE:
(1) Load the structure of interest
(2) Run this code. A GUI will appear.
(3) Set the 'Source:' option to MOE (atoms in the MOE window)
(4) Use the option pull-down to specify the atoms whose solvation
energy you wish to measure
(5) Press 'Compute'. The PBSA solvation energy will be written the
text field in the panel. e.g,

PB Solvation Energy: -8.00656 (kcal/mol)

MDB:
(1) Run this code. A GUI will appear
(2) Set the 'Source:' to MDB
(3) Use the 'VBrowse..' button to select the database
(3) Press 'Compute'. The panel will disappear and the PBSA solvation
energy will be computed for all entries in the database and
written to the field 'PBSA_sol'


TITRATE (i.e.., loop over salts and conditions)
(1) Load this function
(2) Enter

svl>PBSA_Solvation_Energy_MDB_titrate opts

at the SVL command line. Here opts is a tagged vector of options:

opts= [
dbfile: 'database filename',
mfield:'molecule field'
(defaults to first mol field in database)
salts: vector of salt tokens
(default: SALTS_TO_USE = ['NaCl','KCl','MgCl2'];
concs: vector of concentrations
(default: [0.1, 0.05, 0.01];)
temps: vector of temps
(default: 300)
]
(3) The code will loop over all combinations of salts, conc.
and temperature and output the PBSA solvation energies to
the database. The conditions are recorded in the database field
names
--------------------------------------------------------------------

The Excel interface is based on the Apache POI library which has to be downloaded and installed separately. We also recommend having Oracle Java installed.

Installation:
1) Download the Apache POI Library from https://poi.apache.org
2) Extract the downloaded archive file to your hard drive.
3) Add the poi library to your CLASSPATH environment variable.
for example if you extracted the library in c:/javalibs/poi-4.0.1
you would have to add:
C:/javalibs/poi-4.0.1/*;c:/javalibs/poi-4.0.1/ooxml-lib/*;c:/javalibs/poi-4.0.1/lib/*
to CLASSPATH.

4) Extract the archive file you downloaded from the SVL Exchange website.

5) Copy io_excel.svl into $MOE/custom/svl or $HOME/moefiles/svl

6) Copy SVL_Excel.class into $MOE/java

7) In MOE|File|Open you should get the option to Import the XLS file when selecting an excel file.

8) To export an excel file use:
db_ExportXLS[
'$MOE/sample/mol/gmp_inhibitors.mdb',
'$HOME/gmp_inhibitors.xlsx',
'Sheetname',
0,
0
];
where the first number after the sheetname is a flag indicating that only selected entries should be exported.
The second number is for selected fields only.

sdconf and sdconf.bat are tools that can generate small molecule conformations from the operating system command line. They can be combined with other SD files, as explained in the document "SD Pipelining Tools - How to Create a Fragment Database" that can be accessed from (MOE | Help | Tutorials)

Download and unpack the zip file. This contains the SVL script, a PDF document and example files.

(1) Select the protein_pocket_volumes.svl file from a MOE| File | Open panel.

(2) A GUI will appear.

(3) Use the 'Pocket Set A' and 'Pocket Set B' options to select the 'A set' and the 'B set' for the operation. The pocket will be defined around the ligand atoms. To compare proteins with no ligands, first use the Site Finder to create dummy atoms in the binding site. These will then be used to identify the pockets.

(4) Press Create. The pocket volume comparison will proceed and a new graphics object rendering the volume difference will be create and listed in the panel. The volume enclosed by each surface is shown in name of the graphic object.

To install extract the zip file into your MOE installation folder. After restarting MOE you will find a new entry in the configuration (gear wheel) menu (MolPort). Enter the username and password you got from MolPort there. Click OK. Now you can enter a MolPort ID in the text widget of the builder (e.g. MolPort-000-219-650) to get the structure or use the Interface MOE|Extra|MolPort to search or import structures from MolPort. Databases for imported structures contain links to the MolPort website which can be opened from the database browser by clicking (Open in Browser).
If you have questions, please let us know.

To get good QSAR model in MOE, you have to run several commands under the 'Compute' pull-down menu of database viewer, such as Partial Charge, Descriptors, QuaSAR-Contingency, QuaSAR-Model and so on. After you get one model, you have to delete less important descriptors step by step to get better model.

AutoQuaSAR was designed to make this process as easy as possible. By single click, you can get good QSAR model directly either from sd or mdb file with molecule and activity data only. AutoQuaSAR automatically makes working file, calculate descriptors and generate QSAR model by removing less important descriptors step by step to find out the model with best XR2 (XA1, in case of binary model).

Default setting works fine for most cases but AutoQuaSAR also provides lots of options to refine model.

Usage:
load autoqsar.svl
type AutoQSAR [] in the commad line
or just double click autoqsar.svl in 'Open' panel.
Click triangle in 'Oritginal File:' column and select mdb or sd file name in the list of files. If you want to use a file in other directory, select 'Others...' to open File Selection Panel.
For example, select $MOE/sample/mol/blood_brain.mdb for 'Original File' and click 'Fit'.
You will get the result just like this screenshot.

New in ver.2015.12.15
- general bugfixes

New Features in ver.2006.08.02
- Support QuaSAR-Evolution(ver.2006.06.12)
- Added structure display area
- keep all models to allow users to browse/select models

ver.2008.05.14
revised code with manual

Based on a given alignment of a set of small molecules with activities, a grid potential analysis is applied using the principles of CoMFA. Hereby, steric and electrostatic interaction energies of a positively charged carbon probe atom placed on grid points around the aligned molecules are determined and grid potentials are calculated. Then a partial least square model is created by fitting the activities to the grid point potential values. This process starts with a principal component analysis of the descriptors to reduce the number of components in the model to a value set by the user. The model coefficients for the principal components are then mapped back to the descriptors. These contributions can be visualized to identify areas with positive or negative contribution to the model spatially (pre-condition: numerically higher activity values mean higher biological activities).

In 2013 the group of Prof. Rarey (Uni Hamburg) in cooperation with the CADD group of Roche AG published the Torsion Analyzer.
This program accesses the quality (i.e. the strain) of individual bonds of the conformation of a molecule.
The bonds are colored in green/yellow/red based on their individual quality.


INSTALL INSTRUCTIONS:
In the following instructions:
- $HOME needs to be replaced with your home folder (e.g. C:\users\myaccount\)
- folders that do not already exist need to be created.

0) Make sure to use MOE 2018.0101 or newer. Check the Window Title of the MOE window.

1) Get a copy of the XML file which contains the relevant data from the website of the Rarey workgroup:
https://www.zbh.uni-hamburg.de/forschung/amd/datasets/torsion-analyzer-datasets.html
Rename the XML file to rochetor.xml and copy it to $HOME/moefiles/lib/

2) Extract the attached rochetor.zip and copy
rochetor.svl to $HOME/moefiles/svl/
rochetor.menu to $HOME/moefiles/menu/

3) Restart MOE. There will be a red button button "TorAnalyzer" in the RHS. This button will allow to toggle the Torsion Analyzer.


REFERENCES:
Torsion Angle Preferences in Druglike Chemical Space: A Comprehensive Guide
Schaerfer,C., Schulz-Gasch,T., Ehrlich,H.C., Guba,W., Rarey,M., Stahl,M.
J. Med. Chem. 56 (2013) 2016-2028 http://dx.doi.org/10.1021/jm3016816

Torsion Library Reloaded: A New Version of Expert-Derived SMARTS Rules
for Assessing Conformations of Small Molecules.
Guba, W., Meyder, A., Rarey, M., and Hert, J.
J. Chem. Info. Model. (2015) DOI: acs.jcim.5b00522

The authors of the of the original Torsion Analyzer are not affiliated with the authors of this implementation! Please report all bugs to
support@chemcomp.com!

*** update ***
Terminal rotors, such as terminal methy groups, are now filtered. WARNING: This will reduce the number of tosion angles reported by the QuaSAR descriptors!

This script provides a few basic tests to measure the performance of hardware/software running MOE.

Extract the archive to $HOME/moefiles (for current user) or $MOE/custom (for all users).

DESCRIPTION:
This package contains more than 100 small molecule descriptors.
Most of these molecular descriptors are suitable for statistical modelling of molecular properties or activities (e.g. machine learning).

All descriptors use the descriptor plug-in system in MOE. The are divided 3 classes:
2D - conformation-independent
i3D - conformation-dependent
x3D - dependent on conformation and the structure currently in your MOE window

Note: each of the descriptor files contains a number of related descriptors. There are redundant descriptors in case they make sense in more than one file.



INSTALLATION:
1) unzip attached folder
2) Copy the desired descriptor file(s) to $HOME/moefiles/svl, where $HOME is your home directory
(e.g., Windows: C:\Users\walter\moefiles\svl\ or Unix: ~/moefiles/svl/)
3) Restart MOE to load the new file

Notes:
- The directory must be created if it does not already exist.
- $HOME can be determined at the SVL Commands windows. Type:
svl> HOME


USAGE:
The new descriptors are available from the Calculate Descriptors panel
(DBV | Compute | Descriptors | Calculate).

Rogers et al., (ECFP)
J. Chem. Inf. Model. 2010, 50, 742-154;

Differences to the Pipeline Pilot implementation may arise from:
- differences in treating charge-seperation of some moieties (Nitro etc.)
- hypervalent atom treatment (e.g. Sulphur and Phosphorous oxides)
- a different aromaticity perception

update 10-Oct-2018 (ma) E/FCFP now available as MOE descriptors

!!!
NOTE: This latest version (last change: 04-feb-2016 (ma))of the script generates fingerprints that are incompatible to earlier versions
!!!

This is a simple XYZ parser which uses the plug-in system of the MOE Open and Save panels.

Please note that there are multiple XYZ format descriptions. The one used here is:

5

C 0.0000 0.0000 0.0000
H 1.0900 0.0000 0.0000
H -0.3633 0.7267 0.7267
H -0.3633 -0.9926 0.2660
H -0.3633 0.2660 -0.9926

This is mostly an example of how to use the plug-in systems in MOE.

INSTALLATION:
- Copy the file to $HOME/moefiles/svl/ ($HOME is your home folder).
- Restart MOE.
- The XYZ file format should now be available in the Open and Save panel.

Simply unzip the contents to your sublime folder Packages/User/. After restarting Sublime 3 you will have at the bottom right a new language from the drop-down list:

Scientific Vector Language

Sublime should also be able to autodetect SVL files by their extension .svl.

Automatic flexible alignment for all molecules in a database against reference
molecules from the main MOE window. This function uses the FlexAlign []
function. If an output database filename is not supplied then it
will be based on the input filename. For instance if the input database
was in.mdb, the output database would be in_align_all.mdb.

If you set the molecule source to "MOE" then it runs as in the normal
Flexible Alignment panel in MOE.

USAGE

Save this file to your hard disk. For it to be available for use by all
of the users of MOE, create a directory called $MOE/lib/svl/custom/run
and save this file to that directory. eg on Windows and Linux,
this might be: c:\moe\lib\svl\custom\run\align_all.svl
/usr/local/moe/lib/svl/custom/run/align_all.svl

For the script to be available for use by a single user, then it
is best to create a directory in your home directory called svl,
and create a directory called run in that directory. Save this
file into that directory. eg on Windows and Linux, this might be:

c:\Documents and Settings\user_name\svl\run\align_all.svl
/usr/people/user_name/svl/run/align_all.svl

(a) From the MOE | File | Open panel, select this file
and click the "Run" button.

(b) Alternatively, at the SVL command line, type:
run 'align_all.svl'

The underlying SVL function is align_all. To use this, "load" this
file and at the SVL command prompt, type a command like:

svl> align_all []
svl> align_all ['in.mdb', 'out.mdb', options]
svl> align_all ['in.mdb', ', [heavy_rms:0.1, nfail:40]] //refinment only
svl> align_all ['in.mdb', 'out.mdb', [file:']]

- If you use MOE batch, then you have to pass the name of the input MDB file or
the script will fail. eg
moebatch -exec "run['align_all.svl', ['in.mdb',',[nfail:10]]]"

Align_all options are the following:

esel: 0, // use selected entries only
nsolns: 0, // no. of output solutions
mfield: ', // molecule field
dump: 1, // suppress CLI output if 0
fixed: 1, // fix template coordinates
write_template: 1 // 1: write only aligned pose to output (Pose)
// 2: write both in one field (Alignment)
// 3: write pose and template to diff fields (both)

To run the script with the "refine only" option, pass an option with file:'.
For other flexible alignment options see FLEX_DEFAULTS below.


*References*
An Alternative Method for the Alignment of Molecular
Structures: Maximizing Electrostatic and Steric Overlap
Tetrahedron Computer Methodology, Vol. 3, No. 6C, pp615-633, 1990

A New Approach to Probing Conformational Space with
Molecular Mechanics: Random Incremental Pulse Search
J. Am. Chem. Soc. 1989, 111, 4371-4378
D.M.Ferguson D.J.Raber

Calculates CNS MPO from "Moving beyond Rules:
The Development of a Central Nervous System Multiparameter Optimization
(CNS MPO) Approach To Enable Alignment of Druglike Properties;
Wager TT, Hou X, Verhoest PR, Villalobos A.
ACS Chem Neurosci. 2010 Jun 16; 1(6): 435-449"

1. Save this file in $HOME/moefiles/svl
2. Use DBV | Compute | Descriptors | Calculate

The ZINC database is an excellent resource for structural information on small molecules. The script let's you import pre-calculated 3D conformations including relative energies from the DB2 format to a MOE Database. ZINC 2D structures can be downloaded as SDF file and imported using standard MOE tools (please contact support@chemcomp.com if you need further assistance).

Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs).

This script builds models for the linker between two ligands, bound in two binding sites.

There is more information in the Open Access paper:

Pérez-Benito, Laura, et al. "The Size Matters? A Computational Tool to Design Bivalent Ligands." Bioinformatics.
https://doi.org/10.1093/bioinformatics/bty422

Additional information is in the supplementary information for the pap

Usage:
(1) Load this function
(2) At the SVL command line, enter

svl> db_moledit []

(3) Select target MDB file at the panel.

For more information, please see the bubblehelps.
This program may be useful if you want to extract ligand, receptor, or contacting solvent from molecule type field in MDB of molecular dynamics output.

Run PSILO searches from MOE

New in this version:
- Pocket Similarity Search
- Highlight hit atoms and meters for 3D interaction search
- Bugfixes

The following searches are currently implemented:
- Full Text Search
- Ligand Search (via SMILES /MOL)
- Small Molecule Name
- Protein Sequence Search
- Protein Name
- Title (PDB)
- Project Code
- PSILO Family
- Similar Pockets
- CIF Field
- 3D Search

Workflow for Exercise 1 "Compound Library Management" as described in the documentation of the MOE Extensions for KNIME version 2.4.0
See documentation for further details.

Description:
Application to detect duplicate conformations in a MOE database.
Use cabn set RMSD theshold for unique conformations. Option of
eithe selecting of deleting duplicate conformations

Note: ALL entries must not be empty!

Usage:
(A) GUI
If you run this code from a MOE File Open panel(click on "Run SVL"),
a GUI will appear. Use the Browse button to select the input
database. Specify the molecule and (optional) mseq fields
(in the case of databases with multiple conformations).
Molecular symmetry will be taken into account.

Choose the actions for the duplicates - either select in the
database viewer or delete. Press OK to perform teh operation.

(B) Command Line:
The operation can be run using the SVL command:

svl> conf_unique opts

where opts is a tagged vector of options:

opts = [
dbfile:'database_filename',
mfield:'molecule field name'
mseq:'mseq field'
RMSD:'rmsd threshold in A'
action:'Select' or 'Delete' -> action to perform on duplicate

Installation Instructions Updated for MOE 2016

Please refer to the provided manual for details of usage.

To stay up to date about future updates, please subscribe to our newsletter: http://eepurl.com/_PqlD

How to cite Docktite:
Scholz, C.; Knorr, S.; Hamacher, K.; Schmidt, B., DOCKTITE — A Highly Versatile Step-by-Step Workflow for Covalent Docking and Virtual Screening in the Molecular Operating Environment. J. Chem. Inf. Model. 2015, 55, 398-406

This is a workflow for creating QSAR models in MOE. The process works like so:

1. Load (and optionally back up) an MDB of molecules with an activity field
2. Calculate descriptors and specify the name of the activity field
3. Split the database (optional step) into training and test sets. Various methods are available, including k means clustering and selection from binned activity values.
4. Select descriptors. The selection can be pruned by looking for low variance or highly correlated descriptors. In addition the training and test sets can be compared in PC space based on the current descriptor selection.
5. Build the QSAR model. Currently only the PLS and PCR methods are available.
6. Validate the model; a variety of methods, including MOE's standard Leave One Out method, are available.
7. If a test set has been created, the model can now be applied to the test set.

In the zip file you download you will find a PDF explaining how to use the Wizard and the SVL file itself. To start a QSARWizard session, simply save the SVL to your hard-drive and run it e.g. using the MOE | File | Open panel, browse to qsarwizard.svl and click OK.

RECON Descriptors

RECON is an algorithm for the rapid reconstruction of molecular charge
densities and charge density-based electronic properties of molecules,
using atomic charge density fragments precomputed from ab initio
wavefunctions. The method is based on Bader's quantum theory of Atoms
in Molecules.

A library of atomic charge density fragments has been built in a form
that allows for the rapid retrieval of the fragments and molecular
assembly. The RECON algorithm reads in a molecular database, determines
atom types and environments, assigns the closest match from the library
of atom types and combines the densities of the atomic fragments to
compute a large set of new and traditional QSAR descriptors. QSPR and
QSER indices for large pharmaceutical databases or proteins can be
computed within seconds.

This is an SVL implementation of MOE RECON TRAD. A much faster linux
executable with additional features is available by license through the
RPI office of Technology Commercialization http://www.rpitechnology.com

Relevant references:

S. Oloff, S. Zhang, N. Sukumar, C.M. Breneman*, and A. Tropsha*,
"Chemometric Analysis of Ligand Receptor Complementarity: Identifying
Complementary Ligands based on Receptor Information (CoLiBRI)",
Journal of Chemical Information and Modeling, 46 (2), 844-851, 2006.

B.K. Lavine, C.E. Davidson, W. Katt, C.M. Breneman*, "Analysis of
odor structure relationships using electronic Van der Waals surface
property descriptors and genetic algorithms," ACS Symposium Series
(2005), 894 Chemometrics and Cheminformatics, 127-143.

B.K. Lavine, C.E. Davidson, C.M. Breneman*, W. Katt, "Electronic van
der Waals surface property descriptors and genetic algorithms for
developing structure-activity correlations in olfactory databases"
Journal of Chemical Information and Computer Sciences, 43 (6):
1890-1905, Nov-Dec 2003.

C.M. Breneman* and M. Rhem, "A QSPR Analysis of HPLC Column Capacity
Factors for a set of High-Energy Materials Using Electronic Van der
Waals Surface Property Descriptors Computed by the Transferable Atom
Equivalent Method." J. Comp. Chem. 18(2),182-197, 1997.

http://www.drugmining.com/files/RECON/recondoc/WinRecon.html

Installation of RECON TAE Molecule Descriptors

Create the directory $MOE/custom/svl/ (if necessary)
Copy the q_recon_tae.svl file to $MOE/custom/svl/

For instance, if MOE is installed in c:\moe on a Windows PC, this would be
c:\moe\custom\svl\q_recon_tae.svl

If MOE is installed in /usr/local/lib/moe on a Unix system, this would be
/usr/local/lib/moe/custom/svl/q_recon_tae.svl

Copy entire 'reconBin' directory to a nice location (e.g. $MOE)

Set system variable 'recon_moe_root' to the reconBin directory (remember
to include the final '\' or '/' in your variable!) e.g.

recon_moe_root = c:\moe\reconBin\


An alternative to setting a system environment variable is to
replace line 2295 in q_recon_tae.svl with with a line that
points to the reconBin directory. Examples of this could be

local RecRoot = 'c:/moe/reconBin/';
or
local RecRoot = '/usr/local/lib/moe/reconBin/';

Note that paths within MOE use / forward slashes, on both Windows and Unix.

Calculates Central Nervous System Multiparameter Optimized Index from

"Moving beyond rules: the development of a central nervous system multiparameter optimization (CNS MPO) approach to enable alignment of druglike properties."
Wager TT, Hou X, Verhoest PR, Villalobos A.
ACS Chem Neurosci. 2010 Jun 16;1(6):435-49. doi: 10.1021/cn100008c. Epub 2010 Mar 25.

This scripts sends protein sequences to PFAM website and pareses the resulting hits, colors various detected domains and creates MOE 'sets' to access the domains as detected.

This is a basic interface to the ORCA QM package written by F. Neese. It allows to read and write ORCA coordinates (.xyz) and facilitates the creation of input files for reaction coordinate scans and constrained calculations.

USAGE:
MOE | File | Open -> ORCA_coordinates.svl

If you want to start ORCA from the panel the following executables and paths need to be specified in the script:
ORCA_DIR, QUEUE_TMP_DIR, MPI_EXE_DIR, MPI_LIB_DIR

The profiles allow you to create a series of input files with identical header and footer. E.g. calculation of NMR shieldings for different compounds with identical settings.

Press [Show] at the very bottom of the panel to get the input file w/o starting it locally or on a queing system.

For ORCA options etc. please refer to the ORCA documentation.

The output of ORCA calculations can be read by the script (geometric data only!). Scan calculations are imported to an MDB file and can be visualized using MOPAC_rb_energyTable.svl (uploaded to SVL exchange as well).

I do NOT guarantee that this script will work with future releases of ORCA and MOE.
Tested with MOE2015/2016, ORCA3.0.3 in a Linux environment and a SUN grid engine derivative.

UPDATE for MOE2016

All scripts have been tested with MOE2016.08 in a Linux environment.
MOE built-in public domain MOPAC does NOT work due to missing support of mixed cartesian and internal coordinates.

This collection consists of 5 svl scripts:
- MOPAC_rb_createInput.svl
- - create a MOPAC input file for 1D and 2D calculations
- MOPAC_rb_interpretOutput.svl & MOPAC_rb_lib.svl
- - extract structures and energies from MOPAC output file, import them into an .mdb
- MOPAC_rb_energyTable.svl
- - use the extracted structures and energies to visualize an "energy hypersurface"
- MOPAC_rb_applyFixed.svl
- - restore fix-flags when using a MOPAC output geometry as input for a new MOPAC calculation

These scripts are stored in a .zip file.


=== MOPAC_rb_createInput.svl ===
FEATURES:
- creates a MOPAC input file
- fixed atoms get no-optimization flag
- "normal" atoms will keep their cartesian form
- reaction coordinate atom(s) will be transformed into internal coordinates and appended to the "normal" atoms
- for restoration of atom-names and residues a .mmeta-file will be generated which can be used in

BEFORE running this script:
- adjust the constants MOPAC_EXE, MOPAC_LIC and MOPAC_LIB

USAGE:
- prepare the system you want to use as input (all atoms will be used - that is: every atom in Atoms [])
- run this svl-script via 'File' - 'Open'

=== MOPAC_rb_lib ===
- library functions for MOPAC_rb_interpretOutput.sv

=== MOPAC_rb_interpretOutput.svl ===
FEATURES:
- imports structures and energies of MOPAC output file
- if present, uses a .mmeta file to restore atom names and residues

BEFORE running this script:
- adjust the constant MOPAC_RB_LIB to point the location of the library file

USAGE:
- run this svl-script via 'File' - 'Open'


=== MOPAC_rb_energyTable.svl ===
FEATURES:
- visualizes an "energy hypersurface" of a MOPAC grid/scan calculation

USAGE:
- run this svl-script via 'File' - 'Open'
- select the .mdb file (needs to be same format as output of MOPAC_rb_interpretOutput.svl)


MOPAC_rb_applyFixed.svl
In MOE molecular databases (.mdb) are no information whether an atom is fixed or not (fix flag). For every no-optimization flagged atom present in the MOPAC input .dat-file this script looks for an atom of the same element at the same position in the MOE system, and fixes this atom.

FEATURES:
- restores the fix flags of an MOPAC output scene loaded from an .mdb

USAGE:
- load the MOPAC scene from the .mdb
- run this svl-script via 'File' - 'Open'
- select the .dat-file (MOPAC inputfile which has been used to generate the scene)
- PLEASE, check the results: 'Select' - 'Potential' - 'Fixed

INSTALLATION

1) Copy this file to $HOME/moefiles/svl/run/, where $HOME is the user’s home directory (e.g., Windows: C:\Users\walter\moefiles\svl\run\ or Linux/Unix/macOS: ~/moefiles/svl/run/)

2) Start the script using, from the SVL command line or a menu:
run 'plane_gui.svl'

Notes:
- The directory must be created if it does not already exist.
- Alternatively, drag-and-drop this SVL file in MOE window to run.
- $HOME can be determined at the SVL Commands windows. Type:
svl

Installation

1. Unpack the zip file in any directory.
2. Copy the two folders (menu and svl) to your $HOME/moefiles folder.
Now you should have $HOME/moefiles/svl/ted_extensions.svl and
$HOME/moefiles/menu/menu-ted_2016. Please be aware that the menu will over-ride
the default MOE text editor menu and any other modifications to those menu
you might have made.
3. Restart MOE.

Added Features

SVL Search (TED | Find | SVL Search or Ctrl + f in the TED window) - runs the SVL search utility. This allows a more fully featured search either in the current file or in MOE or User directories (as defined in step 4 of the installation process). Searches in file may fail if the file has been modified; please save then search if problems are encountered.

Manual Page (TED | SVL | Manual Page or Ctrl+m) - to enable this function, first type _UPDATE_ in a text editor window, highlight the text and choose TED | SVL | Manual Page. Once the index has been built (you can monitor the indexing in the SVL Commands Window), the manual page for an SVL function can be launched by highlighting the function name and choosing TED | SVL | Manual Page.

Default text blocks: Insert Header (Ctrl+h) and Database Looper Motif will insert a standard SVL file header (as seen in MOE SVL files) and a sample db_NextEntry-style database loop.

Block indent (Ctrl+. indents in, Ctrl+, indents out) - allows easy formatting of SVL. Select a few lines of text and indent it further right or left to make your code easier to read.

Block commenting (Ctrl+I comments out, Ctrl+Shift+I comments back in) -
allows easy commenting in/out of multiple lines of code in a similar manner to block indenting.

Move cursor back one word with Ctrl+b.

Unload the currently viewed SVL file from MOE by TED | SVL | Unload

TED | File | Save, TED | SVL | Save and Load, and TED | SVL | Save and Run all perform a "two stage" save, i.e. the last saved version is copied to a file with the same filename apart from an additional ~ suffix, then the main file is saved.

TED | Edit | Copy doesn't clear the clipboard if no text is selected

TED | Edit | Paste will delete any selected text before pasting at the cursor position.

TED | SVL | Save and Read Menu File will save the current file, reload the standard MOE menus then load the current file with the ReadMenuFile SVL command.

TED | SVL | Execute Command will execute any selected text, or if there is no selected text it will open a command line in the current text editor. This makes it fairly easy to step through a function executing one line at a time. Selected text is automatically expanded to complete lines. If the first line starts with "local", "global", "static", or "const" then that word is removed (this only affects the first line).

TED | SVL | Insert Date inserts today's date at the cursor position

TED | SVL | History Browser launches a listbox containing the contents of your history[]. Commands can be executed by selecting one or more of them and clicking "Execute", or can also be executed by double-clicking. Selected commands can be sent to a MOE text editor by clicking on "Make script".

This application uses a collection of proteins aligned using a common ligand and will highlight pharmacophore features that can be used to maximize
selectivity for each individual protein. An experimental scoring table (provided in the dabatase file imatrix.mdb) is used for this.

The resulting visuals can be used to modify the ligand to increase selectivity, i.e. to improve the activity on the chosen target while reducing the activity on others.

Sample files provided for quick testing.

This script will take as input the MDB from an MD simulation and write out every nth frame as a BMP file.
The user needs to set up the visualization before hand.

Load the protein in MOE and display ribbons. The "run" this script.
By default all the receptor chains with protein ribbons will be
modified. Chain keys can be supplied as an argument to adjust
only specific protein chains.

(a) From the MOE | File | Open panel, select this file
and click the "Run" button.

(b) Alternatively, at the SVL command line, type:
run 'ribbon_sidechains.svl'
run ['ribbon_sidechains.svl', protein_chains]

This implementation of the classic game tetris is published for demonstration purpose only and does NOT contain any scientific relevant methods etc.

For those of you, who need to relax some minutes ... enjoy

(Updated to work with file locking in MOE 2015.1001)

Syntax highlighting for the NEdit text editor for Motif/X11 available from https://sourceforge.net/projects/nedit/

To use this, download this file and save it to ~/.nedit/nedit.rc

Summary:
- Reads a MOE database or a folder containing MOE recognized structures
- Converts all files or database entries to the format give in 'Internal Format'
- Uses a simple mark-up syntax to generate command line required to call third party binary
- Uses the same format to designate the output files expected from third party binary
- collects the results in a MOE database.

It can also create batch scripts (unix only, for now) to do this offline.

The Scenes_Manager is a great way to enable collaboration between MOE users.
This will allow a user to save "scenes" in a fashion similar to other molecular viewers. Saved scences can be shared electronically via a .mse and accompanying .moe file.

It is common for small drug-like compounds and peptides that are constructed from SMILES strings and even SDfiles to have amide bonds where the hydrogen of the nitrogen atom and the oxygen of the carboxyl group result in a dihedral angle of approximately 0 degrees. This function sets all amide dihedral angles to 180 degrees and provides the option for energy minimization (this functionality is not implemented in the MDB version because of the energy minimization function with the MDB framework). The interactive (MOE 3D window) version allows the user to select the amide bond of interest.

The MDB version of this function provides the option to overwrite the original molecule or create a new column (field) for the compounds. The MDB version does indicate how many amide bonds were corrected. It is recommended that after changing the amide dihedral angle(s) within a database of compounds that the MDB energy minimization (MDB | Compute | Energy Minimize) function be performed using the desired molecular force field.

Note: Amide groups within rings are not modified.

The default values are:

No energy minimization for MOE 3D window (nrgOpt:0)
Column (field) with molecule is 'mol' (molField:'mol')
Save corrected compound to a new column 'modMol' (newMolField:'modMol')
Maximum peptide length is set with 'maxPeptideLength' (maxPeptideLength:15)

Instructions

Load this function
load '~/code/svl/utilities/db_SetAmide180.svl'

For a compound in the MOE 3D window (either a single small-organic drug-like compound or a selected molecule)
SetAmide180[]

For a compound in the MOE 3D window with energy minimization of the new conformation
SetAmide180[nrgOpt:1]

For a collection of compound in a MDB saving the new conformation to the originating column of the compound
db_SetAmide180['drugLike_compounds.mdb', [molField:'mol', newMolField:']]

This function counts the occurrence of the 166 MACCS keys for each molecule in a MOE database (mdb). The number of each key's occurrence is written to a user specified CSV file. The option to include the compounds' names and endpoint values in the CSV file is an option from the opts variable.

The default values are:

Column (field) with molecule is 'mol' (molField:'mol')
Include the endpoints if the endpoint column name is provided (endpointField:')
Include the compound's name or SMILES string. Not indicating a compound name file (nameField:') uses the compound's name returned with cName Chains[] or the compound's SMILES string if no name is present.
Column headers for MACCS keys will be returned in lower case letters (uppers:0)

Instructions

Load this function
load '~/code/svl/utilities/MACCSkeyCounts.svl';

Calculate the MACCS Key Counts for a MDB of compounds
MACCSkeyCounts['Selwood_dataset.mdb', 'Selwood_MACCSkeyCounts.csv', [molField:'mol', endpointField:'pEC50', nameField:'name']]

Display a MOE readable structure (*.moe, *.pdb, etc) via internet browser.

Requires jsmol. Please contact support@chemcomp.com if you have any difficulty installing jsmol.

residuesOfInterest = Residues[] |
cNumber rChain Residues[] == 2 or cNumber rChain Residues[] == 4;
atomEnvironment = Atoms[] |
cNumber aChain Atoms[] == 2 or cNumber aChain Atoms[] == 4 or
cNumber aChain Atoms[] == 5 or cNumber aChain Atoms[] == 6;
solvent_exposure_difference[
residuesOfInterest:residuesOfInterest,
atomEnvironment:atomEnvironment,
diffCutoff:0,
colorResidues:0,
printAtoms:1,
outfileSADiff:'solvent_exposure_difference.txt'];

DESCRIPTION:
The buttons contain: a molecular builder with added handling of chirality constraints, a simplified interface for flexible alignment of structures, a simplified interface for Contact Potentials and a pharmacophore query editor.

USAGE:
Save the attached files in your working directory (along with your molecule files) and open the moe-menus.pharmacophores file with the MOE File Open panel. A new set of buttons will appear at the left side of the MOE window.

D. Bondarev

Run Torsion Scan via internal MOPAC engine.


USAGE:

1. Load a desired small molecule structure in MOE.

2. Energy minimiza using MOPAC at AM1 level using
MOE | Compute | Energy Minimize | QM

3. Save this file to your working folder and run it by using
MOE | File | Open

then select 'torscan_mopac.svl' and press 'OK'

The code will run a '1scf' job via MOPAC and will display the zmatrix
created in MOPAC. Heavy atoms in main MOE window will be labelled
accordingly.

4. Select the desired dihedral from the list by clicking on the corresponding
line, then press 'Dihedral'.

ATTENTION: please verify that the selected atoms in main MOE window
correspond to a valid dihedral torsion. Currently we make no attempt
to verify this.


5. Press 'Run MOPAC'.

The script will use the internal MOPAC engine to run a torsion scan and
will plot the results.

The MOPAC input/output files can be displayed by checking 'Open Results
in Text Editor'.

These files are written into temp folder. To find the temp folder, Open
the SVL command window and type TMP. All files are named mopac*.*.

The current version uses external MOPAC such as MOPAC2012 for the actual
torsion scan, since some bugs in MOPAC7 affect certain structures.

Given a structure and a threshold in kcal/mol, the script first detects all "links". A link exists when two neighbouring residues have hydrophobic interaction stronger than the given thresholds.

The script then detects and isolates networks made solely via hydrophobic links as described. Users can include or exclude non-protein residues, multi chain interactions, etc.

An interactive histogram of all links in the system is also displayed, to help users set the hydrophobic threshold value.

This code computes an electrostatic similarity value between a molecule in the MOE window and (each of)a database of compounds previously aligned with the target molecule (for example, resulting from a FlexAlign procedure). It is observed that compounds with higher electrostatic similarity values may show similar biological activities. Thus this procedure may be used to rank docking poses where a known ligand is the target molecule, or can be used for ranking alternate scaffolds from a lead hopping program such as BREED. Electrostatic similarity is computed by solving the PB equation for both target and query compounds and computing a grid based Tanimoto (Tes) between the resulting grids. Thus results are in the range -1->1, with 1 = identical field and -1 = inverse field. This code uses the original PB solver in MOE prior to re-write in 2005.09 release and thus these dependent MOE functions are included in the script to work with versions of MOE released after 2004.

Usage :
(1) Load a molecule into the MOE window and compute MMFF charges. Also calculate MMFF charges for database ligands. nb all code has been tested using MMFF charges so these are recommended for correct use.
(2) Run db_ElectroSim[mdb,molfield] where the mdb contains a set of molecules in 'molfield' (defaults to mol if absent) to compute values for. Procedure outputs scores in field 'Tes' (range -1->1, where 1 = identical field)

Simply open the script if MOE: MOE | File | Open and you will get a graphical panel that allows you to pick the source, target, and number of entries per database.

Upon future execution, use the same output folder, it scans for already downloaded items and only downloads the new ones, and appends them to the last volume.

Syntax highlighting for the jEdit editor available from www.jedit.org

Save the xml file in JEDIT/modes and add the following lines at the end of JEDIT/modes/catalog

<MODE NAME="SVL" FILE="svl.xml"
FILE_NAME_GLOB="*.svl" /

Starting with a force-field energy minimized complex, the program calculates and displays the atomic interaction forces between the ligand and the pocket, as well as the per residue interaction energies of such interactions.

Originally published in Interaction force diagrams: new insight into ligand-receptor binding, H. Shadnia, J.S.Wright, J.M.Anderson Journal of Computer Aided Molecular Design, Vol. 23, No. 3, 1573,185-194, 2009.

Description
Use the EHT pharmacophore scheme to calculate where the Don2, Acc2, Aro and
Hyd features are, bin the distances along the triangle edges and convert to
fingerprint indexes. Hyd-Hyd-Hyd triangles are removed, and rings that are
both Hyd and Aro are only counted as Aro.

Usage
Load this file into MOE, e.g. save to $HOME/svl and restart MOE.
Use the fingerprint tools as normal, search for FP:EHTTAT.

DESCRIPTION:

db_nb_mols_incommon calculates a matrix for the number of compounds each database
has in common with each of the other databases based on either structure or
fingerprints.

db_compare_fp can be used to output the average tanimoto between two databases
based on fingerprints if called with

db_uniqBits_fp calculates a matrix for the number of unique fingerprint bits per
database, i.e.the number of bits not present in both compared databases

USAGE:

1. Save and load this function

2. A) At the SVL command line prompt type the following

svl> db_nb_mols_incommon [mdbs, outfile, opt]

where,

mdbs is a vector of tokens of filenames of the databases to
operate on
outfile is a token representing the filename of the output file
opt is a tagged vector of options
[
fp_flag: flag,
thresh: threshold,
fp_code: fingerprint code,
sim_code: similarity code,
average_Tani: flag,
write_cpds: flag,
cpds_outfile: cpds_outfile
]
- fp_flag: if 1 then similarity is calculated based on fingerprints
default is 0 (based on structure)
- molecules with similarity above the value of threshold are
regarded to be the same

- average_Tani: if 1 the the average similarity between two
databases is calculated instead of the number of molecules which
are regarded to be the same. fp_flag will be set to 1 automatically
(the average similarity is calculated in the following way: take the
maximum value for each row (so the highest similarity for each entry
to any other compound in the comparison database), then the average
value for those row-wise maxima.)

- fp_code and sim_code are the fingerprint and similarity metrics
to be used in the comparison. A list of fingperprints can be found
using the command ph4_FingerprintList []
Default options are fp_flag:0, thresh:0.85 using bit-packed MACCS
keys and the tanimoto coefficient.

- write_cpds: if 1 then another textfile is written containing all
compounds and the databases that contain the compound; if a cpd is
only found in one db, then it is not written to the textfile
- cpds_outfile: name of the output file for the above

2. B)
svl> db_uniqBits_fp [mdbs, outfile, opt]

mdbs is a vector of tokens of filenames of the databases to
operate on
outfile is a token representing the filename of the output file
opt is a tagged vector of options, in this case only the fp_code option,
e.g.: [fp_code:'FP:MACCS'] or [fp_code:'FP:piDAPH4']
A list of fingperprints can be found using the command
ph4_FingerprintList []

EXAMPLE:

db_nb_mols_incommon [['db1.mdb', 'db2.mdb', 'db3.mdb'], 'MatrixOut.txt']

// for all databases that are open in the database viewer
db_nb_mols_incommon [dbv_KeyList[], 'MatrixOut.txt']

// use fingerprint similarity
db_nb_mols_incommon [dbv_KeyList[], 'Matrix.txt',
[fp_code:'FP:MACCS', sim_code:'tanimoto', thresh:0.95, fp_flag:1]]

// output average tanimoto between two databases based on fingerprints
db_compare_fp[['a.mdb', 'b.mdb', 'c.mdb'], 'Matrix.txt',
[fp_code:'FP:MACCS', sim_code:'tanimoto', average_Tani:1]]

db_uniqBits_fp [['db1.mdb', 'db2.mdb', db3.mdb], 'Matrix.txt',
[fp_code:'FP:piDAPH4']]

A demonstration of MOE's url reading and XML parsing functions: given a list of RSS feed URLs (hard coded in the constant URL), read the XML and show the feed's contents as a set of buttons displaying title and first 140 characters of each feed's description. Clicking on the buttons should launch the relevant webpage in your browser.

I've included the SVL Exchange as the first feed; if clicking the button doesn't jump you straight to the item, first log in to the Exchange so the cookie is set in your browser then try again.

MOE can import .png format images, so I might be able to get those into the buttons one day.

Usage:

Run this file, e.g. use the File | Open panel, browse to rss.svl and click OK.

Align protein chains in MOE using external alignment tool MUSCLE (http://www.drive5.com/muscle)

DESCRIPTION:
MOGUL descriptors - get descriptors for the number of rare and unusual bond angles, lengths and torsions.


USAGE:
1. Save and load this function
2. Compute Descriptors:
DBV | Compute | Descriptors | Calculate
3. Search in the descriptor list for mogul

Simultaneously browse corresponding residues of aligned sequences. This utility uses MOE's selection language to define what to browse, what to hide, and what to keep displayed. It should be useful in comparing corresponding structures such as various crystal structures belonging to the same protein family.

Calculate and display twist angle between two planes, each defined via three atoms

Print mutation codes of a given protein chain vs a reference chain. Output in SVL command window and text file.

Export RISM grids from *.RISM files to *.CCP4 format so they can be displayed via MOE Grid Analyzer etc.

Constraint (staple) alignment of protein chains to current state

Builds Kinase Homology Models using input files for template and sequences.

Allows various fusion modes described in the accompanying PDF document.

Graphical panel to run sdfrag from MOE. This generates a database of fragments using six possible methods listed on the panel.

Renders the protein backbone as a graphical object

This script display the force-field potential status of the system in SVL text window and creates graphical text labels. It can be used to verify
whether atoms are fixed, inert, or tethered.

Creates restraints on alpha helix hydrogen bonds

Cycle transparency setting of all or selected atoms. Automatically assigns CTRL+F8 on keyboard to execute itself.

This is a script to connect molecular fragments from a database with a connection point of a scaffold.

For defining the connection point in
the database a SMARTS pattern is used, whereas the first atom in the SMARTS filter is the connection point.

The scaffold has to be loaded in MOE and its connection point atom has to be named "CON".

Only fragments with exact one connection point are written to the output database.
An mseq integer for every entry in the input database is required!

For use in MOE_Batch use a tagged vector. Tags are:
'smfilter' for the SMARTS Filter.
'dbinuput' for the input file.
'dboutput' for the output file.

Allows multi-processor utilization for LowModeMD and Stochastic conformational search.

Instructions:

For single user mode, place these files under $HOME/svl/run and restart MOE. If unsure about location of the $HOME folder, please run MOE and in the SVL command window type HOME and press enter. Please notice that the 'svl' and 'run folders may not exist under $HOME and should be created.

To install for all users, place under $MOE/lib/svl/custom/run and restart MOE.

For a single machine multi-processor mode use the command:
MOE -mpu N

Where N is the number of parallel processes. For best performance this should be set equal to the number of physical processors available. For example for a dual core machine use MOE -mpu 2 . Similarly for a Core i5 or Core i7 machine N should be 4 or 6 respectively. For network multi-processor mode please refer to MOE help pages under 'Installing and Running MOE/smp' .

Next, you can run the conformational search via
MOE | Compute | Conformations | Conformational Search ...

and choose LowModeMD or Stochastic search to run the search.

The script uses a GUI to manage the scenes. To start the GUI, type "Scenes_Gui []" in the command line or you can add a button to the menu.
Once you have a view of the atoms in the graphical window that you want to save, you can add it to the list of the saved ones by clicking "Append" or "Prepend". The program will ask you for a unique name and the you'll see the name in the list.
Clicking on the name will bring the scene back on the screen: visibility, colors, surfaces, ribbons and so on.
If the animation between the scenes is too slow, you can deactivate it or decrease the sleep time between frames in the Options.
Scenes can be deleted, moved, modified and renamed easily with the buttons in the GUI.

When saving the scenes (Save button in the Gui), the script saves two files: one with all the molecules in the system (.moe extension) and another one with the information about the visibility and color of the objects (.mse extension); this means that you can load the scenes only through the GUI by clicking on Load.
This means also that, if you want to share the scenes, you have to send both files.

DESCRIPTION:
Calculates a descriptor, Plane of Best Fit, after Firth et al., Plane of best Fit: A Novel Method to Characterize the Three-Dimensionality of Molecules, JCIM 2012, 52, 2516-2525.

USAGE:

1. Save and load this function, e.g. save to $HOME/svl and restart MOE
2. In the database view menu, choose DBV | Compute | Descriptors | Calculate
3. Look for a descriptor called PBF (type it into the Filter: line)

Description
An interface to allow substructure searching of a MOE database with the substructure defined by an external sketcher rather than SMARTS.

When this tool is run for the first time, the ten slots for stored searches will be empty. Click on New Sketch to launch the sketcher and create a substructure search; once transferred back to MOE, the sketch will appear in the first available slot. Click on the sketch search the selected MOE database for that substructure. If all ten slots are full, the first slot will be emptied to make room. Alternatively, Ctrl-click on a slot to remove that sketch from the history.

Usage
Run this SVL file, e.g. use the MOE | File | Open panel, browse to db_substructuresketch.svl and click OK.

The MOE Protein Properties application calculates a set of key physicochemical protein properties that are of interest in protein engineering applications.

If there is more than one receptor (protein) chain, a protein-protein interface interaction report is generated. For each pair of interacting residues there is a listing of the atom, residue and UID along with the type of interaction (HB, pi-pi etc), the distance and the strength.

To install the MOE/web application, unzip it into your $MOE/web folder and (re)start the MOE/web server. Alternatively unzip it to a different directory, e.g. webapps, then start the MOE/web server with the -app option:

$MOE/bin/moebeb -app webapps

The scripts finds all scaffold in a database, writes them to a separate database and calculates the complexity and cyclicity of the molecules like described in the paper. Also the frequency of the scaffold in the original database is written to the scaffold database.

Usage: load 'sca.svl' then use: db_SCAclassify['databasename',opt] where opt = 0 means that chiral constraints are not taken into account and 1 means they are taken into account. You may also run the script in GUI mode. In this case you will be prompted for a database and you can select how chiral constraints should be handled.

The MOE SA/Report tool uses a file ($MOE/web/sareport/units.txt) to determine both the units (eg. pIC50, IC50, percent) that a column in your data can have, and also regular expression style patterns in the column name that enable MOE to guess its type (eg. *PIC50 implies it is a PIC50). The format of the units.txt file is a little cryptic and difficult to edit by hand. Hence, this module allows the user to define new units or change the regular expression patterns, names or matching order of already defined units.

Place the module in $MOE/lib/svl/custom so it is loaded on startup, and run it without arguments from the SVL command line:
SAreport_range_editor [];

The GUI allows you to define, for each unit being edited, the following information:

* Autotype fields: a limited regular expression syntax consisting of asterisks at the beginning or end (or both) of a match string, and with multiple patterns separated by commas. Eg. *pK*,*pXC50*
* Name: a name for the new unit
* Type: activity (eligible for analysis by SAReport) or property (not eligible). Most new units should be defined as activities (even "properties" like MW or ClogP), unless you don't want to allow SA/Report to analyze them.
* Minimum, Mid-point and Maximum numeric points of the range for this data type. These will correspond to red, yellow and green colors respectively in the standard SA/REport three-color scheme, and values between these three thresholds will be interpolated.
* Multiplier/Offset/Function: scaling, addition and functional (eg. log) transforms applied while interpolating numbers between low/medium/high thresholds. I.e. y = f(ax+b)

PS: this module is referenced in the MOE UGM poster and ACS Philadelphia talk given by Deepak Bandyopadhyay of GSK in 2012.

When run, a menu is created in the Sequence Editor that allows the user to choose the set of residues to be acted upon and the renumbering action (whether it is on the alignment position or the residue ID) to be performed. This is useful for preparing a set of protein-ligand complexes for input into PLIF and Ligand Interactions, as those applications rely on a consistent rUID scheme across all receptors.

Additionally, antibody structures can be renumbered to a canonical scheme.

To use the code, run it, e.g. use the MOE | File | Open panel, browse to resIDfns.svl and click OK, or use this SVL command:

run 'renumber_residues.svl'

Classification of Organic Molecules by Molecular Quantum Numbers
Kong T. Nguyen, Lorenz C. Blum, Ruud van Deursen, and Jean-Louis Reymond
ChemMedChem (2009), 4:1803-1805
DOI: 10.1002/cmdc.200900317

USAGE:
1. Save and load this function
2. The 42 mqn_* Descriptors will be available in MOE

Run the file. A task called 'SurfSync' will show in the Cancel menu (and stay there). If you have a surface around some atoms and then move the atoms, the surface will move with them.

- If only a subset of defining atoms is moved, the surface will try to superpose itself on all the defining atoms but will not change shape. (To change the shape of the surface, you must select it in the surfmap panel and press Apply).

- The surfaces will follow atoms only while the task (named 'SurfSync') is running.

- A new instance of SurfSync will destroy any previous instances and will reset all surfaces according to the current positions of atoms that define them.

DESCRIPTION:

Prediction of three-dimensional molecular structures using information
from infrared spectra
Hemmer, MC and Gasteiger, J
Analytica Chimica Acta 420 (2000) 145-154

NOTE: Different scaling factors can be found in the literature.
This code does not apply any scaling facor, i.e. both, factor
"f" and "B", as sometimes refered to in the literature,
are both set to one.



This code implements RDF descriptors in a variety of flavours:
1.) Distance:
- topological distances: tRDF*
- geometric distances: gRDF*

2.) Properties encoded:
- none/unweighted: ?RDF*u
- atomic mass: ?RDF*m
- atomic volume: ?RDF*v
- electronegativity: ?RDF*e
- SMR increments: ?RDF*s
- SLogP increments: ?RDF*L
- VSA values: ?RDF*a
- partial charge: Q_?RDF*c
- positive partial charge: Q_?RDF*c++
- negative partial charges: Q_?RDF*c--
(NOTE: partial charges need to be calculated
via DBV|Compute|Molecule|Partial Charge)

- E-Hueckel based strength parameters for hydrogen bond donors
or acceptors: ?RDF*ha and ?RDF*hb



USAGE:

1. Save and load this file
2. The descriptors will be available via DBV|Compute|Descriptors...
3. You can find them by typing "tRDF" or "gRDF" in the
filtering textbox

DESCRIPTION:

Moreau/Broto's Autocorrelation vector descriptors
Moreau, G.; Broto, P. NouV. J. Chim. 1980, 4, 359-360.



This code implements ATS descriptors in a variety of flavours:
1.) Distance:
- topological distances: tATS*
- geometric distances: gATS*

2.) Properties encoded:
- none/unweighted: ?ATS*u
- atomic mass: ?ATS*m
- atomic volume: ?ATS*v
- electronegativity: ?ATS*e
- SMR increments: ?ATS*s
- SLogP increments: ?ATS*L
- VSA values: ?ATS*a
- partial charge: Q_?ATS*c
- positive partial charge: Q_?ATS*c++
- negative partial charges: Q_?ATS*c--
(NOTE: partial charges need to be calculated
via DBV|Compute|Molecule|Partial Charge)

- E-Hueckel based strength parameters for hydrogen bond donors
or acceptors: ?ATS*ha and ?ATS*hb

USAGE:

1. Save and load this file
2. The descriptors will be available via DBV|Compute|Descriptors...
3. You can find them by typing "ATS" in the filtering textbox

Calculates the Sterimol descriptors L, B1, B2, B3, B4, B5.
Load file into MOE. The six new descriptors will appear in the Descriptor panel.
Descriptors are calculated for the whole molecule if no attachment point (A0) is found. If an attachment point is found L is measured along the substitution axis and B1-B4 orthogonal to the substitution axis. (see picture)

MOE-FMOutil consists of several tools and DEMO files.
> Structure preparations (Truncation around the ligand, etc.)
> Input-file preparation (Suport: PIEDA, FMO/PCM, FMO-D optimization, partial MP2 PIE etc.)
> Read and validate GAMESS-FMO inputfile
> Import optimize coodinates into MDB file
> Output Analysis (Summation of Pair Interaction Energy for insisted fragments, Visualize Charge transfer, etc.)

USAGE:
Load fmoutil.svl and type FMOutil []

This script draws cartoons of molecular fragments, to provide non-atomic visualization of macromolecules. The molecular fragments and the corresponding cartoons are read from a customizable dictionary. The current dictionary includes some carbohydrates and nucleic acids.

As described in the accompanying documentation, SEBAS is a program for automatically assigning NMR spectra to a given amino acid sequence. Peak files in standard Sparky and/or NMRPipe format are curated with possible interaction with the user. The curated files are used to construct generic spin systems (GS), and then these are assigned to the given sequence (in FASTA format). Possibly multiple assignments are written out as an Excel spreadsheet for easy examination by the user. Several different assignment algorithms are offered, some being suited to extensive examination of all possible assignments of high quality data, and others are better suited to finding any sort of consensus of assignment from noisier data coming from fewer different triple resonance experiments.

Usage: Simply choose MOE | File | Open then double click on
visualize_dipole.svl

For it to be available for use by all of the users of MOE, create a directory called $MOE/lib/svl/custom/run and save this file to that directory. eg on Windows and Linux, this might be:

c:\moe\lib\svl\custom\run\dbimport.svl
/usr/local/moe/lib/svl/custom/run/dbimport.svl

For the script to be available for use by a single user, then it is best to create a directory in your home directory called svl, and create a directory called run in that directory. Save this file into that directory. eg on Windows and Linux, this might be:

c:\Documents and Settings\user_name\svl\run\dbimport.svl
/usr/people/user_name/svl/run/dbimport.svl

DESCRIPTION:

This plugin computes the change in exposed surface area upon ligand binding to the receptor.

USAGE:

Select the ligand atoms and run this code. The results will be printed into the Commands Window: exposed surface area for a 'free' ligand in its current conformation, and the exposed surface are for the ligand in its complex with the receptor.

Description:
Calculate the ROC curve based on a real activity field and a
predicted activity field. The user has the option of ascending ranks
(small is good) or descending ranks (big is good - default) in both
the real and predicted activity fields. If the used specified
a prediction threshold, a confusion matrix of predicted active/
inactive vs. actual active/inactive is calculated based on the
threshold.

Usage:
(1) Run this file (TED | SVL | Save and Run...)
(2) You will be prompted for a database
(3) The ROC Plot panel will open.
(4) Use the option pull-downs to specify the real and
predicted activity fields. Set the thresholds for each if appropriate.
(5) Press 'Calculate' to perform the matrix calculation.
(6) Use the Radio button at the bottom of the panel to
display the results as either absolute numbers or percentages.
(7) Various statistics are also displayed.
(8) Press the 'Dump Plot to Text' to produce a tab-delimited
file that can be imported in to Excel

ROC "receiver operating characteristic curve" see following links

http:www.statsoft.com/textbook/glosr.html#ROC
http:www.anaesthetist.com/mnm/stats/roc/

Given a database with predicted active and true activity,
the function plots an ROC curve. Sample data from one of the
above links is given at the end of this file.

Usage:
(1) Load this function
(2) Open the database with the activities and predictions.
(3) Enter

svl> ROC_Plot[]

at the SVL command line

(4) Choose the predicted active and the true active fields
(5) Set the thresholds so that compounds where

'Activity Field >= Activity Threshold'

'are deemed "active"' and compounds where

'Prediction Fiels >= Prediction Treshold'

are deemed 'predicted active'

(6) Press 'Plot' to plot the ROC curve.

HITS are CORRECT PREDICTIONS
MISSES are INCORRECT PREDICTIONS

The area beneath the curve is also reported

(7) Press 'Report' to produce a text report including
area calculation and TPF vs. FPF data

Start the program with "Run SVL File". The program will ask for a .kont file. All grids contained in the .kont file will be exported as .grid files in MOE grid format into the directory of the .kont file.
The exported grid files can be visualized with MOE's Surfaces and Maps.

Please note: The binary GRID format conversion is still "experimental". If you encounter problems please contact gkirsten @ chemcomp.com

DESCRIPTION:
A low footprint GUI to set various atom properties on a per chain basis, e.g. visibility, render mode. The interface is based on an original idea by Francis Atkinson and Andrew Leach (GSK).

The application allows you to set the visibility for particular chains (or show only one chain, or show all or no chains), set chain names and tags, select individual or all chains, change the rendering mode from line bonds to stick bonds and alter the potential settings (none, fixed, inert). Which hydrogen atoms are shown can also be changed, and carbon atoms can be coloured by chain or by unique chain tags.

The application will remember if you have hidden atoms in a particular chain, e.g. if you are showing only the pocket atoms from a receptor. If you wish to reset the set of hidden atoms for a particular chain, simply hold down the Shift key and click on the visibility toggle button. If you wish to reset the whole system's hidden status, Shift click the "Show" button at the top.

See the bubble help for the action of each button by holding the mouse pointer over it for a few seconds.

Note:
If there are more than a certain number of chains, the buttons are spread over a set of pages, controlled by < and > buttons in the bottom right of the panel. Currently the maximum number of rows, 18, is hard coded within the script in line 88.

You can make the Chain Manager window stay in the same place when it is recreated after creating or removing chains, set the SAVEPOS constant to 1 in line 89.

USAGE:

1. Save this function
2. Run the SVL, e.g. use the MOE | File | Open panel, browse to chainmanager.svl and click OK.

MOE2005.06 does not support installation on Solaris 10 (due to the fact it was released prior to Solaris 10). The included files will be necessary to install MOE on a machine with Solaris 10, as well as register the license manager as a Solaris 10 service.

Files included in MOE_Solaris.tar:
readme.txt - these instructions
install.sh - MOE install script
arch.sh - MOE architecture redirection script
lmgrd.xml - lmgrd service file

Instructions:
1. Install MOE with the included install.sh file with the -i option. For example, if the install.sh file is at the current directory, and the MOE CD is mounted, type:
./install.sh -i /cdrom/moe_2005_06

2. Replace $MOE/bin/arch.sh with the included arch.sh file.

3. Install the license manager as follows:

a. Edit the paths to your MOE/bin, license file, and log file directories on each of the "exec" lines. In this particular file, MOE is installed at /usr/moe, and the log file resides in the root MOE folder.

b. Copy lmgrd.xml to /var/svc/manifest/application/lmgrd.xml (this is not essential, but this is where all of the other service manifest reside)

c. Type "svccfg import /var/svc/manifest/application/lmgrd.xml"

d. Type "svcadm enable application/lmgrd"

e. The service is now registered and started. It will start on each boot. To stop the license manager, type "svcadm disable application/lmgrd", and to reread the license file, type "svcadm refresh application/lmgrd". To restart the license manager, type "svcadm restart application/lmgrd"

OR
a. Copy the "lmgrd_init" file from the SGI architecture included in MOE to /etc/rc3.d/S95lmgrd_init (replacing the appropriate variables in the script).

b. Treat it as you would any other legacy (pre-Solaris 10) service.


Known Issues:
1. lmgrd is run as root. This can be easily changed by adding "username -c " to the beginning of the commands preceded by "exec" in lmgrd.xml

2. The log file is continuously written, so it may get too large. Simply make a copy of the current log, delete it, and restart the license manager.

This application facilitates exploration of the conformations generated via a conformational search method, including LowModeMD.

It scans a database of conformers and displays the global minimum geometry in the MOE window. Then it display a 2D plot of observed values of 'important'
dihedrals. Next, it allows selection of one or more observed dihedrals and displays the values and some statistics in the MOE window, and selects the
corresponding conformations in the database viewer.

Fingerprint using pharmacophore typing with distances calculated over a reduced molecular graph, after JCIM 2006, 46, 208-220 (DOI: 10.1021/ci050457y)

USAGE:

1. Save and load this function e.g. using MOE | File | Open, browse to ph4erg.svl and click "Load SVL".
2. Calculate fingerprints using DBV | Compute | Fingerprints: the new fingerprint will be available, called FP:ErG. Searches can be performed using the DBV | File | Similarity Search tool.

At the moment the amount of fuzzy incrementation has a default of 0.3. To alter this, set the value of a variable called ERG_FUZZY at the command line, for example:

svl> ERG_FUZZY = 0.2

Fuzzy incrementation as implemented here differs slightly from that described in the paper. The paper suggests that for a particular bit that has e.g. three matching paths in the molecule, that bit and its two neighbours would be:

3*FUZZY, 3, 3*FUZZY

In this implementation that has been changed to:

FUZZY, 1 + FUZZY, FUZZY

In addition each bit can only be incremented once, so if bits A and C were matched in the molecule, bit B would only be incremented once. This has been done in order for the similarity scores to more closely match the ones presented in the paper. If the scheme described in the paper is used, similarity scores are too low due to the large value in the main bit.

The correlation between the scores from this SVL and the paper are still not 100% - I am looking into the reason for this. In my tests with thiscode, the similarity scores are always slightly too low.

The collapse of hydrophobic ring systems into a single Hf feature will occur if any two or more ring systems share more than two atoms.

There is no limit to the number of flip-flop atoms (atoms that are both donor and acceptor). Instead, a descriptor is included in the code (accessible via the MOE tool DBV | Compute | Descriptors, called "flipflop". This allows the user to filter their databases manually prior to fingerprint calculation.

_Normalising the fingerprint_
A database can be used to calculate a normalisation vector, that is a set of numbers by which each fingerprint bit will be multiplied when comparing similarities. This allows the effect of overly-represented bits to be reduced. To create a file of normals, use this command:

svl> CalculateErGNormalisation [mdbfile, normalsfile]

Where both mdbfile and normalsfile are tokens containing the filename of the input MDB file and the output text file. If normalsfile is null then the filename will be {mdbfile}_ergnormals.txt.

A file of normals can be "loaded" into MOE using this command:

svl> LoadErGNormals normalsfile

This will then be used when performing similarity searches.

_Visualisation and Reverse Fingerprinting_

Diagnostic tools are available in this code:

1. If you load a molecule into the main MOE window and run this SVL (e.g. use the MOE | File | Open panel, browse to ph4erg.svl and click OK) then the feature types assigned to the molecule will be drawn in the system. In addition the list of features and the distance matrix of features will be printed in the SVL Commands Window. Only the first feature distance matrix is printed in the case that the molecule contains flipflop atoms. If there are flip flop atoms then for each possibility a different graphic object is created. These can be viewed individually using the Graphics Object Manager.

2. With this SVL file loaded into MOE you can see the scores for each property-distance-property triplet with a molecule loaded in MOE using this command:

svl> Reverse_ErG_Fingerprint ph4_ErG_Fingerprint Atoms[]

Fuzzy incrementation is applied, but weighting is not.

This code is also compatible with the MI reverse fingerprinting package, although as of 10th Sept 2009 it does not elaborate flip flop atoms.

3. ErG Bit Descriptors
This will create a new column for each of the bits in the ErG fingerprint (default settings will create 315 new fields). These are then filled with the counts for each bit. Fuzzy incrementation is applied, but only once, i.e. if a bit is matched three times, that bit and its two neighbours would be:

1*FUZZY, 3, 1*FUZZY

To use this tool, enter the following SVL Command:

svl> ErG_Bit_Descriptors mdb

where mdb is a token containing the filename of an MDB, e.g. 'test.mdb'.

4. ErG Bit Descriptor Heat Map
This will load the descriptors calculated using the previous command from the database and create a graphic object showing the counts as a heat map coloured using a black-body palette.

svl> ErG_Bit_HeatMap mdb

The graphic can be saved in various image formats (JPG, PNG, etc) by clicking the Export button in the Heat Map panel.

Wildman & Crippen, J. Chem. Inf. Comput. Sci., 43, 629-636 (2003).

Given the 3xN array of Cartesian coordinates of N atoms plus N associated "weights" or property values of the atoms, calculate a chirality value for the distribution of the weights in space. The weights must be positive.

This corresponds to CIP chirality when the weights are the atomic weights or priorities of the four substituents of a chiral center, and the coordinates are the coordinates of the substituents. If the four substituents are coplanar or if two weights are the same, the calculated chirality is zero. The mirror image of the configuration gives the opposite sign chirality.

The difference is that this chirality measure is a continuous function of the coordinates and the weights of four or more points. Not only does it give the same value when the coordinates are transformed by a translation or proper rotation, but the value is unchanged when the coordinates or the weights are multiplied by a positive constant.

SMARTCyp predicts positions in drug-like molecules that are likely sites for metabolism by cytochromes P450. It acts as a wrapper for the SMARTCyp method within MOE.

SMARTCypMOE processes single molecules loaded in the main MOE window or groups of molecules stored in a MOE molecular database. Molecular databases can be processed interactively from MOE’s Database Viewer or in batch mode via the moebatch program. Reactivity ranks and scores are added to the molecular database and the results can be viewed graphically in the MOE Database Browser.

Using a heme model system,Rydberg et al. computed transition state energies for more than 200 fragments and from these energies and one accessibility-like descriptor a very simple scoring function was deduced which is used to score and rank atoms in a molecule. The model covers all common reactions in P450s and has a good performance for all P450 isoforms except 2C9 and 2D6.

Score = Reactivity – 8 x Accessbility

Where Reactivity is the activation energy in kJ/mol and the accessbility is the span descriptor which measures the relative topological distance to the end of the molecule (which always is between 0.5 and 1).

Prerequisites:
The SMARTCyp executable jar file which can be downloaded at http://www.farma.ku.dk/smartcyp/

The distribution is a zip archive containing the SVL script (cypscore.svl), entries for the moe-menus file (moe-menus-cypscore) and a manual (SMARTCypMOE.pdf), which provides full details on how to install and use the script.

References:
Patrik Rydberg, David E. Gloriam, Jed Zaretzki, Curt Breneman and Lars Olsen (2010) "SMARTCyp: A 2D Method for Prediction of Cytochrome P450-Mediated Drug Metabolism" ACS Med. Chem. Lett., 1, 96-100
Patrik Rydberg, David Gloriam and Lars Olsen (2010) "The SMARTCyp cytochrome P450 metabolism prediction server" Bioinformatics, DOI: 10.1093/bioinformat

Divides numerical data into bins defined by the lower bounds, the upper bonds
and width. Optionally, bins are divided by the number of bins if the width is
not set. The values normalized to the lower bounds of bins will be written in
the input MDB and the lower bounds and the frequencies of bins will be output
to the other MDB.

The histogram of the values of fieldname and their lower bounds
within their bins in Gaussian distribution will be plotted.

Installation:

(a) Place dbv_histogram.svl in $MOE/lib/svl/custom/ directory
and invoke MOE for the general use
(b) Place dbv_histogram.svl in $HOME/svl/ directory
and invoke MOE for the personal use
(c) From the MOE | File | Open panel, select this file
and click the "Load" button for the temporal use.

USAGE:

moe> db_Histogram [filename: '/home/userid/work/filename.mdb', \
nfname: 'SlogP', \
esel: 0, lb: -9, ub: 9, w: 0.5]
('\' stands for the continuation, i.e., enter these in a single line)

or

1. Select a field in MDB on DBV,
2. DBV | CLI | dbv_Histogram [lb: -9, ub: 9, w: 0.5]

where the parameters are:

filename: the mdb filename
If not specified, the filename of the default view will
be used.
nfname: the numeric field name to calculate
The default is the first numerical field name.
esel: calculate only on the selected entries
lb: the value of the lower bound
The default is the minimum value in the field.
ub: the value of the upper bound
The default is the maximum value in the field.
w: the width of the bins
If not specified, (ub - lb) / nbin will be used.
nbin: the number of bins
The default value is 10.
If w is specified, (ub - lb) / w will be used.
openodb: If on, open the output database in database viewer
The default is 1 (on) for GUI
and 0 (off) for batch mode.
drawhis: If on, draw the histogram
The default is 1 (on) for GUI
and 0 (off) for batch mode.

e.g.,

DBV | CLI | dbv_Histogram [nfname: 'SlogP', lb: -9, ub: 9, w: 0.5]

or

1. Select a field,
2. DBV | CLI | dbv_Histogram []
3. Set optional parameters and click OK

if w is given, nbin will be ignored.

The outputs are:
filename.mdb:
l_fieldname: the lower bound of the descriptor in each bin

filename_fieldname_histogram.mdb
l_fieldname: the lower bound of the descriptor in each bin
f_fieldname: the frequencies of the entry in each bin

The histogram between the values of fieldname and their lower bounds
within their bins in Gaussian distribution will be plotted.

There are three main routines you can call :

gamin 'gamess.log'
Update co-ordinates and charges of current Moe molecule with those from the gamess output file. You MUST have the appropriate molecule already read into Moe or the results will be unpredictable.

gampotgui 'gamess.log'
Read the molecular electrostatic potential (MEP) from the gamess output. Contour the MEP on a grid around the molecule using a Gui to control the three contour levels. Works with the 2010 version of Gamess (WinGamess 10).

gamess_input_gui []
Produces a Gui to generate basic Gamess input files. input is a 6-31G** Hartree Fock calculation. You will need to change the basis set manually but you can select an option in the Gui to perform a B3-LYP calculation. The input file is generated when you enter a file name and press return. A .moe file of the current structure is saved at the same time. This is needed when you come to analyse the results.

There is an option to select a Z-matrix (i.e. internal co-ordinate) input file / optimisation rather than the cartesian input normally produced by Moe. This can lead to faster optimisation - I achieved optimisation of 3-Me pyridine in 7 cycles rather than 42 with cartesian, so quite a useful saving.

It is possible to select a partial Z-matrix optimisation in which some degrees of freedom (DoFs) are frozen at their input values. This is handy for calculating strain energies, where you can freeze all torsion angles and allow the bonds and angles to relax. A Gui is generated which allows interactive selection of the frozen DoFs.

The input file for the MEP calculation can also be generated. This option generates a grid of points around the molecule which can be manipulated using a Gui to give the size of grid you need to observe features of interest.

Note : optimisation will move the molecule slightly from the starting point, so you may want to optimise first and then do a second calculation without opt. to calculate your MEP.
Z-matrix optimisation moves the atoms to the co-ordinate origin, so your grid WILL be in the wrong place.

The default version of Gamess limits you to only 100 grid points. You will need to recompile it. Go to the source directory in the unpacked distribution and edit prplib.src. Global replace MXPTPT=100 with MAXPTPT=99999. Then save and follow the compilation instructions. (Easier on linux than windows as you will need compilers, which only come free with good operating systems).

How to use the script:

In MOE command window:
load 'selectNentry.svl'

Individual functions can be used as follows:

selectNtry [num_entry, which_entry, entry_option]
examples:
selectNtry [20,0,0] for selecting every 20th entry of the database
selectNtry [20,1,0] for selecting every 21st entry of the database
selectNtry [20,2,0] for selecting every 22nd entry of the database
selectNtry [5,0,1] for selecting every 5th entry and creating a new group for it


indexgroup [fieldname, indexname]
Where:
fieldname: name of the field that needs to be indexed
indexname: name of the index field that will be created
example: indexgroup ['mseq','Rank']

selectdups [fieldname]
Where:
fieldname: name of the field in which duplicates will be selected
example: selectdups ['CMPDID']
Note: The field should be sorted before running the command

replacentry [fieldname, searchstrng, replacestrng, esel]
Where:
fieldname: name of the field in which the string to be replaced will be searched for.
searchstrng: this string will be searched in the fieldname of the database
replacestrng: this string will replace the search string in all entries of the database
esel: if set to 1 will operate on selected entries otherwise entire database
example: findentry ['HIT_TYPE',"HTS1","HTS2",1]
Note: string for searching and replacement should be in double quotes "

DESCRIPTION
This utility allows import/exporting MOE-Restraints from/to tab-separated-text files. It also compares the existing MOE-restraints vs. the currently loaded
structures and saves the 'unsatisfied' restraints to a text file.

This function allows BREED operations on aligned structures with general formula of A-B-C where B is a ring system. For n input structures, it generates all A(n)-B(n)-C(n) combinations. A demo input file is provided (_rb_demo.moe).

The -protomer switch in sdwash creates the tautomers along with protomers.In case enumerating new tautomers from the input structures is undesirable this program can be used to mark the protomers in the enumerated database, so that they can be separated from the tautomers.

This is a set of scripts that allow the CEPOS InSilico Ltd. programs ParaSurf and ParaFit to be run from within MOE. ParaSurf generates molecular surface properties based on semi-empirical molecular orbital calculations. ParaFit superposes and compares molecules using the spherical harmonic expansions of the molecular surface and properties calculated by ParaSurf. The scripts enable the user to run ParaSurf on single molecules or molecular databases; use ParaSurf descriptors and rotationally invariant fingerprints (RIFs) in MOE QuaSAR; visualize ParaSurf surfaces and properties; align molecules with ParaFit based on their ParaSurf properties; and produce ParaFit 'canonical' alignments.

Prerequisites:
ParaSurf and ParaFit are available from CEPOS InSilico Ltd. (www.ceposinsilico.com). ParaSurf requires the results of semi-empirical molecular orbital calculations, which can be provided either by VAMP or by a public domain version of MOPAC that is available for free download from www.ceposinsilico.com.

The distribution is a zip archive containing the SVL scripts, entries for the moe-menus file and a manual in PDF format that provides full details on how to install and use the scripts.

Radial Distribution Function Descriptors
based upon:
Hemmer, M. C.; Steinhauer, V.; Gasteiger, J.
J. Vibrat. Spect., (1999), 19, 151-164
The Coefficients f and B have been derived from fits against literature data.

This is an independent implementation (D.Bondarev) of the reference below:

J.Chem.Inf.Comput.Sci. 43(4), pp. 1151-1157, 2003
Design and Evaluation of a Molecular Fingerprint Involving the Transformation of Property Descriptor Values into a Binary Classification Scheme. Ling Xue, Jeffrey Godden, Florence Stahura, Jurgen Bajorath.

USAGE:
Load this file with the MOE File Open panel (click OK).

You can then use MOE fingerprint tools as usual.

Compute the MPMFP fingerprint or use any fingerprint tools (including any scripts from this website). The fingerprint code is MPMFP, the similarity metric code is avg-tanimoto.

This is an implementation of the CypScore method devised by Hennemann et al. for the prediction of cytochrome P450 reactive sites in drug-like molecules. The script processes single molecules loaded in the main MOE window or groups of molecules stored in a molecular database. Molecular databases can be processed interactively from MOE’s Database Viewer or in batch mode via the moebatch program. Reactivity scores are added to the molecular database and the results can be viewed graphically in the MOE Database Browser.

For six oxidation reactions, Hennemann et al. derived linear models predicting an atom’s reactivity from the values of descriptors calculated by ParaSurf, with each model producing reactivity predictions on an individual scale. The original paper described a procedure to establish a common reactivity scale for all the models, providing a CypScore value for each atom in the range from 0 (stable) to 100 (reactive), and implemented this procedure for a proprietary data set. Unfortunately, Hennemann et al. do not include sufficient detail to define completely the relationship between the original model values and the scaled scores, and it is not possible to reproduce the missing information as the data from which it was derived is not available publicly. To overcome this, this script employs a common reactivity scale that is a close approximation to the one used by Hennemann et al., derived from secondary information in the original paper through a process described in the accompanying manual. Consequently, the results obtained will approximate closely, but not match exactly, those in the original publication.

Prerequisites:
The CypScore predictions are based on molecular surface descriptors calculated by ParaSurf, which is available from CEPOS Insilico Ltd. (www.ceposinsilico.com). ParaSurf requires the results of semi-empirical molecular orbital calculations, which can be provided either by VAMP or by a public domain version of MOPAC that is available for free download from www.ceposinsilico.com.

The distribution is a zip archive containing the SVL script (cypscore.svl), entries for the moe-menus file (moe-menus-cypscore) and a manual (CypScoreMOE.pdf), which provides full details on how to install and use the script.

Reference:
M. Hennemann, A. Friedl, M. Lobell, J. Keldenich, A. Hillisch, T. Clark and A. H. Goller. CypScore: quantitative prediction of reactivity toward cytochromes P450 based on semiempirical molecular orbital theory. ChemMedChem, 2009, 4

Usage:
(1) Load this function
(2) Enter

svl> GA[]

at the SVL command line

You can get linear model, perform cross validation and
save your models in mdb file format from this panel.

To evaluate new molecules with your model, enter

svl> ModelCalc[]

-----------
TUTORIAL:

'$MOE/sample/mol/blood_brain.mdb' may be used as a tutorial.

(1)Calculate 32 VSA descriptors and TPSA
(2)Run QuaSAR-Evolution with default parameters
Equation file 'ga_*_eq.mdb' and trajectory file 'ga_*_trj.mdb' will open
(3)Plot SMR_VSA4 and SlogP_VSA0 in plot area of trajectory file.
and observe how equations evolve
(4)Calculate TPSA descriptor
(5)Click 'CONTINUE'
(6)Plot TPSA and observe equations
You will see that the combination of SMR_VSA4
and TPSA is very important for good model

Given a sample of (x,y) points where x values are not necessarily uniformly spaced, kappaxyval [x, y] returns 0 if there are just two tight clusters of points, 1 if y is a linear function of x, >1 if y is a nonlinear function of x and/or there is some noise, and 2 if the points are a uniform random scatter inside some rectangular region. This can be handy for automatically detecting linear and nonlinear relations in samples of points, whereas the usual correlation coefficient may be nearly zero for nonlinear relations. This also lends itself to fitting nonlinear and/or noisy samples of points. yfit = fitKxy [x,y] gives a fitted y value for each sampled x value. The yfit values form a relatively smooth curve.

All this is described in: G.M. Crippen, Computational Biology and Chemistry, 33, 357-360 (2009).

Simply run the function, e.g. by choosing MOE | File | Open then double click on torsion_editor.svl. Once the application is running one needs to select two atoms that constitute a rotatable bond. This is done by selecting the first atom then the second atom (by holding down the shift key). Atoms attached to the second selected atom will be the ones that move when the bond is rotated (Alt - left drag).

- to visualize which half of the molecule will rotate once a selection is made, toggle on Show bond arrow. The arrowhead points toward the movable portion

- to reverse which half of the molecule will rotate click on the Reverse button

- to see the current value of the selected torsion toggle on Show meters. Note the values run from 180 degrees to -180 degrees in a clockwise sense when the rotatable atom is in front

- to add a torsion to the list, select a valid set of atoms then click on Create. The atoms that comprise the torsion, the original angle (when the torsion was added to the list) and the current angle is shown. Now whenever this value changes, it will be updated in the list (and the meters as well if Show Meters is toggled on)

- to reselect a torsion in the list, simply click on it and the corresponding atoms will be selected in MOE

- to set a torsion back to its original value, select it in the list and click Reset. (The original torsion angle is value that existed when the torsion was added to the list.)

- to delete a torsion form the list, select it then click on the Delete button

Important Note: This application works by FIXING (MOE | Edit | Potential | Fix) atoms in space to force the required portion of the molecule to rotate. If you decide to minimize the system while the editor is running please make sure that no atoms are selected as the application may create sets of fixed atoms.

1. Save this script to your hard disk.

2. (a) From the MOE | File | Open panel, select the db_mol_setenv.svl
file and double click with the left mouse button.

(b) Alternatively from the SVL command line, type:

run 'db_mol_setenv.svl'

3. To make these changes from MOE/batch, use the following syntax:

$MOE/bin/moebatch -exit -exec "run ['db_mol_setenv.svl', ['dock.mdb',
['ligmol', 'ligmol1'], ['recmol', 'recmol1], ['score1',[]]]]

Score the sequences in the main MOE window by similarity. The similarity matrix is used from the default setting for the MOE system (blosum62 by default), but you can pick a different one from the list. The identity and similarity between two chains is shown in the panel. If you select the "Color Residues" checkbox then the residues will be colored so that similar residues are green and dissimilar residues are red. If you have any chains selected, the colouring and sequence similarity and identity scores are restricted to selected chains.

If you click the "Report" button, then an HTML file is written out with the sequence alignment and the similarity matrix, comparing all the sequences against all the others.

This matrix is not symmetrical - in one half of the matrix the scores are divided by the number of amino acids in one sequence. In the other half of of the matrix, the scores are divided by the length of the other sequence.

If you have any residues selected in the Sequence Editor when you click the "Report" button, then the rest of the residues in the same column will be selected. The similarity matrix in the HTML file will only consider the residues in those columns. In the sequence alignment, the residues which had been selected will be shown in bold and upper case. Unselected residues will be in the
normal font and lower case.

The rPos residue numbers for the whole alignment and the rUID numbers for the first sequence are shown at the top of the alignment.

In the rPos residue numbering the residue number is shown for the residues where the rPos is a multiple of 10. The residue number is positioned so that the unit digit from the number is above the residue. Residues where the rPos number ends in 5 are shown with a : character.

Similarly, the rUID residue numbers are shown for residues where the rUID is a multiple of 10. If there is an rINS insertion character for the residue, then this is shown. Otherwise if the rUID number ends in 5, then this is shown with a : character. If there is another residue it is shown with a - character. If there is a gap in the sequence for that structure, then a space is shown. rUID residue numbers may not be consecutive or unique.

For instance if you had a stretch of amino acids where the rLetter, rPos and rUID values were as below:

rLetter rPos rUID
A 1 18
C 2 19
D 3 19A
E 4 20
F 5 20A
G 6 25
H 7 26
I 8 27
K 9 27A
L 10 27B
M 11 27C
N 12 28
P 13 29
Q 17 30
R 18 31
S 19 32
T 20 33

----:---10----:---20
--20A:--ABC-- 30---
ACDEFGHIKLMNP---QRST

Note:

The rINS for residue 19A is not shown as it is overwritten by the label for residue 20. The rINS value is shown for residue 20A. There are no residues where the rUID is 21, 22, 23 or 24 so these are missing from the rUID numbering line. There are 3 gaps in the alignment, but there are only 2 spaces shown in the rUID numbering line as the label for residue 30 overwrites one of the spaces.

USAGE

1. Save this file to your hard disk.

For it to be available for use by all of the users of MOE, create a directory called $MOE/lib/svl/patch/run and save this file to that directory. eg on Windows and Linux, this might be:

c:\moe\lib\svl\patch\run\seq_sim_monitor.svl
/usr/local/moe/lib/svl/patch/run/seq_sim_monitor.svl

For the script to be available for use by a single user, then it is best to create a directory in your home directory called svl, and create a directory called run in that directory. Save this file into that directory. eg on Windows and Linux, this might be:

c:\Documents and Settings\user_name\svl\run\seq_sim_monitor.svl
/usr/people/user_name/svl/run/seq_sim_monitor.svl

2. Run this file.
(a) From the MOE | File | Open panel, select this file and click the "Run" button.

(b) At the SVL command line, type:

run 'seq_sim_monitor.svl

A Method for the Alignment of Molecular Strucures: Maximizing Electrostatic and Steric Overlap.

1. Load this file. eg
(a) From the MOE | File | Open panel, select this file
and click the "Load" button.
(b) Create an environment variable MOE_SVL_LOAD that points
to a subdirectory such as $MOE/local_svl and then save
this file into that subdirectory.

2. In the Database Viewer select the Compute | Descriptors menu.
The new descriptors will be added to the list. If you type
IND_I in the "Filter" box at the bottom of the panel, then just
these descriptors will be shown.

Select the descriptors the you want to calculate and click the OK
button.

Warning: When unparameterised atoms are found, the electronegativity
will be set to -999 and the covalent atomic radius to 0. A warning will
be printed to the SVL command line.

You can choose which chain you want to select from with the pull-down. The tool then allows you to visualise specific parts of the secondary structure by adjusting the two sliders.
Residues between the sliders (red) are selected, residues outside the sliders (green) are unselected. You can skip forward or backward one residue using the single arrow buttons. The double arrows take you to the start/end of the current structural feature (e.g. back to the start of the current helix).
Once you have identified the feature you are interested in, you can (de)select the residues in the sequence editor using the buttons at the bottom. You can then show/hide residues or use the Synchronise option etc. in the normal way.
Pressing "done" will return the 2ary structure rendering to how it was when you started.

To use the code, load it into Moe and type SSS []

NEWLEAD is a computer program for the automatic
generation of candidate structures. The input for the program
is a set of fragments in the three-dimensional orientation
corresponding to a given pharmacophore model. The treatment
consists in connecting the fragments with spacers assembled
from small chemical entities (atoms, chains or ring moieties).

The results are new structures containing the fragments in
the orientation defined in the input. The program offers
the opportunity of rapidly applying a pharmacophore model
in drug-design by generating automatically a set of chemical
sructures conforming to the model.

This distribution contains several files enabling the use of
NEWLEAD within MOE.

Tschinke, V.; Cohen N.C. The NEWLEAD Program: A New Method for the Design of Candidate Structures from Pharmacophoric
Hypotheses. J. Med.Chem. 1993, 36, 3863-3870.

Binary executable files for NEWLEAD are available for Linux, OS/X and Windows at

http://www.ccl.net/

GI-MOE is a full-featured SVL-based interface for integrating
quantum calculations into MOE using Gaussian 03 as computational back-end. Calculations can be performed on single molecules or on multiple structures (in either .com or MOE .mdb format). Bypassing the cryptic Gaussian command line, various calculations can be performed in GI-MOE including: ground state and transition state calculations (3 methods), checkpoint file restarts, potential scans (with or without geometry optimization), full ONIOM implementation (2-layer and 3-layer), frequency calculations, solvation, and automated RESP charge fitting. Tools for .com file inspection and z-matrix conversion are included. All GI-MOE parameters can be saved and read from a settings file. Calculations can be performed on a local installation of Gaussian or via the GI-MOE remote computing window, which allows for automatic login and submission of GI-MOE jobs to a distributed computing cluster. A complete user manual is provided.

Performs Spectral Clustering, based upon
M.L. Brewer, J. Chem. Inf. Model, 2007 (47), 1727
Code written by M.H. Charlton at Chroma Therapeutics, 2007.

To use this code, load it into MOE and type :
specclus_gui []
Pressing Calculate performs the computation and writes the cluster information to an mdb field '$CLUSTER'.
The code uses precalculated fingerprint fields, but if none is present, it will calculate MACCS FPs for you.
The code needs the tanimoto metric to be applicable to the FP you select.

The code uses two methods of filtering the clusters, as outlined in the paper. The default "threshold" method using a value of 0.85 seems to be quite reliable, although dropping it to 0.65 is often good in the "Diverse Islands" situation (molecules in the clusters are relatively similar, but the clusters are well separated).

Note on timings : the code uses MOE's eigenvector solution function and can get quite time consuming for large sets, potentially scaling with the third power of the number of compounds.
Sample timings 2.66 GHz Xeon, MACCS FPs :
n_mol: 206 500 1000 1500 2000
t(s): 1.1 6.0 24.5 65.4 135.7
fitted line : t=1.14e-08(n^3) + 8.86e-06(n^2) + 4.80e-03(n) - 0.268

Times for larger mdbs depend on the value of the threshold used, but compare well with MOE's standard FP clustering :
~6s for 500 cpds, 25s for 1,000 cpds.

A major problems in energy minimization of the raw PDB files is the lack of hydrogen atoms in the crystallographic coordinates. Software packages including MOE can add the hydrogens in ‘default’ orientations, which lead to many steric clashes. At this point, a regular energy minimization will lead to significant perturbations and movements of hydrogens as well as the heavy atoms.

One way to decrease such perturbations is to try to ‘shift’ them towards less important areas in the protein, i.e. residues which are far from active site. The ‘PDB-Thawing’ program works on this concept.

See the online

Update 12-09-08: fixed navigator button sensitivity

The inspiration for this application is to create an interactive version of the DBV | File | Print panel.

To use this script, save the file to your hard drive then run it, e.g. using MOE | File | Open, browse to db_bigbrowser.svl and click OK.

Once the Big Browser window is open, the basic browsing functionality provided by the arrow buttons at the bottom of the panel is the same as for the standard DBV | File | Browse application:
|< moves to page 1
< moves back a page
> moves forward a page
>| moves to the last page

At the top right, the page position is shown as "n/pages". Clicking that button allows one to jump to a specific page number.

The Subject menu allows one to choose a molecule field.

Define template allows a common substructure to be entered based on atoms in the main MOE window. To use this, click Define Template, then either build the scaffold using the Builder or simply click a molecule to toggle it into MOE. Select the scaffold atoms and click Load in the "Define Template" window. Adjust the depiction (rotation or X/Y flip) and click Apply.

Set Data Fields allows the user to choose a set of data fields to be shown as text. Also a numeric field can be shown as a colour-coded square; use the Activity Indicator drop down menu to select the field.

The Action button allows the user to change what happens when a molecule is clicked in the Big Browser panel. If Toggle MOE is selected, the molecule is either created or destroyed in the main MOE window. If Toggle Sel is chosen, the entry selection state will be toggled. If Mark is chosen, the value in the text field will be written to the field specified in the drop down menu'

The Clear button clears whatever the action is currently set to, i.e. if the action is MOE, all molecules sent to the main MOE windows are deleted; if action is Select the entry selection in the database is cleared. There is no effect for Mark mode.

The Grid drop down menu sets the number of molecules depicted; the grid is defined horizontal x vertical.

The Size menu sets the size of each cell; 1 is smallest, 5 is largest.

Browse Selected sets the browser to show only the entries that were selected when the box was checked. This means that if you are browsing a set of selected molecules and you deselect one of the entries, it will not be removed from the browser. To remove it, you will need to turn off "Browse Selected" then re-enable the option.

The application is not sensitive to changes in the database (e.g. creating or deleting fields or entries); making such changes is liable to result in the application crashing. I aim to address this in a future revision. Please send any bug reports or feature requests to sgrimshaw@chemcomp.com.

Save this file to your hard disk. For it to be available for use by all
of the users of MOE, create a directory called $MOE/lib/svl/custom/run
and save this file to that directory. eg on Windows and Linux,
this might be: c:\moe\lib\svl\custom\run\extract_conformations.svl
/usr/local/moe/lib/svl/custom/run/extract_conformations.svl

For the script to be available for use by a single user, then it
is best to create a directory in your home directory called svl,
and create a directory called run in that directory. Save this
file into that directory. eg on Windows and Linux, this might be:

c:\Documents and Settings\user_name\svl\run\extract_conformations.svl
/usr/people/user_name/svl/run/extract_conformations.svl

(a) From the MOE | File | Open panel, select this file
and click the "Run" button.

(b) Alternatively, at the SVL command line, type:
run 'extract_conformations.svl

Reference:
Oyarzabal, Julen; Howe, Trevor; Alcazar, Jesus; Andres, Jose Ignacio; Alvarez, Rosa M.; Dautzenberg, Frank; Iturrino, Laura; Martinez, Sonia; Van der Linden, Ilse
Novel Approach for Chemotype Hopping Based on Annotated Databases of Chemically Feasible Fragments and Prospective Case study: New Melanin
Concentrating Hormone Antagonists
J. Med. Chem. Articles ASAP March 16, 2009

This set of tools was developed for Inhibox by Daniel Butler in conjunction with excellent UK CCG support.

The user can develop an SQL statement using their favourite (Oracle) cartridge and these scripts will generate a 3D SDF from the submitted SQL statement which is a very handy and efficient way to use the MOE tokens available within a given organisation.

The SVL modules are tied together with a script (MOE_SQL_TO_SDF.sh or MOE_SQL_TO_SDF.bat) which prompts the user for the relevant information and a SMILES column and primary key column must be specified. If the primary key is of type number then this must be converted to a varchar using oracle to_char function (if not oracle then equivalent). If running on linux you may need to run dos2unix on the files for them to work.

Enjoy + I hope it makes a difference!

Daniel Butler.
Jan 2009

POV-Ray is a freeware ray tracing program. See http://www.povray.org for more details. This script writes a POV-Ray scene from MOE. It creates spheres for atoms, cylinders for bonds and triangles for surfaces for MOE graphical objects. The advantage of using a ray tracing program like POV-Ray is that you get much more control over the appearance of the image, with the option to use settings which would be too slow for use in an interactive program. You can also output the images in any size you like (eg for publications, posters or an exhibition stand), with no loss of crispness.

At the moment it does not draw anything for a backbone representation. The trace representation is easy (cylinders between the CA atoms). For the cartoon and tube representations, you can use $MOE/sample/cartoon.svl This creates a graphical object that is then rendered the same as a surface.

The script now handles graphic objects with a dotted, line and solid surfaces.

The position of the camera is not exactly correct, so the perspective view looks a little different to the MOE window. Any suggestions would be gratefully received.

The positions of the lights in the POV-Ray scene are also static. Edit the coordinates and add extra lights as you see fit. You can also modify the textures of the objects with POV-Ray commands (eg to make a surface shiny, smooth, bumpy, etc).

TO DO

Clean up code (error checking, comments, split into functions)

Load positions of camera, lights and scale

Cartoon, Tube and C-Alpha Trace Backbone

Double/triple bond sticks (including SO2, NO2)

Pick up display settings from MOE

Textures

Stereo views

Display labels and meters

USAGE

1. Load a view in MOE which has the surfaces, atoms and bonds displayed as you want in the output. If you want backbone representations, then use $MOE/sample/cartoon.svl.

2. Load this script using one of these methods:

(a) From the File | Open panel, click on this with your left mouse button and click the "load" button to the left hand side of the panel

(b) At the SVL command line type
load 'write_pov.svl'

(c) Set an environment variable MOE_SVL_LOAD that points to a subdirectory like $MOE/local_svl and save this file into that subdirectory. It will then be loaded the next time that MOE starts.

3. At the SVL command line type: Write_POV 'out.pov'

4. Load the out.pov file into POV-Ray and render it.

DESCRIPTION:
This plugin saves an MOE screen picture several times per second into .bmp image files, while you manipulate a molecular system in MOE. It is a MOE video recorder, rather than a 'macro recorder'.
You can then create a video file from these movie frames.

USAGE:
1. MOE | File | Open | "Run SVL"
2. Once you have finished recording, use a third-party program to join the .bmp image frames into a video file. One such program (there are others as well) available at the time of writing, is bmp2avi
http://www.divx-digest.com/software/bmp2avi.html

D.Bondarev

Calculates pairwise rmsd in a database between mols in first two adjacent mol fields or to fields
indicated by the fields parameter. WE ASSUME these are pairs of conformers arbitrarily positioned in space and that these have atoms in the same order.

Usage:
1) Save and Load the file pairwise_rmsd_superpose.svl
2) Type the following at the CLI

svl> pairwise_rmsd_superpose []

This will work on any open database containing two fields of molecules in different conformations. Te explicitly set the database and fields to operate on, type:

svl> pairwise_rmsd_superpose ['myDB.mdb', ['mol_field1', 'mol_field2']]

where 'myDB.mdb' is the database to work on and 'mol_field1' and 'molfield2' are fields in the database containing molecules to operate on.

Output

Overlap

=======

With in-house With it-self

============= ============

25.4% 17.1%



Number of Clusters

==================

Total number of entries: 101040.0

Total number of inhouse entries: 125170.0

Total number of outside entries: 9600.0



Number of cluster from outside libraries
having more than one compound 1034.000
Number of entries that overlap with reference database 2441.0
Number of entries that overlap with itself 1639.0

Number of cluster from outside libraries
having one compound 7043.000



Number of cluster from inhouse libraries
having more than one compound 16156.000
Number of cluster from inhouse libraries
having one compound 78679.000

====================

outside freq_outside

0 92963
1 7043
2 750
3 193
4 44
5 23
6 10
7 7
9 2
10 3
14 1
16 1


====================

reference freq_inhouse

0 6205
1 78679
2 10682
3 2898
4 1059
5 548
6 320
7 151
8 113
9 91
10 59
11 44
12 33
13 29
14 21
15 15
16 17
17 8
18 5
19 6
20 6
21 4
22 6
23 4
24 4
25 2
26 2
27 3
28 1
29 2
30 1
31 1
32 2
34 4
35 1
37 1
38 1
39 1
42 1
43 2
48 2
50 1
56 1
58 1
60 1
61 1
67 1

guides the user through a few simple steps to perform a FlexX computation:

- check for a valid FlexX installation
- load a protein
- protonate properly
(optionally using the new protonate3D)
- launch the FlexX-in-MOE interface
with pre-configured default settings
- determine a reference ligand
for optional RMSD calculation

The results will be ready for post processing
in a MOE mdb database.

db_COMBINE calculates the components of ligand-residue interaction energies as descriptors like COMBINE descriptors (Ortiz, et al., 1995). The interaction energies per residue are divided into van der Waals and electrostatic contributions as described in Arakawa et al. (2008). The Born solvation contributions are also calculated with the constant or distance-independent dielectric constant with the value of 1 if the Generalized Born implicit solvation model is selected as the solvation model. The changes in the internal conformation energies of the ligand and the receptor included in the original COMBINE descriptors are not calculated. Cutoff for the non-bonded interaction energies is set to off to ensure calculating non-distorted energies contributed from distant receptor atoms. Receptor residues are defined as the residues with selected atoms in MOE window if there are selected atoms, or the all residues if no atoms are selected.

Installation:
1. Load db_combine.svl in MOE. Or save it in $MOE/lib/svl/rsi.svl/ or $MOE/lib/svl/patch directory to automatically load it when MOE is invoked.

Usage:
1. Click Window | Potential Setup and select Born as Solvation in the Potential Setup window if you want to include the contribution from the solvation energy of Generalized Born implicit solvation model.
2. Read a receptor into MOE.
3. Select atoms to limit the residues to be calculated. If no atoms are selected, all atoms in MOE window will be taken into account.
4. Open a ligand MDB in DBV.
5. DBV | CLI | db_ COMBINE []

The fields for contributions of van der Waals, electrostatic and, if selected, solvation interaction energies per residue will appear in MDB.

References:
[1] A.R. Ortiz, M.T. Pisabarro, F. Gago, R.C. Wade, Prediction of Drug Binding Affinities by Comparative Binding Energy Analysis, J. Med. Chem. 1995 38, 26812691.
[2] M. Arakawa, K. Hasegawa and K. Funatsu, Tailored scoring function of Trypsin-benzamidine complex using COMBINE descriptors and support vector regression, Chemom. Intell. Lab. Syst. 2008 92, 145-151.

PDF tutorial, step by step, illustrated.

Example files on Way244

SVL utility for template-based superposition and alignment.

PostDock is a new visualization tool for the analysis and comparison of molecular docking results. It processes a docking database and displays an interactive pseudo-3D snapshot of multiple ligand docking poses such that their docking energies and docking poses are visually encoded for rapid visual assessment. The docking energies are represented by a transparency scale whereas the docking poses are encoded by a color scale. The submission includes postdock.svl, a user manual, as well as a reprint of the publication.

The current version is using a revised rule set. In rule 7 the aromaticity definition is used in a more strict manner than in previous versions. Rule 8 and 9 now include heteroatoms connected by double bonds to carbocycles to distinguish between carbo cycles.
The resulting clustering should be improved compared to the original rule set.

Implemenation based on:

Schuffenhauer, Ansgar; Ertl, Peter; Roggo, Silvio; Wetzel, Stefan; Koch, Marcus A.; Waldmann, Herbert
The Scaffold Tree - Visualization of the Scaffold
Universe by Hierarchical Scaffold Classification
J. Chem. Inf. Model. (47) 2007, 47-58

Calculates a hierarchical scaffold tree analysis for a database of compounds. Writes out statistics for frequencies of scaffolds in
actives and inactives.
For cyclic compounds the algorithm is based on the paper mentioned above. For acyclic compounds I have used my own rules.
Rules can be overwritten by a priorityscaffold. This Scaffold will be kept together until no other ring is present.

USAGE: run 'ScaffoldTree.svl' to start the analysis or display previously calculated results

Or load this file and use:
db_ScaffoldTree [database, opt]
where opt is a tagged vector:
[molecules: molfield, activity: activityfield,
threshold: activitythreshold,
operator: activityoperator as token,
prioscaff: SMARTS token]
runs the analysis in batch. operator may be '<' or '>'
if no activity field is specified use 'none' or []
prioscaff is the SMARTS for a prioritzed scaffold
(leave ' if none) only molecules is a required tag!

To extend an existing scaffold tree:
db_ResumeScaffoldTree [database, opt]
opt is defined as for db_ScaffoldTree but has the additional tag 'SCT_Output' for the scaffold database.

See pdf file for more details.

Description:
-similar to the conformational superpose, except that
selected atoms determine a molecular subfragment that is
used to perform the superpose. This is useful when superposing
a set of ligands on to a common substructure

Usage:
(1) First, load a reference ligand into the 3D MOE
window. The ligand should contain the common core
and be in the correct absolute coordinates (if appropriate).
The coordinates of this structure will be used to superpose
the database ligands
(2) Run this function -a GUI will appear. Use the browse button
to select the database of ligands to align.
(3) The alignment may be performed based one of two
criteria for specifying superposition point correspondence:
(i) A SMILES/SMARTS match of the conserved core
(ii) Atoms with unique names throughout the database
(4) (i) Superposing with a SMILES/SMARTS string.
To superpose on a SMILES/SMARTS string of the conserved
core, select the atoms of the consered core and
press the 'Extract' button. A SMILES string of the
selected atoms will be produced. This can be edited to
include wildcards, etc. Then press
'Align by SMARTS'
(ii) Superposing by unique atom name
When the database file is specified, the program will scan
through the first molecule field in the database and
automatically detect 'unique atom names' i.e. (C1, N2),
which are names that occur with only one instance in
any molecule, and occur in every molecule.
The names will be listed in the Alignment Points box.
To superpose on these named atoms, load a template
3D reference structure which has the appropriate
atom names. Select the named atoms in the MOE window
AND in the Alignment Points list. Then press
'Align by Name'

(5) A new field named transformed_mol_field will appear in the
database, with the transformed structures. A comment field
will also be generated, indicating if a the point mapping
and transformation were successful.

This is an interface to the InChI molecule identifier calculator available at:
http://www.iupac.org/inchi/download/InChI-1.zip

When the program is started for the first time you have to browse to the command line version of the InChI tool in the 'Settings & Options' panel to set the path to the executable.

REFERENCE:
Gregori-Puigjane, Elisabet; Mestres, Jordi
SHED: Shannon Entropy Descriptors from Topological Feature Distributions
J. Chem. Inf. Model. 2006, 46, 1615-1622

USAGE: load 'fp_shed.svl'
FP:SHED is added to fingerprints and to the database browser. Molecules in the main window are drawn in black. The current entry is colored in red.
To create a shedplot of all molecules in the main window use SHED_Plot[]. Chain color is used for the plots of different molecules.

Calculates CATS vectors, based upon Schneider et al, Angew Chemie, 38, 2894-2896 with augmentation to included aromatic atom types based upon unpublished work by M.H Charlton, M.L. Brewer and P.N. Mortenson carried out at Evotec.
This code written by M.H. Charlton at Chroma Therapeutics, 2007.

To use the code, load it into MOE and type :
db_cats []
Alternatively, to use a specific molecule field (e.g. mol_wash), type :
db_cats ['mol_wash']
Published version normalises the CATS vectors by the number of heavy atoms. To achieve this, use :
db_cats ['mol_wash', 1]
Doing this will produce 210 fields of type float. Not setting the second argument produces a simple count of the various CATS features (210 * integer).
Failing to provide a molecule field will default to the first such field in the database.
Code uses MOE's definitions of various Ph4 types (e.g. anion, cation) and is therefore best used after washing the molecules to ensure these features are recognised.

Usage>
1. Load 'lightdirection.svl'
2. Type the "LightDirection[]" instruction in CLI

Now light directions can be setted by [MOE|Render|Setup...]. But because we need to input the values of the vector, it is little difficult.

By using this program, you can adjust them smoothly.

Usage:
(1) Open the template structure that is in the
correct reference position. Select the atoms of the
common sub-structure against which you wish to use
to calculate for the RMS deviation.
(2) Open the database of structures you wish to align
(3) run fragment_rmsd.svl
(4) Alternatively, load fragment_rmsd.svl and enter
svl>database_molecule_rmsd ['database.mdb']
(5) A new field 'RMS_Deviation#' will be produced where #
is an integer that increments to generate a unique field name