Activity prediction of substrates in NADH-dependent carbonyl reductase by docking requires catalytic constraints and charge parameterization of catalytic zinc environment

Advanced Insights into Catalytic and Structural Features of the Zinc-Dependent Alcohol Dehydrogenase from Thauera aromatica

The asymmetric reduction of ketones to chiral hydroxyl compounds by alcohol dehydrogenases (ADHs) is an established strategy for the provision of valuable precursors for fine chemicals and pharmaceutics. However, most ADHs favor linear aliphatic and aromatic carbonyl compounds, and suitable biocatalysts with preference for cyclic ketones and diketones are still scarce. Among the few candidates, the alcohol dehydrogenase from Thauera aromatica (ThaADH) stands out with a high activity for the reduction of the cyclic α‐diketone 1,2‐cyclohexanedione to the corresponding α‐hydroxy ketone. This study elucidates catalytic and structural features of the enzyme. ThaADH showed a remarkable thermal and pH stability as well as stability in the presence of polar solvents. A thorough description of the substrate scope combined with the resolution and description of the crystal structure, demonstrated a strong preference of ThaADH for cyclic α‐substituted cyclohexanones, and indicated structural determinants responsible for the unique substrate acceptance.

Structural Studies of β-Diketones and Their Implications on Biological Effects

The paper briefly summarizes methods to determine the structure of β-diketones with emphasis on NMR methods. Density functional calculations are also briefly treated. Emphasis is on the tautomeric equilibria of β-diketones in relation to biological effects. Relevant physical parameters such as acidity and solubility are treated. A series of biologically active molecules are treated with respect to structure (tautomerism). Characteristic molecules or groups of molecules are usnic acids, tetramic and tetronic acids, o-hydroxydibenzoylmethanes, curcumines, lupulones, and hyperforines.

Reversal of Regioselectivity in Zinc-Dependent Medium-Chain Alcohol Dehydrogenase from Rhodococcus erythropolis toward Octanone Derivatives

The zinc-dependent medium-chain alcohol dehydrogenase from Rhodococcus erythropolis (ReADH) is one of the most versatile biocatalysts for the stereoselective reduction of ketones to chiral alcohols. Despite its known broad substrate scope, ReADH only accepts carbonyl substrates with either a methyl or an ethyl group adjacent to the carbonyl moiety; this limits its use in the synthesis of the chiral alcohols that serve as a building blocks for pharmaceuticals. Protein engineering to expand the substrate scope of ReADH toward bulky substitutions next to carbonyl group (ethyl 2-oxo-4-phenylbutyrate) opens up new routes in the synthesis of ethyl-2-hydroxy-4-phenylbutanoate, an important intermediate for anti-hypertension drugs like enalaprilat and lisinopril. We have performed computer-aided engineering of ReADH toward ethyl 2-oxo-4-phenylbutyrate and octanone derivatives. W296, which is located in the small binding pocket of ReADH, sterically restricts the access of ethyl 2-oxo-4-phenylbutyrate, octan-3-one or octan-4-one toward the catalytic zinc ion and thereby limits ReADH activity. Computational analysis was used to identify position W296 and site-saturation mutagenesis (SSM) yielded an improved variant W296A with a 3.6-fold improved activity toward ethyl 2-oxo-4-phenylbutyrate when compared to WT ReADH (ReADH W296A: 17.10 U/mg and ReADH WT: 4.7 U/mg). In addition, the regioselectivity of ReADH W296A is shifted toward octanone substrates. ReADH W296A has a more than 16-fold increased activity toward octan-4-one (ReADH W296A: 0.97 U/mg and ReADH WT: 0.06 U/mg) and a more than 30-fold decreased activity toward octan-2-one (ReADH W296A: 0.23 U/mg and ReADH WT: 7.69 U/mg). Computational and experimental results revealed the role of position W296 in controlling the substrate scope and regiopreference of ReADH for a variety of carbonyl substrates.

Insight into the Highly Conserved and Differentiated Cofactor-Binding Sites of meso-Diaminopimelate Dehydrogenase StDAPDH

meso-Diaminopimelate dehydrogenase ( meso-DAPDH) is a good candidate for one-step synthesis of d-amino acid from 2-keto acids. Our previous research revealed the classification of meso-DAPDH family and showed that type II meso-DAPDH, such as the meso-DAPDH from Symbiobacterium thermophilum (StDAPDH), could catalyze reductive amination. In this article, seven residues of StDAPDH, which are highly conserved in each subfamily but are different between two subfamilies, were targeted to explore the relationships between structure and function. Determination of kinetic parameters showed that the amino acid residues, including P69, K159, V68, S90, V14, and V156, played very important roles in the catalytic function of StDAPDH. Molecular dynamics simulation revealed that these point mutations reduced the productive conformations by the newly formed or eliminated interactions between the residues and ligands. These results strengthen our understanding of the catalytic mechanism and evolution of meso-DAPDH and can aid future endeavors in enzyme engineering.

Computational insight into the phthalocyanine-DNA binding via docking and molecular dynamics simulations

Phthalocyanines are considered as good DNA binders, which makes them promising anti-tumor drug leads. The purpose of this study is to investigate the interactions between DNA and quaternary metallophthalocyanine derivatives (Q-MPc) possessing varying metals (M = Zn, Ni, Cu, Fe, Mg and Ca) by molecular docking since there seems to be a lack of information in the literature regarding this issue. In this direction, Autodock Vina and Molegro Virtual Docker programs were employed. Autodock Vina results reveal that each Q-MPc derivative binds to DNA strongly with similar binding energies and almost identical binding modes. They bind to the grooves of DNA by constituting favorable interactions between phosphate groups of DNA and Q-MPcs. Although changing the metal has no significant effect on binding, presence of quaternary amine substituents increases the binding constant K_b by about 2-fold comparing to the core Pc (ZnPc). Contrary to Autodock Vina, the calculated Molegro Virtual Docker binding scores have been more diverse indicating that the scoring function of Molegro is better in differentiating these metals. Despite the fact that Molegro is superior to Autodock Vina in terms of metal characterization, Autodock Vina and Molegro exhibit similar binding sites for the studied metallophthalocyanines. We propose that Q-MPc derivatives designed in this study are promising anti-tumor lead compounds since they tightly bind to DNA with considerably high K_b values. Cationic substituents and presence of metal have both positive effects on DNA binding which is critical for designing DNA-active drugs. Additional calculations employing molecular dynamics (MD) simulations verified the stability of Q-MPc-DNA complexes which remained in contact after 20 ns via attractive interactions mainly between DNA backbone and the Pc metal center.

Catalytic, Computational, and Evolutionary Analysis of the d-Lactate Dehydrogenases Responsible for d-Lactic Acid Production in Lactic Acid Bacteria

d-Lactate dehydrogenase (d-LDH) catalyzes the reversible reaction pyruvate + NADH + H⁺ ↔ lactate + NAD⁺, which is a principal step in the production of d-lactate in lactic acid bacteria. In this study, we identified and characterized the major d-LDH (d-LDH1) from three d-LDHs in Leuconostoc mesenteroides, which has been extensively used in food processing. A molecular simulation study of d-LDH1 showed that the conformation changes during substrate binding. During catalysis, Tyr101 and Arg235 bind the substrates by hydrogen bonds and His296 acts as a general acid/base for proton transfer. These residues are also highly conserved and have coevolved. Point mutations proved that the substrate binding sites and catalytic site are crucial for enzyme activity. Network and phylogenetic analyses indicated that d-LDH1 and the homologues are widely distributed but are most abundant in bacteria and fungi. This study expands the understanding of the functions, catalytic mechanism, and evolution of d-LDH.

Rational design of Meso-2,3-butanediol dehydrogenase by molecular dynamics simulation and experimental evaluations

Meso-2,3-butanediol dehydrogenase (meso-2,3-BDH) catalyzes NAD⁺ -dependent conversion of meso-2,3-butanediol to acetoin, a crucial external energy storage molecule in fermentive bacteria. In this study, the active tunnel of meso-2,3-BDH was identified. The two short α helixes positioned away from the α4-helix possibly expose the hydrophobic ligand-binding cavity, gating the exit of product and cofactor from the activity pocket. Further MM/GBSA-binding free energy analysis shows that Phe212 and Asn146 function as the key product-release sites. Site-directed mutagenesis experiments targeted to the sites show that the k_cat of Phe212Tyr is enhanced up to (4-8)-fold. The original activity of Asn146Gln is retained, but the activity of Asn146Ala mutation is lost. These results could provide helpful guidance on rational design of short-chain dehydrogenases/reductases.

What's My Substrate? Computational Function Assignment of Candida parapsilosis ADH5 by Genome Database Search, Virtual Screening, and QM/MM Calculations

Zinc-dependent medium chain reductase from Candida parapsilosis can be used in the reduction of carbonyl compounds to pharmacologically important chiral secondary alcohols. To date, the nomenclature of cpADH5 is differing (CPCR2/RCR/SADH) in the literature, and its natural substrate is not known. In this study, we utilized a substrate docking based virtual screening method combined with KEGG, MetaCyc pathway, and Candida genome databases search for the discovery of natural substrates of cpADH5. The virtual screening of 7834 carbonyl compounds from the ZINC database provided 94 aldehydes or methyl/ethyl ketones as putative carbonyl substrates. Out of which, 52 carbonyl substrates of cpADH5 with catalytically active docking pose were identified by employing mechanism based substrate docking protocol. Comparison of the virtual screening results with KEGG, MetaCyc database search, and Candida genome pathway analysis suggest that cpADH5 might be involved in the Ehrlich pathway (reduction of fusel aldehydes in leucine, isoleucine, and valine degradation). Our QM/MM calculations and experimental activity measurements affirmed that butyraldehyde substrates are the potential natural substrates of cpADH5, suggesting a carbonyl reductase role for this enzyme in butyraldehyde reduction in aliphatic amino acid degradation pathways. Phylogenetic tree analysis of known ADHs from Candida albicans shows that cpADH5 is close to caADH5. We therefore propose, according to the experimental substrate identification and sequence similarity, the common name butyraldehyde dehydrogenase cpADH5 for Candida parapsilosis CPCR2/RCR/SADH.