Qi-Shi Du

Virtual Screening for SARS-CoV Protease Based on KZ7088 Pharmacophore Points †

Pharmacophore modeling can provide valuable insight into ligand-receptor interactions. It can also be used in 3D (dimensional) database searching for potentially finding biologically active compounds and providing new research ideas and directions for drug-discovery projects. To stimulate the structure-based drug design against SARS (severe acute respiratory syndrome), a pharmacophore search was conducted over 3.6 millions of compounds based on the atomic coordinates of the complex obtained by docking KZ7088 (a derivative of AG7088) to SARS CoV M(pro) (coronavirus main proteinase), as reportedly recently (Chou, K. C.; Wei, D. Q.; Zhong, W. Z. Biochem. Biophys. Res. Commun. 2003, 308, 148-151). It has been found that, of the 3.6 millions of compounds screened, 0.07% are with the score satisfying five of the six pharmacophore points. Moreover, each of the hit compounds has been evaluated for druggability according to 13 metrics based on physical, chemical, and structural properties. Of the 0.07% compounds thus retrieved, 17% have a perfect score of 1.0; while 23% with one druggable rule violation, 13% two violations, and 47% more than two violations. If the criterion for druggability is set at a maximum allowance of two rule violations, we obtain that only about 0.03% of the compounds screened are worthy of further tests by experiments. These findings will significantly narrow down the search scope for potential compounds, saving substantial time and money. Finally, the featured templates derived from the current study will also be very useful for guiding the design and synthesis of effective drugs for SARS therapy.

Progress in computational approach to drug development against SARS.

Since the outbreak of SARS (severe acute respiratory syndrome) in November 2002 in Southern Chinas Guangdong Province, considerable progress has been made in the development of drugs for SARS therapy. The present mini review is focused on the area of computer-aided drug discovery, i.e., the advances achieved mainly from the approaches of structural bioinformatics, pharmacophore modeling, molecular docking, peptide-cleavage site prediction, and other computational means. It is highlighted that the compounds C(28)H(34)O(4)N(7)Cl, C(21)H(36(O)5)N(6) and C(21)H(36)O(5)N(6) (Wei et al., Amino Acids, 2006, 31: 73-80), as well as KZ7088 (Chou et al. Biochem. Biophys. Res. Commun., 2003, 308: 148-151), a derivative of AG7088, might be the promising candidates for further investigation, and that the octapeptides ATLQAIAS and ATLQAENV, as well as AVLQSGFR, might be converted to effective inhibitors against the SARS enzyme. Meanwhile, how to modify these octapeptides based on the distorted key theory to make them become potent inhibitors is explicitly elucidated. Finally, a brief introduction is given for how to use computer-generated graphs to rapidly diagnose SARS coronavirus.

Molecular modeling and chemical modification for finding peptide inhibitor against severe acute respiratory syndrome coronavirus main proteinase

n Abstractn n Severe acute respiratory syndrome (SARS) is a respiratory disease caused by a newly found virus, called SARS coronavirus. In this study, the cleavage mechanism of the SARS coronavirus main proteinase (Mpro or 3CLpro) on the octapeptide NH2-AVLQ↓SGFR-COOH was investigated using molecular mechanics and quantum mechanics simulations based on the experimental structure of the proteinase. It has been observed that the catalytic dyad (His-41/Cys-145) site between domains I and II attracts the π electron density from the peptide bond Gln–Ser, increasing the positive charge on C(CO) of Gln and the negative charge on N(NH) of Ser, so as to weaken the Gln–Ser peptide bond. The catalytic functional group is the imidazole group of His-41 and the S in Cys-145. Nδ1 on the imidazole ring plays the acid–base catalytic role. Based on the “distorted key theory” [K.C. Chou, Anal. Biochem. 233 (1996) 1–14], the possibility to convert the octapeptide to a competent inhibitor has been studied. It has been found that the chemical bond between Gln and Ser will become much stronger and no longer cleavable by the SARS enzyme after either changing the carbonyl group CO of Gln to CH2 or CF2 or changing the NH of Ser to CH2 or CF2. The octapeptide thus modified might become an effective inhibitor or a potential drug candidate against SARS.n n

Assessment of chemical libraries for their druggability

High throughput virtual screening is acknowledged as the initial means for identifying hit compounds that will be eventually transformed to leads or drug candidates. To improve quality of screening, it is essential to have powerful methods for the analysis of the compound databases. For this purpose, we have developed a novel and practical scoring function to assess the druggability of compounds. The proposed function consists of 12 metrics that take into account physical, chemical and structural properties as well as the presence of undesirable functional groups. We have applied this 12-metric scoring function to 44 different databases that include more than 3.8 million compounds, which are commercially available. The overall quality of each database was evaluated according to the score and rank measured by our 12-metric function. Our findings suggest that, the majority of compounds that do not satisfy druggable rules do so due to high molecular weight, high logP values and the presence of reactive functional groups.

Inhibitor Design for SARS Coronavirus Main Protease Based on “Distorted Key Theory”

In order to find effective peptide inhibitors against SARS CoV M(pro), an analysis was performed for 11 oligo-peptides that can be cleaved by the SARS coronavirus main protease (CoV M(pro), or 3CL(pro)). Flexible molecular alignments of the 11 cleavable peptides have provided useful insights into the chemical properties of their amino acid residues close to the cleavage site. Moreover, it was found through the ligand-receptor docking studies that of the 11 cleavable peptides, NH2-ATLQ / AIAS-COOH and NH2-ATLQ / AENV-COOH had the highest affinity with SARS CoV M(pro). The two octapeptides were selected as initial templates for further chemical modification to make them become effective inhibitors against the SARS enzyme according to the distorted key theory [K. C. Chou, Analytical Biochemistry 233 (1996) 1-14]. The possible chemical modification methods are proposed and examined. The approach developed in this study and the findings thus obtained might stimulate new strategies and provide useful information for drug design against SARS.

Virtual screening for finding natural inhibitor against cathepsin-L for SARS therapy

Summary.Recently Simmons et al. reported a new mechanism for SARS virus entry into target cells, where MDL28170 was identified as an efficient inhibitor of CTSL-meditated substrate cleavage with IC50 of 2.5u2009nmol/l. Based on the molecule fingerprint searching method, 11 natural molecules were found in the Traditional Chinese Medicines Database (TCMD). Molecular simulation indicates that the MOL376 (a compound derived from a Chinese medicine herb with the therapeutic efficacy on the human body such as relieving cough, removing the phlegm, and relieving asthma) has not only the highest binding energy with the receptor but also the good match in geometric conformation. It was observed through docking studies that the van der Waals interactions made substantial contributions to the affinity, and that the receptor active pocket was too large for MDL21870 but more suitable for MOL736. Accordingly, MOL736 might possibly become a promising lead compound for CTSL inhibition for SARS therapy.

Correlations of amino acids in proteins

A correlation analysis among 20 amino acids is performed for four protein structural classes (alpha, beta, alpha/beta, and alpha+beta) in a total of 204 proteins. The correlation relationships among amino acids can be classified into the following four types: (1) strong positive correlation, (2) strong negative correlation, (3) weak correlation, and (4) no correlation. The correlation relationships are different for different proteins and are correlated with the features of their structural classes. The amino acids with the weak correlation relationship can be treated as the independent basis functions for the space where proteins are defined. The amino acids with large correlation coefficients are linear correlative with each other and they are not independent. The strong correlation among amino acids reflects their mutual constrained relationship, as exhibited by their relevant structural features. The information obtained through the correlation analysis is used for predicting protein structural classes and a better prediction quality is obtained than that by the simple geometry distance methods without taking into account the correlation effects.

Polyprotein cleavage mechanism of SARS CoV Mpro and chemical modification of the octapeptide

n Abstractn n The cleavage mechanism of severe acute respiratory syndrome (SARS) coronavirus main proteinase (Mpro or 3CLpro) for the octapeptide AVLQSGFR is studied using molecular mechanics (MM) and quantum mechanics (QM). The catalytic dyad His-41 and Cys-145 in the active pocket between domain I and II seem to polarize the π-electron density of the peptide bond between Gln and Ser in the octapeptide, leading to an increase of positive charge on C(CO) of Gln and negative charge on N(NH) of Ser. The possibility of enhancing the chemical bond between Gln and Ser based on the “distorted key” theory [Anal. Biochem. 233 (1996) 1] is examined. The scissile peptide bond between Gln and Ser is found to be solidified through “hybrid peptide bond” by changing the carbonyl group CO of Gln to CH2 or CF2. This leads to a break of the π-bond system for the peptide bond, making the octapeptide (AVLQSGFR) a “distorted key” and a potential starting system for the design of anti SARS drugs.n n

Application of Bioinformatics in Search for Cleavable Peptides of SARSCoV Mpro and Chemical Modification of Octapeptides

According to the distorted key theory as elaborated in a review article years ago (Chou, K.C.: Analytical Biochemistry, 1996, 233, 1-14), the knowledge of the cleavable peptides by SARS-CoV M(pro) (severe acute respiratory syndrome coronavirus main proteinase) can provide very useful insights on developing drugs against SARS. In view of this, the softwares, ZCURVE_CoV 1.0 and ZCURVE_CoV 2.0 (http://tubic.tju.edu.cn/sars/), developed recently for SARS-Coronavirus are used to analyze the 36 complete SARS-Coronavirus RNA sequences in the gene bank NCBI (http://www.ncbi.nlm.nih.gov/) from different sources for protein coding genes, and to search for the cleavage sites of SARS-CoV M(pro) in polyproteins pp1a and pp1ab. A total of 396 cleavage points are found in the 36 SARS-Coronavirus and 11 cleavable octapeptides abstracted from the 396 cleavage sites. The statistical distributions of amino acids for the cleavable octapeptides at the subsites R4, R3, R2, R1, R1, R2, R3 and R4 are calculated. The cleavage-specific positions are on R2, R1 and R1, and the positions R3 and R4 are featured by some certain specificity for SARS-CoV M(pro). The structural characters of amino acid residues around the cleavage-specific positions are discussed. Two most promising octapeptides, i.e., NH(2)-ATLQ downward arrowAIAS-COOH and NH(2)-ATLQ downward arrowAENV-COOH, are selected to be the candidates for chemical modification, converting into the inhibitors of SARS-CoV M(pro). A possible strategy to convert a cleavable octapeptide by SARS enzyme into a drug candidate against SARS is elucidated.