Browse
Submit Data
Databases
API
Help

Dataset Information

20 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Mechanism of inhibition of SARS-CoV-2 M^pro by N3 peptidyl Michael acceptor explained by QM/MM simulations and design of new derivatives with tunable chemical reactivity.

ABSTRACT: The SARS-CoV-2 main protease (M^pro) is essential for replication of the virus responsible for the COVID-19 pandemic, and one of the main targets for drug design. Here, we simulate the inhibition process of SARS-CoV-2 M^pro with a known Michael acceptor (peptidyl) inhibitor, N3. The free energy landscape for the mechanism of the formation of the covalent enzyme-inhibitor product is computed with QM/MM molecular dynamics methods. The simulations show a two-step mechanism, and give structures and calculated barriers in good agreement with experiment. Using these results and information from our previous investigation on the proteolysis reaction of SARS-CoV-2 M^pro, we design two new, synthetically accessible N3-analogues as potential inhibitors, in which the recognition and warhead motifs are modified. QM/MM modelling of the mechanism of inhibition of M^pro by these novel compounds indicates that both may be promising candidates as drug leads against COVID-19, one as an irreversible inhibitor and one as a potential reversible inhibitor.

SUBMITTER: Arafet K

PROVIDER: S-EPMC8179034 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Determining Michael Acceptor Reactivity from Kinetic, Mechanistic, and Computational Analysis for the Base-catalyzed Thiol-Michael Reaction.

Project description:A combined experimental and computational study of the reactivities of seven commonly used Michael acceptors paired with two thiols within the framework of photobase-catalyzed thiol-Michael reactions is reported. The rate coefficients of the propagation (kP), reverse propagation (k-P), chain-transfer (kCT), and overall reaction (koverall) were experimentally determined and compared with the well-accepted electrophilicity parameters of Mayr and Parr, and DFT-calculated energetics. Both Mayr's and Parr's electrophilicity parameters predict the reactivities of these structurally varying vinyl functional groups well, covering a range of overall reaction rate coefficients from 0.5 to 6.2 s-1. To gain insight into the individual steps, the relative energies have been calculated using DFT for each of the stationary points along this step-growth reaction between ethanethiol and the seven alkenes. The free energies of the individual steps reveal the underlying factors that control the reaction barriers for propagation and chain transfer. Both the propagation and chain transfer steps are under kinetic control. These results serve as a useful guide for Michael acceptor selection to design and predict thiol-Michael-based materials with appropriate kinetic and material properties.

| S-EPMC8409055 | biostudies-literature

An optimized RAD51 inhibitor that disrupts homologous recombination without requiring Michael acceptor reactivity.

Project description:Homologous recombination (HR) is an essential process in cells that provides repair of DNA double-strand breaks and lesions that block DNA replication. RAD51 is an evolutionarily conserved protein that is central to HR. Overexpression of RAD51 protein is common in cancer cells and represents a potential therapeutic target in oncology. We previously described a chemical inhibitor of RAD51, called RI-1 (referred to as compound 1 in this report). The chloromaleimide group of this compound is thought to act as a Michael acceptor and react with the thiol group on C319 of RAD51, using a conjugate addition-elimination mechanism. In order to reduce the likelihood of off-target effects and to improve compound stability in biological systems, we developed an analogue of compound 1 that lacks maleimide-based reactivity but retains RAD51 inhibitory activity. This compound, 1-(3,4-dichlorophenyl)-3-(4-methoxyphenyl)-4-morpholino-1H-pyrrole-2,5-dione, named RI-2 (referred to as compound 7a in this report), appears to bind reversibly to the same site on the RAD51 protein as does compound 1. Like compound 1, compound 7a specifically inhibits HR repair in human cells.

| S-EPMC3619390 | biostudies-other

Improvement of the Chemical Reactivity of Michael Acceptor of Ethacrynic Acid Correlates with Antiproliferative Activities.

Project description:The present study aims to report the design, synthesis, and biological activity of new ethacrynic acid (EA) analogs (6-10) obtained by the double modulation of the carboxylic acid moiety and the aromatic ring with the aim to increase the chemical reactivity of Michael acceptor of EA. All obtained compounds were characterized by 1H and 13C NMR, IR, and high-resolution mass spectrometry. The antiproliferative activity was evaluated in vitro using MMT test, in a first step, against HL60 cell line and in a second step, on a panel of human cancer cell lines such as HCT116, A549, MCF7, PC3, U87-MG, and SKOV3, and normal cell line MRC5 in comparison with positive control doxorubicin. Among all the tested compounds, the product 8 containing a propargyl and a hydroxyl groups, allowing an intramolecular hydrogen bond with the keto group of EA, exhibited a pronounced and selective activity in a nanomolar range against HL60, A549, PC3, and MCF7 with IC50 values of 15, 41.2, 68.7, and 61.5 nM, respectively. Compound 8 also showed a good selectivity index (SI) against HL60 and moderate SI against the other three human cancer cells (A549, PC3, and MCF7). The study of the structure-activity relationship showed that both modifications of the carboxylic group and the introduction of an intramolecular hydrogen bond are highly required to improve the antiproliferative activities. The molecular modeling studies of compound 8 revealed that it favorably binds to the glutathione S-transferase active site, which may explain its interesting anticancer activity. These new compounds have potential to be developed as novel therapeutic agents against various cancer types.

| S-EPMC9865193 | biostudies-literature

QM/MM calculations with deMon2k.

Project description:The density functional code deMon2k employs a fitted density throughout (Auxiliary Density Functional Theory), which offers a great speed advantage without sacrificing necessary accuracy. Powerful Quantum Mechanical/Molecular Mechanical (QM/MM) approaches are reviewed. Following an overview of the basic features of deMon2k that make it efficient while retaining accuracy, three QM/MM implementations are compared and contrasted. In the first, deMon2k is interfaced with the CHARMM MM code (CHARMM-deMon2k); in the second MM is coded directly within the deMon2k software; and in the third the Chemistry in Ruby (Cuby) wrapper is used to drive the calculations. Cuby is also used in the context of constrained-DFT/MM calculations. Each of these implementations is described briefly; pros and cons are discussed and a few recent applications are described briefly. Applications include solvated ions and biomolecules, polyglutamine peptides important in polyQ neurodegenerative diseases, copper monooxygenases and ultra-rapid electron transfer in cryptochromes.

| S-EPMC6272552 | biostudies-literature

QM/QM studies for Michael reaction in coronavirus main protease (3CL Pro).

Project description:Severe acute respiratory syndrome (SARS) is an illness caused by a novel corona virus wherein the main proteinase called 3CL(Pro) has been established as a target for drug design. The mechanism of action involves nucleophilic attack by Cys145 present in the active site on the carbonyl carbon of the scissile amide bond of the substrate and the intermediate product is stabilized by hydrogen bonds with the residues of the oxyanion hole. Based on the X-ray structure of 3CL(Pro) co-crystallized with a trans-alpha,beta-unsaturated ethyl ester (Michael acceptor), a set of QM/QM and QM/MM calculations were performed, yielding three models with increasingly higher the number of atoms. A previous validation step was performed using classical theoretical calculation and PROCHECK software. The Michael reaction studies show an exothermic process with -4.5 kcal/mol. During the reaction pathway, an intermediate is formed by hydrogen and water molecule migration from a histidine residue and solvent, respectively. In addition, similar with experimental results, the complex between N3 and 3CL(Pro) is 578 kcal/mol more stable than N1-3CL(Pro) using Own N-layer Integrated molecular Orbital molecular Mechanics (ONIOM) approach. We suggest 3CL(Pro) inhibitors need small polar groups to decrease the energy barrier for alkylation reaction. These results can be useful for the development of new compounds against SARS.

| S-EPMC7110475 | biostudies-literature

CHARMM-GUI QM/MM Interfacer for a Quantum Mechanical and Molecular Mechanical (QM/MM) Simulation Setup: 1. Semiempirical Methods.

Project description:Quantum mechanical (QM) treatments, when combined with molecular mechanical (MM) force fields, can effectively handle enzyme-catalyzed reactions without significantly increasing the computational cost. In this context, we present CHARMM-GUI QM/MM Interfacer, a web-based cyberinfrastructure designed to streamline the preparation of various QM/MM simulation inputs with ligand modification. The development of QM/MM Interfacer has been achieved through integration with existing CHARMM-GUI modules, such as PDB Reader and Manipulator, Solution Builder, and Membrane Builder. In addition, new functionalities have been developed to facilitate the one-stop preparation of QM/MM systems and enable interactive and intuitive ligand modifications and QM atom selections. QM/MM Interfacer offers support for a range of semiempirical QM methods, including AM1(+/d), PM3(+/PDDG), MNDO(+/d, +/PDDG), PM6, RM1, and SCC-DFTB, tailored for both AMBER and CHARMM. A nontrivial setup related to ligand modification, link-atom insertion, and charge distribution is automatized through intuitive user interfaces. To illustrate the robustness of QM/MM Interfacer, we conducted QM/MM simulations of three enzyme-substrate systems: dihydrofolate reductase, insulin receptor kinase, and oligosaccharyltransferase. In addition, we have created three tutorial videos about building these systems, which can be found at https://www.charmm-gui.org/demo/qmi. QM/MM Interfacer is expected to be a valuable and accessible web-based tool that simplifies and accelerates the setup process for hybrid QM/MM simulations.

| S-EPMC11209942 | biostudies-literature

Revealing the molecular mechanisms of proteolysis of SARS-CoV-2 M<sup>pro</sup> by QM/MM computational methods.

Project description:SARS-CoV-2 Mpro is one of the enzymes essential for the replication process of the virus responsible for the COVID-19 pandemic. This work is focused on exploring its proteolysis reaction by means of QM/MM methods. The resulting free energy landscape of the process provides valuable information on the species appearing along the reaction path and suggests that the mechanism of action of this enzyme, taking place in four steps, slightly differs from that of other cysteine proteases. Our predictions, which are in agreement with some recently published experimental data, can be used to guide the design of COVID-19 antiviral compounds with clinical potential.

| S-EPMC8162313 | biostudies-literature

Comparative Computational Approach To Study Enzyme Reactions Using QM and QM-MM Methods.

Project description:Choline oxidase catalyzes oxidation of choline into glycine betaine through a two-step reaction pathway employing flavin as the cofactor. On the light of kinetic studies, it is proposed that a hydride ion is transferred from α-carbon of choline/hydrated-betaine aldehyde to the N5 position of flavin in the rate-determining step, which is preceded by deprotonation of hydroxyl group of choline/hydrated-betaine aldehyde to one of the possible basic side chains. Using the crystal structure of glycine betaine-choline oxidase complex, we formulated two computational systems to study the hydride-transfer mechanism including main active-site amino acid side chains, flavin cofactor, and choline as a model system. The first system used pure density functional theory calculations, whereas the second approach used a hybrid ONIOM approach consisting of density functional and molecular mechanics calculations. We were able to formulate in silico model active sites to study the hydride-transfer steps by utilizing noncovalent chemical interactions between choline/betaine aldehyde and active-site amino acid chains using an atomistic approach. We evaluated and compared the geometries and energetics of hydride-transfer process using two different systems. We highlighted chemical interactions and studied the effect of protonation state of an active-site histidine base on the energetics of transfer. Furthermore, we evaluated energetics of the second hydride-transfer process as well as hydration of betaine aldehyde.

| S-EPMC6643517 | biostudies-literature

QM/MM Benchmarking of Cyanobacteriochrome Slr1393g3 Absorption Spectra.

Project description:Cyanobacteriochromes are compact and spectrally diverse photoreceptor proteins that are promising candidates for biotechnological applications. Computational studies can contribute to an understanding at a molecular level of their wide spectral tuning and diversity. In this contribution, we benchmark methods to model a 110 nm shift in the UV/Vis absorption spectrum from a red- to a green-absorbing form of the cyanobacteriochrome Slr1393g3. Based on an assessment of semiempirical methods to describe the chromophore geometries of both forms in vacuo, we find that DFTB2+D leads to structures that are the closest to the reference method. The benchmark of the excited state calculations is based on snapshots from quantum mechanics/molecular mechanics molecular dynamics simulations. In our case, the methods RI-ADC(2) and sTD-DFT based on CAM-B3LYP ground state calculations perform the best, whereas no functional can be recommended to simulate the absorption spectra of both forms with time-dependent density functional theory. Furthermore, the difference in absorption for the lowest energy absorption maxima of both forms can already be modelled with optimized structures, but sampling is required to improve the shape of the absorption bands of both forms, in particular for the second band. This benchmark study can guide further computational studies, as it assesses essential components of a protocol to model the spectral tuning of both cyanobacteriochromes and the related phytochromes.

| S-EPMC6540152 | biostudies-literature

QM/MM X-ray refinement of zinc metalloenzymes.

Project description:Zinc metalloenzymes play an important role in biology. However, due to the limitation of molecular force field energy restraints used in X-ray refinement at medium or low resolutions, the precise geometry of the zinc coordination environment can be difficult to distinguish from ambiguous electron density maps. Due to the difficulties involved in defining accurate force fields for metal ions, the QM/MM (quantum-mechanical/molecular-mechanical) method provides an attractive and more general alternative for the study and refinement of metalloprotein active sites. Herein we present three examples that indicate that QM/MM based refinement yields a superior description of the crystal structure based on R and R(free) values and on the inspection of the zinc coordination environment. It is concluded that QM/MM refinement is an useful general tool for the improvement of the metal coordination sphere in metalloenzyme active sites.

| S-EPMC2843799 | biostudies-literature