Revealing the molecular mechanisms of proteolysis of SARS-CoV-2 Mpro by QM/MM computational methods.
Ontology highlight
ABSTRACT:
SUBMITTER: Swiderek K
PROVIDER: S-EPMC8162313 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
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
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
Project description:The catalytic mechanism of the decarboxylation of 5-carboxyvanillate by LigW producing vanillic acid has been studied by using QM cluster and hybrid QM/MM methodologies. In the QM cluster model, the environment of a small QM model is treated with a bulky potential while two QM/MM models studies include partial and full protein with and without explicitly treated water solvent. The studied reaction involves two sequential steps: the protonation of the carbon of the 5-carboxy-vanillate substrate and the decarboxylation of the intermediate from which results deprotonated vanillic acid as product. The structures and energetics obtained by using three structural models and two density functionals are quite consistent to each other. This indicates that the small QM cluster model of the presently considered enzymatic reaction is appropriate enough and the reaction is mainly influenced by the active site.
| S-EPMC6031855 | biostudies-literature
Project description:Mixed quantum mechanics/molecular mechanics (QM/MM) methods offer a valuable computational tool for understanding the electron transfer pathway in protein-substrate interactions and protein-protein complexes. These hybrid methods are capable of solving the Schrödinger equation on a small subset of the protein, the quantum region, describing its electronic structure under the polarization effects of the remainder of the protein. By selectively turning on and off different residues in the quantum region, we are able to obtain the electron pathway for short- and large-range interactions. Here, we summarize recent studies involving the protein-substrate interaction in cytochrome P450 camphor, ascorbate peroxidase and cytochrome c peroxidase, and propose a novel approach for the long-range protein-protein electron transfer. The results on ascorbate peroxidase and cytochrome c peroxidase reveal the importance of the propionate groups in the electron transfer pathway. The long-range protein-protein electron transfer has been studied on the cytochrome c peroxidase-cytochrome c complex. The results indicate the importance of Phe82 and Cys81 on cytochrome c, and of Asn196, Ala194, Ala176 and His175 on cytochrome c peroxidase.
| S-EPMC2706105 | biostudies-literature
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
Project description: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.
| S-EPMC8004986 | biostudies-literature
Project description:To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.
| S-EPMC4688044 | biostudies-literature
Project description:The SARS-CoV-2 main protease (Mpro) 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 Mpro 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 Mpro, 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 Mpro 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.
| S-EPMC8179034 | biostudies-literature
Project description:Programmed ribosomal frameshifting (PRF) serves as an intrinsic translational regulation mechanism employed by some viruses to control the ratio between structural and enzymatic proteins. Most viral mRNAs which use PRF adapt an H-type pseudoknot to stimulate -1 PRF. The relationship between the thermodynamic stability and the frameshifting efficiency of pseudoknots has not been fully understood. Recently, single-molecule force spectroscopy has revealed that the frequency of -1 PRF correlates with the unwinding forces required for disrupting pseudoknots, and that some of the unwinding work dissipates irreversibly due to the torsional restraint of pseudoknots. Complementary to single-molecule techniques, computational modeling provides insights into global motions of the ribosome, whose structural transitions during frameshifting have not yet been elucidated in atomic detail. Taken together, recent advances in biophysical tools may help to develop antiviral therapies that target the ubiquitous -1 PRF mechanism among viruses.
| S-EPMC3321566 | biostudies-literature
Project description:We have made a systematic combined quantum mechanical and molecular mechanical (QM/MM) investigation of possible structures of the N2 bound state of nitrogenase. We assume that N2 is immediately protonated to a N2H2 state, thereby avoiding the problem of determining the position of the protons in the cluster. We have systematically studied both end-on and side-on structures, as well as both HNNH and NNH2 states. Our results indicate that the binding of N2H2 is determined more by interactions and steric clashes with the surrounding protein than by the intrinsic preferences of the ligand and the cluster. The best binding mode with both the TPSS and B3LYP density-functional theory methods has trans-HNNH terminally bound to Fe2. It is stabilised by stacking of the substrate with His-195 and Ser-278. However, several other structures come rather close in energy (within 3-35 kJ/mol) at least in some calculations: The corresponding cis-HNNH structure terminally bound to Fe2 is second best with B3LYP. A structure with HNNH2 terminally bound to Fe6 is second most stable with TPSS (where the third proton is transferred to the substrate from the homocitrate ligand). Structures with trans-HNNH, bound to Fe4 or Fe6, or cis-HNNH bound to Fe6 are also rather stable. Finally, with the TPSS functional, a structure with cis-HNNH side-on binding to the Fe3-Fe4-Fe5-Fe7 face of the cluster is also rather low in energy, but all side-on structures are strongly disfavoured by the B3LYP method.
| S-EPMC7186253 | biostudies-literature