Visualizing the active site oxyanion loop transition upon ensitrelvir binding and transient dimerization of SARS-CoV-2 main protease

N-terminal autoprocessing from its polyprotein precursor enables creating the mature-like stable dimer interface of SARS-CoV-2 main protease (MPro), concomitant with the active site oxyanion loop equilibrium transitioning to the active conformation (E*) and onset of catalytic activity. Through mutagenesis of critical interface residues and evaluating noncovalent inhibitor (ensitrelvir, ESV) facilitated dimerization through its weak binding to MPro, we demonstrate that residues extending from Ser1 through Glu14 are critical for dimerization. Combined mutations G11A, E290A and R298A (MProTM) restrict dimerization even upon binding of ESV to monomeric MProTM with an inhibitor dissociation constant of 7.4 ± 1.6 µM. Contrasting the covalent inhibitor NMV or GC373 binding to monomeric MPro, ESV binding enabled capturing the transition of the oxyanion loop conformations in the absence of a reactive warhead and independent of dimerization. Characterization of complexes by room-temperature X-ray crystallography reveals ESV bound to the E* state of monomeric MPro as well as an intermediate approaching the inactive state (E). It appears that the E* to E equilibrium shift occurs initially from G138-F140 residues, leading to the unwinding of the loop and formation of the 310-helix. Finally, we describe a transient dimer structure of the MPro precursor held together through interactions of residues A5-G11 with distinct states of the active sites, E and E*, likely representing an intermediate in the autoprocessing pathway.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MPro^WT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373) and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decreases in free energy (ΔG) and thermal stability (T_m), respectively, relative to MPro^WT. The H41A mutation results in a 20-fold increase in the dissociation constant (K_d), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and T_m, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MPro^H41A, whereas no significant change is observed in binding to MPro^WT. Structures of the four inhibitor complexes with MPro^1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MPro^WT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MPro^WT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MPro^C145A also results in a significant, but smaller, increase in K_d and decrease in ΔG and T_m, relative to NMV.

Unmasking the Conformational Stability and Inhibitor Binding to SARS-CoV-2 Main Protease Active Site Mutants and Miniprecursor

We recently demonstrated that inhibitor binding reorganizes the oxyanion loop of a monomeric catalytic domain of SARS CoV-2 main protease (MPro) from an unwound (E) to a wound (active, E*) conformation, independent of dimerization. Here we assess the effect of the flanking N-terminal residues, to imitate the MPro precursor prior to its autoprocessing, on conformational equilibria rendering stability and inhibitor binding. Thermal denaturation (T_m) of C145A mutant, unlike H41A, increases by 6.8 °C, relative to wild-type mature dimer. An inactivating H41A mutation to maintain a miniprecursor containing TSAVL[Q or E] of the flanking nsp4 sequence in an intact form [^(-6)MPro^H41A and ^(-6*)MPro^H41A, respectively], and its corresponding mature MPro^H41A were systematically examined. While the H41A mutation exerts negligible effect on T_m and dimer dissociation constant (K_dimer) of MPro^H41A, relative to the wild type MPro, both miniprecursors show a 4-5 °C decrease in T_m and > 85-fold increase in K_dimer as compared to MPro^H41A. The K_d for the binding of the covalent inhibitor GC373 to ^(-6*)MPro^H41A increases ∼12-fold, relative to MPro^H41A, concomitant with its dimerization. While the inhibitor-free dimer exhibits a state in transit from E to E* with a conformational asymmetry of the protomers' oxyanion loops and helical domains, inhibitor binding restores the asymmetry to mature-like oxyanion loop conformations (E*) but not of the helical domains. Disorder of the terminal residues 1-2 and 302-306 observed in both structures suggest that N-terminal autoprocessing is tightly coupled to the E-E* equilibrium and stable dimer formation.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Å resolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Room-temperature X-ray crystallography reveals the oxidation and reactivity of cysteine residues in SARS-CoV-2 3CL Mpro: insights into enzyme mechanism and drug design

The emergence of the novel coronavirus SARS-CoV-2 has resulted in a worldwide pandemic not seen in generations. Creating treatments and vaccines to battle COVID-19, the disease caused by the virus, is of paramount importance in order to stop its spread and save lives. The viral main protease, 3CL Mpro, is indispensable for the replication of SARS-CoV-2 and is therefore an important target for the design of specific protease inhibitors. Detailed knowledge of the structure and function of 3CL Mpro is crucial to guide structure-aided and computational drug-design efforts. Here, the oxidation and reactivity of the cysteine residues of the protease are reported using room-temperature X-ray crystallography, revealing that the catalytic Cys145 can be trapped in the peroxysulfenic acid oxidation state at physiological pH, while the other surface cysteines remain reduced. Only Cys145 and Cys156 react with the alkylating agent N-ethylmaleimide. It is suggested that the zwitterionic Cys145-His45 catalytic dyad is the reactive species that initiates catalysis, rather than Cys145-to-His41 proton transfer via the general acid-base mechanism upon substrate binding. The structures also provide insight into the design of improved 3CL Mpro inhibitors.

Flexibility and Design: Conformational Heterogeneity along the Evolutionary Trajectory of a Redesigned Ubiquitin

Although protein design has been used to introduce new functions, designed variants generally only function as well as natural proteins after rounds of laboratory evolution. One possibility for this pattern is that designed mutants frequently sample nonfunctional conformations. To test this idea, we exploited advances in multiconformer modeling of room-temperature X-ray data collection on redesigned ubiquitin variants selected for increasing binding affinity to the deubiquitinase USP7. Initial core mutations disrupt natural packing and lead to increased flexibility. Additional, experimentally selected mutations quenched conformational heterogeneity through new stabilizing interactions. Stabilizing interactions, such as cation-pi stacking and ordered waters, which are not included in standard protein design energy functions, can create specific interactions that have long-range effects on flexibility across the protein. Our results suggest that increasing flexibility may be a useful strategy to escape local minima during initial directed evolution and protein design steps when creating new functions.