Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors - PubMed (original) (raw)

Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors

Veena Nukoolkarn et al. J Theor Biol. 2008.

Abstract

Since the emergence of the severe acute respiratory syndrome (SARS) to date, neither an effective antiviral drug nor a vaccine against SARS is available. However, it was found that a mixture of two HIV-1 proteinase inhibitors, lopinavir and ritonavir, exhibited some signs of effectiveness against the SARS virus. To understand the fine details of the molecular interactions between these proteinase inhibitors and the SARS virus via complexation, molecular dynamics simulations were carried out for the SARS-CoV 3CL(pro) free enzyme (free SARS) and its complexes with lopinavir (SARS-LPV) and ritonavir (SARS-RTV). The results show that flap closing was clearly observed when the inhibitors bind to the active site of SARS-CoV 3CL(pro). The binding affinities of LPV and RTV to SARS-CoV 3CL(pro) do not show any significant difference. In addition, six hydrogen bonds were detected in the SARS-LPV system, while seven hydrogen bonds were found in SARS-RTV complex.

PubMed Disclaimer

Figures

Fig. 1

Fig. 1

Chemical formula of HIV-1 PR inhibitor, ritonavir and lopinavir, which exhibits signs of effectives against SARS.

Fig. 2

Fig. 2

Structures of (a) SARS-free enzyme, (b) SARS–lopinavir complex and (c) SARS–ritronavir complex.

Fig. 3

Fig. 3

RMSDs of SARS-CoV 3CLpro free enzyme ((a) SARS) and its complexes with lopinavir ((b) SARS–LPV) and ritonavir ((c) SARS–RTV).

Fig. 4

Fig. 4

The probability distributions of the associated interatomic distances between the distances (a) _d_1 and (b) _d_2 from the center of mass of the following residues, H41 to C145 and D48 to Q189, respectively.

Fig. 5

Fig. 5

(a) Definitions and (b) probability distributions of the torsion angles of the two catalytic residues (Tor1: CA-CB-CG-CD of H41 and TorR2: N-CA-CB-S of C145) and two amino acids at the tip of the loop (Tor3: CA-CB-CG-OD of D48 and Tor4: CB-CG-CD-CA of Q189) during 600–2000 ps of the simulations.

Fig. 6

Fig. 6

Radial distribution functions, g(r), centered on the inhibitor atoms (see Fig. 1 for definitions) to oxygen atoms of modeled water of the two complexes, SARS–LPV and SARS–RTV. The plots are classified into two categories where the central atom types for the two inhibitors are topologically (a) equivalent and (b) different.

Similar articles

Cited by

References

    1. Anand K., Palm G.J., Mesters J.R., Siddell S.G., Ziebuhr J., Hilgenfeld R. Structure of coronavirus main proteinase reveals combination of a chymotrypsin fold with an a-extra helical domain. EMBO J. 2002;21:3213–3224. - PMC - PubMed
    1. Anand K., Ziebuhr J., Wadhwani P., Mesters J.R., Hilgenfeld R. Coronavirus main protease (3CLpro) structure: basis for design of anti-SARS drugs. Science. 2003;300:1763–1767. - PubMed
    1. Carlson H.A. Protein flexibility and drug design: how to hit a moving target. Curr. Opin. Chem. Biol. 2002;6:447–452. - PubMed
    1. Carlson H.A., McCammon J.A. Accommodating protein flexibility in computational drug design. Mol. Pharmacol. 2000;57:213–218. - PubMed
    1. Case D.A., Pearlman J.W., Caldwell T.E., Cheatham J., Jr., Wang W.S., Ross C.L., Simmerling T.A., Darden K.M., Merz R.V., Stanton A.L., Cheng J.J., Vincent M., Crowley V., Tsui H., Gohlke R.J., Radmer Y., Duan J., Pitera I., Massova G.L., Seibel U.C., Singh P.K., Kollman P.A. University of California; San Francisco, CA: 2002. AMBER. (Version 7.0 ed.)

Publication types

MeSH terms

Substances

LinkOut - more resources