From quantum-derived principles underlying cysteine reactivity to combating the COVID-19 pandemic

The COVID-19 pandemic poses a challenge in coming up with quick and effective means to counter its cause, the SARS-CoV-2. Here, we show how the key factors governing cysteine reactivity in proteins derived from combined quantum mechanical/continuum calculations led to a novel multi-targeting strategy against SARS-CoV-2, in contrast to developing potent drugs/vaccines against a single viral target such as the spike protein. Specifically, they led to the discovery of reactive cysteines in evolutionary conserved Zn2+-sites in several SARS-CoV-2 proteins that are crucial for viral polypeptide proteolysis as well as viral RNA synthesis, proofreading, and modification. These conserved, reactive cysteines, both free and Zn2+-bound, can be targeted using the same Zn-ejector drug (disulfiram/ebselen), which enables the use of broad-spectrum anti-virals that would otherwise be removed by the virus's proofreading mechanism. Our strategy of targeting multiple, conserved viral proteins that operate at different stages of the virus life cycle using a Zn-ejector drug combined with other broad-spectrum anti-viral drug(s) could enhance the barrier to drug resistance and antiviral effects, as compared to each drug alone. Since these functionally important nonstructural proteins containing reactive cysteines are highly conserved among coronaviruses, our proposed strategy has the potential to tackle future coronaviruses. This article is categorized under:Structure and Mechanism > Reaction Mechanisms and CatalysisStructure and Mechanism > Computational Biochemistry and BiophysicsElectronic Structure Theory > Density Functional Theory.

Promiscuous Catalytic Activity of a Binuclear Metallohydrolase: Peptide and Phosphoester Hydrolyses

In this study, chemical promiscuity of a binuclear metallohydrolase Streptomyces griseus aminopeptidase (SgAP) has been investigated using DFT calculations. SgAP catalyzes two diverse reactions, peptide and phosphoester hydrolyses, using its binuclear (Zn-Zn) core. On the basis of the experimental information, mechanisms of these reactions have been investigated utilizing leucine p-nitro aniline (Leu-pNA) and bis(4-nitrophenyl) phosphate (BNPP) as the substrates. The computed barriers of 16.5 and 16.8 kcal/mol for the most plausible mechanisms proposed by the DFT calculations are in good agreement with the measured values of 13.9 and 18.3 kcal/mol for the Leu-pNA and BNPP hydrolyses, respectively. The former was found to occur through the transfer of two protons, while the latter with only one proton transfer. They are in line with the experimental observations. The cleavage of the peptide bond was the rate-determining process for the Leu-pNA hydrolysis. However, the creation of the nucleophile and its attack on the electrophile phosphorus atom was the rate-determining step for the BNPP hydrolysis. These calculations showed that the chemical nature of the substrate and its binding mode influence the nucleophilicity of the metal bound hydroxyl nucleophile. Additionally, the nucleophilicity was found to be critical for the Leu-pNA hydrolysis, whereas double Lewis acid activation was needed for the BNPP hydrolysis. That could be one of the reasons why peptide hydrolysis can be catalyzed by both mononuclear and binuclear metal cofactors containing hydrolases, while phosphoester hydrolysis is almost exclusively by binuclear metallohydrolases. These results will be helpful in the development of versatile catalysts for chemically distinct hydrolytic reactions.

How the Local Environment of Functional Sites Regulates Protein Function

Proteins form complex biological machineries whose functions in the cell are highly regulated at both the cellular and molecular levels. Cellular regulation of protein functions involves differential gene expressions, post-translation modifications, and signaling cascades. Molecular regulation, on the other hand, involves tuning an optimal local protein environment for the functional site. Precisely how a protein achieves such an optimal environment around a given functional site is not well understood. Herein, by surveying the literature, we first summarize the various reported strategies used by certain proteins to ensure their correct functioning. We then formulate three key physicochemical factors for regulating a protein's functional site, namely, (i) its immediate interactions, (ii) its solvent accessibility, and (iii) its conformational flexibility. We illustrate how these factors are applied to regulate the functions of free/metal-bound Cys and Zn sites in proteins.

Complexes of Zn(II)-Triazoles with CO2 and H2O: Structures, Energetics, and Applications

Using a first-principle methodology, we investigate the stable structures of the nonreactive and reactive clusters formed between Zn²⁺-triazoles ([Zn²⁺-Tz]) clusters and CO₂ and/or H₂O. In sum, we characterized two modes of bonding of [Zn²⁺-Tz] with CO₂/H₂O: the interaction is established through (i) a covalent bond between Zn²⁺ of [Zn²⁺-Tz] and oxygen atoms of CO₂ or H₂O and (ii) hydrogen bonds through N-H or C-H of [Zn²⁺-Tz] and oxygen atoms of H₂O or CO₂, N-H···O. We also identified intramolecular proton transfer processes induced by complexation. Indeed, water drastically changes the shape of the energy profiles of the tautomeric phenomena through strong lowering of the potential barriers to tautomerism. The comparison to [Zn²⁺-Im] subunits formed with Zn²⁺ and imidazole shows that the efficiency of Tz-based compounds for CO₂ capture and uptake is due to the incorporation of more accessible nitrogen donor sites in Tzs compared to imidazoles. Since [Zn²⁺-Tz] clusters are subunits of an organometallic nanoporous materials and Zn-proteins, our data are useful for deriving force fields for macromolecular simulations of these materials. Our work also suggests the consideration of traces of water to better model the CO₂ sequestration and reactivity on macromolecular entities such as pores or active sites.

An efficient protocol for computing the pKa of Zn-bound water

We present an efficient and accurate method for computing absolute pK_w values in Zn²⁺ complexes.