Geoffrey R. Hutchison

Open Babel: An open chemical toolbox

BackgroundA frequent problem in computational modeling is the interconversion of chemical structures between different formats. While standard interchange formats exist (for example, Chemical Markup Language) and de facto standards have arisen (for example, SMILES format), the need to interconvert formats is a continuing problem due to the multitude of different application areas for chemistry data, differences in the data stored by different formats (0D versus 3D, for example), and competition between software along with a lack of vendor-neutral formats.ResultsWe discuss, for the first time, Open Babel, an open-source chemical toolbox that speaks the many languages of chemical data. Open Babel version 2.3 interconverts over 110 formats. The need to represent such a wide variety of chemical and molecular data requires a library that implements a wide range of cheminformatics algorithms, from partial charge assignment and aromaticity detection, to bond order perception and canonicalization. We detail the implementation of Open Babel, describe key advances in the 2.3 release, and outline a variety of uses both in terms of software products and scientific research, including applications far beyond simple format interconversion.ConclusionsOpen Babel presents a solution to the proliferation of multiple chemical file formats. In addition, it provides a variety of useful utilities from conformer searching and 2D depiction, to filtering, batch conversion, and substructure and similarity searching. For developers, it can be used as a programming library to handle chemical data in areas such as organic chemistry, drug design, materials science, and computational chemistry. It is freely available under an open-source license from http://openbabel.org.

Avogadro: an advanced semantic chemical editor, visualization, and analysis platform

BackgroundThe Avogadro project has developed an advanced molecule editor and visualizer designed for cross-platform use in computational chemistry, molecular modeling, bioinformatics, materials science, and related areas. It offers flexible, high quality rendering, and a powerful plugin architecture. Typical uses include building molecular structures, formatting input files, and analyzing output of a wide variety of computational chemistry packages. By using the CML file format as its native document type, Avogadro seeks to enhance the semantic accessibility of chemical data types.ResultsThe work presented here details the Avogadro library, which is a framework providing a code library and application programming interface (API) with three-dimensional visualization capabilities; and has direct applications to research and education in the fields of chemistry, physics, materials science, and biology. The Avogadro application provides a rich graphical interface using dynamically loaded plugins through the library itself. The application and library can each be extended by implementing a plugin module in C++ or Python to explore different visualization techniques, build/manipulate molecular structures, and interact with other programs. We describe some example extensions, one which uses a genetic algorithm to find stable crystal structures, and one which interfaces with the PackMol program to create packed, solvated structures for molecular dynamics simulations. The 1.0 release series of Avogadro is the main focus of the results discussed here.ConclusionsAvogadro offers a semantic chemical builder and platform for visualization and analysis. For users, it offers an easy-to-use builder, integrated support for downloading from common databases such as PubChem and the Protein Data Bank, extracting chemical data from a wide variety of formats, including computational chemistry output, and native, semantic support for the CML file format. For developers, it can be easily extended via a powerful plugin mechanism to support new features in organic chemistry, inorganic complexes, drug design, materials, biomolecules, and simulations. Avogadro is freely available under an open-source license from http://avogadro.openmolecules.net.

The Blue Obelisk-interoperability in chemical informatics.

The Blue Obelisk Movement (http://www.blueobelisk.org/) is the name used by a diverse Internet group promoting reusable chemistry via open source software development, consistent and complimentary chemoinformatics research, open data, and open standards. We outline recent examples of cooperation in the Blue Obelisk group:  a shared dictionary of algorithms and implementations in chemoinformatics algorithms drawing from our various software projects; a shared repository of chemoinformatics data including elemental properties, atomic radii, isotopes, atom typing rules, and so forth; and Web services for the platform-independent use of chemoinformatics programs.

Pybel: a Python wrapper for the OpenBabel cheminformatics toolkit

BackgroundScripting languages such as Python are ideally suited to common programming tasks in cheminformatics such as data analysis and parsing information from files. However, for reasons of efficiency, cheminformatics toolkits such as the OpenBabel toolkit are often implemented in compiled languages such as C++. We describe Pybel, a Python module that provides access to the OpenBabel toolkit.ResultsPybel wraps the direct toolkit bindings to simplify common tasks such as reading and writing molecular files and calculating fingerprints. Extensive use is made of Python iterators to simplify loops such as that over all the molecules in a file. A Pybel Molecule can be easily interconverted to an OpenBabel OBMol to access those methods or attributes not wrapped by Pybel.ConclusionPybel allows cheminformaticians to rapidly develop Python scripts that manipulate chemical information. It is open source, available cross-platform, and offers the power of the OpenBabel toolkit to Python programmers.

De novo design of molecular wires with optimal properties for solar energy conversion

The area of organic photovoltaic materials has elicited great interest in both the scientific and technological communities due to its potential to deliver cheap and highly efficient solar cells [1]. To date, however, such so-called molecular wires have typically yielded energy conversion efficiencies of only ~5-6% despite a theoretical maximum of 13% [2]. We present an approach that uses a genetic algorithm to search the space of synthetically accessible molecular wires for those with optimal electronic structures. This approach combines both cheminformatics (SMILES to 3D using OpenBabel) and computational chemistry (semi-empirical calculations using Gaussian09). Using this method, we have found hundreds of candidates with predicted efficiencies over 8% including many with efficiencies over 10%.

Tuning the Kondo Effect with a Mechanically Controllable Break Junction

We study electron transport through C(60) molecules in the Kondo regime using a mechanically controllable break junction. By varying the electrode spacing, we are able to change both the width and the height of the Kondo resonance, indicating modification of the Kondo temperature and the relative strength of coupling to the two electrodes. The linear conductance as a function of T/T(K) agrees with the scaling function expected for the spin-1/2 Kondo problem. We are also able to tune finite-bias Kondo features which appear at the energy of the first C(60) intracage vibrational mode.

Open Data, Open Source and Open Standards in chemistry: The Blue Obelisk five years on

BackgroundThe Blue Obelisk movement was established in 2005 as a response to the lack of Open Data, Open Standards and Open Source (ODOSOS) in chemistry. It aims to make it easier to carry out chemistry research by promoting interoperability between chemistry software, encouraging cooperation between Open Source developers, and developing community resources and Open Standards.ResultsThis contribution looks back on the work carried out by the Blue Obelisk in the past 5 years and surveys progress and remaining challenges in the areas of Open Data, Open Standards, and Open Source in chemistry.ConclusionsWe show that the Blue Obelisk has been very successful in bringing together researchers and developers with common interests in ODOSOS, leading to development of many useful resources freely available to the chemistry community.

Cinfony – combining Open Source cheminformatics toolkits behind a common interface

BackgroundOpen Source cheminformatics toolkits such as OpenBabel, the CDK and the RDKit share the same core functionality but support different sets of file formats and forcefields, and calculate different fingerprints and descriptors. Despite their complementary features, using these toolkits in the same program is difficult as they are implemented in different languages (C++ versus Java), have different underlying chemical models and have different application programming interfaces (APIs).ResultsWe describe Cinfony, a Python module that presents a common interface to all three of these toolkits, allowing the user to easily combine methods and results from any of the toolkits. In general, the run time of the Cinfony modules is almost as fast as accessing the underlying toolkits directly from C++ or Java, but Cinfony makes it much easier to carry out common tasks in cheminformatics such as reading file formats and calculating descriptors.ConclusionBy providing a simplified interface and improving interoperability, Cinfony makes it easy to combine complementary features of OpenBabel, the CDK and the RDKit.

Turing patterns in a self-replicating mechanism with a self-complementary template

A variety of nonlinear chemical models, such as the Selkov–Schnakenberg, exhibit Turing patterns. The Templator, which is based on a minimal autocatalytic monomer–dimer system, is a newer two-variable model also able to show Turing patterns. Here we find that the dynamic behavior of the Templator is quite similar to other models with cubic nonlinearities. This is demonstrated through a series of computer simulations in two dimensions utilizing the cellular automata approach. The selection of parameter values is based on linear stability analysis, which provides a relatively simple method of predicting Turing pattern formation. The simulations reveal spot, labyrinth, and striped patterns, in agreement with the predictions of the analysis. Other behaviors, such as honeycomb patterns, are also observed. For some parameter values, we study transient spot replication. Our findings strongly suggest that the Templator may belong to the same class of models previously studied by Pearson.

Piezoelectric Effects of Applied Electric Fields on Hydrogen-Bond Interactions: First-Principles Electronic Structure Investigation of Weak Electrostatic Interactions

The piezoelectric properties of 2-methyl-4-nitroaniline crystals were explored qualitatively and quantitatively using an electrostatically embedded many-body (EE-MB) expansion scheme for the correlation energies of a system of monomers within the crystal. The results demonstrate that hydrogen bonding is an inherently piezoelectric interaction, deforming in response to the electrostatic environment. We obtain piezo-coefficients in excellent agreement with the experimental values. This approach reduces computational cost and reproduces the total resolution of the identity (RI)-Møller-Plesset second-order perturbation theory (RI-MP2) energy for the system to within 1.3 × 10(-5)%. Furthermore, the results suggest novel ways to self-assemble piezoelectric solids and suggest that accurate treatment of hydrogen bonds requires precise electrostatic evaluation. Considering the ubiquity of hydrogen bonds across chemistry, materials, and biology, a new electromechanical view of these interactions is required.