Conformational coverage by a genetic algorithm: saturation of conformational space

Combined 4D-QSAR and target-based approaches for the determination of bioactive Isatin derivatives

The hybrid method of the Electron-Conformational Genetic Algorithm (EC-GA) was used to determine the pharmacophore groups and to estimate anticancer activity in isatin derivatives using a robust 4D-QSAR software (EMRE). To build the model, each compound is represented by a set of conformers rather than a single conformation. The Electron Conformational Matrix of Congruity (ECMC) is composed via EMRE software. Electron Conformational Submatrix of Activity (ECSA) was calculated by the comparison of these matrices. Genetic algorithm was used to select important variables to predict theoretical activity. The model with the best seven parameters produced satisfactory results. The E statistics technique was applied to the generated EC-GA model to evaluate the individual contribution of each of the descriptors on biological activity. The r² and q² values of the training set compounds were found to be 0.95 and 0.93, respectively. Because no previous 4D-QSAR studies on isatin derivatives have been conducted, this study is important in the development of new isatin derivatives. In this study, 27 isatin derivatives whose activities were estimated using the hybrid EC-GA method were also investigated through molecular docking and molecular dynamics simulations for their BCL-2 inhibitory activity.

Structural characterisation of natural products by means of quantum chemical calculations of NMR parameters: new insights

In this review, we focus in all aspects of NMR simulation of natural products, from the fundamentals to the new computational toolboxes available, combining advanced quantum chemical calculations with upstream data processing and machine learning.

Automated exploration of the low-energy chemical space with fast quantum chemical methods

We propose and discuss an efficient scheme for the in silico sampling for parts of the molecular chemical space by semiempirical tight-binding methods combined with a meta-dynamics driven search algorithm. The focus of this work is set on the generation of proper thermodynamic ensembles at a quantum chemical level for conformers, but similar procedures for protonation states, tautomerism and non-covalent complex geometries are also discussed. The conformational ensembles consisting of all significantly populated minimum energy structures normally form the basis of further, mostly DFT computational work, such as the calculation of spectra or macroscopic properties. By using basic quantum chemical methods, electronic effects or possible bond breaking/formation are accounted for and a very reasonable initial energetic ranking of the candidate structures is obtained. Due to the huge computational speedup gained by the fast low-cost quantum chemical methods, overall short computation times even for systems with hundreds of atoms (typically drug-sized molecules) are achieved. Furthermore, specialized applications, such as sampling with implicit solvation models or constrained conformational sampling for transition-states, metal-, surface-, or noncovalently bound complexes are discussed, opening many possible applications in modern computational chemistry and drug discovery. The procedures have been implemented in a freely available computer code called CREST, that makes use of the fast and reliable GFNn-xTB methods.

A novel hybrid method named electron conformational genetic algorithm as a 4D QSAR investigation to calculate the biological activity of the tetrahydrodibenzazosines

To understand the structure-activity correlation of a group of tetrahydrodibenzazocines as inhibitors of 17β-hydroxysteroid dehydrogenase type 3, we have performed a combined genetic algorithm (GA) and four-dimensional quantitative structure-activity relationship (4D-QSAR) modeling study. The computed electronic and geometry structure descriptors were regulated as a matrix and named as electron-conformational matrix of contiguity (ECMC). A chemical property-based pharmacophore model was developed for series of tetrahydrodibenzazocines by EMRE software package. GA was employed to choose an optimal combination of parameters. A model has been developed for estimating anticancer activity quantitatively. All QSAR models were established with 40 compounds (training set), then they were considered for selective capability with additional nine compounds (test set). A statistically valid 4D-QSAR ( R training 2 = 0.856 , R test 2 = 0.851 and q² = 0.650) with good external set prediction was obtained.

Conformation Generation: The State of the Art

The generation of conformations for small molecules is a problem of continuing interest in cheminformatics and computational drug discovery. This review will present an overview of methods used to sample conformational space, focusing on those methods designed for organic molecules commonly of interest in drug discovery. Different approaches to both the sampling of conformational space and the scoring of conformational stability will be compared and contrasted, with an emphasis on those methods suitable for conformer sampling of large numbers of drug-like molecules. Particular attention will be devoted to the appropriate utilization of information from experimental solid-state structures in validating and evaluating the performance of these tools. The review will conclude with some areas worthy of further investigation.

Discovery and Optimization of Materials Using Evolutionary Approaches

Materials science is undergoing a revolution, generating valuable new materials such as flexible solar panels, biomaterials and printable tissues, new catalysts, polymers, and porous materials with unprecedented properties. However, the number of potentially accessible materials is immense. Artificial evolutionary methods such as genetic algorithms, which explore large, complex search spaces very efficiently, can be applied to the identification and optimization of novel materials more rapidly than by physical experiments alone. Machine learning models can augment experimental measurements of materials fitness to accelerate identification of useful and novel materials in vast materials composition or property spaces. This review discusses the problems of large materials spaces, the types of evolutionary algorithms employed to identify or optimize materials, and how materials can be represented mathematically as genomes, describes fitness landscapes and mutation operators commonly employed in materials evolution, and provides a comprehensive summary of published research on the use of evolutionary methods to generate new catalysts, phosphors, and a range of other materials. The review identifies the potential for evolutionary methods to revolutionize a wide range of manufacturing, medical, and materials based industries.

Free Energy-Based Conformational Search Algorithm Using the Movable Type Sampling Method

In this article, we extend the movable type (MT) sampling method to molecular conformational searches (MT-CS) on the free energy surface of the molecule in question. Differing from traditional systematic and stochastic searching algorithms, this method uses Boltzmann energy information to facilitate the selection of the best conformations. The generated ensembles provided good coverage of the available conformational space including available crystal structures. Furthermore, our approach directly provides the solvation free energies and the relative gas and aqueous phase free energies for all generated conformers. The method is validated by a thorough analysis of thrombin ligands as well as against structures extracted from both the Protein Data Bank (PDB) and the Cambridge Structural Database (CSD). An in-depth comparison between OMEGA and MT-CS is presented to illustrate the differences between the two conformational searching strategies, i.e., energy-based versus free energy-based searching. These studies demonstrate that our MT-based ligand conformational search algorithm is a powerful approach to delineate the conformational ensembles of molecular species on free energy surfaces.

4D-QSAR study of HEPT derivatives by electron conformational-genetic algorithm method

In this work, the EC-GA method, a hybrid 4D-QSAR approach that combines the electron conformational (EC) and genetic algorithm optimization (GA) methods, was applied in order to explain pharmacophore (Pha) and predict anti-HIV-1 activity by studying 115 compounds in the class of 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio) thymine (HEPT) derivatives as non-nucleoside reverse transcriptase inhibitors (NNRTIs). The series of NNRTIs were partitioned into four training and test sets from which corresponding quantitative structure-activity relationship (QSAR) models were constructed. Analysis of the four QSAR models suggests that the three models generated from the training and test sets used in previous works yielded comparable results with those of previous studies. Model 4, the data set of which was partitioned randomly into two training and test sets with 11 descriptors, including electronical and geometrical parameters, showed good statistics both in the regression (r2(training) )= 0.867, r2test = 0.923) and cross-validation (q (2) = 0.811, q2(ext1) = 0.909, q2(ext2) = 0.909) for the training set of 80 compounds and the test set of 27 compounds. The prediction of the anti-HIV-1 activity of HEPT compounds by means of the EC-GA method allowed for a quantitatively consistent QSAR model. In addition, eight novel compounds never tested experimentally have been designed theoretically using model 4.

Androgen receptor binding affinity: a QSAR evaluation

The multiparameter formulation of the COmmon REactivity PAttern (COREPA) approach has been used to describe the structural requirements for eliciting rat androgen receptor (AR) binding affinity, accounting for molecular flexibility. Chemical affinity for AR binding was related to the distances between nucleophilic sites and structural features describing electronic and hydrophobic interactions between the receptor and ligands. Categorical models were derived for each binding affinity range in terms of specific distances, local (maximal donor delocalizability associated with the oxygen atom of the A ring), global nucleophilicity (partial positive surface areas and energy of the highest occupied molecular orbital) and hydrophobicity (log Kow) of the molecules. An integral screening tool for predicting binding affinity to AR was constructed as a battery of models, each associated with different activity bins. The quality of the screening battery of models was assessed using a high value (0.9) of the Pearson contingency coefficient. The predictability of the model was assessed by testing the model performance on external validation sets. A recently developed technique for selection of potential androgenically active chemicals was used to test the performance of the model in its applicability domain. Some of the selected chemicals were tested for AR transcriptional activation. The experimental results confirmed the theoretical predictions.

4D-QSAR analysis and pharmacophore modeling: electron conformational-genetic algorithm approach for penicillins

4D-QSAR studies were performed on a series of 87 penicillin analogues using the electron conformational-genetic algorithm (EC-GA) method. In this EC-based method, each conformation of the molecular system is described by a matrix (ECMC) with both electron structural parameters and interatomic distances as matrix elements. Multiple comparisons of these matrices within given tolerances for high active and low active penicillin compounds allow one to separate a smaller number of matrix elements (ECSA) which represent the pharmacophore groups. The effect of conformations was investigated building model 1 and 2 based on ensemble of conformers and single conformer, respectively. GA was used to select the most important descriptors and to predict the theoretical activity of the training (74 compounds) and test (13 compounds, commercial penicillins) sets. The model 1 for training and test sets obtained by optimum 12 parameters gave more satisfactory results (R(training)(2)=0.861, SE(training)=0.044, R(test)(2)=0.892, SE(test)=0.099, q(2)=0.702, q(ext1)(2)=0.777 and q(ext2)(2)=0.733) than model 2 (R(training)(2)=0.774, SE(training)=0.056, R(test)(2)=0.840, SE(test)=0.121, q(2)=0.514, q(ext1)(2)=0.641 and q(ext2)(2)=0.570). To estimate the individual influence of each of the molecular descriptors on biological activity, the E statistics technique was applied to the derived EC-GA model.

Bioactive conformational generation of small molecules: a comparative analysis between force-field and multiple empirical criteria based methods

Background

Conformational sampling for small molecules plays an essential role in drug discovery research pipeline. Based on multi-objective evolution algorithm (MOEA), we have developed a conformational generation method called Cyndi in the previous study. In this work, in addition to Tripos force field in the previous version, Cyndi was updated by incorporation of MMFF94 force field to assess the conformational energy more rationally. With two force fields against a larger dataset of 742 bioactive conformations of small ligands extracted from PDB, a comparative analysis was performed between pure force field based method (FFBM) and multiple empirical criteria based method (MECBM) hybrided with different force fields.

Results

Our analysis reveals that incorporating multiple empirical rules can significantly improve the accuracy of conformational generation. MECBM, which takes both empirical and force field criteria as the objective functions, can reproduce about 54% (within 1Å RMSD) of the bioactive conformations in the 742-molecule testset, much higher than that of pure force field method (FFBM, about 37%). On the other hand, MECBM achieved a more complete and efficient sampling of the conformational space because the average size of unique conformations ensemble per molecule is about 6 times larger than that of FFBM, while the time scale for conformational generation is nearly the same as FFBM. Furthermore, as a complementary comparison study between the methods with and without empirical biases, we also tested the performance of the three conformational generation methods in MacroModel in combination with different force fields. Compared with the methods in MacroModel, MECBM is more competitive in retrieving the bioactive conformations in light of accuracy but has much lower computational cost.

Conclusions

By incorporating different energy terms with several empirical criteria, the MECBM method can produce more reasonable conformational ensemble with high accuracy but approximately the same computational cost in comparison with FFBM method. Our analysis also reveals that the performance of conformational generation is irrelevant to the types of force field adopted in characterization of conformational accessibility. Moreover, post energy minimization is not necessary and may even undermine the diversity of conformational ensemble. All the results guide us to explore more empirical criteria like geometric restraints during the conformational process, which may improve the performance of conformational generation in combination with energetic accessibility, regardless of force field types adopted.

Dynamic clustering threshold reduces conformer ensemble size while maintaining a biologically relevant ensemble

Representing the 3D structures of ligands in virtual screenings via multi-conformer ensembles can be computationally intensive, especially for compounds with a large number of rotatable bonds. Thus, reducing the size of multi-conformer databases and the number of query conformers, while simultaneously reproducing the bioactive conformer with good accuracy, is of crucial interest. While clustering and RMSD filtering methods are employed in existing conformer generators, the novelty of this work is the inclusion of a clustering scheme (NMRCLUST) that does not require a user-defined cut-off value. This algorithm simultaneously optimizes the number and the average spread of the clusters. Here we describe and test four inter-dependent approaches for selecting computer-generated conformers, namely: OMEGA, NMRCLUST, RMS filtering and averaged-RMS filtering. The bioactive conformations of 65 selected ligands were extracted from the corresponding protein:ligand complexes from the Protein Data Bank, including eight ligands that adopted dissimilar bound conformations within different receptors. We show that NMRCLUST can be employed to further filter OMEGA-generated conformers while maintaining biological relevance of the ensemble. It was observed that NMRCLUST (containing on average 10 times fewer conformers per compound) performed nearly as well as OMEGA, and both outperformed RMS filtering and averaged-RMS filtering in terms of identifying the bioactive conformations with excellent and good matches (0.5 < RMSD < 1.0 A). Furthermore, we propose thresholds for OMEGA root-mean square filtering depending on the number of rotors in a compound: 0.8, 1.0 and 1.4 for structures with low (1-4), medium (5-9) and high (10-15) numbers of rotatable bonds, respectively. The protocol employed is general and can be applied to reduce the number of conformers in multi-conformer compound collections and alleviate the complexity of downstream data processing in virtual screening experiments.

Mechanism-based common reactivity pattern (COREPA) modelling of aryl hydrocarbon receptor binding affinity

The aryl hydrocarbon receptor is a ligand-activated transcription factor responsive to both natural and synthetic environmental compounds, with the most potent agonist being 2,3,7,8-tetrachlotrodibenzo-p-dioxin. The aim of this work was to develop a categorical COmmon REactivity PAttern (COREPA)-based structure-activity relationship model for predicting aryl hydrocarbon receptor ligands within different binding ranges. The COREPA analysis suggested two different binding mechanisms called dioxin- and biphenyl-like, respectively. The dioxin-like model predicts a mechanism that requires a favourable interaction with a receptor nucleophilic site in the central part of the ligand and with electrophilic sites at both sides of the principal molecular axis, whereas the biphenyl-like model predicted a stacking-type interaction with the aryl hydrocarbon receptor allowing electron charge transfer from the receptor to the ligand. The current model was also adjusted to predict agonistic/antagonistic properties of chemicals. The mechanism of antagonistic properties was related to the possibility that these chemicals have a localized negative charge at the molecule's axis and ultimately bind with the receptor surface through the electron-donating properties of electron-rich groups. The categorization of chemicals as agonists/antagonists was found to correlate with their gene expression. The highest increase in gene expression was elicited by strong agonists, followed by weak agonists producing lower increases in gene expression, whereas all antagonists (and non-aryl hydrocarbon receptor binders) were found to have no effect on gene expression. However, this relationship was found to be quantitative for the chemicals populating the areas with extreme gene expression values only, leaving a wide fuzzy area where the quantitative relationship was unclear. The total concordance of the derived aryl hydrocarbon receptor binding categorical structure-activity relationship model was 82% whereas the Pearson's coefficient was 0.88.

Toward the development of "nano-QSARs": advances and challenges

The most significant achievements and challenges relating to an application of quantitative structure-activity relationship (QSAR) approach in the risk assessment of nanometer-sized materials are highlighted. Recent advances are discussed in the context of "classical" QSAR methodology. The possible ways for the structural characterization of compounds existing at the nanoscale (at least one dimension of 100 nm or less) are briefly reviewed. The applicability of the existing toxicological data for developing QSAR models is evaluated. Finally, the existing models are presented. The need to develop new interpretative descriptors for the nanosystems is also highlighted. It is suggested that, due to high variability in the molecular structures and different mechanisms of toxicity, individual classes of nanoparticles should be modeled separately.

Cyndi: a multi-objective evolution algorithm based method for bioactive molecular conformational generation

BACKGROUND

Conformation generation is a ubiquitous problem in molecule modelling. Many applications require sampling the broad molecular conformational space or perceiving the bioactive conformers to ensure success. Numerous in silico methods have been proposed in an attempt to resolve the problem, ranging from deterministic to non-deterministic and systemic to stochastic ones. In this work, we described an efficient conformation sampling method named Cyndi, which is based on multi-objective evolution algorithm.

RESULTS

The conformational perturbation is subjected to evolutionary operation on the genome encoded with dihedral torsions. Various objectives are designated to render the generated Pareto optimal conformers to be energy-favoured as well as evenly scattered across the conformational space. An optional objective concerning the degree of molecular extension is added to achieve geometrically extended or compact conformations which have been observed to impact the molecular bioactivity (J Comput -Aided Mol Des 2002, 16: 105-112). Testing the performance of Cyndi against a test set consisting of 329 small molecules reveals an average minimum RMSD of 0.864 A to corresponding bioactive conformations, indicating Cyndi is highly competitive against other conformation generation methods. Meanwhile, the high-speed performance (0.49 +/- 0.18 seconds per molecule) renders Cyndi to be a practical toolkit for conformational database preparation and facilitates subsequent pharmacophore mapping or rigid docking. The copy of precompiled executable of Cyndi and the test set molecules in mol2 format are accessible in Additional file 1.

CONCLUSION

On the basis of MOEA algorithm, we present a new, highly efficient conformation generation method, Cyndi, and report the results of validation and performance studies comparing with other four methods. The results reveal that Cyndi is capable of generating geometrically diverse conformers and outperforms other four multiple conformer generators in the case of reproducing the bioactive conformations against 329 structures. The speed advantage indicates Cyndi is a powerful alternative method for extensive conformational sampling and large-scale conformer database preparation.