Revealing the molecular mechanisms of proteolysis of SARS-CoV-2 Mpro by QM/MM computational methods

SARS-CoV-2 M^pro 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.

The molecular mechanism of the proteolysis reaction catalyzed by SARS-CoV-2 M^pro, one of the enzymes essential for the replication process of the virus responsible for the COVID-19 pandemic, is described using computational QM/MM methods.     

On the nature of the transition state in catechol O-methyltransferase. A complementary study based on molecular dynamics and potential energy surface explorations

The way in which enzymes influence the rate of chemical processes is still a question of debate. The protein promotes the catalysis of biochemical processes by lowering the free energy barrier in comparison with the reference uncatalyzed reaction in solution. In this article we are reporting static and dynamic aspects of the enzyme catalysis in a bimolecular reaction, namely a methyl transfer from S-adenosylmethionine to the hydroxylate oxygen of a substituted catechol catalyzed by catechol O-methyltransferase. From QM/MM optimizations, we will first analyze the participation of the environment on the transition vector. The study of molecular dynamics trajectories will allow us to estimate the transmission coefficient from a previously localized transition state as the maximum in the potential of mean force profile. The analysis of the reactive and nonreactive trajectories in the enzyme environment and in solution will also allow studying the geometrical and electronic changes, with special attention to the chemical system movements and the coupling with the environment. The main result, coming from both analyses, is the approximation of the magnesium cation to the nucleophilic and the hydroxyl group of the catecholate as a result of a general movement of the protein, stabilizing in this way the transition state. Consequently, the free energy barrier of the enzyme reaction is dramatically decreased with respect to the reaction in solution.     

Preorganization and reorganization as related factors in enzyme catalysis: the chorismate mutase case

In this paper a deeper insight into the chorismate-to prephenate-rearrangement, catalyzed by Bacillus subtilis chorismate mutase, is provided by means of a combination of statistical quantum mechanics/molecular mechanics simulation methods and hybrid potential energy surface exploration techniques. The main aim of this work is to present an estimation of the preorganization and reorganization terms of the enzyme catalytic rate enhancement. To analyze the first of these, we have studied different conformational equilibria of chorismate in aqueous solution and in the enzyme active site. Our conclusion is that chorismate mutase preferentially binds the reactive conformer of the substrate--that presenting a structure similar to the transition state of the reaction to be catalyzed--with shorter distances between the carbon atoms to be bonded and more diaxial character. With respect to the reorganization effect, an energy decomposition analysis of the potential energies of the reactive reactant and of the reaction transition state in aqueous solution and in the enzyme shows that the enzyme structure is better adapted to the transition structure. This means not only a more negative electrostatic interaction energy with the transition state but also a low enzyme deformation contribution to the energy barrier. Our calculations reveal that the structure of the enzyme is responsible for stabilizing the transition state structure of the reaction, with concomitant selection of the reactive form of the reactants. This is, the same enzymatic pattern that stabilizes the transition structure also promotes those reactant structures closer to the transition structure (i.e., the reactive reactants). In fact, both reorganization and preorganization effects have to be considered as the two faces of the same coin, having a common origin in the effect of the enzyme structure on the energy surface of the substrate.     

On Transition Structures for Hydride Transfer Step in Enzyme Catalysis. A Comparative Study on Models of Glutathione Reductase Derived from Semiempirical, HF, and DFT Methods

As a model of the chemical reactions that take place in the active site of gluthatione reductase, the nature of the molecular mechanism for the hydride transfer step has been characterized by means of accurate quantum chemical characterizations of transition structures. The calculations have been carried out with analytical gradients at AM1 and PM3 semiempirical procedures, ab initio at HF level with 3-21G, 4-31G, 6-31G, and 6-31G basis sets and BP86 and BLYP as density functional methods. The results of this study suggest that the endo relative orientation on the substrate imposed by the active site is optimal in polarizing the C4-Ht bond and situating the system in the neighborhood of the quadratic region of the transition structure associated to the hydride transfer step on potential energy surface. The endo arrangement of the transition structure results in optimal frontier HOMO orbital interaction between NADH and FAD partners. The geometries of the transition structures and the corresponding transition vectors, that contain the fundamental information relating reactive fluctuation patterns, are model independent and weakly dependent on the level of theory used to determine them. A comparison between simple and complex molecular models shows that there is a minimal set of coordinates describing the essentials of hydride transfer step. The analysis of transition vector components suggests that the primary and secondary kinetic isotope effects can be strongly coupled, and this prompted the calculation of deuterium and tritium primary, secondary, and primary and secondary kinetic isotope effects. The results obtained agree well with experimental data and demonstrate this coupling.     

Theoretical modeling of enzyme catalytic power: analysis of "cratic" and electrostatic factors in catechol O-methyltransferase

A comparative theoretical study of a bimolecular reaction in aqueous solution and catalyzed by the enzyme catechol O-methyltransferase (COMT) has been carried out by a combination of two hybrid QM/MM techniques: statistical simulation methods and internal energy minimizations. In contrast to previous studies by other workers, we have located and characterized transition structures for the reaction in the enzyme active site, in water and in a vacuum, and our potential of mean force calculations are based upon reaction coordinates obtained from features of the potential energy surfaces in the condensed media, not from the gas phase. The AM1/CHARMM calculated free energy of activation for the reaction of S-adenosyl methionine (SAM) with catecholate catalyzed by COMT is 15 kcal mol(-1) lower the AM1/TIP3P free-energy barrier for the reaction of the trimethylsulfonium cation with the catecholate anion in water at 300 K, in agreement with previous estimates. The thermodynamically preferred form of the reactants in the uncatalyzed model reaction in water is a solvent-separated ion pair (SSIP). Conversion of the SSIP into a contact ion pair, with a structure resembling that of the Michaelis complex (MC) for the reaction in the COMT active site, is unfavorable by 7 kcal mol(-1), largely due to reorganization of the solvent. We have considered alternative ways to estimate the so-called "cratic" free energy for bringing the reactant species together in the correct orientation for reaction but conclude that direct evaluation of the free energy of association by means of molecular dynamics simulation with a simple standard-state correction is probably the best approach. The latter correction allows for the fact that the size of the unit cell employed with the periodic boundary simulations does not correspond to the standard state concentration of 1 M. Consideration of MC-like species allows a helpful decomposition of the catalytic effect into preorganization and reorganization phases. In the preorganization phase, the substrates are brought together into the MC-like species, either in water or in the enzyme active site. In the reorganization phase, the roles of the enzymic and aqueous environments may be compared directly because reorganization of the substrate is about the same in both cases. Analysis of the electric field along the reaction coordinate demonstrates that in water the TS is destabilized with respect to the MC-like species because the polarity of the solute diminishes and consequently the reaction field is also decreased. In the enzyme, the electric field is mainly a permanent field and consequently there is only a small reorganization of the environment. Therefore, destabilization of the TS is lower than in solution, and the activation barrier is smaller.     

Carbohydrate-based DNA ligands: sugar-oligoamides as a tool to study carbohydrate-nucleic acid interactions

Sugar-oligoamides have been designed and synthesized as structurally simple carbohydrate-based ligands to study carbohydrate-DNA interactions. The general design of the ligands 1-3 has been done as to favor the bound conformation of Distamycin-type gamma-linked covalent dimers which is a hairpin conformation. Indeed, NMR analysis of the sugar-oligoamides in the free state has indicated the presence of a percentage of a hairpin conformation in aqueous solution. The DNA binding activity of compounds 1-3 was confirmed by calf thymus DNA (ct-DNA) NMR titration. Interestingly, the binding of the different sugar-oligoamides seems to be modulated by the sugar configuration. Semiquantitative structural information about the DNA ligand complexes has been derived from NMR data. A competition experiment with Netropsin suggested that the sugar-oligoamide 3 bind to DNA in the minor groove. The NMR titrations of 1-3 with poly(dA-dT) and poly(dG-dC) suggested preferential binding to the ATAT sequence. TR-NOE NMR experiments for the sugar-oligoamide 3-ct-DNA complex both in D(2)O and H(2)O have confirmed the complex formation and given information on the conformation of the ligand in the bound state. The data confirmed that the sugar-oligoamide ligand is a hairpin in the bound state. Even more relevant to our goal, structural information on the conformation around the N-glycosidic linkage has been accessed. Thus, the sugar asymmetric centers pointing to the NH-amide and N-methyl rims of the molecule have been characterized.     

Inclusion of nickel(<Emphasis Type="SmallCaps">II</Emphasis>) and copper(<Emphasis Type="SmallCaps">II</Emphasis>) complexes with aliphatic polyamines in cucurbit[8]uril

The inclusion compound of macrocyclic cavitand cucurbit[8]uril (CB[8]) with the nickel(II) complex containing the tetraazamacrocyclic ligand cyclam, {[Ni(cyclam)]@CB[8]}Cl₂··16H₂O (1), and the inclusion compounds of CB[8] with the copper(II) bis-ethylene-diamine complex, {trans-[Cu(en)₂(H₂O)₂]@CB[8]}Cl₂·{CB[8]}·42H₂O (2a) and {trans-[Cu(en)₂(H₂O)₂]@CB[8]}Cl₂·17H₂O (2b), were synthesized and characterized by X-ray diffraction analysis, IR and ESR spectroscopy, and electrospray mass spectrometry. Guest—host inclusion compounds can be directly synthesized starting from a metal complex and cucurbit[8]uril, as was exemplified by the preparation of compounds 2a and 2b.__________Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 2414–2419, November, 2004.     

Synthesis and structure of the incomplete cuboidal clusters [W3Se4H3(dmpe)3]+, [W3Se4H(3-x)(OH)x(dmpe)3]+ and [W3Se4(OH)3(dmpe)3]+, and the mechanism of the acid-assisted substitution of the coordinated hydrides

The novel incomplete cuboidal cluster [W3Se4H3(dmpe)3](PF6), [1](PF6), has been prepared by reduction of [W3Se4Br3(dmpe)3](PF6) with LiBH4 in THF solution. The trihydroxo complex [W3Se4(OH)3(dmpe)3](PF6), [2](PF6), was obtained by reacting [W3Se4Br3(dmpe)3](PF6) with NaOH in MeCN-H2O solution. The complexes [1](PF6) and [2](PF6) were converted to their BPh4- salts by treatment with NaBPh4. Recrystallisation of [1](BPh4) in the presence of traces of water affords the mixed dihydride hydroxo complex [W3Se4H2(OH)(dmpe)3](BPh4). The crystal structures of [1](BPh4), [2](BPh4) and [W3Se4H2(OH)(dmpe)3](BPh4) have been resolved. Although the [1]+ trihydride does not react with an excess of halide salts, reaction with HX leads to [W3Se4X3(dmpe)3]+ (X = Cl, Br). The kinetics of this reaction has been studied at 25 degrees C in MeCN-H2O solution (1:1, v/v) and found to occur with two consecutive kinetic steps. The first step is independent of the nature and concentration of the X(-) anion but shows a first order dependence on the concentration of acid (k1 = 12.0 mol(-1) dm(3) s(-1)), whereas the second one is independent of the nature and concentration of both the acid and added salts (k2 = 0.024 s(-1)). In contrast, the reaction of [2]+ with acids occurs in a single step with kobs = 0.63 s(-1)(HCl) and 0.17 s(-1)(HBr). These kinetic results are discussed on the basis of the mechanism previously proposed for the reactions of the analogous [W3S4H3(dmpe)3]+ cluster, with special emphasis on the effects caused by the change of S by Se on the rate constants for the different processes involved.     

Interface Amorphization of Two-Dimensional Black Phosphorus upon Treatment with Diazonium Salts

Two-dimensional (2D) black phosphorus (BP) represents one of the most appealing 2D materials due to its electronic, optical, and chemical properties. Many strategies have been pursued to face its environmental instability, covalent functionalization being one of the most promising. However, the extremely low functionalization degrees and the limitations in proving the nature of the covalent functionalization still represent challenges in many of these sheet architectures reported to date. Here we shine light on the structural evolution of 2D-BP upon the addition of electrophilic diazonium salts. We demonstrated the absence of covalent functionalization in both the neutral and the reductive routes, observing in the latter case an unexpected interface conversion of BP to red phosphorus (RP), as characterized by Raman, ³¹P-MAS NMR, and X-ray photoelectron spectroscopies (XPS). Furthermore, thermogravimetric analysis coupled to gas chromatography and mass spectrometry (TG-GC-MS), as well as electron paramagnetic resonance (EPR) gave insights into the potential underlying radical mechanism, suggesting a Sandmeyer-like reaction.