Home > Research > Computation & Modelling

Computation & Modelling

Parallel Computing

A major activity at SCRU is the porting of popular scientific software used in chemistry, biophysics and materials science to Graphical Processing Units (GPU). The SCRU laboratories have developed a Quantum Supercharger Library capable of accelerating legacy codes such as GAMESS-UK, GAMESS-US and NWChem on hardware accelerators such as GPGPUs. Kevin Naidoo who directs SCRU research has been collaborating with Quantum code developers, Paul Sherwood and Martyn Guest in the UK to accelerate the GAMESS-UK code on a GPU cluster.


Kyle D. Fernandes, C. Alicia Renison and Kevin J. Naidoo. Quantum Supercharger Library: Hyper-Parallelism of the Hartree Fock Method. J. Comp. Chem. 2015, 36, 1399-1409.

C. Alicia Renison, Kyle D. Fernandes and Kevin J. Naidoo. Quantum Supercharger Library: Hyper-Parallel Integral Derivatives Algorithms for Ab Initio QM/MM Dynamics. J. Comp. Chem. 2015, 36, 1410-1419

Wilkinson, K.; Sherwood, P.; Guest, M. F.; Naidoo, K. J., Acceleration of the GAMESS-UK Electronic Structure Package on Graphical Processing Units. J. Comp. Chem. 2011, 32, 2313-2318.

Kevin J. Naidoo*, Karl Wilkinson and Kyle Fernandes, Accelerating Scientific Computing Code in FORTRAN: The Quantum Chemistry Project; PGI News October 2011 (

QM/MM Computations

We develop methods for chemical biological applications. Our laboratory designed a new method to cut polar (glycosidic bonds) making it possible to accurately compute complex glycans using a mixed quantum classical model. We call this the Simple Link Atom Hybrid Saccharide (SLASH) method. We have developed semi-empirical quantum methods to model chemical glycobiological events that we termed AM1/d-CB1. We are able to simulate glycosylation reactions more accurately using these methods that are specifically designed for chemical glycobiological applications.


Werner Crous, Martin J. Field, and Kevin J. Naidoo, Simple Link Atom Saccharide Hybrid (SLASH) Treatment for Glycosidic Bonds at the QM/MM Boundary, J. Chem. Theor. Comput. 2014 10 (4), 1727-1738

Krishna Govender, Jiali Gao and Kevin J. Naidoo. AM1/d-CB1: A Semiempirical Model for QM/MM Simulations of Chemical Glycobiology Systems. J. Chem. Theor. Comput. 2014, 10, 4694-4707

Barnett, C. B.; Naidoo, K. J., Ring puckering: A metric for evaluating the accuracy of AM1, PM3, PM3CARB-1 and SCC-DFTB carbohydrate QM/MM simulations. J. Phys. Chem. B. 2010, 114, 17142-17154

Free Energy Computations

We take the view that the Free Energy of a system is a highly informative means with which to study it and design and improve on it. To this end have developed a method capable of calculating the Free Energy from Adaptive Reaction Coordinate Forces (FEARCF) for molecular conformation, association, and chemical reactions{Rajamani, 2003 #1666} in multiple dimensions.


Naidoo, K. J.; Brady, J. W., Calculation of the Ramachandran Potential of Mean Force for a Disaccharide in Aqueous Solution. J. Am. Chem. Soc. 1999, 121, 2244-2252.

Rajamani, R.; Naidoo, K. J.; Gao, J., Implementation of an Adaptive Umbrella Sampling Method for the Calculation of Multidimensional Potential of Mean Force of Chemical Reactions in Solution. J. Comp. Chem. 2003, 24, 1775–1781.

Naidoo, K. J., Multidimensional free energy volumes offer unique insights into reaction mechanisms, molecular conformation and association. Phys. Chem. Chem. Phys. 2012, 14, 9026-9036.

Strümpfer, J.; Naidoo, K. J., Computing free energy hypersurfaces for anisotropic intermolecular associations. J. Comp. Chem. 2010, 31, 308-316.

Coarse Grain Computations

Now that the human genome has been decoded the next step is to map out the complete set of protein structures that are translated from the genome. This is the field of proteomics. The major problem is to quickly solve 3D protein structures this is largely a computational problem. We develop tools to use a coarse grained potential functions embedded in a molecular dynamics calculations to fold and refine protein and glycan conformations. Our specific structural interests are glycosyltransferases, glycosidases and glucokinases.


Naidoo, K. J.; Schnitker, J.; Weeks, J. D., Two-Dimensional Melting Revisited: Molecular Dynamics Simulations Initiated with Optical Microscopy Data., Mol. Phys. 1993, 80, 1-43.

Naidoo, K. J.; Schnitker, J., Melting of two-dimensional colloidal crystals: A simulation study of the Yukawa system. J. Chem. Phys. 1994, 100, 3114-3121.

Molecular Modelling

We apply the computational methods that we develop along with standard molecular dynamics and hybrid quantum classical simulation techniques to two classes of molecules. Firstly, we discover glycosidase reaction mechanisms, from which we develop protocols for rational drug design from carbohydrate scaffolds. A second class of molecules are those of solvents making up the condensed phase. Here we investigate solution properties and design Ionic Liquid properties.


Matthews, R. P.; Venter, G. A.; Naidoo, K. J., Using Solvent Binding and Dielectric Friction to Interpret the Hydration Behavior of Complex Anions. J. Phys. Chem. B 2011, 115, 1045–1055.

Boscaino, A.; Naidoo, K. J., The Extent of Conformational Rigidity Determines Hydration in Non Aromatic HexaCyclic Systems. J. Phys. Chem. B 2011, 115, 2608–2616.

Barnett, C. B.; Wilkinson, K. A.; Naidoo, K. J., Molecular Details from Computational Reaction Dynamics for the Cellobiohydrolase I Glycosylation Reaction. J. Am. Chem. Soc. 2011, 133, 19474-19482.