Dynamic Modelling Reveals 'Hotspots' on the Pathway to Enzyme-Substrate Complex Formation

Transient binding sites at the surface of haloalkane dehalogenase LinB as locations for fine-tuning enzymatic activity

The haloalkane dehalogenase LinB is a well-known enzyme that contains buried active site and is used for many modelling studies. Using classical molecular dynamics simulations of enzymes and substrates, we searched for transient binding sites on the surface of the LinB protein by calculating maps of enzyme-ligand interactions that were then transformed into sparse matrices. All residues considered as functionally important for enzyme performance (e.g., tunnel entrances) were excluded from the analysis to concentrate rather on non-obvious surface residues. From a set of 130 surface residues, twenty-six were proposed as a promising improvement of enzyme performance. Eventually, based on rational selection and filtering out the potentially unstable mutants, a small library of ten mutants was proposed to validate the possibility of fine-tuning the LinB protein. Nearly half of the predicted mutant structures showed improved activity towards the selected substrates, which demonstrates that the proposed approach could be applied to identify non-obvious yet beneficial mutations for enzyme performance especially when obvious locations have already been explored.

Molecular Mechanism of Methanol Inhibition in CALB-Catalyzed Alcoholysis: Analyzing Molecular Dynamics Simulations by a Markov State Model

Lipases are widely used enzymes that catalyze hydrolysis and alcoholysis of fatty acid esters. At high concentrations of small alcohols such as methanol or ethanol, many lipases are inhibited by the substrate. The molecular basis of the inhibition of Candida antarctica lipase B (CALB) by methanol was investigated by unbiased molecular dynamics (MD) simulations, and the substrate binding kinetics was analyzed by Markov state models (MSMs). The modeled fluxes of productive methanol binding at concentrations between 50 mM and 5.5 M were in good agreement with the experimental activity profile of CALB, with a peak at 300 mM. The kinetic and structural analysis uncovered the molecular basis of CALB inhibition. Beyond 300 mM, the kinetic bottleneck results from crowding of methanol in the substrate access channel, which is caused by the gradual formation of methanol patches close to Leu140 (helix α5), Leu278, and Ile285 (helix α10) at a distance of 4-5 Å from the active site. Our findings demonstrate the usefulness of unbiased MD simulations to study enzyme-substrate interactions at realistic substrate concentrations and the feasibility of scale-bridging by an MSM analysis to derive kinetic information.

Fibrinolytic protease from Bacillus cereus S46: Purification, characterization, and evaluation of its in vitro thrombolytic potential

Intravascular thrombosis is a prime cause of cardiac complications worldwide. Microbial fibrinolytic proteases are of clinical significance in thrombosis treatment. The present study discusses the purification and characterization of a protease from Bacillus cereus S46, ascertaining its in vitro thrombolytic activity against a blood clot. By the three-step purification involving precipitation, dialysis, and diethylaminoethyl-cellulose ion-exchange chromatography, a 12.37-fold purification of the enzyme to homogeneity was achieved. The apparent molecular mass of the protease was 30 kDa, as found by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The optimum activity of the enzyme was observed at pH 8.0 and 40°C. The enzyme retained an 82.19% residual activity at pH 8.0 and 40°C for 1 h. The K_m and V_max values of the protease with casein were 0.0027 mM and 9.712 µmol/min, respectively. In an in vitro assay, the purified protease resulted in 97.02% lysis of the blood clot. The fibrinolytic potential of the enzyme, together with its characteristics of being active and stable under near-physiological conditions, may suggest its application as a therapeutic agent.

Substrate Locking Promotes Dimer-Dimer Docking of an Enzyme Antibiotic Target

Protein dynamics manifested through structural flexibility play a central role in the function of biological molecules. Here we explore the substrate-mediated change in protein flexibility of an antibiotic target enzyme, Clostridium botulinum dihydrodipicolinate synthase. We demonstrate that the substrate, pyruvate, stabilizes the more active dimer-of-dimers or tetrameric form. Surprisingly, there is little difference between the crystal structures of apo and substrate-bound enzyme, suggesting protein dynamics may be important. Neutron and small-angle X-ray scattering experiments were used to probe substrate-induced dynamics on the sub-second timescale, but no significant changes were observed. We therefore developed a simple technique, coined protein dynamics-mass spectrometry (ProD-MS), which enables measurement of time-dependent alkylation of cysteine residues. ProD-MS together with X-ray crystallography and analytical ultracentrifugation analyses indicates that pyruvate locks the conformation of the dimer that promotes docking to the more active tetrameric form, offering insight into ligand-mediated stabilization of multimeric enzymes.

Dihydrodipicolinate synthase is absent in fungi

The class I aldolase dihydrodipicolinate synthase (DHDPS) catalyzes the first committed step of the diaminopimelate (DAP) lysine biosynthesis pathway in bacteria, archaea and plants. Despite the existence, in databases, of numerous fungal sequences annotated as DHDPS, its presence in fungi has been the subject of contradictory claims. We report the characterization of DHDPS candidates from fungi. Firstly, the putative DHDPS from Coccidioides immitis (PDB ID: 3QFE) was shown to have negligible enzyme activity. Sequence analysis of 3QFE showed that three out of the seven amino acid residues critical for DHDPS activity are absent; however, exact matches to catalytic residues from two other class I aldolases, 2-keto-3-deoxygluconate aldolase (KDGA), and 4-hydroxy-2-oxoglutarate aldolase (HOGA), were identified. The presence of both KDGA and HOGA activity in 3QFE was confirmed in vitro using enzyme assays, the first report of such dual activity. Subsequent analyses of all publically available fungal sequences revealed that no entry contains all seven residues important for DHDPS function. The candidate with the highest number of identities (6 of 7), KIW77228 from Fonsecaea pedrosoi, was shown to have trace DHDPS activity in vitro, partially restored by substitution of the seventh critical residue, and to be incapable of complementing DHDPS-deficient E. coli cells. Combined with the presence of all seven sequences for the alternative α-aminoadipate (AAA) lysine biosynthesis pathway in C. immitis and F. pedrosoi, we believe that DHDPS and the DAP pathway are absent in fungi, and further, that robust informed methods for annotating genes need to be implemented.

Molecular evolution of an oligomeric biocatalyst functioning in lysine biosynthesis

Dihydrodipicolinate synthase (DHDPS) is critical to the production of lysine through the diaminopimelate (DAP) pathway. Elucidation of the function, regulation and structure of this key class I aldolase has been the focus of considerable study in recent years, given that the dapA gene encoding DHDPS has been found to be essential to bacteria and plants. Allosteric inhibition by lysine is observed for DHDPS from plants and some bacterial species, the latter requiring a histidine or glutamate at position 56 (Escherichia coli numbering) over a basic amino acid. Structurally, two DHDPS monomers form the active site, which binds pyruvate and (S)-aspartate β-semialdehyde, with most dimers further dimerising to form a tetrameric arrangement around a solvent-filled centre cavity. The architecture and behaviour of these dimer-of-dimers is explored in detail, including biophysical studies utilising analytical ultracentrifugation, small-angle X-ray scattering and macromolecular crystallography that show bacterial DHDPS tetramers adopt a head-to-head quaternary structure, compared to the back-to-back arrangement observed for plant DHDPS enzymes. Finally, the potential role of pyruvate in providing substrate-mediated stabilisation of DHDPS is considered.

Comparison of untagged and his-tagged dihydrodipicolinate synthase from the enteric pathogen Vibrio cholerae

Given the emergence of multi drug resistant Vibrio cholerae strains, there is an urgent need to characterize new anti-cholera targets. One such target is the enzyme dihydrodipicolinate synthase (DHDPS; EC 4.3.3.7), which catalyzes the first committed step in the diaminopimelate pathway. This pathway is responsible for the production of two key metabolites in bacteria and plants, namely meso-2,6-diaminopimelate and L-lysine. Here, we report the cloning, expression and purification of untagged and His-tagged recombinant DHDPS from V. cholerae (Vc-DHDPS) and provide comparative structural and kinetic analyses. Structural studies employing circular dichroism spectroscopy and analytical ultracentrifugation demonstrate that the recombinant enzymes are folded and exist as dimers in solution. Kinetic analyses of untagged and His-tagged Vc-DHDPS show that the enzymes are functional with specific activities of 75.6 U/mg and 112 U/mg, K_M (pyruvate) of 0.14 mM and 0.15 mM, K_M (L-aspartate-4-semialdehyde) of 0.08 mM and 0.09 mM, and k_cat of 34 and 46 s^-1, respectively. These results demonstrate there are no significant changes in the structure and function of Vc-DHDPS upon the addition of an N-terminal His tag and, hence, the tagged recombinant product is suitable for future studies, including screening for new inhibitors as potential anti-cholera agents. Additionally, a polyclonal antibody raised against untagged Vc-DHDPS is validated for specifically detecting recombinant and native forms of the enzyme.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.