Skip to main content

Structural insights into inhibition of PRRSV Nsp4 revealed by structure-based virtual screening, molecular dynamics, and MM-PBSA studies



Porcine reproductive and respiratory syndrome respiratory sickness in weaned and growing pigs, as well as sow reproductive failure, and its infection is regarded as one of the most serious swine illnesses worldwide. Given the current lack of an effective treatment, in this study, we identified natural compounds capable of inhibiting non-structural protein 4 (Nsp4) of the virus, which is involved in their replication and pathogenesis.


We screened natural compounds (n = 97,999) obtained from the ZINC database against Nsp4 and selected the top 10 compounds for analysing protein–ligand interactions and physicochemical properties. The five compounds demonstrating strong binding affinity were then subjected to molecular dynamics simulations (100 ns) and binding free energy calculations. Based on analysis, we identified four possible lead compounds that represent potentially effective drug-like inhibitors.


These methods identified that these natural compounds are capable of inhibiting Nsp4 and possibly effective as antiviral therapeutics against PRRSV.


Porcine reproductive and respiratory syndrome virus (PRRSV) infection is an economically important disease in swine and accountable for significant losses to the pork industry worldwide [1]. PRRSV is an enveloped, single-stranded RNA virus of the genus Arterivirus [1,2,3,4] that causes respiratory disease and is responsible for severe reproductive failure in pigs [4]. Generally, the disease is further complicated by secondary infection and results in a high mortality rate [1]. The available treatments, including vaccines, poorly control the disease, owing to the high genetic and antigenic heterogeneity of the virus [5]. Antivirals might be useful in controlling and managing PRRSV, and recent studies have reported the ability of herbal extracts to inhibit PRRSV infection [6,7,8]. Many compounds derived from natural sources such as plants have shown inhibitory activity against viruses and a wide variety of pathogens [9, 10]. Several natural compounds have demonstrated antiviral activity (for instance, immense potential for inhibiting viral replication) and drug-like activity [10]. The identification of potent inhibitors of PRRSV from natural sources is challenging, as it requires compound extraction and evaluation of antiviral activity, both of which require funds and specific experimental equipment. Given the difficulty in screening large numbers of natural compounds, chemo-informatics approaches are useful for identifying lead compounds from databases by targeting essential viral proteins [11, 12].

Open reading frame (ORF)1a and ORF1b comprise ~ 80% of the PRRSV genome and respectively encode pp1a and pp1b polyproteins that are cleaved by viral proteases into non-structural proteins [13]. Papain-like proteases (PL1pro and PL2pro) and a 3C-like serine protease [3CLSP; non-structural protein 4 (Nsp4)] are the viral proteases involved in polyprotein cleavage and required for Arterivirus replication. Additionally, the reported involvement of Nsp4 in interferon inhibition is linked with PRRSV pathogenesis, suggesting it as a promising molecular target for novel therapeutics [14,15,16].

The goal of this study was to identify natural compounds capable of serving as novel inhibitors of PRRSV replication via their targeting of Nsp4. Specifically, we aimed to identify natural compounds through structure-based virtual screening, analyse their physicochemical properties, perform molecular dynamics (MD) simulations, and determine the binding affinities of the potential inhibitors using molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) methods.

Materials and methods

Retrieval and preparation of ligand structures

We retrieved 97,999 compounds from a subset of the ZINC database housing natural compounds [17]. All natural compounds were downloaded in the structure-data file format, and these files were subsequently converted to AutoDock PDBQT [Protein Data Bank (PDB), partial charge (Q), and atom type (T)] files using OpenBabel software ( These files were then used for interaction studies with target protein through structure-based virtual screening and molecular docking using AutoDock vina [18,19,20].

Retrieval and preparation of the target protein structure

The crystal structure of PRRSV Nsp4 (PDB ID: 5Y4L) was retrieved from RCSB-PDB and visualised and analysed using UCSF Chimera-1.15 software [21, 22]. An AutoDock tool was used for the addition of partial atomic charges (Kollman charge), hydrogen atoms, generation of gridbox, and preparation of the Nsp4 structure. The grid box was generated with centre (X = − 4.034, Y = 6.285, Z = 17.57) and size (X = 46, Y = 44, Z = 42) coordinates that were defined in a configuration file (exhaustiveness and energy ranges: 8 and 4, respectively). The prepared structure was saved in the PDBQT file format for molecular docking [23].

Structure-based virtual screening

Virtual screening is a computational method widely used for the identification of lead molecules by docking large numbers of compounds with a molecular target of interest to allow evaluation of the binding free energy of the docked/screened compounds using drug discovery programs [24]. AutoDock vina is a molecular docking and virtual screening program that determines the preferred relative orientation of a ligand during docking or interaction with a molecular target and provides a stable protein–ligand complex structure that exhibits a minimum binding energy [20]. Here, we used AutoDock vina to screen the retrieved natural compounds against PRRSV Nsp4 and generated protein–ligand complexes of the top 10 screened compounds using PyMol ( Two-dimensional (2D) models of the complexes were visualised using Discovery Studio Visualizer ( to determine the amino acid residues involved in the interactions [25].

Drug-likeness analysis

A total of seven principal descriptors were included to evaluate the drug-likeness of the top 10 screened natural compounds. These included molecular weight (MW), logP value, status as a hydrogen-bond donor (HBD), and acceptor (HBA), polar surface area (2D; PSA), polarizability (P), and van der Waals surface area (VWSA). MW, logP, HBD, and HBA were obtained from the ZINC database [17], whereas PSA, P, and VWSA were calculated using the MarvinSketch software (

MD simulation

The MD simulations were performed using GROMACS (v.2018.1; for stability predictions of the Nsp4–ligand complexes [26, 27]. Six systems were generated and subjected to 100 ns MD simulations—one system estimated the stability of Nsp4 and the other five estimated the stability of the Nsp4–ligand complexes. The dynamic nature of the target protein and the docked-ligand complex was predicted in the presence of a solvent. All six systems were solvated in a box using a simple point-charge model. The topology of the ligands was created using ProDRG [28] and that for Nsp4 was generated using the GROMOS9653a6 force field [29]. The systems were neutralised by adding 1 Na+ ion. To eliminate steric hindrance, the steepest energy minimization was used for all systems in order to obtain the maximal force below 1000 kJ/mol/nm. Long-range electrostatic interactions were determined using the particle mesh Ewald (PME) method [30]. For computation of Lennard–Jones and Coulomb interactions, we used a radius cut-off of 1.0 nm; the LINCS algorithm was used to constrain H-bond lengths [31]. All simulations applied a consistent time step of 2 fs. Short-range non-bonded interactions were predicted using a 10-Å cut-off distance, whereas long-range electrostatics were predicted using the PME method with 1.6-Å Fourier grid spacing. Shake algorithms were used to fix all bonds, including H-bonds [32]. After energy minimisation, the systems were equilibrated, followed by position-restraint simulations under NVT and NPT conditions to maintain the volume, temperature, and pressure. Finally, all systems were subjected to a 100 ns MD simulation; coordinates were saved at 2 fs intervals. Root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), H-bonds, and Gibbs free energy landscapes were calculated, and principal component analysis (PCA) was performed to predict correlated motions generated during protein–ligand interactions;‘gmx rms’, ‘gmx rmsf’, ‘gmx gyrate’, ‘gmx sasa’, ‘gmx hbond’, ‘gmx sham’, and ‘gmx covar’ (GROMACSv.2018.1;, respectively, were used for these purposes. The resulting files were analysed and visualised using xmgrace (

Binding-energy calculation

MM-PBSA is a widely used and well-accepted method for calculating the binding free energy of protein–ligand complexes [33]. Here, we used the g_mmpbsa tool ( to calculate the binding energy by integrating high-throughput MD simulation data [34]. The binding energy calculations can be described in the following equation:

$$ \Delta {\mathrm{G}}_{\mathrm{bind}}=\Delta {\mathrm{G}}_{\mathrm{mm}}+\Delta {\mathrm{G}}_{\mathrm{ps}}+\Delta {\mathrm{G}}_{\mathrm{nps}}-\mathrm{T}\Delta \mathrm{S} $$

Here, van der Waals and electrostatic interaction were calculated in molecular mechanics energy (ΔGmm). ΔGps and ΔGnps are the polar and non-polar solvation energies, and TΔS is refer to the entropic contribution where temperature and entropy are denoted by T and S, respectively. The average binding energies of the top five protein–ligand complexes and amino acid residues contributing to the binding activity were calculated by using ‘python’ scripts included in the g_mmpbsa tool.


Identification of lead compounds through structure-based virtual screening

Structure-based virtual screening enables the prediction of optimal interactions between ligands and a macromolecular target for complex formation. The ligands are subsequently sorted according to their binding free energy for the target. This requires the three-dimensional structure of the target, with the compounds obtained from a database and categorised according to their affinity. In the present study, we downloaded a subset of natural compounds (n = 97,999) from the ZINC database for virtual screening against PRRSV Nsp4. We subsequently identified the top 10 compounds sorted according to their minimum binding free energy (range:− 10.0 to − 9.2 kcal/mol) for further analyses (Table 1).

Table 1 Binding free energies of the top 10 screened compounds along with the amino acid residues involved in interactions. The amino acid residues shown in bold are involved in hydrogen-bonding interactions

Analysis and visualisation of the screened compounds

The active site residues of Nsp4 include His39, Asp64, and Ser118, as well as His133 and Ser136 that are reportedly essential for protein activity. The compound showing the optimal binding free energy (− 10.0 kcal/mol; ZINC38167083) demonstrated ligand interactions such as His39-mediated van der Waals interactions and Asp64-mediated Pi-anion interaction. Fig. 1 shows protein–ligand H-bond interactions involving the active site residues His39, Ser118, and Ser136, besides, Asp64 and His133 is interacted with van der Waals and pi-anion interaction with ZINC08877407. The interacting amino acid residues of top 10 compounds are shown in Table 1.

Fig. 1
figure 1

2D representation of the binding interactions of top five screened natural compounds with Nsp4 depicted key amino acid residues contributed in protein-ligand interactions. A ZINC38167083, (B) ZINC16919178, (C) ZINC08792350, (D) ZINC01510656, and (E) ZINC08877407

Assessment of drug likeness through physicochemical property analysis

Physicochemical property analysis is one of the fundamental tasks in any drug discovery program. The top 10 screened compounds were then subjected to analysis of their physicochemical properties according to 7 principal descriptors (MW, logP value, status of HBDs and HBAs, 2D PSA, P, and VWSA). According to a previous study, a good drug should have an MW < 500 Da, an HBD < 5, and an HBA < 10. The MW, logP, HBD, and HBA of the selected compounds met the Lipinski’s rule. Additionally, PSA, P, and VWSA results displayed drug-like behaviour. (Table 2).

Table 2 Physicochemical properties of the top 10 screened compounds

MD simulation analysis

The structure of Nsp4 and top five screened compound-complex with Nsp4 was employed for 100 ns MD simulation study for predicting the dynamic changes during protein-ligand interaction and their nature of stability. The present study included various parameters i.e. RMSD, RMSF, Rg, SASA, H-bond, PCA, Gibbs free energy landscape, and binding free energy calculation.

Structural deviation analysis through RMSD

The RMSD value describes the dynamic behaviour among native structures to a new pose. After a 70 ns of simulation to obtain a stable trajectory, the RMSD values were 0.35, 0.25, 0.29, 0.23, 0.38, and 0.39 nm for Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407, respectively. These data suggest that Nsp4-ZINC08792350 and Nsp4-ZINC38167083 are highly stable complexes relative to the others. Because each Nsp4–compound complex demonstrated stability after the 70 ns simulation, we performed further evaluations on each for last 30 ns trajectory (Fig. 2A).

Fig. 2
figure 2

Stability analysis (A) RMSD values for the Nsp4–compound complexes. Flexibility analysis (B) RMSF values for the Nsp4–compound complexes over the final 30 ns of the simulations. Compactness (C) Rg, and Solvent accessible surface area analysis (D) SASA values for the final 30 ns of the simulations. Black, red, green, blue, orange, and violet colours represent Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407, respectively. E Changes in the number of hydrogen bonds in each respective complex according to data from the final 30 ns of the simulations. Red, green, blue, orange, and violet colour represent Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 respectively

Flexibility analysis through RMSF

Evaluation of the RMSF values used to assess structural rigidity revealed values of 0.08, 0.11, 0.12, 0.10, 0.11, and 0.11 nm for Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407, respectively (Fig. 2B). Higher RMSF values were due to ligand binding, causing alterations in protein geometry. Minimal fluctuations were observed in Nsp4-ZINC08792350 and Nsp4-ZINC38167083 complex compared with that in other complexes.

Radius of gyration (Rg) analysis

Assessment of complex compactness according to Rg calculation revealed values of 1.50, 1.27, 1.43, 1.46, 1.39, and 1.44 nm for Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407, respectively (Fig. 2C). The results indicate that the Nsp4-ZINC38167083 complex showed a more compact structure than the other complexes.

Solvent accessible surface area (SASA) analysis

To identify changes in the solvent-accessible regions of the complexes, we determined SASA values over the course of the final 30 ns of the simulation. Our study revealed values of 95.88, 98.33, 98.98, 97.13, 96.97, and 100.92 nm2 for Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 (Fig. 2D), revealing relatively minimal changes after binding by each of the compounds.

Interaction analysis through hydrogen bonding

Hydrogen bonding is the most important bond for stabilizing protein–ligand interactions. The average number of hydrogen bonds for the complexes Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 over the final 30 ns of the simulations was 0–1 and that for Nsp4-ZINC38167083 and Nsp4-ZINC16919178 was 0–2 and 0–3, respectively (Fig. 2E). Hence, these compounds interacted with Nsp4 and provided a stable complex during protein–ligand interactions.

Principal component analysis (PCA)

In PCA, the sum of the eigenvalues suggests the overall flexibility of a structure under different conditions. Therefore, the first 5 of 50 eigenvectors used to calculate eigenvalues from the final 30 ns of the simulation were used to determine the percentage change in structural movement. The results revealed that these five eigenvectors accounted for 42.85, 63.97, 63.27, 59.14, 64.83, and 71.05% of the motions for Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 respectively (Fig. 3A), suggesting increased movement after the binding of each ligand. Moreover, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, and Nsp4-ZINC01510656 showed less overall motion relative to Nsp4-ZINC08877407. Additionally, generation of a 2D plot for assessing protein dynamics after ligand binding suggested the overall stability (lowcorrelated motions) of Nsp4, Nsp4-ZINC38167083, and Nsp4-ZINC08792350(Fig. 3B), indicating these compounds as possible leads for further evaluation as inhibitors.

Fig. 3
figure 3

Principal component analysis (A) Eigenvalues derived from the final 30 ns of each simulation and used for PCA depicted Eigenvalues vs. first fifty eigenvector. B First two eigenvectors depicted Nsp4 motion in space for all the systems. Black, red, green, blue, orange, and violet colours represent Nsp4, Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 respectively

Gibbs free energy landscape

We then calculated the Gibbs free energy landscape using the first two principal components (PC1 and PC2) in order to visualize the results. Fig. 4 shows the colour-coded plots generated for Nsp4 along with each complex. The lowest free energy values (≤9.08 kJ/mol) were observed for Nsp4-ZINC38167083, suggesting that this complex demonstrated overall thermodynamic stability. The other complexes (Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407) had values of to 11.4 kJ/mol, implying that these complexes have numerous high-energy minima.

Fig. 4
figure 4

The color-coded illustration of the Gibbs free energy landscape plotted using PC1 and PC2. The lower energy systems are represented by the deeper blue color on the contour map. A Nsp4, (B) Nsp4-ZINC38167083, (C) Nsp4-ZINC16919178, (D) Nsp4-ZINC08792350, (E) Nsp4-ZINC01510656, and (F) Nsp4-ZINC08877407

Binding free energy

We then evaluate the binding free energy associated with each ligand through MM-PBSA using the final 10 ns of the simulation, for calculation of van der Waals and electrostatic interactions, Polar solvation, and SASA. The calculated binding free energy for Nsp4-ZINC38167083, Nsp4-ZINC16919178, Nsp4-ZINC08792350, Nsp4-ZINC01510656, and Nsp4-ZINC08877407 was − 124.54, − 128.44, − 159.33, − 122.50, and − 78.19 kJ mol− 1 respectively (Table 3).

Table 3 Average binding free energies of Nsp4 complexes in kJ mol − 1

The investigation of residual binding energy is a key method for identifying residues important to ligand binding. Fig. 5 shows that amino acid residues at positions 5 to 142 contributed significantly to binding of ZINC38167083, ZINC16919178, ZINC08792350, and ZINC01510656, which are the catalytic residues in the active site. Fewer contacts were observed in relation to ZINC08877407 binding, suggesting that ZINC38167083, ZINC16919178, ZINC08792350, and ZINC01510656 represent potential Nsp4 inhibitors.

Fig. 5
figure 5

Plot depicting the amino acid residues of Nsp4 contributing to the binding with natural compounds. Red, green, blue, orange, and violet colours represent ZINC38167083, ZINC16919178, ZINC08792350, ZINC01510656, and ZINC08877407, respectively


PRRSV is a recalcitrant and intricate disease in a pig when working as a cofactor in a porcine respiratory disease complex (PRDC) or primary infectious agent. It was identified as the most frequent virus linked to PRDC [35,36,37,38,39]. Furthermore, PRRSV has been shown to impair the host immune system, which can lead to more serious secondary infections, and chronic disorders [35]. The involvement of Nsp4 in PRRSV replication and pathogenesis is decoded and recommended as one of the key molecular targets for drug development [14]. Therefore, identification of Nsp4 inhibitors is needed to prevent and manage the disease. Natural compounds have made immense contributions in the identification of lead molecule(s) with antiviral potential. It is believed that the disease can be controlled successfully by developing small molecules that can inhibit Nsp4 activity linked with pathogenesis [14]. In the present study, computational approaches are utilized for the identification of possible lead compounds via molecular docking of natural compounds database through structure-based virtual screening followed by downstream analysis. Structure based-virtual screening is a powerful computational approach that is used to investigate important lead molecule(s) from a big set of a compound database based on the lowest binding energy required for stabilizing the protein-ligand complex [40].

From the structure-based virtual screening, we have selected the top ten natural compounds that show interaction with key residues. Further, protein-ligand analysis of the top 5 compounds demonstrated that the ZINC38167083 interacted with Nsp4 and formed one conventional hydrogen bond at position Gly63. Besides, amino acid residues Ser18, His39, Leu41, Thr42, Gly43, Asn44, Thr134, and Thr145 were involved in van der Waals interactions, and Val61 and Ile143 formed alkyl bonds. Ala38 and Phe151 participated in interaction through amide-pi and pi-pi t-shaped bonding. Additionally, Asp64 contributed to interaction through the pi-anion bond. ZINC16919178 bonded with Nsp4 at position Asn11 by one conventional hydrogen bond. In addition, amino acid residues Phe3, Thr5, Ser9, Leu10, Phe26, Ile123, and Gly127 formed van der Waals interactions; Pro78, Tyr92, Leu94, Val99, and Pro101 formed alkyl and pi-alkyl bonds. ZINC08792350 interacted with Nsp4 Thr5, Ser9, Asn11, Val76, Pro78, Ile123, Phe166, and Asp192 through van der Waals interactions; Phe3 and Tyr92 formed pi-pi t-shaped bonding. In addition, the amino acid residues Leu10, Pro101, Leu196 contributed to interaction through pi-alkyl bonding; Arg90 and Val99 formed pi-anion and pi-sigma bonds with ZINC08792350, respectively. ZINC01510656 bonded with Nsp4 at Thr5, Asn11, Val76, Phe26, Leu94, and Gly127 through van der Waals interactions; amino acid residues Leu10, Pro78, and IIe123 formed pi-alkyl bonds, and Ser9, Tyr92, Val99 contributed in interaction by carbon-hydrogen bonding, pi-pi t-shaped, and pi-sigma interactions, respectively. ZINC08877407 formed conventional hydrogen bonds with Nsp4 at position His39, Ser118, ASP117, and Ser136; amino acid residues Gly63, Asp64, Ala114, Cys115, Gly116, Thr134, Lys138, and Thr145 contributed to protein-ligand interaction through van der Waals. Additionally, Gly135 formed a carbon-hydrogen bond, and His133, Ile143, and Phe151 interacted through pi-anion, pi-alkyl, and pi-sigma bonding, respectively. Medicinal chemists have traditionally been interested in noncovalent interactions that are indicative of attraction, directed intermolecular forces in their search for the “glue” that keeps ligand and their molecular target together. In recent years, with the rapid increase in the number of solved biomolecular structures and the performance enhancement of computational methods, it is now possible to provide a more thorough understanding of protein-ligand interaction [41]. Therefore, based on the results, it was concluded that the screened compounds can inhibit the virulence activity of Nsp4 [14]. Besides, the results of physicochemical properties prediction suggest that the screened compounds demonstrated good drug-like behavior and could be considered for further analysis [42, 43]. Therefore, 100 ns MD simulation analysis was conducted for Nsp4 and top 5 natural compounds i.e. ZINC38167083, ZINC16919178, ZINC08792350, ZINC01510656, and ZINC08877407, respectively with Nsp4 to evaluate the dynamic behavior of protein and protein-ligand complexes. It is recognized as a powerful approach for predicting the conformational stability of macromolecules before and after ligand binding, besides the simulated data can be utilized for calculation of real binding energy of small molecules concerning time along with a contribution of binding amino acid residues present in the macromolecular target [44]. Several structural parameters were calculated, including RMSD, RMSF, Rg, SASA, H-bonding, PCA, and gibbs free energy [45,46,47,48]. The RMSD value indicated that all of the complexes were stable and creating an equilibrated trajectory for further investigation. As a result, we determined RMSF, Rg, SASA, PCA, and Gibbs free energy to determine the nature of each system subjected for MD simulation. Drug selectivity, metabolization, and stability all require H-bonds. To better understand the H-bond and its contributions to the overall stability of each system, an H-bond analysis of natural compounds-Nsp4 complexes were calculated. The hydrogen bonding study indicates that all of the Nsp4-complexes are stable and made bonding with essential catalytic residues [49]. The overall analysis revealed that each complex was stabilizing after 70 ns indicating better interaction with Nsp4 in terms of stability that is required for its inhibition. Further, MM-PBSA binding free energy and residual binding energy were calculated to assess the binding affinities of natural compounds with Nsp4. For determining the binding free energy of protein–ligand complexes by using MD simulation trajectory, it is a frequently used and well-accepted method [33, 50]. The strength of the binding contacts between the ligand and the target protein is measured by ligand binding affinity, which is directly linked to ligand potency. In the field of drug discovery, its evaluation is crucial. Furthermore, in favorable reactions, the free energy is negative. So, lowering the binding energy improves interactions, and low binding energy corresponds to the high binding affinity of protein-ligand complexes [51]. MM-PBSA analysis demonstrated that the compounds ZINC38167083, ZINC16919178, ZINC08792350, and ZINC01510656 can act as a potential lead for inhibition of Nsp4 [52, 53]. Whereas, ZINC08877407 was not recommended as a lead because their binding energy was found to be higher as compared to other compounds.

In past years, the identification of lead compounds for drug development take much time and cost as well as required good infrastructure experimental facilities [11, 54]. Due to advances in structural biology, computer science, and bioinformatics, it becomes easy to find out putative molecule(s) by a screening of a big database that has a strong affinity with the target for experimental evaluation [24, 55]. It saves the cost and time of the scientific community. Most of the medicines available in the market are from a natural source or it is a derivative of naturally occurring molecules [11, 24]. Natural compounds have immense potential to inhibit virus and pathogenic proteins and act as antiviral drugs [56,57,58]. The results presented in this work are, therefore, informative for understanding the antiviral potential of suggested compounds as therapeutics for PRRSV. It might be also useful for the prevention of pigs and other animals from different viral diseases [59, 60].


PRRSV infection is a main concern for the global swine industry, and there is a need to identify novel and effective therapeutic agents. Given the importance of Nsp4 in PRRSV replication and pathogenesis, we employed computational and MD approaches to screen and identify natural compounds as novel inhibitors of Nsp4 activity. The results identified four possible lead compounds that represent potentially effective drug-like inhibitors for application as antiviral therapeutics. Further studies are warranted to confirm these findings through experimental and clinical evaluations in order to promote future management of PRRSV infection.

Availability of data and materials

All data generated or analysed during this study are included in the manuscript.





3C-like serine protease


Hydrogen-bond acceptor


Hydrogen-bond donor


Molecular dynamics


Molecular mechanics Poisson–Boltzmann surface area


Molecular weight


Non-structural protein 4


Open reading frame


Papain-like protease


Particle mesh Ewald


Porcine reproductive and respiratory syndrome virus


Polar surface area


van der Waals surface area


  1. Huang C, Bernard D, Zhu J, Dash RC, Chu A, Knupp A, et al. Small molecules block the interaction between porcine reproductive and respiratory syndrome virus and CD163 receptor and the infection of pig cells. Virol J. 2020;17(1):1–11.

    Article  Google Scholar 

  2. Cavanagh D. Nidovirales: a new order comprising Coronaviridae and Arteriviridae. Arch Virol. 1997;142:629–33.

    Google Scholar 

  3. Benfield DA, Nelson E, Collins JE, Harris L, Goyal SM, Robison D, et al. Characterization of swine infertility and respiratory syndrome (SIRS) virus (isolate ATCC VR-2332). J Vet Diagn Investig. 1992;4(2):127–33.

    Article  Google Scholar 

  4. Music N, Gagnon CA. The role of porcine reproductive and respiratory syndrome (PRRS) virus structural and non-structural proteins in virus pathogenesis. Anim Health Res Rev. 2010;11(2):135–63.

    Article  Google Scholar 

  5. Meng XJ. Heterogeneity of porcine reproductive and respiratory syndrome virus: implications for current vaccine efficacy and future vaccine development. Vet Microbiol. 2000;74(4):309–29.

    Article  Google Scholar 

  6. Sun N, Li E, Wang Z, Zhao J, Wang S, He J, et al. Sodium tanshinoneIIA sulfonate inhibits porcine reproductive and respiratory syndrome virus via suppressing N gene expression and blocking virus-induced apoptosis. Antivir Ther. 2014;19:89–95.

    Article  Google Scholar 

  7. Zhou DH, Wang X, Yang M, Shi X, Huang W, Feng Q. Combination of low concentration of (−)-epigallocatechin gallate (EGCG) and curcumin strongly suppresses the growth of non-small cell lung cancer in vitro and in vivo through causing cell cycle arrest. Int J Mol Sci. 2013;14(6):12023–36.

    Article  Google Scholar 

  8. Du T, Nan Y, Xiao S, Zhao Q, Zhou EM. Antiviral strategies against PRRSV infection. Trends Microbiol. 2017;25(12):968–79.

    Article  Google Scholar 

  9. Evans AB, Dong P, Loyd H, Zhang J, Kraus GA, Carpenter S. Identification and characterization of small molecule inhibitors of porcine reproductive and respiratory syndrome virus. Antivir Res. 2017;146:28–35.

    Article  Google Scholar 

  10. Zitterl-Eglseer K, Marschik T. Antiviral medicinal plants of veterinary importance: a literature review. Planta Med. 2020;86(15):1058–72.

    Article  Google Scholar 

  11. Agnihotry S, Pathak RK, Srivastav A, Shukla PK, Gautam B. Molecular docking and structure-based drug design. In: Computer-aided drug design. Singapore: Springer; 2020. p. 115–31.

    Chapter  Google Scholar 

  12. Singh DB, Pathak RK. Computational approaches in drug designing and their applications. In: Experimental protocols in biotechnology. New York, NY: Humana; 2020. p. 95–117.

    Chapter  Google Scholar 

  13. Ziebuhr J, Snijder EJ, Gorbalenya AE. Virus-encoded proteinases and proteolytic processing in the Nidovirales. Microbiol. 2000;81(4):853–79.

    Article  Google Scholar 

  14. Tian X, Lu G, Gao F, Peng H, Feng Y, Ma G, et al. Structure and cleavage specificity of the chymotrypsin-like serine protease (3CLSP/nsp4) of porcine reproductive and respiratory syndrome virus (PRRSV). J Mol Biol. 2009;392(4):977–93.

    Article  Google Scholar 

  15. Shi Y, Lei Y, Ye G, Sun L, Fang L, Xiao S, et al. Identification of two antiviral inhibitors targeting 3C-like serine/3C-like protease of porcine reproductive and respiratory syndrome virus and porcine epidemic diarrhea virus. Vet Microbiol. 2018;213:114–22.

    Article  Google Scholar 

  16. An TQ, Li JN, Su CM, Yoo D. Molecular and cellular mechanisms for PRRSV pathogenesis and host response to infection. Virus Res. 2020;286:197980.

    Article  Google Scholar 

  17. Irwin JJ, Sterling T, Mysinger MM, Bolstad ES, Coleman RG. ZINC: a free tool to discover chemistry for biology. J Chem Inf Model. 2012;52(7):1757–68.

    Article  Google Scholar 

  18. Dalby A, Nourse JG, Hounshell WD, Gushurst AK, Grier DL, Leland BA, et al. Description of several chemical structure file formats used by computer programs developed at molecular design limited. J Chem Inf Comput Sci. 1992;32(3):244–55.

    Article  Google Scholar 

  19. O'Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open babel: an open chemical toolbox. J Chem. 2011;3(1):1–14.

    Article  Google Scholar 

  20. Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.

    Article  Google Scholar 

  21. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The protein data bank. Nucleic Acids Res. 2000;28(1):235–42.

    Article  Google Scholar 

  22. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12.

    Article  Google Scholar 

  23. Goodsell DS, Morris GM, Olson AJ. Automated docking of flexible ligands: applications of AutoDock. J Mol Recognit. 1996;9(1):1–5.<1::AID-JMR241>3.0.CO;2-6.

    Article  Google Scholar 

  24. Pathak RK, Singh DB, Sagar M, Baunthiyal M, Kumar A. Computational approaches in drug discovery and design. In: Computer-aided drug design. Singapore: Springer; 2020. p. 1–21.

    Chapter  Google Scholar 

  25. Jejurikar BL, Rohane SH. Drug designing in discovery studio. Asian J Res Chem. 2021;14:135–8.

    Google Scholar 

  26. Pronk S, Páll S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013;29(7):845–54.

    Article  Google Scholar 

  27. Abraham MJ, Murtola T, Schulz R, Páll S, Smith JC, Hess B, et al. GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 2015;1:19–25.

    Article  Google Scholar 

  28. Schüttelkopf AW, Van Aalten DM. PRODRG: a tool for high-throughput crystallography of protein–ligand complexes. Acta Crystallogr D Biol Crystallogr. 2004;60(8):1355–63.

    Article  Google Scholar 

  29. Oostenbrink C, Villa A, Mark AE, Van Gunsteren WF. A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6. J Comput Chem. 2004;25(13):1656–76.

    Article  Google Scholar 

  30. Darden T, York D, Pedersen L. Particle mesh Ewald: an N· log (N) method for Ewald sums in large systems. J Chem Phys. 1993;98(12):10089–92.

    Article  Google Scholar 

  31. Hess B, Bekker H, Berendsen HJ, Fraaije JG. LINCS: a linear constraint solver for molecular simulations. J Comput Chem. 1997;18(12):1463–72.<1463::AID-JCC4>3.0.CO;2-H.

    Article  Google Scholar 

  32. Ryckaert JP, Ciccotti G, Berendsen HJ. Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977;23(3):327–41.

    Article  Google Scholar 

  33. Genheden S, Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discovery. 2015;10(5):449–61.

    Article  Google Scholar 

  34. Kumari R, Kumar R. Open source drug discovery consortium and Lynn, A. g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations. J Chem Inf Model. 2014;54(7):1951–62.

    Article  Google Scholar 

  35. Lim B, Kim S, Lim KS, Jeong CG, Kim SC, Lee SM, et al. Integrated time-serial transcriptome networks reveal common innate and tissue-specific adaptive immune responses to PRRSV infection. Vet Res. 2020;51(1):1–8.

    Article  Google Scholar 

  36. Brockmeier SL, Palmer MV, Bolin SR. Effects of intranasal inoculation of porcine reproductive and respiratory syndrome virus Bordetella bronchiseptica, or a combination of both organisms in pigs. Am J Vet Res. 2000;61(8):892–9.

    Article  Google Scholar 

  37. Allan G, McNeilly F, Kennedy S, Meehan B, Ellis J, Krakowka S. Immunostimulation, PCV-2 [porcine circovirus] and PMWS [porcine wasting syndrome]. Vet Rec. 2000;147(6):170–1.

    Google Scholar 

  38. Thacker EL, Halbur PG, Ross RF, Thanawongnuwech R, Thacker BJ. Mycoplasma hyopneumoniae potentiation of porcine reproductive and respiratory syndrome virus-induced pneumonia. J Clin Microbiol. 1999;37(3):620–7.

    Article  Google Scholar 

  39. Brockmeier SL, Halbur PG, Thacker EL. Porcine respiratory disease complex. In: Brogden KA, Guthmiller JM, editors. Polymicrobial diseases. Chapter 13. Washington: ASM press; 2002.

    Google Scholar 

  40. Maia EHB, Assis LC, de Oliveira TA, da Silva AM, Taranto AG. Structure-based virtual screening: from classical to artificial intelligence. Front Chem. 2020;8:343.

    Article  Google Scholar 

  41. Zhou P, Huang J, Tian F. Specific noncovalent interactions at protein-ligand interface: implications for rational drug design. Curr Med Chem. 2012;19(2):226–38.

    Article  Google Scholar 

  42. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 1997;23(1-3):3–25.

    Article  Google Scholar 

  43. Jiao Y, Peng J, Ye X, Hu H, Gan L, Yang J, et al. Study on pharmacological properties and cell absorption metabolism of novel daidzein napsylates. Royal Society open. 2021;science 8(1):201475.

  44. Shukla R, Singh TR. High-throughput screening of natural compounds and inhibition of a major therapeutic target HsGSK-3 β for Alzheimer’s disease using computational approaches. J Genet Eng Biotechnol. 2021;19(1):1–17.

    Article  Google Scholar 

  45. Pathak RK, Gupta A, Shukla R, Baunthiyal M. Identification of new drug-like compounds from millets as xanthine oxidoreductase inhibitors for treatment of hyperuricemia: a molecular docking and simulation study. Comput Biol Chem. 2018;76:32–41.

    Article  Google Scholar 

  46. Rai SK, Pathak RK, Singh DB, Bhatt A, Baunthiyal M. Chemo-informatics guided study of natural inhibitors targeting rho GTPase: a lead for treatment of glaucoma. In Silico Pharm. 2021;9(1):1–11.

    Article  Google Scholar 

  47. Shukla R, Singh TR. Identification of small molecules against cyclin dependent kinase-5 using chemoinformatics approach for Alzheimer’s disease and other tauopathies. J Biomol Struct Dyn. 2020:1–13.

  48. Karthick V, Nagasundaram N, Doss CGP, Chakraborty C, Siva R, Lu A, et al. Virtual screening of the inhibitors targeting at the viral protein 40 of Ebola virus. Infect Dis Pov. 2016;5(1):1–10.

    Article  Google Scholar 

  49. Agrawal S, Govind Kumar V, Gundampati RK, Moradi M, Kumar TK. Characterization of the structural forces governing the reversibility of the thermal unfolding of the human acidic fibroblast growth factor. Sci Rep. 2021;11(1):1–3. 15579.

    Article  Google Scholar 

  50. Kumar VG, Agrawal S, Kumar TK, Moradi M. Binding affinity estimation from restrained umbrella sampling simulations. bioRxiv. 2021.

  51. Wan S, Bhati AP, Zasada SJ, Coveney PV. Rapid, accurate, precise and reproducible ligand–protein binding free energy prediction. Interface Focus. 2020;10(6):20200007.

    Article  Google Scholar 

  52. Poli G, Granchi C, Rizzolio F, Tuccinardi T. Application of MM-PBSA methods in virtual screening. Molecules. 2020;25(8):1971.

    Article  Google Scholar 

  53. Kuhn B, Gerber P, Schulz-Gasch T, Stahl M. Validation and use of the MM-PBSA approach for drug discovery. J Med Chem. 2005;48(12):4040–8.

    Article  Google Scholar 

  54. Eisenstein M. Active machine learning helps drug hunters tackle biology. Nat Biotechnol. 2020;38(5):512–5.

    Article  Google Scholar 

  55. Schenone M, Dančík V, Wagner BK, Clemons PA. Target identification and mechanism of action in chemical biology and drug discovery. Nat Chem Biol. 2013;9(4):232–53 53.

    Article  Google Scholar 

  56. Omrani M, Keshavarz M, Nejad Ebrahimi S, Mehrabi M, McGaw LJ, Ali Abdalla M, et al. Potential natural products against respiratory viruses: a perspective to develop anti-COVID-19 medicines. Front Pharmacol. 2020;11:2115.

    Google Scholar 

  57. Martinez JP, Sasse F, Brönstrup M, Diez J, Meyerhans A. Antiviral drug discovery: broad-spectrum drugs from nature. Nat Prod Rep. 2015;32(1):29–48.

    Article  Google Scholar 

  58. Bhattacharya R, Dev K, Sourirajan A. Antiviral activity of bioactive phytocompounds against coronavirus: an update. J Virol Methods. 2021;23:114070.

    Article  Google Scholar 

  59. Li WG, Dai FY, Cheng YX, Yin GF, Bi JL, Li DP. Identification of porcine reproductive and respiratory syndrome virus inhibitors through an oriented screening on natural products. Chem Res Chin Univ. 2013;29(2):290–3.

    Article  Google Scholar 

  60. Gómez-García M, Puente H, Argüello H, Mencía-Ares Ó, Rubio P, Carvajal A. In vitro assessment of antiviral effect of natural compounds on porcine epidemic diarrhea coronavirus. Front Vet Sci. 2021;8.

Download references


The authors thank Chung-Ang University, Anseong-si, the Republic of Korea for providing High-Performance Computing (HPC) and other necessary facilities.


This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2018R1A6A1A03025159).

Author information

Authors and Affiliations



JMK designed the experiments and supervised the research; RKP performed experiments, analyse results and wrote the manuscript. YJS helped in analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jun-Mo Kim.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pathak, R.K., Seo, YJ. & Kim, JM. Structural insights into inhibition of PRRSV Nsp4 revealed by structure-based virtual screening, molecular dynamics, and MM-PBSA studies. J Biol Eng 16, 4 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Swine
  • Nsp4
  • Molecular dynamics
  • Protein–ligand interaction