The prediction of main protease SARS-CoV-2 inhibition based on models of enzyme-inhibitor complexes

A set of linear regression equations predicting the IC50 values for SARS-CoV-2 main protease inhibitors was analyzed. For 180 competitive inhibitors, we have simulated the molecular dynamics of enzyme-inhibitor complexes with known structures or modeled using molecular docking. In the docking procedure, the selection of final poses was restricted by similarity to known structural analogs. The values of the energy contributions obtained by means of calculation of the free energy change of the enzyme-inhibitor complex performed by two variants of the MMPBSA (MMGBSA) method and a number of physicochemical characteristics of the inhibitors were used as independent variables. During the learning process, indicator variables were used for inhibitor subsets obtained from various literature sources to compensate the existing systematic deviations from the target value. A leave one out and leave 20% out cross validation procedures were used to evaluate the prediction quality. For the total logarithmic range width of 3.71, the mean error in predicting the lg(IC50) value was 0.45 log units. The stability of the prediction depending on the variability of the complex in molecular dynamics was investigated.