Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease

The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.     

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.     

Identification of essential acidic residues of outer membrane protease OmpT supports a novel active site

Escherichia coli outer membrane protease OmpT has previously been classified as a serine protease with Ser(99) and His(212) as active site residues. The recently solved X-ray structure of the enzyme was inconsistent with this classification, and the involvement of a nucleophilic water molecule was proposed. Here, we substituted all conserved aspartate and glutamate residues by alanines and measured the residual enzymatic activities of the variants. Our results support the involvement of a nucleophilic water molecule that is activated by the Asp(210)/His(212) catalytic dyad. Activity is also strongly dependent on Asp(83) and Asp(85). Both may function in binding of the water molecule and/or oxyanion stabilization. The proposed mechanism implies a novel proteolytic catalytic site.     

Catalytic sites of medicinal leech enzyme destabilase-lysozyme (Mldl). Structure-functional correlation

Based on three-dimensional model of the bifunctional enzyme Destabilase-Lysozyme (mlDL-Ds2) in complex with trimer of N-acetylglucosoamine (NAG)3 the functional role of the stereochemically based group of amino acids (Glu14, Asp26, Ser 29, Ser31, Lys38, His92), in manifestation of glycosidase and isopeptidase activities has been elucidated. By method of site-directed mutagenesis it has been shown that mlDL glycosidase active site includes catalytic Glu14 and Asp26, and isopeptidase site functions as Ser/Lys dyad presented by catalytic residues Lys38 and Ser29. Thus, among the invertebrate lysozymes mlDL presents first example of the bifunctional enzyme with identified position of the isopeptidase active site and localization of the corresponding catalytic residues.     

pH‐dependent conformational dynamics of beta‐secretase 1: A molecular dynamics study

Beta‐secretase 1 (BACE‐1) is an aspartyl protease implicated in the overproduction of β‐amyloid fibrils responsible for Alzheimer disease. The process of β‐amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH‐dependent dynamics of BACE‐1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β‐amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double‐protonated, mono‐protonated, and double‐deprotonated states over the known pH range of higher activity. These differences include Tyr 71‐aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE‐1 activity.     

First principles computational study of the active site of arginase

Ab initio density functional theory (DFT) methods were used to investigate the structural features of the active site of the binuclear enzyme rat liver arginase. Special emphasis was placed on the crucial role of the second shell ligand interactions. These interactions were systematically studied by performing calculations on models of varying size. It was determined that a water molecule, and not hydroxide, is the bridging exogenous ligand. The carboxylate ligands facilitate the close approach of the Mn (II) ions by attenuating the metal-metal electrostatic repulsion. Of the two metals, Mn(A) was shown to carry a larger positive charge. Analysis of the electronic properties of the active site revealed that orbitals involving the terminal Asp234 residue, as well as the flexible micro-1,1 bridging Asp232, lie at high energies, suggesting weaker coordination. This is reflected in certain structural variability present in our models and is also consistent with recent experimental findings. Finally, implications of our findings for the biological function of the enzyme are delineated.     

Catalytic sites of medicinal leech enzyme destabilase-lysozyme (mlDL): Structure-function relationship

Based on the three-dimensional model of the bifunctional enzyme destabilase-lysozyme of the medicinal leech (mlDL) in complex with trimer of N-acetylglucosamine (NAG)₃ by site-directed mutagenesis method, the functional role of the group of amino acids (Glu14, Asp26, Ser29, Ser31, Lys38, His92) in manifestation of lysozyme (glycosidase, muramidase) and isopeptidase activities has been investigated by site-directed mutagenesis. The results obtained go well with hypothesis, that lysozyme active site of mlDL includes catalytic Glu14 and Asp26 residues, and isopeptidase site functions as Ser/Lys catalytic dyad presented by catalytic residues Ser29 and Lys38. Thus, among the invertebrate lysozymes, mlDL presents the first example of a bifunctional enzyme with identified position of the isopeptidase active site and localization of the corresponding catalytic residues.     

Insights from atomic-resolution X-ray structures of chemically synthesized HIV-1 protease in complex with inhibitors

The human immunodeficiency virus 1 (HIV-1) protease (PR) is an aspartyl protease essential for HIV-1 viral infectivity. HIV-1 PR has one catalytic site formed by the homodimeric enzyme. We chemically synthesized fully active HIV-1 PR using modern ligation methods. When complexed with the classic substrate-derived inhibitors JG-365 and MVT-101, the synthetic HIV-1 PR formed crystals that diffracted to 1.04- and 1.2-A resolution, respectively. These atomic-resolution structures revealed additional structural details of the HIV-1 PR's interactions with its active site ligands. Heptapeptide inhibitor JG-365, which has a hydroxyethylamine moiety in place of the scissile bond, binds in two equivalent antiparallel orientations within the catalytic groove, whereas the reduced isostere hexapeptide MVT-101 binds in a single orientation. When JG-365 was converted into the natural peptide substrate for molecular dynamic simulations, we found putative catalytically competent reactant states for both lytic water and direct nucleophilic attack mechanisms. Moreover, free energy perturbation calculations indicated that the insertion of catalytic water into the catalytic site is an energetically favorable process.     

A molecular dynamics study of the structural stability of HIV-1 protease under physiological conditions: the role of Na+ ions in stabilizing the active site

HIV-1 protease is most active under weakly acidic conditions (pH 3.5-6.5), when the catalytic Asp25 and Asp25' residues share 1 proton. At neutral pH, this proton is lost and the stability of the structure is reduced. Here we present an investigation of the effect of pH on the dynamics of HIV-1 protease using MD simulation techniques. MD simulations of the solvated HIV-1 protease with the Asp25/25' residues monoprotonated and deprotonated have been performed. In addition we investigated the effect of the inclusion of Na(+) and Cl(-) ions to mimic physiological salt conditions. The simulations of the monoprotonated form and deprotonated form including Na(+) show very similar behavior. In both cases the protein remained stable in the compact, "self-blocked" conformation in which the active site is blocked by the tips of the flaps. In the deprotonated system a Na(+) ion binds tightly to the catalytic dyad shielding the repulsion between the COO(-) groups. Ab initio calculations also suggest the geometry of the active site with the Na(+) bound closely resembles that of the monoprotonated case. In the simulations of the deprotonated form (without Na(+) ions), a water molecule bound between the Asp25 Asp25' side-chains. This disrupted the dimerization interface and eventually led to a fully open conformation.     

Energetic and mechanistic studies of glucoamylase using molecular recognition of maltose OH groups coupled with site-directed mutagenesis

Molecular recognition using a series of deoxygenated maltose analogues was used to determine the substrate transition-state binding energy profiles of 10 single-residue mutants at the active site of glucoamylase from Aspergillus niger. The individual contribution of each substrate hydroxyl group to transition-state stabilization with the wild type and each mutant GA was determined from the relation Delta(DeltaG()) = -RT ln[(k(cat)/K(M))(x)/(k(cat)/K(M))(y)], where x represents either a mutant enzyme or substrate analogue and y the wild-type enzyme or parent substrate. The resulting binding energy profiles indicate that disrupting an active site hydrogen bond between enzyme and substrate, as identified in crystal structures, not only sharply reduces or eliminates the energy contributed from that particular hydrogen bond but also perturbs binding contributions from other substrate hydroxyl groups. Replacing the active site acidic groups, Asp55, Glu180, or Asp309, with the corresponding amides, and the neutral Trp178 with the basic Arg, all substantially reduced the binding energy contribution of the 4'- and 6'-OH groups of maltose at subsite -1, even though both Glu180 and Asp309 are localized at subsite 1. In contrast, the substitution, Asp176 --> Asn, located near subsites -1 and 1, did not substantially perturb any of the individual hydroxyl group binding energies. Similarly, the substitutions Tyr116 --> Ala, Ser119 --> Tyr, or Trp120 --> Phe also did not substantially alter the energy profiles even though Trp120 has a critical role in directing conformational changes necessary for activity. Since the mutations at Trp120 and Asp176 reduced k(cat) values by 50- and 12-fold, respectively, a large effect on k(cat) is not necessarily accompanied by changes in hydroxyl group binding energy contributions. Two substitutions, Asn182 --> Ala and Tyr306 --> Phe, had significant though small effects on interactions with 3- and 4'-OH, respectively. Binding interactions between the enzyme and the glucosyl group in subsite -1, particularly with the 4'- and 6'-OH groups, play an important role in substrate binding, while subsite 1 interactions may play a more important role in product release.     

Exploring the drug resistance mechanism of active site,non-active site mutations and their cooperative effects in CRF01_AE HIV-1 protease: molecular dynamics simulations and free energy calculations

Human immunodeficiency virus type 1 protease is essential for virus replication and maturation and has been considered as one of the important drug target for the antiretroviral treatment of HIV infection. The majority of HIV infections are caused due to non-B subtypes in developing countries. Subtype AE is spreading rapidly and infecting huge population worldwide. Understanding the interdependence of active and non-active site mutations in conferring drug resistance is crucial for the development effective inhibitors in subtype AE protease. In this work, we have investigated the mechanism of resistance against indinavir (IDV) due to therapy selected active site mutation V82F, non-active site mutations PF82V and their cooperative effects PV82F in subtype AE-protease using molecular dynamics simulations and binding free energy calculations. The simulations suggested all the three complexes lead to decrease in binding affinity of IDV, whereas the PF82V complex resulted in an enhanced binding affinity compared to V82F and PV82F complexes. Large positional deviation of IDV was observed in V82F complex. The preservation of hydrogen bonds of IDV with active site Asp25/Asp25′ and flap residue Ile50/50′ via a water molecule is crucial for effective binding. Owing to the close contact of 80s loop with Ile50′ and Asp25, the alteration between residues Thr80 and Val82, further induces conformational change thereby resulting in loss of interactions between IDV and the residues in the active site cavity, leading to drug resistance. Our present study shed light on the effect of active, non-active site mutations and their cooperative effects in AE protease.

Communicated by Ramaswamy H. Sarma     

Suppression of the back proton-transfer from Asp85 to the retinal Schiff base in bacteriorhodopsin: a theoretical analysis of structural elements

The transfer of a proton from the retinal Schiff base to the nearby Asp85 protein group is an essential step in the directional proton-pumping by bacteriorhodopsin. To avoid the wasteful back reprotonation of the Schiff base from Asp85, the protein must ensure that, following Schiff base deprotonation, the energy barrier for back proton-transfer from Asp85 to the Schiff base is larger than that for proton-transfer from the Schiff base to Asp85. Here, three structural elements that may contribute to suppressing the back proton-transfer from Asp85 to the Schiff base are investigated: (i) retinal twisting; (ii) hydrogen-bonding distances in the active site; and (iii) the number and location of internal water molecules. The impact of the pattern of bond twisting on the retinal deprotonation energy is dissected by performing an extensive set of quantum-mechanical calculations. Structural rearrangements in the active site, such as changes of the Thr89:Asp85 distance and relocation of water molecules hydrogen-bonding to the Asp85 acceptor group, may participate in the mechanism which ensures that following the transfer of the Schiff base proton to Asp85 the protein proceeds with the subsequent photocycle steps, and not with back proton transfer from Asp85 to the Schiff base.     

On the multiple functional roles of the active site histidine in catalysis and allosteric regulation of Escherichia coli glucosamine 6-phosphate deaminase

The active site of glucosamine-6-phosphate deaminase (EC 3.5.99.6, formerly 5.3.1.10) from Escherichia coli was first characterized on the basis of the crystallographic structure of the enzyme bound to the competitive inhibitor 2-amino-2-deoxy-glucitol 6-phosphate. The structure corresponds to the R allosteric state of the enzyme; it shows the side-chain of His143 in close proximity to the O5 atom of the inhibitor. This arrangement suggests that His143 could have a role in the catalysis of the ring-opening step of glucosamine 6-phosphate whose alpha-anomer is the true substrate. The imidazole group of this active-site histidine contacts the carboxy groups from Glu148 and Asp141, via its Ndelta1 atom [Oliva et al. (1995) Structure 3, 1323-1332]. These interactions change in the T state because the side chain of Glu148 moves toward the allosteric site, leaving at the active site the dyad Asp141-His143 [Horjales et al. (1999) Structure 7, 527-536]. In this research, a dual approach using site-directed mutagenesis and controlled chemical modification of histidine residues has been used to investigate the role of the active-site histidine. Our results support a multifunctional role of His143; in the forward reaction, it is involved in the catalysis of the ring-opening step of the substrate, glucosamine 6-P. In the reverse reaction, the substrate fructose 6-P binds in its open chain, carbonylic form. The role of His143 in the binding of both glucosamine 6-P and reaction intermediates in their extended-chain forms was demonstrated by binding experiments using the reaction intermediate analogue, 2-amino-2-deoxy-D-glucitol 6-phosphate. His143 was also shown to be a critical residue for the conformational coupling between active and allosteric sites. From the pH dependence of the reactivity of the active site histidine to diethyl dicarbonate, we observed a pK(a) change of 1.2 units to the acid side when the enzyme undergoes the allosteric T to R transition during which the side chain of Glu148 moves toward the active site. The kinetic study of the Glu148-Gln mutant deaminase shows that the loss of the carboxy group and its replacement with the corresponding amide modifies the k(cat) versus pH profile of the enzyme, suggesting that the catalytic step requiring the participation of His143 has become rate-limiting. This, in turn, indicates that the interaction Glu148-His143 in the wild-type enzyme in the R state contributes to make the enzyme functional over a wide pH range.     

An inhibitor binding pocket distinct from the catalytic active site on human beta-APP cleaving enzyme

Beta-APP cleaving enzyme (BACE) is responsible for the first of two proteolytic cleavages of the APP protein that together lead to the generation of the Alzheimer's disease-associated Abeta peptide. It is widely believed that halting the production of Abeta peptide, by inhibition of BACE, is an attractive therapeutic modality for the treatment of Alzheimer's disease. BACE is an aspartyl protease, and there is significant effort in the pharmaceutical community to apply traditional design methods to the development of active site-directed inhibitors of this enzyme. We report here the discovery of a ligand binding pocket within the catalytic domain of BACE that is distinct from the enzymatic active site (i.e., an exosite). Peptides, initially identified from combinatorial phage peptide libraries, contain the sequence YPYF(I/L)P(L/I) and bind specifically to this exosite, even in the presence of saturating concentrations of active site-directed inhibitors. Binding of peptides to the BACE exosite leads to a concentration-dependent inhibition of proteolysis for APP-related, protein-based substrates of BACE. The discovery of this exosite opens new opportunities for the identification and development of novel and potentially selective small molecule inhibitors of BACE that act through exosite, rather than active site, binding interactions.     

Protein engineering of uridine phosphorylase from Escherichia coli K-12. II. Comparative study of hybrid and mutant forms of uridine phosphorylases]

Genes for hybrid uridine phosphorylases (UPases) consisting of fragments of amino acid sequences of UPases from Escherichia coli and Salmonella typhimurium were constructed. Producing strains of the corresponding proteins were genetically engineered. Mutant forms of the E. coli K-12 UPase were produced by site-directed mutagenesis. A comparative study of the enzyme properties of the mutant and hybrid forms of bacterial UPases was performed. It was shown that Asp27 unlike Asp5 and Asp29 residues of the E. coli UPase forms part of the active site of the protein. A scheme of the involvement of Asp27 in the binding of inorganic phosphate is proposed.     

Calculation of substrate binding affinities for a bacterial GH78 rhamnosidase through molecular dynamics simulations

Ram2 from Pediococcus acidilactici is a rhamnosidase from the glycoside hydrolase family 78. It shows remarkable selectivity for rutinose rather than para-nitrophenyl-alpha-l-rhamnopyranoside (p-NPR). Molecular dynamics simulations were performed using a homology model of this enzyme, in complex with both substrates. Free energy calculations lead to predicted binding affinities of −34.4 and −30.6 kJ mol⁻¹ respectively, agreeing well with an experimentally estimated relative free energy of 5.4 kJ mol⁻¹. Further, the most relevant binding poses could be determined. While p-NPR preferably orients its rhamnose moiety toward the active site, rutinose interacts most strongly with its glucose moiety. A detailed hydrogen bond analysis confirms previously implicated residues in the active site (Asp217, Asp222, Trp226, Asp229 and Glu488) and quantifies the importance of individual residues for the binding. The most important amino acids are Asp229 and Phe339 which are involved in many interactions during the simulations. While Phe339 was observed in more simulations, Asp229 was involved in more persistent interactions (forming an average of at least 2 hydrogen bonds during the simulation). These analyses directly suggest mutations that could be used in a further experimental characterization of the enzyme. This study shows once more the strength of computer simulations to rationalize and guide experiments at an atomic level.     

Molecular docking based virtual screening of natural compounds as potential BACE1 inhibitors: 3D QSAR pharmacophore mapping and molecular dynamics analysis

Beta-site APP cleaving enzyme1 (BACE1) catalyzes the rate determining step in the generation of Aβ peptide and is widely considered as a potential therapeutic drug target for Alzheimer’s disease (AD). Active site of BACE1 contains catalytic aspartic (Asp) dyad and flap. Asp dyad cleaves the substrate amyloid precursor protein with the help of flap. Currently, there are no marketed drugs available against BACE1 and existing inhibitors are mostly pseudopeptide or synthetic derivatives. There is a need to search for a potent inhibitor with natural scaffold interacting with flap and Asp dyad. This study screens the natural database InterBioScreen, followed by three-dimensional (3D) QSAR pharmacophore modeling, mapping, in silico ADME/T predictions to find the potential BACE1 inhibitors. Further, molecular dynamics of selected inhibitors were performed to observe the dynamic structure of protein after ligand binding. All conformations and the residues of binding region were stable but the flap adopted a closed conformation after binding with the ligand. Bond oligosaccharide interacted with the flap as well as catalytic dyad via hydrogen bond throughout the simulation. This led to stabilize the flap in closed conformation and restricted the entry of substrate. Carbohydrates have been earlier used in the treatment of AD because of their low toxicity, high efficiency, good biocompatibility, and easy permeability through the blood–brain barrier. Our finding will be helpful in identify the potential leads to design novel BACE1 inhibitors for AD therapy.     

Nucleotide substrate recognition by UDP-N-acetylglucosamine acyltransferase (LpxA) in the first step of lipid A biosynthesis

Lipid A is an integral component of the lipopolysaccharide (LPS) that forms the selective and protective outer monolayer of Gram-negative bacteria, and is essential for bacterial growth and viability. UDP-N-acetylglucosamine acyltransferase (LpxA) initiates lipid A biosynthesis by catalyzing the transfer of R-3-hydroxymyristic acid from acyl carrier protein to the 3'-hydroxyl group of UDP-GlcNAc. The enzyme is a homotrimer, and previous studies suggested that the active site lies within a positively charged cleft formed at the subunit-subunit interface. The crystal structure of Escherichia coli LpxA in complex with UDP-GlcNAc reveals details of the substrate-binding site, with prominent hydrophilic interactions between highly conserved clusters of residues (Asn198, Glu200, Arg204 and Arg205) with UDP, and (Asp74, His125, His144 and Gln161) with the GlcNAc moiety. These interactions serve to bind and orient the substrate for catalysis. The crystallographic model supports previous results, which suggest that acylation occurs via nucleophilic attack of deprotonated UDP-GlcNAc on the acyl donor in a general base-catalyzed mechanism involving a catalytic dyad of His125 and Asp126. His125, the general base, interacts with the 3'-hydroxyl group of UDP-GlcNAc to generate the nucleophile. The Asp126 side-chain accepts a hydrogen bond from His125 and helps orient the general base to participate in catalysis. Comparisons with an LpxA:peptide inhibitor complex indicate that the peptide competes with both nucleotide and acyl carrier protein substrates.     

Conformation, stability, and active-site cysteine titrations of Escherichia coli D26A thioredoxin probed by Raman spectroscopy.

The active-site cysteines (Cys 32 and Cys 35) of Escherichia coli thioredoxin are oxidized to a disulfide bridge when the protein mediates substrate reduction. In reduced thioredoxin, Cys 32 and Cys 35 are characterized by abnormally low pKa values. A conserved side chain, Asp 26, which is sterically accessible to the active site, is also essential to oxidoreductase activity. pKa values governing cysteine thiol-thiolate equilibria in the mutant thioredoxin, D26A, have been determined by direct Raman spectrophotometric measurement of sulfhydryl ionizations. The results indicate that, in D26A thioredoxin, both sulfhydryls titrate with apparent pKa values of 7.5+/-0.2, close to values measured previously for wild-type thioredoxin. Sulfhydryl Raman markers of D26A and wild-type thioredoxin also exhibit similar band shapes, consistent with minimal differences in respective cysteine side-chain conformations and sulfhydryl interactions. The results imply that neither the Cys 32 nor Cys 35 SH donor is hydrogen bonded directly to Asp 26 in the wild-type protein. Additionally, the thioredoxin main-chain conformation is largely conserved with D26A mutation. Conversely, the mutation perturbs Raman bands diagnostic of tryptophan (Trp 28 and Trp 31) orientations and leads to differences in their pH dependencies, implying local conformational differences near the active site. We conclude that, although the carboxyl side chain of Asp 26 neither interacts directly with active-site cysteines nor is responsible for their abnormally low pKa values, the aspartate side chain may play a role in determining the conformation of the enzyme active site.     

Investigation into mechanism of orotidine 5′-monophosphate decarboxylase enzyme by MM-PBSA/MM-GBSA and molecular docking

In this work, we studied the binding affinity of orotidine 5′-monophosphate (OMP) and 6-hydroxy-UMP (BMP) for Saccharomyces cerevisiae orotidine 5′-monophosphate decarboxylase (OMPDC) enzyme by using Molecular Mechanics-Poisson–Boltzmann Surface Area (MM-PBSA) and the Molecular Mechanics-Generalised Born Surface Area (MM-GBSA) calculations. In all simulations, Asp91, which is an important residue in the enzyme active site, was considered in both anionic (present in the native form of the enzyme) and neutral states. A series of 10-ns molecular dynamics simulations were performed for the four OMPDC–ligand complexes, two ligand-free enzymes and two free ligands, followed by MM-PBSA/MM-GBSA calculations on the collected snapshots, and molecular docking calculations using the free enzymes and ligands. The results of MM-PBSA/MM-GBSA calculations indicate that all of the OMPDC–ligand complexes form favourable systems in water, which is in agreement with corresponding experimental data. The results of the MM-PBSA and molecular docking methods also showed that OMPDC–BMP complexes, transition state analogue and inhibitor of the OMPDC enzyme have the highest binding affinities. The fact that in the native anionic state BMP shows a higher binding affinity compared with the substrate suggests the contribution of a transition state stabilisation mechanism in the debatable catalytic mechanism of the OMPDC enzyme.