Macs in Chemistry

Insanely Great Science

computational chemistry

Mathematica version 12 released


I'm sure many will have noticed the release of Mathematica version 12, the website contains details of the many additional new features, but I happened to notice there were a couple of Chemistry features including:

The Molecule is a symbolic representation of a chemical species and is a fully computable first-class member of the Wolfram Language. More details here.


Open Forcefield 0.2


The 0.2.0 release of the Open Force Field Toolkit, featuring RDKit support and the new-and-improved SMIRNOFF v0.2 force field spec has been announced.

We're excited to announce the public release of the Open Force Field toolkit version 0.2.0! Most notably, this release adds the ability to assign SMIRNOFF parameters and AM1-BCC charges with a completely open-source backend, adding support for the RDKit and AmberTools via a new ToolkitWrapper infrastructure that can be extended in the future to support additional cheminformatics toolkits. The OpenEye Toolkit will continue to be supported, as well as used internally our parameter-fitting pipelines in the short term. We're extremely grateful to the long list of contributors that have made this release possible, especially Shuzhe Wang from the Riniker group for piloting much of the RDKit functionality.


QUBEKit: QUantum BEspoke FF toolKit


Just saw an interesting paper "QUBEKit: Automating the Derivation of Force Field Parameters from Quantum Mechanics" DOI.

QUBEKit is python based force field derivation toolkit that allows users to derive accurate molecular mechanics parameters directly from quantum mechanical calculations.

Code is available on GitHub QUBEKit, and there is a user tutorial on the Wiki Page.


  • Anaconda3
  • Biochemical and Organic Simulation System (BOSS)
  • OpenMM
  • Gaussian09
  • Matlab 2017

Python modules used:

  • numpy
  • argparse
  • collections
  • colorama
  • matplotlib


Py-ChemShell 2019 released


Just got this notice.

We are pleased to announce Py-ChemShell 2019 (v19.0), the first beta release of the Python-based version of ChemShell. Py-ChemShell 2019 offers new interfaces to ORCA and DL_POLY 4, a complete task-farmed parallelisation framework (including parallel finite difference gradients), RESP charge fitting procedures, and case studies for problems in materials modelling.

Py-ChemShell can be downloaded free of charge under the open source GNU LGPL v3 licence from this site.

ChemShell is a scriptable computational chemistry environment for multiscale modelling. While it supports standard quantum chemical or force field calculations, its main strength lies in hybrid QM/MM calculations. The concept is to leave the time-consuming energy evaluation to external specialised codes, while ChemShell takes over the communication and data handling.

ChemShell provides interfaces to a variety of QM and MM codes, including:

  • NWChem
  • FHI-aims
  • MNDO
  • Orca
  • Molpro
  • Gaussian
  • DMol3
  • Q-Chem
  • GULP


MOE 2019.1 released


The 2019.01 release of Chemical Computing Group's Molecular Operating Environment software includes a variety of new features, enhancements

Full release notes are here.


  • Calculate and Analyze pH-Dependent Protein Properties
  • MOEsaic Session Sharing and Project Customization
  • Determine Conformation Population from NMR NOE Data
  • Predict Relative Binding Energies with AMBER Thermodynamic Integration

Worth noting there is an updated Version of Flexera License Manager.

MOE now uses an updated version of the Flexera license manager. The license manager server components lmgrd, chemcompd, and lmutil have all been updated to version Note that older versions of MOE will continue to run with updated license manager servers.


LICHEM: Layered Interacting CHEmical Models


An update to LICHEM: Layered Interacting CHEmical Models has been published DOI

LICHEM is an open-source (GPLv3) interface between QM and MM software so that QM/MM calculations can be performed with polarizable and frozen electron density force fields. Functionality is also present for standard point-charge based force fields, pure MM, and pure QM calculations.

Available from GitHub

Note, On OSX machines, the SEDI, TEX, BIB, and CXXFLAGS variables will need to be modified.


ORCA Version 4.1 released


With the release of ORCA 4.1, they have moved our forum and download site to a new server at the Max Planck Institute fuer Kohlenforschung, where the ORCA team now has its home base. Now at

ORCAis an ab initio quantum chemistry program package that contains modern electronic structure methods including density functional theory, many-body perturbation, coupled cluster, multireference methods, and semi-empirical quantum chemistry methods. Its main field of application is larger molecules, transition metal complexes, and their spectroscopic properties. DOI.

List of new features for ORCA 4.1:


  • B97M-V, wB97M-V, wB97X-V plus various D3 variants of B97 functionals
  • Simple input keywords for DSD-BLYP, DSD-PBEP86, and DSD-PBEB95
  • CPCM analytic Hessian
  • DLPNO-double hybrid DFT including gradient
  • SymRelax option in %method

Semiempirical methods

  • XTB method of Grimme et al.

Coupled cluster

  • Iterative solution of the full (T) equations for DLPNO-CCSD(T)
  • Open shell DLPNO-CCSD density and spin density matrices
  • Full DLPNO-MP2 gradient
  • CIM (Cluster in molecules) Implementation with MP2, CCSD(T), DLPNO-MP2 and DLPNO-CCSD(T)
  • IP and EA coupled cluster methods and their DLPNO variants
  • STEOM-CCSD for open shells
  • SOC between bt-PNO-STEOM and STEOM states
  • Improved Multilevel implementation including multilevel DLPNO-IP
  • F12-Triples scaling for RHF canonical CCSD(T) based on the CCSD/ CCSD-F12 ratio


  • New CASSCF SuperCIPT converger is reliable and efficient.
  • New options for final orbitals to find partner orbitals for the chosen active space e.g. bonding / anti-bonding partners.
  • MC-RPA (Multiconfigurational random phase approximation)
  • ◦ AO driven integral direct for calculations on larger molecules
  • ◦ Fock matrix -> conventional, direct, RIJ/COSX
  • ◦ MPI parallel
  • ◦ NTOs for visualizing transitions
  • Checking stability of state specific CASSCF wave functions by orca_mcrpa
  • Dynamic correlation dressed (DCD-CAS) method with inclusion of relativistic effects (SOC, spin-spin, magnetic fields)
  • CASSCF RIJCOSX allows two separate auxiliary basis sets
  • CASCI/NEVPT2 protocol for XAS and RIXS


  • Nudge elastic band method to locate transition states
  • Enabled 3-dimensional relaxed potential energy surface scan
  • Improvement of redundandant internal coordinate generation
  • Faster and more smooth convergence for 3-dimensional systems and embedded cluster models
  • Intrinsic reaction coordinate (IRC) following
  • Swart model Hessian (good for weak interactions)

Molecular Dynamics

  • MD simulations can now use Cartesian, distance, angle, and dihedral angle constraints.
  • The MD module now features cells of several geometries (cube, orthorhombic, parallelepiped, sphere, ellipsoid), which can help to keep the system inside of a well-defined volume.
  • The cells can be defined as elastic, such that their size adapts to the system. This enables to run simulations under constant pressure.
  • Ability to define regions (subsets of atoms) enables applications such as thermostating different parts of the system to different temperatures (cold solute in hot solvent, temperature gradients, ...)
  • Trajectories can now be written in XYZ and PDF file format.
  • A restart file is written in every simulation step. Simulations can be restarted to seamlessly continue.
  • The energy drift of the simulation is now displayed in every step.
  • The MD module now works with a broader range of methods (semiempirics, ECPs, QM/MM).
  • Fixed a bug in the time integration of the equations of motion which compromised energy conservation.

Spectroscopic properties

  • orca_pnmr module tool to calculate paramagnetic NMR spectra
  • NMR chemical shifts with RI-MP2 and double hybrid DFT including GIAO’s, spin-component scaling and CPCM
  • NMR Spin-Spin coupling in calculations with DFT/HF
  • NMR wth ZORA
  • Maximoff-Scuseria correction for the kinetic energy density in GIAO-based calculations with meta-GGA functionals
  • Exact and gauge invariant transition moments and approximate decomposition into dipole, quadrupole etc terms in all modules.
  • PNO-ROCIS method for more efficient X-ray absorption calculations
  • IP-ROCISD for high spin ROHF references
  • TD-DFT:
  • Transient spectra (excited state absorption) for CIS/TDA
  • Triplet gradients (with RIJ, COSX and all) for all cases.
  • Spin orbit coupling (including CPCM) and gradients
  • Root following scheme for optimization
  • Slow term to correct energy of relaxed excited state
  • Full TD-DFT with double hybrids
  • ESD module to calculate spectroscopic properties
  • Vibrationally resolved absorption spectra including Duschinsky rotation and/or vibronic coupling.
  • Fluorescence and Phosphorescence rates with same options.
  • Resonance Raman spectra with the same options
  • works with CIS/TDDFT, ROCIS, CASSCF and EOM/STEOM.
  • Seven different schemes for obtaining an excited state PES and five different choices of coordinate systems

Analysis tools:

  • Open Shell LED
  • Dispersion interaction Density plots
  • LED for DLPNO-MP2
  • LED for the frozen state
  • Update of AIM interface
  • NBO 7 compatibility (i4)
  • Miscellaneous
  • Compound method (Infrastructure, plus W2.2, W1, G2(MP2), G2(MP2-SVP), G2(MP2-SV) methods)
  • Property file (additional properties, plus new infrastructure)
  • Decomposition of correlation energy for canonical RHF CCSD energies to singlet - triple pairs
  • Additional EP2 extrapolation schemes using RI-MP2 and DLPNO-MP2 methods as cheap methods (request from forum)
  • Lanthanide new def2 basis sets
  • def2-XVP/C auxiliary basis sets for Ce-Lu by Chmela and Harding.
  • Robust Second order optimizer for localized orbitals
  • Added a few basis sets.


DIRAC18 released


The DIRAC program computes molecular properties using relativistic quantum chemical methods. It is named after P.A.M. Dirac, the father of relativistic electronic structure theory.

I can be downloaded from the zenodo repository.

New features are described here.

DIRAC, a relativistic ab initio electronic structure program, Release DIRAC18 (2018), written by T. Saue, L. Visscher, H. J. Aa. Jensen, and R. Bast, with contributions from V. Bakken, K. G. Dyall, S. Dubillard, U. Ekström, E. Eliav, T. Enevoldsen, E. Faßhauer, T. Fleig, O. Fossgaard, A. S. P. Gomes, E. D. Hedegård, T. Helgaker, J. Henriksson, M. Iliaš, Ch. R. Jacob, S. Knecht, S. Komorovský, O. Kullie, J. K. Lærdahl, C. V. Larsen, Y. S. Lee, H. S. Nataraj, M. K. Nayak, P. Norman, G. Olejniczak, J. Olsen, J. M. H. Olsen, Y. C. Park, J. K. Pedersen, M. Pernpointner, R. Di Remigio, K. Ruud, P. Sałek, B. Schimmelpfennig, A. Shee, J. Sikkema, A. J. Thorvaldsen, J. Thyssen, J. van Stralen, S. Villaume, O. Visser, T. Winther, and S. Yamamoto (available at, see also


Installing Osprey 3.0 under Mac OS X


A recent publication described OSPREY 3.0: Open-Source Protein Redesign for You, with Powerful New Feature DOI.

We present Osprey 3.0, a new and greatly improved release of the osprey protein design software. Osprey 3.0 features a convenient new Python interface, which greatly improves its ease of use. It is over two orders of magnitude faster than previous versions of osprey when running the same algorithms on the same hardware. Moreover, osprey 3.0 includes several new algorithms, which introduce substantial speedups as well as improved biophysical modeling. It also includes GPU support, which provides an additional speedup of over an order of magnitude. Like previous versions of osprey, osprey 3.0 offers a unique package of advantages over other design software, including provable design algorithms that account for continuous flexibility during design and model conformational entropy. Finally, we show here empirically that osprey 3.0 accurately predicts the effect of mutations on protein–protein binding.

Osprey 3.0 is available at as free and open‐source software GPLv2.

The source code is available on GitHub

Unfortunately the installation instructions do not include Mac OSX but there are instructions for "Debian-like Linux" which seemed promising. With the invaluable help of Nathan Guerin I was able to get OSPREY installed.



ConstruQt API


Just got details of an interesting service

ChemAlive ( would like to offer ConstruQt, its core molecular design tool based on quantum mechanics (QM), for trial.

Currently you can:

  • Transforms list of SMILES or InChI molecular designations into state-of-the-art 3D molecular structures in SD format
  • Manages the conformational space of the molecules with a robust shape searching algorithm
  • Generates all reasonable tautomeric forms of the molecule and prioritizes them by energy
  • Generates all diastereomeric forms of the molecules and differentiates them by energy
  • All molecules are stored in our unique database architecture making the calculations easily augmented and carried through to other processes

The last bullet point is worth noting, so don't submit anything confidential.


IGMPlot release 2.4


The new IGMPlot release 2.4, is available for download at . It provides chemists with a visual analysis of covalent and non-covalent interactions

Detailed installation notes are in the documentation (page 5).

IGMPlot is written in C++. It has been installed and tested on several platforms: computational centers (linux), MacOS, Windows10, and several compilers and versions (GNU, Intel, PGI), it can be compiled with or without OpenMP support

On MacOs machines, a sequential version of IGMPlot can be obtained with the Clang compiler. In the Makefile choose the options:

  • CppCompilerFamily=GNU
  • CppCompilerVersion=5andabove o OpenMP=NO
  • CC=g++

On MacOS machines, to leverage OpenMP multicore execution, you must install a gcc (g++) version different from the one provided within the compiler front end “Clang” which until now has not built-in support for OpenMP. You might install gcc with the command: ‘brew install gcc -- without-multilib’ (see for instance openmp-on-mac-os-x-10-11). This way, the compiler might be installed somewhere like /usr/local/Cellar/gcc/7.1.0/bin/g++-7. In this example, make sure the g++-7 command be available with your PATH and adjust the IGMPlot makefile accordingly (changing the g++ command with g++-7 for instance).

This link might also be useful OpenMP under MacOSX.



The SAMPL6 Blind Prediction Challenge for Computational Chemistry


Now on GitHub

SAMPL6 Part II will include a octanol-water log P prediction challenge and will be followed by a joint D3R/SAMPL workshop in San Diego, Aug 22-23, 2019, immediately before the San Diego ACS National Meeting. A special issue or special section of JCAMD will be organized to disseminate the results of this challenge.


Open Force Field Consortium


The Open Force Field Consortium, an academic-industry collaboration designed to improve small molecule force fields used to guide pharmaceutical drug discovery.

The Consortium will develop an extensible, open source toolkit for constructing, applying, and evaluating force fields; produce and curate public datasets necessary to build high-accuracy biomolecular force fields; and apply these tools and datasets to generate improved force fields. Academic and industry partners work together to ensure its success.


Turbomole Update


There is a new Turbomole release

TURBOMOLE has been developed to provide a fast and stable code to treat molecules for industrial application. With the TUBROMOLE implementation of RI-DFT, one of the fastest DFT methods will be available at your fingertips.

TURBOMOLE V7.3 (July 2018) New features:

  • PNO-CCSD(T0) and PNO-CCSD(T) energies for closed-shell systems [1]
  • new DFT-D4 dispersion correction based on xTB [2]
  • modernized NMR (with RI-J, COSMO, meta-GGAs, low-order scaling HF-exchange, SMP parallelization) [3]
  • VCD spectra using COSMO
  • periodic DFT with larger basis sets (treatment of linear dependency)
  • two-photon absorption cross sections and analytic frequency-dependent hyperpolarizabilities with TDDFT/TDHF [4]
  • X2C gradients for 1- and 2-component DFT, full X2C and DLU-X2C [5]
  • vibronic absorption/emission spectra (new module: radless) [6]
  • CC2 vertical excited states with COSMO [7]
  • NTO (natural transition orbitals) for TDDFT
  • RI-GW based on dRPA (very fast GW and BSE) [8]


  • GW and Bethe-Salpeter based on fast dRPA
  • support of RI-J and linear scaling HF exchange in NMR calculations
  • PNO-MP2 closed shell energy calculations significantly more efficient


  • new scripts for parallel execution which recognize the most frequently used queuing systems

TmoleX (4.4) now supports:

  • DFT-D4 dispersion correction
  • X2C relativistic two-component treatment for spin-orbit coupling terms, and new X2C basis sets
  • Fukui indices and functions (calculation and visualization)
  • movie exports to mp4 file format
  • B97-3c functional


OMEGA v3.0.1 released


OpenEye have just announced the release of OMEGA v3.0.1 This upgrade fixes several bugs and adds a number of internal improvements.

Major bug fixes

  • A bug that caused memory leaks in OMEGA classic, dense, pose, and rocs modes, has been fixed. Previously, a substantial memory leak was experienced when running OMEGA on a large database.
  • OMEGA macrocycle no longer uses excessive memory for molecules with terminal heavy atoms.

OMEGA performs rapid conformational expansion of drug-like molecules, yielding a throughput of tens of thousands of compounds per day per processor. OMEGA is very effective at reproducing bioactive conformations, and provides an optimal balance between speed and performance when used on large compound databases.


ACS awards for Computers in Chemistry


Nominations are now open for the Computers in Chemistry division of the ACS awards.

More details here


Schrödinger Software Release 2018-2


Schrödinger have announced a major update their software suite.

Full details are here


BiosolveIT update SeeSAR and more.


BioSolveIT have announced significant changes and improvements in SeeSAR resulting in another major release to version 8. The biggest change is that they now provide full protein visualization support. While the focus of the tool is for the most part still on the defined binding site, you can now...: see the whole protein in all its glory! As always, a major update means that HYDE scores must be re-calculated to stay in line with the changes made in the underlying structures. We certainly believe that these enhancements are well worth it:

  • improved alignment
  • full protein support in the seqence view
  • search&find specific amino acids, waters or other protein components
  • all protein visualization controls bundled
  • enhanced pharmacophore handling
  • fragment growing for covalent binders

For details see:

They also have two new tools:

REALSpaceNavigator is the world's largest, ultra-fast searchable chemical space developed in collaboration with Enamine Ltd. It comprises roughly 3.8 billion compounds today, which will be delivered on demand in less than 4 weeks with an exceptional success rate of 80% and above.

PepSee is a software tool for interactive, visual compound prioritization as well as the design of next-generation peptide therapeutics. Peptide design ideally supports a multi-parameter optimization to maximize the likelihood of success. PepSee visualizes the relevant parameters at hand, side by side with the sequence data. Color-coded display stimulates SAR exploration. The main features of PepSee comprise:

  • comfortable sequence & data import (from Excel, FASTA, PLN, Text, even PDF)
  • automated as well as manual sequence alignment
  • various data coloring and plotting options
  • organizing and annotating your compounds
  • interactive design of novel peptides

Chemical Information and Computer Applications Group (CICAG) website


The new RSC CICAG website is now live why not have a look and provide suggestions and feedback.


The Chemical Information and Computer Applications Group (CICAG) is one of the RSC’s many member-led Interest Groups, which exist to benefit RSC members and the wider chemical science community.

Also provides links to the social media feeds (Twitter, LinkedIn etc.)


deMon2k code version 5 released


deMon (density of Montréal) is a software package for density functional theory (DFT) calculations. It uses the linear combination of Gaussian-type orbital (LCGTO) approach for the self-consistent solution of the Kohn-Sham (KS) DFT equations. The calculation of the four-center electron repulsion integrals is avoided by introducing an auxiliary function basis for the variational fitting of the Coulomb potential.

The user guide provides installation instructions and requires a Fortran compiler, BASH and MPI.


Amber 18 and AmberTools 18released


Amber is a suite of biomolecular simulation programs. It began in the late 1970's, and is maintained by an active develpment community

Amber 18 ajor new features include:

  • Free energy calculations on GPUs
  • GPU support for 12-6-4 ion potentials
  • Domain decomposition for CPU-parallelism
  • Nudged elastic band calculations for pmemd (CPU and partial GPU implementation)
  • Constant redox potential calculations, to supplement constant pH simulations
  • Support and significant performance improvements for the latest Maxwell, Pascal and Volta GPUs from NVIDIA.
  • New pmemd.gem code for advanced force fields, including AMOEB

AmberTools 18 new features include

  • CUDA-enabled pbsa solver; extensions for membrane modeling with PB *lambda-dynamics method for constant pH simulations *packmol_memgen tool for building lipids and bilayers *New ("middle") integration algorithms in sander *Build tools based on CMake *Continued updates and extensions to cpptraj: *ability to obtain energies from snapshots of PME simulations *Pairlist and other speedups *improved scripting abilities

Instructions for installing Amber under Mac OSX are here

You will need to install gfortran, whilst you can download the binary it might be worth considering using Homebrew as described here


NWChem updated


Just catching up.

NWChem 6.8 is now available on Github

NWChem provides many methods for computing the properties of molecular and periodic systems using standard quantum mechanical descriptions of the electronic wavefunction or density. Its classical molecular dynamics capabilities provide for the simulation of macromolecules and solutions, including the computation of free energies using a variety of force fields. These approaches may be combined to perform mixed quantum-mechanics and molecular-mechanics simulations.

Instructions for compiling NWChem on various platforms including Mac OSX


SAMSON 0.7.0 is available


SAMSON has been updated with a number of cool features, I particularly like the embedded Jupyter console.

SAMSON is a platform for computational nanoscience.

Python scripting is now available! Most of the SAMSON API is exposed in Python, and a Jupyter console embedded in SAMSON allows you to create models and run simulations, generate movies, perform analysis and reporting, etc., directly from scripts.


What’s more, Python makes it even easier to integrate and pipeline SAMSON and SAMSON Elements with well-known packages from diverse fields, e.g. TensorFlow, PyRosetta, RDKit, ASE, etc., to name a few.


MOE update 2018.01 released


The latest update to Chemical Computing Group's Molecular Operating Environment (MOE) software includes a variety of new features, enhancements

Windows XP (finally!) and macOS 10.6 have been removed from the list of officially supported platforms. Supported Windows platforms are Vista/7/8/10, and the minimum supported macOS is 10.7 (Lion).

Amber14:EHT Forcefield. The Amber14 parameter set is now supported in MOE. The new parameters consist of improvements to nucleic acids; otherwise, protein and small molecule parameters (and charges) are unchanged. The forcefield can be selected in the MOE | Footer.

TCR-MHC Protein Complex Database. A new MOE Project database containing T-Cell Receptor (TCR) – Major Histocompatibility Complex (MHC) x-ray structures has been added to MOE. The database can be accessed with MOE | Protein | Search | TCR-MHC | TCR-MHC which will launch the MOE Project Search panel.

Several applications have been parallelized to run in the moe -mpu environment:

  • Descriptor calculations with the SVL function QuaSAR_DescriptorMDB.
  • Energy minimization in the Database Viewer DBV | Compute | Molecule | Energy Minimize.
  • Conformational search using MDB input files in MOE | Compute | Conformations | Search.
  • Rotamer library generation with DBV | Compute | Build Rotamer Library.
  • Project database creation with the SVL run file dbupdate.svl and the scripts $MOE/bin/projupdate and $MOE/bin/projupdate.bat.

I plan to review the latest version of MOE in the near future.


Predicting the Conformational Energy of Small Molecules


An interesting publication in JCIM, Atom Types Independent Molecular Mechanics Method for Predicting the Conformational Energy of Small Molecules, DOI.

We report herein our effort to incorporate lone pairs into our model to extend its applicability domain to any saturated small molecules. The developed model H-TEQ 2 has been validated on a wide variety of molecules from polyaromatic molecules to carbohydrates and molecules with high heteroatoms/carbon ratios.




YANK is a GPU-accelerated Python framework for exploring algorithms for alchemical free energy calculations.


  • Modular Python framework to facilitate development and testing of new algorithms
  • GPU-accelerated via the OpenMM toolkit
  • Alchemical free energy calculations in both explicit and implicit solvent
  • Hamiltonian exchange among alchemical intermediates with Gibbs sampling framework
  • General Markov chain Monte Carlo framework for exploring enhanced sampling methods
  • Built-in equilibration detection and convergence diagnostics
  • Support for AMBER prmtop/inpcrd files
  • Support for absolute binding free energy calculations
  • Support for transfer free energies (such as hydration or partition free energies)

Install using conda

$ conda config --add channels omnia --add channels conda-forge
$ conda install yank

conda will install dependencies from binary packages automatically, including difficult-to-install packages such as OpenMM, numpy, and scipy. YANK runs on Python 3.5, and Python 3.6


Python and compchem


Python seems to becoming the lingua franca for scientific scripting/progamming and it is perhaps not surprising that we now see increasing support for computational chemistry.

Chemtools is a set of modules that is intended to help with more advanced computations using common electronic structure methods/ programs. Currently the is some limited support for Gamess-US and MolPro program packages but other codes can be easily interfaced. It requires:

  • Python works with Python 2.7.x and 3.x
  • numba
  • numpy
  • mendeleev
  • scipy
  • setuptools

Chemtools is NOT hosted on pypi yet but in can be installed by pip from the bibbucket repository with:

pip install

Pygamess is a GAMESS wrapper for Python, it requires:

  • Python 2.6 or later (not support 3.x)
  • RDKit

It can be installed using pip

pip install pygamess


single point calculation with RDKit

from pygamess import Gamess
from rdkit import Chem
from rdkit.Chem import AllChem
m = Chem.MolFromSmiles("CC")
m = Chem.AddHs(m)
g = Gamess()
nm =

PyQuante: Python Quantum Chemistry, an open-source suite of programs for developing quantum chemistry methods, it currently supports

  • Hartree-Fock: Restricted closed-shell HF and unrestricted open-shell HF;
  • DFT: LDA (SVWN, Xalpha) and GGA (BLYP) functionals;
  • Optimized-effective potential DFT;
  • Two electron integrals computed using Huzinaga, Rys, or Head-Gordon/Pople techniques; C and Python interfaces to all of these programs;
  • MINDO/3 semiempirical energies and forces;
  • CI-Singles excited states;
  • DIIS convergence acceleration;
  • Second-order Moller-Plesset (MP2) perturbation theory.

cclib is an open source library, written in Python, for parsing and interpreting the results of computational chemistry packages. The goals of cclib are centered around the reuse of data obtained from these programs and contained in output files, specifically

  • ADF (versions 2007 and 2013)
  • DALTON (versions 2013 and 2015)
  • Firefly, formerly known as PC GAMESS (version 8.0)
  • GAMESS (US) (version 2012)
  • GAMESS-UK (version 7.0)
  • Gaussian (versions 03 and 09)
  • Jaguar (versions 7.0 and 8.3)
  • Molpro (versions 2006 and 2012)
  • NWChem (versions 6.0 and 6.5)
  • ORCA (versions 2.9 and 3.0)
  • Psi (versions 3.4 and 4.0)
  • Q-Chem (version 4.2)

FragBuilder a tool to create, setup and analyse QM calculations on peptides. DOI.


And of course there is OpenBabel that can be used create input files for a variety of computational chemistry packages.

If I've missed anything please feel free to let me know.




LigParGen is a web-based service that provides force field (FF) parameters for organic molecules or ligands, offered by the Jorgensen group.

It is available here.

LigParGen provides bond, angle, dihedral, and Lennard-Jones OPLS-AA parameters with 1.14CM1A or 1.14CM1A-LBCC partial atomic charges. Server provides parameter and topology files for commonly used molecular dynamics and Monte Carlo packages OpenMM, Gromacs, NAMD, CHARMM, LAMMPS, CNS/X-PLOR, Q, DESMOND, BOSS and MCPRO. Also, the PQR file is generated. Supported input formats: SMILES, MOL and PDB.Maximum ligand size allowed is 200 atoms.

It is also possible to install LigParGen locally on Mac or Linux machines using Anaconda as described here

More details here

LigParGen web server: an automatic OPLS-AA parameter generator for organic ligands, DOI


Creating a Bioconda recipe


A little while back I mentioned BioConda. You can read more details in this publication "Bioconda: A sustainable and comprehensive software distribution for the life sciences", DOI. Conda is a platform- and language-independent package manager that sports easy distribution, installation and version management of software.

The conda package manager has recently made installing software a vastly more streamlined process. Conda is a combination of other package managers you may have encountered, such as pip, CPAN, CRAN, Bioconductor, apt-get, and homebrew. Conda is both language- and OS-agnostic, and can be used to install C/C++, Fortran, Go, R, Python, Java etc

The bioconda channel is a Conda channel providing bioinformatics related packages for Linux and Mac OS. Looking through the packages it is clear there it already contains a number of chemistry packages. These include: Updated 24 November 2017

  • OpenBabel
  • Rdkit
  • Opsin
  • chemfp
  • gromacs
  • osra
  • Autodock Vina
  • openmg
  • align-it
  • strip-it
  • shape-it
  • np-likeness-scorer
  • Smina

Bioconda offers a collection of over 3100 software tools, which are continuously maintained, updated, and extended by a growing global community of more than 330 contributors. Rather than try to duplicate this effort for a "Chemconda" it seems more efficient to encourage chemists to contribute to Bioconda. If you do package a chemistry application for Bioconda please let me know and I'll publicise it on my blog and add it to the list above. To start things rolling I've added to Bioconda and I've written a page describing how to create a bioconda recipe.

Link to page Creating a Bioconda recipe


Atomic Simulation Environment


The Atomic Simulation Environment (ASE) is a set of tools and Python modules for setting up, manipulating, running, visualizing and analyzing atomistic simulations. The code is freely available under the GNU LGPL license.

It requires

  • Python 2.7, 3.4-3.6
  • NumPy


  • SciPy
  • For ase.gui: Matplotlib (2D Plotting)

It can be installed using PIP

pip install --upgrade --user ase

Full details of the MacOSX installation are here.


SAMPL6 comp chem challenge update


Version 1.3: Add pKa prediction challenge instructions, input files, submission template files, update on the future plans of logD challenge.

SAMPL6 pKa Challenge Instructions

Challenge timeframe: Oct 25, 2017 to Jan 10, 2018

This challenge consists of predicting microscopic and macroscopic acid dissociation constants(pKa)s of 24 small organic molecules. These fragment-like small molecules are selected for their similarity to kinase inhibitors and for experimental tractability. Our aim is to evaluate how well current pKa prediction methods perform with drug fragment-like molecules through blind predictions.

Three formats of pKa prediction results will be evaluated:

  • microscopic pKa values and related microstates
  • microstate populations as a function of pH
  • macroscopic pKa values


The SAMPL6 Blind Prediction Challenge for Computational Chemistry


SAMPL6 includes challenges based on aqueous host-guest binding data (binding free energies and, optionally, binding enthalpies) for three different host molecules; and on physical properties (distribution coefficients and possibly solubilities), for a set of fragment-like molecules. The host-guest systems are useful to test simulation methods, force fields, and solvent models, in the context of binding, without posing the setup issues and computational burden of protein simulations. The physical properties offer efficient tests of force field accuracy when detailed simulations are used, and can also test pKa prediction methods, continuum solvation models, and knowledge-based prediction methods. SAMPL6 will also introduce a new challenge component, the “SAMPLing challenge”, in which computational methods will be evaluated on how efficiently their calculations approach well-converged reference results generated by the organizers. Participants will be provided with machine readable setup files for the molecular systems, including force field setups, along with recommended cutoffs and treatments of long-ranged interactions. The SAMPLing challenge is expected to include one or more cases from each challenge component (host-guest binding on each system; log D calculation).

If you would like to participate in the challenge join up here.


IM-UFF: extending the Universal Force Field for interactive molecular modeling


An interesting development for those working in nano materials.

We have completed the development of IM-UFF (Interactive Modeling - UFF), an extension of UFF that combines the possibility to significantly modify molecular structures (as with reactive force fields) with a broad diversity of supported systems thanks to the universality of UFF. Such an extension lets the user easily build and edit molecular systems interactively while being guided by physically-based inter-atomic forces. This approach introduces weighted atom types and weighted bonds, used to update topologies and atom parameterizations at every time step of a simulation. IM-UFF has been evaluated on a large set of benchmarks and is proposed as a self-contained implementation integrated in a new module for the SAMSON software platform for computational nanoscience.

This contribution has been submitted to the Journal of Molecular Modeling.

SAMSON is a novel software platform for computational nanoscience. Rapidly build models of nanotubes, proteins, and complex nanosystems. Run interactive simulations to simulate chemical reactions, bend graphene sheets, (un)fold proteins. SAMSON’s generic architecture makes it suitable for material science, life science, physics, electronics, chemistry, and even education. SAMSON is developed by the NANO-D group at INRIA, and means “Software for Adaptive Modeling and Simulation Of Nanosystems”.



Grand Challenge 3


I've written a couple of docking/virtual screening workflows using SMINA, a freely available tool DOI. There are a number of other alternatives and it is very difficult to get good comparisons which is why the Grand Challenges are useful.

The Grand Challenge 3 (GC3) is a blinded prediction challenge for the computational chemistry community, with components addressing pose-prediction, affinity ranking, and free energy calculations. GC3 is based on six different protein targets, Cathepsin S and five different kinases, and is separated into five subchallenges, some of which involve multiple protein targets. Only Cathepsin S is associated with cocrystal structures, so the kinase components of this challenge focus on affinity ranking and/or free energy predictions. Three of the datasets, Cathepsin S, JAK2 and TIE2, include a free energy prediction component.

This is an ideal opportunity to test novel algorithms on a carefully curated dataset.

Computational practices often employ a number of computational algorithms and dataset preparation steps for meaningful results. D3R will provide a forum for deposition, dissemination and discussion of such workflows through this website. Workflows will represent methods used successfully in the blinded challenges and methods donated by our pharmaceutical and/or academic collaborators. GitHub.

In Stage 1 September 1 - October 1, your predicted poses for the 24 ligands, in a coordinate system aligned with the S04-bound Cathepsin S structure provided in the inputs. Your predicted affinities, or affinity rankings, for all 136 compounds and/or your predicted absolute or relative binding affinities (in kcal/mol) for the free energy subset of 33 compounds. When Stage 1 closes, we will release the crystallographic poses of the 24 ligands.

In Stage 2 October 1 - December 1, your predictions of the affinity rankings of all 136 compounds and/or absolute or relative binding affinities (in kcal/mol) for the free energy subset of 33 compounds.

Full details and registration.


Facio GUI for CompChem


Facio is a free GUI for computational chemistry softwares (TINKER, MSMS, Firefly, Gamess, MOPAC and Gaussian).

Current Version of Facio : 20.1.2 (released on August 17 2017) for Mac (Mac OS X10.6 or later and macOS Sierra)

Pasted Graphic


Gabedit 2.5.0 released


The new version (2.5.0) of Gabedit has been released and is available for download at

Gabedit is a graphical user interface to computational chemistry packages like Gamess-US, Gaussian, Molcas, Molpro, MPQC, OpenMopac, Orca, PCGamess and Q-Chem. It can display a variety of calculation results including support for most major molecular file formats. The advanced "Molecule Builder" allows to rapidly sketch in molecules and examine them in 3D. Graphics can be exported to various formats, including animations.

Here is a list of the significant changes between 2.5.0 and 2.4.8:

  • Minor bugs fixed.
  • New tools for VASP :

    • read geometries (Optimization or M. Dynamic) from VASP OUTCAR file
    • read geometry from VASP POSCAR file
    • Create VASP POSCAR file
    • Read dielectric function from a VASP xml file and compute optic properties : the refractive index n(w), the extinction coefficient k(w), the absorption coefficient alpha(w), the reflectivity R(w), the energy loss spectrum L(w), and the optical conductivity sigma(w).

      • read data from vasprun.xml and plot DOS, pDOS and, Bands structures
  • Gabedit can now read the hessian from .hess orca file. After reading of the hessian, Gabedit computes frequencies, modes and effective masses.
  • Export in CChemI : update
  • Tv accepted (used by Gaussian and Mopac for periodic system). Using Tv, Gabedit can generate other cells.
  • deMon2k is now supported (Thanks to Dennis Salahub, Mauricio Chagas da Silva, Jonathan Kung and Morteza Chehelamirani for their suggestions, corrections, comments,...)
  • Gabedit can now compute the anharmonic spectrum by QM/MMFF94 method using iGVPT2 program. Gabedit can read the harmonic and anharmonic spectra from an iGVPT2 output file.
  • Energy, geometry optimization, MD, MD Conformations search by MMFF94, MMFF94s, UFF and Ghemical potentials are now supported by Gabedit via Open Babel.
  • Energy, geometry optimization, MD, MD Conformations search using your own program (potential) are now supported. DFTB+ is supported via this new tool.

Instructions for compiling under Mac OSX are here


Schrödinger Updated


The Schrödinger small molecule discovery suite has ben updated. This looks to be a substantial update and is described in the video below.

Supported MacOS X 10.12 and Mac OS X 10.9 - 10.11

3D Support, Supported: Interlaced stereo via Zalman 3D Monitors


CRYSTAL14 v1.0.4 released


Crystal14 Version 1.0.4 is mainly a bugfix release. All users are encouraged to upgrade to v1.0.4 as soon as possible.

This is the list of options which have been fixed in the v1.0.4 version of CRYSTAL14:

  • Electronic band structure for open-shell systems (asymmetric k-points only)
  • Piezoelectricity of 1D and 2D cases (non-periodic directions fixed)
  • Piezoelectricity for open-shell systems
  • BREAKELAST option restored
  • ATOMINSE+SUPERCELL combination fixed
  • Cell gradients for charged systems
  • Range-separated hybrids (conflict with bipolar approximation fixed)
  • Some static limits for CLUSTER option have been relaxed
  • Small fixes (normalization, k-points coordinates) in phonon bands and density-of-states
  • CUBE format for 3D plots of charge density and electrostatic potential fixed
  • Various fixes in restart options of vibration frequencies + IR and Raman intensities calculation
  • HJS exchange hole for open shell systems. This affects HSE06, HSEsol, HISS, LC-wPBE, LC-wPBEsol and LC-wBLYP functionals.
  • Fixes in the POTC option


Psi4 1.1: An Open-Source Electronic Structure Program


A recent paper describes Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability DOI

Psi4 is an ab initio electronic structure program providing methods such as Hartree–Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools.

Psi4 1.1 can be downloaded from here with versions supporting Python 2.7, 3.5 and 3.6.

Note the installation instructions for Mac: Install XCode via the App Store, Make sure you open XCode and accept the license agreement after you install.


Conformer generation


The generation of multiple conformations is an important step in a number of operations from input to ab initio calculations to providing input files for docking studies. A recent paper compared seven freely available conformer ensemble generators: Balloon (two different algorithms), the RDKit standard conformer ensemble generator, the Experimental-Torsion basic Knowledge Distance Geometry (ETKDG) algorithm, Confab, Frog2 and Multiconf-DOCK DOI, and also provided a dataset of ligand conformations taken from the PDB.

A recent twitter discussion involving Greg Landrum and David Koes prompted Greg to publish a blog post describing conformation generation within RDKit. The post compares using distance geometry to select diverse conformations versus an approach that combines the distance geometry approach with experimental torsion-angle preferences obtained from small-molecule crystallographic data (ETKDG). He also looks at the impact of force-field minimisation.

A really interesting read with code provided.


ADF Modeling Suite 2017


Just looking at the release notes for ADF2017 and support pages and a couple of things caught my eye.

The python distribution shipped with ADF was upgraded to python 3.5. Among others, new and updated modules include the iPython interpreter for easier development of python codes, a number of useful packages such as numpy 1.11.3 / scipy 0.18.1, ASE 3.13.0, matplotlib and RDKit 2016.09.

The Mac OSX version no longer requires XQuartz and is much faster in visualizing large systems.


ADF Modeling Suite 2017


I just got this message

We are proud to announce the 2017 release of the ADF Modeling Suite, with excellent contributions from our collaborators and the continued efforts of the SCM team in Amsterdam.

Exciting new features include

  • Many new NLO properties (TPA, THG, ...): Hu, Autschbach & Jensen
  • Constrained DFT with excited states: Ramos & Pavanello
  • LFDFT for d-d and d-f transitions: Ramanantoanina & Daul
  • CV-DFT for singlet-triplet excitations: Krykunov, Senn, Park & Seidu
  • Faster periodic response with TD(C)DFT, including 2D systems: Raupach
  • VCD analysis tools: Nicu

Reactivity & Analysis-

  • Latest xc functionals (SCAN, MN15-L, ....): interface to libxc 3.0
  • Special points, fat bands and improved pDOS analysis
  • GUI support for NEGF with BAND (Thijssen group, includes self-consistent NEGF, gate & bias potential, spin transport) and post-SCF DFTB-NEGF (Heine group)
  • FDE + local COSMO: Goez & Neugebauer
  • Reactivity descriptors from conceptual DFT and QTAIM: Tognetti & Joubert
  • Geometry optimization with SpinFlip in QUILD: Swart
  • Spin-polarization and l-dependency for DFTB: Melix, Oliveira, Rueger, Heine
  • Much faster periodic DFTB(+D) optimizations, latest parameters freely available
  • eReaxFF including explicit electrons: based on Islam, Verstraelen & van Duin
  • Controllable mass-scaling for force bias Monte Carlo ReaxFF: Bal & Neyts
  • Improved VLE, LLE, IDAC, kOW with reparameterized COSMO-SAC: Chen & Lin

GUI & Builders-

  • Quantum ESPRESSO: GUI interface & binaries
  • MOF builder and UFF4MOFsII: Coupry, Addicoat, Heine
  • Much faster visualization of large and periodic systems
  • Set up and visualize 'molecule gun' calculations with ReaxFF

For a more comprehensive list and details see:


ORCA 4.0 released


The latest update to ORCA has just been released.

The program ORCA is a modern electronic structure program package written by F. Neese, with contributions from many current and former coworkers and several collaborating groups. The binaries of ORCA are available free of charge for academic users for a variety of platforms. ORCA is a flexible, efficient and easy-to-use general purpose tool for quantum chemistry with specific emphasis on spectroscopic properties of open-shell molecules. It features a wide variety of standard quantum chemical methods ranging from semiempirical methods to DFT to single- and multireference correlated ab initio methods. It can also treat environmental and relativistic effects. Due to the user-friendly style, ORCA is considered to be a helpful tool not only for computational chemists, but also for chemists, physicists and biologists that are interested in developing the full information content of their experimental data with help of calculations.

New Features of Version 4.0:

New Methods:

  • Linear scaling DLPNO-CCSD(T) open shell. New restricted open-shell formulation
  • Linear scaling DLPNO-MP2 (RHF and UHF)
  • Linear scaling DLPNO-MP2-F12 (RHF)
  • Linear scaling DLPNO-CCSD(T) (the 2013 implementation is still available)
  • Linear scaling DLPNO-CCSD(T) local energy decomposition scheme
  • Linear scaling DLPNO-CCSD closed shell density
  • Linear scaling cluster in molecule (CIM): MP2, CCSD(T), DLPNO-CCSD(T)
  • Linear scaling DLPNO-NEVPT2
  • NEVPT2-F12
  • Updated interface to BLOCK 1.0
  • Closed shell EOM-CCSD energies
  • Closed shell STEOM-CCSD energies
  • Partial PNO-EOM-CCSD method for excited states
  • Partial PNO-STEOM-CCSD method for excited states
  • Mukherjee Mk-LPNO-MRCCSD(T)
  • Powerful iterative configuration expansion (ICE-CI) approximation to Full-CI
  • ICE-CI for large active space CASSCF calculations
  • MREOM-CCSD (also with SOC)
  • Fully internally contracted MRCI
  • Full TD-DFT energies and gradient for hybrid functionals
  • Super-fast approximate TD-DFT: sTDA/sTDDFT of Grimme and co-workers
  • PBEh-3c method of Grimme and co-workers

SCF, DFT and Hessian:

  • Large performance improvements for calculations with four center integrals
  • Improved performance with RI-J with conventionally stored integrals
  • Gradient for range separated hybrids
  • Gradient for range double hybrid functionals with meta GGAs
  • Gradient for range double hybrid functionals with range separated functionals
  • Gradient for RI-JK
  • Frequencies for range separated functionals
  • Stability analysis and automatic search for broken symmetry states
  • Local spin analysis
  • Fractional occupation number analysis (FOD) for detection of MR character

MDCI module:

  • All improvements for DLPNO methods as listed above
  • Closed shell EOM-CCSD energies
  • Closed shell STEOM-CCSD energies
  • Automatic closed shell STEOM-CCSD active space selection
  • EOM-CCSD(2) and STEOM-CCSD(2) approximations
  • EOM-CCSD transition moments
  • EOM/STEOM-CCSD core level excited states
  • ADC(2) and CC(2) methods (initial implementation)
  • Improved automatic frozen core handling
  • Core-correlation in automatic basis set extrapolation

AUTOCI module:

  • FIC-MRCI, CEPA/0 variant and DDCI3


  • Detailed tutorial showing CASSCF/NEVPT2 usage
  • Accelerated CI (ACCCI) a more efficient CI step for multi-root calculations
  • Automatic implementation of AbInitio ligand-field theory
  • Simplified generation of double-shell orbitals
  • Active space protection scheme and improved warnings
  • ICE-CI as CI solver for larger active spaces
  • Partially Contracted NEVPT2 with and without RI
  • Updated interface to BLOCK 1.0
  • DMRG-NEVPT2 for active spaces up to 20 orbitals
  • Magnetization and magnetic susceptibility
  • Printing of the wavefunction in terms of CSFs and spin-determinants
  • MREOM-CCSD (also with SOC)
  • Local spin analysis for CASSCF
  • Fragment decomposition of the spin-spin interaction
  • Cumulant approximation for NEVPT2
  • ACCCI as CIStep for FIC and DLPNO-NEVPT2
  • Explicitly correlated RI-FIC-NEVPT2 (NEVPT2-F12)


  • Full TD-DFT for hybrid functionals
  • Gradient for full TD-DFT with hybrid functionals
  • TD-DFT/TDA gradient with range separated functionals
  • ROCIS magnetic properties (hyperfine, g-tensor, ZFS tensor, MCD)
  • ROCIS-RIXS spectra
  • PNO-ROCIS for spectacular performance improvements
  • Super-fast approximate TD-DFT: sTDA/sTDDFT
  • Natural transition orbitals in TD-DFT and ROCIS


  • GIAO implementation for NMR chemical shifts. Various aproximations (RIJOCOSX, RIJK)
  • New Handling of basis set names. Now fully consistent with TurboMole def2-defaults (including ECPs) SARC basis sets separately available
  • New reading of basis sets and ECPs together
  • New correlation consistent basis sets added
  • New SARC basis sets for the lanthanides; good for correlated calculations
  • New ANO-RCC basis sets added
  • Improved frozen core handling in correlation calculations
  • Improved automatic auxiliary basis set generation
  • Corrections for low-frequency modes in thermochemistry
  • New and improved NBO interface
  • CPCM and improved SMD solvent models
  • Intrinsic atomic orbital (IAO) and bond orbital implementation
  • Improved performance in Boys localization
  • Updated and improved mapspc program
  • Atomic Mean Field (AMFI) spin-orbit coupling operators
  • EPRNMR works with range separated hybrid functionals
  • New molecular dynamics module


Conformation generation


Most drug-like molecules contain a number of rotatable bonds and prediction of bioactivities, docking etc. require an understanding of conformation. Whilst systematic methods can in theory explore all conformational space, however as the number of rotatable bonds increases a systematic search becomes prohibitive both in terms of the computational cost in generating conformations but also the time taken to process all the generated conformations (e.g. dock into protein). Thus there is great interest in rapid means to generate ensembles of representative conformations.


A recent paper DOI helps to address the problem by compiling a high quality dataset of structures generated using ligands from the Protein Data Bank.

The datasets were applied to benchmarking seven freely available conformer ensemble generators: Balloon (two different algorithms), the RDKit standard conformer ensemble generator, the Experimental-Torsion basic Knowledge Distance Geometry (ETKDG) algorithm, Confab, Frog2 and Multiconf-DOCK. Substantial differences in the performance of the individual algorithms were observed, with RDKit and ETKDG generally achieving a favorable balance of accuracy, ensemble size and runtime.

The dataset is freely available for download


MOE updated


Chemical Computing Group have announced and update to MOE. The MOE 2016.0802 update contains a number of updates to the biomolecule modelling including improved hydrogen bond detection, and addition of a number of unnatural amino acids.

There have also been improvements to MOE/Web MOE/web. The MOE/web version compatibility check has been broadened. MOE/web license waiting has been improved. HTTPS authentication proxy server support has been improved.


Tinker Updated


The TINKER molecular modeling software is a complete and general package for molecular mechanics and dynamics, with some special features for biopolymers. TINKER has the ability to use any of several common parameter sets, such as Amber (ff94, ff96, ff98, ff99, ff99SB), CHARMM (19, 22, 22/CMAP), Allinger MM (MM2-1991 and MM3-2000), OPLS (OPLS-UA, OPLS-AA), Merck Molecular Force Field (MMFF), Liam Dang's polarizable model, and the AMOEBA (2004, 2009, 2013) polarizable atomic multipole force field.

The TINKER package contains a variety of interesting algorithms such as: flexible implementation of atomic multipole-based electrostatics with explicit dipole polarizability, various continuum solvation treatments including several generalized Born (GB/SA) models, generalized Kirkwood implicit solvation for AMOEBA, an interface to APBS for Poisson-Boltzmann calculations, efficient truncated Newton (TNCG) local optimization, surface areas and volumes with derivatives, free energy calculations via the Bennett Acceptance Ratio (BAR) method, normal mode vibrational analysis, minimization in Cartesian, torsional or rigid body space, symplectic RESPA multiple time step integration for molecular dynamics, velocity Verlet stochastic dynamics, pairwise neighbor lists and splined spherical energy cutoff methods, particle mesh Ewald (PME) summation for partial charges and polarizable multipoles, a novel reaction field treatment of long range electrostatics, fast distance geometry metrization with better sampling than standard methods, Elber's reaction path algorithm, potential smoothing and search (PSS) methods for global optimization, Monte Carlo Minimization (MCM) for efficient potential surface scanning, tools for fitting charge, multipole and polarization models to QM-based electrostatic potentials and more....

TINKER 8 is a major new release of the Ponder Lab tool set for molecular mechanics and dynamics calculations. An important change in this new version is the switch from old-style common blocks to Fortran modules. Use of modules and greatly increased use of dynamic memory allocation means TINKER can now support very large molecular systems. TINKER 8 also implements improved OpenMP parallelization throughout many parts of the code. Additional big improvements include parallel neighbor list building and updating, and big reduction in iteration needed to converge AMOEBA polarization via an efficient PCG solver. Other changes from the previous TINKER version include new and updated force field parameter sets and numerous minor additions and bug fixes, many of them suggested by users of the package. Please note that as with prior new releases, version 8 is neither backward nor forward compatible with earlier versions of TINKER. In particular, older versions of parameter files should not be used with TINKER 8 executables and vice versa.


Publishing computational notebooks with Binder


I've now written a couple of Jupyter notebooks and one of the issues that has come up is how to share the notebooks in a way that ensures the results will be reproducible in an environment when updates to components occur regularly.

Binder is a collection of tools for building and executing version-controlled computational environments that contain code, data, and interactive front ends, like Jupyter notebooks. It's 100% open source.

At a high level, Binder is designed to make the following workflow as easy as possible

  • Users specify a GitHub repository
  • Repository contents are used to build Docker images
  • Deploy containers on-demand in the browser on a cluster running Kubernetes

Common use cases include:

  • sharing scientific work
  • sharing journalism
  • running tutorials and demos with minimal setup
  • teaching courses


If you want to find out more have a look at this blog post by the developers.


Maestro 11 training


A nice image from the Maestro 11 training session run by Schrödinger, it looks like there are a few more training sessions coming up also.


Maestro 11 is the portal to all of Schrödinger's computational technology – far more than just a user interface, Maestro 11 also helps researchers organize and analyze data.


Added PUPIL to alphabetical listing


PUPIL,Program for User Package Interface and Linking, is a software environment - the program - that allows developers to link quickly and efficiently together multiple pieces of software in a fully automated multi-scale simulation. More specifically, it supports QM/MM MD simulations where the user might choose among any of the different MD engines and QM engines, which are connected to PUPIL as external programs through a tiny specific interface. One of the main advantages here is that the user can use most of the functionalities that may have those external programs interfaced without the necessity to be implemented again on independent interfaces. In fact, this simulation interface concentrates all the common code involved in the coupling terms of the QM/MM approach.


ICM version 3.8-5


MolSoft have announced the release of ICM version 3.8-5.

  • Generate a 2D Interaction Diagram of a ligand with the binding pocket. The image is annotated with hydrogen bonds and interacting residues.
  • 3D ligand editor is a powerful tool for the interactive design of new lead compounds in 3D
  • ICMJS is a JavaScript/HTML5 viewer for 3D Molecular Graphics which does not require any plugin or installation.
  • Support for MMTF format. The Macromolecular Transmission Format (MMTF)
  • Support for Mac retina display
  • Add docking restraints by selecting atoms in the receptor
  • Updates to protein modelling, bioinformatics and cheminformatics

Full release notes are here


MOE 2016.08 released


Chemical Computing Group have just announced an update to MOE. This release has fixed a couple of Mac OSX 10.12 (Sierra) issues but also brings a host of new features.

  • MOEsaic: Web-Application for Ligand Analytics
  • Spectral Analysis for Structure Determination
  • Enhanced Protein Patch Analyzer
  • Integrated Antibody Project Database and Antibody Homology Modeler
  • Small Footprint MOE to Facilitate Large Scale Deployments
  • Physical and Virtual Rendering of Structures

A more detailed description of the new and enhanced features in MOE 2016.08 can be found at


APBS 1.5 released


The latest version of APBS includes several notable features and bug fixes. This release includes the addition of Poisson-Boltzmann Analytical-Method (PB-AM), Poisson-Boltzmann Semi-Analytical Method (PB-SAM) and the Treecode-Accelerated Boundary Integral Poisson-Boltzmann method (TABI). Additionally, we have made improvements to the build system and the system tests, as well as miscellaneous bug fixes.

APBS & PDB2PQR: Electrostatic and solvation properties from complex molecules. Solve the Poisson-Boltzmann and related equations to calculate solvation energies and electrostatic properties for analysis and visualization

APBS 1.5 changes

Binary releases may be found on GitHub and on SourceForge. New Features

  • Poisson-Boltzmann Analytical Method (PBAM, see Lotan & Head-Gordon) and Semi-Analytical Method (PBSAM, see Yap & Head-Gordon) integrated with APBS.
  • PBSAM is currently only available in the Linux and OS X distributions.
  • Examples are located with the APBS examples in the pbam/ and pbsam/ directories.
  • More information and documentation may be found in the PBAM and PBSAM sections of the APBS-PDB2PQR website.
  • Tree-Code Accelerated Boundary Integral Poisson-Boltzmann Method (TABI-PB) integrated with APBS.(See Geng & Krasny)
  • Examples are located with the APBS examples in the bem/, bem-pKa/, and bem-binding-energies/ folders
  • Included NanoShaper alternative to MSMS.
  • More information and documentation may be found in the Contributions section of the APBS-PDB2PQR website
  • Added binary DX format support to the appropriate APBS tools.
  • Test suite amended and expanded.
  • Removed hard-coded limitation to number of grid points used to determine surface accessibility.

Known Bugs / Limitations

  • PBSAM not building in windows due to C standard restrictions in the Microsoft compiler implementation.

Full details here


QUACPAC v1.7.0 released


Even the best charge models are useless if protonation states are wrong. QUACPAC attempts to offer everything necessary to do charges correctly. It includes pKa and tautomer enumeration in order to get correct protonation states, partial charges using multiple models that cover a range of speed and accuracy, and electrostatic potential map construction and storage.

This update contains amount other things

  • Significant improvements have been made to the reasonable tautomer algorithm that affect its aliphatic and non-aromatic resonance portions.
  • Inclusion of AM1BCC ELF10, a new method for applying partial charges to a molecule
  • TIP3P water charges are now assigned when using Amber charge sets on molecules containing waters.

Support for Mac OS X 10.10 and 10.11 has been added.
Mac OS X 10.7 and 10.8 are no longer supported.

Full release notes are here


BioExcel: An Introduction


An overview of the BioExcel Project.

BioExcel is based on improving three aspects of biomolecular research. Firstly, improving the performance and scalability of the most commonly used software, such as GROMACS (, HADDOCK ( and CPMD (, to take advantage of next-gen HPC systems and the expected increase in the amount of data produced. It’s also important to improve how easy it is for users to access and use these types of software. Not all researchers have experience in efficiently handling data and software. BioExcel aims to provide customisable workflow environments, which will allow relatively novice HPC/HTC users take advantage of the analysis software provided in ways that suit their specific research. In addition to this, hands-on training and public webinars are already underway, aiming to teach researchers best practices and how to best utilise the software and resources available.




Facio is a free GUI for computational chemistry softwares (TINKER, MSMS, Firefly, Gamess, MOPAC and Gaussian).



GPView Wave Function Analysis and Visualization


The GPView program is a C++ package for wave function analysis and visualization.

It is developed and maintained by Tian Shi and Ping Wang Ref](

In this manuscript, we will introduce a recently developed program GPView, which can be used for wave function analysis and visualization. The wave function analysis module can calculate and generate 3D cubes for various types of molecular orbitals and electron density related with electronic excited states, such as natural orbitals, natural transition orbitals, natural difference orbitals, hole-particle density, detachment-attachment density and transition density. The visualization module of GPView can display molecular and electronic (iso-surfaces) structures. It is also able to animate single trajectories of molecular dynamics and non-adiabatic excited state molecular dynamics using the data stored in existing files. There are also other utilities help to extract and process the output of quantum chemistry calculations. The GPView provides full graphic user interface (GUI) which makes it very easy to use.



Amber16 and AmberTools16 released


AmberTools consists of several independently developed packages that work well by themselves, and with Amber itself. The suite can also be used to carry out complete molecular dynamics simulations, with either explicit water or generalized Born solvent models.

The AmberTools suite is free of charge, and its components are mostly released under the GNU General Public License (GPL). A few components are included that are in the public domain or which have other, open-source, licenses. The sander program now has the LGPL license. AmberTools is distributed in source code format, and must be compiled in order to be used. You will need C, C++ and Fortran90 compilers

The Amber16 package builds on AmberTools16 by adding the pmemd program, which resembles the sander (molecular dynamics) code in AmberTools, but provides (much) better performance on multiple CPUs, and dramatic speed improvements on GPUs. Major new features include:

Semi-Isotropic Pressure Scaling (GPU)
Charmm VDW Force Switch (CPU, GPU)
Enhanced NMR Restraint support + R^6 averaging support (GPU)
Gaussian Accelerated Molecular Dynamics (CPU, GPU)
Support for external electric fields (CPU)
Expanded umbrella sampling support (GPU)
Constant pH supported with replica exchange along pH coordinate (GPU)
Support for gas phase MD (igb=6) (CPU, GPU)
Support and significant performance improvements for the latest Kepler, Maxwell and Pascal GPUs from NVIDIA.


PDB2PQR 2.1 released


A new version of PDB2PQR has been released.

APBS (Adaptive Poisson-Boltzmann Solver) and PDB2PQR are software packages designed to help you analyze the solvation properties of small and macro-molecules such as proteins, nucleic acids, and other complex systems

  • Added alternate method to do visualization using 3dmol.
  • Replaced the Monte Carlo method for generating titration curves with graph cut. See (If you prefer the Monte Carlo method, please use
  • Added compile options to allow for arbitrary flags to be added. Helps work around some platforms where scons does not detect the needed settings correctly.
  • Added a check before calculating pKa's for large interactions energies.

OSX binaries require OSX 10.6 or newer. The OSX binary is 64-bit.

Dolinsky TJ, Czodrowski P, Li H, Nielsen JE, Jensen JH, Klebe G, Baker NA. PDB2PQR: Expanding and upgrading automated preparation of biomolecular structures for molecular simulations. Nucleic Acids Res, 35, W522-5, 2007. DOI

Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup, execution, and analysis of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res, 32, W665-W667, 2004. DOI


ADF Modeling Suite 2016


The ADF modelling suite has been a popular modelling package used in many areas of chemistry and materials science. SCM have recently announced an important update to ADF

Software for Chemistry & Materials (SCM) is an Amsterdam-based computational chemistry software company. Originally spinning out from the Vrije Universiteit as Scientific Computing & Modelling NV in 1995, the SCM team supports and develops the ADF Modeling Suite, centered around the flagship program Amsterdam Density Functional (ADF), which was originally developed in the 1970s in the theoretical chemistry department.

Key new features::

  • New XC functionals: range-separated hybrid HSE06, long-range corrected hybrids and new meta-GGAs
  • SM12 solvation model and Constrained DFT in ADF
  • Spectroscopy: excitations from CV(n)-DFT, fast TDDFT+TB and sTDDFT, X-ray emission, surface-enhanced ROA, new kernel for periodic TDDFT, vibrationally resolved spectra from TDDFTB, hole states in BAND
  • Improved robustness for SCC-DFTB, COSMO, periodic optimizations, new HF exchange scheme
  • Analysis: unrestricted and periodic energy decomposition (pEDA), Fukui functions, Natural Transition Orbitals
  • Reactivity, PES: Automated reaction pathways for ReaxFF trajectories, new Transition State search options, MECP, analytic lattice gradients
  • GUI: job chaining, multiple spectra, orbital interaction visualization
  • Scripting: ASE interface for all codes, scripting support for COSMO-RS, job chaining, extensions to FlexMD multi-scale dynamics
  • Parameters, databases: Quasinano15 including repulsive potentials for light elements, latest 3OB parameters, new parameter sets for ReaxFF, ionic liquids database COSMO-RS

Full release notes are here




iSpartan has been updated.

What's New in Version 1.4.6 Improved stability for iOS9 for Spartan'14 Parallel Suite (Serve)r users Enhanced charge labels in Sketch mode Increased electrostatic potential surfaces integrity

iSpartan is a versatile app for molecular modeling on the iPad, iPhone, and iPod Touch. The app allows any chemist fast and easy access to computational methods that have proven reliable for a large range of molecules.



OpenEye Toolkit Updated


OpenEye have announced the release of OpenEye Toolkits v2016.Feb. These libraries include the usual support for C++, Python, C# and Java.

The update address several key features.


OpenEye toolkits are used in web services that require protection from malicious users. The most obvious attack vector against the OpenEye toolkits is file format parsing since scientific file formats are complex and often underdefined and there is the potential for embedded malicious code. This update closes a number of potential vulnerabilities.

FastROCS TK: Database Loading Performance

An interesting development is that physical memory limits on GPU's mean that for loading larger libraries the loading of the dataset actually takes longer than the actual search. This release addresses that issue.

OEMedChem TK

This also contains first official release of OEMedChem TK, in particular access to matched molecular pairs.

This 2016.Feb release no longer support OSX 10.8 , but support has been added for OSX 10.11. This 2016.Feb release supports Python 3.5 for the following platforms: OSX 10.10, OSX 10.11, Ubuntu 12, Ubuntu 14, RedHat 6, and RedHat 7

Full release notes are here ….




Added NCIPLOT to the alphabetical listing of applications. NCI (Non-Covalent Interactions) is a visualization index based on the density and its derivatives. It enables identification of non-covalent interactions. NCIPLOT is available for download.

Our approach reveals the underlying chemistry that compliments the covalent structure. It provides a rich representation of van der Waals interactions, hydrogen bonds, and steric repulsion in small molecules, molecular complexes, and solids. Most importantly, the method, requiring only knowledge of the atomic coordinates, is efficient and applicable to large systems, such as proteins or DNA. Across these applications, a view of nonbonded interactions emerges as continuous surfaces rather than close contacts between atom pairs, offering rich insight into the design of new and improved ligands.


Erin R. Johnson, Shahar Keinan, Paula Mori-Sanchez, Julia Contreras-Garcia, Aron J. Cohen, and Weitao Yang, J. Am. Chem. Soc. 2010, 132, pp 6498-6506.
J. Contreras-Garcia, E. R. Johnson, S. Keinan, R. Chaudret, J-P. Piquemal, D. N. Beratan, and W. Yang. J. Chem. Theory Comput. 2011, 7, pp 625-632.




WebMO version 16

WebMO version 16 has been officially released and is available for FREE download at

New features in WebMO version 16 include:

  • Support for both iOS and Android WebMO free apps
  • Support for the newest versions of ORCA, PSI, and Q-Chem
  • Support for hydrogen bonds in the WebMO editor
  • Improved support for building/cleanup of fused aromatic rings
  • Incorporation of modern HTML5 user interface elements including: drag-and-drop of jobs between folders, fewer pop-up windows, live updating of output files while jobs run, and expandable results page
  • Support for unlimited number of user folders and sub-folders (Pro)
  • Enhanced support for resource specification with external batch queues (Enterprise)
  • Various bug fixes

Important changes in the previous version include:

  • Java security warnings are reduced/eliminated with new signed applet
  • Support for the newest versions of GAMESS, MOPAC
  • Lookup feature to build structures by name (type "aspirin" or "vancomycin")
  • Links to external databases to lookup molecular properties, spectra, and data
  • Huckel molecular orbital calculations from within Build Molecule page (Pro)

WebMO is the most popular interface to computational chemistry programs, with over 22,000 unique WebMO licenses issued to date!


WebMO has released a new Android app and updated its popular Apple iOS app!

The WebMO apps for portable devices are FREE.  

The Android and iOS WebMO apps provide equivalent features including:

  • Standalone molecular editing, optimization, symmetry, and orbital calculations
  • Integrated web lookup of molecular properties, spectra, and data
  • Access to WebMO servers (Free, Pro, and Enterprise versions of WebMO 16 and higher)

The WebMO apps supports all chemistry courses by calculating and displaying:

  • VSEPR structures for General Chemistry
  • 3-D structures to replace model kits for Organic Chemistry
  • Data and properties for Analytical Chemistry
  • Point group and symmetry elements for Inorganic Chemistry
  • Small proteins and bimolecules for Biochemistry
  • Molecular orbitals and vibrations for Physical Chemistry
  • WebMO calculations for undergraduate and graduate level Research

Thousands of students and scientists use the WebMO apps on their smartphones and tablets each month.  You can join them today for free!


More computational services are moving away from local hardware and into the "cloud".  Cloud computing requires zero upfront capital outlay, eliminates hardware maintenance, and can dramatically lower the cost of computing.

WebMO now supports installation on cloud computing platforms.  In less than 10 minutes, you can use your web browser to create a virtual machine, run the WebMO SITC ("Server in the Cloud") script, and have a functioning WebMO 16 server.

The cost of a WebMO cloud server can be amazingly low: <$1 per day for a WebMO server capable of running jobs on standard ab initio chemistry engines, and just $0.12 per day for a WebMO server running MOPAC!  And when you don't need it, you just turn it off.

Complete documentation is available at

Since Google Compute Engines are free for the first 60 days, so there is no cost or risk to try cloud computing today.


BCL::Conf: small molecule conformational sampling using a knowledge based rotamer library


Sampling conformational space is a key requirement for several areas of ligand design in small molecule drug discovery. A recent paper BCL::Conf: small molecule conformational sampling using a knowledge based rotamer library DOI describes a new conformational search method.

The project homepage provides a download, supports Mac OSX 10.4 or higher.


CCP4 released


The current version is CCP4 7.0 (07 January 2016). The new packages include:
SHELX suite: co distribution for academic users
CCP4I2: new ccp4 interface
DIALS: data processing and integration
ARCIMBOLDO-LITE: molecular replacement pipeline

CCP4 exists to produce and support a world-leading, integrated suite of programs that allows researchers to determine macromolecular structures by X-ray crystallography, and other biophysical techniques.

The latest version supports Mac OSX and Linux and requires X11. Source code is also available.


An early Christmas present from Chemical Computing Group


Chemical Computing Group have just released an up date to MOE, version 2015.10 includes:-

Protein-Protein Docking

  • Generate docked poses using FFT followed by all atom minimization
  • Define receptor and ligand sites to focus docking
  • Automatically detect antibody CDR sites

Integrated Alignment, Consensus and Superposition in the Sequence Editor

  • Manipulate multimeric protein sequences using split side-by-side Sequence Editor panes
  • Use dendrograms to visualize pairwise similarity, identity and RMSD relationships
  • Select residues based on plotted values using resizable sequence editor plots

Distributed Pharmacophore Searching

  • Run pharmacophore searches on a cluster directly from MOE GUI
  • Perform fast corporate database searches
  • Access multiple databases stored on a central server

Covalent Docking and Electron Density Docking

  • Use reaction-based organic transformations to covalently docking
  • Minimize ligand strain energy while maximizing ligand fit to electron density
  • Run docking through an enhanced streamlined scenario-based interface

Extended Hückel Descriptors and pKa Model

  • Compute molecular properties such as logP, logS and molar refractivity
  • Determine populations of ligand protonation states at a given pH
  • Calculate the pKa and pKb of small molecules

13C NMR Analysis

  • Apply QM conformation refinement to calculate 13C NMR shielding
  • Convert computed shieldings and predict 13C NMR chemical shifts
  • Compare computed chemical shifts to experimental shifts for structure determination

I'll write a review in the New Year.


BCL::Conf: small molecule conformational sampling using a knowledge based rotamer library


I just thought I'd flag a recent paper on sampling conformation space DOI.

The interaction of a small molecule with a protein target depends on its ability to adopt a three-dimensional structure that is complementary. Therefore, complete and rapid prediction of the conformational space a small molecule can sample is critical for both structure- and ligand-based drug discovery algorithms such as small molecule docking or three-dimensional quantitative structure–activity relationships.

The software is available from the MeilerLab home page


OpenEye Toolkits v2015.October released


OpenEye have announced the release of OpenEye Toolkits v2015.October. These libraries include the usual support for C++, Python, C# and Java.

New Features

  • FastROCS TK was added to the OpenEye toolkits collection
  • Molecule reading performance improvement in OEChem TK
  • The capabilities of the OEBio-Fragment Network have been expanded
  • 213 new ring templates have been added to the OEChem TK built-in ring dictionary

The full release notes give more details

In particular note the 2015.Oct release is the last to support Mac OSX 10.8 so time to upgrade if you have not already done so.


Computational chemistry guides & tools


The Medicines for Malaria Venture have an interesting page in which they are accumulating a list of computational tools and guides describing the use of the tools to address key issues within the drug discovery process.

Tools were chosen to address common needs expressed by medicinal and computational chemists working in the not-for-profit area. Recognising that this is a global effort, we have selected software packages on the basis of being free for all users.

The guides are either text descriptions or webcasts showing the tool in action. To date they include DataWarrior, KNIME, YASARA, ChEMBL and PK Tool.


Dalton 2015 entry updated


The Dalton2015 suite consists of two separate executables, Dalton and LSDalton. The Dalton code is a powerful tool for a wide range of molecular properties at different levels of theory, whereas LSDalton is a linear-scaling HF and DFT code suitable for large molecular systems, now also with some CCSD capabilities.

Together, the two programs provide an extensive functionality for the calculations of molecular properties at the HF, DFT, MCSCF, and CC levels of theory. Many of these properties are only available in the Dalton2015 suite.

Dalton can be run on a variety of systems running the UNIX operating system. The current release of the program supports Linux, Cray, SGI, and MacOS using GNU or Intel compilers (we plan to publish patches for PGI and XL compilers).. The program is written in FORTRAN 77, FORTRAN 90 and C, with machine dependencies isolated using C preprocessor directives. All floating-point computations are performed in 64-bit precision, but if 32-bit integers are available the code will take advantage of this to reduce storage requirements in some sections.


wxMacMolPlt updated


The popular GUI for GAMESS MacMolPlt has now been renamed and moved to GitHub, wxMacMolPlt is a cross-platform (Mac OS X, Linux and Windows) GUI for preparing, submitting and visualizing input and output for the GAMESS quantum chemistry package. Features include a graphical molecule builder, GAMESS input generation, animation of output and visualization of molecules, normal modes, orbitals and other properties


wxMacMolPlt can be downloaded from here

Bode, B. M. and Gordon, M. S. J. Mol. Graphics Mod., 16, 1998, 133-138. DOI


SeeSAR version 3.2 released


SeeSAR 3.2 can now be used in cases where there is no protein. For example if for the analysis of a ligand-based virtual screen


Release Notes version 3.2 2015-07-24

Utilizing SeeSAR without a protein

  • SeeSAR has grown from a single purpose affinity assessment tool to a multi purpose 3D structural viewer for compound design and prioritization. With this update it is now also possible to use all these nice features for 3D ligand-based projects. You may, for example, visually inspect small molecules alignments. You may filter a hit list by means of calculated and external properties...
  • Big data booster, With version 3.0 we first equipped SeeSAR with database functionality. Version 3.2 comes with a load of performance enhancements that speed up the calculation by up to a factor of 5, now utilizing all available CPUs in your computer even more efficiently.
  • 3D graphics enhancement, In this version we updated the graphics support and SeeSAR is now compatible with more graphics cards than ever before. Especially the compatibility with integrated graphics cards (the type most frequently found in laptop computers) - which used to be the major trouble makers - has been greatly improved.

Kinetiscope, interactive method for the accurate simulation of chemical reactions.


Kinetiscope is a scientific software tool that provides the bench scientist with an easy-to-use, rapid, interactive method for the accurate simulation of chemical reactions.


The program package is completely self-contained, and requires no programming or extensive user training. This lets you become immediately productive,

Kinetiscope comes with a User's Manual, a set of tutorials and an extensive library of example simulations that show you the kinds of systems that can be studied with it and that illustrate techniques for handling various types of physical situations. The library includes simulations of gas phase, solution phase and solid state reactions such as co- and terpolymerization ... radical chain-initiated polymerization (including a sample spreadsheet for extracting molecular weight distributions) ... kinetic resolution of enantiomeric mixtures ... chemistry in supercritical media ... pH-dependent model enzyme kinetics ... thermogravimetric analysis ... temperature programmed desorption ... smog chemistry ... silane chemistry in a chemical vapor deposition reactor ... model batch and flow catalytic reactors ... curing of polymers with significant volume shrinkage ... synthetic protocols for preparation of a photosensitizer ... chemical oscillators ... electrochemical reactions studied by cyclic voltammetry ... photochemical reactions from a pulse light source ... pharmacokinetics of drug dosing ... and imaging chemistry in photoresist materials.

Kinetiscope is a 64-bit application and requires Mac OS X 10.6 or above.

  • D. L. Bunker, B. Garrett, T. Kliendienst and G.S. Long III, Combustion and Flame 23, 373 (1974).
  • D. T. Gillespie, Journal of Computational Physics 22, 403 (1976).

Workshop on Resources for Computational Drug Discovery


I just thought I'd mention this upcoming workshop, as ever there is a top class line up.

EMBL-EBI/Wellcome Trust Workshop on Resources for Computational Drug Discovery 2-6 November 2015 Wellcome Trust Genome Campus, Hinxton, Cambridge, UK

The draft agenda is here


SeeSAR Updated


SeeSAR has been updated to version 3.1, the release notes highlight two significant new features.

SeeSAR is a software tool for interactive, visual compound prioritization as well as compound evolution.

  • Working with "big data" With this update we lifted the limit of handling only a maximum of 5000 poses in SeeSAR. We know that a lot of people like to do their compound analysis and prioritization after virtual screening campaigns also with much bigger sets. It is not likely that you will look at more than a couple of hundred poses, however, since the filtering (see also below) is extremely efficient, it provides quite an attractive opportunity to load all your data (not just the top x) and do your prioritization with all properties at hand right here in SeeSAR.
  • Enhanced filtering Behind the scenes SeeSAR knows so much more about your compounds than what is displayed in the table. The basic stuff like no. of acceptors and donors, rotatable bonds, etc. to do the usual Lipinski-type filtering is of course available, but also more elaborate stuff like the number of hydrogen bonds formed or the number of torsions that lie outside the statistical "norm". All of these are now available for filtering to help you optimally trim down your data to find the really interesting part.

NOTE! SeeSAR project files from older versions are incompatible and cannot be loaded. By default SeeSAR puts a new version in a separate location. The recommendation is to export your data from the old project file with the old version and import it into the latest SeeSAR. This is a one-time effort, which allows you to benefit from the features of the most up-to-date version.


Findings and SZMAP updated


Findings 1.2.3 the electronic notebook has just been released. Already using Findings? Menu Findings Check for updates...

OpenEye have just announced the release of SZMAP v1.2.1. This is a critical upgrade that fixes several bugs.

  • The Water Orientation VIDA Extension now saves orientation probability data and ensemble energy, vdw and order data on probe molecules as SDData. When probes generated by this extension are written to an .sdf file, these values will be written as well.
  • The Water Orientation VIDA Extension (version 1.1.3) now displays the correct water orientations for results from szmap version 1.2.0. It is also more compatible with 2D depictions in the 3D display.
  • gameplan no longer crashes when all stabilization test sites clash.
  • pch option -nonsymmetrized_charges now works to control AM1BCC partial charges. The default has been changed to true so as not to alter the default behavior.
  • pch no longer treats CYS residues bonded through the sidechain to a non-CYS as anionic.
  • pch no longer duplicates non-standard groups containing a metal (such as cofactors) and groups they are bonded to.
  • The hydrogen charge on the standard water probe is now correctly listed as +0.327 not +0.237 in SZMAP Theory.

Cresset update Torch and Forge


The update of Forge, a computational chemistry workbench for ligand-based design, includes over 170 new or improved features. Of particular note is Activity Atlas a new component enables you to summarize the SAR for a series into a visual 3D model that can be used to aid new molecule design. Forge V10.4 can now connect to an external web service, through a REST interface, to import external properties and data computed or retrieved by such web services as additional columns in the Molecules Table.

The latest version of Torch now includes Multi-parameter Optimization (MPO) options including condensing many activity and physicochemical properties into a single score representing the fit to the project profile, this has been coupled to an improved radial plot visualisation and tile display.


OpenEye toolkits updated


OpenEye has announced the release of OpenEye Toolkits v2015.June. These libraries include the usual support for C++, Python, C# and Java and are now available for download.

New Features Highlights:

  • PDB Splitting in OEBio TK
  • PAINS (Pan Assay Interference Compounds) filter in OEMolProp TK
  • Matched molecular pair improvements in OEMedChem TK
  • Custom ring template dictionaries in OEChem TK
  • Anaconda support for easier Python toolkit installation



The recently updated AMBER tools now includes software to analyse the structure and thermodynamics of water at protein and other surfaces. In addition to 3D-RISM, AMBER Tools now includes a new method, called Grid Inhomogeneous Solvation Theory, uses the same underlying theory as WaterMap but generates 3D grids of water properties, rather than focusing on discrete hydration sites.

The method is described in this paper, Grid inhomogeneous solvation theory: Hydration structure and thermodynamics of the miniature receptor cucurbit[7]uril DOI.

Solvation and desolvation plays a critical role in ligand binding but can be difficult to determine computationally.

AmberTools consists of several independently developed packages that work well by themselves, and with Amber itself. The suite can also be used to carry out complete molecular dynamics simulations, with either explicit water or generalized Born solvent models. AmberTools is distributed in source code format, and must be compiled in order to be used. You will need C, C++ and Fortran90 compilers.


Gamess-US on a Raspbery Pi


GAMESS is a program for ab initio molecular quantum chemistry. Briefly, GAMESS can compute SCF wavefunctions ranging from RHF, ROHF, UHF, GVB, and MCSCF. In an interesting development it is now possible to compile Gamess-US on the $40 Raspberry Pi, full instructions are available here. This provides a low cost system for demonstrating quantum chemistry.

To aid construction of the chemistry input WebMO-14.0, a web-based interface for computational chemistry programs, now has an App for iPhone/iPad. This allows you to compose e.g. a Gamess computation on your iPad or iPhone and send it to your Raspberry Pi (in the local network or on the web with a forwarded IP) to run the job and get the results back for visualizing and printing.



MacMolPlt 7.6


I just got this message regarding MacMolPlt a graphics program for plotting 3-D molecular structures and normal modes (vibrations).

I have just finished posting the final binaries for MacMolPlt 7.6. In addition to the code changes I have had to move the home site and binaries download sites due to the impending shutdown of Google code.

This version includes the following changes:

  • Energy plot window now accepts the same keystrokes as the main display window to change frames (left, right, home, end)
  • Fix a parsing issue with ROHF GAMESS log files
  • Account for an $EFRAG group change to allow all fragment atoms rather than only the first three.
  • Fixed a crash in the frequency window (when no normal modes are present)
  • Fixed the positioning of lone pairs in the builder when the coordination number plus the number of lone pairs is 5.
  • Fixed a couple of parsing issues with MolDen and Molekel (mkl) files.
  • Added GIF and TIFF export image formats (requires wxWidgets 2.9 or newer).

I have also reworked the 3D orbital generation code to significantly improve performance. Please report anything that doesn't look quite right.

Binaries and source files are available at:


Lucas, Molecule viewer and editor


Lucas is a novel program for graphical display and editing of molecular systems. The program allows fast and easy building and/or editing different molecular structures, up to several thousands of atoms large. Luscus is able to visualise dipole moments, normal modes, molecular orbitals, electron densities and electrostatic potentials. In addition, simple geometrical objects can be rendered in order to reveal a geometrical feature or a physical quantity. The program is developed as a graphical interface for the MOLCAS program package, however its adaptive nature makes possible to use luscus with other computational program packages and chemical formats. All data files are opened via simple plug-ins which makes easy to implement a new file format in luscus. The easiness of editing molecular geometries makes luscus suitable for teaching students chemical concepts and molecular modelling.


Journal of Cheminformatics 2015, 7:16 [DOI](}

The source code is available on Sourceforge


Scientific Applications under Yosemite


I just thought I'd like to thank all those who contributed to the Scientific Applications under Yosemite web page, many users and developers contacted me either via email or in the comments section and they certainly added information about applications that I don't have access to.

To date the page has been viewed well over 10,000 times with readers from 188 different countries. Viewers spent an average of just under two minutes on the page and it still attracts 800 pages views a month.

Given that 75% of the visitors to the site are now using Yosemite I suspect most scientists have now made the transition and I won't be updating the page any more. Once again thanks for the contributions.


MOE updated


MOE2014.0901 Update is now available. MOE is a fully integrated molecular modelling and drug discovery software package.

MOE 2014.0901 updates:


Protein Builder

  • Option for AMBER residue name
  • Append/prepend multiple residue sequence specified by single-letter names Builder:
  • Added H’s inherit color if there is a consistent coloring in the residue

sddesc: New -smi:p option causes field headers to be written to the output ASCII file

Bug Fixes:

  • MOESVLRUNPATH now properly honored
  • Combinatorial Builder now honors different attachment point locations on the same R-group
  • Database Save As one entry per file mode now properly generates unique filenames
  • Dock Template Forcing batch file now correctly generated
  • Saved views in .moe files now properly restored
  • Auto-save when Database Viewer display attributes are changed can now be disabled to prevent changes to the database file modification date when only the display is changed and not the database content
  • SVL function Deprotonate now works properly
  • Various MOE Project and Project Database Update bugs
  • Various minor bug fixes

There are reviews of MOE available here

Moe:- Molecular modeling
Moe Update (Jan 2009):- Molecular modeling
Review of MOE (2009.10 release):- Molecular modeling
Moe Update (December 2010.10 release):- Molecular modeling
Moe Update (December 2011 release):- Molecular modeling
Moe Update (December 2012 release):- Molecular modeling


MOE updated


The molecular modelling platform MOE 2014.09 has been updated, some of the new features in MOE include:

  • MOE Project for Organizing SBDD Data
  • Focused Protein and Antibody Libraries - Virtual Phage Display
  • Quantum Mechanical Refinement of Conformations and Energy Minimization
  • Template Forced Docking and Molecular Superposition
  • Non-natural Amino Acid Support for Protein and Peptide Design
  • Specialized Protein Family Databases and Search Interface

MOE is a software system designed to support Cheminformatics, Molecular Modelling, Bioinformatics, Virtual Screening, Structure-based-design and can be used to build new applications based on SVL (Scientific Vector Language).

There are several reviews of the previous versions of MOE here


Ocura (aka StarVue)


Previously known as StarVue, the latest release of Ocura, version 6.0, is now available. This latest release will enable you to open and view files from the latest releases of StarDrop 6.0 and Sentira 1.0.

Ocura is a desktop application is specifically designed for scientists who want a simple way to load a set of molecule structures and easily browse through the data.


POSIT updated


OpenEye ihave announce the release of POSIT v3.1, the component of the OEDocking suite devoted to pose prediction.

This update includes:

  • The HYBRID and FRED algorithms have been incorporated into POSIT, the appropriate method is determined by analyzing the ligand to pose against the input receptors.
  • Multiprocessing has been enabled through the use of MPI, to speed calculations.
  • POSIT now supports a list of receptors files or .lst file as input. This overcomes command-line limitations for the number of receptors that can be used simultaneously.
  • Added a MEDIOCRE result rating for results between 33% and 50% probability.
  • Command line parameters have been simplified and updated to be compatible with the OEDocking Suite of tools.

POSIT is designed for the posing problem in lead optimization, i.e. how best to leverage project information from previous protein-ligand structures to predict the pose of a new ligand. It does this by assessing the similarity of the new ligand to known bound structures. Performance degrades as similarity decreases and so at some point it is worth searching more exhaustively.


SeeSAR updated


SeeSAR 1.5 has been released. SeeSAR it is intended as an interactive tool for designing/improving ligands for drug discovery.

The latest release covers two major topics: 1. A series of features that make the editing more swift and easier. To this end they introduced hot-keys, context menus and drawing a bond by drag&drop. 2. Often times people use SeeSAR for visual inspection e.g. after docking. Now normally you'll have multiple poses per compound. For a better overview the Table now allows you to collapse all poses to just one line per compound.


Furthermore you can set a bookmark to indicate what you like and export only the ones on the wish list.

There is a review of SeeSAR here.


FORECASTER, docking tool


I just got this email

Thank you for your collaboration in helping us to test the beta version of the FORECASTER Suite 2014. From your feedback and bug reports, we have now released the final version of the Suite. The files were updated and posted on the download page. Please send us any bugs that you might have not yet reported.

The FITTED docking tool was initially been developed as a suite of three programs: SMART (used to prepare the small molecules for docking), PROCESS (used to prepare the protein files for docking) and the docking program FITTED. More recently, these three programs together with several others have been integrated into a single package, namely FORECASTER.

More information can be found here


2014 release of the ADF modeling suite


The ADF modeling suite consists of the GUI, the powerful DFT codes ADF (molecules) and BAND (surfaces, bulk), the semi-empirical DFTB and MOPAC2012 modules, ReaxFF, and COSMO-RS. The binaries for the entire suite work out of the box, in parallel, on all popular platforms (Windows, Mac, Linux).

A summary of new features and improvements in the 2014 release:

  • lower-memory, better parallel SCF in ADF
  • significant speed-ups in DFTB and the periodic DFT code BAND (AO-based)
  • many-body dispersion functionals (Tkatchenko et al.)
  • conformer search, support for multiple configurations, spectra averaging
  • TD-DFTB and DFTB-NEGF, with electronic parameters for 87 elements
  • ReaxFF force field optimizer, Grand Canonical Monte Carlo
  • COSMO-SAC 2013-ADF parameters

KiSTheIP added to alphabetical listing


KiSThelP is a cross-platform free open-source program developed to estimate molecular and reaction properties from electronic structure data. To date, three computational chemistry software formats are supported (Gaussian, GAMESS, NWChem). Some key features are:

  • gas-phase molecular thermodynamic properties (offering hindered rotor treatment)
  • thermal equilibrium constants
  • transition state theory rate coefficients (TST, VTST) including one-dimensional tunnelling effects (Wigner and Eckart)
  • RRKM rate constants, for elementary reactions with well-defined barriers.

KiSThelP is intended as a working tool both for the general public and also for more expert users. It provides graphical front-end capabilities designed to facilitate calculations and interpreting results. KiSThelP enables to change input data and simulation parameters directly through the GUI and to visually probe how it affects results. Users can access results in the form of graphs and tables. The graphical tool offers customizing of 2D-plots, exporting images and data files.


Sentira Review


Sentira is a new chemical data visualisation tool from Optibrium. The focus is on ease of use data visualisation and as such is probably targeted at the bench scientist rather than a specialist computational scientist. It supports a selection of plotting and SAR tools.

I’ve written a review of my first impressions.

There is also a list of data visualisation applications here.




Chargemol program performs atomic population analysis to determine DDEC net atomic charges, atomic spin moments, and effective bond orders. Because the DDEC net atomic charges are simultaneously optimized to reproduce atomic chemical states and the electrostatic potential surrounding a material, they are well-suited for constructing force-fields used in atomistic simulations (e.g., classical molecular dynamics or monte carlo simulations) and for quantifying electron transfer between atoms in complex materials and during chemical reactions

The DDEC method is described in the following publication and references therein.

Thomas A. Manz and David S. Sholl, "Improved Atoms-in-Molecule Charge Partitioning Functional for Simultaneously Reproducing the Electrostatic Potential and Chemical States in Periodic and Non-Periodic Materials", J. Chem. Theory Comput., Vol. 8 (2012) 2844-2867. DOI

The program can be run using either Matlab or Fortran source codes, which yield identical numbers. The Fortran code is parallelized with OpenMP and runs much faster than the Matlab code.


Quantum Expresso Updated


Quantum Expresso 5.1 is available for download.

QUANTUM ESPRESSO is an integrated suite of Open-Source computer codes for electronic-structure calculations and materials modelling at the nanoscale. It is based on density-functional theory, plane waves, and pseudopotentials.


New additions and updates


McQSAR:  A Multiconformational Quantitative Structure−Activity Relationship Engine Driven by Genetic Algorithms

McQSAR, an extension to the traditional GA approach to derive QSARs. McQSAR is able to use descriptors for multiple representations per compound, such as different conformers, tautomers, or protonation forms. Test runs show that the algorithm converges to a set of representations that describe the binding mode of the set of input molecules to a reasonable resolution provided that suitable descriptors based on the three-dimensional structure are used.

Mikko J. Vainio and Mark S. Johnson (2005) McQSAR: A Multiconformational Quantitative Structure-Activity Relationship Engine Driven by Genetic Algorithms. J. Chem. Inf. Model. 45, 1953-1961 DOI.

The recently updated Balloon creates 3D atomic coordinates from molecular connectivity via distance geometry and confomer ensembles using a multi-objective genetic algorithm. The input can be SMILES, SDF or MOL2 format. Output is SDF or MOL2. Flexibility of aliphatic rings and stereochemistry about double bonds and tetrahedral chiral atoms is handled.


Installing ACPC on a Mac


One of the advantages of using a Mac for science is that you can often make use of the UNIX underpinnings of Mac OSX to access programs written for Linux.

A recent publication in Journal of Cheminformatics caught my eye, screening of molecules using electrostatics is usually a very time-consuming process, but this publication describes an interesting and very quick way to screen molecules.

A rotation-translation invariant molecular descriptor of partial charges and its use in ligand-based virtual screening Francois Berenger, Arnout Voet, Xiao Yin Lee and Kam YJ Zhang Journal of Cheminformatics 2014, 6:23 doi

I’ve written instructions for how to install ACPC under Mac OSX.


A Review of Forge V10.2 on the New MacPro


Now that I have my new MacPro I thought it might be interesting to try out a couple of the software packages that I’ve previously reviewed. ForgeV10 allows the scientist to use Cresset’s proprietary electrostatic and physicochemical fields to align, score and compare diverse molecules. It allows the user to build field based pharmacophores to understand structure activity and then use the template to undertake a virtual screen to identify novel scaffolds. I’ve previously reviewed ForgeV10 and as it was formally known FieldAlign so I’m going to focus on the support for multiple processors and a few of the new features.

Read the review here


There is a compilation of software reviews here


Amber 14 and Amber 14 Tools released


Version 14 of the Amber software suite has been released (There was no "unlucky" Amber13.)

  • Force fields: Amber has two new fixed-charge protein force fields, ff14SB and ff14ipq, a new modular lipid force field, Lipid14, and updates to nucleic acid and carbohydrate force fields.
  • Improved options for self-guided Langevin dynamics and accelerated molecular dynamics, to enchance sampling along soft degrees of freedom.
  • A completely reorganized Reference Manual
  • QM/MM calculations can interface with a variety of external quantum chemistry programs, expanding the types of quantum models available
  • More features from sander have been added to pmemd for both CPU and GPU platforms, including performance improvements, and support for extra points, multi-dimension replica exchange, a Monte Carlo barostat, ScaledMD, Jarzynski sampling, explicit solvent constant pH, GBSA, and hydrogen mass repartitioning. Support is also included for the latest Kepler, Titan and GTX7xx GPUs.
  • Expanded methods are available for free energy calculations that change Hamiltonian models, including better procedures for appearing and disappearing atoms, and tighter integration with replica-exchange simulations, and a new absolute free energy method.
  • New facilities are present for using electron density maps (e.g. from cryo EM/ET experiments) as constraints, and to support rigid (or partially flexible) groups in simulations.

Amber Tools have also been updated.

Among the new features in AmberTools14:

  • The sander module, our workhorse simulation program, is now a part of AmberTools;
  • Greatly expanded and improved cpptraj program for analyzing trajectories;
  • new documentation and tools for inspecting and modifying Amber parameter files;
  • Improved workflow for setting up and analyzing simulations;
  • new capability for semi-empirical Born-Oppenheimer molecular dynamics;
  • EMIL: a new absolute free energy method using TI;
  • New Free Energy Workflow (FEW) tool automates free energy calculations (LIE, TI, and MM/PBSA-type calculations);
  • Completely reorganized Reference Manual

Porting of BUDE (Bristol University Docking Engine) to OpenCL.


A recently publication “High Performance in silico Virtual Drug Screening on Many-Core Processors” DOI describes porting BUDE (Bristol University Docking Engine) to OpenCL.

Our highly optimized OpenCL implementation of BUDE sustains 1.43 TFLOP/s on a single NVIDIA GTX 680 GPU, or 46% of peak performance. BUDE also exploits OpenCL to deliver effective performance portability across a broad spectrum of different computer architectures from different vendors, includ- ing GPUs from NVIDIA and AMD, Intel’s Xeon Phi and multi-core CPUs with SIMD instruction sets.

BUDE is now one the fastest HPC applications ever developed and nicely demonstrates the portability of OpenCL across different architectures.

There is a list of GPU accelerated applications here.


Academic access to CrunchYard HPC


CrunchYard is pleased to announce their academic offering for HPC in the cloud. The academic offering allows anyone with a valid academic e-mail to instantly access the online HPC facility.

Group accounts can be created for your research group, where simulation credits are shared. Credits are valid for 2 months after purchase.

Academic accounts include access to over 300 true CPU cores, and the following simulation software: FEKO CP2K NWCHEM GAMESS GROMACS LAMMPS

Other licensed codes are available on request (such as CPMD, STAR-CCM+ etc.)


Asteris and RCSB PDB updated


The official mobile app for viewing PDB structures RCSB PDB has been updated to include the April molecule of the month.

The molecule designer app Asteris has been updated to version 1.0.1 with a number of bug fixes and performance improvements.



GTKDynamo:- Molecular modeling extension for PyMOL


GTKDynamo is free/open source software which, together with pDynamo, transforms PyMOL into a powerful interface for molecular modeling. The interface has been designed to facilitate determining reaction pathways in biological systems, specially using hybrid QC/MM (or QM/MM) methods.

Some capabilities include:

  • Pure QC simulations - ab initio and SMO.
  • Pure MM simulations - using AMBER, CHARMM or OPLS force fields.
  • Hybrid QCMM simulations.
  • Single point calculations.
  • Energy minimization.
  • Molecular dynamics.
  • Reaction coordinate scanning.
  • Umbrella sampling.
  • Reaction path calculations - using NEB.

GTKDynamo is available for download for linux and Mac platforms . Please, make sure that you have installed:

  • pDynamo
  • Matplotlib
  • Numpy / Pylab
  • Pymol 1.x
  • pyGTK
  • ORCA, ab initio calculations.


J. F. R. Bachega, L. F. S. M. Timmers, L. Assirati, L. B. Bachega, M. J. Field, T. Wymore. J. Comput. Chem. 2013, 34, 2190-2196. DOI:


QMForge 2.3.2 released


QMForge is a program used to analyze the results of quantum chemistry (DFT) calculations. Gaussian 98/03/09, ADF, GAMESS (US), GAMESS (UK), PC-GAMESS, Jaguar, and ORCA files are supported. The following analyses are available:

  • The support of several QM output formats including Gaussian, ADF, GAMESS (and its various derivatives), Jaguar, Molpro, and ORCA,
  • Population analyses such as Mulliken, Lowdin, C-squared, and Overlap on user-defined sets ("fragments") of basis functions,
  • Gross Population Analysis to compare Mulliken and Lowdin populations and spin densities for each atom and orbital,
  • Fragment Analysis to interpret the contributions of fragment MOs to molecular MOs,
  • Frenking’s Charge Decomposition Analysis,
  • Calculation of Mayer's bond orders,
  • Visualization of each step in a geometry optimization, with the ability to save any of those structures as XYZ or PDB files,
  • Plots of Convergence and Energy in the Geometry Optimization tab,
  • A simple, yet powerful, XYZ editor with tools that allow translations and rotations to align bonds to specific axes,
  • Animation of the normal modes of a frequency calculation and the ability to save these as animated GIFs,
  • Plots of IR and/or Raman spectra in the Frequency Tab, and
  • A plot of electronic transitions and a easy-to-read list of the corresponding orbital excitations from TDDFT/CIS calculations.

QMForge has been created using the following dependencies:

  • The Python scripting language (2.7.5),
  • The Qt4 toolkit and its Python extensions PyQt4,
  • NumPy (1.7.1),
  • PyQwt (5.2.0), and
  • The cclib computational library (v1.2b).
  • OpenBabel and its Python bindings
  • The Python Imaging Library (PIL)
  • simplejson


QMForge has only been extensively tested on Mac OS X Snow Leopard and Mavericks

Tenderholt, Adam L. QMForge, Version 2.3.2,
Tenderholt, Adam L. "QMForge: A Program to Analyze Quantum Chemistry Calculations", Version 2.3.2,


Schrodinger Small Molecule Drug Discovery Suite Updated


The Schrodinger Small Molecule Drug Discovery Suite was updated over the weekend, this is a major update that brings in a host of new features and improvements.

Maestro Graphical Interface

Improved flexible ligand superposition Additional graphics settings
Real-time antialiasing Real-time ambient occlusion, outlines, and cartoon shading effects Multivariate ranking in the Project Table
Simultaneously maximize or minimize up to four property values, and rank entries based on the optimization Date Created and Date Modified fields automatically generated in the Project Table Workspace responsiveness of atom labels is up to 2.5x faster Click and drag to rearrange atom, measurement, and adjustment labels in the Workspace Support for bond labels Installed scripts and Tools menu items now searchable in the Task Tree Significant improvements to the Property Calculation interface in the project facility
Simultaneously calculate multiple properties Additional 2D properties now available: AlogP, #Hbond acceptors, #HBond donors, #rotatable bonds, polar surface area, molar refractivity, and polarizability

Ligand Docking

Ligand efficiencies are now calculated from the DockingScore instead of the GlideScore Generate per-residue interaction energies in Virtual Screening Workflow (VSW) for visualization New server mode in Glide Ligand Designer enables near real-time interactive docking (Glide Ligand Designer Script)

Pharmacophore Modeling

Performance improvements to Phase database operations, including faster deletion and insertion of ligands Automatic restart of Phase database subjobs

Field-Based QSAR

Use QM-calculated fields in 3D QSAR (command line only; phasefqsar script)
fqsar script generates Jaguar input files for computing QM electrostatic fields for use in 3D QSAR

Molecular Dynamics

Monitor secondary structure elements over the course of the trajectory (Simulation Interactions Diagram; SID)

Quantum Mechanics

New interface to compute thermodynamic properties for reactions New faster TDDFT algorithm and graphical interface Compute Raman intensities Several improvements to the results script Jaguar pKa displays the computed pKa as an atom label by default Heat of formation graphical interface now supports bromine and iodine Improved numerical stability of the 1st and 2nd derivatives of the D3 correction Increased utility of script
Script acts on a group of isomers and skips structures with unique stoichiometries

Protein X-Ray Refinement

Optionally set hydrogen B-values

Workflows & Pipelining

Includes the latest version of KNIME (v2.9)
Many new features including a Send Email node and ability to save workflows under different names; see for a complete list of new features Use any Glide simulation option in the Glide Ligand Docking node Employ a specific template in the Prime Build Homology Modeling node Import ungrouped structures to PyMOL from Run PyMOL node

Job Control

Improved fault tolerance Improved handling of suspended jobs in queueing systems

There are also updates to the Biologics Suite and the Materials Science Suite.


MolSoft Releases ICM Version 3.8


MolSoft have announced that a new version of ICM is now available for download from the support site . A description of the key new features can be found on the news page and release notes. To help you get to know the new features we will be holding a free webinar next week (2/11) - we hope you can join us, please register here.

Some of the key new features include:

  • Anaglyph Stereo
  • MolSkin - high quality surface graphics
  • Movies from Slides
  • SCARE - induced fit docking
  • Fragment screening
  • MolScreen - >360 high quality fingerprint and 3D pharmacophore models
  • Blast search direct from the GUI
  • ToxScore - new score for drug reactivity and toxicity

QMForge 2.3.1 released


QMForge 2.3.1, a cross-platform, open-source program for interpreting and analyzing the results of QM calculations has just been released

QMForge 2.3.1 builds upon the previous versions with the addition of the following features:

  • Plots of Convergence and Energy in the Geometry Optimization tab,
  • Plots of IR and/or Raman spectra in the Frequency Tab,
  • Ability to save the normal modes of a Frequency calculation as animated GIFs,
  • Gross Population Analysis to compare Mulliken and Lowdin populations and spin densities for each atom and orbital

Other notable features include

  • The support of several QM output formats including Gaussian, ADF, GAMESS (and its various derivatives), Jaguar, and ORCA,
  • Population analyses such as Mulliken, Lowdin, C-squared, and Overlap on user-defined sets ("fragments") of basis functions,
  • Fragment Analysis to interpret the contributions of fragment MOs to molecular MOs,
  • Charge Decomposition Analysis,
  • Mayer's bond orders,
  • Visualization of every step in a geometry optimization, with the ability to save any of those structures as XYZ or PDB files,
  • A simple, yet powerful, XYZ editor with tools that allow translations and rotations to align bonds to specific axes,
  • Animation of the normal modes of a frequency calculation, and
  • A plot of electronic transitions and a easy-to-read list of the corresponding orbital excitations from TDDFT/CIS calculations.

A review of FAst MEtabolizer (FAME)


Whilst much computational work is undertaken to support, library design, virtual screening, hit selection and affinity optimisation the reality is that the most challenging issues to resolve in drug discovery often revolve around absorption, distribution, metabolism and excretion (ADME). Whilst we can measure the levels of parent drug in various medium tracking metabolic fate can often be a considerably more difficult proposition requiring significant resources. For this reason prediction of sites of metabolism has become the subject of current interest.

FAME DOI is a collection of random forest models trained on a comprehensive and highly diverse data set of 20,000 small molecules annotated with their experimentally determined sites of metabolism taken from multiple species (rat, dog and human). In addition dedicated models are available to predict sites of metabolism of phase I and II processes.


FAME offers a high performance prediction of sites of metabolism mediated by a wide variety of mechanisms.

The full review is available here

There is a list of software reviews here.


MOE Update


Chemical Computing Group have just announced an update for MOE.

A patch update is now available for MOE 2013.08. This patch contains a series of important updates for better performance, s

In the MOE 2013.0801 patch update:

  • System Manager browsing speed-up for large systems
  • System Manager tag/group expansion only on Ctrl-click during browsing
  • Mac OS X NVIDIA graphics driver bug workaround
  • Mac OS X real-time ray tracing now supported by default on the new Intel Iris and Intel Iris Pro graphics cards
  • Bug fixes

Computational Chemistry Freeware


I recently came across this brilliant collection of software from Michel Petitjean

  • ARMS: Spatial Alignment with the RMS (Root Mean Square) method. (fixed pairwise correspondence)
  • ASV: Analytical calculation of van der Waals surfaces and volumes. (or any union of spheres)
  • CCCPP: Computes Cavites, Channels, Pores and Pockets in proteins.
  • CSR: The Combined SDM/RMS Algorithm for spatial alignment of two molecules. (pairwise correspondence computed)
  • CYL: Minimal radius enclosing cylinder. Minimal radius circumscribed cylinder.
  • DIVCF: Selects by clustering major conformations of a molecule in a set of its conformers.
  • DOG: Docking Geometrically two molecules. (fixed pairwise correspondence)
  • GRD: Computation of the Radius and Diameter of a molecular graph. (computes also the topological shape index)
  • MCG: Optimal Partition (classification): numerical variables and non-euclidean spaces. The number of classes is computed.
  • POP: Optimal Partition (classification): categorical variables. The number of classes is computed.
  • POSE: Computes the RMSD between two ligand poses. No rotation translation is performed.
  • QCM: Quantitative Chirality Measure of a conformer (graph automorphisms enumeration included)
  • RADI: Computation of the Radius and Diameter of a spatial set. (computes also various other geometrical parameters)
  • VIRAPOPS: A forward simulator dedicated to rapidly evolved viral populations.

Binaries are available for MacOSX and Linux.




I noticed that CPMD a parallelized plane wave / pseudopotential implementation of Density Functional Theory, particularly designed for ab-initio molecular dynamics. is now available on Crunchyard expanding the list of available Computational Chemistry packages. The following packages are also available. CP2K, LAMMPS, GAMESS, GROMACS , NWCHEM


CPMD added to alphabetical listing

I’ve added CPMD to the alphabetical listing of applications.

CPMD is an ab initio electronic structure and molecular dynamics (MD) program using a plane wave/pseudopotential implementation of density functional theory (DFT). It is mainly targeted at Car-Parrinello MD simulations, but also supports geometry optimizations, Born-Oppenheimer MD, path integral MD, response functions, QM/MM, excited states and calculation of some electronic properties.

Full installation instructions are available on the website together with the user manual and examples.

The examples can run effectively on an Intel Core i5 (2.53 GHz, OS/X) with 4 GB of available memory (most of tests require ca. 1GB, few more). Nonetheless, references have been generated on an IBM Blade power7 (with 8 mpi tasks and 1 OMP task per CPMD run) with 64 GB of available memory. Estimates of the execution time and memory requirements will be given based on this more performant setup.

CPMD capabilities

  • Works with norm conserving or ultrasoft pseudopotentials
  • LDA, LSD and the most popular gradient correction schemes; free energy density functional implementation
  • Isolated systems and system with periodic boundary conditions; k-points
  • Molecular and crystal symmetry
  • Wavefunction optimization: direct minimization and diagonalization
  • Geometry optimization: local optimization and simulated annealing
  • Molecular dynamics: constant energy, constant temperature and constant pressure
  • Path integral MD
  • Response functions
  • Excited states
  • Many electronic properties
  • Time-dependent DFT (excitations, molecular dynamics in excited states)
  • Coarse-grained non-Markovian metadynamics

CPMD is free for non-profit organisations.


SYBYL-X 2.1.1 released


SYBYL-X 2.1.1 is now available, the focus of this release is to extend the capabilities available via the standalone PYTHON interface to 3D-QSAR, which was introduced in SYBYL-X 2.1 earlier this year. The PYTHON API allows 3D-QSAR models (CoMFA, CoMSIA, and Topomer CoMFA) to be created and used for predictions outside of SYBYL-X.

  1. Hologram QSAR (HQSAR) is a now available via Python.  HQSAR has been successfully applied to generate predictive global QSAR models for on- and off-target effects and models for important ADME related properties; the HQSAR method employs 2D-substructural counts as descriptors.

  2. Similarity computations and similarity searches (UNITY 2D fingerprints) are now accessible via Python to support various workflows, such as lead expansion, lead hopping, and cluster analysis.


SZYBKI has been updated


OpenEye have just released an update to SZYBKI with a host of new features. SZYBKI is used to optimise the three dimensional structure of molecules prior to their use in other programs. SZYBKI also refines portions of a protein structure and optimize ligands within a protein active site, making it useful in conjunction with docking programs.

  • New utility program called FreeForm is available for Szybki users. It provides two distinct functionalities: evaluation of the solvation free energy of the input molecules and free energies of solution conformations. Please Note: FreeForm is not available under 32-bit Windows because of the high memory requirements of this application.
  • A new forcefield for protein-ligand interations is available upon selecting a new option -ff followed by AmberMMFF94 or AmberMMFF94S. It is a combination of MMFF94 (or MMFF94S) with Amber. In this combined force field, MMFF94 (or MMFF94S) is used to describe the intramolecular interactions of the ligand and the Amber force field is used for the VdW and Coulomb interactions between ligand and protein. Currently, this force field can be used only for ligands inside rigid proteins.
  • Entropy estimation based on analytical MMFF Hessian is extended for ligands bound in a rigid protein using the option -entropy.
  • Constraining torsion potential in the form: f8ebfd663514a377a56dfefaeb661230a663dc0f, where 81251f7be3a3c36e3c6602884accebe2d6ab3292 is the user specified force constant and 6b57c4a2f949d360120e9eec70ed44786c71b408 is the reference torsion dihedral angle is available with the use of the input flag -tor_constr.
  • New option -optMethod is introduced for optimization method selection. The possible choices of optimizers type are BFGS, conjugate gradient, steepest descent and mixtures of steepest descent preoptimization followed by BFGS of conjugent gradient. Option -conj is no longer supported.
  • Optimization of molecular systems with large number of degrees of freedom (>= 500) is by default done with the conjugate gradient method unless specifically requested with the option -optMethod.
  • New option -optGeometry is introduced for selection of coordinate system to be used during optimization. It replaces no longer supported options -optcart, -opttorsions, -solid and -noopt. It can be also used to optimize hydrogen atoms positions only, instead of using for that purpose -fixsmarts.
  • Molecules from the input file which failed during processing are by default written to the separate molecular file. See the description of new option -keepFailures.
  • Starting from this release flags -fixsmarts and -harmsmarts are followed by the name of the text file containing a single line with a SMARTS pattern used to fix or constrain atoms. Input of SMARTS strings on the command line for those flags are not used anymore because some SMARTS strings may contain special characters which may be interpreted incorrectly on some platforms.

MOE 2013.08 Released


Chemical Computing Group have just announced the release of the latest update of MOE (Molecular Operating Environment).

General Updates

Extended Hückel Theory for Pharmacophore Discovery

  • Apply EHT strengths to identify and discover weak and strong interactions
  • Annotate non-standard interactions: halogen and CH bond donors
  • Score pharmacophore hits using sum of interaction energies

de novo Loop and Linker Modeler

  • Search and browse for de novo and knowledge-based loop candidates
  • Generate multiple loop conformations and score loop-loop interactions • Investigate linkers for fusion proteins and dual variable domains

Protein Alignments and Superpositions

  • Superpose protein structures independent of sequence
  • Apply new threading methodology for sequence to structure alignments
  • Align sequences and superpose structures using STOVCA criteria

Core System Enhancements (with 64-bit support) The default versions of MOE for Linux, Windows, and Mac OSX are now all 64-bit. The current release also includes 32-bit versions for each architecture.

  • Create additional data grouping level in System Manager
  • Handle large number of protein:ligand complexes in real-time
  • Colour database viewer text and cells by data values for enhanced analysis, Empty cells can be specially colored

Interface to Mogul from CCDC

  • Access a knowledge-based library of small molecule crystallographic data
  • Visualize histograms and statistics fo rbond angles, dihedrals and torsions • Ensure ligand conformation is consistent with the CSD

Solvent Analysis using 3D-RISM Enhancements

  • Calculate water densities on full protein or protein:protein interface
  • Accurate placement and stability prediction of water sites
  • Use3 D-RISM densities to validate water position in crystal structures

chEMBL Library. MOE's SD Pipeline Command Tools were used to generate fragments from Release 14 of the ChEMBL database of bioactive drug-like small molecules. The resulting database of conformations, chemblr14_frag.mdb, for approximately 830,000 fragments, suitable for Scaffold Replacement or combinatorial chemistry methodologies

Mac OS X Enhancements.

  1. OS X Dock Icon Enhancements. A new Dock menu, accessed by right-clicking (or control-clicking) the MOE Dock icon, allows for standard Mac OS X operations – Hide, Show, and Quit – across all open iterations of MOE. The Dock Menu also contains menu items in order for users to view the About panel bundled in MOE, as well as access the Preferences panel. Additionally, the Dock indicator light under the MOE Dock icon remains persistent after having launched a new instance of MOE, until all opened MOE instances have been closed. Users continue to be able to drag and drop compatible MOE documents (i.e. .pdb, .svl, etc.) onto the Desktop and/or Dock icons for direct opening in MOE.

  2. New Preferences Menu. A new Preferences panel allows for certain parameters to be set and/or changed graphically. These parameters include the defaults key and path to MOE, the current configuration of MOE (32- or 64-bit, or Auto-select), as well as a selection of common command line tags, with a custom section for inputting other command line options.

  3. Output Continuity. stdout and stderr received from MOE have been standardized to be displayed in a native panel in a scrollable text box.


Scripting Vortex 16


OCHEM is a free open access site of annotated models and chemical data. OCHEM contains 1831772 experimental records for about 477 properties collected from 12457 sources you are free to upload your own data and also build predictive models using existing or your own data.

There are also a number of already built models that the public can access, these include

  • Ames test
  • CYP1A2 inhibition
  • LogP and Solubility

You can run predictions on OCHEM using simple REST-like web services, these vortex scripts submit tasks to the various models and then retrieve the resulting prediction.


Turning up the heat with Torch and Forge


Anyone involved in a drug discovery programme will be aware of the challenge presented by trying to visualise and explore structure-activity relationships (SAR), in particular visualising questions like :-

“What is the largest change that can be made whilst maintaining activity?”

Activity Miner from Cresset is a new tool designed to rapidly interrogate and decipher SAR in both Torch and Forge. Activity Miner is intended to help identify key elements of the SAR by starting from a set of aligned molecules and then automatically comparing them to each other.

Activity Miner1

More details are here


ROCS Updated


OpenEye have just announced that the virtual screening tool ROCS v3.2 has been released.

Several noteworthy features have been added to this version including a -subrocs option that can drastically improve substructure alignments. Also included is an application rocs-report that uses our 2D depiction technology to make pdf reports of hitlists displayed with 2D similarity, shape and color overlaps, as well as property histograms. Substantial upgrades have been made to vROCS. An improved sketcher now highlights unspecified stereochemistry in atoms and bonds in query structures, and requires the user to correct any unspecified stereochemistry.

ROCS is available for download here.


StarDrop 5.4

StarDrop was recently updated to version 5.4, this brings an update to the virtual library design module and scaffold based design, there have also been improvements to the plotting and data visualisation.

There are now seven optional plugins with three exciting new options.

Derek Nexus™ - Knowledge based toxicity prediction The new Derek Nexus module for StarDrop provides Lhasa Limited's world-leading technology for knowledge-based prediction of key toxicities. Using data from published and donated (unpublished) sources, Derek Nexus identifies structure-toxicity relationships that alert you to the potential for your compounds to cause toxicity. The Derek Nexus module provides predictions of the likelihood of a compound causing toxicity in over 40 endpoints, including mutagenicity, hepatotoxicity and cardiotoxicity.

BIOSTER™ - A world of chemistry experience BIOSTER is developed and updated in collaboration with Digital Chemistry and is available as an optional extension to StarDrop's Nova module. This combination enables you to quickly and easily search the comprehensive BIOSTER database to identify transformations that are relevant to your compounds. These can be automatically applied to generate novel structures with a high likelihood of biological activity and synthetic accessibility, prioritised against the property profile you require for your project. BIOSTER brings the collective experience of the chemistry community to help you to discover new active analogues of your compounds based on the tried and tested principle of isosterism. The BIOSTER module contains a unique compilation of over 20,000 precedented bioisosteric transformations, manually curated from the literature by Dr István Ujváry, complete with references to the original publications in which they are described.

torch3D™ The renamed torch3D module, using Cresset’s unique Field technology to understand and apply 3D Structure Activity Relationship (SAR), has been updated to include the latest version of Cresset’s XED force field providing insight into compounds’ 3D structures, biological activities and interactions.

These certainly significantly expand the potential utility of StarDrop, but note that these are not part of the standard install and may require additional licensing.