Theoretical Analysis of the Reaction Mechanism of Biotin Carboxylase

Human Acetyl-CoA Carboxylase 1 Is an Isomerase: Carboxyl Transfer Is Activated by Catalytic Effect of Isomerization

Obesity and its related diseases such as cancer and diabetes are leading life-threatening issues in the modern world. Thus, new drugs toward obesity and obesity-caused diseases are highly desired. Human acetyl-CoA carboxylase 1 (hACC1) in charge of the rate-limiting step of the human fatty acid synthesis was recognized as an attractive target for rational drug design. The fundamental reaction mechanism and nature of the transition state of hACC1 remain unclear. In this study, combined quantum mechanics and molecular mechanics (QM/MM), molecular dynamics (MD), and free-energy simulations were performed to investigate the catalytic mechanism of the hACC1-catalyzed carboxyl-transfer reaction. Our computational results show a three-step mechanism for carboxyl transferase (CT)-catalyzed reaction, including isomerization of carboxybiotin, proton-transfer from acetyl-CoA to carboxybiotin, and carboxylation of acetyl-CoA enolate. Interestingly, isomerization of carboxybiotin is the rate-limiting step of the entire reaction pathway, indicating hACC1 has the catalytic effect of isomerization and thus might be an isomerase also. The activation free-energy barrier of carboxyl-transfer catalyzed by hACC1 was calculated to be 16.4 kcal/mol, in excellent agreement with the experimental result (16.7 kcal/mol). The obtained reaction mechanism together with the nature of the transition state provides helpful knowledge not only for future investigation of other ACCs but also for rational design of hACC1 inhibitors, such as TS analogue. The catalytic effect of hACC1 isomerization is discussed.

Four amino acids define the CO2 binding pocket of enoyl-CoA carboxylases/reductases

Carboxylases capture and convert CO₂, which makes them key enzymes in photosynthesis and the global carbon cycle. However, the question how enzymes bind atmospheric CO₂ is still unsolved. We studied enoyl-CoA carboxylases/reductases (Ecrs), the fastest CO₂-fixing enzymes in nature, using structural biology, biochemistry, and advanced computational methods. Ecrs create a highly specific CO₂-binding pocket with 4 amino acids at the active site. The pocket controls the fate of the gaseous molecule during catalysis and shields the catalytic center from oxygen and water. This exquisite control makes Ecrs highly efficient carboxylases outcompeting RuBisCO, the key enzyme of photosynthesis, by an order of magnitude. Our findings define the atomic framework for the future development of CO₂-converting catalysts in biology and chemistry.

Carboxylases are biocatalysts that capture and convert carbon dioxide (CO₂) under mild conditions and atmospheric concentrations at a scale of more than 400 Gt annually. However, how these enzymes bind and control the gaseous CO₂ molecule during catalysis is only poorly understood. One of the most efficient classes of carboxylating enzymes are enoyl-CoA carboxylases/reductases (Ecrs), which outcompete the plant enzyme RuBisCO in catalytic efficiency and fidelity by more than an order of magnitude. Here we investigated the interactions of CO₂ within the active site of Ecr from Kitasatospora setae. Combining experimental biochemistry, protein crystallography, and advanced computer simulations we show that 4 amino acids, N81, F170, E171, and H365, are required to create a highly efficient CO₂-fixing enzyme. Together, these 4 residues anchor and position the CO₂ molecule for the attack by a reactive enolate created during the catalytic cycle. Notably, a highly ordered water molecule plays an important role in an active site that is otherwise carefully shielded from water, which is detrimental to CO₂ fixation. Altogether, our study reveals unprecedented molecular details of selective CO₂ binding and C–C-bond formation during the catalytic cycle of nature’s most efficient CO₂-fixing enzyme. This knowledge provides the basis for the future development of catalytic frameworks for the capture and conversion of CO₂ in biology and chemistry.

Behavior of the ATP grasp domain of biotin carboxylase monomers and dimers studied using molecular dynamics simulations

The enzyme biotin carboxylase (BC) uses adenosine triphosphate (ATP) to carboxylate biotin and is involved in fatty acid synthesis. Structural evidence suggests that the B domain of BC undergoes a large hinge motion of ∼45° when binding and releasing substrates. Escherichia coli BC can function as a natural homodimer and as a mutant monomer. Using molecular dynamics simulations, we evaluate the free energy profile along a closure angle of the B domain of E. coli BC for three cases: a monomer without bound Mg(2)ATP, a monomer with bound Mg(2)ATP, and a homodimer with bound Mg(2)ATP in one subunit. The simulation results show that a closed state is the most probable for the monomer with or without bound Mg(2)ATP. For the dimer with Mg(2)ATP in one of its subunits, communication between the two subunits was observed. Specifically, in the dimer, the opening of the subunit without Mg(2)ATP caused the other subunit to open, and hysteresis was observed upon reclosing it. The most stable state of the dimer is one in which the B domain of both subunits is closed; however, the open state for the B domain without Mg(2)ATP is only approximately 2k(B)T higher in free energy than the closed state. A simple diffusion model indicates that the mean times for opening and closing of the B domain in the monomer with and without Mg(2)ATP are much smaller than the overall reaction time, which is on the order of seconds.

Quantum Chemical Modeling of Enzymatic Reactions: The Case of Decarboxylation

We present a systematic study of the decarboxylation step of the enzyme aspartate decarboxylase with the purpose of assessing the quantum chemical cluster approach for modeling this important class of decarboxylase enzymes. Active site models ranging in size from 27 to 220 atoms are designed, and the barrier and reaction energy of this step are evaluated. To model the enzyme surrounding, homogeneous polarizable medium techniques are used with several dielectric constants. The main conclusion is that when the active site model reaches a certain size, the solvation effects from the surroundings saturate. Similar results have previously been obtained from systematic studies of other classes of enzymes, suggesting that they are of a quite general nature.

Umbrella sampling simulations of biotin carboxylase: is a structure with an open ATP grasp domain stable in solution?

Biotin carboxylase is a homodimer that utilizes ATP to carboxylate biotin. Studies of the enzyme using X-ray crystallography revealed a prominent conformational change upon binding ATP. To determine the importance of this closing motion, the potential of mean force with the closure angle as a reaction coordinate was calculated using molecular dynamics simulations and umbrella sampling for a monomer of Escherichia coli biotin carboxylase in water with restraints to simulate attachment to a surface. The result suggests that the most stable state for the enzyme is a closed state different from both the ATP-bound and open state X-ray crystallography structures. There is also a significant motion of a region near the dimer interface not predicted by considering only open and closed configurations, which may have implications for the dynamics and activity of the dimer.