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

Free Energy Profile of Domain Movement in Ligand-Free Citrate Synthase

Citrate synthase plays a fundamental role in the metabolic cycle of the cell. Its catalytic mechanism is complex involving the binding of two substrates that cause a domain movement. In this paper, we used classical molecular dynamics simulations and umbrella-sampling simulations to determine the potential of mean force along a reaction coordinate for the domain movement in ligand-free citrate synthase from pig ( Sus scrofa). The results show that, at 293 K, the closed-domain conformation has a ∼4 k_b T higher energy than the open-domain conformation. In a simple two-state model, this difference means that the enzyme spends 98% of the time in the open-domain conformation ready to receive the substrate, oxaloacetate, rather than the closed-domain conformation where the binding site would be inaccessible to the substrate. Given that experimental evidence indicates that the binding of oxaloacetate induces at least partial closure, this would imply an induced-fit mechanism which we argue is applicable to all enzymes with a functional domain movement for reasons of catalytic efficiency. A barrier of 4 k_b T gives an estimation of the mean first passage time in the range 1-10 μs.

Crystal structure of the essential biotin-dependent carboxylase AccA3 fromMycobacterium tuberculosis

Biotin‐dependent acetyl‐CoA carboxylases catalyze the committed step in type II fatty acid biosynthesis, the main route for production of membrane phospholipids in bacteria, and are considered a key target for antibacterial drug discovery. Here we describe the first structure of AccA3, an essential component of the acetyl‐CoA carboxylase system in Mycobacterium tuberculosis (MTb). The structure, sequence comparisons, and modeling of ligand‐bound states reveal that the ATP cosubstrate‐binding site shows distinct differences compared to other bacterial and eukaryotic biotin carboxylases, including all human homologs. This suggests the possibility to design MTb AccA3 subtype‐specific inhibitors.

Database

Coordinates and structure factors have been deposited in the Protein Data Bank with the accession number 5MLK.

Maurice M, Cleland WW. Novel insights into the biotin carboxylase domain reactions of pyruvate carboxylase from Rhizobium etli

The catalytic mechanism of the MgATP-dependent carboxylation of biotin in the biotin carboxylase domain of pyruvate carboxylase from R. etli (RePC) is common to the biotin-dependent carboxylases. The current site-directed mutagenesis study has clarified the catalytic functions of several residues proposed to be pivotal in MgATP-binding and cleavage (Glu218 and Lys245), HCO(3)(-) deprotonation (Glu305 and Arg301), and biotin enolization (Arg353). The E218A mutant was inactive for any reaction involving the BC domain and the E218Q mutant exhibited a 75-fold decrease in k(cat) for both pyruvate carboxylation and the full reverse reaction. The E305A mutant also showed a 75- and 80-fold decrease in k(cat) for both pyruvate carboxylation and the full reverse reaction, respectively. While Glu305 appears to be the active site base which deprotonates HCO(3)(-), Lys245, Glu218, and Arg301 are proposed to contribute to catalysis through substrate binding interactions. The reactions of the biotin carboxylase and carboxyl transferase domains were uncoupled in the R353M-catalyzed reactions, indicating that Arg353 may not only facilitate the formation of the biotin enolate but also assist in coordinating catalysis between the two spatially distinct active sites. The 2.5- and 4-fold increase in k(cat) for the full reverse reaction with the R353K and R353M mutants, respectively, suggests that mutation of Arg353 allows carboxybiotin increased access to the biotin carboxylase domain active site. The proposed chemical mechanism is initiated by the deprotonation of HCO(3)(-) by Glu305 and concurrent nucleophilic attack on the γ-phosphate of MgATP. The trianionic carboxyphosphate intermediate formed reversibly decomposes in the active site to CO(2) and PO(4)(3-). PO(4)(3-) then acts as the base to deprotonate the tethered biotin at the N(1)-position. Stabilized by interactions between the ureido oxygen and Arg353, the biotin-enolate reacts with CO(2) to give carboxybiotin. The formation of a distinct salt bridge between Arg353 and Glu248 is proposed to aid in partially precluding carboxybiotin from reentering the biotin carboxylase active site, thus preventing its premature decarboxylation prior to the binding of a carboxyl acceptor in the carboxyl transferase domain.

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.