Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Electronic properties of 4,4 ',5,5 '-tetramethyl-2,2 '-biphosphinine (tmbp) in the redox series fac-[Mn(Br)(CO)(3)(tmbp)], [Mn(CO)(3)(tmbp)](2), and [Mn(CO)(3)(tmbp)](-): crystallographic, spectroelectrochemical, and DFT computational study

Stepwise electrochemical reduction of the complex fac-[Mn(Br)(CO)(3)(tmbp)] (tmbp = 4,4',5,5'-tetramethyl-2,2'-biphosphinine) produces the dimer [Mn(CO)(3)(tmbp)](2) and the five-coordinate anion [Mn(CO)(3)(tmbp)](-). All three members of the redox series have been characterized by single-crystal X-ray diffraction. The crystallographic data provide valuable insight into the localization of the added electrons on the (carbonyl)manganese and tmbp centers. In particular, the formulation of the two-electron-reduced anion as [Mn-0(CO)(3)(tmbp(-))](-) also agrees with the analysis of its IR nu(CO) wavenumbers and with the results of density functional theoretical (DFT) MO calculations on this compound. The strongly delocalized pi-bonding in the anion stabilizes its five-coordinate geometry and results in the appearance of several mixed Mn-to-tmbp charge-transfer/IL(tmbp) transitions in the near-UV-vis spectral region. A thorough voltammetric and UV-vis/IR spectroelectrochemical study of the reduction path provided evidence for a direct formation of [Mn(CO)(3)(tmbp)](-) via a two-electron ECE mechanism involving the [Mn(CO)(3)(tmbp)](.) radical transient. At ambient temperature [Mn(CO)(3)(tmbp)](-) reacts rapidly with nonreduced fac-[Mn(Br)(CO)(3)(tmbp)] to produce [Mn(CO)(3)(tmbp)](2). Comparison with the analogous 2,2'-bipyridine complexes has revealed striking similarity in the bonding properties and reactivity, despite the stronger pi-acceptor character of the tmbp ligand.

Titanocene anticancer complexes and their binding mode of action to human serum albumin: a computational study

Due to the pivotal role played by human serum albumin (HSA) in the transport and cytotoxicity of titanocene complexes, a docking study has been performed on a selected set of titanocene complexes to aid in the current understanding of the potential mode of action of these titanocenes upon binding HSA. Analysis of the docking results has revealed potential binding at the known drug binding sites in HSA and has provided some explanation for the specificity and subsequent cytotoxicity of these titanocenes. Additionally, a new alternative binding site for these titanocenes has been postulated.

Computational modes and grid imprinting on five quasi-uniform spherical C-grids

Currently, most operational forecasting models use latitude-longitude grids, whose convergence of meridians towards the poles limits parallel scaling. Quasi-uniform grids might avoid this limitation. Thuburn et al, JCP, 2009 and Ringler et al, JCP, 2010 have developed a method for arbitrarily-structured, orthogonal C-grids (TRiSK), which has many of the desirable properties of the C-grid on latitude-longitude grids but which works on a variety of quasi-uniform grids. Here, five quasi-uniform, orthogonal grids of the sphere are investigated using TRiSK to solve the shallow-water equations. We demonstrate some of the advantages and disadvantages of the hexagonal and triangular icosahedra, a Voronoi-ised cubed sphere, a Voronoi-ised skipped latitude-longitude grid and a grid of kites in comparison to a full latitude-longitude grid. We will show that the hexagonal-icosahedron gives the most accurate results (for least computational cost). All of the grids suffer from spurious computational modes; this is especially true of the kite grid, despite it having exactly twice as many velocity degrees of freedom as height degrees of freedom. However, the computational modes are easiest to control on the hexagonal icosahedron since they consist of vorticity oscillations on the dual grid which can be controlled using a diffusive advection scheme for potential vorticity.