Team:SCUT-China/Modeling/Simulation

From 2014.igem.org

Simulation

Reactions and Mechanism

1.Mechanism

The model for enzyme action, first suggested by Brown and Henri but later established more thoroughly Michaelis and Menten, suggests the binding of free enzyme to the substrate forming a enzyme-substrate(ES) complex. This complex undergoes a transformation, releasing product and free enzyme. The free enzyme is then available for another round of binding to new substrate.

In the single substrate enzyme catalyzed reactions, the mechanism is often written as:

And the Michaelis-Menten equation is:

(2) Kinetics of multi substrates enzyme catalyzed reactions

We apply two mechanisms of kinetics of multi substrates enzyme catalyzed reactions to our modeling section, one is Theorell-Chance bi-bi (T-C bi-bi), the other is Ping-pong bi-bi.

Sequential bi-bi reaction means that all substrates must bind to the enzyme before any reaction takes place. Now we take two substrates as an example. In sequential bi-bi reaction, a ternary complex is formed when both substrates bind to the enzyme. But the complex is not steady and two products release after new covalent bonds are formed and old covalent bonds are broken in the complex.

Sequential bi-bi reaction can be random or ordered. The random sequential bi-bi means that any substrate can bind first to the enzyme and any product can release first. The ordered sequential bi-bi means that the substrates bind to the enzyme and the products release in a specific order.

T-C bi-bi reaction is a kind of the ordered sequential bi-bi reactions. Now we take two substrates as an example again. The characteristic of T-C bi-bi is that two products can release in a flash because of the extreme instability of the ternary complex. The allosteric process is not obvious.

Ping-pong bi-bi reaction is a double-displacement reaction. Substrate A bind first to the enzyme followed by product P release. Typically, product P is a fragment of the original substrate A. The rest of the substrate is covalently attached to the enzyme E, which we now designate as F. The second Substrate, B, binds and reacts with enzyme F. Substrate B form a covalent adduct with the covalent fragment of A in the enzyme F to form product Q. The enzyme is finally restored to its initial form, E.

In Ping-pong bi-bi reactions, the mechanism is often written as:

And the equation is:




2.Reactions

Although each domain has its independent function, some reactions are catalyzed by several domains together in the process of polyketide’s synthesis, such as activation of propionyl groups, combination of the extended unit and polyketide chain and so on. Therefore, in analysis of reactions of their mechanism, the introduction will be started from the perspective of the equations of chemical reactions.

(1)Acyl carrier protein, ACP

ACP domain is a carrier. Its function is anchoring the polyketide chains needed to extend. The synthesis of polyketides is engaged on these carriers. The polyketide chains are activated by AT domain, and then are anchored on ACP domains. Until the catalysis of TE domain, polyketide chains are released from ACP domain and become ripe polyketides.

(2)Loading Part

The leading role of the loading part is Acyltransfer-ase, AT domain. In this part, we will introduce the function of AT domain by analyzing the chemical reactions between AT domain, ACP domain and polyketide chains.

The function of AT domain is activating the extended units, propionyl groups which will be combined with the polyketide chains.

In Loading, Proplonyl-CoA reacts with the oxhydryl on the specific binding sites of AT domains and binds to AT domain, and then CoA-SHs are displaced. At this time, proplonyl groups have been activated. They are transferred onto the ACP domains in Loading and anchored there.

In Module 1 or Module 2, AT domains anchor the proplonyl groups of (2S)Methylmalonyl-CoA to the ACP domains in Module 1 or Module 2 by the same mechanism. The proplonyl groups anchored wait for the catalysis of KS domain in the next step.

The concrete equation of chemical reaction is :

After the analysis of equation above, we conclude that this enzyme catalyzed reaction is match with Ping-pong bi-bi mechanism.

The graph of mechanism is written as:

(3)KS Part

The leading role of the KS part is Ketocaylsynthase, KS domain. Its last step is activation of propionyl groups. AT domains in Module 1 or Module 2 activate the propionyl groups in (2S)Methylmalonyl-CoA (details in Loading Part). Afterwards, KS domains combine the extended units anchored on the ACP domains belonging to the Modules which the KS domains belong to with the polyketide chains anchored on the ACP domains belonging to the former Modules (or Loadings). Then new polyketide chains form and are anchored the ACP domains belonging to the Modules which the KS domains belong to.

The function of KS domain is combining the propionyl group with acetyl group at the end of the polyketide chain by forming the C—C bond. In KS domain catalyzed reactions, polyketide chains are transferred from ACP doamains belonging to the former Module (or Loading) onto KS domains, and then occur decarboxylic reaction with the extended units anchored on the ACP domains belonging to the Modules which the KS domains belong to. As a result, the polyketide chains are elongated.

The concrete equation of chemical reaction is :

(4)KR Part

The process of KR domain catalyzed reactions is a little more complex. Its leading role is KR domain. The function of KR domain is deoxidizing the extended unit and enabling it to form the β-hydroxyl ester bond. In KR domain catalyzed reactions, the carbonyl groups of polyketide chains anchored on the ACP domains are combined with the specific binding sites of KR domains through weak chemical bonds. Then the electron cloud in the carbonyl groups is moving towards oxygen atoms so that the carbon atoms have net positive charge. At this time, NADPHs have a chance to bind to the polyketide chains and form instable ternary complexes. In order to reduce the inner energy of this ternary complexes and run towards stable, NADPHs provide hydride ions and KR domains provide protons for the complexes and convert the carbonyl groups to theβ-hydroxyl esters. Finally, the KR domain residues which have lost protons are deoxidized to the initial KR domains by other NADPHs.

The concrete equation of chemical reaction is :

Because KR domains catalyzed reactions are involved several steps, it is inpossible for them to be match with any simple mechanism of enzyme catalyzed reactions. After the analysis of equation above, we define that these enzyme catalyzed reactions may be match with [T-C bi-bi + Ping-pong bi-bi] mechanism. In other words, KR domains catalyzed reactions may belong to Ping-pong bi-bi mechanism which include T-C bi-bi mechanism in its former half part.

The graph of mechanism is written as:

(5)DH Part

The leading role of DH Part is DH domain. The function of DH domain is dehydrating the extended unit and enabling it to form the α, β- enol ester bond. In DH domain catalyzed reactions, the oxhydryls in the extended units and the hydrogen atoms binding to the neighbor carbon atoms are removed by DH domains together and form . After catalysis, α, β-enol ester bonds are formed in the polyketide chains.

The concrete equation of chemical reaction is :

After the analysis of equation above, we conclude that this enzyme catalyzed reaction is match with the mechanism of single substrate enzyme catalyzed reactions.

The graph of mechanism is written as:

(6)ER Part

The leading role of ER Part is ER domain. The function of ER domain is deoxidizing the extended unit and enabling it to form the saturated methylene. The mechanism of ER domain catalyzed reactions is similar to KR domain’s. In ER domain catalyzed reactions, NADPHs provide hydride ions and ER domains provide protons for the complexes and deoxidize carbon-carbon double bond to form the saturated methylene.

The concrete equation of chemical reaction is :

After the analysis of equation above, we conclude that the mechanism of this enzyme catalyzed reaction is similar to KR domain’s.

The graph of mechanism is written as:

(7)TE Part

The leading role of TE Part is TE domain. The function of TE domain is removing the polyketide chain from PKS. In TE domain catalyzed reactions, polyketide chains anchored on the ACP domains react with the oxhydryl on the pecific binding sites of TE domains. Then the polyketide chains are transferred from ACP domains to TE domains. However, polyketide chains are instable when binding to TE domains and they release from TE domains, that is, release from the whole PKS. Finally, they become ripe polyketides.

The concrete equation of chemical reaction is :

After the analysis of equation above, we conclude that this enzyme catalyzed reaction is similar to the single substrate enzyme catalyzed reactions. We confirm that we can apply the mechanism of single substrate enzyme catalyzed reactions to TE domain catalyzed reactions by using King-Altman Method.

The graph of mechanism is written as:






Kinetic Equations and Simplified Model Analysis

In the establishment of kinetic models, we employed the kinetic equations educed from single substrate enzyme catalyzed reactions and multi substrate enzyme catalyzed reactions to each part (details in Reactions and Mechanism). Kinetic equations of enzyme catalyzed reactions (including single substrate’s and multi substrate’s) are educed from Law of Mass Action initially. We wrote out the related kinetic equations according to the mechanism in each part.

1.Loading Part

The reaction mechanism of Loading Part is Ping-pong bi-bi. The kinetic equation is :

KmAis the Michaelis Constant of substrate A when [B] is saturated.

KmBis the Michaelis Constant of substrate B when [A] is saturated.

Acoording to prerequisite of Michaelis-Menten equation, the concentration of substrates are much more than the concentration of enzyme, that Is, [A]>>[B]. Therefore, this kinetic model can be simplified. The simplified equation is:

After perfection, the final equation is:

[A0]is the concentration of DEBS 1.
[B] is the concentration of product in AT domain catalyzed reactions.


2. KS Part
(1) AT Reaction

The reaction mechanism of AT Reaction is same as AT Part’s-------Ping-pong bi-bi. Therefore, the kinetic equation is the same.

(2) KS Reaction

The reaction mechanism of KS Reaction is Ping-pong bi-bi. The kinetic equation is :

KmAis the Michaelis Constant of substrate A when [B] is saturated.
KmBis the Michaelis Constant of substrate B when [A] is saturated.

In KS Reaction, both two substrates (Proplonyl-CoA and (2S)Methylmalonyl-CoA) are anchored on ACP domains. And what's more, ACP domains anchoring substrates belong to the same DEBS 1. We can easily draw a conclusion that the concentration of substrate A equals to substrate B’s, that is, [A]=[B]. Therefore, this kinetic model can be simplified. The simplified equation is:

After perfection, the final equation is:

[A0]is the concentration of DEBS 1.
[F] is the concentration of product in enzyme catalyzed reactions in last step.
[D] is the concentration of product in KS domain catalyzed reactions.

3. KR Part

The reaction mechanism of KR Reaction is Ping-pong bi-bi + T-C bi-bi. It is a little more complex so that the related kinetic equation can’t be written out. Therefore, we simplify the reaction mechanism directly.

Because NADPH and NADP+ is dynamic-balanced in vivo of E.coli and their concentrations are constant, we needn’t consider the amount of NADPH participating in the reactions and the output of NADP+ in the analysis simplifying reaction mechanism. And the graph of mechanism is written as:

However, the amount of NADPH does have a certain influence on KR domain catalyzed reactions, so we should regard NADPH as a kind of inhibitor of KR domain, that is, KR domains will lose efficacy after they provide protons and NADPHs don’t compensate immediately. Then description is similar to enzyme inhibitor’s. Therefore, the final simplified equation is:

[A0]is the concentration of DEBS 1.
[H] is the concentration of NADPH and KR domain.

KIis the dissociation constant of NADPH and KR domain.
[J] is the concentration of product in enzyme catalyzed reactions in last step.
[Q] is the concentration of product in KS domain catalyzed reactions.

4. DH Part

The reaction mechanism of DH Part is single substrate enzyme catalyzed reactions. The kinetic equation is :

[A0]is the concentration of DEBS 1.
[V] is the concentration of product in enzyme catalyzed reactions in last step.
[Y] is the concentration of product in DH domain catalyzed reactions.

5. ER Part

The reaction mechanism of ER Part is same as KR Part’s. Therefore, the kinetic equation is the same.

6. TE Part

The reaction mechanism of TE Part is single substrate enzyme catalyzed reactions. The kinetic equation is :

[A0]is the concentration of DEBS 1.
[W] is the concentration of product in enzyme catalyzed reactions in last step.
[Z] is the concentration of product in TE domain catalyzed reactions.






Reference

[1] Ana Pereda et al. : The loading domain of the erythromycin polyketide synthase is not essential for erythromycin biosynthesis in Saccharopolyspora erythraea, Microbiology (1 998),144, 543-553
[2] David E. Cane: Programming of Erythromycin Biosynthesis by a Modular Polyketide Synthase, J. Biol. Chem. 2010, 285:27517-27523
[3] Alice Y. Chen et al. : Structure-Based Dissociation of a Type I Polyketide Synthase Module, chembiol.2007.05.015
[4] Kira J. Weissman et al. : Evidence for a Protein-Protein Interaction Motif on an Acyl Carrier Protein Domain from a Modular Polyketide Synthase, Chemistry & Biology 13, 625–636, June 2006
[5] Rembert Pieper et al. : Erythromycin Biosynthesis: Kinetic Studies on a Fully Active Modular Polyketide Synthase Using Natural and Unnatural Substrates, Biochemistry 1996, 35, 2054-2060
[6] R.A. Azizyan et al. : Computational Modeling of Kinetics of the Bisubstrate Enzymatic Reaction with Ping-Pong Mechanism, Biolog. Journal of Armenia, 2 (64), 2012
[7] James T. Wu et al. : Stability of NADPH: Effect of Various Factors on the Kinetics of Degradation, CLIN. CHEM. 32/2, 314-319 (1986)