Background:Severe acute respiratory syndrome coronavirus (SARS-CoVs) have emerged as a global health threat, which had caused a high rate of mortality. There is an urgent need to find effective drugs against these viruses.Objective:This study aims to predict the activity of unsymmetrical aromatic disulfides by constructing a QSAR model, and to design new compounds according to the structural and physicochemical attributes responsible for higher activity towards SARS-CoVs main protease.Methods:All molecules were constructed in ChemOffice software and molecular descriptors were calculated by CODESSA software. A regression-based linear heuristic method was established by changing descriptors datasets and calculating predicted IC50values of compounds. Then, some new compounds were designed according to molecular descriptors from the heuristic method model. The compounds with predicted values smaller than a set point were constantly screened out. Finally, the properties analysis and molecular docking were conducted to further understand the structure-activity relationships of these finalized compounds.Results:The heuristic method explored the various descriptors responsible for bioactivity and gained the best linear model with R20.87. The success of the model fully passed the testing set validation, proving that the model has both high statistical significance and excellent predictive ability. A total of 5 compounds with ideal predicted IC50were found from the 96 newly designed derivatives and their properties analyze was carried out. Molecular docking experiments were conducted for the optimal compound 31a, which has the best compound activity with good target protein binding capability.Conclusion:The heuristic method was quite reliable for predicting IC50values of unsymmetrical aromatic disulfides. The present research provides meaningful guidance for further exploration of the highly active inhibitors for SARS-CoVs.