Pharmacokinetics Profiler (PhaKinPro): Model Development, Validation, and Implementation as a Web Tool for Triaging Compounds with Undesired Pharmacokinetics Profiles

Computational models that predict pharmacokinetic properties are critical to deprioritize drug candidates that emerge as hits in high-throughput screening campaigns. We collected, curated, and integrated a database of compounds tested in 12 major end points comprising over 10,000 unique molecules. We then employed these data to build and validate binary quantitative structure-activity relationship (QSAR) models. All trained models achieved a correct classification rate above 0.60 and a positive predictive value above 0.50. To illustrate their utility in drug discovery, we used these models to predict the pharmacokinetic properties for drugs in the NCATS Inxight Drugs database. In addition, we employed the developed models to predict the pharmacokinetic properties of all compounds in the DrugBank. All models described in this paper have been integrated and made publicly available via the PhaKinPro Web-portal that can be accessed at https://phakinpro.mml.unc.edu/.

Development, Validation, and Application of a Human Reproductive Toxicity Prediction Model Based on Adverse Outcome Pathway

A growing number of environmental contaminants have been proved to have reproductive toxicity to males and females. However, the unclear toxicological mechanism of reproductive toxicants limits the development of virtual screening methods. By consolidating androgen (AR)-/estrogen receptors (ERs)-mediated adverse outcome pathways (AOPs) with more than 8000 chemical substances, we uncovered relationships between chemical features, a series of pathway-related effects, and reproductive apical outcomes─changes in sex organ weights. An AOP-based computational model named RepTox was developed and evaluated to predict and characterize chemicals' reproductive toxicity for males and females. Results showed that RepTox has three outstanding advantages. (I) Compared with the traditional models (37 and 81% accuracy, respectively), AOP significantly improved the predictive robustness of RepTox (96.3% accuracy). (II) Compared with the application domain (AD) of models based on small in vivo datasets, AOP expanded the ADs of RepTox by 1.65-fold for male and 3.77-fold for female, respectively. (III) RepTox implied that hydrophobicity, cyclopentanol substructure, and several topological indices (e.g., hydrogen-bond acceptors) were important, unbiased features associated with reproductive toxicants. Finally, RepTox was applied to the inventory of existing chemical substances of China and identified 2100 and 7281 potential toxicants to the male and female reproductive systems, respectively.

Simplex representation of molecular structure as universal QSAR/QSPR tool

We review the development and application of the Simplex approach for the solution of various QSAR/QSPR problems. The general concept of the simplex method and its varieties are described. The advantages of utilizing this methodology, especially for the interpretation of QSAR/QSPR models, are presented in comparison to other fragmentary methods of molecular structure representation. The utility of SiRMS is demonstrated not only in the standard QSAR/QSPR applications, but also for mixtures, polymers, materials, and other complex systems. In addition to many different types of biological activity (antiviral, antimicrobial, antitumor, psychotropic, analgesic, etc.), toxicity and bioavailability, the review examines the simulation of important properties, such as water solubility, lipophilicity, as well as luminescence, and thermodynamic properties (melting and boiling temperatures, critical parameters, etc.). This review focuses on the stereochemical description of molecules within the simplex approach and details the possibilities of universal molecular stereo-analysis and stereochemical configuration description, along with stereo-isomerization mechanism and molecular fragment "topography" identification.

Two Decades of 4D-QSAR: A Dying Art or Staging a Comeback?

A key question confronting computational chemists concerns the preferable ligand geometry that fits complementarily into the receptor pocket. Typically, the postulated ‘bioactive’ 3D ligand conformation is constructed as a ‘sophisticated guess’ (unnecessarily geometry-optimized) mirroring the pharmacophore hypothesis—sometimes based on an erroneous prerequisite. Hence, 4D-QSAR scheme and its ‘dialects’ have been practically implemented as higher level of model abstraction that allows the examination of the multiple molecular conformation, orientation and protonation representation, respectively. Nearly a quarter of a century has passed since the eminent work of Hopfinger appeared on the stage; therefore the natural question occurs whether 4D-QSAR approach is still appealing to the scientific community? With no intention to be comprehensive, a review of the current state of art in the field of receptor-independent (RI) and receptor-dependent (RD) 4D-QSAR methodology is provided with a brief examination of the ‘mainstream’ algorithms. In fact, a myriad of 4D-QSAR methods have been implemented and applied practically for a diverse range of molecules. It seems that, 4D-QSAR approach has been experiencing a promising renaissance of interests that might be fuelled by the rising power of the graphics processing unit (GPU) clusters applied to full-atom MD-based simulations of the protein-ligand complexes.

Pred-Skin: A Web Portal for Accurate Prediction of Human Skin Sensitizers

Safety assessment is an essential component of the regulatory acceptance of industrial chemicals. Previously, we have developed a model to predict the skin sensitization potential of chemicals for two assays, the human patch test and murine local lymph node assay, and implemented this model in a web portal. Here, we report on the substantially revised and expanded freely available web tool, Pred-Skin version 3.0. This up-to-date version of Pred-Skin incorporates multiple quantitative structure-activity relationship (QSAR) models developed with in vitro, in chemico, and mice and human in vivo data, integrated into a consensus naïve Bayes model that predicts human effects. Individual QSAR models were generated using skin sensitization data derived from human repeat insult patch tests, human maximization tests, and mouse local lymph node assays. In addition, data for three validated alternative methods, the direct peptide reactivity assay, KeratinoSens, and the human cell line activation test, were employed as well. Models were developed using open-source tools and rigorously validated according to the best practices of QSAR modeling. Predictions obtained from these models were then used to build a naïve Bayes model for predicting human skin sensitization with the following external prediction accuracy: correct classification rate (89%), sensitivity (94%), positive predicted value (91%), specificity (84%), and negative predicted value (89%). As an additional assessment of model performance, we identified 11 cosmetic ingredients known to cause skin sensitization but were not included in our training set, and nine of them were accurately predicted as sensitizers by our models. Pred-Skin can be used as a reliable alternative to animal tests for predicting human skin sensitization.

4D- quantitative structure-activity relationship modeling: making a comeback

Introduction: Predictive Quantitative Structure-Activity Relationship (QSAR) modeling has become an essential methodology for rapidly assessing various properties of chemicals. The vast majority of these QSAR models utilize numerical descriptors derived from the two- and/or three-dimensional structures of molecules. However, the conformation-dependent characteristics of flexible molecules and their dynamic interactions with biological target(s) is/are not encoded by these descriptors, leading to limited prediction performances and reduced interpretability. 2D/3D QSAR models are successful for virtual screening, but typically suffer at lead optimization stages. That is why conformation-dependent 4D-QSAR modeling methods were developed two decades ago. However, these methods have always suffered from the associated computational cost. Recently, 4D-QSAR has been experiencing a significant come-back due to rapid advances in GPU-accelerated molecular dynamic simulations and modern machine learning techniques. Areas covered: Herein, the authors briefly review the literature regarding 4D-QSAR modeling and describe its modern workflow called MD-QSAR. Challenges and current limitations are also highlighted. Expert opinion: The development of hyper-predictive MD-QSAR models could represent a disruptive technology for analyzing, understanding, and optimizing dynamic protein-ligand interactions with countless applications for drug discovery and chemical toxicity assessment. Therefore, there has never been a better time and relevance for molecular modeling teams to engage in hyper-predictive MD-QSAR modeling.

Consensus QSAR modelling of SIRT1 activators using simplex representation of molecular structure

Hierarchical QSAR technology (HiT QSAR) was used for consensus QSAR modelling of 65 SIRT1 activators. Simplex representation of molecular structure (SiRMS) has been used for descriptor generation. The predictive QSAR models were developed using the partial least squares (PLS) method. The QSAR models were built up according to OECD principles. One hundred rounds of Y-scrambling were performed for each selected model to exclude chance correlations. A successful consensus model (r² = 0.830, [Formula: see text] = 0.754) was obtained from the five best QSAR models. Leverage, ellipsoid and local tree domain of applicability (DA) approaches have been used for evaluation of the quality of predictions. Molecular fragments responsible for an increase and decrease of the activation properties have been determined by mechanistic interpretation of the developed QSAR model.

xMaP-An Interpretable Alignment-Free Four-Dimensional Quantitative Structure-Activity Relationship Technique Based on Molecular Surface Properties and Conformer Ensembles

A novel alignment-free molecular descriptor called xMaP (flexible MaP descriptor) is introduced. The descriptor is the advancement of the previously published translationally and rotationally invariant three-dimensional (3D) descriptor MaP (mapping property distributions onto the molecular surface) to the fourth dimension (4D). In addition to MaP, xMaP is independent of the chosen starting conformation of the encoded molecules and is therefore entirely alignment-free. This is achieved by using ensembles of conformers, which are generated by conformational searches. This step of the procedure is similar to Hopfinger's 4D quantitative structure-activity relationship (QSAR). A five-step procedure is used to compute the xMaP descriptor. First, a conformational search for each molecule is carried out. Next, for each of the conformers an approximation to the molecular surface with equally distributed surface points is computed. Third, molecular properties are projected onto this surface. Fourth, areas of identical properties are clustered to so-called patches. Fifth, the spatial distribution of the patches is converted into an alignment-free descriptor that is based on the entire conformer ensemble. The resulting descriptor can be interpreted by superimposing the most important descriptor variables and the molecules of the data set. The most important descriptor variables are identified with chemometric regression tools. The novel descriptor was applied to several benchmark data sets and was compared to other descriptors and QSAR techniques comprising a binary fingerprint, a topological pharmacophore descriptor (Cats2D), and the field-based 3D-QSAR technique GRID/PLS which is alignment-dependent. The use of conformer ensembles renders xMaP very robust. It turns out that xMaP performs very well on (almost) all data sets and that the statistical results are comparable to GRID/PLS. In addition to that, xMaP can also be used to efficiently visualize the derived quantitative structure-activity relationships.

Databases and QSAR for Cancer Research

In this review, we take a survey of bioinformatics databases and quantitative structure-activity relationship studies reported in published literature. Databases from the most general to special cancer-related ones have been included. Most commonly used methods of structure-based analysis of molecules have been reviewed, along with some case studies where they have been used in cancer research. This article is expected to be of use for general bioinformatics researchers interested in cancer and will also provide an update to those who have been actively pursuing this field of research.

Alarms about structural alerts

Structural alerts are widely accepted in chemical toxicology and regulatory decision support as a simple and transparent means to flag potential chemical hazards or group compounds into categories for read-across. However, there has been a growing concern that alerts disproportionally flag too many chemicals as toxic, which questions their reliability as toxicity markers. Conversely, the rigorously developed and properly validated statistical QSAR models can accurately and reliably predict the toxicity of a chemical; however, their use in regulatory toxicology has been hampered by the lack of transparency and interpretability. We demonstrate that contrary to the common perception of QSAR models as "black boxes" they can be used to identify statistically significant chemical substructures (QSAR-based alerts) that influence toxicity. We show through several case studies, however, that the mere presence of structural alerts in a chemical, irrespective of the derivation method (expert-based or QSAR-based), should be perceived only as hypotheses of possible toxicological effect. We propose a new approach that synergistically integrates structural alerts and rigorously validated QSAR models for a more transparent and accurate safety assessment of new chemicals.

QSAR modeling: where have you been? Where are you going to?

Quantitative structure-activity relationship modeling is one of the major computational tools employed in medicinal chemistry. However, throughout its entire history it has drawn both praise and criticism concerning its reliability, limitations, successes, and failures. In this paper, we discuss (i) the development and evolution of QSAR; (ii) the current trends, unsolved problems, and pressing challenges; and (iii) several novel and emerging applications of QSAR modeling. Throughout this discussion, we provide guidelines for QSAR development, validation, and application, which are summarized in best practices for building rigorously validated and externally predictive QSAR models. We hope that this Perspective will help communications between computational and experimental chemists toward collaborative development and use of QSAR models. We also believe that the guidelines presented here will help journal editors and reviewers apply more stringent scientific standards to manuscripts reporting new QSAR studies, as well as encourage the use of high quality, validated QSARs for regulatory decision making.

Does rational selection of training and test sets improve the outcome of QSAR modeling?

Prior to using a quantitative structure activity relationship (QSAR) model for external predictions, its predictive power should be established and validated. In the absence of a true external data set, the best way to validate the predictive ability of a model is to perform its statistical external validation. In statistical external validation, the overall data set is divided into training and test sets. Commonly, this splitting is performed using random division. Rational splitting methods can divide data sets into training and test sets in an intelligent fashion. The purpose of this study was to determine whether rational division methods lead to more predictive models compared to random division. A special data splitting procedure was used to facilitate the comparison between random and rational division methods. For each toxicity end point, the overall data set was divided into a modeling set (80% of the overall set) and an external evaluation set (20% of the overall set) using random division. The modeling set was then subdivided into a training set (80% of the modeling set) and a test set (20% of the modeling set) using rational division methods and by using random division. The Kennard-Stone, minimal test set dissimilarity, and sphere exclusion algorithms were used as the rational division methods. The hierarchical clustering, random forest, and k-nearest neighbor (kNN) methods were used to develop QSAR models based on the training sets. For kNN QSAR, multiple training and test sets were generated, and multiple QSAR models were built. The results of this study indicate that models based on rational division methods generate better statistical results for the test sets than models based on random division, but the predictive power of both types of models are comparable.

Existing and Developing Approaches for QSAR Analysis of Mixtures

This review is devoted to the critical analysis of advantages and disadvantages of existing mixture descriptors and their usage in various QSAR/QSPR tasks. We describe good practices for the QSAR modeling of mixtures, data sources for mixtures, a discussion of various mixture descriptors and their application, recommendations about proper external validation specific for mixture QSAR modeling, and future perspectives of this field. The biggest problem in QSAR of mixtures is the lack of reliable data about the mixtures' properties. Various mixture descriptors are used for the modeling of different endpoints. However, these descriptors have certain disadvantages, such as applicability only to 1 : 1 binary mixtures, and additive nature. The field of QSAR of mixtures is still under development, and existing efforts could be considered as a foundation for future approaches and studies. The usage of non-additive mixture descriptors, which are sensitive to interaction effects, in combination with best practices of QSAR model development (e.g., thorough data collection and curation, rigorous external validation, etc.) will significantly improve the quality of QSAR studies of mixtures.

QSAR analysis of [(biphenyloxy)propyl]isoxazoles: agents against coxsackievirus B3

BACKGROUND

Antiviral drugs are urgently needed for the treatment of acute and chronic diseases caused by enteroviruses such as coxsackievirus B3 (CVB3). The main goal of this study is quantitative structure-activity relationship (QSAR) analysis of anti-CVB3 activity (clinical CVB3 isolate 97927 [log IC50, µM]) and investigation of the selectivity of 25 ([biphenyloxy]propyl)isoxazoles, followed by computer-aided design and virtual screening of novel active compounds.

DISCUSSION

The 2D QSAR obtained models are quite satisfactory (R(2) = 0.84-0.99, Q(2) = 0.76-0.92, R(2)(ext) = 0.62-0.79). Compounds with high antiviral activity and selectivity have to contain 5-trifluoromethyl-[1,2,4]oxadiazole or 2,4-difluorophenyl fragments. Insertion of 2,5-dimethylbenzene, napthyl and especially biphenyl substituents into investigated compounds substantially decreases both their antiviral activity and selectivity. Several compounds were proposed as a result of design and virtual screening. A high level of activity of 2-methoxy-1-phenyl-1H-imidazo[4,5-c]pyridine (sm428) was confirmed experimentally.

CONCLUSION

Simplex representation of molecular structure allows successful QSAR analysis of anti-CVB3 activity of ([biphenyloxy]propyl)isoxazole derivatives. Two possible ways of battling CVB3 are considered as a future perspective.

QSAR analysis of the toxicity of nitroaromatics in Tetrahymena pyriformis: structural factors and possible modes of action

The Hierarchical Technology for Quantitative Structure-Activity Relationships (HiT QSAR) was applied to 95 diverse nitroaromatic compounds (including some widely known explosives) tested for their toxicity (50% inhibition growth concentration, IGC₅₀) against the ciliate Tetrahymena pyriformis. The dataset was divided into subsets according to putative mechanisms of toxicity. The Classification and Regression Trees (CART) approach implemented within HiT QSAR has been used for prediction of mechanism of toxicity for new compounds. The resulting models were shown to have ~80% accuracy for external datasets indicating that the mechanistic dataset division was sensible. The Partial Least Squares (PLS) statistical approach was then used to develop 2D QSAR models. Validated PLS models were explored to: (1) elucidate the effects of different substituents in nitroaromatic compounds on toxicity; (2) differentiate compounds by probable mechanisms of toxicity based on their structural descriptors; and (3) analyse the role of various physical-chemical factors responsible for compounds' toxicity. Models were interpreted in terms of molecular fragments promoting or interfering with toxicity. It was also shown that mutual influence of substituents in benzene ring plays the determining role in toxicity variation. Although chemical mechanism based models were statistically significant and externally predictive (r²(ext) = 0.64 for the external set of 63 nitroaromatics identified after all calculations have been completed), they were also shown to have limited coverage (57% for modelling and 76% for external set).

Per aspera ad astra: application of Simplex QSAR approach in antiviral research

This review explores the application of the Simplex representation of molecular structure (SiRMS) QSAR approach in antiviral research. We provide an introduction to and description of SiRMS, its application in antiviral research and future directions of development of the Simplex approach and the whole QSAR field. In the Simplex approach every molecule is represented as a system of different simplexes (tetratomic fragments with fixed composition, structure, chirality and symmetry). The main advantages of SiRMS are consideration of the different physical–chemical properties of atoms, high adequacy and good interpretability of models obtained and clear procedures for molecular design. The reliability of developed QSAR models as predictive virtual screening tools and their ability to serve as the basis of directed drug design was validated by subsequent synthetic and biological experiments. The SiRMS approach is realized as the complex of the computer program ‘HiT QSAR’, which is available on request.

QSPR Analysis of Peroxidase Substrates Reactivity

Quantitative structure-property relationship (QSPR) analysis of phenol derivatives reactivity in the horseradish peroxidase catalyzed oxidative reactions was carried out. The statistic models, which describe the substituted phenols reactivity (Кm-1, Vmax) quite adequately, were obtained by multiple linear regression and partial least squares (PLS) methods. The electronic parameters of molecules, their lipophylicity, molecular refraction, and form parameters were used as descriptors for molecular structure. The obtained models allow to predict the reactivity of the new phenolic substrates with satisfactory reliability.

Virtual Screening and Molecular Design Based on Hierarchical Qsar Technology

This chapter is devoted to the hierarchical QSAR technology (HiT QSAR) based on simplex representation of molecular structure (SiRMS) and its application to different QSAR/QSPR tasks. The essence of this technology is a sequential solution (with the use of the information obtained on the previous steps) of the QSAR paradigm by a series of enhanced models based on molecular structure description (in a specific order from 1D to 4D). Actually, it’s a system of permanently improved solutions. Different approaches for domain applicability estimation are implemented in HiT QSAR. In the SiRMS approach every molecule is represented as a system of different simplexes (tetratomic fragments with fixed composition, structure, chirality, and symmetry). The level of simplex descriptors detailed increases consecutively from the 1D to 4D representation of the molecular structure. The advantages of the approach presented are an ability to solve QSAR/QSPR tasks for mixtures of compounds, the absence of the “molecular alignment” problem, consideration of different physical–chemical properties of atoms (e.g., charge, lipophilicity), and the high adequacy and good interpretability of obtained models and clear ways for molecular design. The efficiency of HiT QSAR was demonstrated by its comparison with the most popular modern QSAR approaches on two representative examination sets. The examples of successful application of the HiT QSAR for various QSAR/QSPR investigations on the different levels (1D–4D) of the molecular structure description are also highlighted. The reliability of developed QSAR models as the predictive virtual screening tools and their ability to serve as the basis of directed drug design was validated by subsequent synthetic, biological, etc. experiments. The HiT QSAR is realized as the suite of computer programs termed the “HiT QSAR” software that so includes powerful statistical capabilities and a number of useful utilities.

[Artificial ribonucleases: quantitative analysis of the structure-activity relationship and new insight into the strategy of design of highly efficient RNase mimetics]

The dependence of hydrolytic activity of artificial ribonucleases toward an HIV-I RNA fragment, a 21-mer oligonucleotide, and tRNA Asp on the structure of the RNase mimetic was analyzed. The quantitative structure-activity relationship (QSAR task) was determined by the method of simplex representation of the molecular structure where the amounts of four-atom fragments (simplexes) of fixed structure, symmetry, and chirality served as descriptors. Not only the types of atoms participating in simplexes but also their physicochemical properties (e.g., partial charges, lipophilicities, etc.) were taken into account. This allowed the estimation of the relative role of various factors affecting the interaction of molecules under study with the corresponding biological target. The 2D QSAR models obtained by the method of projection to latent structures have quite satisfactory statistical indices (R2 = 0.82-0.96; Q2 = 0.73-0.89), which help predict the activities of new compounds. The electrostatic properties of ribonuclease atoms were shown to contribute significantly to the manifestation of the hydrolytic activity of ribonucleases in the case of the 21-mer oligonucleotide and tRNA. In addition, the structural fragments that most greatly contribute to the alteration of the hydrolytic activity of RNases were identified. The models obtained were used for the virtual screening and molecular design of new highly efficient RNase mimetics.

The effects of characteristics of substituents on toxicity of the nitroaromatics: HiT QSAR study

The present study applies the Hierarchical Technology for Quantitative Structure-Activity Relationships (HiT QSAR) for (i) evaluation of the influence of the characteristics of 28 nitroaromatic compounds (some of which belong to a widely known class of explosives) as to their toxicity; (ii) prediction of toxicity for new nitroaromatic derivatives; (iii) analysis of the effects of substituents in nitroaromatic compounds on their toxicity in vivo. The 50% lethal dose concentration for rats (LD50) was used to develop the QSAR models based on simplex representation of molecular structure. The preliminary 1D QSAR results show that even the information on the composition of molecules reveals the main tendencies of changes in toxicity. The statistic characteristics for partial least squares 2D QSAR models are quite satisfactory (R2 = 0.96-0.98; Q2 = 0.91-0.93; R2 (test) = 0.89-0.92), which allows us to carry out the prediction of activity for 41 novel compounds designed by the application of new combinations of substituents represented in the training set. The comprehensive analysis of toxicity changes as a function of substituent position and nature was carried out. Molecular fragments that promote and interfere with toxicity were defined on the basis of the obtained models. It was shown that the mutual influence of substituents in the benzene ring plays a crucial role regarding toxicity. The influence of different substituents on toxicity can be mediated via different C-H fragments of the aromatic ring.

Quantitative structure-activity relationship studies of [(biphenyloxy)propyl]isoxazole derivatives. Inhibitors of human rhinovirus 2 replication

The 50% cytotoxic concentration (CC50) in HeLa cells, the 50% inhibitory concentration (IC50) against human rhinovirus 2 (HRV-2), and the selectivity index (SI = CC50/IC50) of [(biphenyloxy)propyl]isoxazole derivatives were used to develop quantitative structure-activity relationships (QSAR) based on simplex representation of molecular structure. Statistic characteristics for partial least-squares models are quite satisfactory (R2 = 0.838 - 0.918; Q2 = 0.695 - 0.87) for prediction of CC50, IC50, and SI values and permit the virtual screening and molecular design of new compounds with strong anti-HRV-2 activity. The quality of prognosis for designed compounds was additionally estimated by analysis of domain applicability for each QSAR model. A hypothesis to the effect that terminal benzene substituents must have negative electrostatic potential and definite length (approximately 5.5-5.6 A) to possess strong antiviral activity has been suggested. The quality of developed analysis, i.e., high level of antiviral action of three new designed compounds, has been confirmed experimentally.

Identification of individual structural fragments of N,N'-(bis-5-nitropyrimidyl)dispirotripiperazine derivatives for cytotoxicity and antiherpetic activity allows the prediction of new highly active compounds

OBJECTIVES

The objectives of this study were (i) to apply computer-based technologies to evaluate the structure of 48 N,N'-(bis-5-nitropyrimidyl)dispirotripiperazines which belong to a new class of highly active antiviral compounds binding to cell surface heparan sulphates, (ii) to understand the chemical- biological interactions governing their activities, and (iii) to design new compounds with strong antiviral activity.

METHODS

The logarithm of 50% cytotoxic concentration (CC(50)) in GMK cells, of 50% inhibitory concentration (IC(50)) against herpes simplex virus type 1, and of selectivity index (SI = CC(50)/IC(50)) was used to develop quantitative structure-activity relationships (QSARs) based on simplex representation of molecular structure. The QSAR model was applied to design new compounds. Two of these compounds were synthesized, physico-chemically characterized and tested for cytotoxicity and antiviral activity.

RESULTS

Statistic characteristics for partial least squares models allow the prediction of CC(50), IC(50) and SI values. The QSAR results demonstrate a high impact of individual structural fragments for antiviral activity. Molecular fragments that promote and interfere with antiviral activity were defined on the basis of the obtained models. Electrostatic factors (38%) and hydrophobicity (34%) were the most important determinants of antiherpetic activity. Using the established method, new potential dispirotripiperazine derivatives were computationally designed. Two of these computationally designed compounds were synthesized. The biological test results confirm the computationally predicted values of these compounds.

CONCLUSIONS

The established QSAR model is suitable for the design of new antiherpetic compounds and prediction of their activity.

Hierarchic system of QSAR models (1D–4D) on the base of simplex representation of molecular structure

In this work, a hierarchic system of QSAR models from 1D to 4D is considered on the basis of the simplex representation of molecular structure (SiRMS). The essence of this system is that the QSAR problem is solved sequentially in a series of the improved models of the description of molecular structure. Thus, at each subsequent stage of a hierarchic system, the QSAR problem is not solved ab ovo, but rather the information obtained from the previous step is used. Actually, we deal with a system of solutions defined more exactly. In the SiRMS approach, a molecule is represented as a system of different simplex descriptors (tetratomic fragments with fixed composition, structure, chirality and symmetry). The level of simplex-descriptor detail increases consecutively from 1D to 4D representations of molecular structure. It enables us to determine the fragments of structure that promote or interfere with the given biological activity easily. Molecular design of compounds with a given level of activity is possible on the basis of SiRMS. The efficiency of the method is demonstrated for the example of the analysis of substituted piperazines affinity for the 5-HT1A receptor.