Enzyme Flexibility Research Articles

Molecular insights and inhibitory dynamics of flavonoids in targeting Pim-1 kinase for cancer therapy.

Pim-1 kinase, a serine/threonine kinase, is often overexpressed in various cancers, contributing to disease progression and poor prognosis. In this study, we explored the potential of flavonoids as inhibitors of Pim-1 kinase using a combination of molecular docking and steered molecular dynamics (SMD) simulations. Our docking studies revealed two main binding orientations for the flavonoid molecules. The SMD simulations showed that the binding mode with higher pulling forces was linked to stronger inhibitory activity, with a strong positive correlation (R 2 ≈ 0.92) between pulling forces and IC50 values. Quercetin stood out as the most potent inhibitor, showing a pulling force of about 820pN and an IC_(5) 0 of less than 6µM. Further dynamic simulations indicated that quercetin's hydroxyl groups at the C3, C-5 and C-7 positions formed stable hydrogen bonds with key residues GLU-121, Leu-44 and Val-126, respectively enhancing its binding stability and effectiveness. Our results emphasized the critical role of the hydroxyl group at the C-3 position, which plays a pivotal function in effectively anchoring these molecules in the active site of Pim-1 kinase. Principal component analysis (PCA) of Pim-1 kinase's conformational changes revealed that potent inhibitors like quercetin, galangin, and kaempferol significantly restricted the enzyme's flexibility, suggesting potential inhibitory effect. These findings provide insights into the structural interactions between flavonoids and Pim-1 kinase, offering a foundation for future experimental investigations. However, further studies, including in vitro and in vivo validation, are necessary to assess the pharmacological relevance and specificity of flavonoids in cancer therapy.

A curated list of targeted optimized promiscuous ketoreductases (TOP-K).

Enzymes are either specific or promiscuous catalysts in nature. The latter is portrayed by protein families like CYP450Es, Aldo-ketoreductases and short/medium-chain dehydrogenases which participate in detoxification or secondary metabolite production. However, enzymes are evolutionarily 'blind' to an ever-increasing synthetic substrate library. Industries and laboratories have circumvented this by high-throughput screening or site-specific engineering to synthesize the product of interest. However, this paradigm entails cost and time-intensive one-enzyme, one-substrate catalysis model. One of the superfamilies regularly used for chiral alcohol synthesis are short-chain dehydrogenases/reductases (SDRs). Our objective is to determine a superset of promiscuous SDRs that can catalyze multiple ketones. They are typically classified into shorter 'Classical' and longer 'Extended' type ketoreductases. However, current analysis of modelled SDRs reveals a length-independent conserved N-terminus Rossmann-fold and a variable substrate-binding C-terminus substrate-binding region for both categories. The latter is recognized to influence the enzyme's flexibility and substrate promiscuity and we hypothesize these properties are directly linked with each other. We tested this by catalyzing ketone intermediates with the essential and specific enzyme: FabG_E, as well as non-essential SDRs such as UcpA and IdnO. The experimental results confirmed this biochemical-biophysical association, making it an interesting filter for ascertaining promiscuous enzymes. Hence, we created a dataset of physicochemical properties derived from the protein sequences and employed machine learning algorithms to examine potential candidates. This resulted in 24 targeted optimized ketoreductases (TOP-K) from 81 014 members. The experimental validation of select TOP-Ks demonstrated the correlation between the C-terminal lid-loop structure, enzyme flexibility and turnover rate on pro-pharmaceutical substrates.

Rational design of a highly efficient catalytic system for the production of PAPS from ATP and its application in the synthesis of chondroitin sulfate

The compound 3'-phosphoadenosine-5'-phosphosulfate (PAPS) serves as a sulfate group donor in the production of valuable sulfated compounds. However, elevated costs and low conversion efficiency limit the industrial applicability of PAPS. Here, we designed and constructed an efficient and controllable catalytic system for the conversion of adenosine triphosphate (ATP) (disodium salt) into PAPS without inhibition from by-products. In vitro and in vivo testing in Escherichia coli identified adenosine-5'-phosphosulfate kinase from Penicillium chrysogenum (PcAPSK) as the rate-limiting enzyme. Based on analysis of the catalytic steps and molecular dynamics simulations, a mechanism-guided "ADP expulsion" strategy was developed to generate an improved PcAPSK variant (L7), with a specific activity of 48.94 U·mg-1 and 73.27-fold higher catalytic efficiency (kcat/Km) that of the wild-type enzyme. The improvement was attained chiefly by reducing the ADP-binding affinity of PcAPSK, as well as by changing the enzyme's flexibility and lid structure to a more open conformation. By introducing PcAPSK L7 in an in vivo catalytic system, 73.59 mM (37.32 g·L-1 ) PAPS was produced from 150 mM ATP in 18.5 h using a 3-L bioreactor, and achieved titer is the highest reported to date and corresponds to a 98.13% conversion rate. Then, the PAPS catalytic system was combined with the chondroitin 4-sulfotransferase using a one-pot method. Finally, chondroitin sulfate was transformed from chondroitin at a conversion rate of 98.75%. This strategy has great potential for scale biosynthesis of PAPS and chondroitin sulfate.

An Atomistic Understanding of Allosteric Inhibition of Glutamate Racemase: a Dampening of Native Activation Dynamics.

Glutamate racemases (GR) are members of the family of bacterial enzymes known as cofactor-independent racemases and epimerases and catalyze the stereoinversion of glutamate. D-amino acids are universally important for the proper construction of viable bacterial cell walls, and thus have been repeatedly validated as attractive targets for novel antimicrobial drug design. Significant aspects of the mechanism of this challenging stereoinversion remain unknown. The current study employs a combination of MD and QM/MM computational approaches to show that the GR from H. pylori must proceed via a pre-activation step, which is dependent on the enzyme's flexibility. This mechanism is starkly different from previously proposed mechanisms. These findings have immediate pharmaceutical relevance, as the H. pylori GR enzyme is a very attractive allosteric drug target. The results presented in this study offer a distinctly novel understanding of how AstraZeneca's lead series of inhibitors cripple the H. pylori GR's native motions, via prevention of this critical chemical pre-activation step. Our experimental studies, using SPR, fluorescence and NMR WaterLOGSY, show that H. pylori GR is not inhibited by the uncompetitive mechanism originally put forward by Lundqvist et al.. The current study supports a deep connection between native enzyme motions and chemical reactivity, which has strong relevance to the field of allosteric drug discovery.

Mechanism of Inhibition of Glycoside Hydrolases Investigated by Molecular Dynamics Simulations

Biofuels are a well-known alternative to the largely used fossil-derived fuels. One of the main challenges for the production of the second-generation biofuels, produced by enzymatic degradation of agricultural waste, is the inhibition of the glycoside hydrolases enzymes by non-suitable substrates or by their own product. Employing molecular dynamics (MD) simulations, we investigated how two enzymes of the glycoside hydrolases family are affected by a non-suitable substrate (Man5B) or by their own product (CEL7A). Combined with previous experimental results, our simulations reveal that the reduction of Man5B enzymatic efficiency in the presence of a poor substrate (gluco-oligosaccharides) is associated with a loss of the enzyme's flexibility, the latter being required to bind new substrate, while the presence of a more suitable substrate (manno-oligosaccharides) does not pose this problem (1). MD simulations of the glycoside hydrolase enzyme CEL7A interacting with cellulose substrate of different lengths in position −1 to −7 offers a detailed view of its enzymatic mechanism and also on how the enzyme dynamics is affected by substrate. Employing steered molecular dynamics (SMD) we were able to analyze the difference in free-energy when a substrate is moving through CEL7A's catalytic tunnel under different conditions. All together our results are indicating that the inhibition of the CEL7A by its product is directly related to the cellobiose (enzymatic product) that was just cleaved and does not leave the exit pocket because of the already high concentration of cellobiose in solution.

Implementation of Traditional Medicine into Structure-Based Approaches for Antiviral Lead Discovery

The search for new drug leads against viruses remains an area of active investigations. In the last decades the influenza neuraminidase (iNA) has been extensively exploited by computational approaches in the search for small molecule inhibitors of this key enzyme responsible for viral replication and spread. These efforts fueled the development of currently worldwide licensed iNA inhibitors. However, the increase of drug-resistant influenza virus strains demonstrates the urgent need for innovative and effective antiviral agents. It is anticipated that the recent findings of the enzyme's flexibility and thereby extended binding pocket will help to overcome these hindrances.1

Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches

Based on integrative computational hybrid approaches that combined statistical coupling analysis (SCA), molecular dynamics (MD), and normal mode analysis (NMA), evolutionarily coupled residues involved in functionally relevant motion in the adenylate kinase protein family were identified. The hybrids identified four top-ranking site pairs that belong to a conserved hydrogen bond network that is involved in the enzyme's flexibility. A second group of top-ranking site pairs was identified in critical regions for functional dynamics, such as those related to enzymatic turnover. The high consistency of the results obtained by SCA with NMA (SCA.NMA) and by SCA.MD hybrid analyses suggests that suitable replacement of the matrix of cross-correlation analysis of atomic fluctuations (derived by using NMA) with those based on MD contributes to the identification of such sites by means of a fast computational calculation. The analysis presented here strongly supports the hypothesis that evolutionary forces, such as coevolution at the sequence level, have promoted functional dynamic properties of the adenylate kinase protein family. Finally, these hybrid approaches can be used to identify, at the residue level, protein motion coordination patterns not previously observed, such as in hinge regions.