PROFICS: A bacterial selection system for directed evolution of proteases.

Proteases serve as important tools in biotechnology and as valuable drugs or drug targets. Efficient protein engineering methods to study and modulate protease properties are thus of great interest for a plethora of applications. We established PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection), a bacterial selection system, which enables the optimization of proteases for biotechnology, therapeutics or diagnosis in a simple overnight process. During the PROFICS process, proteases are selected for their ability to specifically cut a tag from a reporter enzyme and leave a native N-terminus. Precise and efficient cleavage after the recognition sequence reverses the phenotype of an Escherichia coli knockout strain deficient in an essential enzyme of pyrimidine synthesis. A toolbox was generated to select for proteases with different preferences for P1' residues (the residue immediately following the cleavage site). The functionality of PROFICS is demonstrated with viral proteases and human caspase-2. PROFICS improved caspase-2 activity up to 25-fold after only one round of mutation and selection. Additionally, we found a significantly improved tolerance for all P1' residues caused by a mutation in a substrate interaction site. We showed that this improved activity enables cells containing the new variant to outgrow cells containing all other mutants, facilitating its straightforward selection. Apart from optimizing enzymatic activity and P1' tolerance, PROFICS can be used to reprogram specificities, erase off-target activity, optimize expression via tags/codon usage, or even to screen for potential drug-resistance-conferring mutations in therapeutic targets such as viral proteases in an unbiased manner.

Efficient In Silico Saturation Mutagenesis of a Member of the Caspase Protease Family.

Rational-design methods have proven to be a valuable toolkit in the field of protein design. Numerical approaches such as free-energy calculations or QM/MM methods are fit to widen the understanding of a protein-sequence space but require large amounts of computational time and power. Here, we apply an efficient method for free-energy calculations that combines the one-step perturbation (OSP) with the third-power-fitting (TPF) approach. It is fit to calculate full free energies of binding from three different end states only. The nonpolar contribution to the free energies are calculated for a set of chosen amino acids from a single simulation of a judiciously chosen reference state. The electrostatic contributions, on the other hand, are predicted from simulations of the neutral and charged end states of the individual amino acids. We used this method to perform in silico saturation mutagenesis of two sites in human Caspase-2. We calculated relative binding free energies toward two different substrates that differ in their P1' site and in their affinity toward the unmutated protease. Although being approximate, our approach showed very good agreement upon validation against experimental data. 76% of the predicted relative free energies of amino acid mutations was found to be true positives or true negatives. We observed that this method is fit to discriminate amino acid mutations because the rate of false negatives is very low (<1.5%). The approach works better for a substrate with medium/low affinity with a Matthews correlation coefficient (MCC) of 0.63, whereas for a substrate with very low affinity, the MCC was 0.38. In all cases, the combined TPF + OSP approach outperformed the linear interaction energy method.

Enhancing the promiscuity of a member of the Caspase protease family by rational design.

The N-terminal cleavage of fusion tags to restore the native N-terminus of recombinant proteins is a challenging task and up to today, protocols need to be optimized for different proteins individually. Within this work, we present a novel protease that was designed in-silico to yield enhanced promiscuity toward different N-terminal amino acids. Two mutations in the active-site amino acids of human Caspase-2 were determined to increase the recognition of branched amino-acids, which show only poor binding capabilities in the unmutated protease. These mutations were determined by sequential and structural comparisons of Caspase-2 and Caspase-3 and their effect was additionally predicted using free-energy calculations. The two mutants proposed in the in-silico studies were expressed and in-vitro experiments confirmed the simulation results. Both mutants showed not only enhanced activities toward branched amino acids, but also smaller, unbranched amino acids. We believe that the created mutants constitute an important step toward generalized procedures to restore original N-termini of recombinant fusion proteins.

Study of key residues in MERS-CoV and SARS-CoV-2 main proteases for resistance against clinically applied inhibitors nirmatrelvir and ensitrelvir.

The Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is an epidemic, zoonotically emerging pathogen initially reported in Saudi Arabia in 2012. MERS-CoV has the potential to mutate or recombine with other coronaviruses, thus acquiring the ability to efficiently spread among humans and become pandemic. Its high mortality rate of up to 35% and the absence of effective targeted therapies call for the development of antiviral drugs for this pathogen. Since the beginning of the SARS-CoV-2 pandemic, extensive research has focused on identifying protease inhibitors for the treatment of SARS-CoV-2. Our intention was therefore to assess whether these protease inhibitors are viable options for combating MERS-CoV. To that end, we used previously established protease assays to quantify inhibition of SARS-CoV-2, MERS-CoV and other main proteases. Nirmatrelvir inhibited several of these proteases, whereas ensitrelvir was less broadly active. To simulate nirmatrelvir's clinical use against MERS-CoV and subsequent resistance development, we applied a safe, surrogate virus-based system. Using the surrogate virus, we previously selected hallmark mutations of SARS-CoV-2-Mpro, such as T21I, M49L, S144A, E166A/K/V and L167F. In the current study, we selected a pool of MERS-CoV-Mpro mutants, characterized the resistance and modelled the steric effect of catalytic site mutants S142G, S142R, S147Y and A171S.

Highly specific SARS-CoV-2 main protease (Mpro) mutations against the clinical antiviral ensitrelvir selected in a safe, VSV-based system.

In the SARS-CoV-2 pandemic, the so far two most effective approved antivirals are the protease inhibitors nirmatrelvir, in combination with ritonavir (Paxlovid) and ensitrelvir (Xocova). However, antivirals and indeed all antimicrobial drugs are sooner or later challenged by resistance mutations. Studying such mutations is essential for treatment decisions and pandemic preparedness. At the same time, generating resistant viruses to assess mutants is controversial, especially with pathogens of pandemic potential like SARS-CoV-2. To circumvent gain-of-function research with non-attenuated SARS-CoV-2, a previously developed safe system based on a chimeric vesicular stomatitis virus dependent on the SARS-CoV-2 main protease (VSV-Mpro) was used to select mutations against ensitrelvir. Ensitrelvir is clinically especially relevant due to its single-substance formulation, avoiding drug-drug interactions by the co-formulated CYP3A4 inhibitor ritonavir in Paxlovid. By treating VSV-Mpro with ensitrelvir, several highly-specific resistant mutants against this inhibitor were selected, while being still fully or largely susceptible to nirmatrelvir. We then confirmed several ensitrelvir-specific mutants in gold standard enzymatic assays and SARS-CoV-2 replicons. These findings indicate that the two inhibitors can have distinct viral resistance profiles, which could determine treatment decisions.

A comprehensive study of SARS-CoV-2 main protease (Mpro) inhibitor-resistant mutants selected in a VSV-based system.

Nirmatrelvir was the first protease inhibitor (PI) specifically developed against the SARS-CoV-2 main protease (3CLpro/Mpro) and licensed for clinical use. As SARS-CoV-2 continues to spread, variants resistant to nirmatrelvir and other currently available treatments are likely to arise. This study aimed to identify and characterize mutations that confer resistance to nirmatrelvir. To safely generate Mpro resistance mutations, we passaged a previously developed, chimeric vesicular stomatitis virus (VSV-Mpro) with increasing, yet suboptimal concentrations of nirmatrelvir. Using Wuhan-1 and Omicron Mpro variants, we selected a large set of mutants. Some mutations are frequently present in GISAID, suggesting their relevance in SARS-CoV-2. The resistance phenotype of a subset of mutations was characterized against clinically available PIs (nirmatrelvir and ensitrelvir) with cell-based and biochemical assays. Moreover, we showed the putative molecular mechanism of resistance based on in silico molecular modelling. These findings have implications on the development of future generation Mpro inhibitors, will help to understand SARS-CoV-2 protease-inhibitor-resistance mechanisms and show the relevance of specific mutations in the clinic, thereby informing treatment decisions.

CASPON platform technology: Ultrafast circularly permuted caspase-2 cleaves tagged fusion proteins before all 20 natural amino acids at the N-terminus.

Fusion protein technologies improve the expression and purification of recombinant proteins, but the removal of the tags involved requires specific proteases. The circularly permuted caspase-2 (cpCasp2) with its specific cleavage site, efficiently generates the untagged protein. While cleavage with cpCasp2 is possible before all 20 proteinogenic amino acids, cleavage before valine, leucine, isoleucine, aspartate and glutamate suffers from slow, and before proline extremely slow, turnover. To make the platform fusion protein process even more general such that any protein with an authentic N-terminus can be produced with high efficiency, the bacterial selection system PROFICS (PRotease Optimization via Fusion-Inhibited Carbamoyltransferase-based Selection) was used to evolve cpCasp2 into a variant with a catalytic turnover two orders of magnitude higher and the ability to cleave before any amino acid. The high specificity and the stability of the original circularly permuted protease was fully retained in this mutant, while the high manufacturability was mostly retained, albeit with decreased soluble titer. Four point-mutations are responsible for this change in activity, two of which are located in or near the binding pocket of the active site. This variant was named CASPON enzyme and is a major component of the CASPase-based fusiON (CASPON) platform technology. Applicability for the production of recombinant proteins was demonstrated by enzymatic removal of the CASPON tag from five model proteins. The CASPON tag enables high soluble expressions, affinity purification and good accessibility for cleavage. The five industry-relevant proteins of interest were FGF2, TNF, GH, GCSF and PTH.