[Structure motif guided mining of MHET hydrolase and development of a two-enzyme cascade for plastics depolymerization at mild temperature].

The utilization of polyethylene terephthalate (PET) has caused significant and prolonged ecological repercussions. Enzymatic degradation is an environmentally friendly approach to addressing PET contamination. Hydrolysis of mono(2-hydroxyethyl) terephthalate (MHET), a competitively inhibited intermediate in PET degradation, is catalyzed by MHET degrading enzymes. Herein, we employed bioinformatic methods that combined with sequence and structural information to discover an MHET hydrolase, BurkMHETase. Enzymatic characterization showed that the enzyme was relatively stable at pH 7.5-10.0 and 30-45 ℃. The kinetic parameters kcat and Km on MHET were (24.2±0.5)/s and (1.8±0.2) µmol/L, respectively, which were similar to that of the well-known IsMHETase with higher substrate affinity. BurkMHETase coupled with PET degradation enzymes improved the degradation of PET films. Structural analysis and mutation experiments indicated that BurkMHETase may have evolved specific structural features to hydrolyze MHET. For MHET degrading enzymes, aromatic amino acids at position 495 and the synergistic interactions between active sites or distal amino acids appear to be required for MHET hydrolytic activity. Therefore, BurkMHETase may have substantial potential in a dual-enzyme PET degradation system while the bioinformatic methods can be used to broaden the scope of applicable MHETase enzymes.

Unveiling the repressive mechanism of a PPS-like regulator (PspR) in polyhydroxyalkanoates biosynthesis network.

Poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV) is a type of polyhydroxyalkanoates (PHA) that exhibits numerous outstanding properties and is naturally synthesized and elaborately regulated in various microorganisms. However, the regulatory mechanism involving the specific regulator PhaR in Haloferax mediterranei, a major PHBV production model among Haloarchaea, is not well understood. In our previous study, we showed that deletion of the phosphoenolpyruvate (PEP) synthetase-like (pps-like) gene activates the cryptic phaC genes in H. mediterranei, resulting in enhanced PHBV accumulation. In this study, we demonstrated the specific function of the PPS-like protein as a negative regulator of phaR gene expression and PHBV synthesis. Chromatin immunoprecipitation (ChIP), in situ fluorescence reporting system, and in vitro electrophoretic mobility shift assay (EMSA) showed that the PPS-like protein can bind to the promoter region of phaRP. Computational modeling revealed a high structural similarity between the rifampin phosphotransferase (RPH) protein and the PPS-like protein, which has a conserved ATP-binding domain, a His domain, and a predicted DNA-binding domain. Key residues within this unique DNA-binding domain were subsequently validated through point mutation and functional evaluations. Based on these findings, we concluded that PPS-like protein, which we now renamed as PspR, has evolved into a repressor capable of regulating the key regulator PhaR, and thereby modulating PHBV synthesis. This regulatory network (PspR-PhaR) for PHA biosynthesis is likely widespread among haloarchaea, providing a novel approach to manipulate haloarchaea as a production platform for high-yielding PHA. KEY POINTS: • The repressive mechanism of a novel inhibitor PspR in the PHBV biosynthesis was demonstrated • PspR is widespread among the PHA accumulating haloarchaea • It is the first report of functional conversion from an enzyme to a trans-acting regulator in haloarchaea.

Computational redesign of a hydrolase for nearly complete PET depolymerization at industrially relevant high-solids loading.

Biotechnological plastic recycling has emerged as a suitable option for addressing the pollution crisis. A major breakthrough in the biodegradation of poly(ethylene terephthalate) (PET) is achieved by using a LCC variant, which permits 90% conversion at an industrial level. Despite the achievements, its applications have been hampered by the remaining 10% of nonbiodegradable PET. Herein, we address current challenges by employing a computational strategy to engineer a hydrolase from the bacterium HR29. The redesigned variant, TurboPETase, outperforms other well-known PET hydrolases. Nearly complete depolymerization is accomplished in 8 h at a solids loading of 200 g kg-1. Kinetic and structural analysis suggest that the improved performance may be attributed to a more flexible PET-binding groove that facilitates the targeting of more specific attack sites. Collectively, our results constitute a significant advance in understanding and engineering of industrially applicable polyester hydrolases, and provide guidance for further efforts on other polymer types.

Traceless enzymatic protein synthesis without ligation sites constraint.

Protein synthesis and semisynthesis offer immense promise for life sciences and have impacted pharmaceutical innovation. The absence of a generally applicable method for traceless peptide conjugation with a flexible choice of junction sites remains a bottleneck for accessing many important synthetic targets, however. Here we introduce the PALME (protein activation and ligation with multiple enzymes) platform designed for sequence-unconstrained synthesis and modification of biomacromolecules. The upstream activating modules accept and process easily accessible synthetic peptides and recombinant proteins, avoiding the challenges associated with preparation and manipulation of activated peptide substrates. Cooperatively, the downstream coupling module provides comprehensive solutions for sequential peptide condensation, cyclization and protein N/C-terminal or internal functionalization. The practical utility of this methodology is demonstrated by synthesizing a series of bioactive targets ranging from pharmaceutical ingredients to synthetically challenging proteins. The modular PALME platform exhibits unprecedentedly broad accessibility for traceless protein synthesis and functionalization, and holds enormous potential to extend the scope of protein chemistry and synthetic biology.

GRAPE, a greedy accumulated strategy for computational protein engineering.

Nature harbors fascinating enzymatic catalysts with high efficiency, chemo-, regio- and stereoselectivity. However, the insufficient stability of the enzymes often prevents their widespread utilization for industrial processes. Not content with the finite repertoire of naturally occurring enzymes, protein engineering holds promises to extend the applications of the improved enzymes with desired physical and catalytic properties. Herein, we devised a computational strategy (greedy accumulated strategy for protein engineering, GRAPE) to enhance the thermostability of enzymes. Through scanning of all point mutations of the structural and evolutionary consensus analysis, a library containing fewer than 100 mutations was established for characterization. After preliminary experimental verification, effective mutations are clustered in a multidimensional physical property space and then accumulated via the greedy algorithm to produce the final designed enzyme. Using the recently reported IsPETase from Ideonella sakaiensis that decomposes PET under ambient temperatures as a starting point, we adopted the GRAPE strategy to come up with a DuraPETase (TM=77°C, raised by 31°C) which showed drastically enhanced degradation performance (300-fold) on semicrystalline PET films at 40°C.

Bioretrosynthesis of Functionalized N-Heterocycles from Glucose via One-Pot Tandem Collaborations of Designed Microbes.

The design of multistrain systems has markedly expanded the prospects of using long biosynthetic pathways to produce natural compounds. However, the cooperative use of artificially engineered microbes to synthesize xenobiotic chemicals from renewable carbohydrates is still in its infancy. Here, a microbial system is developed for the production of high-added-value N-heterocycles directly from glucose. Based on a retrosynthetic analysis, eleven genes are selected, systematically modulated, and overexpressed in three Escherichia coli strains to construct an artificial pathway to produce 5-methyl-2-pyrazinecarboxylic acid, a key intermediate in the production of the important pharmaceuticals Glipizide and Acipimox. Via one-pot tandem collaborations, the designed microbes remarkably realize high-level production of 5-methyl-2-pyrazinecarboxylic acid (6.2 ± 0.1 g L-1) and its precursor 2,5-dimethylpyrazine (7.9 ± 0.7 g L-1). This study is the first application of cooperative microbes for the total biosynthesis of functionalized N-heterocycles and provides new insight into integrating bioretrosynthetic principles with synthetic biology to perform complex syntheses.

Characterization and efficient production of a thermostable, halostable and organic solvent-stable cellulase from an oil reservoir.

The manufacture of biofuels from cellulose is regarded as one of practicable strategies to meet increasing energy demand and alleviate environmental issues. Cellulases, which play an important role in the production of second-generation biofuels, are expected to be highly thermostable, halostable and organic solvent-stable to adapt to the harsh conditions in practical application. Here we cloned and characterized a novel cellulase (MaCel) from Mahella australiensis 50-1 BON, an anaerobic thermophile isolated from an oil reservoir. MaCel exhibited excellent thermostability, halostability as well as organic solvent stability, and could be efficiently produced in a yield of 1.7 × 106 U/L in 15 h with inexpensive culture medium. These results indicate that MaCel may be a suitable candidate for industrial applications, illustrating the potential benefits of enzymes from oil reservoir extremophiles in the manufacture of biofuels.

Reductase of Mutanobactin Synthetase Triggers Sequential C-C Macrocyclization, C-S Bond Formation, and C-C Bond Cleavage.

Mutanobactins (MUBs) and their congeners that contain a macrocycle and/or a thiazepane ring are lipopeptides from Streptococcus mutans, a major causative agent of dental caries. Here we show that the C-terminal reductase domain of MubD releases the lipohexapeptide intermediates in an aldehyde form, which enables a spontaneous C-C macrocyclization. In the presence of a thiol group, the macrocyclized MUBs can further undergo spontaneous C-S bond formation and C-C bond cleavage to generate diverse MUB congeners.

Improving the System Performance of the Asymmetric Biosynthesis of d-Pantoic Acid by Using Artificially Self-Assembled Enzymes in Escherichia coli.

d-Pantoic acid (d-PA) is an important chiral precursor of a broad range of biologically active compounds. The asymmetric synthesis of d-PA through reductase coupling with NADPH regeneration systems is highly promising, but the process is restricted by expensive cofactor consumption and low cofactor recycling frequency. Here, an effective construction of self-assembled ketopantoic acid reductase and glucose dehydrogenase via protein-peptide interaction of PDZ domain and PDZ ligand was established. The self-assembled enzymes exhibited highly ordered two-dimensional threadlike macromolecular structures with improved cofactor regeneration. Furthermore, the bioconversion with whole-cell catalysis showed that the robustness and efficiency of the system with self-assembled enzymes were significantly higher than those of the unassembled enzymes. This study provides a strategy for the effective asymmetric biosynthesis of d-PA with a trace amount of cofactor and shows potential for industrial applications.

Biochemical and structural characterization of a highly active branched-chain amino acid aminotransferase from Pseudomonas sp. for efficient biosynthesis of chiral amino acids.

Aminotransferases (ATs) are important biocatalysts for the synthesis of chiral amines because of their capability of introducing amino group into ketones or keto acids as well as their high enantioselectivity, high regioselectivity. Among all ATs, branched-chain amino acid aminotransferase (BCAT) can use branched-chain amino acids (BCAAs) as substrate, including L-valine, L-leucine, and L-isoleucine, with α-ketoglutarate to form the corresponding α-keto acids and L-glutamate. Alternatively, BCATs have been used for the biosynthesis of unnatural amino acids, such as L-tert-leucine and L-norvaline. In the present study, the BCAT from Pseudomonas sp. (PsBCAT) was cloned and expressed in Escherichia coli for biochemical and structural analyses. The optimal reaction temperature and pH of PsBCAT were 40 °C and 8.5, respectively. PsBCAT exhibited a comparatively broader substrate spectrum and showed remarkably high activity with bulked aliphatic L-amino acids (kcat up to 220 s-1). Additionally, PsBCAT had activities with aromatic L-amino acids, L-histidine, L-lysine, and L-threonine. This substrate promiscuity is unique for the BCAT family and could prove useful in industrial applications. To analyze the catalytic mechanism of PsBCAT with the broad substrate spectrum, the crystal structure of PsBCAT was also determined. Based on the determined crystal structure, we found some differences in the organization of the substrate binding cavity, which may influence the substrate specificity of the enzyme. Finally, conjugated with the ornithine aminotransferase (OrnAT) to shift the reaction equilibrium towards the product formation, the coupled system was applied to the asymmetric synthesis of L-tert-leucine and L-norvaline. In summary, the structural and functional characteristics of PsBCAT were analyzed in detail, and this information will be conducive to industrial production of enantiopure chiral amino acids by aminotransferase.

Thermostability improvement of the glucose oxidase from Aspergillus niger for efficient gluconic acid production via computational design.

Gluconic acid (GA) and its alkali salts are extensively used in the food, feed, beverage, textile, pharmaceutical and construction industries. However, the cost-effective and eco-friendly production of GA remains a challenge. The biocatalytic process involving the conversion of glucose to GA is catalysed by glucose oxidase (GOD), in which the catalytic efficiency is highly dependent on the GOD stability. In this study, we used in silico design to enhance the stability of glucose oxidase from Aspergillus niger. A combination of the best mutations increased the apparent melting temperature by 8.5 °C and significantly enhanced thermostability and thermoactivation. The variant also showed an increased optimal temperature without compromising the catalytic activity at lower temperatures. Moreover, the combined variant showed higher tolerance at pH 6.0 and 7.0, at which the wild-type enzyme rapidly deactivated. For GA production, an approximate 2-fold higher GA production yield was obtained, in which an almost complete conversion of 324 g/L d-glucose to GA was achieved within 18 h. Collectively, this work provides novel and efficient approaches for improving GOD thermostability, and the obtained variant constructed by the computational strategy can be used as an efficient biocatalyst for GA production at industrially viable conditions.

Functional characterization of tryptophan437 at subsite +2 in pullulanase from Bacillus subtilis str. 168.

Pullulanase, a typical debranching enzyme, specifically hydrolyzes α-1,6 glycosidic linkages in pullulan, starch and so on. There is accumulated knowledge about the catalysis of pullulanase, but functional roles of acceptor subsites get less attention. According to the crystal structure of pullulanase from Bacillus subtilis str. 168, stacking interaction between tryptophan437 and glycosyl unit of α-CD was found. Trp437 located in conserved region III was highly conserved in pullulanases. But functional role of this residue remains unclear. Site-directed mutagenesis was used to determine the function of Trp437. Replacement of Trp437 with glycine, phenylalanine, proline, arginine resulted in mutants with undetectable hydrolysis activity and reverse hydrolysis activity. The secondary structure of mutated proteins showed no difference from the wild-type pullulanase, but they lost the capability to bind ß-CD based on the ITC measurement. Molecular docking was performed to investigate the affinity of proteins to ligands including maltotriose, 62-α-D-maltotriosyl-maltotriose and α-CD, showing the affinity of mutants to ligands became weaker compared to that of PulA. And the hydrogen bond between Trp437 and Glu525 in PulA only was found in mutant W437R. These results show that Trp437 at subsite +2 plays crucial role in catalysis of PulA, and it has little tolerance to change.

Molecular dynamics investigations of structural and functional changes in Bcl-2 induced by the novel antagonist BDA-366.

Apoptosis is a fundamental biological phenomenon, in which anti- or proapoptotic proteins of the Bcl-2 family regulate a committed step. Overexpression of Bcl-2, the prototypical antiapoptotic protein in this family, is associated with therapy resistance in various human cancers. Accordingly, Bcl-2 inhibitors intended for cancer therapy have been developed, typically against the BH3 domain. Recent experimental evidences have shown that the antiapoptotic function of Bcl-2 is not immutable, and that BDA-366, a novel antagonist of the BH4 domain, converts Bcl-2 from a survival molecule to an inducer of cell death. In this study, the underlying mechanisms of this functional conversion were investigated by accelerated molecular dynamics simulation. Results revealed that Pro127 and Trp30 in the BH4 domain rotate to stabilize BDA-366 via π-π interactions, and trigger a series of significant conformational changes of the α3 helix. This rearrangement blocks the hydrophobic binding site (HBS) in the BH3 domain and further prevents binding of BH3-only proteins, which consequently allows the BH3-only proteins to activate the proapoptotic proteins. Analysis of binding free energy confirmed that BDA-366 cross-inhibits BH3-only proteins, implying negative cooperative effects across separate binding sites. The newly identified blocked conformation of the HBS along with the open to closed transition pathway revealed by this study advances the understanding of the Bcl-2 transition from antiapoptotic to proapoptotic function, and yielded new structural insights for novel drug design against the BH4 domain. Communicated by Ramaswamy H. Sarma.

Cryptococcus neoformans sexual reproduction is controlled by a quorum sensing peptide.

Bacterial quorum sensing is a well-characterized communication system that governs a large variety of collective behaviours. By comparison, quorum sensing regulation in eukaryotic microbes remains poorly understood, especially its functional role in eukaryote-specific behaviours, such as sexual reproduction. Cryptococcus neoformans is a prevalent fungal pathogen that has two defined sexual cycles (bisexual and unisexual) and is a model organism for studying sexual reproduction in fungi. Here, we show that the quorum sensing peptide Qsp1 serves as an important signalling molecule for both forms of sexual reproduction. Qsp1 orchestrates various differentiation and molecular processes, including meiosis, the hallmark of sexual reproduction. It activates bisexual mating, at least in part through the control of pheromone, a signal necessary for bisexual activation. Notably, Qsp1 also plays a major role in the intercellular regulation of unisexual initiation and coordination, in which pheromone is not strictly required. Through a multi-layered genetic screening approach, we identified the atypical zinc finger regulator Cqs2 as an important component of the Qsp1 signalling cascade during both bisexual and unisexual reproduction. The absence of Cqs2 eliminates the Qsp1-stimulated mating response. Together, these findings extend the range of behaviours governed by quorum sensing to sexual development and meiosis.

Computational redesign of enzymes for regio- and enantioselective hydroamination.

Introduction of innovative biocatalytic processes offers great promise for applications in green chemistry. However, owing to limited catalytic performance, the enzymes harvested from nature's biodiversity often need to be improved for their desired functions by time-consuming iterative rounds of laboratory evolution. Here we describe the use of structure-based computational enzyme design to convert Bacillus sp. YM55-1 aspartase, an enzyme with a very narrow substrate scope, to a set of complementary hydroamination biocatalysts. The redesigned enzymes catalyze asymmetric addition of ammonia to substituted acrylates, affording enantiopure aliphatic, polar and aromatic ß-amino acids that are valuable building blocks for the synthesis of pharmaceuticals and bioactive compounds. Without a requirement for further optimization by laboratory evolution, the redesigned enzymes exhibit substrate tolerance up to a concentration of 300 g/L, conversion up to 99%, ß-regioselectivity >99% and product enantiomeric excess >99%. The results highlight the use of computational design to rapidly adapt an enzyme to industrially viable reactions.

Engineering improved thermostability of the GH11 xylanase from Neocallimastix patriciarum via computational library design.

Xylanases, which cleave the ß-1,4-glycosidic bond between xylose residues to release xylooligosaccharides (XOS), are widely used as food additives, animal feeds, and pulp bleaching agents. However, the thermally unstable nature of xylanases would hamper their industrial application. In this study, we used in silico design in a glycoside hydrolase family (GH) 11 xylanase to stabilize the enzyme. A combination of the best mutations increased the apparent melting temperature by 14 °C and significantly enhanced thermostability and thermoactivation. The variant also showed an upward-shifted optimal temperature for catalysis without compromising its activity at low temperatures. Moreover, a 10-fold higher XOS production yield was obtained at 70 °C, which compensated the low yield obtained with the wild-type enzyme. Collectively, the variant constructed by the computational strategy can be used as an efficient biocatalyst for XOS production at industrially viable conditions.

Molecular dynamic modeling of CYP51B in complex with azole inhibitors.

Cytochrome P450 14α-sterol demethylase (CYP51), the key enzyme in sterol biosynthesis pathway, is an important target protein of cholesterol-lowering agents, antifungal drugs, and herbicides. CYP51B enzyme is one of the CYP51 family members. In the present study, we have chosen four representative inhibitors of CYP51B: diniconazole (Din), fluconazole (Flu), tebuconazole (Teb), and voriconazole (Vor), and launched to investigate the binding features of CYP51B-inhibitor and gating characteristics via molecular docking and molecular dynamics (MD) simulations. The results show that the ranking of binding affinities among these inhibitors is Din > Teb > Vor > Flu. Din shows the strongest binding affinity, whereas Flu shows the weakest binding affinity. More importantly, based on the calculated binding free energy results, we hypothesize that the nonpolar interactions are the most important contributors, and three key residues (Thr77, Ala258, and Lys454) play crucial role in protein-inhibitor binding. Besides, residue Phe180 may play a molecular switch role in the process of the CYP51B-Teb and CYP51B-Vor binding. Additionally, Tunnel analysis results show that the major tunnel of Din, Flu, and Teb is located between helix K, FG loop, and ß4 hairpin (Tunnel II).The top ranked possible tunnel (Tunnel II) corresponds to Vor exits through helix K, F and helix J. This study further revealed the CYP51B relevant structural characteristics at the atomic level as well as provided a basis for rational design of new and more efficacious antifungal agents.

Molecular dynamics investigations of membrane-bound CYP2C19 polymorphisms reveal distinct mechanisms for peripheral variants by long-range effects on the enzymatic activity.

Increasing sophistication in methods used to account for human polymorphisms in susceptibility to drug metabolism has been one of the success stories in the prevention of adverse drug reactions. Genetic polymorphisms in drug-metabolizing enzymes can affect enzyme activity and cause differences in treatment response or drug toxicity. CYP2C19 is one of the most highly polymorphic CYP enzymes and acts on 10-15% of drugs in current clinical use. Despite the number of experimental analyses carried out for this system, the detailed structural basis for altered catalytic properties of polymorphic CYP2C19 variants remains unexplained at the atomic level. To this end, we have investigated the mutation effects on structural characteristics and tunnel geometry upon single point mutations to elucidate the underlying molecular mechanism for the enzymatic activity deficiencies by using the fully atomistic molecular dynamics (MD) simulations in their native, membrane-bound cellular environment. The obtained results demonstrate how significant sequence divergence causes heterogeneous variations, and further affects the shape and chemical properties of the substrate binding site. Principal component analysis (PCA) results combined with free energy calculations have revealed distinct mechanisms for different peripheral variants, implying a more complicated process for the decrease/loss of enzymatic activity upon the introduction of point mutations in CYP2C19 rather than simply structural changes of the region where the mutation is located. Overall, our present study provides important insights into the current pharmacogenetic knowledge of human drug-metabolizing CYP2C19 to understand the large inter-individual variability in drug clearance. The knowledge of heterogeneous variations in structural features could guide future experimental and computational work on efficient and safe drug treatment with better pharmacokinetic properties based on the common variant alleles of CYP genes, which varies among different ethnic populations.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms.

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.

Molecular dynamics investigations of BioH protein substrate specificity for biotin synthesis.

BioH, an enzyme of biotin synthesis, plays an important role in fatty acid synthesis which assembles the pimelate moiety. Pimeloyl-acyl carrier protein (ACP) methyl ester, which is long known to be a biotin precursor, is the physiological substrate of BioH. Azelayl methyl ester, which has a longer chain than pimeloyl methyl ester, conjugated to ACP is also indeed accepted by BioH with very low rate of hydrolysis. To date, the substrate specificity for BioH and the molecular origin for the experimentally observed rate changes of hydrolysis by the chain elongation have remained elusive. To this end, we have investigated chain elongation effects on the structures by using the fully atomistic molecular dynamics simulations combined with binding free energy calculations. The results indicate that the substrate specificity is determined by BioH together with ACP. The added two methylenes would increase the structural flexibility by protein motions at the interface of ACP and BioH, instead of making steric clashes with the side chains of the BioH hydrophobic cavity. On the other hand, the slower hydrolysis of azelayl substrate is suggested to be associated with the loose of contacts between BioH and ACP, and with the lost electrostatic interactions of two ionic/hydrogen bonding networks at the interface of the two proteins. The present study provides important insights into the structure-function relationships of the complex of BioH with pimeloyl-ACP methyl ester, which could contribute to further understanding about the mechanism of the biotin synthetic pathway, including the catalytic role of BioH.