Molecular dynamics simulation of the last step of a catalytic cycle: product release from the active site of the enzyme chorismate mutase from Mycobacterium tuberculosis

Geometric Remodeling of Nitrilase Active Pocket Based on ALF-Scanning Strategy To Enhance Aromatic Nitrile Substrate Preference and Catalytic Efficiency

Nitrilase can catalyze nitrile compounds to generate corresponding carboxylic acids. Nitrilases as promiscuous enzymes can catalyze a variety of nitrile substrates, such as aliphatic nitriles, aromatic nitriles, etc. However, researchers tend to prefer enzymes with high substrate specificity and high catalytic efficiency. In this study, we developed an active pocket remodeling (ALF-scanning) based on modulating the geometry of the nitrilase active pocket to alter substrate preference and improve catalytic efficiency. Using this strategy, combined with site-directed saturation mutagenesis, we successfully obtained 4 mutants with strong aromatic nitrile preference and high catalytic activity, W170G, V198L, M197F, and F202M, respectively. To explore the synergistic relationship of these 4 mutations, we constructed 6 double-combination mutants and 4 triple-combination mutants. By combining mutations, we obtained the synergistically enhanced mutant V198L/W170G, which has a significant preference for aromatic nitrile substrates. Compared with the wild type, its specific activities for 4 aromatic nitrile substrates are increased to 11.10-, 12.10-, 26.25-, and 2.55-fold, respectively. By mechanistic dissection, we found that V198L/W170G introduced a stronger substrate-residue π-alkyl interaction in the active pocket and obtained a larger substrate cavity (225.66 Å3 to 307.58 Å3), making aromatic nitrile substrates more accessible to be catalyzed by the active center. Finally, we conducted experiments to rationally design the substrate preference of 3 other nitrilases based on the substrate preference mechanism and also obtained the corresponding aromatic nitrile substrate preference mutants of these three nitrilases and these mutants with greatly improved catalytic efficiency. Notably, the substrate range of SmNit is widened. IMPORTANCE In this study, the active pocket was largely remodeled based on the ALF-scanning strategy we developed. It is believed that ALF-scanning not only could be employed for substrate preference modification but might also play a role in protein engineering of other enzymatic properties, such as substrate region selectivity and substrate spectrum. In addition, the mechanism of aromatic nitrile substrate adaptation we found is widely applicable to other nitrilases in nature. To a large extent, it could provide a theoretical basis for the rational design of other industrial enzymes.

Computational investigations of allostery in aromatic amino acid biosynthetic enzymes

Allostery, in which binding of ligands to remote sites causes a functional change in the active sites, is a fascinating phenomenon observed in enzymes. Allostery can occur either with or without significant conformational changes in the enzymes, and the molecular basis of its mechanism can be difficult to decipher using only experimental techniques. Computational tools for analyzing enzyme sequences, structures, and dynamics can provide insights into the allosteric mechanism at the atomic level. Combining computational and experimental methods offers a powerful strategy for the study of enzyme allostery. The aromatic amino acid biosynthesis pathway is essential in microorganisms and plants. Multiple enzymes involved in this pathway are sensitive to feedback regulation by pathway end products and are known to use allostery to control their activities. To date, four enzymes in the aromatic amino acid biosynthesis pathway have been computationally investigated for their allosteric mechanisms, including 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase, anthranilate synthase, chorismate mutase, and tryptophan synthase. Here we review the computational studies and findings on the allosteric mechanisms of these four enzymes. Results from these studies demonstrate the capability of computational tools and encourage future computational investigations of allostery in other enzymes of this pathway.

Structural insights into E. coli porphobilinogen deaminase during synthesis and exit of 1-hydroxymethylbilane

Porphobilinogen deaminase (PBGD) catalyzes the formation of 1-hydroxymethylbilane (HMB), a crucial intermediate in tetrapyrrole biosynthesis, through a step-wise polymerization of four molecules of porphobilinogen (PBG), using a unique dipyrromethane (DPM) cofactor. Structural and biochemical studies have suggested residues with catalytic importance, but their specific role in the mechanism and the dynamic behavior of the protein with respect to the growing pyrrole chain remains unknown. Molecular dynamics simulations of the protein through the different stages of pyrrole chain elongation suggested that the compactness of the overall protein decreases progressively with addition of each pyrrole ring. Essential dynamics showed that domains move apart while the cofactor turn region moves towards the second domain, thus creating space for the pyrrole rings added at each stage. Residues of the flexible active site loop play a significant role in its modulation. Steered molecular dynamics was performed to predict the exit mechanism of HMB from PBGD at the end of the catalytic cycle. Based on the force profile and minimal structural changes the proposed path for the exit of HMB is through the space between the domains flanking the active site loop. Residues reported as catalytically important, also play an important role in the exit of HMB. Further, upon removal of HMB, the structure of PBGD gradually relaxes to resemble its initial stage structure, indicating its readiness to resume a new catalytic cycle.

Heme is the prosthetic group at the core of the oxygen carrier metalloprotein hemoglobin. Heme consists of a tetrapyrrole called porphyrin bound to an iron ion. It is synthesized by the heme biosynthetic pathway, which is common to all eukaryotes and most prokaryotes. Porphobilinogen deaminase, an enzyme in the heme biosynthetic pathway, catalyzes the formation of a linear tetrapyrrole product, 1-hydroxymethylbilane, from four units of porphobilinogen. In this study we carried out molecular dynamics simulations to understand the structural changes that the enzyme undergoes while catalyzing this reaction. There are three segments to the study: 1) understanding the changes in the enzyme when the porphobilinogen units get attached to the dipyrromethane cofactor, thereby forming a polypyrrole chain; 2) exit of the product from the active site of the enzyme via steered molecular dynamics; and 3) the relaxation of the enzyme to the initial stage to resume its catalytic cycle. Molecular dynamics simulations of the protein through the different stages of pyrrole chain elongation gives insight into the motions of domains, active site loop and role of conserved active site residues in facilitating the accommodation of the polypyrrole chain. In addition to this, we propose a possible exit path for the product and demonstrate the relaxation of the enzyme after the exit of the product to resume the catalytic cycle.