The stability of binding ligands on free protein and the protein-ligand complex was analyzed using molecular dynamics by comparing the 50 ns MD trajectories.
3.4.1 Root Mean Square Deviation (RMSD)
The root-mean-square-deviation (RMSD) values provide deviation at Å level with time duration nanosecond (ns) level for protein and ligand, indicating the complex (compound) stability as a whole. The protein RMSD was relatively stable throughout the simulation, with insubstantial fluctuations ranging from 0.5 Å- 1 Å, which indicates that VP28 has not undergone significant conformational changes. The RMSD plot (Fig. 6a) indicates that the stigmasta-4,22-dien-3-one -VP28 (2ED6) complex stabilized at 10 ns with respect to the reference time as 0 ns after commencing simulation. The fluctuations were observed at 45 ns and was stabilized. The variations were in the order of 1 Å indicating the complex has not undergone conformational changes. The overall RMSD plot of stigmasta-4,22-dien-3-one -VP28 complex suggest that ligand is bound stably to VP28 and has not diffused away from the protein. The RMSD plot (Fig. 6b) indicates that stigmasta-4,22-dien-3-one -VP26 (2EDM) complex stabilized at 20 ns with respect to the reference time as 0 ns after commencing simulation. The fluctuation in RMSD values for protein remain within 3.0 Å during the simulation period. The ligand fit to proteins oscillate within 1.5 Å up to 40 ns. This indicates that the complex is very much stable. The RMSD plot (Fig. 6c) indicates that stigmasta-4,22-dien-3-one -VP24 (5HLJ) complex and the protein converged at 15 ns and the RMSD values remain within 3 Å and diverged at 35 ns. The complex stabilized at 38 ns and variations were in the order of 2.5 Å indicating the complex has not undergone conformational changes. The overall RMSD plot of stigmasta-4,22-dien-3-one -VP24 complex suggest that ligand is bound stably to VP24 and has not diffused away from the protein.
3.4.2 Root Mean Square Fluctuation (RMSF)
The root mean square fluctuation (RMSF) represents the fluctuation of the protein and ligand at the atomic and residual levels. In this study, we have only analysed the protein fluctuation at the residual level in terms of residue number and interaction with stigmasta-4,22-dien-3-one. The residues of VP28 have performed steadily that only a few residues have fluctuated beyond 2 Å and those residues are THR32, VAL33, THR34, LYS35 GLY200 and THR201. The ligand molecule interacted with the protein 15 times involving VAL140, TYR193, ASN51, THR135, VAL140, PHE158, TYR193, HIS195, ASN51, LEU52, VAL140, PHE158, TYR193, VAL194 and HIS195 (Fig. 7a). The VP26 has performed stable during the simulation period and the ligand atoms interacted with the protein 18 times involving VAL24, ALA44, TYR46, PRO72, LEU73, THR124, GL125, VAL126, ASP127, ASN193, VAL194, ILE195, ASP196, ILE197, LYS198, ASP199, GLU200 and ILE201 (Fig. 7b). The VP24 residues fluctuated slightly more than 2 Å and those residues are PHE203, LYS204, VA205, GLY206, GLU207 and LYS208. The ligand molecule interacted with the protein 10 times THR108, LEU109, ASN110, VAL111, ASN112, GLU118, ASP136, VAL138, PHE156, VAL157, LEU158, TYR159, PHE160, LYS161 and THR178 (Fig. 7c). In the overall protein-RMSF analysis some residues showed slight fluctuations and the interactions were not observed to be specific to the range of fluctuated amino acids and the width of the line (green colour) indicates the duration of the interaction’s duration. The stigmasta-4,22-dien-3-one fit appropriately in the pocket of the proteins and destabilizes its core functionalities to inhibit the viral growth.