Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms

C ucumis melo compounds: A new avenue for ALR-2 inhibition in diabetes mellitus

Diabetes mellitus (DM) is a prominent contributor to morbidity and mortality in developed nations, primarily attributable to vascular complications such as atherothrombosis occurring in the coronary arteries. Aldose reductase (ALR2), the main enzyme in the polyol pathway, catalyzes the conversion of glucose to sorbitol, leading to a significant buildup of reactive oxygen species in different tissues. It is therefore a prime candidate for therapeutic targeting, and extensive study is currently underway to discover novel natural compounds that can inhibit it. Cucumis melo (C. melo) has a long history as a lipid-lowering ethanopharmaceutical plant. In this study, compounds derived from C. melo were computationally evaluated as possible lead candidates. Various computational filtering methods were employed to assess the drug-like properties and ADMET (absorption, distribution, metabolism, excretion, and toxicity) profiles of the compounds. The compounds were subsequently addressed to analysis of their interactions, molecular docking, and molecular dynamics simulation studies. When compared to the conventional therapeutic compounds, three compounds exhibited enhanced binding affinity and intra-molecular residue interactions, resulting in increased stability and specificity. Consequently, four potent inhibitors, namely PubChem CIDs 119205, 65373, 6184, and 332427, have been identified. These inhibitors exhibit promising potential as pharmacological targets for the advancement of novel ALR-2 inhibitors.

Novel Insights into the Catalytic Mechanism of Collagenolysis by Zn(II)-Dependent Matrix Metalloproteinase-1

Collagen hydrolysis, catalyzed by Zn(II)-dependent matrix metalloproteinases (MMPs), is a critical physiological process. Despite previous computational investigations into the catalytic mechanisms of MMP-mediated collagenolysis, a significant knowledge gap in understanding remains regarding the influence of conformational sampling and entropic contributions at physiological temperature on enzymatic collagenolysis. In our comprehensive multilevel computational study, employing quantum mechanics/molecular mechanics (QM/MM) metadynamics (MetD) simulations, we aimed to bridge this gap and provide valuable insights into the catalytic mechanism of MMP-1. Specifically, we compared the full enzyme-substrate complex in solution, clusters in solution, and gas-phase to elucidate insights into MMP-1-catalyzed collagenolysis. Our findings reveal significant differences in the catalytic mechanism when considering thermal effects and the dynamic evolution of the system, contrasting with conventional static potential energy surface QM/MM reaction path studies. Notably, we observed a significant stabilization of the critical tetrahedral intermediate, attributed to contributions from conformational flexibility and entropy. Moreover, we found that protonation of the scissile bond nitrogen occurs via proton transfer from a Zn(II)-coordinated hydroxide rather than from a solvent water molecule. Following C-N bond cleavage, the C-terminus remains coordinated to the catalytic Zn(II), while the N-terminus forms a hydrogen bond with a solvent water molecule. Subsequently, the release of the C-terminus is facilitated by the coordination of a water molecule. Our study underscores the pivotal role of protein conformational dynamics at physiological temperature in stabilizing the transition state of the rate-limiting step and key intermediates, compared to the corresponding reaction in solution. These fundamental insights into the mechanism of collagen degradation provide valuable guidance for the development of MMP-1-specific inhibitors.

Tetrahydrobiopterin from cyanide-degrading bacterium Bacillus pumilus strain SVD06 induces apoptosis in human lung adenocarcinoma cell (A549)

Tetrahydrobiopterin (BH4) is an essential biological cofactor and a derivative of pterin which is considered potent anticancer agents. In continuation of our previous study on the identification of BH4 from cyanide-degrading Bacillus pumilus, the present study focuses on evaluating the anticancer properties of BH4 on A549, a human lung adenocarcinoma. Anticancer activity analysis shows that BH4 inhibited A549 cell growth after 24 h of incubation with 0.02 mg/mL. In acridine orange/ethidium bromide staining, BH4-treated A549 cells showed apoptotic morphology. BH4 treatment caused cell cycle arrest at G0/G1 phase compared to control cells. BH4 augmented p53 expression in alveolar cancer cells by downregulating MDM2 levels. There was downregulation of casp-3 and upregulation of iNOS gene in BH4-treated A549 cells. Further, docking studies indicated that BH4 had significant interactions with the above proteins affirming the apoptosis mechanism. Thus, BH4 could be considered a potential anticancer drug.

Absorption wavelength (TD-DFT) and adsorption of metal chalcogen clusters with methyl nicotinate: Structural, electronic, IRI, SERS, pharmacological and antiviral studies (HIV and omicron)

The DFT B3LYP-LAND2DZ technique is used to examine interactions of Methyl nicotinate with copper selenide and zinc selenide clusters. The existence of reactive sites is determined using ESP maps and Fukui data. The energy variations between HOMO and LUMO are utilised to calculate various energy parameters. The Atoms in Molecules and ELF (Electron Localisation Function) maps are employed to investigate the topology of the molecule. The Interaction Region Indicator is used to determine the existence of non-covalent zones in the molecule. The UV-Vis spectrum using the TD-DFT method and DOS graphs are used to obtain the theoretical determination of electronic transition and properties. Structural analysis of the compound is obtained using theoretical IR spectra. To explore the adsorption of copper selenide and zinc selenide clusters on the Methyl nicotinate, the adsorption energy and theoretical SERS spectra are employed. Furthermore, pharmacological investigations are carried out to confirm the drug's non-toxicity. The compound's antiviral efficacy against HIV and Omicron is demonstrated via protein-ligand docking.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

Mechanism of reactivation of the peroxidase catalytic activity of human cyclooxygenases by reducing cosubstrate quercetin

Our earlier studies show that the peroxidase activity of cyclooxygenase 1 and 2 (COX-1 and COX-2) can be reactivated in vitro and in vivo by the presence of certain naturally-occurring flavonoids such as quercetin and myricetin, which serve as reducing cosubstrates. These compounds can activate COX at nanomolar concentrations. In the present study, quercetin is used as a representative model compound to investigate the chemical mechanism by which the peroxidase activity of human COX-1 and COX-2 is reactivated after each catalytic cycle. Molecular docking and quantum mechanics calculations are carried out to probe the interactions of quercetin with the peroxidase sites of COX-1/2 and the reactivation mechanism. It is found that some of the partially-ionized states of quercetin can bind tightly and closely inside the peroxidase active sites of the COX enzymes and directly interact with heme Fe ion. While quercetin contains several phenolic hydroxyl groups, it is found that only the C-3'-OH group can effectively donate an electron for the reduction of heme because it not only can bind closely and tightly inside the peroxidase sites of COX-1/2, but it can also facilely donate an electron to heme Fe ion. This investigation provides a mechanistic explanation for the chemical process by which quercetin reactivates COX-1/2 peroxidases. This knowledge would aid in the rational design of drugs that can selectively target the peroxidase sites of COX-1/2 either as activators or inhibitors.

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Design and SAR Analysis of Covalent Inhibitors Driven by Hybrid QM/MM Simulations

Quantum mechanics/molecular mechanics (QM/MM) hybrid technique is emerging as a reliable computational method to investigate and characterize chemical reactions occurring in enzymes. From a drug discovery perspective, a thorough understanding of enzyme catalysis appears pivotal to assist the design of inhibitors able to covalently bind one of the residues belonging to the enzyme catalytic machinery. Thanks to the current advances in computer power, and the availability of more efficient algorithms for QM-based simulations, the use of QM/MM methodology is becoming a viable option in the field of covalent inhibitor design. In the present review, we summarized our experience in the field of QM/MM simulations applied to drug design problems which involved the optimization of agents working on two well-known drug targets, namely fatty acid amide hydrolase (FAAH) and epidermal growth factor receptor (EGFR). In this context, QM/MM simulations gave valuable information in terms of geometry (i.e., of transition states and metastable intermediates) and reaction energetics that allowed to correctly predict inhibitor binding orientation and substituent effect on enzyme inhibition. What is more, enzyme reaction modelling with QM/MM provided insights that were translated into the synthesis of new covalent inhibitor featured by a unique combination of intrinsic reactivity, on-target activity, and selectivity.