High-Resolution Substrate Specificity Profiling of SARS-CoV-2 Mpro; Comparison to SARS-CoV Mpro

The SARS-CoV-2 Mpro protease from COVID-19 cleaves the pp1a and pp2b polyproteins at 11 sites during viral maturation and is the target of Nirmatrelvir, one of the two components of the frontline treatment sold as Paxlovid. We used the YESS 2.0 platform, combining protease and substrate expression in the yeast endoplasmic reticulum with fluorescence-activated cell sorting and next-generation sequencing, to carry out the high-resolution substrate specificity profiling of SARS-CoV-2 Mpro as well as the related SARS-CoV Mpro from SARS 2003. Even at such a high level of resolution, the substrate specificity profiles of both enzymes are essentially identical. The population of cleaved substrates isolated in our sorts is so deep, the relative catalytic efficiencies of the different cleavage sites on the SARS-CoV-2 polyproteins pp1a and pp2b are qualitatively predicted. These results not only demonstrated the precise and reproducible nature of the YESS 2.0/NGS approach to protease substrate specificity profiling but also should be useful in the design of next generation SARS-CoV-2 Mpro inhibitors, and by analogy, SARS-CoV Mpro inhibitors as well.

Chemical and linguistic considerations for encoding Chinese characters: an embodiment using chain-end degradable sequence-defined oligourethanes created by consecutive solid phase click chemistry

Sequence-defined polymers (SDPs) are currently being investigated for use as information storage media. As the number of monomers in the SDPs increases, with a corresponding increase in mathematical base, the use of tandem-MS for de novo sequencing becomes more challenging. In contrast, chain-end degradation routines are truly de novo, potentially allowing very large mathematical bases for encoding. While alphabetic scripts have a few dozen symbols, logographic scripts, such as Chinese, can have several thousand symbols. Using a new in situ consecutive click reaction approach on an oligourethane backbone for writing, and a previously reported chain-end degradation routine for reading, we encoded/decoded a confucius proverb written in Chinese characters using two encoding schemes: Unicode and Zhèng Mă. Unicode is an internationally standardized arbitrary string of hexadecimal (base-16) symbols which efficiently encodes uniquely identifiable symbols but requires complete fidelity of transmission, or context-based inferential strategies to be interpreted. The Zhèng Mă approach encodes with a base-26 system using the visual characteristics and internal composition of Chinese characters themselves, which leads to greater ambiguity of encoded strings, but more robust retrievability of information from partial or corrupted encodings. The application of information-encoded oligourethanes to two different encoding systems allowed us to establish their flexibility and versatility for data storage. We found the oligourethanes immensely adaptable to both encoding schemes for Chinese characters, and we highlight the expected tradeoff between the efficiency and uniqueness of Unicode encoding on the one hand, and the fidelity to a scripts' particular visual characteristics on the other.

Aggregation of Charge Acceptors on Nanocrystal Surfaces Alters Rates of Photoinduced Electron Transfer

Semiconductor nanocrystals (NCs) interfaced with molecular ligands that function as charge and energy acceptors are an emerging platform for the design of light-harvesting, photon-upconverting, and photocatalytic materials. However, NC systems explored for these applications often feature high concentrations of bound acceptor ligands, which can lead to ligand-ligand interactions that may alter each system's ability to undergo charge and energy transfer. Here, we demonstrate that aggregation of acceptor ligands impacts the rate of photoinduced NC-to-ligand charge transfer between lead(II) sulfide (PbS) NCs and perylenediimide (PDI) electron acceptors. As the concentration of PDI acceptors is increased, we find the average electron transfer rate from PbS to PDI ligands decreases by nearly an order of magnitude. The electron transfer rate slowdown with increasing PDI concentration correlates strongly with the appearance of PDI aggregates in steady-state absorption spectra. Electronic structure calculations and molecular dynamics (MD) simulations suggest PDI aggregation slows the rate of electron transfer by reducing orbital overlap between PbS charge donors and PDI charge acceptors. While we find aggregation slows electron transfer in this system, the computational models we employ predict ligand aggregation could also be used to speed electron transfer by producing delocalized states that exhibit improved NC-molecule electronic coupling and energy alignment with NC conduction band states. Our results demonstrate that ligand aggregation can alter rates of photoinduced electron transfer between NCs and organic acceptor ligands and should be considered when designing hybrid NC:molecule systems for charge separation.

Evaluating the Effect of Dye-Dye Interactions of Xanthene-Based Fluorophores in the Fluorosequencing of Peptides

A peptide sequencing scheme utilizing fluorescence microscopy and Edman degradation to determine the amino acid position in fluorophore-labeled peptides was recently reported, referred to as fluorosequencing. It was observed that multiple fluorophores covalently linked to a peptide scaffold resulted in a decrease in the anticipated fluorescence output and worsened the single-molecule fluorescence analysis. In this study, we report an improvement in the photophysical properties of fluorophore-labeled peptides by incorporating long and flexible (PEG)₁₀ linkers at the peptide attachment points. Long linkers to the fluorophores were installed using copper-catalyzed azide-alkyne cycloaddition conditions. The photophysical properties of these peptides were analyzed in solution and immobilized on a microscope slide at the single-molecule level under peptide fluorosequencing conditions. Solution-phase fluorescence analysis showed improvements in both quantum yield and fluorescence lifetime with the long linkers. While on the solid support, photometry measurements showed significant increases in fluorescence brightness and 20 to 60% improvements in the ability to determine the amino acid position with fluorosequencing. This spatial distancing strategy demonstrates improvements in the peptide sequencing platform and provides a general approach for improving the photophysical properties in fluorophore-labeled macromolecules.

Using Spectator Ligands to Enhance Nanocrystal-to-Molecule Electron Transfer

Semiconductor nanocrystals (NCs) have emerged as promising photocatalysts. However, NCs are often functionalized with complex ligand shells that contain not only charge acceptors but also other "spectator ligands" that control NC solubility and affinity for target reactants. Here, we show that spectator ligands are not passive observers of photoinduced charge transfer but rather play an active role in this process. We find the rate of electron transfer from quantum-confined PbS NCs to perylenediimide acceptors can be varied by over a factor of 4 simply by coordinating cinnamate ligands with distinct dipole moments to NC surfaces. Theoretical calculations indicate this rate variation stems from both ligand-induced changes in the free energy for charge transfer and electrostatic interactions that alter perylenediimide electron acceptor orientation on NC surfaces. Our work shows NC-to-molecule charge transfer can be fine-tuned through ligand shell design, giving researchers an additional handle for enhancing NC photocatalysis.

Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes

The molecular composition and binding epitopes of the immunoglobulin G (IgG) antibodies that circulate in blood plasma after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are unknown. Proteomic deconvolution of the IgG repertoire to the spike glycoprotein in convalescent subjects revealed that the response is directed predominantly (>80%) against epitopes residing outside the receptor binding domain (RBD). In one subject, just four IgG lineages accounted for 93.5% of the response, including an amino (N)-terminal domain (NTD)-directed antibody that was protective against lethal viral challenge. Genetic, structural, and functional characterization of a multidonor class of "public" antibodies revealed an NTD epitope that is recurrently mutated among emerging SARS-CoV-2 variants of concern. These data show that "public" NTD-directed and other non-RBD plasma antibodies are prevalent and have implications for SARS-CoV-2 protection and antibody escape.

YESS 2.0, a Tunable Platform for Enzyme Evolution, Yields Highly Active TEV Protease Variants

Here we describe YESS 2.0, a highly versatile version of the yeast endoplasmic sequestration screening (YESS) system suitable for engineering and characterizing protein/peptide modifying enzymes such as proteases with desired new activities. By incorporating features that modulate gene transcription as well as substrate and enzyme spatial sequestration, YESS 2.0 achieves a significantly higher operational and dynamic range compared with the original YESS. To showcase the new advantages of YESS 2.0, we improved an already efficient TEV protease variant (TEV-EAV) to obtain a variant (eTEV) with a 2.25-fold higher catalytic efficiency, derived almost entirely from an increase in turnover rate (k_cat). In our analysis, eTEV specifically digests a fusion protein in 2 h at a low 1:200 enzyme to substrate ratio. Structural modeling indicates that the increase in catalytic efficiency of eTEV is likely due to an enhanced interaction between the catalytic Cys151 with the P1 substrate residue (Gln). Furthermore, the modeling showed that the ENLYFQS peptide substrate is buried to a larger extent in the active site of eTEV compared with WT TEV. The new eTEV variant is functionally the fastest TEV variant reported to date and could potentially improve efficiency in any TEV application.

Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes in COVID-19 convalescent plasma

Although humoral immunity is essential for control of SARS-CoV-2, the molecular composition, binding epitopes and effector functions of the immunoglobulin G (IgG) antibodies that circulate in blood plasma following infection are unknown. Proteomic deconvolution of the circulating IgG repertoire (Ig-Seq ¹ ) to the spike ectodomain (S-ECD ² ) in four convalescent study subjects revealed that the plasma response is oligoclonal and directed predominantly (>80%) to S-ECD epitopes that lie outside the receptor binding domain (RBD). When comparing antibodies directed to either the RBD, the N-terminal domain (NTD) or the S2 subunit (S2) in one subject, just four IgG lineages (1 anti-S2, 2 anti-NTD and 1 anti-RBD) accounted for 93.5% of the repertoire. Although the anti-RBD and one of the anti-NTD antibodies were equally potently neutralizing in vitro , we nonetheless found that the anti-NTD antibody was sufficient for protection to lethal viral challenge, either alone or in combination as a cocktail where it dominated the effect of the other plasma antibodies. We identified in vivo protective plasma anti-NTD antibodies in 3/4 subjects analyzed and discovered a shared class of antibodies targeting the NTD that utilize unmutated or near-germline IGHV1-24, the most electronegative IGHV gene in the human genome. Structural analysis revealed that binding to NTD is dominated by interactions with the heavy chain, accounting for 89% of the entire interfacial area, with germline residues uniquely encoded by IGHV1-24 contributing 20% (149 Å ² ). Together with recent reports of germline IGHV1-24 antibodies isolated by B-cell cloning ^3,4 our data reveal a class of shared IgG antibodies that are readily observed in convalescent plasma and underscore the role of NTD-directed antibodies in protection against SARS-CoV-2 infection.

Quantitative Analysis of the Substrate Specificity of Human Rhinovirus 3C Protease and Exploration of Its Substrate Recognition Mechanisms

Human rhinovirus 3C protease (HRV 3C-P) is a high-value commercial cysteine protease that could specifically recognize the short peptide sequence of LEVLFQ↓GP. In here, a strategy based on our previous Yeast Endoplasmic Reticulum Sequestration Screening (YESS) approach was developed in Saccharomyces cerevisiae, a model microorganism, to fully characterize the substrate specificity of a typical human virus protease, HRV 3C-P, in a quantitative and fast manner. Our results demonstrated that HRV 3C-P had very high specificity at P1 and P1' positions, only recognizing Gln/Glu at the P1 position and Gly/Ala/Cys/Ser at the P1' position, respectively. Comparably, it exhibited efficient recognition of most residues at the P2' position, except Trp. Further biochemical characterization through site mutagenesis, enzyme structural modeling, and comparison with other 3C proteases indicated that the S1 pocket of HRV 3C-P was constituted by neutral and basic amino acids, in which His160 and Thr141 specifically interacted with Gln or Glu residues at the substrate P1 position. Additionally, the stringent S1' pocket determined its unique property of only accommodating residues without or with short side chains. Based on our characterization, LEVLFQ↓GM was identified as a more favorable substrate than the original LEVLFQ↓GP at high temperature, which might be caused by the conversion of random coils to β-turns in HRV 3C-P along with the temperature increase. Our studies prompted a further understanding of the substrate specificity and recognition mechanism of HRV 3C-P. Besides, the YESS-PSSC combined with the enzyme modeling strategy in this study provides a general strategy for deciphering the substrate specificities of proteases.

Prompting Fab Yeast Surface Display Efficiency by ER Retention and Molecular Chaperon Co-expression

For antibody discovery and engineering, yeast surface display (YSD) of antigen-binding fragments (Fabs) and coupled fluorescence activated cell sorting (FACS) provide intact paratopic conformations and quantitative analysis at the monoclonal level, and thus holding great promises for numerous applications. Using anti-TNFα mAbs Infliximab, Adalimumab, and its variants as model Fabs, this study systematically characterized complementary approaches for the optimization of Fab YSD. Results suggested that by using divergent promoter GAL1-GAL10 and endoplasmic reticulum (ER) signal peptides for co-expression of light chain and heavy chain-Aga2 fusion, assembled Fabs were functionally displayed on yeast cell surface with sigmoidal binding responses toward TNFα. Co-expression of a Hsp70 family molecular chaperone Kar2p and/or protein-disulfide isomerase (Pdi1p) significantly improved efficiency of functional display (defined as the ratio of cells displaying functional Fab over cells displaying assembled Fab). Moreover, fusing ER retention sequences (ERSs) with light chain also enhanced Fab display quality at the expense of display quantity, and the degree of improvements was correlated with the strength of ERSs and was more significant for Infliximab than Adalimumab. The feasibility of affinity maturation was further demonstrated by isolating a high affinity Fab clone from 1:10³ or 1:10⁵ spiked libraries.

Rapid Screen for Tyrosine Kinase Inhibitor Resistance Mutations and Substrate Specificity

We present a rapid and high-throughput yeast and flow cytometry based method for predicting kinase inhibitor resistance mutations and determining kinase peptide substrate specificity. Despite the widespread success of targeted kinase inhibitors as cancer therapeutics, resistance mutations arising within the kinase domain of an oncogenic target present a major impediment to sustained treatment efficacy. Our method, which is based on the previously reported YESS system, recapitulated all validated BCR-ABL1 mutations leading to clinical resistance to the second-generation inhibitor dasatinib, in addition to identifying numerous other mutations which have been previously observed in patients, but not yet validated as drivers of resistance. Further, we were able to demonstrate that the newer inhibitor ponatinib is effective against the majority of known single resistance mutations, but ineffective at inhibiting many compound mutants. These results are consistent with preliminary clinical and in vitro reports, indicating that mutations providing resistance to ponatinib are significantly less common; therefore, predicting ponatinib will be less susceptible to clinical resistance relative to dasatinib. Using the same yeast-based method, but with random substrate libraries, we were able to identify consensus peptide substrate preferences for the SRC and LYN kinases. ABL1 lacked an obvious consensus sequence, so a machine learning algorithm utilizing amino acid covariances was developed which accurately predicts ABL1 kinase peptide substrates.

Exploiting the interactions of aromatic units for folding and assembly in aqueous environments

A variety of non-covalent interactions (including hydrogen bonding, ionic interactions, metal coordination and desolvation/solvation) have been utilized to organize oligomers into well-defined structures. Herein is described a survey of aromatic foldamers that capitalize on electrostatic complementarity of substituted aromatic units to drive folding and assembly in aqueous environments. A brief description of recent advances in the understanding of aromatic interactions is provided, followed by examples of foldamers that exploit interactions between aromatic units to drive their assembly in predictable fashion. The history of our aromatic foldamers is traced from the first structure designed to fold into a pleated structure in an aqueous environment to a heteroduplex system more related to nucleic acids. Taken together, the results demonstrate that electrostatic complementarity of aromatic units provides a versatile framework for driving predictable folding and assembly in aqueous environments.

Characterization of aromatic residue-controlled protein retention in the endoplasmic reticulum of Saccharomyces cerevisiae

An endoplasmic reticulum (ER) retention sequence (ERS) is a characteristic short sequence that mediates protein retention in the ER of eukaryotic cells. However, little is known about the detailed molecular mechanism involved in ERS-mediated protein ER retention. Using a new surface display-based fluorescence technique that effectively quantifies ERS-promoted protein ER retention within Saccharomyces cerevisiae cells, we performed comprehensive ERS analyses. We found that the length, type of amino acid residue, and additional residues at positions -5 and -6 of the C-terminal HDEL motif all determined the retention of ERS in the yeast ER. Moreover, the biochemical results guided by structure simulation revealed that aromatic residues (Phe-54, Trp-56, and other aromatic residues facing the ER lumen) in both the ERS (at positions -6 and -4) and its receptor, Erd2, jointly determined their interaction with each other. Our studies also revealed that this aromatic residue interaction might lead to the discriminative recognition of HDEL or KDEL as ERS in yeast or human cells, respectively. Our findings expand the understanding of ERS-mediated residence of proteins in the ER and may guide future research into protein folding, modification, and translocation affected by ER retention.

Profiling Protease Specificity: Combining Yeast ER Sequestration Screening (YESS) with Next Generation Sequencing

An enzyme engineering technology involving yeast endoplasmic reticulum (ER) sequestration screening (YESS) has been recently developed. Here, a new method is established, in which the YESS platform is combined with NextGen sequencing (NGS) to enable a comprehensive survey of protease specificity. In this approach, a combinatorial substrate library is targeted to the yeast ER and transported through the secretory pathway, interacting with any protease(s) residing in the ER. Multicolor FACS screening is used to isolate cells labeled with fluorophore-conjugated antibodies, followed by NGS to profile the cleaved substrates. The YESS-NGS method was successfully applied to profile the sequence specificity of the wild-type and an engineered variant of the tobacco etch mosaic virus protease. Proteolysis in the yeast secretory pathway was also mapped for the first time in vivo revealing a major cleavage pattern of Ali/Leu-X-Lys/Arg-Arg. Here Ali is any small aliphatic residue, but especially Leu. This pattern was verified to be due to the well-known endogenous protease Kex2 after comparison to a newly generated Kex2 knockout strain as well as cleavage of peptides with recombinant Kex2 in vitro. This information is particularly important for those using yeast display technology, as library members with Ali/Leu-X-Lys/Arg-Arg patterns are likely being removed from screens via Kex2 cleavage without the researcher's knowledge.

Yeast Endoplasmic Reticulum Sequestration Screening for the Engineering of Proteases from Libraries Expressed in Yeast

There is significant interest in engineering proteases with desired proteolytic properties. We describe a high-throughput fluorescence-activated cell sorting (FACS) assay for detecting altered proteolytic activity of protease in yeast, at the single cell level. This assay relies on coupling yeast endoplasmic reticulum (ER) retention, yeast surface display, and FACS analysis. The method described here allows facile screening of large libraries, and of either protease or substrate variants, including the screening of protease libraries against substrate libraries. We demonstrate the application of this technique in the screening of libraries of Tobacco Etch Virus protease (TEV-P) for altered proteolytic activities. In addition, the generality of this method is also validated by other proteases such as human granzyme K and the hepatitis C virus protease, and the human Abelson tyrosine kinase.

NDI and DAN DNA: nucleic acid-directed assembly of NDI and DAN

Two novel DNA base surrogate phosphoramidites 1 and 2, based upon relatively electron-rich 1,5-dialkoxynaphthalene (DAN) and relatively electron-deficient 1,4,5,8-naphthalenetetracarboxylic diimide (NDI), respectively, were designed, synthesized, and incorporated into DNA oligonucleotide strands. The DAN and NDI artificial DNA bases were inserted within a three-base-pair region within the interior of a 12-mer oligonucleotide duplex in various sequential arrangements and investigated with CD spectroscopy and UV melting curve analysis. The CD spectra of the modified duplexes indicated B-form DNA topology. Melting curve analyses revealed trends in DNA duplex stability that correlate with the known association of DAN and NDI moieties in aqueous solution as well as the known favorable interactions between NDI and natural DNA base pairs. This demonstrates that DNA duplex stability and specificity can be driven by the electrostatic complementarity between DAN and NDI. In the most favorable case, an NDI-DAN-NDI arrangement in the middle of the DNA duplex was found to be approximately as stabilizing as three A-T base pairs.

Time-Dependent Solid State Polymorphism of a Series of Donor-Acceptor Dyads

In order to exploit the use of favorable electrostatic interactions between aromatic units in directing the assembly of donor-acceptor (D-A) dyads, the present work examines the ability of conjugated aromatic D-A dyads with symmetric side chains to exhibit solid-state polymorphism as a function of time during the solid formation process. Four such dyads were synthesized and their packing in the solid-state from either slower (10-20 days) or faster (1-2 days) evaporation from solvent was investigated using single crystal X-ray analysis and powder X-ray diffraction. Two of the dyads exhibited tail-to-tail (A-A) packing upon slower evaporation from solvent and head-to-tail (D-A) packing upon faster evaporation from solvent. A combination of single crystal analysis and XRD patterns were used to create models wherein a packing model for the other two dyads is proposed. Our findings suggest that while side chain interactions in asymmetric aromatic dyads can play an important role in enforcing segregated D-A dyad assembly, slowly evaporating symmetrically substituted aromatic dyads allows for favorable electrostatic interactions between the aromatic moieties to facilitate the organization of the dyads in the solid-state.

Threading polyintercalators with extremely slow dissociation rates and extended DNA binding sites

The development of small molecules that bind DNA sequence specifically has the potential to modulate gene expression in a general way. One mode of DNA binding is intercalation, or the insertion of molecules between DNA base pairs. We have developed a modular polyintercalation system in which intercalating naphthalene diimide (NDI) units are connected by flexible linkers that alternate between the minor and major grooves of DNA when bound. We recently reported a threading tetraintercalator with a dissociation half-life of 16 days, the longest reported to date, from its preferred 14 bp binding site. Herein, three new tetraintercalator derivatives were synthesized with one, two, and three additional methylene units in the central major groove-binding linker. These molecules displayed dissociation half-lives of 57, 27, and 18 days, respectively, from the 14 bp site. The optimal major groove-binding linker was used in the design of an NDI hexaintercalator that was analyzed by gel-shift assays, DNase I footprinting, and UV-vis spectroscopy. The hexaintercalator bound its entire 22 bp binding site, the longest reported specific binding site for a synthetic, non-nucleic acid-based DNA binding molecule, but with a significantly faster dissociation rate compared to the tetraintercalators.

More than meets the eye: conformational switching of a stacked dialkoxynaphthalene-naphthalenetetracarboxylic diimide (DAN-NDI) foldamer to an NDI-NDI fibril aggregate

The thermally induced conformational switching of a stacked dialkxoynaphthalene-naphthalenetetracarboxylic diimide (DAN-NDI) amphiphilic foldamer to an NDI-NDI fibril aggregate is described. The aggregated fibril structures were explored by UV/Vis, circular dichroism (CD), atomic-force microscopy (AFM), and TEM techniques. Our findings indicate that the aromatic DAN-NDI interactions of the original foldamer undergoes transformation to a fibrillar assembly with aromatic NDI-NDI stacked interactions. These structural insights could help inform new molecular designs and increase our understanding of fibrillar assembly and aggregation process in aqueous solution.

Commercial proteases: present and future

This review presents a brief overview of the general categories of commercially used proteases, and critically surveys the successful strategies currently being used to improve the properties of proteases for various commercial purposes. We describe the broad application of proteases in laundry detergents, food processing, and the leather industry. The review also introduces the expanding development of proteases as a class of therapeutic agents, as well as highlighting recent progress in the field of protease engineering. The potential commercial applications of proteases are rapidly growing as recent technological advances are producing proteases with novel properties and substrate specificities.

Subtle recognition of 14-base pair DNA sequences via threading polyintercalation

Small molecules that bind DNA in a sequence-specific manner could act as antibiotic, antiviral, or anticancer agents because of their potential ability to manipulate gene expression. Our laboratory has developed threading polyintercalators based on 1,4,5,8-naphthalene diimide (NDI) units connected in a head-to-tail fashion by flexible peptide linkers. Previously, a threading tetraintercalator composed of alternating minor-major-minor groove-binding modules was shown to bind specifically to a 14 bp DNA sequence with a dissociation half-life of 16 days [Holman, G. G., et al. (2011) Nat. Chem. 3, 875-881]. Herein are described new NDI-based tetraintercalators with a different major groove-binding module and a reversed N to C directionality of one of the minor groove-binding modules. DNase I footprinting and kinetic analyses revealed that these new tetraintercalators are able to discriminate, by as much as 30-fold, 14 bp DNA binding sites that differ by 1 or 2 bp. Relative affinities were found to correlate strongly with dissociation rates, while overall C(2) symmetry in the DNA-binding molecule appeared to contribute to enhanced association rates.