1. Introduction
White spot syndrome is a kind of syndrome caused by white spot syndrome virus (WSSV). It is one of the main diseases of shrimp epidemic outbreak, with high mortality and serious harm to shrimp breeding industry. It is an important aquatic animal infectious disease to be reported by OIE [1]. Since the outbreak of the disease in 1992, it caused a large area of death of prawns, and has still brought great economic losses to China’s prawn breeding industry. The prevention and treatment of white spot syndrome is of great significance to the healthy cultivation of shrimp in China [2].
WSSV is the first marine biological DNA virus to obtain the whole genome in the world. The genome size of WSSV is about 300 KB. It is predicted that WSSV encodes about 180 proteins, among which more than 50 proteins encoded by WSSV virus have been identified as structural proteins [3]. Among these structural proteins, envelope protein is an important determinant of viral infection and pathogenicity [4], which is considered to be the first molecule in direct contact with the host [5] [6], so it plays a key role in cell targeting and host defense. Although the exact mechanism of WSSV virus entry into shrimp is still unknown, it is speculated that the virus entry into host cells is through the interaction between virus envelope protein and host cell receptor [7]. Since the virus enters the host cell through the interaction between the virus envelope protein and the host cell receptor, it is very important for the discovery of interaction blockers, especially for the development of antiviral drugs [8].
Among all of the known structural proteins of WSSV, VP28 is the main envelope [9] protein, encoded by the reading frame wsv421 [10]. The crystal structure of VP28 showed that VP28 is a trimmer with triple symmetry. Confirmed by immune-electron microscopy, VP28 was distributed on the outer surface of the virus [11]. In the process of infection, VP28, as a virus binding protein, binds to Rab7 protein in prawn cells and mediates the fusion of virus envelope and host cell membrane, thus helping the virus enter the host cell [12] [13], which plays a key role in the early stage of WSSV infection of prawn [14]. In the study of WSSV molecular mechanism, it was found that Rab7 (PmRab7) protein in Penaeusmonodon and PcRab7 in Penaeuschinensis can be specifically combined with VP28. However, the interaction mechanism of PcRab7 and VP28 remains unknown.
Although analysis of the crystal tertiary structure of protein complex is the most direct way to clarify the interaction mechanism, due to the limitation of experimental conditions the crystal tertiary structure of Rab7 and VP28 has not been obtained. The molecular mechanism of the interaction has not been fully elucidated. The development of structural modeling and molecular simulation technology and the improvement of simulation accuracy provide a new way and method for interaction research. In this research, we determined the interaction interface and interaction sites by using the methods of molecular simulation. This study will contribute to a deep understanding of the infection mechanism of WSSV on the host.
2. Methods
The target sequence of PcRab7 (accession NO. AEF33797) was downloaded from protein database NCBI. Meanwhile, a sequence alignment for PcRab7 against other sequences with available tertiary structure was performed on the online server BLAST (http://blast.ncbi.nlm.nih.gov). The automated sequence alignment and analysis of the template and target was carried out using the program of Espript2.2 [15]. The tertiary structure of PcRab7 was built by I-TASSER [16], which is a hierarchical approach to protein structure prediction. The predicted structures were saved in the PDB format. Furthermore, we select PROCHECK [17] to evaluate the quality of the generated models. Select PcRab7 and VP28 as the single structure and docked by employed the protein-protein docking program ZDOCK3.0.2 [18]. After docking, 2000 structures of PcRab7-VP28 complex were generated and sorted according to binding energy, the complex with the lowest binding free energy was as the initial for further analysis.
3. Results and Discussion
Through BLAST template search based on the protein data bank (PDB), it can be seen that the sequence consistency of PcRab7 and four Rab7 proteins with known crystal tertiary structure is between 86% - 91%, and the sequence coincidence rate is more than 78% (Figure 1). The differences between PcRab7 and the other four homologous proteins are mainly concentrated in the two ends of the protein, which indicates that the structure constructed with these four homologous structures as templates is reliable.
Using 1VG0 as the template, the three-level structure of PcRab7 was constructed (Figure 2(A)). According to the Ramachandran plot test (Figure 2(B)), 98.75% of the regions in the structure are within the acceptable range.
The constructed structure of PcRab7 was used for molecular docking with protein-protein docking software ZDOCK 3.0.2, and 2000 complex structures were obtained. The RMSD values of these 2000 complexes were calculated as 3.0 Å. All the structures in the largest cluster were sorted according to the free energy, and the conformation with the lowest binding free energy was selected as the conformation of PcRab7-VP28 complex (Figure 3). From Figure 3 we can easily found that, PcRab7 binds to VP28 through the second β-sheet, which is composed of residues L73-D86. In the previous research, Western blot and Dynabeads magnetic beads showed that the Rab7 protein of Penaeuschinensis (PcRab7) could bind to VP28 in the same way as the Rab7 protein of Penaeusmonodon (PmRab7). Furthermore, Western blot and FIA were used to further
Figure 1. Multiple sequence alignment of PcRab7 with 1VG0, 3LAW, 1T91 and 1VG1. The blue spiral and arrow at the top represent the α-helix and β-sheet. This figure was generated by Espript2.2.
Figure 2. PcRab7 modelled structure and its rationality test. (A) The cartoon image of the tertiary structure of PcRab7 is generated by the molecular graphics software PyMOL and displayed in rainbow color (N-terminal is blue, C-terminal is red); (B) According to the Ramachandran plot test of PcRab7, the image was generated by DS4.2, in which the red triangle represents the residues in the unacceptable range.
Figure 3. Complex of PcRab7 and VP28. VP28 on the left and PcRab7 on the right. Binding zone in PcRab7 was shown in yellow, α-helix, β-sheets and loops were shown in red, blue and gray respectively.
analyze the specific action fragments of PcRab7 and VP28, and the preliminary conclusion was that the L73-T111 and P157-D192 (L41-T79 and P125-D160 in PmRab7) regions of PcRab7 had binding activity with VP28. Results of region prediction analysis based on docking study are consistent with those based on experiment. Through analysis of the interaction interface of the complex structure, we can see that E81, F77 and D76 of PcRab7 bind to VP28 by Pi-Pi and Pi-Alkey interaction (Figure 4). Therefore, residues E81, F77 and D76 were speculated to be the key site affecting the binding PcRab7 and VP28.
Figure 4. Binding sites of PcRab7 and VP28. Residues in PcRab7 were shown in cartoon stick and those in VP28 were shown in stick.
4. Conclusion
It is not difficult to find that the interaction mechanism of PcRab7 and VP28 is important and urgent in both theoretical research and application needs. To solve this problem, this research first modeled the tertiary structure of PcRab7 and then docked it to VP28 by employing the method of protein-protein docking. Furthermore, we determined the interaction interface is in the zone of L73-D86, which is consistent with the previous experimental results. Meanwhile the interaction sites of complex of PcRab7 and VP28 in PcRab7 are E81, F77 and D76. These results will contribute to a deep understanding of the infection mechanism of WSSV on the host, and provide a target and theoretical reference for further drug screening and design.
Acknowledgements
This work is supported by the Central Public-interest Scientific Institution Basal Research Fund, YSFRI, CAFS (NO. 20603022019023 and 20603022017006).