Exploring the selectivity of cytochrome P450 for enhanced novel anticancer agent synthesis.

Cytochrome P450 monooxygenases are an extensive and unique class of enzymes, which can regio- and stereo-selectively functionalise hydrocarbons by way of oxidation reactions. These enzymes are naturally occurring but have also been extensively applied in a synthesis context, where they are used as efficient biocatalysts. Recently, a biosynthetic pathway where a cytochrome P450 monooxygenase catalyses a critical step of the pathway was uncovered, leading to the production of a number of products that display high antitumour potency. In this work, we use computational techniques to gain insight into the factors that determine the relative yields of the different products. We use conformational search algorithms to understand the substrate stereochemistry. On a machine-learned 3D protein structure, we use molecular docking to obtain a library of favourable poses for substrate-protein interaction. With molecular dynamics, we investigate the most favourable poses for reactivity on a molecular level, allowing us to investigate which protein-substrate interactions favour a given product and thus gain insight into the product selectivity.

Gaseous inhibition of the transsulfuration pathway by cystathionine ß-synthase.

The transsulfuration pathway plays a key role in mammals for maintaining the balance between cysteine and homocysteine, whose concentrations are critical in several biochemical processes. Human cystathionine ß-synthase is a heme-containing, pyridoxal 5'-phosphate (PLP)-dependent enzyme found in this pathway. The heme group does not participate directly in catalysis, but has a regulatory function, whereby CO or NO binding inhibits the PLP-dependent reactions. In this study, we explore the detailed structural changes responsible for inhibition using quantum chemical calculations to validate the experimentally observed bonding patterns associated with heme CO and NO binding and molecular dynamics simulations to explore the medium-range structural changes triggered by gas binding and propagating to the PLP active site, which is more than 20 Å distant from the heme group. Our results support a previously proposed mechanical signaling model, whereby the cysteine decoordination associated with gas ligand binding leads to breaking of a hydrogen bond with an arginine residue on a neighbouring helix. In turn, this leads to a shift in position of the helix, and hence also of the PLP cofactor, ultimately disrupting a key hydrogen bond that stabilizes the PLP in its catalytically active form.

Exploration of biochemical reactivity with a QM/MM growing string method.

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

The Mechanism of Biochemical NO-Sensing: Insights from Computational Chemistry.

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO-receptor in humans is soluble guanylyl cyclase (sGC) - a histidine-ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO-signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high-level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H-NOX protein from Shewanella oneidensis (So H-NOX) homologous to the NO-binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six-coordinate intermediate heme-NO complex, followed by relatively facile His decoordination yielding a five-coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.

Soil biotransformation of thiodiglycol, the hydrolysis product of mustard gas: understanding the factors governing remediation of mustard gas contaminated soil.

Thiodiglycol (TDG) is both the precursor for chemical synthesis of mustard gas and the product of mustard gas hydrolysis. TDG can also react with intermediates of mustard gas degradation to form more toxic and/or persistent aggregates, or reverse the pathway of mustard gas degradation. The persistence of TDG have been observed in soils and in the groundwater at sites contaminated by mustard gas 60 years ago. The biotransformation of TDG has been demonstrated in three soils not previously exposed to the chemical. TDG biotransformation occurred via the oxidative pathway with an optimum rate at pH 8.25. In contrast with bacteria isolated from historically contaminated soil, which could degrade TDG individually, a consortium of three bacterial strains isolated from the soil never contaminated by mustard gas was able to grow on TDG in minimal medium and in hydrolysate derived from an historical mustard gas bomb. Exposure to TDG had little impacts on the soil microbial physiology or on community structure. Therefore, the persistency of TDG in soils historically contaminated by mustard gas might be attributed to the toxicity of mustard gas to microorganisms and the impact to soil chemistry during the hydrolysis. TDG biodegradation may form part of a remediation strategy for mustard gas contaminated sites, and may be enhanced by pH adjustment and aeration.