Browse
Submit Data
Databases
API
Help

Dataset Information

15 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

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

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

Impact of Warhead Modulations on the Covalent Inhibition of SARS-CoV-2 M<sup>pro</sup> Explored by QM/MM Simulations.

Project description:The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus-2, SARS-CoV-2, shows the need for effective antiviral treatments. Here, we present a simulation study of the inhibition of the SARS-CoV-2 main protease (Mpro), a cysteine hydrolase essential for the life cycle of the virus. The free energy landscape for the mechanism of the inhibition process is explored by QM/MM umbrella sampling and free energy perturbation simulations at the M06-2X/MM level of theory for two proposed peptidyl covalent inhibitors that share the same recognition motif but feature distinct cysteine-targeting warheads. Regardless of the intrinsic reactivity of the modeled inhibitors, namely a Michael acceptor and a hydroxymethyl ketone activated carbonyl, our results confirm that the inhibitory process takes place by means of a two-step mechanism, in which the formation of an ion pair C145/H41 dyad precedes the protein-inhibitor covalent bond formation. The nature of this second step is strongly dependent on the functional groups in the warhead: while the nucleophilic attack of the C145 sulfur atom on the Cα of the double bond of the Michael acceptor takes place concertedly with the proton transfer from H41 to Cβ, in the compound with an activated carbonyl, the sulfur attacks the carbonyl carbon concomitant with a proton transfer from H41 to the carbonyl oxygen via the hydroxyl group. An analysis of the free energy profiles, structures along the reaction path, and interactions between the inhibitors and the different pockets of the active site on the protein shows a measurable effect of the warhead on the kinetics and thermodynamics of the process. These results and QM/MM methods can be used as a guide to select warheads to design efficient irreversible and reversible inhibitors of SARS-CoV-2 Mpro.

| S-EPMC8751016 | 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

Projected hybrid orbitals: a general QM/MM method.

Project description:A projected hybrid orbital (PHO) method was described to model the covalent boundary in a hybrid quantum mechanical and molecular mechanical (QM/MM) system. The PHO approach can be used in ab initio wave function theory and in density functional theory with any basis set without introducing system-dependent parameters. In this method, a secondary basis set on the boundary atom is introduced to formulate a set of hybrid atomic orbtials. The primary basis set on the boundary atom used for the QM subsystem is projected onto the secondary basis to yield a representation that provides a good approximation to the electron-withdrawing power of the primary basis set to balance electronic interactions between QM and MM subsystems. The PHO method has been tested on a range of molecules and properties. Comparison with results obtained from QM calculations on the entire system shows that the present PHO method is a robust and balanced QM/MM scheme that preserves the structural and electronic properties of the QM region.

| S-EPMC4306490 | biostudies-literature