In silico identification of potential inhibitors against main protease of SARS-CoV-2 6LU7 from Andrographis panniculata via molecular docking,binding energy calculations and molecular dynamics simulation studies

BackgroundThe ongoing global outbreak of new corona virus (SARS-CoV-2) has been recognized as global public health concern since it causes high morbidity and mortality every day. Due to the rapid spreading and re-emerging, we need to find a potent drug against SARS-CoV-2. Synthetic drugs, such as hydroxychloroquine, remdisivir have paid more attention and the effects of these drugs are still under investigation, due to their severe side effects. Therefore, the aim of the present study was performed to identify the potential inhibitor against main protease SARS-CoV-2 6LU7.ObjectiveIn this study, RO5, ADME properties, molecular dynamic simulations and free binding energy prediction were mainly investigated.ResultsThe molecular docking study findings revealed that andrographolide had higher binding affinity among the selected natural diterpenoids compared to co-crystal native ligand inhibitor N3. The persistent inhibition of Ki for diterpenoids was analogous. Furthermore, the simulations of molecular dynamics and free binding energy findings have shown that andrographolide possesses a large amount of dynamic properties such as stability, flexibility and binding energy.ConclusionIn conclusion, findings of the current study suggest that selected diterpenoids were predicted to be the significant phytonutrient-based inhibitor against SARS-CoV-2 6LU7 (M^pro). However, preclinical and clinical trials are needed for the further scientific validation before use.     

Flexibility between the Protease and Helicase Domains of the Dengue Virus NS3 Protein Conferred by the Linker Region and Its Functional Implications

The dengue virus (DENV) NS3 protein is essential for viral polyprotein processing and RNA replication. It contains an N-terminal serine protease region (residues 1–168) joined to an RNA helicase (residues 180–618) by an 11-amino acid linker (169–179). The structure at 3.15 Å of the soluble NS3 protein from DENV4 covalently attached to 18 residues of the NS2B cofactor region (NS2B₁₈NS3) revealed an elongated molecule with the protease domain abutting subdomains I and II of the helicase (Luo, D., Xu, T., Hunke, C., Grüber, G., Vasudevan, S. G., and Lescar, J. (2008) J. Virol. 82, 173–183). Unexpectedly, using similar crystal growth conditions, we observed an alternative conformation where the protease domain has rotated by ∼161° with respect to the helicase domain. We report this new crystal structure bound to ADP-Mn²⁺ refined to a resolution of 2.2 Å. The biological significance for interdomain flexibility conferred by the linker region was probed by either inserting a Gly residue between Glu¹⁷³ and Pro¹⁷⁴ or replacing Pro¹⁷⁴ with a Gly residue. Both mutations resulted in significantly lower ATPase and helicase activities. We next increased flexibility in the linker by introducing a Pro¹⁷⁶ to Gly mutation in a DENV2 replicon system. A 70% reduction in luciferase reporter signal and a similar reduction in the level of viral RNA synthesis were observed. Our results indicate that the linker region has evolved to an optimum length to confer flexibility to the NS3 protein that is required both for polyprotein processing and RNA replication.     

Influence of proline on the thermostability of the active site and membrane arrangement of transmembrane proteins

Proline residues play a fundamental and subtle role in the dynamics, structure, and function in many membrane proteins. Temperature derivative spectroscopy and differential scanning calorimetry have been used to determine the effect of proline substitution in the structural stability of the active site and transmembrane arrangement of bacteriorhodopsin. We have analyzed the Pro-to-Ala mutation for the helix-embedded prolines Pro⁵⁰, Pro⁹¹, and Pro¹⁸⁶ in the native membrane environment. This information has been complemented with the analysis of the respective crystallographic structures by the FoldX force field. Differential scanning calorimetry allowed us to determine distorted membrane arrangement for P50A and P186A. The protein stability was severely affected for P186A and P91A. In the case of Pro⁹¹, a single point mutation is capable of strongly slowing down the conformational diffusion along the denaturation coordinate, becoming a barrier-free downhill process above 371 K. Temperature derivative spectroscopy, applied for first time to study thermal stability of proteins, has been used to monitor the stability of the active site of bacteriorhodopsin. The mutation of Pro⁹¹ and Pro¹⁸⁶ showed the most striking effects on the retinal binding pocket. These residues are the Pro in closer contact to the active site (activation energies for retinal release of 60.1 and 76.8 kcal/mol, respectively, compared to 115.8 kcal/mol for WT). FoldX analysis of the protein crystal structures indicates that the Pro-to-Ala mutations have both local and long-range effects on the structural stability of residues involved in the architecture of the protein and the active site and in the proton pumping function. Thus, this study provides a complete overview of the substitution effect of helix-embedded prolines in the thermodynamic and dynamic stability of a membrane protein, also related to its structure and function.     

D614G mutation in the SARS-CoV-2 spike protein enhances viral fitness by desensitizing it to temperature-dependent denaturation

The D614G mutation in the spike protein of SARS-CoV-2 alters the fitness of the virus, leading to the dominant form observed in the COVID-19 pandemic. However, the molecular basis of the mechanism by which this mutation enhances fitness is not clear. Here we demonstrated by cryo-electron microscopy that the D614G mutation resulted in increased propensity of multiple receptor-binding domains (RBDs) in an upward conformation poised for host receptor binding. Multiple substates within the one RBD-up or two RBD-up conformational space were determined. According to negative staining electron microscopy, differential scanning calorimetry, and differential scanning fluorimetry, the most significant impact of the mutation lies in its ability to eliminate the unusual cold-induced unfolding characteristics and to significantly increase the thermal stability under physiological pH. The D614G spike variant also exhibited exceptional long-term stability when stored at 37 °C for up to 2 months. Our findings shed light on how the D614G mutation enhances the infectivity of SARS-CoV-2 through a stabilizing mutation and suggest an approach for better design of spike protein-based conjugates for vaccine development.     

Investigation of changes in protein stability and substrate affinity of 3CL-protease of SARS-CoV-2 caused by mutations

3CL^pro of SARS-CoV-2 is one of the enzymes required for the replication process of the virus responsible for the COVID-19 pandemic. In this study, changes in protein stability and substrate affinity caused by mutations were investigated to stir the development of potent inhibitors. Sequence data of samples were obtained from the NCBI Virus database. Mutation analyses were performed with RDP4 and MegaX. 3CL^pro tertiary models were created using Robetta. Molecular docking for peptidomimetic substrate and inhibitor ligand was done with Autodock v4.2 and Haddock v2.4. Protein stability analysis was performed using mCSM stability and DynaMut2. Twenty-four missense mutations in 3CL^pro were identified in this study. Changes in the 3CL^pro structure induced by the mutations Met49Thr, Leu167Ser, and Val202Ala resulted in significant levels of instability (-2.029,-2.612,-2.177 kcal.mol^-1, respectively). The lowest interaction energy for substrate was -58.7 kcal.mol^-1 and -62.6 kcal.mol^-1 in wild-type and mutant, respectively. The lowest docking energy for ligand was -6.19 and -9.52 kcal.mol^-1 for wild-type and mutant, respectively. This study reports for the first time that mutations cause increased substrate affinity of 3CL^pro from SARS-CoV-2. This research provides important data for the development of potent peptidomimetic inhibitors for the treatment of COVID-19. Keywords: 3CL-protease, mutation analysis, protein stability, SARS-CoV-2 genome, substrate affinity     

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203–205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral ‘RG’ motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2’s continued adaptation to human infection.     

Secondary mutations M36I and A71V in the human immunodeficiency virus type 1 protease can provide an advantage for the emergence of the primary mutation D30N

Development of resistance mutations in enzymatic targets of human immunodeficiency virus 1 (HIV-1) hampers the ability to provide adequate therapy. Of special interest is the effect mutations outside the active site of HIV-1 protease have on inhibitor binding and virus viability. We engineered protease mutants containing the active site mutation D30N alone and with the nonactive site polymorphisms M36I and/or A71V. We determined the K(i) values for the inhibitors nelfinavir, ritonavir, indinavir, KNI272, and AG1776 as well as the catalytic efficiency of the mutants. Single and double mutation combinations exhibited a decrease in catalytic efficiency, while the triple mutant displayed catalytic efficiency greater than that of the wild type. Variants containing M36I or A71V alone did not display a significant change in binding affinities to the inhibitors tested. The variant containing mutation D30N displayed a 2-6-fold increase in K(i) for all inhibitors tested, with nelfinavir showing the greatest increase. The double mutants containing a combination of mutations D30N, M36I, and A71V displayed -0.5-fold to +6-fold changes in the K(i) of all inhibitors tested, with ritonavir and nelfinavir most affected. Only the triple mutant showed a significant increase (>10-fold) in K(i) for inhibitor nelfinavir, ritonavir, or AG-1776 displaying 22-, 19-, or 15-fold increases, respectively. Our study shows that the M36I and A71V mutations provide a greater level of inhibitor cross-resistance combined with active site mutation D30N. M36I and A71V, when present as natural polymorphisms, could aid the virus in developing active site mutations to escape inhibitor binding while maintaining catalytic efficiency.     

Selection of Drug-Resistant Feline Immunodeficiency Virus (FIV) Encoding FIV/HIV Chimeric Protease in the Presence of HIV-Specific Protease Inhibitors

An infectious chimeric feline immunodeficiency virus (FIV)/HIV strain carrying six HIV-like protease (PR) mutations (I37V/N55M/V59I/I98S/Q99V/P100N) was subjected to selection in culture against the PR inhibitor lopinavir (LPV), darunavir (DRV), or TL-3. LPV selection resulted in the sequential emergence of V99A (strain S-1X), I59V (strain S-2X), and I108V (strain S-3X) mutations, followed by V37I (strain S-4X). Mutant PRs were analyzed in vitro, and an isogenic virus producing each mutant PR was analyzed in culture for LPV sensitivity, yielding results consistent with the original selection. The 50% inhibitory concentrations (IC₅₀s) for S-1X, S-2X, S-3X, and S-4X were 95, 643, 627, and 1,543 nM, respectively. The primary resistance mutations, V99⁸²A, I59⁵⁰V, and V37³²I, are consistent with the resistance pattern developed by HIV-1 under similar selection conditions. While resistance to LPV emerged readily, similar PR mutations causing resistance to either DRV or TL-3 failed to emerge after passage for more than a year. However, a G37D mutation in the nucleocapsid (NC) was observed in both selections and an isogenic G37D mutant replicated in the presence of 100 nM DRV or TL-3, whereas parental chimeric FIV could not. An additional mutation, L92V, near the PR active site in the folded structure recently emerged during TL-3 selection. The L92V mutant PR exhibited an IC₅₀ of 50 nM, compared to 35 nM for 6s-98S PR, and processed the NC-p2 junction more efficiently, consistent with increased viral fitness. These findings emphasize the role of mutations outside the active site of PR in increasing viral resistance to active-site inhibitors and suggest additional targets for inhibitor development.     

A computational method to characterize the missense mutations in the catalytic domain of GAA protein causing Pompe disease

Pompe disease is an autosomal recessive lysosomal storage disease caused by acid α-glucosidase (GAA) deficiency, resulting in intralysosomal accumulation of glycogen, including cardiac, skeletal, and smooth muscle cells. The GAA gene is located on chromosome 17 (17q25.3), the GAA protein consists of 952 amino acids; of which 378 amino acids (347-726) falls within the catalytic domain of the protein and comprises of active sites (518 and 521) and binding sites (404, 600, 616, and 674). In this study, we used several computational tools to classify the missense mutations in the catalytic domain of GAA for their pathogenicity and stability. Eight missense mutations (R437C, G478R, N573H, Y575S, G605D, V642D, L705P, and L712P) were predicted to be pathogenic and destabilizing to the protein structure. These mutations were further subjected to phenotyping analysis using SNPeffect 4.0 to predict the chaperone binding sites and structural stability of the protein. The mutations R437C and G478R were found to compromise the chaperone-binding activity with GAA. Molecular docking analysis revealed that the G478R mutation to be more significant and hinders binding to the DNJ (Miglustat) compared with the R437C. Further molecular dynamic analysis for the two mutations demonstrated that the G478R mutation was acquired higher deviation, fluctuation, and lower compactness with decreased intramolecular hydrogen bonds compared to the mutant R437C. These data are expected to serve as a platform for drug design against Pompe disease and will serve as an ultimate tool for variant classification and interpretations.     

The SARS-CoV-2 and other human coronavirus spike proteins are fine-tuned towards temperature and proteases of the human airways

The high transmissibility of SARS-CoV-2 is related to abundant replication in the upper airways, which is not observed for the other highly pathogenic coronaviruses SARS-CoV and MERS-CoV. We here reveal features of the coronavirus spike (S) protein, which optimize the virus towards the human respiratory tract. First, the S proteins exhibit an intrinsic temperature preference, corresponding with the temperature of the upper or lower airways. Pseudoviruses bearing the SARS-CoV-2 spike (SARS-2-S) were more infectious when produced at 33°C instead of 37°C, a property shared with the S protein of HCoV-229E, a common cold coronavirus. In contrast, the S proteins of SARS-CoV and MERS-CoV favored 37°C, in accordance with virus preference for the lower airways. Next, SARS-2-S-driven entry was efficiently activated by not only TMPRSS2, but also the TMPRSS13 protease, thus broadening the cell tropism of SARS-CoV-2. Both proteases proved relevant in the context of authentic virus replication. TMPRSS13 appeared an effective spike activator for the virulent coronaviruses but not the low pathogenic HCoV-229E virus. Activation of SARS-2-S by these surface proteases requires processing of the S1/S2 cleavage loop, in which both the furin recognition motif and extended loop length proved critical. Conversely, entry of loop deletion mutants is significantly increased in cathepsin-rich cells. Finally, we demonstrate that the D614G mutation increases SARS-CoV-2 stability, particularly at 37°C, and, enhances its use of the cathepsin L pathway. This indicates a link between S protein stability and usage of this alternative route for virus entry. Since these spike properties may promote virus spread, they potentially explain why the spike-G614 variant has replaced the early D614 variant to become globally predominant. Collectively, our findings reveal adaptive mechanisms whereby the coronavirus spike protein is adjusted to match the temperature and protease conditions of the airways, to enhance virus transmission and pathology.     

Dynamic allostery highlights the evolutionary differences between the CoV-1 and CoV-2 main proteases

The SARS-CoV-2 coronavirus has become one of the most immediate and widely studied systems since its identification and subsequent global outbreak from 2019 to 2022. In an effort to understand the biophysical changes as a result of mutations, the mechanistic details of multiple different proteins within the SARS-CoV-2 virus have been studied and compared with SARS-CoV-1. Focusing on the main protease (mPro), we explored the long-range dynamics using the Dynamic Coupling Index (DCI) to investigate the dynamic coupling between the catalytic site residues and the rest of the protein, both inter- and intrachain, for the CoV-1 and CoV-2 mPro. We found that there is significant cross-chain coupling between these active sites and specific distal residues in the CoV-2 mPro not present in CoV-1. The enhanced long-distance interactions, particularly between the two chains, suggest subsequently enhanced cooperativity for CoV-2. A further comparative analysis of the dynamic flexibility using the dynamic flexibility index (DFI) between the CoV-1 and CoV-2 mPros shows that the inhibitor binding near active sites induces change in flexibility to a distal region of the protein, opposite in behavior between the two systems; this region becomes more flexible upon inhibitor binding in CoV-1, while it becomes less flexible in the CoV-2 mPro. Upon inspection, we show that, on average, the dynamic flexibility of the sites substituted from CoV-1 to CoV-2 changes significantly less than the average calculated across all residues within the structure, indicating that the differences in behaviors between the two systems is likely the result of allosteric influence, in which the new substitutions in CoV-2 induce flexibility and dynamic changes elsewhere in the structure.     

Human Immunodeficiency Virus Type 1 Protease Genotypes and In Vitro Protease Inhibitor Susceptibilities of Isolates from Individuals Who Were Switched to Other Protease Inhibitors after Long-Term Saquinavir Treatment

An understanding of the mechanisms of virologic cross-resistance between human immunodeficiency virus type 1 protease inhibitors is important for the establishment of effective treatment strategies for patients who no longer respond to their initial protease inhibitor. Protease gene sequencing results from patients treated with saquinavir showed significant increases in the frequency of the G48V protease mutation in patients receiving higher doses of the drug. In addition, all six patients who developed the G48V mutation during saquinavir therapy developed the V82A mutation either on continued saquinavir or after a switch to nelfinavir or indinavir. In vitro susceptibility assays showed that all 13 isolates with reduced susceptibilities to two or more protease inhibitors had either the G48V or L90M mutation, along with an average of six other protease mutations. Reduced susceptibility to nelfinavir was found in 14 isolates, but only 1 possessed the D30N mutation. These results suggest that mutations selected in vivo by initial saquinavir therapy may provide more cross-resistance to the other protease inhibitors than has been previously reported.     

Structural and kinetic analyses of the protease from an amprenavir-resistant human immunodeficiency virus type 1 mutant rendered resistant to saquinavir and resensitized to amprenavir

Recent drug regimens have had much success in the treatment of human immunodeficiency virus (HIV)-infected individuals; however, the incidence of resistance to such drugs has become a problem that is likely to increase in importance with long-term therapy of this chronic illness. An analysis and understanding of the molecular interactions between the drug(s) and the mutated viral target(s) is crucial for further progress in the field of AIDS therapy. The protease inhibitor amprenavir (APV) generates a signature set of HIV type 1 (HIV-1) protease mutations associated with in vitro resistance (M46I/L, I47V, and I50V [triple mutant]). Passage of the triple-mutant APV-resistant HIV-1 strain in MT4 cells, in the presence of increasing concentrations of saquinavir (SQV), gave rise to a new variant containing M46I, G48V, I50V, and I84L mutations in the protease and a resulting phenotype that was resistant to SQV and, unexpectedly, resensitized to APV. This phenotype was consistent with a subsequent kinetic analysis of the mutant protease, together with X-ray crystallographic analysis and computational modeling which elucidated the structural basis of these observations. The switch in protease inhibitor sensitivities resulted from (i) the I50V mutation, which reduced the area of contact with APV and SQV; (ii) the compensating I84L mutation, which improved hydrophobic packing with APV; and (iii) the G-to-V mutation at residue 48, which introduced steric repulsion with the P3 group of SQV. This analysis establishes the fine detail necessary for understanding the loss of protease binding for SQV in the quadruple mutant and gain in binding for APV, demonstrating the powerful combination of virology, molecular biology, enzymology, and protein structural and modeling studies in the elucidation and understanding of viral drug resistance.     

Crystal Structure of the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Papain-like Protease Bound to Ubiquitin Facilitates Targeted Disruption of Deubiquitinating Activity to Demonstrate Its Role in Innate Immune Suppression

Middle East respiratory syndrome coronavirus (MERS-CoV) is a newly emerging human pathogen that was first isolated in 2012. MERS-CoV replication depends in part on a virus-encoded papain-like protease (PL^pro) that cleaves the viral replicase polyproteins at three sites releasing non-structural protein 1 (nsp1), nsp2, and nsp3. In addition to this replicative function, MERS-CoV PL^pro was recently shown to be a deubiquitinating enzyme (DUB) and to possess deISGylating activity, as previously reported for other coronaviral PL^pro domains, including that of severe acute respiratory syndrome coronavirus. These activities have been suggested to suppress host antiviral responses during infection. To understand the molecular basis for ubiquitin (Ub) recognition and deconjugation by MERS-CoV PL^pro, we determined its crystal structure in complex with Ub. Guided by this structure, mutations were introduced into PL^pro to specifically disrupt Ub binding without affecting viral polyprotein cleavage, as determined using an in trans nsp3↓4 cleavage assay. Having developed a strategy to selectively disable PL^pro DUB activity, we were able to specifically examine the effects of this activity on the innate immune response. Whereas the wild-type PL^pro domain was found to suppress IFN-β promoter activation, PL^pro variants specifically lacking DUB activity were no longer able to do so. These findings directly implicate the DUB function of PL^pro, and not its proteolytic activity per se, in the inhibition of IFN-β promoter activity. The ability to decouple the DUB activity of PL^pro from its role in viral polyprotein processing now provides an approach to further dissect the role(s) of PL^pro as a viral DUB during MERS-CoV infection.     

RyR2 disease mutations at the C-terminal domain intersubunit interface alter closed-state stability and channel activation

Ryanodine receptors (RyRs) are ion channels that mediate the release of Ca²⁺ from the sarcoplasmic reticulum/endoplasmic reticulum, mutations of which are implicated in a number of human diseases. The adjacent C-terminal domains (CTDs) of cardiac RyR (RyR2) interact with each other to form a ring-like tetrameric structure with the intersubunit interface undergoing dynamic changes during channel gating. This mobile CTD intersubunit interface harbors many disease-associated mutations. However, the mechanisms of action of these mutations and the role of CTD in channel function are not well understood. Here, we assessed the impact of CTD disease-associated mutations P4902S, P4902L, E4950K, and G4955E on Ca²⁺− and caffeine-mediated activation of RyR2. The G4955E mutation dramatically increased both the Ca²⁺-independent basal activity and Ca²⁺-dependent activation of [³H]ryanodine binding to RyR2. The P4902S and E4950K mutations also increased Ca²⁺ activation but had no effect on the basal activity of RyR2. All four disease mutations increased caffeine-mediated activation of RyR2 and reduced the threshold for activation and termination of spontaneous Ca²⁺ release. G4955D dramatically increased the basal activity of RyR2, whereas G4955K mutation markedly suppressed channel activity. Similarly, substitution of P4902 with a negatively charged residue (P4902D), but not a positively charged residue (P4902K), also dramatically increased the basal activity of RyR2. These data suggest that electrostatic interactions are involved in stabilizing the CTD intersubunit interface and that the G4955E disease mutation disrupts this interface, and thus the stability of the closed state. Our studies shed new insights into the mechanisms of action of RyR2 CTD disease mutations.     

Ca2+-dependent mechanism of membrane insertion and destabilization by the SARS-CoV-2 fusion peptide

Cell penetration after recognition of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus by the ACE2 receptor and the fusion of its viral envelope membrane with cellular membranes are the early steps of infectivity. A region of the Spike protein of the virus, identified as the “fusion peptide” (FP), is liberated at its N-terminal site by a specific cleavage occurring in concert with the interaction of the receptor-binding domain of the Spike. Studies have shown that penetration is enhanced by the required binding of Ca²⁺ ions to the FPs of coronaviruses, but the mechanisms of membrane insertion and destabilization remain unclear. We have predicted the preferred positions of Ca²⁺ binding to the SARS-CoV-2-FP, the role of Ca²⁺ ions in mediating peptide-membrane interactions, the preferred mode of insertion of the Ca²⁺-bound SARS-CoV-2-FP, and consequent effects on the lipid bilayer from extensive atomistic molecular dynamics simulations and trajectory analyses. In a systematic sampling of the interactions of the Ca²⁺-bound peptide models with lipid membranes, SARS-CoV-2-FP penetrated the bilayer and disrupted its organization only in two modes involving different structural domains. In one, the hydrophobic residues F833/I834 from the middle region of the peptide are inserted. In the other, more prevalent mode, the penetration involves residues L822/F823 from the LLF motif, which is conserved in CoV-2-like viruses, and is achieved by the binding of Ca²⁺ ions to the D830/D839 and E819/D820 residue pairs. FP penetration is shown to modify the molecular organization in specific areas of the bilayer, and the extent of membrane binding of the SARS-CoV-2 FP is significantly reduced in the absence of Ca²⁺ ions. These findings provide novel mechanistic insights regarding the role of Ca²⁺ in mediating SARS-CoV-2 fusion and provide a detailed structural platform to aid the ongoing efforts in rational design of compounds to inhibit SARS-CoV-2 cell entry.     

Structural and functional characterization of MERS coronavirus papain-like protease

Backgrounds

A new highly pathogenic human coronavirus (CoV), Middle East respiratory syndrome coronavirus (MERS-CoV), has emerged in Jeddah and Saudi Arabia and quickly spread to some European countries since September 2012. Until 15 May 2014, it has infected at least 572 people with a fatality rate of about 30% globally. Studies to understand the virus and to develop antiviral drugs or therapy are necessary and urgent. In the present study, MERS-CoV papain-like protease (PL^pro) is expressed, and its structural and functional consequences are elucidated.

Results

Circular dichroism and Tyr/Trp fluorescence analyses indicated that the secondary and tertiary structure of MERS-CoV PL^pro is well organized and folded. Analytical ultracentrifugation analyses demonstrated that MERS-CoV PL^pro is a monomer in solution. The steady-state kinetic and deubiquitination activity assays indicated that MERS-CoV PL^pro exhibits potent deubiquitination activity but lower proteolytic activity, compared with SARS-CoV PL^pro. A natural mutation, Leu105, is the major reason for this difference.

Conclusions

Overall, MERS-CoV PL^pro bound by an endogenous metal ion shows a folded structure and potent proteolytic and deubiquitination activity. These findings provide important insights into the structural and functional properties of coronaviral PL^pro family, which is applicable to develop strategies inhibiting PL^pro against highly pathogenic coronaviruses.     

N-Terminal Finger Stabilizes the S1 Pocket for the Reversible Feline Drug GC376 in the SARS-CoV-2 Mpro Dimer

The main protease (M^pro, also known as 3CL protease) of SARS-CoV-2 is a high priority drug target in the development of antivirals to combat COVID-19 infections. A feline coronavirus antiviral drug, GC376, has been shown to be effective in inhibiting the SARS-CoV-2 main protease and live virus growth. As this drug moves into clinical trials, further characterization of GC376 with the main protease of coronaviruses is required to gain insight into the drug’s properties, such as reversibility and broad specificity. Reversibility is an important factor for therapeutic proteolytic inhibitors to prevent toxicity due to off-target effects. Here we demonstrate that GC376 has nanomolar K_i values with the M^pro from both SARS-CoV-2 and SARS-CoV strains. Restoring enzymatic activity after inhibition by GC376 demonstrates reversible binding with both proteases. In addition, the stability and thermodynamic parameters of both proteases were studied to shed light on physical chemical properties of these viral enzymes, revealing higher stability for SARS-CoV-2 M^pro. The comparison of a new X-ray crystal structure of M^pro from SARS-CoV complexed with GC376 reveals similar molecular mechanism of inhibition compared to SARS-CoV-2 M^pro, and gives insight into the broad specificity properties of this drug. In both structures, we observe domain swapping of the N-termini in the dimer of the M^pro, which facilitates coordination of the drug’s P1 position. These results validate that GC376 is a drug with an off-rate suitable for clinical trials.     

Engineered disulfide reveals structural dynamics of locked SARS-CoV-2 spike

The spike (S) protein of SARS-CoV-2 has been observed in three distinct pre-fusion conformations: locked, closed and open. Of these, the function of the locked conformation remains poorly understood. Here we engineered a SARS-CoV-2 S protein construct “S-R/x3” to arrest SARS-CoV-2 spikes in the locked conformation by a disulfide bond. Using this construct we determined high-resolution structures confirming that the x3 disulfide bond has the ability to stabilize the otherwise transient locked conformations. Structural analyses reveal that wild-type SARS-CoV-2 spike can adopt two distinct locked-1 and locked-2 conformations. For the D614G spike, based on which all variants of concern were evolved, only the locked-2 conformation was observed. Analysis of the structures suggests that rigidified domain D in the locked conformations interacts with the hinge to domain C and thereby restrains RBD movement. Structural change in domain D correlates with spike conformational change. We propose that the locked-1 and locked-2 conformations of S are present in the acidic high-lipid cellular compartments during virus assembly and egress. In this model, release of the virion into the neutral pH extracellular space would favour transition to the closed or open conformations. The dynamics of this transition can be altered by mutations that modulate domain D structure, as is the case for the D614G mutation, leading to changes in viral fitness. The S-R/x3 construct provides a tool for the further structural and functional characterization of the locked conformations of S, as well as how sequence changes might alter S assembly and regulation of receptor binding domain dynamics.     

新型冠状病毒抗体和小分子抑制剂的研究现状

世界卫生组织已宣布新型冠状病毒感染（coronavirus disease 2019,COVID-19）的爆发为全球大流行。中和抗体和小分子抑制剂在预防及治疗COVID-19中发挥重要作用。尽管已开发出了多种中和抗体以及疫苗,但是随着病原体严重急性呼吸综合征冠状病毒2(severe acute respiratory syndrome coronavirus 2,SARS-CoV-2)的不断变异,现有的抗体及疫苗面临巨大的挑战。小分子抑制剂主要通过干扰病毒与宿主的结合以及病毒自身的复制达到消灭病毒以及抑制病毒感染的作用,并且对SARS-CoV-2突变株具有广谱抑制作用,是当前研究的热点。近年来国内外学者对SARS-CoV-2的小分子抑制剂做了大量的研究工作,本文根据中和抗体识别的抗原表位以及小分子抑制剂的作用机制分别对用于预防及治疗COVID-19的中和抗体和小分子抑制剂进行综述,讨论其研究现状,并展望小分子抑制剂的应用前景,以期为该领域的进一步研究提供参考。