Virtual screening for R-groups, including predicted pIC50 contributions, within large structural databases, using Topomer CoMFA

Structural Insights on Hyp-Gly-Containing Peptides as Antiplatelet Compounds through Topomer CoMFA and CoMSIA Analysis

Increasing evidence has shown collagen hydrolysate involves a variety of bioactivities. In our previous study, multiple antiplatelet peptides containing Hyp/Pro-Gly were identified in collagen hydrolysates from Salmo salar and silver carp skin and exhibited anti-thrombosis activity without bleeding risks in vivo. However, the relationship between structure and activity remains unknown. We performed 3D-QSAR studies on 23 Hyp/Pro-Gly-containing peptides in which 13 peptides were reported before. CoMFA, Topomer CoMFA and CoMSIA analyses were used to generate the QSAR models. Topomer CoMFA analysis showed a q² value of 0.710, an r² value of 0.826, an r²_pred value of 0.930, and the results showed that Hyp instead of Pro was more important for improving the antiplatelet activity. CoMSIA analysis showed a q² value of 0.461, an r² value of 0.999, and an r²_pred value of 0.999. Compared with the electrostatic field and hydrogen bond donor field, the steric field, hydrophobic field and hydrogen bond receptor field have great influence on the activity of antiplatelet peptides. The predicted peptide EOGE exhibited antiplatelet activity induced by ADP, and inhibited thrombus formation (300 μmol/kg bw) without bleeding risks. Combined results of these studies indicate that OG-containing peptides had a potential to be developed into an effective specific medical food in the prevention of thrombotic diseases.

Insights into the structural requirements of triazole derivatives as promising DPP IV inhibitors: computational investigations

Diabetes is one of the leading causes of death globally as per World Health Organization 2019. To cope up with side effects of current diabetes therapy, researchers have found several novel targets for the treatment of diabetes. Currently, dipeptidyl peptidase IV (DPP IV) has emerged as a target in modulating the diabetes physiology. In the present work, various 3D-Quantitative structure activity relationship (QSAR) techniques namely comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis, topomer CoMFA and molecular hologram QSAR are used to explore the structural requirements of triazole derivatives as DPP IV inhibitors. Different models generated by 3D QSAR studies had acceptable statistical values for further prediction of molecules. From the contour maps of QSAR results, important structural features are deduced. Substitutions on N₁ and N₂ of triazole ring with H-bond donor group enhances the biological activity. Aliphatic side chain, less bulky group, H-bond donor group and -COOH group on N₃ of triazole ring are vital for the DPP IV inhibition. Moreover, electron withdrawing side chain on the triazole ring improves the biological activity. Further, novel triazole derivatives were designed and docking results of these compounds proved the efficiency of the developed 3D QSAR model. In future, results of this study may provide promising DPP IV inhibitors for the treatment of diabetes. Communicated by Ramaswamy H. Sarma.

Topomer CoMFA and HQSAR Study on Benzimidazole Derivative as NS5B Polymerase Inhibitor

Background:

In recent years, the number of people infected with the hepatitis C virus (HCV) is increasing rapidly. This has become a major threat to global health, therefore, new anti- HCV drugs are urgently needed. HCV NS5B polymerase is an RNA-dependent RNA polymerase (RdRp), which plays an important role in virus replication, and can effectively prevent the replication of HCV sub-genomic RNA in daughter cells. It is considered a very promising HCV therapeutic target for the design of anti-HCV drugs.

Methods:

In order to explore the relationship between the structure of benzimidazole derivative and its inhibitory activity on NS5B polymerase, holographic quantitative structure-activity relationship (HQSAR) and Topomer comparative molecular field analysis (CoMFA) were used to establish benzimidazole QSAR model of derivative inhibitors.

Results:

The results show that for the Topomer CoMFA model, the cross-validation coefficient q2 value is 0.883, and the non-cross-validation coefficient r2 value is 0.975. The model is reasonable, reliable, and has a good predictive ability. For the HQSAR model, the cross-validated q2 value is 0.922, and the uncross-validated r2 value is 0.971, indicating that the model data fit well and has a high predictive ability. Through the analysis of the contour map and color code diagram, 40 new benzimidazole inhibitor molecules were designed, and all of them have higher activity than template molecules, and the new molecules have significant interaction sites with protein 3SKE.

Conclusion:

The 3D-QSAR model established by Topomer CoMFA and HQSAR has good prediction results and the statistical verification is valid. The newly designed molecules and docking results provide theoretical guidance for the synthesis of new NS5B polymerase inhibitors and for the identification of key residues that the inhibitors bind to NS5B, which helps to better understand their inhibitory mechanism. These findings are helpful for the development of new anti-HCV drugs.

QSAR analysis of 3-pyrimidin-4-yl-oxazolidin-2-one derivatives isocitrate dehydrogenase inhibitors using Topomer CoMFA and HQSAR methods

A series of mIDH1 inhibitors derived from 3-pyrimidine-4-oxazolidin-2-ketone derivatives were studied by QSAR model to explore the key factors that inhibit mIDH1 activity. The generated model was cross-verified and non-cross-verified by Topomer CoMFA and HQSAR methods; the independent test set was verified by PLS method; the Topomer search technology was used for virtual screening and molecular design; and the Surflex-Dock method and ADMET technology were used for molecular docking, pharmacology and toxicity prediction of the designed drug molecules. The Topomer CoMFA and HQSAR cross-validation coefficients q² are 0.783 and 0.784, respectively, and the non-cross-validation coefficients r² are 0.978 and 0.934, respectively. Ten new drug molecules have been designed using Topomer search technology. The results of molecular docking and ADMET show that the newly designed drug molecules are effective. The docking situation, pharmacology and toxicity prediction results are good. The model can be used to predict the bioactivity of the same type of new compounds and their derivatives. The prediction results of molecular design, molecular docking and ADMET can provide some ideas for the design and development of novel mIDH1 inhibitor anticancer drugs, and provide certain theoretical basis of the experimental verification of new compounds in the future. Newly designed molecules after docking with corresponding proteins in the PDB library, it can explore the targets of drug molecules acting with large proteins and the related force, which is very helpful for the design of new drugs and the mechanism of drug action.

A computational approach for designing novel SARS-CoV-2 Mpro inhibitors: combined QSAR, molecular docking, and molecular dynamics simulation techniques

The promising compound T₂₁ for treating COVID-19 at the active site of SARS-CoV-2 M^pro.

Using Topomer Comparative Molecular Field Analysis to Elucidate Activity Differences of Aminomethylenethiophene Derivatives as Lysyl Oxidase Inhibitors: Implications for Rational Design of Antimetastatic Agents for Cancer Therapy

Topomer comparative molecular field analysis (topomer CoMFA) is applied to the quantitative structure-activity relationship (QSAR) study of aminomethylenethiophene (AMT) derivatives as lysyl oxidase (LOX) inhibitors. A total of thirty-six AMT derivatives were selected to build the QSAR model. The established topomer CoMFA model has the non-cross-validated correlation coefficient (r2) of 0.912 and the leave-one-out correlation coefficient (q2) of 0.540, which is statistically significant. The theoretically predicted anti-LOX potency agrees well with the experimentally observed inhibitory activity, proving the reasonable predictive ability of the QSAR model. The effect of molecular field information on the LOX inhibition of substituted aminomethylenethiophene was discussed in detail. The structural modification of the aminomethylenethiophene scaffold was carried out, and novel AMT derivatives with theoretically decent LOX inhibition were proposed. The topomer CoMFA modeling could provide a quantitative perspective into the structure-activity relationship of AMT derivatives and potentially speed up the rational design of LOX inhibitors as antimetastatic agents for cancer therapy.

Modeling MEK4 Kinase Inhibitors through Perturbed Electrostatic Potential Charges

MEK4, mitogen-activated protein kinase kinase 4, is overexpressed and induces metastasis in advanced prostate cancer lesions. However, the value of MEK4 as an oncology target has not been pharmacologically validated because selective chemical probes targeting MEK4 have not been developed. With advances in both computer and biological high-throughput screening, selective chemical entities can be discovered. Structure-based quantitative structure-activity relationship (QSAR) modeling often fails to generate accurate models due to poor alignment of training sets containing highly diverse compounds. Here we describe a highly predictive, nonalignment based robust QSAR model based on a data set of strikingly diverse MEK4 inhibitors. We computed the electrostatic potential (ESP) charges using a density functional theory (DFT) formalism of the donor and acceptor atoms of the ligands and hinge residues. Novel descriptors were then generated from the perturbation of the charge densities of the donor and acceptor atoms and were used to model a diverse set of 84 compounds, from which we built a robust predictive model.

Two- and three-dimensional QSAR studies on hURAT1 inhibitors with flexible linkers: topomer CoMFA and HQSAR

hURAT1 (human urate transporter 1) is a successful target for hyperuricemia. Recently, the modification work on hURAT1 inhibitors showed that the flexible linkers would benefit biological activity. The study aimed to investigate the contribution of the linkers and give modification strategies on this kind of structures based on QSAR models (HQSAR and topomer CoMFA). The most effective HQSAR and topomer CoMFA models were generated by applying the training set containing 63 compounds, with the cross-validated q² values of 0.869/0.818 and the non-cross-validated correlation coefficients r² of 0.951/0.978, respectively. The Y-randomization test was applied to ensure the robustness of the models. The external predictive correlation coefficient (r_pred²) grounded on the external test set (21 compounds) of two models was 0.910 and 0.907, respectively. In addition, the models were validated by Golbraikh-Tropsha and Roy methods, as well as other statistical metrics. The results showed that both models were reliable. Topomer CoMFA steric/electrostatic contours and HQSAR atomic contribution maps illustrated the structural features which governed their inhibitory potency. The dependable results could provide important insights to guide the designing of more potential hURAT1 inhibitors.

Antitumor evaluation and 3D-QSAR studies of a new series of the spiropyrroloquinoline isoindolinone/aza-isoindolinone derivatives by comparative molecular field analysis (CoMFA)

In current study, antitumor activity of two series of the newly synthesized spiropyrroloquinoline isoindolinone and spiropyrroloquinoline aza-isoindolinone scaffolds was evaluated against three human breast normal and cancer cell lines (MCF-10A, MCF-7 and SK-BR-3) and compared with cytotoxicity values of doxorubicin and colchicine as the standard drugs. It was found that several compounds were endowed with cytotoxicity in the low micromolar range. Among these two series, compounds 6i, 6j, 6k and 7l, 7m, 7n, 7o containing 3-ethyl-1H-indole moiety were found to be highly effective against both cancer cell lines ranging from [Formula: see text] to [Formula: see text] in comparison with the corresponding analogs. Compared with human cancer cells, the most potent compounds did not show high cytotoxicity against human breast normal MCF-10A cells. Generally, most of the evaluated compounds 6a-l and 7a-o series showed more antitumor activity against SK-BR-3 than MCF-7 cells. Moreover, comparative molecular field analysis (CoMFA) as a popular tools of three-dimensional quantitative structure-activity relationship (3D-QSAR) studies was carried out on 27 spiropyrroloquinolineisoindolinone and spiropyrroloquinolineaza-isoindolinone derivatives with antitumor activity against on SK-BR-3 cells. The obtained CoMFA models showed statistically excellent performance, which also possessed good predictive ability for an external test set. The results confirm the important effect of molecular steric and electrostatic interactions of these compounds on in vitro cytotoxicity against SK-BR-3.

Combined HQSAR, topomer CoMFA, homology modeling and docking studies on triazole derivatives as SGLT2 inhibitors

Aim: Sodium–glucose cotransporter 2 (SGLT2) is a promising target for diabetes therapy. We aimed to develop computational approaches to identify structural features for more potential SGLT2 inhibitors. Materials & methods: In this work, 46 triazole derivatives as SGLT2 inhibitors were studied using a combination of several approaches, including hologram quantitative structure–activity relationships (HQSAR), topomer comparative molecular field analysis (CoMFA), homology modeling, and molecular docking. HQSAR and topomer CoMFA were used to construct models. Molecular docking was conducted to investigate the interaction of triazole derivatives and homology modeling of SGLT2, as well as to validate the results of the HQSAR and topomer CoMFA models. Results: The most effective HQSAR and topomer CoMFA models exhibited noncross-validated correlation coefficients of 0.928 and 0.891 for the training set, respectively. External predictions were made successfully on a test set and then compared with previously reported models. The graphical results of HQSAR and topomer CoMFA were proven to be consistent with the binding mode of the inhibitors and SGLT2 from molecular docking. Conclusion: The models and docking provided important insights into the design of potent inhibitors for SGLT2.

TRPV1 antagonism by piperazinyl-aryl compounds: A Topomer-CoMFA study and its use in virtual screening for identification of novel antagonists

Transient Receptor Potential Vanilloid, member 1 (TRPV1), is a non-selective cation channel belonging to the transient receptor potential (TRP) family of ion channels. It occurs in the peripheral and central nervous system, activated by a variety of exogenous and endogenous stimuli, thus playing a key role in transmission of pain. This has been a target for chronic pain since more than a decade and a number of antagonists that progressed into clinical trials have failed due to the unexpected side effect of core body temperature rise, thus halting progress in this field. Of late, there has been an upsurge in research on this target, with the rat TRPV1 structure being determined, many new antagonists discovered that are temperature-neutral and many new therapeutic avenues being discovered for TRPV1, including diseases of respiratory and digestive systems, skin and bladder. Towards identifying diverse compounds to decipher the role of this target in various indications, here we report a 3D-QSAR model built using the new topomer-CoMFA methodology on a series of piperazinyl-aryl TRPV1 antagonists and the use of this model, along with a pharmacophore model and the shape of one of the potent compounds of this series, to virtually screen a subset of the ZINC database to find novel and diverse hits. These can serve as starting points to develop modality-selective antagonists for chronic pain and to elucidate the critical role of TRPV1 in the various new therapeutic areas.

Molecular modeling studies of [6,6,5] Tricyclic Fused Oxazolidinones as FXa inhibitors using 3D-QSAR, Topomer CoMFA, molecular docking and molecular dynamics simulations

Coagulation factor Xa (Factor Xa, FXa) is a particularly promising target for novel anticoagulant therapy. The first oral factor Xa inhibitor has been approved in the EU and Canada in 2008. In this work, 38 [6,6,5] Tricyclic Fused Oxazolidinones were studied using a combination of molecular modeling techniques including three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, molecular dynamics and Topomer CoMFA (comparative molecular field analysis) were used to build 3D-QSAR models. The results show that the best CoMFA model has q(2)=0.511 and r(2)=0.984, the best CoMSIA (comparative molecular similarity indices analysis) model has q(2)=0.700 and r(2)=0.993 and the Topomer CoMFA analysis has q(2)=0.377 and r(2)=0.886. The results indicated the steric, hydrophobic, H-acceptor and electrostatic fields play key roles in models. Molecular docking and molecular dynamics explored the binding relationship of the ligand and the receptor protein.

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.

Design Some New Type-I c-met Inhibitors Based on Molecular Docking and Topomer CoMFA Research

In this paper, a specific design strategy targeting c-met kinase was reported based on docking modeling and topomer comparative molecular field analysis (Topomer CoMFA). A novel U-shape conformation which is distinct from the literature was demonstrated by molecular docking among 68 U-shape c-met inhibitors. According to the docking results, two Topomer CoMFA models with high predictive ability were established based on the two fragment rule. The results from both docking and topomer CoMFA showed that the π-π stacking interaction with Tyr1230 and the hydrogen bond with hinge region play an important role in inhibitory activity. Furthermore, the flexible linker and the adjacent solvent group would be favorable to stabilize the conformation and to enhance the two interactions mentioned above. Based on our patent, 14 new compounds were designed by our design strategy. The binding mode exhibited as expected and their activities were predicted by topomer CoMFA model. The preliminary biological tests showed most of them have potent activity to c-met kinase. Our study would provide guidelines to design some new U-shaped c-met inhibitors with new scaffolds and optimize the current molecules.

Template CoMFA: the 3D-QSAR Grail?

Template CoMFA, a novel alignment methodology for training or test set structures in 3D-QSAR, is introduced. Its two most significant advantages are its complete automation and its ability to derive a single combined model from multiple structural series affecting a biological target. Its only two inputs are one or more "template" structures having 3D coordinates that share some Cartesian space, as may result from X-ray crystallography or pharmacophoric hypothesis, and one or more connectivity-only SAR tables associated with a common target. Template CoMFA also overcomes the major disadvantages of both existing 3D-QSAR alignment methodologies, specifically the tedium and subjectivity of familiar ad hoc approaches, and the awkwardness, occasional physicochemical heresies, and structural scope limitations of the purely topomer approach. The template CoMFA algorithms are described, and two of its application classes are presented. The first class, general models of binding to factor Xa and P38 map kinase, uses crystallographic structures as templates, with the encouraging result that the statistical qualities of each of these two combined models are equivalent to those of their constituent individual series models. The second, 15 data sets originally collected for validation of topomer CoMFA, with arbitrary structures as templates, confirms that the modeling power of template CoMFA resembles that of its predecessors.

X-ray crystallographic structures as a source of ligand alignment in 3D-QSAR

Three-dimensional quantitative structure-activity relationships (3D-QSAR) analyses are methods correlating a pharmacological property with a mathematical representation of a molecular property distribution around three-dimensional molecular models for a set of congeners. 3D-QSAR methods are known to be highly sensitive to ligand conformation and alignment method. The current study collects 32 unique positions of congeneric ligands co-crystallized with the binding domain of AMPA receptors and aligns them using protein coordinates. Thus, it allows for a unique opportunity to consider a ligands' orientation aligned by their mode of binding in a native molecular target. Comparative molecular field analysis (CoMFA) models were generated for this alignment and compared with the results of analogous modeling using standard structure-based alignment or obtained in docking simulations of the ligands' molecules. In comparison with classically derived models, the model based on X-ray crystallographic studies showed much better performance and statistical significance. Although the 3D-QSAR methods are mainly employed when crystallographic information is limited, the current study underscores the importance that the selection of inappropriate molecular conformations and alignment methods can lead to generation of erroneous models and false conclusions.

R-group template CoMFA combines benefits of "ad hoc" and topomer alignments using 3D-QSAR for lead optimization

Template CoMFA methodologies extend topomer CoMFA by allowing user-designated templates, for example the experimental receptor-bound conformation of a prototypical ligand, to help determine the alignment of training and test set structures for 3D-QSAR. The algorithms that generate its new structural modality, template-constrained topomers, are described. Template CoMFA's resolution of certain topomer CoMFA concerns, by providing user control of topological consistency and structural acceptability, is demonstrated for sixteen 3D-QSAR training sets, in particular the Selwood dataset.

The inevitable QSAR renaissance

QSAR approaches, including recent advances in 3D-QSAR, are advantageous during the lead optimization phase of drug discovery and complementary with bioinformatics and growing data accessibility. Hints for future QSAR practitioners are also offered.

Methods for combinatorial and parallel library design

Diversity has historically played a critical role in design of combinatorial libraries, screening sets and corporate collections for lead discovery. Large library design dominated the field in the 1990s with methods ranging anywhere from purely arbitrary through property based reagent selection to product based approaches. In recent years, however, there has been a downward trend in library size. This was due to increased information about the desirable targets gleaned from the genomics revolution and to the ever growing availability of target protein structures from crystallography and homology modeling. Creation of libraries directed toward families of receptors such as GPCRs, kinases, nuclear hormone receptors, proteases, etc., replaced the generation of libraries based primarily on diversity while single target focused library design has remained an important objective. Concurrently, computing grids and cpu clusters have facilitated the development of structure based tools that screen hundreds of thousands of molecules. Smaller "smarter" combinatorial and focused parallel libraries replaced those early un-focused large libraries in the twenty-first century drug design paradigm. While diversity still plays a role in lead discovery, the focus of current library design methods has shifted to receptor based methods, scaffold hopping/bio-isostere searching, and a much needed emphasis on synthetic feasibility. Methods such as "privileged substructures based design" and pharmacophore based design still are important methods for parallel and small combinatorial library design. This chapter discusses some of the possible design methods and presents examples where they are available.

Computational medicinal chemistry in fragment-based drug discovery: what, how and when

The use of fragment-based drug discovery (FBDD) has increased in the last decade due to the encouraging results obtained to date. In this scenario, computational approaches, together with experimental information, play an important role to guide and speed up the process. By default, FBDD is generally considered as a constructive approach. However, such additive behavior is not always present, therefore, simple fragment maturation will not always deliver the expected results. In this review, computational approaches utilized in FBDD are reported together with real case studies, where applicability domains are exemplified, in order to analyze them, and then, maximize their performance and reliability. Thus, a proper use of these computational tools can minimize misleading conclusions, keeping the credit on FBDD strategy, as well as achieve higher impact in the drug-discovery process. FBDD goes one step beyond a simple constructive approach. A broad set of computational tools: docking, R group quantitative structure–activity relationship, fragmentation tools, fragments management tools, patents analysis and fragment-hopping, for example, can be utilized in FBDD, providing a clear positive impact if they are utilized in the proper scenario – what, how and when. An initial assessment of additive/non-additive behavior is a critical point to define the most convenient approach for fragments elaboration.

Rethinking 3D-QSAR

The average error of pIC50 prediction reported for 140 structures in make-and-test applications of topomer CoMFA by four discovery organizations is 0.5. This remarkable accuracy can be understood to result from a topomer pose’s goal of generating field differences only at lattice intersections adjacent to intended structural change.

Recent advances in fragment-based QSAR and multi-dimensional QSAR methods

This paper provides an overview of recently developed two dimensional (2D) fragment-based QSAR methods as well as other multi-dimensional approaches. In particular, we present recent fragment-based QSAR methods such as fragment-similarity-based QSAR (FS-QSAR), fragment-based QSAR (FB-QSAR), Hologram QSAR (HQSAR), and top priority fragment QSAR in addition to 3D- and nD-QSAR methods such as comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA), Topomer CoMFA, self-organizing molecular field analysis (SOMFA), comparative molecular moment analysis (COMMA), autocorrelation of molecular surfaces properties (AMSP), weighted holistic invariant molecular (WHIM) descriptor-based QSAR (WHIM), grid-independent descriptors (GRIND)-based QSAR, 4D-QSAR, 5D-QSAR and 6D-QSAR methods.

Modeling approaches for ligand-based 3D similarity

3D ligand-based similarity approaches are widely used in the early phases of drug discovery for tasks such as hit finding by virtual screening or compound design with quantitative structure–activity relationships. Here in we review widely used software for performing such tasks. Some techniques are based on relatively mature technology, shape-based similarity for instance. Typically, these methods remained in the realm of the expert user, the experienced modeler. However, advances in implementation and speed have improved usability and allow these methods to be applied to databases comprising millions of compounds. There are now many reports of such methods impacting drug-discovery projects. As such, the medicinal chemistry community has become the intended market for some of these new tools, yet they may consider the wide array and choice of approaches somewhat disconcerting. Each method has subtle differences and is better suited to certain tasks than others. In this article we review some of the widely used computational methods via application, provide straightforward background on the underlying theory and provide examples for the interested reader to pursue in more detail. In the new era of preclinical drug discovery there will be ever more pressure to move faster and more efficiently, and computational approaches based on 3D ligand similarity will play an increasing role in in this process.