Engineering OleT Enzyme to Improve Biofuel Yield

Introduction

Global energy demand is reported to increase significantly over the last one year. The Global Energy and CO2 Status Report that was released in March 2018 indicated that energy requirements all over the world increased by a margin of 2.1% in 2017, which was more than double the requirements of the previous year (International Energy Agency 2018). This increase was attributed to a boost in economic growth all over the world. Approximately 70% of the energy demands were met by gas, coal and oil, whereas renewable energy sources accounted for the rest of the demand. The rapid increase in global fuel consumption has led to many concerns regarding dwindling fuel reserves.

This problem is made worse by volatile geopolitical factors, the inconsistent cost of fossil fuel and increasing apprehensions about national energy security (Covert et al. 2016). Moreover, the use of fossil fuel causes serious environmental concerns due to the emission of greenhouse gases (Nicoletti et al. 2015). The International Energy Agency (2018) reports that the increase in global energy use was associated with the emission of vast quantities of carbon dioxide gas, which is a known greenhouse gas. Carbon dioxide gas emissions increased by 1.7%, whereas carbon emissions rose to 32.5 gigatonnes, which was the highest levels ever recorded.

Thus, these findings show that fossil fuel substitutes should be developed to protect the environment from the adverse effects of large quantities of carbon dioxide gas. Biological materials can be used to generate renewable sources of energy. Additionally, fuels from these sources do not cause adverse environmental consequences. Biofuels have been around for many years. Peanut oil was used to operate diesel engines in the early 20th century (Guo, Song & Buhain 2015).

However, the discovery of cheaper fossil fuels shifted attention to their use thereby causing biofuels to be forgotten. Bioethanol and biodiesel are the most commonly used biofuels in the current market (Pulidindi & Gedanken 2015). Bioethanol is made by the fermentation of carbohydrates in a process that is similar to that of beer brewing. In recent years, bioethanol has been used as an additive to fossil fuels to minimise the production of carbon monoxide and other emissions that lead to the formation of smog (Shahir et al. 2015; Mofijur et al. 2016). For example, in the United States, ethanol from corn is added to gasoline (Chum et al. 2016).

Additionally, flexible-fuel vehicles that can operate on mixtures of gasoline and high proportions of ethanol are now available. Countries such as Brazil have been using ethanol from sugarcane as a biofuel and have cars whose engines can be powered using pure ethanol (Chum et al. 2016; Martins et al. 2018). Similarly, Europe produces biodiesel from palm oil for use in motor vehicle engines (Mathews & Ardiyanto 2015). Biodiesel is produced by mixing alcohol (in most cases methanol) with organic substances such as animal fat, and vegetable oil. The most common use of biodiesel is to decrease vehicle emissions. However, it can still be used as a renewable substitute fuel in its pure form in engines that use diesel.

Other forms of biofuel encompass methanol and reformulated gasoline constituents. Methanol, which is also referred to as wood alcohol, is mainly manufactured from methane or biomass. Gasification, which is the vaporization of biomass at elevated temperatures followed by the elimination of impurities and passing of the gas over a catalyst, is the most common approach that is used to produce natural gas. Examples of reformulated gasoline constituents include methyl tertiary butyl ether and ethyl tertiary butyl ether. These substitutes produce reduced amounts of pollutants as opposed to conventional fuels.

Nevertheless, the best biofuels are bio-hydrocarbons, particularly intermediate to long-chain fatty alkanes or alkenes (Liu et al. 2014). Aliphatic hydrocarbons also serve other purposes such as polymers, detergents and lubricants in industrial and domestic applications. These biofuels are preferred because of their close semblance to petroleum-based fuels with respect to their chemical structures and physical attributes. Therefore, recent research efforts have focused on the biosynthetic pathways of biofuels from different microorganisms as well as their potential in the production of aliphatic hydrocarbons (Liao et al. 2016; Zargar et al. 2017). One benefit of using microorganisms to synthesize fatty alkanes or alkenes is that microbial systems are can be scaled appropriately based on the quantities of the desired products. Consequently, microbial systems are cost-effective.

The production of biofuels in microbial systems is a biological process that involves biological reactions, which require enzymes. Enzymes are biological catalysts that speed up the rate of biochemical reactions. Living organisms produce enzymes naturally to expedite the conversion of substrates to useful products that play various roles in different physiological functions. A typical enzyme reaction involves the attachment of the substrate to the reaction centre of the enzyme, which results in the creation of an enzyme-substrate complex (Cornish-Bowden 2014).

Thereafter, an intricate series of chemical reactions lead to the rearrangement of chemical bonds in the substrate, which leads to the formation of a product. Naturally occurring enzymes are highly efficient in their various roles at biological temperatures and pressures by reducing the rate of reactions substantially. Then again, experimental conditions differ from physiological ones. These differences interfere with the catalysis of the enzyme outside biological systems.

As a result, rigorous research is in progress to elucidate the properties and machinery of biological enzymes to facilitate the exploitation of biocatalysts in the chemical industry. However, such endeavours have been frustrated by the nature of enzyme-catalysed reactions, which are too fast for investigators to obtain accurate details that would be beneficial for comprehending the machinery and reaction pathways. Therefore, novel computational techniques have been developed leading to the creation of accurate models that can envisage the reactions of biocatalysts. Thus, it is now possible to examine large biochemical systems.

Investigations concerning the production of biofuel in microbial systems have shown the involvement of a special group of enzymes known as cytochrome P450 alongside heme enzymes. Cytochrome P450 is an important group of enzymes whose main purpose is to metabolise drugs in the human body (Liu et al. 2014). These enzymes are ubiquitous in living organisms, including plants, fungi, bacteria, and eukaryotic organisms (Nelson 2018).

Different subclasses of cytochrome P450 have been characterized based on their involvement in different biological functions. A heme enzyme reacts with fatty acid substrates in a reaction that involves a decarboxylation step mediated by a cytochrome P450 peroxygenase. This scheme has the capacity to synthesise biofuel (Grant et al. 2015). However, it liberates several by-products that must be removed from the desired product before it can be used. Currently, most investigative efforts report 50% decarboxylation and 50% hydroxylation products in the production of biofuels by cytochrome P450 (Faponle et al. 2016). Therefore, there is a need to engineer the protein to enhance the yield of biofuel.

The use of quantum mechanics or molecular mechanics (QM/MM) techniques has facilitated the characterization of biochemical systems by researchers. Quantum-mechanical (QM) approaches fall into ab initio methods that entail the computation of a system devoid of prior knowledge (de Visser et al. 2014). Overall QM methods work by working out the Schrödinger equation. As much as QM methods can expedite the accurate computation of systems, their requirements increase based on the number of atoms or electrons in the system. Therefore, these methods are unbefitting for investigations involving enzymatic systems (de Visser 2009).

In contrast, alternative methods such as molecular mechanics (MM) that can deal with systems with a higher number of atoms are unable to expound bond-breaking or formation processes. As a result, innovative techniques that amalgamate the rapidity of MM with the exactitude of QM have been developed, leading to computational procedures that can give details of enzymatic systems. One such method is based on the density functional theory (DFT) alongside high computational power, which has led to increased interest in theoretical chemistry processes.

DFT computations are usually executed to provide backing for experimental biochemical and inorganic investigations, facilitate the interpretation of data and give directions for new experiments. This dissertation aimed at reviewing the use of the OleT enzyme in biofuel production to determine existing gaps, engineering the enzyme to improve its yield and reporting the findings of the investigation.

Project Objectives

The aim of the project is to engineer the OleT enzyme to improve its biofuel yield.

The objectives include:

  • Using computational methods to design and develop mutants that give a higher yield of biofuel products.
  • Discussing the outcomes of the computational endeavours and their application in practical settings.
  • Discussing the expected personal, academic, industrial, societal, and economic benefits of the project.

Literature Review

Production of Biofuels

To this point, four bacterial biosynthetic pathways that convert free fatty acids into hydrocarbons have been identified and exemplified. The first pathway, which was identified in cyanobacteria, is made up of two enzymes: an acyl-acyl carrier protein (acyl-ACP) reductase and an aldehyde decarbonylase (Choi & Lee 2013; Wang, Liu & Lu 2013). These two enzymes convert fatty acyl-ACPs to alkanes (Liu et al. 2014).

The second pathway comprises a cluster of three genes in Micrococcus luteus, which catalyse the head to tail linkage of two fatty acyl-coenzyme A (acyl-CoA) molecules to form alkenes with internal double bonds (Chen et al. 2015). The third reaction path uses a unique P450 decarboxylase OleTJE, which was isolated from Jeotgalicoccus sp. ATCC 8456 (Grant et al. 2015; Hsieh et al. 2017; Matthews et al. 2017a). OleTJE catalyses the direct removal of carbon dioxide molecules (decarboxylation) from long-chain fatty acids to produce α-olefins (Liu et al. 2014). This reaction occurs in the presence of H2O2.

The fourth reaction uses a type I polyketide synthase isolated from Synechococcus sp. PCC 7002 (Lin et al. 2015). This enzyme uses five consecutive steps to transform fatty acyl-ACPs into α-olefins (Lennen & Pfleger 2013; Mendez-Perez, Herman & Pfleger 2014). Out of the four methods, the simplest reaction is the one-step decarboxylation of fatty acids by OleTJE P450. In addition, this enzyme can use different molecules as substrates, for example, free fatty acids as substrates in the place of fatty acid thioesters. The ability of OleTJE P450 to use free fatty acids as substrates is considered an advantage in metabolic engineering because free fatty acids are abundant in nature (Faponle et al. 2016).

In addition, free fatty acids are easily manipulated in microorganisms such as E. coli, which is one of the most highly developed microbial cell industrial units (Chen et al. 2013; Marisch et al. 2013; Rosano & Ceccarelli 2014; Wang & Zhu 2018). Hence, OleTJE P450 fatty acid decarboxylative mechanism holds immense potential for engineering into a biological α-alkene-producing scheme.

Cytochrome P450

Cytochrome P450 is a superfamily of heme-containing external monooxygenases (Faponle et al. 2016). The name of this group of enzymes is attributed to their categorisation as hemoproteins (proteins that contain heme as a cofactor) and exhibition of distinct spectral features. P stands for pigment, whereas the term 450 signifies the rare absorption peak of the reduced CO-bound compound at a wavelength of 450 nm. Cytochrome P450s play a crucial role in the metabolism of drugs in humans, particularly in phase 1 reactions where drug molecules undergo various reactions to render them water-soluble (Faponle et al. 2016).

These enzymes are also found in other living organisms, for example, eukaryotic animals, plants, fungi and bacteria. The cytochrome P450 enzyme group is very diverse. Molecular investigations indicated that more than 20,000 distinct sequences of the enzyme had been identified as of January 2015 (Belcher et al. 2015). In humans, the site where cytochrome P450s are found predominantly is in the liver where they catalyse the breakdown of toxic substances (Zanger & Schwab 2013; Rendic & Guengerich 2014; Ji et al. 2015).

However, these enzymes also play a role in other physiological processes such as the biosynthesis of hormones, degradation of vitamin D, metabolism of alkaloids, steroids, eicosanoids, and fatty acids, the formation of bile acid, as well as the instigation of procarcinogens (Nebert, Wikvall & Miller 2013). These processes demonstrate the diverse roles of cytochrome P450s in the body and their versatility in the activation of substrates. In plants, cytochrome P450s play a central role in the metabolism of herbicides and secondary metabolites (Renault et al. 2014; Saika et al. 2014; Zhai et al. 2014).

The identification of more than 20,000 forms of cytochrome P450 as well as the diverse physiological functions of the enzyme implies that it has different isoenzymes (isozymes) in organisms that use it (Gricman, Vogel & Pleiss 2014). Isozymes are enzymes that exist in different molecular forms within the same organism. The substrate-binding pocket of P450 differs in shape and size in its various isozymes. The functional group of cytochrome P450 enzymes is made up of an iron (IV)–oxo heme cation radical moiety that is referred to as Compound I (Faponle et al. 2016). The characterization of Compound I can provide valuable data regarding the mechanism of cytochrome P450.

However, experimental exploration of compound I is convoluted by its transitory nature, which constrains the availability of information about its catalytic machinery. Theoretical and experimental investigations reveal that the most common reaction mechanism executed by the cytochrome P450 enzymes is aliphatic hydroxylation of substrates. Empirical studies have shown that aliphatic hydroxylation is a stepwise sequel that forms alcohol compounds as the main products.

Furthermore, this mechanism has been substantiated by computational modelling, which has revealed that the oxidation state consists of two spin states (lying close to each other) that yield two-state reactivity configurations, each having its own rate constant (de Visser et al. 2014). The two-state reactivity prototype also projects that the lifespan of the radical intermediary could lead to a split of the potential energy. Cytochromes P450 take part in catalytic reactions, for instance, dehalogenation, peroxidation, hydroxylation, sulfoxidation, S-dealkylation, desulfuration, and epoxidation, among many others (Zhao et al. 2014).

During these reactions, the enzyme needs to interact with its corresponding electron donors. This interaction is crucial because it confers enzyme specificity, which ensures adequate catalysis and acceptable reaction rates while guarding the system against side reactions.

The cytochrome P450 enzyme superfamily is among the most flexible naturally occurring biocatalysts. These enzymes use oxygen to catalyse their reactions. Electron-moving proteins transfer electrons to the enzyme to accelerate the activation of oxygen, whereas the electron carriers NADH or NADPH smooth the progress of substrate hydroxylation. Most cytochrome P450 variants are bound by membranes and have different hydrophobic substrates where they catalyse monooxygenation reactions. Some of the substrates involved in these reactions include biological compounds such as steroids, fatty acids, and prostaglandins.

On the other hand, extraneous compounds that take part in monooxygenation reactions include organic solvents, ethanol, anaesthetics and alkyl aryl hydrocarbon products. In monooxygenase reactions, one or more redox affiliate proteins are used to convey two electrons from NADH or NADPH to the heme iron in the reactive centre of the enzyme for dioxygen initiation. Subsequently, one atom of oxygen is inserted into their substrates during catalysis. This reaction mechanism has led to the naming of the enzymes as P450 monooxygenases. A typical cytochrome P450 catalytic cycle is shown below in Figure 1.

The cytochrome P450 catalytic cycle
Figure 1: The cytochrome P450 catalytic cycle (Belcher et al. 2014, p. 6526).

However, an exceptional cytochrome P450 family, the CYP152 family, includes members that rely on H2O2 solely as the only electron and oxygen donor. These members are known as P450 peroxygenases and they include P450SPα, P450BSβ, and OleTJE (Matthews et al. 2017a). The peroxygenase activity of P450 enzymes is usually considered an advantage in practical settings due to the cost differences between H2O2 and NADPH. H2O2 is significantly cheaper than NADPH and redox proteins.

The applications of the peroxygenase activity of cytochrome P450 enzymes include catalytic reactions involved in the biosynthesis of biological molecules as well as enzyme additives in laundry cleansing agents (de Visser 2009). Therefore, there is an increased interest in taking advantage of the oxidative role of cytochromes P450 for industrial purposes.

A significant reaction mechanism in the chemistry of cytochrome P450 enzymes is the desaturation of aliphatic compounds to produce olefins (Munro et al. 2018). Desaturation reactions bear numerous similarities to aliphatic hydroxylation because both processes start with the elimination of a hydrogen atom from the substrate (Ji et al. 2015). Nonetheless, the reaction subsequently divides into two possible product routes the first of which presents OH rebound to form alcohol products, whereas the other path leads to the elimination of a second hydrogen atom to produce an olefin and water (Matthews et al. 2017a).

The first reaction pathway is referred to as aliphatic hydroxylation while the second path is known as desaturation. Desaturation reactions are of immense interest in the production of biofuels from organic substances (Ghosh et al. 2016; Hendry et al. 2016; Carvalho et al. 2017).

Even though comprehensive scientific investigations into the chemistry of P450 have been conducted for decades, many gaps exist in the understanding of its reaction mechanism. For example, there is inadequate information regarding the intricate details of its hydroxylation and desaturation mechanisms as well as possible ways of bioengineering these reaction pathways to produce a specific product.

Cytochrome P450 Peroxygenase OleTJE

Peroxygenases are cytochromes P450 that make use of the peroxide shunt to catalyse the hydroxylation of long chain fatty acids to produce an alcohol (Beller, Lee & Katz 2015; Wang et al. 2017). Adding hydrogen peroxide or any other organic peroxide to the reaction medium facilitates the production of compound 0 or compound I without the need for expensive redox partners such as NADPH and NADH. The CYP152 family contains the most important cytochrome P450 peroxygenases (Hsieh & Makris 2016; Matthews et al. 2017b). A notable peroxygenase is cytochrome P450 peroxygenase OleTJE, which was isolated from Jeotgalicoccus sp. ATCC 8456 (Fu et al. 2015; Lee et al. 2015). OleTJE catalyses the unusual decarboxylation of long-chain fatty acids leading to the formation of α-alkenes in a reaction that uses H2O2 as the only electron and oxygen donor. An alternative reaction mechanism entails the use of an enzymatic redox cascade in addition to molecular oxygen and the recycling of NAPDH as shown in Figure 2.

Oxidative decarboxylation of fatty acids with OleT.
Figure 2: Oxidative decarboxylation of fatty acids with OleT. Path A uses H2O2 via the peroxide shunt, whereas path B uses O2 and NAD(P)H (Dennig et al. 2015, p. 8819).

Even though the decarboxylase activity of OleTJE proceeds with relative ease, the industrial manufacture of cheap α-alkene biofuels on a large scale while depending on the H2O2-dependent enzymatic system is not feasible. The limitations of this reaction are attributed to the need for large quantities of hydrogen peroxide, which would be costly. In addition, long-term exposure of biocatalysts to high concentrations of H2O2 is detrimental to their integrity and would lead to the inactivation of the enzymes. For these reasons, the H2O2-independent action of OleTJE is ideal for the commercial production of α-alkenes in microbial systems (Liu et al. 2014).

Aliphatic α-alkenes are valuable because of their potential use as biofuel substitutes for fossil fuels, manufacture of lubricants, detergents, and polymers (Pai et al. 2017). Therefore, studies on fatty acid decarboxylation by OleTJE are important and findings from such studies can inform the industrial manufacture of biogenic α-alkenes, which are renewable and environmentally friendly fuel sources.

H2O2-Independent Activity of OleTJE

H2O2-independent decarboxylation of fatty acids by OleTJE occurs in vivo in the presence of NADPH and O2 (Hsieh et al. 2017). Under these conditions, OleTJE enzyme efficiently decarboxylates long-chain fatty acids (fatty acids with carbon atoms ranging from 12 to 20) to produce terminal olefins. This process requires OleTJE to work alongside a fused P450 reductase domain RhFRED from Rhodococcus sp., CamA, CamB, or a separate flavodoxin/flavodoxin reductase from E. coli (Zhou et al. 2018).

The uncovering of the H2O2-free activity of OleTJE queries the monooxygenase-like machinery of the peroxygenase and directs future metabolic engineering efforts to enhance redox partners for efficient production of α-alkenes by OleTJE. However, the reaction liberates appreciable quantities of α- and β-hydroxylation by-products whose origin is not clear (Faponle et al. 2016). The peroxide independent activity of OleTJE shows potential in the production of biofuels because it does not pose a threatening environment for microbial systems as opposed to the peroxide-dependent pathway that destroys the integrity of microbial systems.

α- and β-Hydroxylation by-products of OleTJE

QM/MM inquiries into the bifurcation pathways have facilitated the identification and elucidation of the formation of three by-products of OleTJE. These studies have also made it possible for researchers to demonstrate how the enzyme can be engineered for peak desaturase activity. QM/MM investigations show that the availability of a substrate as well as the polarity of the solvent that is present in the binding cleft of the enzyme destabilises the OH-rebound path and permits a thermodynamically disadvantageous decarboxylation reaction through kinetic means as shown in Figure 3 (Faponle et al. 2016).

Reaction products obtained from fatty-acid activation by P450 OleTJE
Figure 3: Reaction products obtained from fatty-acid activation by P450 OleTJE (Faponle et al. 2016, p.5479).

Substrate Specificity of OleTJE Systems in Peroxide-Independent Pathways

The enzyme-substrate affinity of different OleTJE systems provides an indication of substrate preference of various OleTJE systems to improve the comprehension of the in vivo mechanism of fatty acid decarboxylases (Grant, Mitchell & Makris 2016; Xu et al. 2017). This information also points out the metabolic engineering of fatty acid biosynthesis in the right track. A number of straight-chain saturated fatty acids with even numbers of carbon atoms in the hydrocarbon chain have been used to resolve the substrate specificity of different OleTJE systems. Liu et al. (2014) used long-chain fatty acids ranging from C8 to C20 to establish the substrate specificity of OleTJE and OleTJE-RhFRED.

The proportion of the substrate that was converted into the expected α-alkene product was used as an indication of the enzyme system’s substrate preference. Liu et al. (2014) noted that OleTJE had the highest preference for myristic acid (C14) for olefin production with approximately 97% of the fatty acid being transformed to terminal olefins. The proportions of lauric acid (C12) and palmitic acid (C16) that were converted to terminal alkenes were lower than the expected levels. OleTJE could only transform small quantities of stearic acid (C18) and arachidic acid (C20) to alkenes (Figure 4). However, the enzyme could not catalyse the decarboxylation of capric acid (C10) or caprylic acid (C8).

The substrate preference spectrum for OleTJE (A) and OleTJE-RhFRED (B)
Figure 4: The substrate preference spectrum for OleTJE (A) and OleTJE-RhFRED (B) (Liu et al. 2014, p. 5).

On the other hand, it was demonstrated that the OleTJE-RhFRED system had a higher affinity for fatty acids with relatively short chain lengths. This finding implied that the presence of a P450-reductase induced a small structural modification in the active site of OleTJE. The highest activity for OleTJE-RhFRED was shown for lauric acid at a conversion rate of 83.8%. Nevertheless, no activity was observed for arachidic, capric, and caprylic acids (Li et al. 2014).

Another outstanding observation was that the system consisting of OleTJE-RhFRED, NADPH, and O2 demonstrated a lower activity than OleTJE and H2O2 for all fatty acids tested except lauric acid. The observed phenomenon indicates that OleTJE could have evolved over time to augment its peroxygenase activity.

On the other hand, Dennig et al. (2015) established the enzymatic oxidative decarboxylation of short-chain fatty acids ranging from C4 to C9 using OleT, O2 as the oxidant, and NAD(P)H as the electron donor to produce the corresponding terminal C3 to C8 alkenes. Product titres of up to 0.93 gL-1 and total turnover numbers exceeding 2000 were given off. The development of a well-organized electron-transfer series using putidaredoxin CamAB together with the reutilisation of NAD(P)H without regard for formate, glucose, or phosphite made it possible for these yields to be attained (Dennig et al. 2015).

Even though CamAB could oxidize NAD(P)H without an appropriate electron acceptor, no H2O2 was detected after the full oxidation of 1 mm NADH, which indicated the likelihood of a direct electron displacement from CamAB to OleT. This observation also ruled out the possibility of inactivation by hydrogen peroxide. Consequently, the absence of catalase did not affect the throughput of the system.

Astoundingly, fatty acids with short chains (fatty acids with carbon atoms ranging from C4 to C11) were decarboxylated, which resulted in the production of 0.04 to 2.45 mm of terminal alkenes. The existing temperatures and the number of carbon atoms in the fatty acid chains had an effect on the overall yield of the terminal alkenes. The highest activity of the enzyme was recorded for a C18 fatty acid at room temperature, whereas C12 yielded the best outcome at 48oC (3.26 mm of 1-undecene).

On the whole, high conversion rates were recorded for fatty acids ranging from C10 to C16 at temperatures of 48oC. The CamAB-OleTJE system facilitated the production of propene from butyric acid in only one step. This observation disproved the previous hypothesis that OleT had a preference for long-chain fatty acids compared to short-chain fatty acids (Dennig et al. 2015). It was presumed that long-chain fatty acids interacted with the substrate pocket, whereas smaller substrates were held in the binding compartment as evident in correlated enzymes, for example, CYPBSß.

The low conversion rates of fatty acids with intermediate chain lengths (C10 and C11) was attributed to poor binding of the substrate on the active site of the enzyme (Dennig et al. 2015; Hammerer, Winkler & Kroutil 2018). Decarboxylation was the principal reaction, with approximately 86 to 99% selectivity in fatty acids C18 to C22.

The formation of α- and β-hydroxylated products was influenced by chain length and reaction temperature (Faponle et al. 2016). There was 62% α-/β-hydroxylation when fatty acids with 9 carbon atoms reacted at temperatures of 48oC. The observed harmony between CamAB and OleT was associated with similarities in the synchronisation of the adjacent cysteine thiolate moieties in the ferrous states of OleT and P450Cam, which allowed the attachment of oxygen to the heme iron.

Most P450 enzymes of bacterial origin depend on class I electron-transfer schemes (Dennig et al. 2015; Liang et al. 2018), which explains the enhanced catalytic performance of OleT together with CamAB, as opposed to systems illustrated in the past such as OleT–RhFRed and Fdr/Fdx. These findings indicate that metabolic engineering efforts directed at the peroxide independent activity of OleTJE should consider pairing the enzyme with CamAB to take advantage of different lengths of fatty acid substrates.

Product Distribution in OleTJE Catalysis

Environmental effects may influence the branching and regioselectivity inclinations of enzymes. Therefore, Faponle et al. (2016) used QM/MM to study the full OleTJE enzyme and investigate its bifurcation pathway and product distribution. The main objectives of the study were to identify the source of P450 OleTJE product formation and predict the features that break up the decarboxylation from hydroxylation pathways.

The findings revealed that were no significant changes in geometry, spin state collation, and comparative energies. Geometry studies showed that there were no substantial disparities between hydrogen-atom extraction impediments. In all instances, the hydrogen atom was nearly equidistant from the donor and acceptor atoms, notwithstanding the type of the QM region studied. It was evident that hydrogen atom removal impediments by Cpd I from the Cα and Cβ loci of the substrate were close in energy on the doublet or quartet spin states. All four impediments fell within a range of 3.2 (BS1) or 4.0 (BS2) kilocalories per mole for Sn400 (Faponle et al. 2016).

The observed hydrogen atom extraction impediments corroborate empirical findings of a combination of reaction outcomes from Cα and Cβ hydrogen atom removal. The minor disparity in the entropies of the two reaction courses suggests that the chemical reactions that lead to the elimination of a hydrogen atom between Cα and Cβ are competitive. However, the ratio of the Cα and Cβ pathways as depicted by the comparative energies of the transitional state moieties 4,2TSHA,α and 4,2TSHA,b were strongly influenced by environmental factors, for example, the associations between hydrogen bonds (Faponle et al. 2016).

The experiment by Faponle et al. (2016) showed that hydrogen atom extraction was the rate-limiting stage in the reaction mechanism. Therefore, the authors approximated a large kinetic isotope effect and computed an estimate for the substitution of hydrogen atoms by deuterium (molecular mass of 11). Geometry shots were also taken to determine possible hydrogen atom removal from the Cγ locus of the substrate (Sn400). It was noted that the reaction energy exceeded the reactant complex by more than 30 kcal mol-1 hence this mechanism was overlooked (Faponle et al. 2016).

The system slowed down following the removal of hydrogen atoms to form a radical intermediate (IH,α or IH,β) with an iron(IV)–hydroxo group and a reactive substrate on either the Cα or Cβ locus of the substrate. The transitional species 4,2IH,α only reacted through radical rebound past the rebound transition state 4,2TSreb,α thereby leading to the formation of α-hydroxo arachidonic acid upshots (4,2 P OH,α). The rebound barriers recorded during the investigation were larger and higher in energy than is the norm for P450 hydroxylation reactions. Under normal circumstances, the low-spin

pathway does not have any impediments. A decarboxylation reaction cannot take place following the removal of Cα hydrogen; it requires the additional removal of a hydrogen atom from Cβ to Cα, which requires high energy.

The competitive pathways that resulted in decarboxylation and hydroxylation products (PD,β and POH,β, correspondingly) via transition states TSD,β and TSreb,β came from the radical intermediate 4,2IH,β, as estimated by snapshots Sn300 and Sn400. The high-spin pathway through 4 IH,β gave rebound and decarboxylation hindrances that were almost similar in energy, hence could be regarded as competitive.

On the other hand, the low-spin path produced quality decarboxylation products from 2TSH,β directly without stable radical intermediates. These findings indicated that the rebound barriers were bigger than conventionally found in ideal complexes, which was attributed to numerous hydrogen-bonding connections to the iron–hydroxo group. These restrictions increased the obstacles to the rotation and rebound pathways. In addition, the hydroxylation process was dominant following Cα-H removal while decarboxylation was favoured following the extraction of the Cβ-H atom on the doublet spin state.

On the other hand, the surface of the quartet spin state was expected to liberate a mixture of products. This observation implies that the reaction should be modified to promote Cβ-H extraction and impede Cα-H removal to promote decarboxylation. The elevated rebound barrier was linked to hydrogen-bonding interactions because of the water channel in P450 OleTJE in addition to the polarity of the nearby carboxylate-Arg245 salt bridge (Faponle et al. 2016). The importance of water channels in P450 OleTJE was also demonstrated.

Bifurcation Pathways of OleTJE

Faponle et al. (2016) used a valence-bond (VB) diagram to elucidate the bifurcation pathways of P450 OleTJE (Figure 5). The illustration in Figure 5 shows that the wave function of the reactant (YI) merges with an energised state in the geometry of the yields, which is Yreb* for the alcohol products and YD* for the decarboxylation output. The wave function of the product joins the energised excited state in the reactant geometry (4,2IH,β) in the same way.

The extent of the reaction inhibition is a function of the excitation energy needed to transform reactant to product in the reactant geometry. Thus, the bifurcation machinery is determined by the comparative magnitude of the excitation energies Greb or Gdecarb (Faponle et al. 2016). Therefore, the microelectronic discrepancies between the wave functions in the ground and excited states regulate the barrier inhibitions and rate constants for the rebound and decarboxylation reactions.

Evaluating the VB structures in the geometry of 4,2IH,β with emphasis on ground and excited states clarifies the parameters that govern the bifurcation variables. For OH rebound and the development of hydroxylation outcomes, the πyz/π*yz set of orbitals was divided into two atomic orbitals, 3d yz,Fe and 2py,O. The division of orbitals process expended an estimated Eπ/π*yz in energy from the system, which was estimated at 37.2 kcal mol-1. The 2py,O orbital coupled with ΦSub, a substrate radical, to create the C-O bond orbital (σC-O) in a process that consumed a bond severance energy BDECO of 87.2 kcal mol-1.

On the other hand, since the 2py,O had a double occupancy in 4,2IH,β, it was imperative to upgrade one electron to a free low-lying orbital. The two available orbitals were 3dxz in the doublet spin state and the σ*z2 orbital in the foursome spin state, Eexc. The excitation energy was approximated at 75.3 kcal mol-1. Thus the value of Greb was predicted as 25.3 kcal mol-1. As a rule, the height of the barricade is usually 33.33% of the promotion gap, which implied that a rebound barrier of approximately 8 kcal mol-1 was applicable.

A valence-bond illustration of the electronic attributes of decarboxylation versus OH-rebound reactions.
Figure 5: A valence-bond illustration of the electronic attributes of decarboxylation versus OH-rebound reactions. The lines depict chemical bonds, whereas the dots show valence electrons (Faponle et al. 2016, p. 5482).
Computed bifurcation barriers of different metal intermediates
Figure 6: Computed bifurcation barriers of different metal intermediates (Faponle et al. 2016, p. 5481).

However, a different scenario was witnessed in the decarboxylation reaction where two processes occurred simultaneously. The C-Cα bond of the substrate disintegrated (BDECC) while the leaving CO2 group was oxidised to CO2. The second process denoted the electron affinity of CO2 (EACO2). As a result, the electron duo in the C-Cα orbital (σCCα) was split, after which the electron on Cα coupled up with the reactive group on Cβ to create a π bond (π CαCβ) possessing energy Eπ.

The additional electron was taken up by the iron–hydroxo compound, which was reduced by an amount of energy that corresponded to its electron affinity (EAFeOH). The electron affinity of the iron complex was estimated at 56.5 kcal mol-1, whereas the strength of the C-C bond in the substrate, BDECC, was found to be 7.8 kcal mol-1 (Faponle et al. 2016). Conversely, the singlet to triplet energy gap was used to compute the energy of the π bond, which yielded a value of 73.7 kcal mol-1.

A promotion gap for decarboxylation (Gdecarb) of 3.4 kcal mol-1 was computed using the standard electron affinity of CO2 (‑0.60 eV). The findings pointed towards a minor reaction barrier that was estimated to be less than 1 kcal mol-1. Figure 6 indicates the calculated bifurcation impediments of different radical intermediates.

Factors Affecting Decarboxylation Reactions of OleTJE

It is necessary to ascertain the factors that influence a reaction path that leads to the formation of the desired product before designing an enzyme to overproduce a given product. Investigations indicate that the decarboxylation pathway of OleTJE is mainly influenced by substrate factors (Xu et al. 2017). These aspects include the intensity of the C-Cα bond, the energy required to create the π bond between Cα-Cβ, and the electron affinity of CO2 (Faponle et al. 2016; Fang et al. 2017).

The effect of the oxidant on the decarboxylation reaction is only felt through the electron affinity of the iron (IV)–hydroxo moiety. Consequently, boosting the regioselectivity of decarboxylation as over the hydroxylation activity should focus on weakening the hydroxylation path by engineering the substrate-binding compartment or using alternate substrates. The computations presented by Faponle et al. (2016) show that the substrate-binding compartment of P450 OleTJE is highly polar and has numerous hydrogen-bonding interactions, which are crucial to the disruption of the radical rebound reactions that lead to the formation of α- and β-hydroxo fatty acids.

In addition, electron-extracting moieties attached to the Cα or Cβ locus would elevate the energy required to form the C-O bond, which would increase the radical rebound barriers. In contrast, the decarboxylation reactions would be boosted further by using substrates with an alkene bond between Cγ and Cδ. Such substrates would promote the π conjugation with the Cα-Cβ π bond, thereby leading to increased Eπ energy and stabilisation of the decarboxylation process.

Theoretical Background of QM Simulations

Quantum mechanics proposes that the wave function can provide all attributes of a molecular system, which in turn solves the Schrodinger equation (McMahon 2013; Mohebbi, Abbaszadeh & Dehghan 2013; Edmonds 2016). The overall energy of a molecule is depicted as an eigen function of the Hamilton operator (H) as shown in equation 1 below.

Theoretical Background of QM Simulations

The Hamilton operator encompasses the kinetic and potential energy influences of all the electrons and nuclei within a molecule (Omnes 2018). This relationship is summarized by equations 2 and 3 where h denotes Planck’s constant, me is the mass of the electron, e is the charge of an element, ε0 symbolises the permittivity in vacuum, and rij depict the expanse between. On the other hand, the spaces between two electrons and two nuclei are portrayed by rij and RIJ, correspondingly.

Consequently, the Hamiltonian in the second equation is divided into four constituents that exemplify critical components the first of which is the electrons’ kinetic energy. The other components include the attraction force between the electron and nucleus, the repulsive forces between electrons and the repulsion force between two nuclei.

Equation

Therefore, elucidations of the Schrodinger equation are influenced by the total number of electrons and nuclei in a chemical system (Jensen 2017). These solutions increase in complexity as the size of the system increases (Harrison & Valavanis 2016; Le Roy 2017). The complexity of the Schrodinger equation implies that the equation can only be solved precisely for the hydrogen atom and ion. Estimations may be needed for larger chemical systems, which has led to the development of numerous computational chemistry approaches that encompass varying features and levels of precision (Alolyan & Simos 2015; Ning & Simos 2015; Grimme & Schreiner 2018).

On the whole, meticulous benchmarking and calibration are needed to produce accurate results, particularly for systems involving transition metals (Genheden et al. 2017; Mardirossian & Head-Gordon 2017). In such systems, the stability and reproducibility of techniques are ascertained by comparing findings to experiential data or advanced techniques.

Hartree–Fock (HF) theory provides the simplest solution to the Schrodinger equation (Koch & Holthausen 2015; Yanai et al. 2015). The assumption made is that every electron moves in the mean field of other electrons. Therefore, one-electron wave functions are expanded into a group of known orbitals referred to as the basis set. Nonetheless, solving the Schrodinger equation does not imply that the computation is over but only provides the energy of the starting structure (de Visser et al. 2014). The geometry undergoes modifications followed by new computations. Subsequently, the energy is minimised with regard to all geometric coordinates until the first derivatives of the energy reach a specific value. At this point, an optimised geometry is attained.

The HF technique provides reasonable geometries. However, its energies are insufficient to replicate chemical reaction machinery and experimental barrier heights. Coupled-cluster calculations using singles, doubles, and perturbatively treated triples (CCSD(T)) techniques are regarded the gold-standard in QM (Sparta & Neese 2014; Calbo et al. 2015; Liakos & Neese 2015). The main disadvantage of these methods is that their computational demands make them relevant for simple systems with 50 atoms or less. Biomimetic complexes with 100 or more atoms require simpler computational methods (de Visser et al. 2014).

The density functional theory (DFT) was advanced in the 1960s as a substitute for the wave function approaches such as the B3LYP method (Obot, Macdonald & Gasem 2015). The electronic energy is computed as an unknown function of the electron density. The creation of functionals that rely on electron density gradient led to major leaps in the DFT methodology and the development of semi-empirical methods. However, recent studies have highlighted the weaknesses of the B3LYP method, particularly with systems that include transition metals, thereby leading to the use of alternative DFT methods (Maxwell, Pendás & Popelier 2016; Medvedev et al. 2017).

However, there is no general and reliable method that has been recommended. Therefore, researchers use different methods. Furthermore, HF and conventional DFT methods cannot explain dispersion interactions. A recent study has shown that correcting a simple empirical dispersion leads to a significant improvement in the energies computed using the B3LYP method. Therefore, it is necessary to validate DFT computations using experimental data.

QM computations provide the electronic energy of the system under investigation. However, these energies can be likened to experimental values by taking into account the effects of solvents, vibrations and entropy. The extent of occupation of each individual vibrational energy level is determined by a Boltzmann distribution at a specific temperature. The specific partition functions for the translational, gyratory and vibrational entropy constituents can be used to compute the overall entropy of the energy levels. Usually, the partial derivative of the free energy to the temperature at an invariable pressure is used to compute the entropy of a molecule. Quantum chemical software packages have these entropy derivations to estimate the free activation energy of a molecule, which can then be equated with experimentally derived values.

Computational chemistry approaches liberate energies that denote gas-phase data. Therefore, computational outcomes can only be compared with experiential findings following the incorporation of solvents into the computations. Therefore, it is necessary to add specific solvent molecules to the chemical model and re-optimise the structures. Nonetheless, introducing many atoms into the system complicates it and calls for complex computations that may not be feasible because solvent molecules exist in different conformations. Therefore, the most commonly used method entails the addition of solvent adjustments to the energetics with the aid of an implicit solvent model (Jensen 2017).

The QM/MM approach also makes it possible to incorporate the effects of the surrounding into QM computations. The guiding principle of QM/MM techniques is to split the chemical system into two or more distinct shells, which are referred to as the outer and inner cores as shown in Figure 7.

The splitting of an enzymatic system into an inner core (QM region) and an outer core (MM region)
Figure 7: The splitting of an enzymatic system into an inner core (QM region) and an outer core (MM region) (de Visser et al. 2014, p. 267).

The computations of the inner core are implemented using precise QM or DFT techniques, whereas the calculations of the outer core that holds the rest of the protein and solvent are done using computationally inexpensive methods such as molecular mechanics. The inner core is of utmost importance because it holds the essential features of the substrate and the enzyme active site. The interactions of the inner core with other proteins hold it in position. These interactions maintain the conformation of the enzyme active site. Consequently, it is possible to blend precision of QM methods with the rapidity of MM methods.

The literature review has provided useful insights that guide the engineering of OleTJE to enhance its production of biofuels. The following insights have revealed that the enhancement of olefin product distributions requires site-directed mutations that impede α- and β-hydroxylation as well as altering the substrate-binding compartment to limit hydroxylation rebound reactions (Li et al. 2015). It is also evident that the choice of substrate, reaction temperature and OleTJE systems or combination of systems (Dennig et al. 2015) is critical in obtaining optimal terminal olefin yields from fatty acids. These pointers informed the computational methods in the methodology section of the project.

Methodology

Computational modelling of biosystems lays a basis for investigators to set up experimental inquiries because it excludes redundant variables that would otherwise be a waste of resources. The fact that quantum mechanics/molecular mechanics (QM/MM) techniques use computed reaction machinery and energy settings to provide details regarding the starting point of biological and chemical reactions, as well as their by-products, has increased the popularity of these methods over the years (de Visser et al. 2014).

As a result, researchers can collect data regarding factors that have an effect on enzyme reactions, rates and specificities by examining the different states of the molecules involved in enzymatic reactions. On the other hand, cluster replicas, for example, the Density Functional Theory (DFT) facilitate voluminous computations in a short span by probing precise features of a single protein (de Visser 2009; Jones 2015; Melander et al. 2016).

Despite cluster models providing information about the intrinsic activity of enzymes’ active sites, there is inadequate information concerning intricate details such as key bonding patterns and lengths. This section of the project provides a detailed procedure on the use of general QM/MM optimization protocols to promote the decarboxylation mechanisms of cytochrome P450 OleTJE in the production of terminal olefins. The proposed elements of improvement as highlighted in the literature review were incorporated into the general QM procedure (Figure 8) to determine the most effective manipulations that would yield favourable outcomes.

QM/MM Methodology

The QM/MM method was initiated by finding heavy-atom coordinates of the protein (OleTJE enzyme) from the protein databank library. These coordinates were available in a pdb format. As a result, the entire structure required editing as well as further adjustments before any simulations could be run. Examples of relevant modifications that were done to the initial pdb structure entailed the addition of hydrogen atoms, missing chains, substrates and amino acid residues. Another crucial modification was the removal of crystal additive buffers to leave the protein structure in its native state.

The active site of the structure was also considered during the modifications. The resultant molecule was then solvated. These procedures made it possible to stabilize the system to encourage the correct folding of proteins. Snapshots that were required to execute the QM/MM simulation were then produced through a molecular dynamics (MD) simulation. This process led to the split up of the biochemical system into inner and outer cores.

The computations of the inner core were completed using ab initio or DFT techniques that are reported to have high levels of accuracy (Lewars 2016). On the other hand, outer core calculations were achieved using MM approaches. These procedures have been validated using comparisons of experimental and theoretical data in various studies such as those conducted by Zeida et al. (2014), Wang et al. (2014) and Quesne et al. (2016). A sequential illustration of the specific steps of the QM/MM methodology is indicated in Figure 8.

A flowchart for the general set-up of QM/MM calculations from the starting structure in pdb format to the production of QM/MM snapshots
Figure 8: A flowchart for the general set-up of QM/MM calculations from the starting structure in pdb format to the production of QM/MM snapshots (Quesne et al. 2016, p. 2565).

Selection of Starting Structure in pdb Format

The accuracy and precision of QM/MM studies are dependent on the suitability of the starting materials. Therefore, researchers ought to be cautious when choosing the initial crystal structure from the protein databank. The ideal starting structure was not readily available from the protein databank. Most pdb files in the pdb database contained information indicating whether the structure was incomplete and whether any residues were missing.

As a result, modifications were unavoidable. However, the most comprehensive structure that contained many of the atoms in their correct positions was chosen to cut down the time required for modifications. Thus, the student ascertained the nature of the missing attributes (determined whether they were protein side-chains, loop areas, substrates, binding proteins, co-factors, or reducing partners) and selected the file with the closest semblance to the preferred structure.

This procedure also minimised the likelihood of errors during the modification of the structure. Another useful factor that was considered during the selection of the starting crystal structures was the resolution of the pdb files. Files with high resolution were preferred to those of low resolution because it was easier to identify distinct atoms as well as their location in the crystal structure. Therefore, pdb files with a minimum resolution of 2 Å were selected for the subsequent analyses. In addition, the investigator chose files with Rfree values of 0.3 or lower to ensure that only reliable structures were selected.

Modification of the Starting Structure

QM/MM computations of fast reactions provide the most informative data on enzymatic reactions. However, the pace of these reactions and their corresponding short-lived intermediates complicate the process of finding their crystal coordinates. This problem is usually circumvented by modifying a structure that is closely related to the target reaction. Therefore, a previous structure in a resting state was taken and modified appropriately.

Pertinent modifications of the starting structure were informed by the findings and recommendations of Quesne et al. (2016). These changes included involved substrate modification and alteration of the enzyme’s active site. Faponle et al. (2016) recommended that in order to augment the regioselectivity of decarboxylation reactions of OleTJE as opposed to hydroxylation reactions, it was necessary to lay emphasis on destabilising the hydroxylation path by engineering the substrate-binding compartment.

Another viable option was to use alternate substrates, which would bind the active site without leading to product formation thus inhibiting the formation of hydroxylation products. The computations presented by Faponle et al. (2016) presented adequate calculations demonstrating the high polarity of the substrate-binding compartment of P450 OleTJE, which was capable of forming numerous hydrogen bonds.

These bonds were instrumental in the interruption of the radical rebound reactions leading to the formation of α- and β-hydroxo fatty acids. Therefore, it was conjectured that an additional increase in the polarity of P450 OleTJE’s active site would hold back the reactions that form hydroxylation products. Another amendment was appending electron-extracting moieties to the Cα or Cβ locus to elevate the energy needed to form the C-O bond and increase the radical rebound barriers.

Substrate Docking Protocols

In conventional enzymatic reactions, the binding of a substrate to an enzyme prompts a conformational modification in the three-dimensional organization of the resultant enzyme-substrate complex. However, this reaction proceeds at a pace that is too fast to allow the crystallisation of the analogous structure. Advances in technology, particularly molecular and computational docking procedures have made it possible to predict the orientations and loci of binding substrates. Docking procedures exploit basic data about the substrate and protein surface, for example, the types and configurations of atoms, atomic modifications and polarity.

An online server, SwissDock, which is available at www.swissdock.ch, was used for this step. These data were likened to find the most appropriate lock-and-key pairs. The enzyme and substrates were tested in each probable orientation after which the energies for each alignment were worked out via MM methods. The possible enzyme-substrate complex was pinpointed by picking out the arrangement with the least energy. The Autodock software package was used in this process (Quesne et al. 2016).

The Addition of Hydrogen Atoms and Protein Groups to the Structure

The execution of geometrical optimisations was followed by the addition of hydrogen atoms to the protein structure in the pdb file. A number of online programs and packages were beneficial in enhancing the detectability of the correct positions of hydrogen atoms. The addition of hydrogen atoms was determined by the oxidation states, pKa systems and hybridisation of the atoms.

Selection of Molecular Modelling Replicas

The energy bonds were minimized using the right field force prior to the conducting of QM/MM calculations. Energy minimisation made it possible to achieve correct bond lengths and hybridisation structures. Several field forces were available. However, since this project involved an enzyme, e-biochemical field force such as GROMOS, CHARMM, or AMBER were ideal for the process. The force field calculations were directed by the following equations.

Equations

The first equation shows that MM energy (EMM) is the aggregate of bonded and non-bonded interactions. In the above equations, d expresses the bond length, ϴ designates the angles of valence bonding and x symbolises torsional dihedral angles. The dihedral attribute is a cosine function that considers the torsional multiplicity (n) and phase (δ) with respect to a bonding constant (kx). The non-bonding energy between two atoms A and B is distributed between two interactions: a Lennard-Jones expression and a Coulomb interaction as shown in equation 6.

The Lennard-Jones expression consists of the rAB term that identifies interatomic expanse, whereas ϵAB and σAB portray the depth, breadth and site of the inactive well. On the other hand, the (qAqB) term symbolises the operational atomic charge of the two atoms A and B, whereas ϵ0 designates permittivity in a vacuum. The value of qAqB increases as the distance between the two atoms and ϵ0 decreases.

Solvation of the Protein Structure

One inadequacy noted in the crystallised protein structures in pdb format was the lack of some of the water molecules that should normally be available. This insufficiency was attributed to the fact that the removal of water from protein suspensions was a pre-requisite to the crystallisation of protein structures (Jindal & Warshel 2016). Therefore, it was mandatory to add water to the initial starting structure to re-solvate it and make a replica of a realistic enzyme model at room temperature (Liu et al. 2015; Piana et al. 2015). Re-solvation also made better the in vivo and in vitro representation of polarity (Bellissent-Funel et al. 2015).

The major limitation of available resolvating methods was that it was difficult to add water molecules to the innermost regions of the protein. In addition, water molecules have the tendency to leave behind pits in the boundary of the protein, which gets in the way with the calculations (Jindal & Warshel 2016).

Therefore, it was necessary to conduct several re-solvation attempts of the system until the protein structure could not hold any more water molecules. Molecular, dynamic equilibration methods were used to transfer water molecules into the pits. One method required conducting a long MD simulation while exerting pressure from the perimeters of the system to compel the conveyance of water into the centre of the protein. The solvation process was followed by a confirmation of the positions of the additional water molecules to ascertain that no hydrogen bonds had been formed using atoms that were expected to play a role in the subsequent reaction pathway.

Heating, Equilibrating and Conducting Molecular Dynamics Simulation on the Solvated Structure

MD simulations can only be conducted on a system whose temperatures have been raised to the expected reaction temperatures to equilibrate the complex. Consequently, the temperature of the system was elevated in stages at small intervals ranging from 0oK to 298oK. During this process, the applied force field was used to minimise the hydrogen atoms and solvent molecules. When the desired temperature was achieved, the protein backbone became flexible, which equilibrated the system. At that point, MD simulations were run for 500 to 1000 picoseconds to generate snapshots of the protein structure, reaction intermediates and products for the required QM/MM computations. Different temperatures were investigated to determine the impact of temperature on the decarboxylation mechanism of OleTJE. An average of 6 simulations were done for each parameter.

Selection of Snapshots from the MD Simulation

It was necessary to choose several snapshots from the MD simulation because each snapshot was responsible for different computations in connection with the reaction rate, barrier height and product formation. Furthermore, alterations in the patterns of hydrogen bonds within a single water molecule were likely to change the overall energy by approximately 4 kcal mol-1 (Hayashi et al. 2017).

Therefore, taking a number of snapshots made it possible to examine different conformations of the protein. Numerous MD snapshots were taken from various points down the MD path graph. Each recreation was taken as the starting point for the creation of reaction profiles. In addition, snapshots from diverse protein milieus had varying parameters, for example, local minima of comparative energies, spin-state gaps and ordering. Therefore, the mean of the findings was obtained when assessing reaction schemes from more than one snapshot.

Selection of the QM Region, QM Methods and Basis Set

The choice and size of the QM region have an effect on the energetics of a reaction path and the electronic configuration of the active enzyme (Kulik et al. 2016; Karelina & Kulik 2017). Selecting a small QM region can lead to wrong descriptions of the polarisation of the system (Jindal & Warshel 2016). On the contrary, extremely large QM regions can result in excessive polarisation. Incongruences in bond lengths when relating QM regions of varying sizes are responsible for these observations. Therefore, it was necessary to examine several QM regions of assorted sizes to find out consistencies and select the correct regions.

Electrostatic and Mechanical Embedding Methods

Calculations of energies in the QM region should take into account the charges in the MM zone. Consequently, electrostatic or mechanical-embedding systems were helpful in this endeavour. The mechanical-embedding procedures determined van der Waals forces and electrostatic influences of the MM atoms in the QM zone (Jindal & Warshel 2016). The QM and MM energies were separated using geometric parameters such as the expanse of van der Waals interactions. Thereafter, the QM/MM expressions were assimilated directly into the force field estimations.

Calculation of QM/MM Energies

QM/MM energies were calculated by completing geometrical optimisations on the QM and MM zones. Two methods of optimisation were possible: adiabatic and diabatic (Barbatti et al. 2014; Duarte et al. 2015; Cabeza de Vaca, Acebes & Guallar 2016). The adiabatic geometry optimisation approach entailed relaxing the macro element exchanges inside the QM core as well as the micro-component interactions brought about by the neighbouring MM environment.

In contrast, the diabatic style included the freezing of the QM core region in the course of optimisation of the MM environment and vice versa. Nevertheless, linker atoms interacted in the boundary of the two regions, which meant that the total QM/MM energy was not contributed by the QM and MM energies only. Hence, it was necessary to consider the impact of these interactions during the computations. These coupling terms are described by the following equation.

EQM/MM= EQM +EMM+ E border …….…………………………………………7

QM/MM computations can be achieved by assimilating an MM into a QM software package. In this project, the ORCA software was used for the computations. ORCA is a contemporary ab initio and DFT package that can execute geometry optimisations and to forecast many spectroscopic strictures at various stages of theory. ORCA integrates advanced ab initio quantum chemical techniques that apply the configuration interaction and coupled cluster techniques into its functionalities (Metz et al. 2014; Dohm, Spohr & Korth 2017). These calculations can be executed using additive and subtractive procedures as shown in equations 8 and 9.

E Whole QM/MM = E Whole MM + E QM (I+L) – EMM (I+L) ………………………..8

E Whole QM/MM = E Environment MM + E QM (I+L) + E Coupling QM/MM ………………9

The subtractive approach, which was used in the project, involved taking advantage of linker atoms (L) to conceal the inner subsystems. The EQM (I + L) term represented the sum of the QM region (I) and linker atoms (L) at the QM level. Conversely, the EMM (I + L) term denoted a similar calculation at the MM level. However, it was important to subtract the MM energy of the inner subsystem from the energy of the entire system (E whole MM) to preclude double tallying of energies. The main benefit of the subtractive method was that it did not need a QM/MM joining term (Duarte et al. 2015; Řezáč 2016; Roßbach & Ochsenfeld 2017; Cao & Ryde 2018). The findings of the QM/MM processes are presented and explained in the next chapter.

Soundness, Rigour and Limitations of the Methods

To ensure the reliability of the entire procedure and the outcomes, the starting structure was obtained from the protein databank library available at www.rcsb.org. This resource is recognized among the best sources for published crystal structure data where published enzymatic crystal structure coordinates are kept and systematised.

The accuracy of the subsequent procedures was ensured by scrutinizing the initial starting structure to ensure that no errors were present. Additional modifications were also done to this structure as described in the methods section. The rationale or these modifications was to ensure that all essential components were in place before cutting down restrictions on the protein’s backbone to forestall improper folding (Quesne et al. 2016). Another precaution that was taken to ensure the rigour and soundness of the method was taking into account the energy changes of the QM and MM regions separately before considering the energy of the inclusive QM/MM. This procedure made it possible to identify any changes in the overall system that were not reflected in the individual QM and MM energies.

The procedure used in the entire experiment followed the schematic representation indicated in Figure 7. This procedure has been validated in various QM/MM experiments conducted by Zeida et al. (2014) and Wang et al. (2014). Zeida et al. (2014) used experimental means to show demonstrate the oxidative reaction mechanism of the enzyme peroxiredoxin from Mycobacterium tuberculosis AlphE. Through QM/MM techniques, it was established that the enzyme reduced hydrogen peroxide at a rate constant of 8.2×10-4 m-1 s-1, which was consistent to a rate-determining activation barrier of 5.7 kcal mol-1(Zeida et al. 2014).

Nonetheless, it was observed that this reaction was too fast to permit the implementation of accurate kinetic studies using flow methodology. Therefore, QM/MM methods obtained a theoretical activation energy of 4.8 kcal mol1. This difference was attributed to the impact of the environment on the energy profile. However, the experimental value had a close resemblance to the theoretical value, which validated the protocol.

On the other hand, Wang et al. (2014) used a similar approach to look into the enzyme machinery in the repair of exocyclic etheno bridge damage in adenine DNA bases. The authors used QM/MM computations to investigate the repair machinery involved in the mending of the etheno-bridged adenine. This repair process was conducted by AlkB enzymes, particularly iron(IV)-oxo species. Previous experiential studies that employed mass spectrometry and in crystallo separation reported that the etheno-bridge adenine was rectified by the generation of an epoxide that underwent subsequent modifications to form a glycol.

Theoretical computations replicated the energy barrier for the rate-limiting step, which was also observed in experimental trials. This step was also influenced by the pH of the reaction. A protonated enzyme/substrate complex from pdb structures was developed, solvated and exposed to an MD run prior to the selection of three snapshots for QM/MM computations. Wang et al. (2014) a theoretical energy barrier for the rate-determining step as 18.6 kcal mol1, whereas the experimental rate-limiting barrier was shown to be 21 kcal mol1. These findings also indicate the dependability of the QM/MM technique.

Reference List

Alolyan, I & Simos, TE 2015, ‘A high algebraic order multistage explicit four-step method with vanished phase-lag and its first, second, third, fourth and fifth derivatives for the numerical solution of the Schrödinger equation’, Journal of Mathematical Chemistry, vol. 53, no. 8, pp. 1915-1942.

Barbatti, M, Ruckenbauer, M, Plasser, F, Pittner, J, Granucci, G, Persico, M & Lischka, H 2014, ‘Newton‐X: a surface‐hopping program for nonadiabatic molecular dynamics’, Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 4, no. 1, pp. 26-33.

Belcher, J, McLean, KJ, Matthews, S, Woodward, LS, Fisher, K, Rigby, SE, Nelson, DR, Potts, D, Baynham, MT, Parker, DA & Leys, D 2014, ‘Structure and biochemical properties of the alkene producing cytochrome P450 OleTJE (CYP152L1) from the Jeotgalicoccus sp. 8456 bacterium’, Journal of Biological Chemistry, vol. 289, no. 10, pp. 6535-6550.

Beller, HR, Lee, TS & Katz, L 2015, ‘Natural products as biofuels and bio-based chemicals: fatty acids and isoprenoids’, Natural Product Reports, vol. 32, no. 10, pp. 1508-1526.

Bellissent-Funel, MC. Hassanali, A, Havenith, M, Henchman, R, Pohl, P, Sterpone, F, van der Spoel, D, Xu, Y & Garcia, AE 2016, ‘Water determines the structure and dynamics of proteins’, Chemical Reviews, vol. 116, no. 13, pp. 7673-7697.

Cabeza de Vaca, I, Acebes, S & Guallar, V 2016, ‘Ecoupling server: a tool to compute and analyze electronic couplings’, Journal of Computational Chemistry, vol. 37, no. 18, pp. 1740-1745.

Calbo, J, Orti, E, Sancho-García, JC & Arago, J 2015, ‘Accurate treatment of large supramolecular complexes by double-hybrid density functionals coupled with nonlocal van der Waals corrections’, Journal of Chemical Theory and Computation, vol. 11, no. 3, pp. 932-939.

Cao, L & Ryde, U 2018, ‘On the difference between additive and subtractive QM/MM calculations’, Frontiers in Chemistry, vol. 6, pp. 1-15.

Carvalho, AKF, da Conceição, LRV, Silva, JPV, Perez, VH & de Castro, HF 2017, ‘Biodiesel production from Mucor circinelloides using ethanol and heteropolyacid in one and two-step transesterification’, Fuel, vol. 202, pp. 503-511.

Chen, B, Lee, DY & Chang, MW 2015, ‘Combinatorial metabolic engineering of Saccharomyces cerevisiae for terminal alkene production’, Metabolic Engineering, vol. 31, pp. 53-61.

Chen, X, Zhou, L, Tian, K, Kumar, A, Singh, S, Prior, BA & Wang, Z 2013, ‘Metabolic engineering of Escherichia coli: a sustainable industrial platform for bio-based chemical production’, Biotechnology Advances, vol. 31, no. 8, pp. 1200-1223.

Choi, YJ & Lee, SY 2013, ‘Microbial production of short-chain alkanes, Nature, vol. 502, no. 7472, pp. 571-574.

Chum, HL, Warner, E, Seabra, JE & Macedo, IC 2014, ‘A comparison of commercial ethanol production systems from Brazilian sugarcane and US corn’, Biofuels, Bioproducts and Biorefining, vol. 8, no. 2, pp. 205-223.

Cornish-Bowden, A 2014, Principles of enzyme kinetics, Elsevier, New York, NY.

Covert, T, Greenstone, M & Knittel, CR 2016, ‘Will we ever stop using fossil fuels?’ Journal of Economic Perspectives, vol. 30, no. 1, pp. 117-138.

de Visser, SP 2009, ‘Introduction to quantum behaviour – a primer,’ in RK Allemann & NS Scrutton (eds), Quantum tunnelling in enzyme-catalysed reactions, Royal Society of Chemistry, Cambridge, UK, pp. 18-35.

de Visser, SP, Quesne, MG, Martin, B, Comba, P & Ryde, U 2014, ‘Computational modelling of oxygenation processes in enzymes and biomimetic model complexes’, Chemical Communications, vol. 50, no. 3, pp. 262-282.

Dennig, A, Kuhn, M, Tassoti, S, Thiessenhusen, A, Gilch, S, Bülter, T, Haas, T, Hall, M & Faber, K 2015, ‘Oxidative decarboxylation of short‐chain fatty acids to 1‐alkenes’, Angewandte Chemie International Edition, vol. 54, no. 30, pp. 8819-8822.

Dohm, S, Spohr, E & Korth, M 2017, ‘Developing adaptive QM/MM computer simulations for electrochemistry’, Journal of Computational Chemistry, vol. 38, no. 1, pp. 51-58.

Duarte, F, Amrein, BA, Blaha-Nelson, D & Kamerlin, SC 2015, ‘Recent advances in QM/MM free energy calculations using reference potentials’, Biochimica et Biophysica Acta (BBA)-General Subjects, vol. 1850, no. 5, pp. 954-965.

Edmonds, AR 2016, Angular momentum in quantum mechanics, Princeton University Press, Princeton.

Fang, B, Xu, H, Liu, Y, Qi, F, Zhang, W, Chen, H, Wang, C, Wang, Y, Yang, W & Li, S 2017, ‘Mutagenesis and redox partners analysis of the P450 fatty acid decarboxylase OleT JE’, Scientific Reports, vol. 7, p. 44258.

Faponle, AS, Quesne, MG & de Visser, SP 2016, ‘Origin of the regioselective fatty‐acid hydroxylation versus decarboxylation by a cytochrome P450 Peroxygenase: what drives the reaction to biofuel production?’ Chemistry-A European Journal, vol. 22, no. 16, pp. 5478-5483.

Fu, WJ, Chi, Z, Ma, ZC, Zhou, HX, Liu, GL, Lee, CF & Chi, ZM 2015, ‘Hydrocarbons, the advanced biofuels produced by different organisms, the evidence that alkanes in petroleum can be renewable’, Applied Microbiology and Biotechnology, vol. 99, no. 18, pp. 7481-7494.

Genheden, S, Reymer, A, Saenz-Méndez, P & Eriksson, LA 2017, ‘Computational chemistry and molecular modelling basics’, Computational Tools for Chemical Biology, vol. 3, pp. 1-38.

Ghosh, A, Khanra, S, Mondal, M, Halder, G, Tiwari, ON, Saini, S, Bhowmick, TK & Gayen, K 2016, ‘Progress toward isolation of strains and genetically engineered strains of microalgae for production of biofuel and other value added chemicals: a review’, Energy Conversion and Management, vol. 113, pp. 104-118.

Grant, JL, Hsieh, CH & Makris, TM 2015, ‘Decarboxylation of fatty acids to terminal alkenes by cytochrome P450 compound I’, Journal of the American Chemical Society, vol. 137, no. 15, pp. 4940-4943.

Grant, JL, Mitchell, ME & Makris, TM 2016, ‘Catalytic strategy for carbon− carbon bond scission by the cytochrome P450 OleT’, Proceedings of the National Academy of Sciences, vol. 113, no. 36, pp. 10049-10054.

Gricman, L, Vogel, C & Pleiss, J 2014, ‘Conservation analysis of class‐specific positions in cytochrome P450 monooxygenases: functional and structural relevance’, Proteins: Structure, Function, and Bioinformatics, vol. 82, no. 3, pp. 491-504.

Grimme, S & Schreiner, PR 2018, ‘Computational chemistry: the fate of current methods and future challenges’, Angewandte Chemie International Edition, vol. 57, no. 16, pp. 4170-4176.

Guo, M, Song, W & Buhain, J 2015, ‘Bioenergy and biofuels: history, status, and perspective,’ Renewable and Sustainable Energy Reviews, vol. 42, pp. 712-725.

Hammerer, L, Winkler, CK & Kroutil, W 2018, ‘Regioselective biocatalytic hydroxylation of fatty acids by cytochrome P450s’, Catalysis Letters, vol. 148, no. 3, pp. 787-812.

Harrison, P & Valavanis, A 2016, Quantum wells, wires and dots: theoretical and computational physics of semiconductor nanostructures, John Wiley & Sons, Hoboken, NJ.

Hayashi, S, Uchida, Y, Hasegawa, T, Higashi, M, Kosugi, T & Kamiya, M 2017, ‘QM/MM geometry optimization on extensive free-energy surfaces for examination of enzymatic reactions and design of novel functional properties of proteins’, Annual Review of Physical Chemistry, vol. 68, pp. 135-154.

Hendry, JI, Prasannan, CB, Joshi, A, Dasgupta, S & Wangikar, PP 2016, ‘Metabolic model of Synechococcus sp. PCC 7002: prediction of flux distribution and network modification for enhanced biofuel production’, Bioresource Technology, vol. 213, pp. 190-197.

Hsieh, CH & Makris, TM 2016, ‘Expanding the substrate scope and reactivity of cytochrome P450 OleT’, Biochemical and Biophysical Research Communications, vol. 476, no. 4, pp. 462-466.

Hsieh, CH, Huang, X, Amaya, JA, Rutland, CD, Keys, CL, Groves, JT, Austin, RN & Makris, TM 2017, ‘The enigmatic P450 decarboxylase OleT is capable of, but evolved to frustrate, oxygen rebound chemistry’, Biochemistry, vol. 56, no. 26, pp. 3347-3357.

Hsieh, CH, Huang, X, Amaya, JA, Rutland, CD, Keys, CL, Groves, JT, Austin, RN & Makris, TM 2017, ‘The enigmatic P450 decarboxylase OleT is capable of, but evolved to frustrate, oxygen rebound chemistry’, Biochemistry, vol. 56, no. 26, pp. 3347-3357.

International Energy Agency 2018, Global energy demand grew by 2.1% in 2017, and carbon emissions rose for the first time since 2014. Web.

Jensen, F 2017, Introduction to computational chemistry, John Wiley & Sons, Hoboken, NJ.

Ji, L, Faponle, AS, Quesne, MG, Sainna, MA, Zhang, J, Franke, A, Kumar, D, van Eldik, R, Liu, W & de Visser, SP 2015, ‘Drug metabolism by cytochrome P450 enzymes: what distinguishes the pathways leading to substrate hydroxylation over desaturation?’ Chemistry-A European Journal, vol. 21, no. 25, pp. 9083-9092.

Jindal, G & Warshel, A 2016, ‘Exploring the dependence of QM/MM calculations of enzyme catalysis on the size of the QM region’, The Journal of Physical Chemistry B, vol. 120, no. 37, pp. 9913-9921.

Jones, RO 2015, ‘Density functional theory: its origins, rise to prominence, and future’, Reviews of Modern Physics, vol. 87, no. 3, pp. 897-923.

Karelina, M & Kulik, HJ 2017, ‘Systematic quantum mechanical region determination in QM/MM simulation’, Journal of Chemical Theory and Computation, vol. 13, no. 2, pp. 563-576.

Koch, W & Holthausen, MC 2015, A chemist’s guide to density functional theory, John Wiley & Sons, Hoboken, NJ.

Kulik, HJ, Zhang, J, Klinman, JP & Martinez, TJ 2016, ‘How large should the QM region be in QM/MM calculations? The case of catechol O-methyltransferase’, The Journal of Physical Chemistry B, vol. 120, no. 44, pp. 11381-11394.

Le Roy, RJ 2017, ‘LEVEL: A computer program for solving the radial Schrödinger equation for bound and quasibound levels’, Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 186, pp. 167-178.

Lennen, RM & Pfleger, BF 2013, Microbial production of fatty acid-derived fuels and chemicals, Current Opinion in Biotechnology, vol.24, no. 6, pp. 1044-1053.

Lewars, EG 2016, Computational chemistry: introduction to the theory and applications of molecular and quantum mechanics, Springer, New York.

Liakos, DG & Neese, F 2015, ‘Domain based pair natural orbital coupled cluster studies on linear and folded alkane chains’, Journal of Chemical Theory and Computation, vol. 11, no. 5, pp. 2137-2143.

Liang, Y, Wei, J, Qiu, X & Jiao, N 2018, ‘Homogeneous oxygenase catalysis’, Chemical Reviews, vol. 118, no. 10, pp. 4912-4945.

Liao, JC, Mi, L, Pontrelli, S & Luo, S 2016, ‘Fuelling the future: microbial engineering for the production of sustainable biofuels’, Nature Reviews Microbiology, vol. 14, no. 5, pp. 288-304.

Lin, FM, Marsh, ENG & Lin, XN 2015, ‘Recent progress in hydrocarbon biofuel synthesis: pathways and enzymes’, Chinese Chemical Letters, vol. 26, no. 4, pp. 431-434.

Liu, F, Luehr, N, Kulik, HJ & Martínez, TJ 2015, ‘Quantum chemistry for solvated molecules on graphical processing units using polarizable continuum models, Journal of Chemical Theory and Computation, vol. 11, no. 7, pp. 3131-3144.

Liu, Y, Wang, C, Yan, J, Zhang, W, Guan, W, Lu, X & Li, S 2014, ‘Hydrogen peroxide-independent production of α-alkenes by OleT JE P450 fatty acid decarboxylase’, Biotechnology for Biofuels, vol. 7, no. 1, p. 28.

Mardirossian, N & Head-Gordon, M 2017, ‘Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals’, Molecular Physics, vol. 115, no. 19, pp. 2315-2372.

Marisch, K, Bayer, K, Scharl, T, Mairhofer, J, Krempl, PM, Hummel, K, Razzazi-Fazeli, E & Striedner, G 2013, ‘A comparative analysis of industrial Escherichia coli K–12 and B strains in high-glucose batch cultivations on process-, transcriptome-and proteome level’, PloS One, vol. 8, no. 8, pp. 1-16.

Martins, AL, Wanke, P, Chen, Z & Zhang, N 2018, ‘Ethanol production in Brazil: an assessment of main drivers with MCMC generalized linear mixed models’, Resources Conservation and Recycling, vol. 132, pp. 16-27.

Mathews, J & Ardiyanto, A 2015, ‘Estimation of greenhouse gas emissions for palm oil biodiesel production: a review and case study within the Council Directives 2009/28/EC of the European Parliament’, Journal of Oil Palm Environment & Health, vol. 6, pp. 25-41.

Matthews, S, Belcher, JD, Tee, KL, Girvan, HM, McLean, KJ, Rigby, SE, Levy, CW, Leys, D, Parker, DA, Blankley, RT & Munro, AW 2017a, ‘Catalytic determinants of alkene production by the cytochrome P450 peroxygenase OleTJE’, Journal of Biological Chemistry, vol. 292, no. 12, pp. 5128-5143.

Matthews, S, Tee, KL, Rattray, NJ, McLean, KJ, Leys, D, Parker, DA, Blankley, RT & Munro, AW 2017b, ‘Production of alkenes and novel secondary products by P450 OleTJE using novel H2O2‐generating fusion protein systems’, FEBS Letters, vol. 591, no. 5, pp. 737-750.

Maxwell, P, Pendás, ÁM & Popelier, PL 2016, ‘Extension of the interacting quantum atoms (IQA) approach to B3LYP level density functional theory (DFT)’, Physical Chemistry Chemical Physics, vol. 18, no. 31, pp. 20986-21000.

McMahon, D 2013, Quantum mechanics demystified (2nd ed.), McGraw-Hill, New York, NY.

Medvedev, MG, Bushmarinov, IS, Sun, J, Perdew, JP & Lyssenko, KA 2017, ‘Density functional theory is straying from the path toward the exact functional’, Science, vol. 355, no. 6320, pp. 49-52.

Melander, M, Jónsson, EO, Mortensen, JJ, Vegge, T & García Lastra, JM 2016, ‘Implementation of constrained DFT for computing charge transfer rates within the projector augmented wave method’, Journal of Chemical Theory and Computation, vol. 12, no. 11, pp. 5367-5378.

Mendez-Perez, D, Herman, NA & Pfleger, BF 2014, ‘A desaturase gene involved in formation of 1, 14-nonadecadiene in Synechococcus sp. strain PCC 7002’, Applied and Environmental Microbiology, vol. 80, no. 19, pp. 6073-6079.

Metz, S, Kästner, J, Sokol, AA, Keal, TW & Sherwood, P 2014, ‘C hem S hell—a modular software package for QM/MM simulations’, Wiley Interdisciplinary Reviews: Computational Molecular Science, vol. 4, no. 2, pp. 101-110.

Mofijur, M, Rasul, MG, Hyde, J, Azad, AK, Mamat, R & Bhuiya, MMK 2016, ‘Role of biofuel and their binary (diesel–biodiesel) and ternary (ethanol–biodiesel–diesel) blends on internal combustion engines emission reduction’, Renewable and Sustainable Energy Reviews, vol. 53, pp. 65-278.

Mohebbi, A, Abbaszadeh, M & Dehghan, M 2013, ‘The use of a meshless technique based on collocation and radial basis functions for solving the time fractional nonlinear Schrödinger equation arising in quantum mechanics’, Engineering Analysis with Boundary Elements, vol. 37, no. 2, pp. 475-485.

Munro, AW, McLean, KJ, Grant, JL & Makris, TM 2018, ‘Structure and function of the cytochrome P450 peroxygenase enzymes’, Biochemical Society Transactions, vol. 46, no. 1, pp. 183-196.

Nebert, DW, Wikvall, K & Miller, WL 2013, ‘Human cytochromes P450 in health and disease’, Philosophical Transactions of the Royal Society B, vol. 368, no. 1612, pp. 1-21.

Nelson, DR 2018, ‘Cytochrome P450 diversity in the tree of life’, Biochimica et Biophysica Acta (BBA)-Proteins and Proteomics, vol. 1866, no. 1, pp. 141-154.

Nicoletti, G, Arcuri, N, Nicoletti, G & Bruno, R 2015, ‘A technical and environmental comparison between hydrogen and some fossil fuels’, Energy Conversion and Management, vol. 89, pp. 205-213.

Ning, H & Simos, TE 2015, ‘A low computational cost eight algebraic order hybrid method with vanished phase-lag and its first, second, third and fourth derivatives for the approximate solution of the Schrödinger equation’, Journal of Mathematical Chemistry, vol. 53, no. 6, pp. 1295-1312.

Obot, IB, Macdonald, DD & Gasem, ZM 2015, ‘Density functional theory (DFT) as a powerful tool for designing new organic corrosion inhibitors. Part 1: an overview’, Corrosion Science, vol. 99, pp. 1-30.

Omnes, R 2018, The interpretation of quantum mechanics, Princeton University Press, Princeton.

Pai, ZP, Selivanova, NV, Oleneva, PV, Berdnikova, PV & Beskopyl’nyi, AM 2017, ‘Catalytic oxidation of α-alkenes with hydrogen peroxide to carboxylic acids in the presence of peroxopolyoxotungstate complexes’, Catalysis Communications, vol. 88, pp. 45-49.

Piana, S, Donchev, AG, Robustelli, P & Shaw, DE 2015, ‘Water dispersion interactions strongly influence simulated structural properties of disordered protein states’, The Journal of Physical Chemistry B, vol. 119, no. 16, pp. 5113-5123.

Pulidindi, IN & Gedanken, A 2015, ‘Employing novel techniques (microwave and sonochemistry) in the synthesis of biodiesel and bioethanol’, in Z Fang, RL Smith & Qi, X (eds), Production of biofuels and chemicals with ultrasound, Springer, Dordrecht, pp. 159-185.

Quesne, MG, Borowski, T & de Visser, SP 2016, ‘Quantum mechanics/molecular mechanics modelling of enzymatic processes: caveats and breakthroughs,’ Chemistry–A European Journal, vol. 22, no. 8, pp. 2562-2581.

Renault, H, Bassard, JE, Hamberger, B & Werck-Reichhart, D 2014, ‘Cytochrome P450-mediated metabolic engineering: current progress and future challenges’, Current Opinion in Plant Biology, vol. 19, pp. 27-34.

Rendic, S & Guengerich, FP 2014, ‘Survey of human oxidoreductases and cytochrome P450 enzymes involved in the metabolism of xenobiotic and natural chemicals’, Chemical Research in Toxicology, vol. 28, no. 1, pp. 38-42.

Řezáč, J 2016, ‘Cuby: an integrative framework for computational chemistry’, Journal of Computational Chemistry, vol. 37, no. 13, pp. 1230-1237.

Rosano, GL & Ceccarelli, EA 2014, ‘Recombinant protein expression in Escherichia coli: advances and challenges’, Frontiers in Microbiology, vol. 5, pp. 1-17.

Roßbach, S & Ochsenfeld, C 2017, ‘Influence of coupling and embedding schemes on QM size convergence in QM/MM approaches for the example of a proton transfer in DNA’, Journal of Chemical Theory and Computation, vol. 13, no. 3, pp. 1102-1107.

Saika, H, Horita, J, Taguchi-Shiobara, F, Nonaka, S, Ayako, NY, Iwakami, S, Hori, K, Matsumoto, T, Tanaka, T, Itoh, T & Yano, M 2014, ‘A novel rice cytochrome P450 gene, CYP72A31, confers tolerance to acetolactate synthase-inhibiting herbicides in rice and Arabidopsis’, Plant Physiology, vol. 16, no. 33, pp. 1232-1240.

Shahir, SA, Masjuki, HH, Kalam, MA, Imran, A & Ashraful, AM 2015, ‘Performance and emission assessment of diesel–biodiesel–ethanol/bioethanol blend as a fuel in diesel engines: a review’, Renewable and Sustainable Energy Reviews, vol. 48, pp. 62-78.

Sparta, M & Neese, F 2014, ‘Chemical applications carried out by local pair natural orbital based coupled-cluster methods’, Chemical Society Reviews, vol. 43, no. 14, pp. 5032-5041.

Wang, B, Usharani, D, Li, C & Shaik, S 2014, ‘Theory uncovers an unusual mechanism of DNA repair of a lesioned adenine by AlkB enzymes,’ Journal of the American Chemical Society, vol. 136, no. 39, pp. 13895-13901.

Wang, J & Zhu, K 2018, ‘Microbial production of alka (e) ne biofuels’, Current Opinion in Biotechnology, vol. 50, pp. 11-18.

Wang, W, Liu, X & Lu, X 2013, ‘Engineering cyanobacteria to improve photosynthetic production of alka(e)nes’, Biotechnology for Biofuels, vol. 6, no. 1, pp. 1-9.

Wang, Y, Lan, D, Durrani, R & Hollmann, F 2017, ‘Peroxygenases en route to becoming dream catalysts. What are the opportunities and challenges?’ Current Opinion in Chemical Biology, vol. 37, pp. 1-9.

Xu, H, Ning, L, Yang, W, Fang, B, Wang, C, Wang, Y, Xu, J, Collin, S, Laeuffer, F, Fourage, L & Li, S 2017, ‘In vitro oxidative decarboxylation of free fatty acids to terminal alkenes by two new P450 peroxygenases’, Biotechnology for Biofuels, vol. 10, no. 1, p. 208.

Xu, H, Ning, L, Yang, W, Fang, B, Wang, C, Wang, Y, Xu, J, Collin, S, Laeuffer, F, Fourage, L & Li, S 2017, ‘In vitro oxidative decarboxylation of free fatty acids to terminal alkenes by two new P450 peroxygenases’, Biotechnology for Biofuels, vol. 10, no. 1, pp. 1-15.

Yanai, T, Fann, GI, Beylkin, G & Harrison, RJ 2015, ‘Multiresolution quantum chemistry in multiwavelet bases: excited states from time-dependent Hartree–Fock and density functional theory via linear response’, Physical Chemistry Chemical Physics, vol. 17, no. 47, pp. 31405-31416.

Zanger, UM & Schwab, M 2013, ‘Cytochrome P450 enzymes in drug metabolism: regulation of gene expression, enzyme activities, and impact of genetic variation’, Pharmacology & Therapeutics, vol. 138, no. 1, pp. 103-141.

Zargar, A, Bailey, CB, Haushalter, RW, Eiben, CB, Katz, L & Keasling, JD 2017, ‘Leveraging microbial biosynthetic pathways for the generation of ‘drop-in’ biofuels’, Current Opinion in Biotechnology, vol. 45, pp. 156-163.

Zeida, A, Reyes, AM, Lebrero, MC, Radi, R, Trujillo, M & Estrin, DA 2014, ‘The extraordinary catalytic ability of peroxiredoxins: a combined experimental and QM/MM study on the fast thiol oxidation step,’ Chemical Communications, vol. 50, no. 70, pp. 10070-10073.

Zhao, YJ, Cheng, QQ, Su, P, Chen, X, Wang, XJ, Gao, W & Huang, LQ 2014, ‘Research progress relating to the role of cytochrome P450 in the biosynthesis of terpenoids in medicinal plants’, Applied Microbiology and Biotechnology, vol. 98, no. 6, pp. 2371-2383.

Zhou, YJ, Hu, Y, Zhu, Z, Siewers, V & Nielsen, J 2018, ‘Engineering 1-alkene biosynthesis and secretion by dynamic regulation in yeast’, ACS Synthetic Biology, vol. 7, no. 2, pp. 584-590.