key: cord-255815-5d9bqji0
authors: Malik, Ajamaluddin; Alsenaidy, Mohammad A.
title: MERS‐CoV papain-like protease (PL(pro)): expression, purification, and spectroscopic/thermodynamic characterization
date: 2017-05-30
journal: 3 Biotech
DOI: 10.1007/s13205-017-0744-3
sha: 
doc_id: 255815
cord_uid: 5d9bqji0

Within a decade, MERS-CoV emerged with nearly four times higher case fatality rate than an earlier outbreak of SARS-CoV and spread out in 27 countries in short span of time. As an emerging virus, combating it requires an in-depth understanding of its molecular machinery. Therefore, conformational characterization studies of coronavirus proteins are necessary to advance our knowledge of the matter for the development of antiviral therapies. In this study, MERS-CoV papain-like protease (PL(pro)) was recombinantly expressed and purified. Thermal folding pathway and thermodynamic properties were characterized using dynamic multimode spectroscopy (DMS) and thermal shift assay. DMS study showed that the PL(pro) undergoes a single thermal transition and follows a pathway of two-state folding with T (m) and van’t Hoff enthalpy values of 54.4 ± 0.1 °C and 317.1 ± 3.9 kJ/mol, respectively. An orthogonal technique based on intrinsic tryptophan fluorescence also showed that MERS-CoV PL(pro) undergoes a single thermal transition and unfolds via a pathway of two-state folding with a T (m) value of 51.4 °C. Our findings provide significant understandings of the thermodynamic and structural properties of MERS-CoV PL(pro).

Frequent fatal coronavirus outbreaks in humans and animals have caused serious concerns in the healthcare sector, scientific community, and animal husbandry. First human outbreak of severe acute respiratory syndrome coronavirus (SARS-CoV) in 2002 caused life-threatening atypical pneumonia in more than 8000 people in 26 countries with a case fatality rate (CFR) of 10% (Pillaiyar et al. 2015; Al-Tawfiq et al. 2014; WHO 2016b) . Another lethal coronavirus outbreak emerged in the Arab peninsula countries in 2012 which was caused by what is now known as the Middle East respiratory syndrome coronavirus (MERS-CoV) . From September 2012 to December 2016, 1864 laboratory-confirmed MERS-CoV cases of infection in 27 countries with 659 mortalities (nearly four times higher CFR than SARS-CoV) have been reported (Al-Tawfiq et al. 2016; WHO 2016a) . Coronavirus survivors after acute infections suffer from many health issues and require long-term medical assistance (Han et al. 2003; Ong et al. 2005; Chan et al. 2003; Leow et al. 2005; Siu 2016; Cha et al. 2016) . In addition to SARS-CoV and MERS-CoV, at least four other pathological coronaviruses (HCoV-OC43, HCoV-229E, HCoV-HKU1, and HCoV-NL63) are continuously circulating in humans causing relatively mild respiratory conditions that may in some instances escalate to severe pathological illnesses (Mackay et al. 2012; Carbajo-Lozoya et al. 2012; Simon et al. 2007) . Coronaviruses have also caused deadly diseases in animals, leading to huge economic losses in the animal husbandry sector (Vlasova et al. 2014; Lee and Lee 2014; Sun et al. 2016) . Moreover, high mutation and recombination rates in coronaviruses allow them to cross species barriers and adapt to new hosts more easily (Denison et al. 2011; Lau and Chan 2015) .

Viral proteases are essential for pathogenesis and virulence. Like all coronaviruses, MERS-CoV contains two cysteine proteases (main protease and papain-like protease) which processes viral nonstructural polypeptides (Kilianski et al. 2013; Hilgenfeld 2014) . MERS-CoV main protease (M pro , also called the 3C-like protease, 3CL pro ) cleaves at eleven sites, while MERS-CoV papain-like protease (PL pro ) cuts at three sites on the nonstructural polypeptides and releases mature nonstructural proteins (Hilgenfeld 2014) . Thus, MERS-CoV proteases make up a suitable target for antiviral therapies.

MERS-CoV open-reading frame 1 (ORF1) encodes two large polyproteins (pp1a and pp1b). MERS-CoV PL pro domain is encoded on the pp1a proteins (residue 1484-1800) (Yang et al. 2014; Hilgenfeld 2014; Kilianski et al. 2013) . Like other coronaviruses, MERS-CoV PL pro contains a catalytic triad and exhibits similar proteolytic, deubiquitination, and ISG15-linked ISGylation properties (Lin et al. 2014; Chen et al. 2007; Clementz et al. 2010; Yang et al. 2014; Zheng et al. 2008) .

In this study, we expressed and purified MERS-CoV PL pro . Thermal stability was studied by thermal shift assay using intrinsic fluorescence and Dynamic Multimode Spectroscopy (DMS). MERS-CoV PL pro was found to unfold via a single thermal transition and follows a pathway of two-state folding. This study will not only help in the understanding of the folding and stability of MERS-CoV PL pro but also could help shed some light on other deubiquitinating enzymes with the similar folding scaffold.

The ORF of MERS-CoV PL pro (1484-1800 polyprotein residues, GenBank accession number NC_019843.2) was cloned into pET28a plasmid under T7 promoter as published before (Lin et al. 2014) . The codon was optimized (GenScript, USA) and cloned between NcoI and XhoI sites which was in frame of C-terminal His tag present on the vector. E. coli BL21 (DE3) pLysS was used for the expression of recombinant protein. Low-molecular weight protein markers, prepacked Ni-NTA, and Superdex 75 columns were from Amersham Biosciences (United Kingdom). Chicken egg lysozyme was from USB Corporation, USA. Benzonase, ANS, and kanamycin from Sigma. IPTG was purchased from Bio Basic, Canada. All other chemicals used in this study were of reagent grade. Cary 60 spectrometer and Cary Eclipse spectrofluorometer were from Agilent technologies, USA. AKTA purification system was from Amersham Biosciences (United Kingdom) and SDS-PAGE assembly from Bio-Rad (USA). Thermomixer and benchtop cooling centrifuge were from Eppendorf, Germany. Innova 44R Shaking incubator was from New Brunswick, Germany. Chirascan-Plus spectropolarimeter was from Applied photophysics, United Kingdom.

Expression and purification of MERS-CoV PL pro in E. coli BL21 (DE3) pLysS E. coli BL21 (DE3) pLysS harboring pET28a-MPL pro was used for expression of MERS-CoV PL pro . Protein expression and soluble protein extraction were performed as described in Lin et al. (2014) . Purification of MERS-CoV PL pro was performed with minor modification of an earlier published protocol (Lin et al. 2014) . Briefly, 1 mM DTT was used throughout unless described. Cleared crude lysate was passed through a 1-mL Ni-NTA column pre-equilibrated with 20 mM Tris, pH 8.5, 500 mM NaCl, 10 mM imidazole, and 1 mM DTT and washed with 20 CV equilibration buffer. Bound protein was eluted with a linear gradient of 0-50% Buffer B (equilibration buffer containing 500 mM Imidazole) at 1 mL/min flow rate on AKTA purification system. The purity of eluted fractions was analyzed on SDS-PAGE.

The prepacked Superdex 75 equilibrated with 20 mM Tris, pH 8.5, 100 mM NaCl, and 1 mM DTT was calibrated with proteins of known molecular weight. Subsequently, Ni-NTA-purified MERS-CoV PL pro was further purified using Superdex 75 column. The purity of eluted fractions was analyzed on SDS-PAGE. Highly pure fractions were pooled, aliquoted, and stored at -80°C.

Before analysis, the frozen aliquots were thawed and centrifuged at 13,000 rpm for 15 min at 4°C. Protein concentration was determined spectrophotometrically at 280 nm using a molar extinction coefficient of 42,400 M -1 cm -1 .

Tryptophan fluorescence spectra (50 lg/mL) of MERS-CoV PL pro were recorded using a Cary Eclipse Fluorescence Spectrophotometer in a 10-mm-path length cuvette. To measure tryptophan fluorescence, MERS-CoV PL pro sample was excited at 295 nm and emission spectra were collected between 305 and 400 nm (5 nm excitation and 5 nm emission bandwidth). Temperature melting studies of MERS-CoV PL pro were performed at 10°C increments as well as a linear increase of 1°C/min. MERS-CoV PL pro incubated at different temperatures from 20 to 90°C in Peltier-controlled Cary eclipse fluorometer. The temperature of the protein samples was monitored using an internal temperature probe. When the desired temperature was reached, the sample was allowed to equilibrate for 2 min before tryptophan fluorescence spectra were measured. The maximum fluorescence intensity (I max ) and maximum fluorescence wavelength (k max ) were plotted with respect to temperature. In a similar experiment, MERS-CoV PL pro was gradually heated from 20 to 80°C at a rate of 1°C/min during which tryptophan fluorescence was measured by exciting at 295 nm and collecting at 330 and 350 nm to obtain the temperature melting curve.

To study the secondary structure of MERS-CoV PL pro in terms of conformational and thermal stability, dynamic multimode spectroscopy was applied. The measurement was performed using Chirascan-Plus spectrophotometer, calibrated with (1S)-(?)-10-camphorsulfonic acid. In this study, 0.2 mg/mL of MERS-CoV PL pro was gradually heated from 20 to 94°C at 1°C/mL rate. Internal thermal probe was inserted in the 0.1-cm-path length cuvette to precisely monitor the actual temperature of the samples. Far-UV CD spectra from 200 to 250 nm were recorded at each temperature. Thermal transition data were processed using manufacturer's Global 3 software.

The steady-state kinetics of MERS-CoV PL pro was measured as described in Lin et al. (2014) , with slight modification. Briefly, 50 lM fluorogenic peptidyl substrate, Dabcyl-FRLKGGAPIKGV-Edans (GenScript), was mixed with different concentrations of MERS-CoV PL pro (6.8-0.1 lM, in twofold serial dilution) using 50 mM phosphate at pH 6.5 as a buffer at room temperature. The fluorescence signal was measured for 30 min at 60-s intervals in a Hidex Chameleon plate reader using 340 nm (excitation) and 535 nm (emission filter) with 20% gain.

Expression and purification of recombinant MERS-CoV PL pro In this study, MERS-CoV PL pro was overexpressed in E. coli BL21 (DE3) pLysS. MERS-CoV PL pro was purified in two-step chromatography as described in Lin et al. (2014) . Ni-NTA elute contained minor impurities (data not shown). When Ni-NTA elute passed through gel filtration column, one major symmetrical sharp peak was obtained (Fig. 1a) . SDS-PAGE analysis of the pooled fractions showed the yield of highly pure MERS-CoV PL pro (Fig. 1b) . We obtained nearly 10 mg of MERS-CoV PL pro from a 1-L shake flask culture. In an earlier study, the yield of soluble MERS-CoV PL pro was strain dependent and the best yield (*52 mg purified protein from 1 L shake flask culture) was obtained with E. coli BL21 (DE3) STAR strain (Lin et al. 2014) . The difference in the yield of MERS-CoV PL pro was due to the different strain's genetic background. 

Intrinsic fluorescence spectroscopy is a highly sensitive tool, which provides information about the microenvironment of Trp and Tyr residues within proteins. Tyrosine emission maximum is less sensitive to its local environment compared to tryptophan. Indole ring in the tryptophan residues undergoes two isoenergetic transitions which causes polarity sensitivity, while tyrosine undergoes through a single electronic state (Ghisaidoobe and Chung 2014) . During the course of protein unfolding, the polarity of fluorophore microenvironment changes, which in turn leads to changes in the maximum fluorescence intensity (I max ) as well as maximum fluorescence wavelength (k max ). Therefore, intrinsic tryptophan fluorescence emission is sensitive to the tertiary structure and detects subtle protein conformational changes in solution. It has been frequently used for the characterization of protein's conformational changes under different stress conditions (Kumar et al. 2005; Xiao et al. 2015) .

A previous study analyzing the quaternary structure of MERS-CoV PL pro using analytical ultracentrifugation technique has shown that PL pro is found in the monomeric state (Lin et al. 2014) . MERS-CoV PL pro contains ten tyrosine and five tryptophan residues (Fig. 2a) . MERS-CoV PL pro consists of two domains: N-terminal ubiquitinlike (Ubl) domain and a catalytic core domain. The N-terminal Ubl domain consists of 62 residues and contained one a-helix, one 3 10 -helix and five b-strands. The substrate-binding region is solvent exposed and comprised of the right-hand scaffold (Lei et al. 2014) . Three tryptophan residues (W93, 243 and 303) are buried, while W187 and W190 are surface exposed (Fig. 2b) . Figure 3a shows the decrease in intrinsic fluorescence of MERS-CoV PL pro with increasing temperature. At low temperature (20°C), highest fluorescent intensity (I max ) was observed with maximum fluorescence wavelength (k max ) at 343 nm, indicating an overall localization of all five tryptophans in the partially hydrophobic environment. As the temperature increases gradually from 20 to 90°C, I max decreased with major transitions occurring between 50 and 60°C (Fig. 3a,  b) . Initially, k max increased from 343 to 344 nm when the temperature was increased from 20 to 30°C. As the temperature was increased further, a blue-shift in k max with increasing temperature was observed, indicating conformational rearrangements in different temperature regimes. The major blue-shift in k max was found between 50 and 60°C (Fig. 3c) . The 14-nm blue-shift of tryptophan fluorescence spectrum indicated that the microenvironment of the tryptophan residues is becoming more hydrophobic during thermal denaturation. MERS-CoVPL pro thermal unfolding is different from most other proteins. Commonly, protein unfolding fluorescence spectra are characterized by a long wavelength shift ''red-shift.'' But some proteins, Fig. 2 a Sequence of C-terminal His-tagged MERS-CoV PL pro showing ten Tyr and five Trp residues, which are highlighted in green and blue, respectively. Potential internal protease site is highlighted in red and His-tagged residues are underlined. b Three-dimensional structure of MERS-CoV PL pro (PDB ID: 4R3D) with the side chains of the five Trp residues being numbered and colored in blue including MERS-CoV PL pro , exhibit blue-shift upon denaturation (Slutskaya et al. 2015; Duy and Fitter 2006; Pattanaik et al. 2003) . Pig pancreatic a-amylase showed red-shifted fluorescence spectra when chemically unfolded and showed blue-shifted spectra during thermal unfolding (Duy and Fitter 2006) . Equine lysozyme first exhibits a blue-shift transition at lower temperature and red-shift above 50°C (Morozova et al. 1991) .

To obtain temperature melting curve, MERS-CoVPL pro was gradually heated from 20 to 80°C at 1°C/min and the ratio of 330/350 nm tryptophan fluorescence was plotted with respect to temperature (Fig. 4) . Data were fitted according to the equation f = y 0 ? a/(1 ? exp(-(xx 0 )/ b)) with an r 2 value of 0.9910. Our results showed that MERS-CoV PL pro was moderately stable and unfolds via a single transition with a T m value of 51.4°C.

Since MERS-CoV PL pro is a protease, it may undergo autolysis during thermal shift assay. It contains one potential autolysis site (LKGG) as shown in Fig. 2a . To evaluate the extent of autolysis as well as the extent of irreversible thermal unfolding and aggregation, MERS-CoV PL pro (0.2 mg/mL) was incubated on a thermomixer from 20 to 70°C with 10°C intervals. Six samples were gradually heated and equilibrated at the respective temperatures for 3 min. When the desired temperatures were attained, the samples were removed, kept on ice, and centrifuged at 13,000 rpm for 15 min to remove any forming aggregates. Equal volumes of the supernatant were analyzed (Fig. 5) . If MERS-CoV PL pro would undergo autolysis during the course of thermal incubation, we would expect to see the appearance of two or more bands of MERS-CoV PL pro fragments on SDS-PAGE gels or at least a decrease in the band intensity of MERS-CoV PL pro . We found that the intensity of MERS-CoV PL pro band was apparently unchanged, indicating that autolysis was not occurring during the thermal shift assays applied in this study. In another scenario, irreversible unfolding aggregates may form during thermal shift assay. If this is the case, then aggregated protein will be pelleted down upon centrifugation and band intensity would have decreased in the supernatant sample. Our result showed that the band intensity of the supernatant samples incubated from 20 to 70°C was apparently unchanged (Fig. 5) , indicating Fig. 3 Thermally induced structural changes in MERS-CoV PL pro as monitored by the intrinsic tryptophan fluorescence spectroscopy. a To monitor tryptophan fluorescence at different temperatures, MERS-CoV PL pro was slowly heated and allowed to equilibrate for 2 min at the respective temperatures. The sample was excited at 295 nm and the emission spectra were collected from 305 to 400 nm (5 nm excitation and 5 nm emission bandwidth). b Effect of temperature on the tryptophan emission intensity showing the decrease of I max with increasing temperature. Major transition occurred between 50 and 60°C. c Effect of temperature on the tryptophan maximum emission wavelength (k max ). Blue-shift was observed during thermal unfolding of MERS-CoV PL pro . Initially, red-shift was observed when the temperature was increased from 20 to 30°C. Further increase in temperature leads to blue-shift with major changes occurring between 50 and 60°C no or insignificant aggregation occurring during thermal shift assays.

Information about thermal stability, unfolding pathway, and secondary structure of MERS-CoV PL pro was obtained using an orthogonal method. We employed DMS, a newly developed information-rich experimental technique, to obtain spectroscopic and thermodynamic data of melting temperature (T m ) and van't Hoff enthalpy (DH VH ) (John and Weeks 2000; Greenfield 2006; Al-Ahmady et al. 2012) . Temperature-induced secondary structural changes in MERS-CoV PL pro in the far-UV region were monitored to study thermodynamic parameters. Moreover, DMS also characterized thermal unfolding pathway and identified the number of folding intermediate species and their relative concentrations during the unfolding process (Malik et al. 2015 (Malik et al. , 2016 . MERS-CoV PL pro was gradually heated from 20 to 94°C at 1°C/min increments and far-UV CD spectra were recorded from 200 to 250 nm. Far-UV CD spectra at several wavelengths were plotted as a function of temperature (Fig. 6a) . MERS-CoV PL pro underwent a single thermal transition as a function of temperature, suggesting two-state folding. Similar transition was also observed in the thermal shift assay using fluorescence spectroscopy. In an earlier study, secondary structure content was calculated and it was found that b-sheet structure (31%) is dominant in MERS-CoV PL pro (Lin et al. 2014) . Similarly, our results showed that MERS-CoV PL pro presented a single negative minimum at *218 nm suggesting a predominant b-sheet structure as well. When MERS-CoV PL pro was gradually heated, secondary structure was lost and became irregularly disorder structure (Fig. 6b) . The relative concentrations of folded and unfolded species as a function of temperature are shown in Fig. 6c . The thermal melting point (T m ) and van't Hoff enthalpy (DH VH ) of MERS-CoV PL pro calculated using Global 3 analysis software were 54.4 ± 0.1°C and 317.1 ± 3.9 kJ/mol, respectively. The thermal melting points (T m ) of MERS-CoV PL pro calculated by thermal shift assay and DMS were found to be in close proximity, reflecting tertiary-secondary structure unfolding events, respectively. In a recent study, papainlike protease of Murine Coronavirus also underwent a single thermal transition with a T m value of *46°C (Mielech et al. 2015) . A three-dimensional model of the thermal transitions of MERS-CoV PL pro was generated using Global 3 analysis software (Fig. 6d) . A single transition is clearly evident at lower wavelengths in the far-UV spectra region.

MERS-CoV PL pro was expressed in soluble state in E. coli and purified to homogeneity in two-step chromatography. Two orthogonal techniques were used for studying unfolding pathway that include the utilization of DMS and thermal shift assay. The results showed that MERS-CoV PL pro unfolds via a single thermal transition and follows a two-state unfolding pathway. Thermal shift assay calculated a T m value of 51.4°C and DMS method calculated a T m value of 54.4 ± 0.1°C, in agreement with sequential unfolding events of tertiary and secondary structures, respectively. Similar folding behavior and thermal melting point were also observed in papain-like protease of Murine 

Lipid-peptide vesicle nanoscale hybrids for triggered drug release by mild hyperthermia in vitro and in vivo

Travel implications of emerging coronaviruses: SARS and MERS-CoV

Middle East respiratory syndrome coronavirus: current situation and travelassociated concerns

Replication of human coronaviruses SARS-CoV, HCoV-NL63 and HCoV-229E is inhibited by the drug FK506

A case report of a middle east respiratory syndrome survivor with kidney biopsy results

SARS: prognosis, outcome and sequelae

Proteolytic processing and deubiquitinating activity of papain-like proteases of human coronavirus NL63

Deubiquitinating and interferon antagonism activities of coronavirus papain-like proteases

Coronaviruses: an RNA proofreading machine regulates replication fidelity and diversity

How aggregation and conformational scrambling of unfolded states govern fluorescence emission spectra

Intrinsic tryptophan fluorescence in the detection and analysis of proteins: a focus on Forster resonance energy transfer techniques

Using circular dichroism collected as a function of temperature to determine the thermodynamics of protein unfolding and binding interactions

A follow-up study of 69 discharged SARS patients

From SARS to MERS: crystallographic studies on coronaviral proteases enable antiviral drug design

van't Hoff enthalpies without baselines

Assessing activity and inhibition of Middle East respiratory syndrome coronavirus papain-like and 3C-like proteases using luciferasebased biosensors

Second derivative tryptophan fluorescence spectroscopy as a tool to characterize partially unfolded intermediates of proteins

Coronaviruses: emerging and re-emerging pathogens in humans and animals

Outbreak-related porcine epidemic diarrhea virus strains similar to US strains

Crystal structure of the papain-like protease of MERS coronavirus reveals unusual, potentially druggable active-site features

Hypocortisolism in survivors of severe acute respiratory syndrome (SARS)

Structural and functional characterization of MERS coronavirus papain-like protease

Co-circulation of four human coronaviruses (HCoVs) in Queensland children with acute respiratory tract illnesses in 2004

Spectroscopic and thermodynamic properties of recombinant heat shock protein A6 from Camelus dromedarius

Structural and thermodynamic properties of kappa class glutathione transferase from Camelus dromedarius

Murine coronavirus ubiquitin-like domain is important for papain-like protease stability and viral pathogenesis

Stability of equine lysozyme. I. Thermal unfolding behaviour

1-year pulmonary function and health status in survivors of severe acute respiratory syndrome

Unusual fluorescence of W168 in Plasmodium falciparum triosephosphate isomerase, probed by single-tryptophan mutants

Middle East respiratory syndrome-coronavirus (MERS-CoV): an updated overview and pharmacotherapeutics

Acute life threatening event (ALTE) in an infant with human coronavirus HCoV-229E infection

Coping with future epidemics: Tai chi practice as an overcoming strategy used by survivors of severe acute respiratory syndrome (SARS) in post-SARS Hong Kong

Heat-induced conformational changes of TET peptidase from crenarchaeon Desulfurococcus kamchatkensis

Epidemiology and vaccine of porcine epidemic diarrhea virus in China: a minireview

Distinct characteristics and complex evolution of PEDV strains

Middle East respiratory syndrome coronavirus

WHO (2016b) SARS (Severe Acute Respiratory Syndrome

Heat-induced unfolding of apo-CP43 studied by fluorescence spectroscopy and CD spectroscopy

Proteolytic processing, deubiquitinase and interferon antagonist activities of Middle East respiratory syndrome coronavirus papain-like protease

PLP2, a potent deubiquitinase from murine hepatitis virus, strongly inhibits cellular type I interferon production

Acknowledgements The authors extend their appreciation to the Deanship of Scientific Research at King Saud University for funding this work through the Research Project No R5-16-02-05.

Conflict of interest To the best of our knowledge, no conflict of interest, financial or others, exists. All authors are fully aware of this submission.