S1 Subsite Research Articles

Identification of potential inhibitors of cholinergic and β-secretase enzymes from phytochemicals derived from Gongronema latifolium Benth leaf: an integrated computational analysis.

Neurodegenerative disorders (NDDs) are associated with increased activities of the brain acetylcholinesterase (AChE), butyrylcholinesterase (BChE) and β-secretase enzyme (BACE1). Inhibition of these enzymes affords therapeutic option for managing NDDs such as Alzheimer's disease (AD) and Parkinson's disease (PD). Although, Gongronema latifolium Benth (GL) has been widely documented in ethnopharmacological and scientific reports for the management of NDDs, there is paucity of information on its underlying mechanism and neurotherapeutic constituents. Herein, 152 previously reported Gongronema latifolium derived-phytochemicals (GLDP) were screened against hAChE, hBChE and hBACE-1 using molecular docking, molecular dynamics (MD) simulations, free energy of binding calculations and cluster analysis. The result of the computational analysis identified silymarin, alpha-amyrin and teraxeron with the highest binding energies (-12.3, -11.2, -10.5 Kcal/mol) for hAChE, hBChE and hBACE-1 respectively as compared with those of the reference inhibitors (-12.3, -9.8 and - 9.4 for donepezil, propidium and aminoquinoline compound respectively). These best docked phytochemicals were found to be orientated in the hydrophobic gorge where they interacted with the choline-binding pocket in the A-site and P-site of the cholinesterase and subsites S1, S3, S3' and flip (67-75) residues of the pocket of the BACE-1. The best docked phytochemicals complexed with the target proteins were stable in a 100 ns molecular dynamic simulation. The interactions with the catalytic residues were preserved during the simulation as observed from the MMGBSA decomposition and cluster analyses. The presence of these phytocompounds most notably silymarin, which demonstrated dual high binding tendencies to both cholinesterases, were identified as potential neurotherapeutics subject to further investigation.

Discovery of Chalcone-Based Hybrid Structures as High Affinity and Site-Specific Inhibitors against SARS-CoV-2: A Comprehensive Structural Analysis Based on Various Host-Based and Viral Targets.

Previous studies indicated that natural-based chalcones have significant inhibitory effects on the coronavirus enzymes 3CLpro and PLpro as well as modulation of some host-based antiviral targets (HBATs). In this study, a comprehensive computational and structural study was performed to investigate the affinity of our compound library consisting of 757 chalcone-based structures (CHA-1 to CHA-757) for inhibiting the 3CLpro and PLpro enzymes and against twelve selected host-based targets. Our results indicated that CHA-12 (VUF 4819) is the most potent and multi-target inhibitor in our chemical library over all viral and host-based targets. Correspondingly, CHA-384 and its congeners containing ureide moieties were found to be potent and selective 3CLpro inhibitors, and benzotriazole moiety in CHA-37 was found to be a main fragment for inhibiting the 3CLpro and PLpro. Surprisingly, our results indicate that the ureide and sulfonamide moieties are integral fragments for the optimum 3CLpro inhibition while occupying the S1 and S3 subsites, which is fully consistent with recent reports on the site-specific 3CLpro inhibitors. Finding the multi-target inhibitor CHA-12, previously reported as an LTD4 antagonist for the treatment of inflammatory pulmonary diseases, prompted us to suggest it as a concomitant agent for relieving respiratory symptoms and suppressing COVID-19 infection.

Structure-activity relationships of dihydropyrimidone inhibitors against native and auto-processed human neutrophil elastase

BackgroundHuman neutrophil elastase (HNE) is a key driver of systemic and cardiopulmonary inflammation. Recent studies have established the existence of a pathologically active auto-processed form of HNE with reduced binding affinity against small molecule inhibitors. MethodAutoDock Vina v1.2.0 and Cresset Forge v10 software were used to develop a 3D-QSAR model for a series of 47 DHPI inhibitors. Molecular Dynamics (MD) simulations were carried out using AMBER v18 to study the structure and dynamics of sc (single-chain HNE) and tcHNE (two-chain HNE). MMPBSA binding free energies of the previously reported clinical candidate BAY 85–8501 and the highly active BAY-8040 were calculated with sc and tcHNE. ResultsThe DHPI inhibitors occupy the S1 and S2 subsites of scHNE. The robust 3D-QSAR model showed acceptable predictive and descriptive capability with regression coefficient of r2 = 0.995 and cross-validation regression coefficient q2 = 0.579 for the training set. The key descriptors of shape, hydrophobics and electrostatics were mapped to the inhibitory activity. In auto-processed tcHNE, the S1 subsite undergoes widening and disruption. All the DHPI inhibitors docked with the broadened S1′-S2′ subsites of tcHNE with lower AutoDock binding affinities. The MMPBSA binding free energy of BAY-8040 with tcHNE reduced in comparison with scHNE while the clinical candidate BAY 85–8501 dissociated during MD. Thus, BAY-8040 may have lower inhibitory activity against tcHNE whereas the clinical candidate BAY 85–8501 is likely to be inactive. ConclusionSAR insights gained from this study will aid the future development of inhibitors active against both forms of HNE.

Exploration of imatinib and nilotinib-derived templates as the P2-Ligand for HIV-1 protease inhibitors: Design, synthesis, protein X-ray structural studies, and biological evaluation

Structure-based design, synthesis, X-ray structural studies, and biological evaluation of a new series of potent HIV-1 protease inhibitors are described. These inhibitors contain various pyridyl-pyrimidine, aryl thiazole or alkylthiazole derivatives as the P2 ligands in combination with darunavir-like hydroxyethylamine sulfonamide isosteres. These heterocyclic ligands are inherent to kinase inhibitor drugs, such as nilotinib and imatinib. These ligands are designed to make hydrogen bonding interactions with the backbone atoms in the S2 subsite of HIV-1 protease. Various benzoic acid derivatives have been synthesized and incorporation of these ligands provided potent inhibitors that exhibited subnanomolar level protease inhibitory activity and low nanomolar level antiviral activity. Two high resolution X-ray structures of inhibitor-bound HIV-1 protease were determined. These structures provided important ligand-binding site interactions for further optimization of this class of protease inhibitors.

Structure-guided design of direct-acting antivirals that exploit the gem-dimethyl effect and potently inhibit 3CL proteases of severe acute respiratory syndrome Coronavirus-2 (SARS-CoV-2) and middle east respiratory syndrome coronavirus (MERS-CoV)

The high morbidity and mortality associated with SARS-CoV-2 infection, the etiological agent of COVID-19, has had a major impact on global public health. Significant progress has been made in the development of an array of vaccines and biologics, however, the emergence of SARS-CoV-2 variants and breakthrough infections are an ongoing major concern. Furthermore, there is an existing paucity of small-molecule host and virus-directed therapeutics and prophylactics that can be used to counter the spread of SARS-CoV-2, and any emerging and re-emerging coronaviruses. We describe herein our efforts to address this urgent need by focusing on the structure-guided design of potent broad-spectrum inhibitors of SARS-CoV-2 3C-like protease (3CLpro or Main protease), an enzyme essential for viral replication. The inhibitors exploit the directional effects associated with the presence of a gem-dimethyl group that allow the inhibitors to optimally interact with the S4 subsite of the enzyme. Several compounds were found to potently inhibit SARS-CoV-2 and MERS-CoV 3CL proteases in biochemical and cell-based assays. Specifically, the EC50 values of aldehyde 1c and its corresponding bisulfite adduct 1d against SARS-CoV-2 were found to be 12 and 10 nM, respectively, and their CC50 values were >50 μM. Furthermore, deuteration of these compounds yielded compounds 2c/2d with EC50 values 11 and 12 nM, respectively. Replacement of the aldehyde warhead with a nitrile (CN) or an α-ketoamide warhead or its corresponding bisulfite adduct yielded compounds 1g, 1eand1f with EC50 values 60, 50 and 70 nM, respectively. High-resolution cocrystal structures have identified the structural determinants associated with the binding of the inhibitors to the active site of the enzyme and, furthermore, have illuminated the mechanism of action of the inhibitors. Overall, the high Safety Index (SI) (SI=CC50/EC50) displayed by these compounds suggests that they are well-suited to conducting further preclinical studies.

Computational insights for predicting the binding and selectivity of peptidomimetic plasmepsin IV inhibitors against cathepsin D.

Plasmepsins (Plms) are aspartic proteases involved in the degradation of human hemoglobin by P. falciparum and are essential for the survival and growth of the parasite. Therefore, Plm enzymes are reported as an important antimalarial drug target. Herein, we have applied molecular docking, molecular dynamics (MD) simulations, and binding free energy with the Linear Interaction Energy (LIE) approach to investigate the binding of peptidomimetic PlmIV inhibitors with a particular focus on understanding their selectivity against the human Asp protease cathepsin D (CatD). The residual decomposition analysis results suggest that amino acid differences in the subsite S3 of PlmIV and CatD are responsible for the higher selectivity of the 5a inhibitor. These findings yield excellent agreement with experimental binding data and provide new details regarding van der Waals and electrostatic interactions of subsite residues as well as structural properties of the PlmIV and CatD systems.

Discovery of novel thioquinazoline-N-aryl-acetamide/N-arylacetohydrazide hybrids as anti-SARS-CoV-2 agents: Synthesis, in vitro biological evaluation, and molecular docking studies

In the current investigation, two novel series of (tetrahydro)thioquinazoline-N-arylacetamides and (tetrahydro)thioquinazoline-N-arylacetohydrazides were designed, synthesized and investigated for their antiviral activity against SARS-CoV-2. The thioquinazoline-N-arylacetamide 17g as well as the tetrahydrothioquinazoline-N-arylacetohydrazides 18c and 18f showed potent antiviral activity with IC50 of 21.4, 38.45 and 26.4 µM, respectively. In addition, 18c and 18f demonstrated potential selectivity toward the SARS-CoV-2 over the host cells with SI of 10.67 and 16.04, respectively. Further evaluation of the mechanism of action of the three derivatives 17g, 18c, and 18f displayed that they can inhibit the virus at the adsorption as well as at the replication stages, in addition to their virucidal properties. In addition, 17g, 18c, and 18f demonstrated satisfactory physicochemical properties as well as drug-likeness properties to be further optimized for the discovery of novel antiviral agents. The docking simulation on Mpro binding site predicted the binding pattern of the target compounds rationalizing their differential activity based on their hydrophobic interaction and fitting in the hydrophobic S2 subsite of the binding site

Peptide inhibitors of angiotensin-I converting enzyme based on angiotensin (1–7) with selectivity for the C-terminal domain

The renin-angiotensin system (RAS) is a key regulator of human arterial pressure. Several of its effects are modulated by angiotensin II, an octapeptide originating from the action of angiotensin-I converting enzyme (ACE) on the decapeptide angiotensin-I. ACE possess two active sites (nACE and cACE) that have their own kinetic and substrate specificities. ACE inhibitors are widely used as the first-line treatment for hypertension and other heart-related diseases, but because they inactivate both ACE domains, their use is associated with serious side effects. Thus, the search for domain-specific ACE inhibitors has been the focus of intense research. Angiotensin (1–7), a peptide that also belongs to the RAS, acts as a substrate of nACE and an inhibitor of cACE. We have synthetized 15 derivatives of Ang (1–7), sequentially removing the N-terminal amino acids and modifying peptides extremities, to find molecules with improved selectivity and inhibition properties. Ac-Ang (2–7)-NH2 is a good ACE inhibitor, resistant to cleavage and with improved cACE selectivity. Molecular dynamics simulations provided a model for this peptide’s selectivity, due to Val3 and Tyr4 interactions with ACE subsites. Val3 has an important interaction with the S3 subsite, since its removal greatly reduced peptide-enzyme interactions. Taken together, our findings support ongoing studies using insights from the binding of Ac-Ang (2–7)-NH2 to develop effective cACE inhibitors.

Novel covalent and non-covalent complex-based pharmacophore models of SARS-CoV-2 main protease (Mpro) elucidated by microsecond MD simulations

As the world enters its second year of the pandemic caused by SARS-CoV-2, intense efforts have been directed to develop an effective diagnosis, prevention, and treatment strategies. One promising drug target to design COVID-19 treatments is the SARS-CoV-2 Mpro. To date, a comparative understanding of Mpro dynamic stereoelectronic interactions with either covalent or non-covalent inhibitors (depending on their interaction with a pocket called S1’ or oxyanion hole) has not been still achieved. In this study, we seek to fill this knowledge gap using a cascade in silico protocol of docking, molecular dynamics simulations, and MM/PBSA in order to elucidate pharmacophore models for both types of inhibitors. After docking and MD analysis, a set of complex-based pharmacophore models was elucidated for covalent and non-covalent categories making use of the residue bonding point feature. The highest ranked models exhibited ROC-AUC values of 0.93 and 0.73, respectively for each category. Interestingly, we observed that the active site region of Mpro protein–ligand complex undergoes large conformational changes, especially within the S2 and S4 subsites. The results reported in this article may be helpful in virtual screening (VS) campaigns to guide the design and discovery of novel small-molecule therapeutic agents against SARS-CoV-2 Mpro protein.

Structure-based design of anti-mycobacterial drug leads that target the mycolic acid transporter MmpL3

New anti-tubercular agents are urgently needed to address the emerging threat of drug resistance to human tuberculosis. Here, we have used structure-assisted methods to develop compounds that target mycobacterial membrane protein large 3 (MmpL3). MmpL3 is essential for the transport of mycolic acids, an important cell-wall component of mycobacteria. We prepared compounds that potently inhibit the growth of Mycobacterium tuberculosis (Mtb) and other mycobacteria in cell culture. The cryoelectron microscopy (cryo-EM) structure of mycobacterial MmpL3 in complex with one of these compounds (ST004) was determined using lipid nanodiscs at an overall resolution of 3.36Å. The structure reveals the binding mode of ST004 to MmpL3, with the S4 and S5 subsites of the inhibitor-binding pocket in the proton translocation channel playing vital roles. These data are a promising starting point for the development of anti-tuberculosis drugs that target MmpL3.

Exploring structural features of potent dipeptidyl peptidase IV (DPP-IV) inhibitory peptides derived from tilapia (Oreochromis niloticus) skin gelatin by an integrated approach of multivariate analysis and Gly-Pro-based peptide library

There are abundant dipeptidyl-peptidase IV (DPP-IV) inhibitory peptide motifs in collagen sequences due to high content of Pro. However, the structural features of the most potent DPP-IV inhibitory peptides were not fully elucidated. In this study, peptides from tilapia skin gelatin hydrolysates with different DPP-IV inhibitory activities were analyzed by UPLC-MS/MS and multivariate analysis, and a Gly-Pro-type peptide library was constructed to elucidate their structural features. Results showed that peptide length had a dominant effect on the DPP-IV inhibition of collagen-derived peptides. Moreover, Gly-Pro-type peptides (e.g., GPA- and GPI- types) containing 4 ∼ 9 residues showed a potent DPP-IV inhibition, the IC50 values of which were 2.15 ∼ 10.43 times lower than that of Gly-Pro-Xaa tripeptides. More importantly, different proteases had discrepancy in releasing these peptides, among which papain could release them to a greater extent due to its strong preference for Arg in the S1 subsite and Pro in the S3 subsite.

Identification of repurposing therapeutics toward SARS-CoV-2 main protease by virtual screening.

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CLpro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CLpro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CLpro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CLpro inhibition. Among these compounds, lapatinib showed high efficiency of 3CLpro inhibition (IC50 value of 35 ± 1 μM and Ki of 23 ± 1 μM). The binding behavior of lapatinib against 3CLpro was elucidated by molecular dynamics simulations. This drug could well bind with 3CLpro residues in the five subsites S1’, S1, S2, S3, and S4. Moreover, lapatinib’s key chemical pharmacophore features toward SAR-CoV-2 3CLpro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CLpro.

Novel trypanocidal thiophen-chalcone cruzain inhibitors: structure- and ligand-based studies.

Background: Chagas disease is a neglected tropical disease that affects millions of people worldwide and for which no effective treatment is available. Materials & methods: 17chalcones were synthesized, for which the inhibition of cruzain and trypanocidal activity were investigated. Results: Chalcone C8 showed the highest cruzain inhibitory (IC50=0.536μm) and trypanocidal activity (IC50=0.990μm). Molecular docking studies showed interactions involving Asp161 and the thiophen group interacting with the S2 subsite. Furthermore, quantitative structure-activity relationship (q2=0.786; r2=0.953) and density functional theory studies were carried out, and a correlation between the lowest unoccupied molecular orbital surface and trypanocidal activity was observed. Conclusion: These results demonstrate that these chalcones are worthwhile hits to be further optimized in Chagas disease drug discovery programs.

Inhibitory activity of xanthoangelol isolated from Ashitaba (Angelica keiskei Koidzumi) towards α-glucosidase and dipeptidyl peptidase-IV: in silico and in vitro studies

In Indonesia, the sap of Angelica keiskei Koidzumi has been utilized traditionally as a blood-sugar reducer, nonetheless, its molecular mechanism still needs to be studied. This study aimed to isolate xanthoangelol (XA) from the yellow sap of A. keiskei planted in Mount Rinjani, Indonesia, and to investigate its mechanism by in silico and in vitro methods towards α-glucosidase and dipeptidyl peptidase-IV (DPP-IV). The dried yellow sap was macerated using ethanol, subjected to liquid-liquid extraction using a different polarity of solvents, further gradient-eluted with column chromatography. The isolated compound, formed as yellow crystals, melting point 114–114.4 °C, λmax 368 nm, m/z 393.20 [M + H]+, was confirmed as XA. Acarbose, an α-glucosidase inhibitor, and sitagliptin, a DPP-IV inhibitor, respectively, were employed as the reference drugs for both the in silico and in vitro studies. XA interacts with essential amino acid residues 232–237 in the N-terminal N-loop of α-glucosidase by forming a hydrogen bond with Ala234, a salt-bridge with Asp232, and 9 hydrophobic interactions (binding energy -7.81 kcal/mol; Ki = 1.99 μM). These binding modes resemble those of acarbose. Moreover, XA forms hydrogen bonds with Glu205 and Glu206 in the subsite S2 and π-π interaction with Phe357 in the extensive subsite S2 of DPP-IV (binding energy -8.34 kcal/mol; Ki = 0.873 μM), which are similar to those of sitagliptin. XA inhibits both α-glucosidase (IC50 XA = 14.45 μM; IC50 acarbose = 207 μM) and DPP-IV (IC50 XA = 10.49 μM; IC50 sitagliptin = 0.87 μM). Taken together, XA isolated from the yellow sap of A. keiskei Koidzumi might possess the potential to be further developed as an inhibitor of α-glucosidase and DPP-IV.

Arginine-containing dipeptides decrease affinity of gut trypsins and compromise soybean pest development

The design and production of molecules capable of mimicking the binding or/and functional sites of proteins inhibitors represent a promising strategy for the exploration and modulation of gut trypsin function in insect pests, specifically Lepidoptera. Here, for the first time, we characterized the trypsin activity present in the gut, performance and development of Anticarsia gemmatalis (Lepidoptera: Noctuidae) larvae when exposed to arginine-containing dipeptides. In silico assessment showed that arginine-containing dipeptides have a greater affinity for the active site of A. gemmatalis trypsins than lysine-containing peptides due to the presence of the double-charged guanidinium group that enhances the interaction at the S1 subsite of trypsins. Furthermore, the inhibitory and anti-insect potential of the peptides was demonstrated through kinetic and larval life cycle parameters, respectively. These dipeptides showed structural stability, binding to the active site, corroborated in vitro (competitive inhibition), and significant reduction of trypsin enzyme activity in the gut, survival, and weight of the A. gemmatalis larvae. Our findings reinforce the idea that small peptides are promising candidates for lepidopteran pest management. The optimization of DI2 and DI1 peptides, enhancing uptake and affinity to trypsins, may turn the use of these molecules feasible in agriculture.

Repurposing of existing antibiotics for the treatment of diabetes mellitus.

The online version contains supplementary material available at 10.1007/s40203-021-00118-6.

Targeting the S2 Subsite Enables the Structure-Based Discovery of Novel Highly Selective Factor XIa Inhibitors.

FXIa inhibition has been a promising strategy for treating thrombotic diseases. Up to date, many small-molecule FXIa inhibitors have been identified; however, most of them exhibit undesirable selectivity over the homologous plasma kallikrein (PKal). By employing structure-based drug design strategies, we identified many novel selective FXIa inhibitors that have extra interactions with the S2 subsite of FXIa. Among them, compound 35 displayed good inhibitory activity against FXIa and high selectivity over PKal and even several other serine proteases. Additionally, 35 showed significant anticoagulant activity toward the intrinsic pathway without affecting the extrinsic pathway. In vivo, 35 exhibited significant antithrombotic activity without increasing the bleeding risk and obvious toxicity in mice, demonstrating that it could be a promising candidate for further research. This study first demonstrates the importance of the S2 subsite of FXIa, paving the way to design highly selective FXIa inhibitors for clinical uses.

In silico screening-based discovery of novel covalent inhibitors of the SARS-CoV-2 3CL protease

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) 3CL protease (3CLpro) has been regarded as an extremely promising antiviral target for the treatment of coronavirus disease 2019 (COVID-19). Here, we carried out a virtual screening based on commercial compounds database to find novel covalent non-peptidomimetic inhibitors of this protease. It allowed us to identify 3 hit compounds with potential covalent binding modes, which were evaluated through an enzymatic activity assay of the SARS-CoV-2 3CLpro. Moreover, an X-ray crystal structure of the SARS-CoV-2 3CLpro in complex with compound 8, the most potent hit with an IC50 value of 8.50 μM, confirmed the covalent binding of the predicted warhead to the catalytic residue C145, as well as portrayed interactions of the compound with S1’ and S2 subsites at the ligand binding pocket. Overall, the present work not merely provided an experiment-validated covalent hit targeting the SARS-CoV-2 3CLpro, but also displayed a prime example to seeking new covalent small molecules by a feasible and effective computational approach.

Discovery of C-12 dithiocarbamate andrographolide analogues as inhibitors of SARS-CoV-2 main protease: In vitro and in silico studies

A global crisis of coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has impacted millions of people’s lives throughout the world. In parallel to vaccine development, identifying potential antiviral agents against SARS-CoV-2 has become an urgent need to combat COVID-19. One of the most attractive drug targets for discovering anti-SARS-CoV-2 agents is the main protease (Mpro), which plays a pivotal role in the viral life cycle. This study aimed to elucidate a series of twenty-one 12-dithiocarbamate-14-deoxyandrographolide analogues as SARS-CoV-2 Mpro inhibitors using in vitro and in silico studies. These compounds were initially screened for the inhibitory activity toward SARS-CoV-2 Mpro by in vitro enzyme-based assay. We found that compounds 3k, 3l, 3m and 3t showed promising inhibitory activity against SARS-CoV-2 Mpro with >50% inhibition at 10 μM. Afterward, the binding mode of each compound in the active site of SARS-CoV-2 Mpro was explored by molecular docking. The optimum docked complexes were then chosen and subjected to molecular dynamic (MD) simulations. The MD results suggested that all studied complexes were stable along the simulation time, and most of the compounds could fit well with the SARS-CoV-2 Mpro active site, particularly at S1, S2 and S4 subsites. The per-residue decomposition free energy calculations indicated that the hot-spot residues essential for ligand binding were T25, H41, C44, S46, M49, C145, H163, M165, E166, L167, D187, R188, Q189 and T190. Therefore, the obtained information from the combined experimental and computational techniques could lead to further optimization of more specific and potent andrographolide analogues toward SARS-CoV-2 Mpro.

An intermolecular salt bridge linking substrate binding and P1 substrate specificity switch of arterivirus 3C-like proteases

Equine arteritis virus (EAV) and porcine reproductive and respiratory syndrome virus (PRRSV) represent two members of the family Arteriviridae and pose a major threat to the equine- and swine-breeding industries throughout the world. Previously, we and others demonstrated that PRRSV 3C-like protease (3CLpro) had very high glutamic acid (Glu)-specificity at the P1 position (P1-Glu). Comparably, EAV 3CLpro exhibited recognition of both Glu and glutamine (Gln) at the P1 position. However, the underlying mechanisms of the P1 substrate specificity shift of arterivirus 3CLpro remain unclear. We systematically screened the specific amino acids in the S1 subsite of arterivirus 3CLpro using a cyclized luciferase-based biosensor and identified Gly116, His133 and Ser136 (using PRRSV 3CLpro numbering) are important for recognition of P1-Glu, whereas Ser136 is nonessential for recognition of P1-Gln. Molecular dynamics simulations and biochemical experiments highlighted that the PRRSV 3CLpro and EAV 3CLpro formed distinct S1 subsites for the P1 substrate specificity switch. Mechanistically, a specific intermolecular salt bridge between PRRSV 3CLpro and substrate P1-Glu (Lys138/P1-Glu) are invaluable for high Glu-specificity at the P1 position, and the exchange of K138T (salt bridge interruption, from PRRSV to EAV) shifted the specificity of PRRSV 3CLpro toward P1-Gln. In turn, the T139K exchange of EAV 3CLpro showed a noticeable shift in substrate specificity, such that substrates containing P1-Glu are likely to be recognized more efficiently. These findings identify an evolutionarily accessible mechanism for disrupting or reorganizing salt bridge with only a single mutation of arterivirus 3CLpro to trigger a substrate specificity switch.