Differential specificity of SARS-CoV-2 main protease variants on peptide versus protein-based substrates

The SARS-CoV-2 main protease (Mpro ) holds significant importance as a biological target in combating coronaviruses due to its importance in virus replication. Considering the emergence of novel SARS-CoV-2 variants and the mutations observed in the Mpro sequence, we hypothesized that these mutations may have a potential impact on the protease's specificity. To test this, we expressed Mpro corresponding to the original strain and variants Beta1, Beta2, and Omicron and analyzed their activity on protein-based and peptide substrates. Although we observed differential activity on the protein-based substrate, there was very little difference when analyzed on the peptide substrate. We conclude that mutations on the Mpro sequence, despite having a minor effect on a peptide substrate cleavage, did not change the catalytic site environment enough to build resistance to inhibition. Therefore, we propose that inhibitors initially designed for the Mpro of the original strain will be effective in all the variants. Thus, Mpro is likely to continue to be a target of therapeutic interest as mutations in its sequence are rare and, as we show here, have a minor effect on the protease's recognition of peptide-based molecules.

Physiological models to study the effect of molecular crowding on multi-drug bound proteins: insights from SARS-CoV-2 main protease

Molecular Dynamics simulations are often used in drug design. However, such simulations do not account for the physiological environment of the receptor; hence overlook its impact on biomolecular interactions. To address this lacuna, we identified three objectives to pursue - develop models of physiological environment, study a drug-receptor complex in such environments, and identify methods to analyze these complicated simulations. Two novel physiological models were developed and studied. The first, called 'm10', comprises of 10 of the most abundant cytoplasmic metabolites at physiological concentrations. The second, called 'phy', supplements m10 with an additional crowder protein to elicit macromolecular crowding effect. The main protease (M^pro) of SARS-CoV-2, being essential for viral replication, is an attractive drug target for COVID-19. Hence, we chose M^pro docked with multiple drugs as our model drug-receptor system. With a plethora of compounds, physiological systems can be exceedingly large and complex. A novel Spark-based software (SparkTraj) was developed to rapidly analyze non-specific contacts and water interactions. Our study shows that crowding enhances the difference in the dynamics of apo- vs drug-bound complexes. Metabolites, at times as a cluster, were seen interacting with the protease, drugs, and binding sites in drug-free receptor. Except one that crawled to an adjacent pocket in phy, the drugs remained in their respective pockets in all simulations. Given these observations, we hope that the models and approach presented here would help the optimization, evaluation, and selection of potential drugs. Generic biomolecular dynamics could also benefit from such models and tools.Communicated by Ramaswamy H. Sarma.

A genetically encoded BRET-based SARS-CoV-2 Mpro protease activity sensor

The main protease, Mpro, is critical for SARS-CoV-2 replication and an appealing target for designing anti-SARS-CoV-2 agents. Therefore, there is a demand for the development of improved sensors to monitor its activity. Here, we report a pair of genetically encoded, bioluminescence resonance energy transfer (BRET)-based sensors for detecting Mpro proteolytic activity in live cells as well as in vitro. The sensors were generated by sandwiching peptides containing the Mpro N-terminal autocleavage sites, either AVLQSGFR (short) or KTSAVLQSGFRKME (long), in between the mNeonGreen and NanoLuc proteins. Co-expression of the sensors with Mpro in live cells resulted in their cleavage while mutation of the critical C145 residue (C145A) in Mpro completely abrogated their cleavage. Additionally, the sensors recapitulated the inhibition of Mpro by the well-characterized pharmacological agent GC376. Further, in vitro assays with the BRET-based Mpro sensors revealed a molecular crowding-mediated increase in the rate of Mpro activity and a decrease in the inhibitory potential of GC376. The sensors developed here will find direct utility in studies related to drug discovery targeting the SARS-CoV-2 Mpro and functional genomics application to determine the effect of sequence variation in Mpro.

Celastrol: A lead compound that inhibits SARS-CoV-2 replication, the activity of viral and human cysteine proteases, and virus-induced IL-6 secretion

The global emergence of coronavirus disease 2019 (COVID‐19) has caused substantial human casualties. Clinical manifestations of this disease vary from asymptomatic to lethal, and the symptomatic form can be associated with cytokine storm and hyperinflammation. In face of the urgent demand for effective drugs to treat COVID‐19, we have searched for candidate compounds using in silico approach followed by experimental validation. Here we identified celastrol, a pentacyclic triterpene isolated from Tripterygium wilfordii Hook F, as one of the best compounds out of 39 drug candidates. Celastrol reverted the gene expression signature from severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2)‐infected cells and irreversibly inhibited the recombinant forms of the viral and human cysteine proteases involved in virus invasion, such as M^pro (main protease), PL^pro (papain‐like protease), and recombinant human cathepsin L. Celastrol suppressed SARS‐CoV‐2 replication in human and monkey cell lines and decreased interleukin‐6 (IL‐6) secretion in the SARS‐CoV‐2‐infected human cell line. Celastrol acted in a concentration‐dependent manner, with undetectable signs of cytotoxicity, and inhibited in vitro replication of the parental and SARS‐CoV‐2 variant. Therefore, celastrol is a promising lead compound to develop new drug candidates to face COVID‐19 due to its ability to suppress SARS‐CoV‐2 replication and IL‐6 production in infected cells.

Impact of Inert Crowders on Host-Guest Recognition Process

Biological environments typically contain high concentrations (300-400 mg/mL) of different macromolecules at volume fractions as large as 30%-40%. Biomolecular recognition processes, a ubiquitous biological phenomena, occurring in such crowded heterogeneous media would differ significantly compared to the dilute buffer solutions. Here we quantify the potential impact of inert crowders on prototypical host-guest recognition process by explicit-solvent molecular dynamics (MD) simulations in atomic resolution. We demonstrate that the crowders, when smaller in size, would facilitate the binding process of the guest molecule by decreasing the free energy barrier for binding via excluded volume effect and desolvation of the host receptor. However, the extent of crowder-induced stabilization of a host-guest complex is found to be significantly higher when the guest molecule is sterically constricted to approach the host along a centrosymmetric direction, compared to its unrestricted, freely diffusive movement. A kinetic analysis of the recognition process reveals that the origin of a relatively stronger crowder impact during constricted movement of guest molecule lies in the relatively enhanced residence time of the guest inside the host by crowders. Together, our results suggest that the extent of impact of crowding on recognition processes would be contingent upon the presence or absence of constriction on ligand movement.

Protein-complex stability in cells and in vitro under crowded conditions

Biology is beginning to appreciate the effects of the crowded and complex intracellular environment on the equilibrium thermodynamics and kinetics of protein folding. The next logical step involves the interactions between proteins. We review quantitative, wet-experiment based efforts aimed at understanding how and why high concentrations of small molecules, synthetic polymers, biologically relevant cosolutes and the interior of living cells affect the energetics of protein-protein interactions. We then address popular theories used to explain the effects and suggest expeditious paths for a more methodical integration of experiment and simulation.

What macromolecular crowding can do to a protein

The intracellular environment represents an extremely crowded milieu, with a limited amount of free water and an almost complete lack of unoccupied space. Obviously, slightly salted aqueous solutions containing low concentrations of a biomolecule of interest are too simplistic to mimic the “real life” situation, where the biomolecule of interest scrambles and wades through the tightly packed crowd. In laboratory practice, such macromolecular crowding is typically mimicked by concentrated solutions of various polymers that serve as model “crowding agents”. Studies under these conditions revealed that macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, protein-protein interactions, protein-nucleic acid interactions, and pathological aggregation. The goal of this review is to systematically analyze currently available experimental data on the variety of effects of macromolecular crowding on a protein molecule. The review covers more than 320 papers and therefore represents one of the most comprehensive compendia of the current knowledge in this exciting area.

Evidence for substrate binding-induced zwitterion formation in the catalytic Cys-His dyad of the SARS-CoV main protease

The coronavirus main protease (M(pro)) represents an attractive drug target for antiviral therapy of coronavirus (CoV) infections, including severe acute respiratory syndrome (SARS). The SARS-CoV M(pro) and related CoV proteases have several distinct features, such as an uncharged Cys-His catalytic dyad embedded in a chymotrypsin-like protease fold, that clearly separate these enzymes from archetypical cysteine proteases. To further characterize the catalytic system of CoV main proteases and to obtain information about improved inhibitors, we performed comprehensive simulations of the proton-transfer reactions in the SARS-CoV M(pro) active site that lead to the Cys(-)/His(+) zwitterionic state required for efficient proteolytic activity. Our simulations, comprising the free enzyme as well as substrate-enzyme and inhibitor-enzyme complexes, lead us to predict that zwitterion formation is fostered by substrate binding but not inhibitor binding. This indicates that M(pro) employs a substrate-induced catalytic mechanism that further enhances its substrate specificity. Our computational data are in line with available experimental results, such as X-ray geometries, measured pKa values, mutagenesis experiments, and the measured differences between the kinetic parameters of substrates and inhibitors. The data also provide an atomistic picture of the formerly postulated electrostatic trigger involved in SARS-CoV M(pro) activity. Finally, they provide information on how a specific microenvironment may finely tune the activity of M(pro) toward specific viral protein substrates, which is known to be required for efficient viral replication. Our simulations also indicate that the low inhibition potencies of known covalently interacting inhibitors may, at least in part, be attributed to insufficient fostering of the proton-transfer reaction. These findings suggest ways to achieve improved inhibitors.

P-I class metalloproteinase from Bothrops moojeni venom is a post-proline cleaving peptidase with kininogenase activity: insights into substrate selectivity and kinetic behavior

Snake venom metalloproteinases (SVMPs) belonging to P-I class are able to hydrolyze extracellular matrix proteins and coagulation factors triggering local and systemic reactions by multiple molecular mechanisms that are not fully understood. BmooMPα-I, a P-I class SMVP from Bothrops moojeni venom, was active upon neuro- and vaso-active peptides including angiotensin I, bradykinin, neurotensin, oxytocin and substance P. Interestingly, BmooMPα-I showed a strong bias towards hydrolysis after proline residues, which is unusual for most of characterized peptidases. Moreover, the enzyme showed kininogenase activity similar to that observed in plasma and cells by kallikrein. FRET peptide assays indicated a relative promiscuity at its S2-S'2 subsites, with proline determining the scissile bond. This unusual post-proline cleaving activity was confirmed by the efficient hydrolysis of the synthetic combinatorial library MCA-GXXPXXQ-EDDnp, described as resistant for canonical peptidases, only after Pro residues. Structural analysis of the tripeptide LPL complexed with BmooMPα-I, generated by molecular dynamics simulations, assisted in defining the subsites and provided the structural basis for subsite preferences such as the restriction of basic residues at the S2 subsite due to repulsive electrostatic effects and the steric impediment for large aliphatic or aromatic side chains at the S1 subsite. These new functional and structural findings provided a further understanding of the molecular mechanisms governing the physiological effects of this important class of enzymes in envenomation process.

Chimeric exchange of coronavirus nsp5 proteases (3CLpro) identifies common and divergent regulatory determinants of protease activity

Human coronaviruses (CoVs) such as severe acute respiratory syndrome CoV (SARS-CoV) and Middle East respiratory syndrome CoV (MERS-CoV) cause epidemics of severe human respiratory disease. A conserved step of CoV replication is the translation and processing of replicase polyproteins containing 16 nonstructural protein domains (nsp's 1 to 16). The CoV nsp5 protease (3CLpro; Mpro) processes nsp's at 11 cleavage sites and is essential for virus replication. CoV nsp5 has a conserved 3-domain structure and catalytic residues. However, the intra- and intermolecular determinants of nsp5 activity and their conservation across divergent CoVs are unknown, in part due to challenges in cultivating many human and zoonotic CoVs. To test for conservation of nsp5 structure-function determinants, we engineered chimeric betacoronavirus murine hepatitis virus (MHV) genomes encoding nsp5 proteases of human and bat alphacoronaviruses and betacoronaviruses. Exchange of nsp5 proteases from HCoV-HKU1 and HCoV-OC43, which share the same genogroup, genogroup 2a, with MHV, allowed for immediate viral recovery with efficient replication albeit with impaired fitness in direct competition with wild-type MHV. Introduction of MHV nsp5 temperature-sensitive mutations into chimeric HKU1 and OC43 nsp5 proteases resulted in clear differences in viability and temperature-sensitive phenotypes compared with MHV nsp5. These data indicate tight genetic linkage and coevolution between nsp5 protease and the genomic background and identify differences in intramolecular networks regulating nsp5 function. Our results also provide evidence that chimeric viruses within coronavirus genogroups can be used to test nsp5 determinants of function and inhibition in common isogenic backgrounds and cell types.

Formation of protein complexes in crowded environments--from in vitro to in vivo

Traditionally, biochemical studies are performed in dilute homogenous solutions, which are very different from the dense mixture of molecules found in cells. Thus, the physiological relevance of these studies is in question. This recognition motivated scientists to formulate the effect of crowded solutions in general, and excluded volume in particular, on biochemical processes. Using polymers or proteins as crowders, it was shown that while crowding tends to significantly enhance the formation of complexes containing many subunits, dimerizations are only mildly affected. Computer simulations, together with experimental evidence, indicate soft interactions and diffusion as critical factors that operate in a concerted manner with excluded volume to modulate protein binding. Yet, these approaches do not truly mimic the cellular environment. In vivo studies may overcome this shortfall. The few studies conducted thus far suggest that in cells, binding and folding occur at rates close to those determined in dilute solutions. Obtaining quantitative biochemical information on reactions inside living cells is currently a main challenge of the field, as the complexity of the intracellular milieu was what motivated crowding research to begin with.

Structured crowding and its effects on enzyme catalysis

Macromolecular crowding decreases the diffusion rate, shifts the equilibrium of protein-protein and protein-substrate interactions, and changes protein conformational dynamics. Collectively, these effects contribute to enzyme catalysis. Here we describe how crowding may bias the conformational change and dynamics of enzyme populations and in this way affect catalysis. Crowding effects have been studied using artificial crowding agents and in vivo-like environments. These studies revealed a correlation between protein dynamics and function in the crowded environment. We suggest that crowded environments be classified into uniform crowding and structured crowding. Uniform crowding represents random crowding conditions created by synthetic particles with a narrow size distribution. Structured crowding refers to the highly coordinated cellular environment, where proteins and other macromolecules are clustered and organized. In structured crowded environments the perturbation of protein thermal stability may be lower; however, it may still be able to modulate functions effectively and dynamically. Dynamic, allosteric enzymes could be more sensitive to cellular perturbations if their free energy landscape is flatter around the native state; on the other hand, if their free energy landscape is rougher, with high kinetic barriers separating deep minima, they could be more robust. Above all, cells are structured; and this holds both for the cytosol and for the membrane environment. The crowded environment is organized, which limits the search, and the crowders are not necessarily inert. More likely, they too transmit allosteric effects, and as such play important functional roles. Overall, structured cellular crowding may lead to higher enzyme efficiency and specificity.