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.

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.

Sub-kHz lasing of a CaF₂ whispering gallery mode resonator stabilized fiber ring laser

We utilize a high quality calcium fluoride whispering-gallery-mode resonator to passively stabilize a simple erbium doped fiber ring laser with an emission frequency of 196THz (wavelength 1530nm) to an instantaneous linewidth below 650Hz. This corresponds to a relative stability of 3.3 × 10(-12) over 16μs. In order to characterize the linewidth we use two identical self-built lasers and a commercial laser to determine the individual lasing linewidth via the three-cornered-hat method. We further show that the lasers are finely tunable throughout the erbium gain region.

CaF2 whispering-gallery-mode-resonator stabilized-narrow-linewidth laser

A fiber laser is stabilized by introducing a calcium fluoride (CaF(2)) whispering-gallery-mode resonator as a filtering element in a ring cavity. It is set up using a semiconductor optical amplifier as a gain medium. The resonator is critically coupled through prisms, and used as a filtering element to suppress the laser linewidth. A three-cornered-hat method is used and shows a stability of 10(-11) after 10 micros. Using the self-heterodyne beat technique, the linewidth is determined to be 13 kHz. This implies an enhancement factor of 10(3) with respect to the passive cavity linewidth.

Whispering-gallery-mode-resonator-stabilized narrow-linewidth fiber loop laser

We demonstrate a narrow-line fiber loop laser using erbium-doped fiber as the gain material, stabilized by using a microsphere as a transmissive frequency selective element. Stable lasing with a linewidth of 170 kHz is observed, limited by the experimental spectral resolution. A linear increase in output power and a redshift of the lasing mode were also observed with increasing pump power. Its potential applications are discussed.

Atmospheric transfer of optical and radio frequency clock signals

The phase instability induced during the transfer of radio frequency and optical clock signals through the turbulent atmosphere was measured in a rooftop experiment. Radio frequency intensity modulation of a laser to transmit signals over 100 m results in an Allan deviation of 1.31x10(-10) at 1 s. Optical transfer is more accurate at 1.68x10(-13) at 1 s. As a consequence, fiber links are more suitable for the transfer of optical frequencies over very long distances while free space transmission might find applications in short distances of less than 1 km.