In Silico Screening of Inhibitors of the Venezuelan Equine Encephalitis Virus Nonstructural Protein 2 Cysteine Protease.

The Venezuelan equine encephalitis virus (VEEV) nonstructural protein 2 (nsP2) cysteine protease (EC 3.4.22.B79) is essential for viral replication. High throughput in silico/in vitro screening using a focused set of known cysteine protease inhibitors identified two epoxysuccinyl prodrugs, E64d and CA074 methyl ester (CA074me) and a reversible oxindole inhibitor. Here, we determined the X-ray crystal structure of the CA074-inhibited nsP2 protease and compared it with our E64d-inhibited structure. We found that the two inhibitors occupy different locations in the protease. We designed hybrid inhibitors with improved potency. Virus yield reduction assays confirmed that the viral titer was reduced by >5 logs with CA074me. Cell-based assays showed reductions in viral replication for CHIKV, VEEV, and WEEV, and weaker inhibition of EEEV by the hybrid inhibitors. The most potent was NCGC00488909-01 which had an EC50 of 1.76 µM in VEEV-Trd-infected cells; the second most potent was NCGC00484087 with an EC50 = 7.90 µM. Other compounds from the NCATS libraries such as the H1 antihistamine oxatomide (>5-log reduction), emetine, amsacrine an intercalator (NCGC0015113), MLS003116111-01, NCGC00247785-13, and MLS00699295-01 were found to effectively reduce VEEV viral replication in plaque assays. Kinetic methods demonstrated time-dependent inhibition by the hybrid inhibitors of the protease with NCGC00488909-01 (Ki = 3 µM) and NCGC00484087 (Ki = 5 µM). Rates of inactivation by CA074 in the presence of 6 mM CaCl2, MnCl2, or MgCl2 were measured with varying concentrations of inhibitor, Mg2+ and Mn2+ slightly enhanced inhibitor binding (3 to 6-fold). CA074 inhibited not only the VEEV nsP2 protease but also that of CHIKV and WEEV.

Automated SSHHPS Analysis Predicts a Potential Host Protein Target Common to Several Neuroinvasive (+)ssRNA Viruses.

Within the viral genome, short stretches of homologous host pathogen sequences (SSHHPS) span the protease cleavage sites. To identify host proteins that may be cleaved during infection, we searched the human proteome for viral protease cleavage sites (~20 amino acids). We developed a sequence-to-symptom tool, automating the search and pairing process. We used the viral protein sequence, PHI-BLAST, and UniProt database for gene ontologies and disease relationships. We applied the tool to nine neuroinvasive viruses: Venezuelan and Eastern Equine encephalitis virus (VEEV, EEEV); severe acute respiratory syndrome (SARS, SARS-CoV-2); Middle East respiratory syndrome (MERS); EV-71; Japanese encephalitis virus (JEV); West Nile (WNV); and Zika (ZIKV). A comparison of the hits identified a protein common to all nine viruses called ADGRA2 (GPR124). ADGRA2 was a predicted hit of the 3CL main protease and papain-like protease (PLpro) of SARS-CoV-2. ADGRA2 is an adhesion G protein-coupled receptor and a key endothelial regulator of brain-specific angiogenesis. It is a Wnt7A/Wnt7B specific coactivator of beta-catenin signaling and is essential for blood-brain barrier (BBB) integrity in central nervous system (CNS) diseases. We show the cleavage of the predicted sequences in MYOM1, VWF by the SARS-CoV-2 PLpro; DNAH8 (dynein) by the MERS PLpro; ADGRA2 by the alphaviral VEEV nsP2 protease; and POT1 by the SARS-CoV-2 and MERS PLpro.

Engineering an Fc-Fusion of a Capsule Degrading Enzyme for the Treatment of Anthrax.

Polymers of d-glutamic acid (PDGA) form the capsule of the highly virulent Ames strain of B. anthracis. PDGA is antiphagocytic and weakly immunogenic; it enables the bacteria to evade the innate immune responses. CapD is an enzyme that catalyzes the covalent anchoring of PDGA. CapD is an Ntn-amido hydrolase that utilizes an internal Thr-352 as its nucleophile and general acid and base. An internal cleavage produces a free N-terminal Thr-352 and a short and long polypeptide chain. The chains were circularly permuted (CP) to move Thr-352 to the N-terminus of the polypeptide. We previously showed that a branched PEG-CapDS334C-CP could protect mice (80% survival) against a 5 LD50 challenge with B. anthracis Ames without the use of antibiotics, monoclonals, or vaccines. In attempts to improve the in vivo circulation time of CapD and enhance its avidity to its polymeric substrate, an Fc-domain of a mouse IgG1 was fused to CapDS334C-CP and the linker length and sequence were optimized. The resulting construct, Fc-CapDS334C-CP, then was pegylated with a linear 2 kDa mPEG at S334C to produce mPEG-Fc-CapDS334C-CP. Interestingly, the fusion of the Fc-domain and incorporation of the S334C mutation imparted acid stability, but slightly reduced the kcat (∼ 2-fold lower). In vivo, the measured protein concentration in sera was higher for the Fc-fusion constructs compared to the mPEG-Fc-CapDS334C-CP. However, the exposure calculated from measured sera enzymatic activity was higher for the mPEG-CapDS334C-CP. The pegylated Fc-fusion was less active than the PEG-CapDS334C-CP, but detectable in sera at 24 h by immunoblot. Here we describe the engineering of a soluble, active, pegylated Fc-fusion of B. anthracis CapD (mPEG-Fc-CapD-CP) with activity in vitro, in serum, and on encapsulated bacteria.

Treatment of experimental anthrax with pegylated circularly permuted capsule depolymerase.

Anthrax is considered one of the most dangerous bioweapon agents, and concern about multidrug-resistant strains has led to the development of alternative therapeutic approaches that target the antiphagocytic capsule, an essential virulence determinant of Bacillus anthracis, the causative agent. Capsule depolymerase is a γ-glutamyltransferase that anchors the capsule to the cell wall of B. anthracis. Encapsulated strains of B. anthracis can be treated with recombinant capsule depolymerase to enzymatically remove the capsule and promote phagocytosis and killing by human neutrophils. Here, we show that pegylation improved the pharmacokinetic and therapeutic properties of a previously described variant of capsule depolymerase, CapD-CP, when delivered 24 hours after exposure every 8 hours for 2 days for the treatment of mice infected with B. anthracis. Mice infected with 382 LD50 of B. anthracis spores from a nontoxigenic encapsulated strain were completely protected (10 of 10) after treatment with the pegylated PEG-CapD-CPS334C, whereas 10% of control mice (1 of 10) survived with control treatment using bovine serum albumin (P < 0.0001, log-rank analysis). Treatment of mice infected with five LD50 of a fully virulent toxigenic, encapsulated B. anthracis strain with PEG-CapD-CPS334C protected 80% (8 of 10) of the animals, whereas 20% of controls (2 of 10) survived (P = 0.0125, log-rank analysis). This strategy renders B. anthracis susceptible to innate immune responses and does not rely on antibiotics. These findings suggest that enzyme-catalyzed removal of the capsule may be a potential therapeutic strategy for the treatment of multidrug- or vaccine-resistant anthrax and other bacterial infections.

The SARS-CoV-2 SSHHPS Recognized by the Papain-like Protease.

Viral proteases are highly specific and recognize conserved cleavage site sequences of ∼6-8 amino acids. Short stretches of homologous host-pathogen sequences (SSHHPS) can be found spanning the viral protease cleavage sites. We hypothesized that these sequences corresponded to specific host protein targets since >40 host proteins have been shown to be cleaved by Group IV viral proteases and one Group VI viral protease. Using PHI-BLAST and the viral protease cleavage site sequences, we searched the human proteome for host targets and analyzed the hit results. Although the polyprotein and host proteins related to the suppression of the innate immune responses may be the primary targets of these viral proteases, we identified other cleavable host proteins. These proteins appear to be related to the virus-induced phenotype associated with Group IV viruses, suggesting that information about viral pathogenesis may be extractable directly from the viral genome sequence. Here we identify sequences cleaved by the SARS-CoV-2 papain-like protease (PLpro) in vitro within human MYH7 and MYH6 (two cardiac myosins linked to several cardiomyopathies), FOXP3 (an X-linked Treg cell transcription factor), ErbB4 (HER4), and vitamin-K-dependent plasma protein S (PROS1), an anticoagulation protein that prevents blood clots. Zinc inhibited the cleavage of these host sequences in vitro. Other patterns emerged from multispecies sequence alignments of the cleavage sites, which may have implications for the selection of animal models and zoonosis. SSHHPS/nsP is an example of a sequence-specific post-translational silencing mechanism.

Stabilization of a Broadly Neutralizing Anti-Chikungunya Virus Single Domain Antibody.

A single domain antibody (clone CC3) previously found to neutralize a vaccine strain of the chikungunya virus (PRNT50 = 2. 5 ng/mL) was found to be broadly neutralizing. Clone CC3 is not only able to neutralize a wild-type (WT) strain of chikungunya virus (CHIKV), but also neutralizes WT strains of Mayaro virus (MAYV) and Ross River virus (RRV); both arthralgic, Old World alphaviruses. Interestingly, CC3 also demonstrated a degree of neutralizing activity against the New World alphavirus, Venezuelan equine encephalitis virus (VEEV); albeit both the vaccine strain, TC-83, and the parental, WT Trinidad donkey strain had PRNT50 values ~1,000-fold higher than that of CHIKV. However, no neutralization activity was observed with Western equine encephalitis virus (WEEV). Ten CC3 variants designed to possess a range of isoelectric points, both higher and lower, were constructed. This approach successfully identified several lower pI mutants which possessed improved thermal stabilities by as much as 10°C over the original CC3 (Tm = 62°C), and excellent refolding abilities while maintaining their capacity to bind and neutralize CHIKV.

Structural and kinetic evidence of aging after organophosphate inhibition of human Cathepsin A.

Human Cathepsin A (CatA) is a lysosomal serine carboxypeptidase of the renin-angiotensin system (RAS) and is structurally similar to acetylcholinesterase (AChE). CatA can remove the C-terminal amino acids of endothelin I, angiotensin I, Substance P, oxytocin, and bradykinin, and can deamidate neurokinin A. Proteomic studies identified CatA and its homologue, SCPEP1, as potential targets of organophosphates (OP). CatA could be stably inhibited by low µM to high nM concentrations of racemic sarin (GB), soman (GD), cyclosarin (GF), VX, and VR within minutes to hours at pH 7. Cyclosarin was the most potent with a kinetically measured dissociation constant (KI) of 2 µM followed by VR (KI = 2.8 µM). Bimolecular rate constants for inhibition by cyclosarin and VR were 1.3 × 103 M-1sec-1 and 1.2 × 103 M-1sec-1, respectively, and were approximately 3-orders of magnitude lower than those of human AChE indicating slower reactivity. Notably, both AChE and CatA bound diisopropylfluorophosphate (DFP) comparably and had KIDFP = 13 µM and 11 µM, respectively. At low pH, greater than 85% of the enzyme spontaneously reactivated after OP inhibition, conditions under which OP-adducts of cholinesterases irreversibly age. At pH 6.5 CatA remained stably inhibited by GB and GF and <10% of the enzyme spontaneously reactivated after 200 h. A crystal structure of DFP-inhibited CatA was determined and contained an aged adduct. Similar to AChE, CatA appears to have a "backdoor" for product release. CatA has not been shown previously to age. These results may have implications for: OP-associated inflammation; cardiovascular effects; and the dysregulation of RAS enzymes by OP.

Codon harmonization reduces amino acid misincorporation in bacterially expressed P. falciparum proteins and improves their immunogenicity.

Codon usage frequency influences protein structure and function. The frequency with which codons are used potentially impacts primary, secondary and tertiary protein structure. Poor expression, loss of function, insolubility, or truncation can result from species-specific differences in codon usage. "Codon harmonization" more closely aligns native codon usage frequencies with those of the expression host particularly within putative inter-domain segments where slower rates of translation may play a role in protein folding. Heterologous expression of Plasmodium falciparum genes in Escherichia coli has been a challenge due to their AT-rich codon bias and the highly repetitive DNA sequences. Here, codon harmonization was applied to the malarial antigen, CelTOS (Cell-traversal protein for ookinetes and sporozoites). CelTOS is a highly conserved P. falciparum protein involved in cellular traversal through mosquito and vertebrate host cells. It reversibly refolds after thermal denaturation making it a desirable malarial vaccine candidate. Protein expressed in E. coli from a codon harmonized sequence of P. falciparum CelTOS (CH-PfCelTOS) was compared with protein expressed from the native codon sequence (N-PfCelTOS) to assess the impact of codon usage on protein expression levels, solubility, yield, stability, structural integrity, recognition with CelTOS-specific mAbs and immunogenicity in mice. While the translated proteins were expected to be identical, the translated products produced from the codon-harmonized sequence differed in helical content and showed a smaller distribution of polypeptides in mass spectra indicating lower heterogeneity of the codon harmonized version and fewer amino acid misincorporations. Substitutions of hydrophobic-to-hydrophobic amino acid were observed more commonly than any other. CH-PfCelTOS induced significantly higher antibody levels compared with N-PfCelTOS; however, no significant differences in either IFN-γ or IL-4 cellular responses were detected between the two antigens.

Sequence Tolerance of a Single-Domain Antibody with a High Thermal Stability: Comparison of Computational and Experimental Fitness Profiles.

The sequence fitness of a llama single-domain antibody with an unusually high thermal stability is explored by a combined computational and experimental study. Starting with the X-ray crystallographic structure, RosettaBackrub simulations were applied to model sequence-structure tolerance profiles and identify key substitution sites. From the model calculations, an experimental site-directed mutagenesis was used to produce a panel of mutants, and their melting temperatures were determined by thermal denaturation. The results reveal a sequence fitness of an excess stability of approximately 12 °C, a value taken from a decrease in the melting temperature of an electrostatic charge-reversal substitution in the CRD3 without a deleterious effect on the binding affinity to the antigen. The tolerance for the disruption of antigen recognition without loss in the thermal stability was demonstrated by the introduction of a proline in place of a tyrosine in the CDR2, producing a mutant that eliminated binding. To further assist the sequence design and the selection of engineered single-domain antibodies, an assessment of different computational strategies is provided of their accuracy in the detection of substitution "hot spots" in the sequence tolerance landscape.

Paired Carboxylic Acids in Enzymes and Their Role in Selective Substrate Binding, Catalysis, and Unusually Shifted pKa Values.

Cathepsin A (CatA, EC 3.4.16.5, UniProtKB P10619 ) is a human lysosomal carboxypeptidase. Counterintuitively, crystal structures of CatA and its homologues show a cluster of Glu and Asp residues binding the C-terminal carboxylic acid of the product or inhibitor. Each of these enzymes functions in an acidic environment and contains a highly conserved pair of Glu residues with side chain carboxyl group oxygens that are approximately 2.3-2.6 Šapart. In small molecules, carboxyl groups separated by ∼3 Šcan overcome the repulsive interaction by protonation of one of the two groups. The pKa of one group increases (pKa ∼ 11) and can be as much as ∼6 pH units higher than the paired group. Consequently, at low and neutral pH, one carboxylate can carry a net negative charge while the other can remain protonated and neutral. In CatA, E69 and E149 form a Glu pair that is important to catalysis as evidenced by the 56-fold decrease in kcat/Km in the E69Q/E149Q variant. Here, we have measured the pH dependencies of log(kcat), log(Km), and log(kcat/Km) for wild type CatA and its variants and have compared the measured pKa with calculated values. We propose a substrate-assisted mechanism in which the high pKa of E149 (>8.5) favors the binding of the carboxylate form of the substrate and promotes the abstraction of the proton from H429 of the catalytic triad effectively decreasing its pKa in a low-pH environment. We also identify a similar motif consisting of a pair of histidines in S-formylglutathione hydrolase.

Proteolytic cleavage of host proteins by the Group IV viral proteases of Venezuelan equine encephalitis virus and Zika virus.

The alphaviral nonstructural protein 2 (nsP2) cysteine proteases (EC 3.4.22.-) are essential for the proteolytic processing of the nonstructural (ns) polyprotein and are validated drug targets. A common secondary role of these proteases is to antagonize the effects of interferon (IFN). After delineating the cleavage site motif of the Venezuelan equine encephalitis virus (VEEV) nsP2 cysteine protease, we searched the human genome to identify host protein substrates. Here we identify a new host substrate of the VEEV nsP2 protease, human TRIM14, a component of the mitochondrial antiviral-signaling protein (MAVS) signalosome. Short stretches of homologous host-pathogen protein sequences (SSHHPS) are present in the nonstructural polyprotein and TRIM14. A 25-residue cyan-yellow fluorescent protein TRIM14 substrate was cleaved in vitro by the VEEV nsP2 protease and the cleavage site was confirmed by tandem mass spectrometry. A TRIM14 cleavage product also was found in VEEV-infected cell lysates. At least ten other Group IV (+)ssRNA viral proteases have been shown to cleave host proteins involved in generating the innate immune responses against viruses, suggesting that the integration of these short host protein sequences into the viral protease cleavage sites may represent an embedded mechanism of IFN antagonism. This interference mechanism shows several parallels with those of CRISPR/Cas9 and RNAi/RISC, but with a protease recognizing a protein sequence common to both the host and pathogen. The short host sequences embedded within the viral genome appear to be analogous to the short phage sequences found in a host's CRISPR spacer sequences. To test this algorithm, we applied it to another Group IV virus, Zika virus (ZIKV), and identified cleavage sites within human SFRP1 (secreted frizzled related protein 1), a retinal Gs alpha subunit, NT5M, and Forkhead box protein G1 (FOXG1) in vitro. Proteolytic cleavage of these proteins suggests a possible link between the protease and the virus-induced phenotype of ZIKV. The algorithm may have value for selecting cell lines and animal models that recapitulate virus-induced phenotypes, predicting host-range and susceptibility, selecting oncolytic viruses, identifying biomarkers, and de-risking live virus vaccines. Inhibitors of the proteases that utilize this mechanism may both inhibit viral replication and alleviate suppression of the innate immune responses.

Analysis of Group IV Viral SSHHPS Using In Vitro and In Silico Methods.

Alphaviral enzymes are synthesized in a single polypeptide. The nonstructural polyprotein (nsP) is processed by its nsP2 cysteine protease to produce active enzymes essential for viral replication. Viral proteases are highly specific and recognize conserved cleavage site motif sequences (~6-8 amino acids). In several Group IV viruses, the nsP protease(s) cleavage site motif sequences can be found in specific host proteins involved in generating the innate immune responses and, in some cases, the targeted proteins appear to be linked to the virus-induced phenotype. These viruses utilize short stretches of homologous host-pathogen protein sequences (SSHHPS) for targeted destruction of host proteins. To identify SSHHPS the viral protease cleavage site motif sequences can be inputted into BLAST and the host genome(s) can be searched. Cleavage initially can be tested using the purified nsP viral protease and fluorescence resonance energy transfer (FRET) substrates made in E. coli. The FRET substrates contain cyan and yellow fluorescent protein and the cleavage site sequence (CFP-sequence-YFP). This protease assay can be used continuously in a plate reader or discontinuously in SDS-PAGE gels. Models of the bound peptide substrates can be generated in silico to guide substrate selection and mutagenesis studies. CFP/YFP substrates have also been utilized to identify protease inhibitors. These in vitro and in silico methods can be used in combination with cell-based assays to determine if the targeted host protein affects viral replication.

Mutation of Asn-475 in the Venezuelan Equine Encephalitis Virus nsP2 Cysteine Protease Leads to a Self-Inhibited State.

The alphaviral nsP2 cysteine protease of the Venezuelan equine encephalitis virus (VEEV) is a validated antiviral drug target. Clan CN proteases contain a cysteine protease domain that is intimately packed with an S-adenosyl-l-methionine-dependent RNA methyltransferase (SAM MTase) domain. Within a cleft formed at the interface of these two domains, the peptide substrate is thought to bind. The nucleophilic cysteine can be found within a conserved motif, 475NVCWAK480, which differs from that of papain (22CGSCWAFS29). Mutation of the motif residue, N475, to alanine unexpectedly produced a self-inhibited state in which the N-terminal residues flipped into the substrate-binding cleft. Notably, the N-terminal segment was not hydrolyzed-consistent with a catalytically incompetent state. The N475A mutation resulted in a 70-fold decrease in kcat/Km. A side chain-substrate interaction was predicted by the structure; the S701A mutation led to a 17-fold increase in Km. An Asn at the n-2 position relative to the Cys was also found in the coronaviral papain-like proteases/deubiquitinases (PLpro) of the SARS and MERS viruses, and in several papain-like human ubiquitin specific proteases (USP). The large conformational change in the N475A variant suggests that Asn-475 plays an important role in stabilizing the N-terminal residues and in orienting the carbonyl during nucleophilic attack but does not directly hydrogen bond the oxyanion. The state trapped in crystallo is an unusual result of site-directed mutagenesis but reveals the role of this highly conserved Asn and identifies key substrate-binding contacts that may be exploited by peptide-like inhibitors.

Label free checkerboard assay to determine overlapping epitopes of Ebola virus VP-40 antibodies using surface plasmon resonance.

Immunoassay formats, in which antibodies provide sensitivity and specificity, are often utilized to provide rapid and simple diagnostic tests. Surface plasmon resonance is frequently used to evaluate the suitability of antibodies by determining binding kinetics to agents or surrogate antigens. We used SPR to evaluate a number of commercial monoclonal antibodies as well as single domain antibodies produced in-house. All the antibodies targeted the Ebola virus viral protein 40 (VP40). We determined the ability of each antibody to bind to immobilized VP40, and ensured they did not bind Ebola glycoprotein or the nucleoprotein. A subset of the monoclonal antibodies was immobilized to characterize antigen capture in solution. It can be advantageous to utilize antibodies that recognize distinct epitopes when choosing reagents for detection and diagnostic assays. We determined the uniqueness of the epitope recognized by the anti-VP40 antibodies using a checkerboard format that exploits the 6×6 array of interactions monitored by the Bio-Rad ProteOn XPR36 SPR instrument. The results demonstrate the utility of surface plasmon resonance to characterize monoclonal and recombinant antibodies. Additionally, the analysis presented here enabled the identification of pairs of anti-VP40 antibodies which could potentially be utilized in sandwich type immunoassays for the detection of Ebola virus.

Stability of isolated antibody-antigen complexes as a predictive tool for selecting toxin neutralizing antibodies.

Ricin is an A-B ribosome inactivating protein (RIP) toxin composed of an A-chain subunit (RTA) that contains a catalytic N-glycosidase and a B-chain (RTB) lectin domain that binds cell surface glycans. Ricin exploits retrograde transport to enter into the Golgi and the endoplasmic reticulum, and then dislocates into the cytoplasm where it can reach its substrate, the rRNA. A subset of isolated antibodies (Abs) raised against the RTA subunit protect against ricin intoxication, and RTA-based vaccine immunogens have been shown to provide long-lasting protective immunity against the holotoxin. Anti-RTA Abs are unlikely to cross a membrane and reach the cytoplasm to inhibit the enzymatic activity of the A-chain. Moreover, there is not a strict correlation between the apparent binding affinity (Ka) of anti-RTA Abs and their ability to successfully neutralize ricin toxicity. Some anti-RTA antibodies are toxin-neutralizing, whereas others are not. We hypothesize that neutralizing anti-RTA Abs may interfere selectively with conformational change(s) or partial unfolding required for toxin internalization. To test this hypothesis, we measured the melting temperatures (Tm) of neutralizing single-domain Ab (sdAb)-antigen (Ag) complexes relative to the Tm of the free antigen (Tm-shift = Tmcomplex - TmAg), and observed increases in the Tmcomplex of 9-20 degrees. In contrast, non-neutralizing sdAb-Ag complexes shifted the TmComplex by only 6-7 degrees. A strong linear correlation (r2 = 0.992) was observed between the magnitude of the Tm-shift and the viability of living cells treated with the sdAb and ricin holotoxin. The Tm-shift of the sdAb-Ag complex provided a quantitative biophysical parameter that could be used to predict and rank-order the toxin-neutralizing activities of Abs. We determined the first structure of an sdAb-RTA1-33/44-198 complex, and examined other sdAb-RTA complexes. We found that neutralizing sdAb bound to regions involved in the early stages of unfolding. These Abs likely interfere with steps preceding or following endocytosis that require conformational changes. This method may have utility for the characterization or rapid screening of other Ab that act to prevent conformational changes or unfolding as part of their mechanism of action.

Kinetic, Mutational, and Structural Studies of the Venezuelan Equine Encephalitis Virus Nonstructural Protein 2 Cysteine Protease.

The Venezuelan equine encephalitis virus (VEEV) nonstructural protein 2 (nsP2) cysteine protease (EC 3.4.22.-) is essential for viral replication and is involved in the cytopathic effects (CPE) of the virus. The VEEV nsP2 protease is a member of MEROPS Clan CN and characteristically contains a papain-like protease linked to an S-adenosyl-l-methionine-dependent RNA methyltransferase (SAM MTase) domain. The protease contains an alternative active site motif, (475)NVCWAK(480), which differs from papain's (CGS(25)CWAFS), and the enzyme lacks a transition state-stabilizing residue homologous to Gln-19 in papain. To understand the roles of conserved residues in catalysis, we determined the structure of the free enzyme and the first structure of an inhibitor-bound alphaviral protease. The peptide-like E64d inhibitor was found to bind beneath a ß-hairpin at the interface of the SAM MTase and protease domains. His-546 adopted a conformation that differed from that found in the free enzyme; one or both of the conformers may assist in leaving group departure of either the amine or Cys thiolate during the catalytic cycle. Interestingly, E64c (200 µM), the carboxylic acid form of the E64d ester, did not inhibit the nsP2 protease. To identify key residues involved in substrate binding, a number of mutants were analyzed. Mutation of the motif residue, N475A, led to a 24-fold reduction in kcat/Km, and the conformation of this residue did not change after inhibition. N475 forms a hydrogen bond with R662 in the SAM MTase domain, and the R662A and R662K mutations both led to 16-fold decreases in kcat/Km. N475 forms the base of the P1 binding site and likely orients the substrate for nucleophilic attack or plays a role in product release. An Asn homologous to N475 is similarly found in coronaviral papain-like proteases (PLpro) of the Severe Acute Respiratory Syndrome (SARS) virus and Middle East Respiratory Syndrome (MERS) virus. Mutation of another motif residue, K480A, led to a 9-fold decrease in kcat and kcat/Km. K480 likely enhances the nucleophilicity of the Cys. Consistent with our substrate-bound models, the SAM MTase domain K706A mutation increased Km 4.5-fold to 500 µM. Within the ß-hairpin, the N545A mutation slightly but not significantly increased kcat and Km. The structures and identified active site residues may facilitate the discovery of protease inhibitors with antiviral activity.

Comparison of Replica Exchange Simulations of a Kinetically Trapped Protein Conformational State and its Native Form.

Recently an X-ray crystallographic structure of a single-domain antibody was reported with the protein chain trapped in a rare homodimeric form. One of the conformers appears to exhibit a misfolded region, and thus presumably the configurational stability is less favorable. To investigate whether simulation methods can detect any difference between the conformers and buttress the notion that one conformation is trapped on a pathway that incurs lower activation energy to unfold, adaptive temperature-based replica exchange simulations were applied to each chain to model conformational transitions. Simulation results found that the observed crystallographic difference between the two chains in the complementarity determining region CDR2 induces a stark distinction in conformational populations on the energy landscape. An appraisal of the energetic difference between the CDR2 conformations at 300 K revealed a localized order-disorder free-energy transition of roughly equivalent to two peptide hydrogen bonds in solution. It was also found that interconversion between the conformers is slower than the rate to unfold and that near an unfolding transition temperature one conformer retained a greater fraction of native-like contacts and energy over a longer time span before fully populating the denatured state, thus verifying the coexistence of a metastable conformation in the crystallographic assembly.

Can template-based protein models guide the design of sequence fitness for enhanced thermal stability of single domain antibodies?

We investigate the practical use of comparative (template-based) protein models in replica-exchange simulations of single-domain antibody (sdAb) chains to evaluate if the models can correctly predict in rank order the thermal susceptibility to unfold relative to experimental melting temperatures. The baseline model system is the recently determined crystallographic structure of a llama sdAb (denoted as A3), which exhibits an unusually high thermal stability. An evaluation of the simulation results for the A3 comparative model and crystal structure shows that, despite the overall low Cα root-mean-square deviation between the two structures, the model contains misfolded regions that yields a thermal profile of unraveling at a lower temperature. Yet comparison of the simulations of four different comparative models for sdAb A3, C8, A3C8 and E9, where A3C8 is a design of swapping the sequence of the complementarity determining regions of C8 onto the A3 framework, discriminated among the sequences to detect the highest and lowest experimental melting transition temperatures. Further structural analysis of A3 for selected alanine substitutions by a combined computational and experimental study found unexpectedly that the comparative model performed admirably in recognizing substitution 'hot spots' when using a support-vector machine algorithm.

A conformational change in the peripheral anionic site of Torpedo californica acetylcholinesterase induced by a bis-imidazolium oxime.

As part of ongoing efforts to design improved nerve agent antidotes, two X-ray crystal structures of Torpedo californica acetylcholinesterase (TcAChE) bound to the bis-pyridinium oxime, Ortho-7, or its experimental bis-imidazolium analogue, 2BIM-7, were determined. Bis-oximes contain two oxime groups connected by a hydrophobic linker. One oxime group of Ortho-7 binds at the entrance to the active-site gorge near Trp279, and the second binds at the bottom near Trp84 and Phe330. In the Ortho-7-TcAChE complex the oxime at the bottom of the gorge was directed towards the nucleophilic Ser200. In contrast, the oxime group of 2BIM-7 was rotated away from Ser200 and the oxime at the entrance induced a significant conformational change in the peripheral anionic site (PAS) residue Trp279. The conformational change alters the surface of the PAS and positions the imidazolium oxime of 2BIM-7 further from Ser200. The relatively weaker binding and poorer reactivation of VX-inhibited, tabun-inhibited or sarin-inhibited human acetylcholinesterase by 2BIM-7 compared with Ortho-7 may in part be owing to the unproductively bound states caught in crystallo. Overall, the reactivation efficiency of 2BIM-7 was comparable to that of 2-pyridine aldoxime methyl chloride (2-PAM), but unlike 2-PAM the bis-imidazolium oxime lacks a fixed charge, which may affect its membrane permeability.

A novel Vibrio beta-glucosidase (LamN) that hydrolyzes the algal storage polysaccharide laminarin.

The metabolic versatility, tractability and rapid growth potential of the Vibrio spp. have made them increasingly attractive systems for investigating carbon cycling in the marine environment. In this study, an in silico subtractive proteomic strategy was used to identify a novel 101 kDa GH3 family ß-glucosidase (LamN) that was found in bioluminescent Vibrio campbellii strains capable of utilizing the algal storage glucan laminarin. A heterologous overexpression system verified the sequence-predicted function of LamN as it enabled the growth of Escherichia coli on laminarin as a sole carbon source. Quantitative reverse transcription PCR analyses revealed that V. campbellii grown on laminarin demonstrated a 4- to 314-fold induction of lamN gene expression when compared to the same strains grown on glucose or glycerol. Corresponding tandem mass spectrometric analyses detected LamN protein expression only in cells grown on laminarin. Heterologous expression, purification and biochemical characterization identified LamN as a heat stable laminarinase with ß-1,3, ß-1,4 and ß-1,6 glucosidase activity. Collectively, these data identify an enzyme that may allow V. campbellii to exploit some of the most abundant polysaccharides associated with deteriorating phytoplankton blooms and provide support for the potential involvement of V. campbellii in the formation of bioluminescent milky seas.