Structure-based design and protein X-ray analysis of a protein kinase inhibitor

Heterocycles in Breast Cancer Treatment: The Use of Pyrazole Derivatives

Among the aromatic heterocycle rings, pyrazole -a five-membered ring with two adjacent nitrogen atoms in its structure has been postulated as a potent candidate in the pharmacological context. This moiety is an interesting therapeutic target covering a broad spectrum of biological activities due to its presence in many natural substances. Hence, the potential of the pyrazole derivatives as antitumor agents has been explored in many investigations, showing promising results in some cases. In this sense, breast cancer, which is already the leading cause of cancer mortality in women in some countries, has been the topic selected for this review, which covers a range of different research from the earliest studies published in 2003 to the most recent ones in 2021.

Predicting the Prevalence of Alternative Warfarin Tautomers in Solution

Warfarin, a widely used oral anticoagulant, is prescribed as a racemic mixture. Each enantiomer of neutral Warfarin can exist in 20 possible tautomeric states leading to complex pharmacokinetics and uncertainty as to the relevant species under different conditions. Here, the ability of alternative computational approaches to predict the preferred tautomeric form(s) of neutral Warfarin in different solvents is examined. It is shown that varying the method used to estimate the heat of formation in vacuum (direct or via homodesmic reactions), whether entropic corrections were included, and the method used to estimate the free enthalpy of solvation (i.e., PCM, COSMO, or SMD implicit models or explicit solvent) lead to large differences in the predicted rank and relative populations of the tautomers. In this case, only a combination of the enthalpy of formation using homodesmic reactions and explicit solvent to estimate the free enthalpy of solvation yielded results compatible with the available experimental data. The work also suggests that a small but significant subset of the possible Warfarin tautomers are likely to be physiologically relevant.

Contribution of indazolinone tautomers to kinase activity

The design and synthesis of indazolinone containing kinase inhibitors are reported. Regioisomers that showed profound potency variation in previously-reported isoindolinone and aminoindazole systems were surprisingly found to have similar potencies in the case of the indazolinone chemical series. An interpretation using differential hinge hydrogen bonding and tautomeric equilibrium of indazolinone ring system is supported by quantum mechanics calculations. The equipotent inhibition of a representative kinase (KDR) by regioisomeric indazolinones 4 and 5 is clear evidence that in case of the indazolinone hinge, both tautomers are equally favored, and should be considered in design of inhibitors.

X-ray crystallography: Assessment and validation of protein-small molecule complexes for drug discovery

INTRODUCTION: Crystallography is the key initial component for structure-based and fragment-based drug design and can often generate leads that can be developed into high potency drugs. Therefore, huge sums of money are committed based on the outcome of crystallography experiments and their interpretation. AREAS COVERED: This review discusses how to evaluate the correctness of an X-ray structure, focusing on the validation of small molecule-protein complexes. Various types of inaccuracies found within the PDB are identified and the ramifications of these errors are discussed. The reader will gain an understanding of the key parameters that need to be inspected before a structure can be used in drug discovery efforts, as well as an appreciation of the difficulties of correctly interpreting electron density for small molecules. The reader will also be introduced to methods for validating small molecules within the context of a macromolecular structure. EXPERT OPINION: One of the reasons that ligand identification and positioning, within a macromolecular crystal structure, is so difficult is that the quality of small molecules widely varies in the PDB. For this reason, the PDB can not always be considered a reliable repository of structural information pertaining to small molecules, and this makes the derivation of general principles that govern small molecule-protein interactions more difficult.

Inhibitors of SARS-3CLpro: virtual screening, biological evaluation, and molecular dynamics simulation studies

SARS-CoV from the coronaviridae family has been identified as the etiological agent of Severe Acute Respiratory Syndrome (SARS), a highly contagious upper respiratory disease that reached epidemic status in 2002. SARS-3CL(pro), a cysteine protease indispensible to the viral life cycle, has been identified as one of the key therapeutic targets against SARS. A combined ligand and structure-based virtual screening was carried out against the Asinex Platinum collection. Multiple low micromolar inhibitors of the enzyme were identified through this search, one of which also showed activity against SARS-CoV in a whole cell CPE assay. Furthermore, multinanosecond explicit solvent simulations were carried out using the docking poses of the identified hits to study the overall stability of the binding site interactions as well as identify important changes in the interaction profile that were not apparent from the docking study. Cumulative analysis of the evaluated compounds and the simulation studies led to the identification of certain protein-ligand interaction patterns which would be useful in further structure based design efforts.

Challenges in the determination of the binding modes of non-standard ligands in X-ray crystal complexes

Despite its central role in structure based drug design the determination of the binding mode (position, orientation and conformation in addition to protonation and tautomeric states) of small heteromolecular ligands in protein:ligand complexes based on medium resolution X-ray diffraction data is highly challenging. In this perspective we demonstrate how a combination of molecular dynamics simulations and free energy (FE) calculations can be used to correct and identify thermodynamically stable binding modes of ligands in X-ray crystal complexes. The consequences of inappropriate ligand structure, force field and the absence of electrostatics during X-ray refinement are highlighted. The implications of such uncertainties and errors for the validation of virtual screening and fragment-based drug design based on high throughput X-ray crystallography are discussed with possible solutions and guidelines.

Tautomer preference in PDB complexes and its impact on structure-based drug discovery

Tautomer enrichment is a key step of ligand preparation prior to virtual screening. In this paper, we have investigated how tautomer preference in various media (water, gas phase, and crystal) compares to tautomer preference at the active site of the protein by analyzing the different possible H-bonding contacts for a set of 13 tautomeric structures. In addition, we have explored the impact of four different protocols for the enumeration of tautomers in virtual screening by using Flap, Glide, and Gold as docking tools on seven targets of the DUD data set. Excluding targets in which the binding does not involve tautomeric atoms (HSP90, p38, and VEGFR2), we found that the average receiver operating characteristic curve enrichment at 10% was 0.25 (Gold), 0.24 (Glide), and 0.50 (Flap) by considering only tautomers predicted to be unstable in water versus 0.41 (Gold), 0.56 (Glide), 0.51 (Flap) by limiting the enumeration process only to the predicted most stable tautomer. The inclusion of all tautomers (stable and unstable) yielded slightly poorer results than considering only the most stable form in water.

Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution

Generating the appropriate protonation states of drug-like molecules in solution is important for success in both ligand- and structure-based virtual screening. Screening collections of millions of compounds requires a method for determining tautomers and their energies that is sufficiently rapid, accurate, and comprehensive. To maximise enrichment, the lowest energy tautomers must be determined from heterogeneous input, without over-enumerating unfavourable states. While computationally expensive, the density functional theory (DFT) method M06-2X/aug-cc-pVTZ(-f) [PB-SCRF] provides accurate energies for enumerated model tautomeric systems. The empirical Hammett-Taft methodology can very rapidly extrapolate substituent effects from model systems to drug-like molecules via the relationship between pK(T) and pK(a). Combining the two complementary approaches transforms the tautomer problem from a scientific challenge to one of engineering scale-up, and avoids issues that arise due to the very limited number of measured pK(T) values, especially for the complicated heterocycles often favoured by medicinal chemists for their novelty and versatility. Several hundreds of pre-calculated tautomer energies and substituent pK(a) effects are tabulated in databases for use in structural adjustment by the program Epik, which treats tautomers as a subset of the larger problem of the protonation states in aqueous ensembles and their energy penalties. Accuracy and coverage is continually improved and expanded by parameterizing new systems of interest using DFT and experimental data. Recommendations are made for how to best incorporate tautomers in molecular design and virtual screening workflows.

Let's not forget tautomers

A compound exhibits tautomerism if it can be represented by two structures that are related by an intramolecular movement of hydrogen from one atom to another. The different tautomers of a molecule usually have different molecular fingerprints, hydrophobicities and pKa's as well as different 3D shape and electrostatic properties; additionally, proteins frequently preferentially bind a tautomer that is present in low abundance in water. As a result, the proper treatment of molecules that can tautomerize, approximately 25% of a database, is a challenge for every aspect of computer-aided molecular design. Library design that focuses on molecular similarity or diversity might inadvertently include similar molecules that happen to be encoded as different tautomers. Physical property measurements might not establish the properties of individual tautomers with the result that algorithms based on these measurements may be less accurate for molecules that can tautomerize-this problem influences the accuracy of filtering for library design and also traditional QSAR. Any 2D or 3D QSAR analysis must involve the decision of if or how to adjust the observed Ki or IC50 for the tautomerization equilibria. QSARs and recursive partitioning methods also involve the decision as to which tautomer(s) to use to calculate the molecular descriptors. Docking virtual screening must involve the decision as to which tautomers to include in the docking and how to account for tautomerization in the scoring. All of these decisions are more difficult because there is no extensive database of measured tautomeric ratios in both water and non-aqueous solvents and there is no consensus as to the best computational method to calculate tautomeric ratios in different environments.

Lessons for fragment library design: analysis of output from multiple screening campaigns

Over the past 8 years, we have developed, refined and applied a fragment based discovery approach to a range of protein targets. Here we report computational analyses of various aspects of our fragment library and the results obtained for fragment screening. We reinforce the finding of others that the experimentally observed hit rate for screening fragments can be related to a computationally defined druggability index for the target. In general, the physicochemical properties of the fragment hits display the same profile as the library, as is expected for a truly diverse library which probes the relevant chemical space. An analysis of the fragment hits against various protein classes has shown that the physicochemical properties of the fragments are complementary to the properties of the target binding site. The effectiveness of some fragments appears to be achieved by an appropriate mix of pharmacophore features and enhanced aromaticity, with hydrophobic interactions playing an important role. The analysis emphasizes that it is possible to identify small fragments that are specific for different binding sites. To conclude, we discuss how the results could inform further development and improvement of our fragment library.

Binding of small-molecule ligands to proteins: "what you see" is not always "what you get"

We review insights from computational studies of affinities of ligands binding to proteins. The power of structural biology is in translating knowledge of protein structures into insights about their forces, binding, and mechanisms. However, the complementary power of computer modeling is in showing "the rest of the story" (i.e., how motions and ensembles and alternative conformers and the entropies and forces that cannot be seen in single molecular structures also contribute to binding affinities). Upon binding to a protein, a ligand can bind in multiple orientations; the protein or ligand can be deformed by the binding event; waters, ions, or cofactors can have unexpected involvement; and conformational or solvation entropies can sometimes play large and otherwise unpredictable roles. Computer modeling is helping to elucidate these factors.

Molecular modeling and dynamics simulation of human cyclin-dependent kinase 3 complexed with inhibitors

The complex CDK3-cyclin is involved in the control of the progression of G0. While the mechanisms governing early and late G1 progression are well understood, very little is known about the G0-G1 transition. Human CDK3 is closely related to cyclin-dependent kinase 2 (CDK2). Since there is no crystallographic structure of human CDK3, this work describes for the first time a molecular model of human CDK3 complexed with several inhibitors. Comparison of the binary complexes with different inhibitors strongly indicates that those inhibitors should inhibit CDK3 as well as CDK2.

TrixX: structure-based molecule indexing for large-scale virtual screening in sublinear time

Structure-based virtual screening today is basically organized as a sequential process where the molecules of a screening library are evaluated for instance with respect to their fit with a biological target. In this paper, we present a novel structure-based screening paradigm avoiding sequential searching and therefore enabling sublinear runtime behavior. We implemented the novel paradigm in the virtual screening tool TrixX and successfully applied it in screening experiments on four targets from relevant therapeutic areas. With the screening paradigm implemented in TrixX, we propose some important extensions and modifications to traditional virtual screening approaches: Instead of processing all compounds in the screening library sequentially, TrixX first analyzes the geometric and physicochemical binding site characteristics and then draws compounds with matching features from a compound catalog. The catalog organizes the compounds by their physicochemical and geometric features making use of relational database technology with indexed tables in order to support efficient queries for compounds with specific features. A key element of the compound catalog is a highly selective geometric descriptor that carries information on the type of functional groups of the compound, their Euclidian distance, the preferred interaction direction of each functional group, and the location of steric bulk around the triangle. In a re-docking experiment with 200 protein-ligand complexes, we could show that TrixX is able to correctly predict the location of ligand functional groups in co-crystallized complexes. In a retrospective virtual screening experiment for four different targets, the enrichment factors of TrixX are comparable to the enrichment factors of FlexX and FlexX-Scan. With computing times clearly below one second per compound, TrixX counts among the fastest virtual screening tools currently available and is nearly two orders of magnitude faster than standard FlexX.

Hit-to-lead optimization of 1,4-dihydroindeno[1,2-c]pyrazoles as a novel class of KDR kinase inhibitors

A series of 1,4-dihydroindeno[1,2-c]pyrazoles was prepared and evaluated for their enzymatic inhibition of KDR kinase. Computer modeling studies revealed the importance of attaching a basic side chain in predicting the binding mode of those compounds. Further investigation of structure-activity relationships led to 19, a lead compound with an acceptable selectivity profile, activity in whole cells, and good oral efficacy in an estradiol-induced murine uterine edema model of VEGF activity.

Selectivity and potency of cyclin-dependent kinase inhibitors

Members of the cyclin-dependent kinase (CDK) family play key roles in various cellular processes. There are 11 members of the CDK family known till now. CDKs are activated by forming noncovalent complexes with cyclins such as A-, B-, C-, D- (D1, D2, and D3), and E-type cyclins. Each isozyme of this family is responsible for particular aspects (cell signaling, transcription, etc) of the cell cycle, and some of the CDK isozymes are specific to certain kinds of tissues. Aberrant expression and overexpression of these kinases are evidenced in many disease conditions. Inhibition of isozymes of CDKs specifically can yield beneficiary treatment modalities with minimum side effects. More than 80 3-dimensional structures of CDK2, CDK5, and CDK6 complexed with inhibitors have been published. This review provides an understanding of the structural aspects of CDK isozymes and binding modes of various known CDK inhibitors so that these kinases can be better targeted for drug discovery and design. The amino acid residues that constitute the cyclin binding region, the substrate binding region, and the area around the adenosine triphosphate (ATP) binding site have been compared for CDK isozymes. Those amino acids at the ATP binding site that could be used to improve the potency and subtype specificity have been described.

Scaffold hopping

The aim of scaffold hopping is to discover structurally novel compounds starting from known active compounds by modifying the central core structure of the molecule. Scaffold hopping is a central task of modern medicinal chemistry requiring a multitude of techniques, which are discussed in this article. Their application has led to several molecules with chemically completely different core structures, and yet binding to the same receptor. Computational approaches for scaffold hopping highlight the challenges of the field that are still unsolved.:

Design and Characterization of Libraries of Molecular Fragments for Use in NMR Screening against Protein Targets

We have designed four generations of a low molecular weight fragment library for use in NMR-based screening against protein targets. The library initially contained 723 fragments which were selected manually from the Available Chemicals Directory. A series of in silico filters and property calculations were developed to automate the selection process, allowing a larger database of 1.79 M available compounds to be searched for a further 357 compounds that were added to the library. A kinase binding pharmacophore was then derived to select 174 kinase-focused fragments. Finally, an additional 61 fragments were selected to increase the number of different pharmacophores represented within the library. All of the fragments added to the library passed quality checks to ensure they were suitable for the screening protocol, with appropriate solubility, purity, chemical stability, and unambiguous NMR spectrum. The successive generations of libraries have been characterized through analysis of structural properties (molecular weight, lipophilicity, polar surface area, number of rotatable bonds, and hydrogen-bonding potential) and by analyzing their pharmacophoric complexity. These calculations have been used to compare the fragment libraries with a drug-like reference set of compounds and a set of molecules that bind to protein active sites. In addition, an analysis of the overall results of screening the library against the ATP binding site of two protein targets (HSP90 and CDK2) reveals different patterns of fragment binding, demonstrating that the approach can find selective compounds that discriminate between related binding sites.

Imidazo[1,2- a ]pyridines. Part 2 : SAR and optimisation of a potent and selective class of cyclin-dependent kinase inhibitors

Exploration of SAR and optimisation of the imidazo[1,2-a]pyridine CDK inhibitors has lead to the discovery of novel, potent and selective inhibitors of the cyclin-dependent kinase CDK2. Understanding of SAR has identified positions of substitution, which allow modification of physical properties and offer the potential for in vivo optimisation.

Identification of a new chemical class of potent angiogenesis inhibitors based on conformational considerations and database searching

The vascular endothelial growth factor (VEGF) tyrosine kinase receptors KDR and Flt-1 are targets of current interest in anticancer drug research. PTK787/ZK222584 is a potent inhibitor of these enzymes in clinical evaluation as an antiangiogenic agent. The development of a hypothesis concerning the bioactive conformation of this compound has led to the discovery of a new class of potent inhibitors of KDR and Flt-1, the anthranilamides. This could be achieved with a limited experimental effort, which only involved the testing of one archive compound and the synthesis and testing of one appropriate analogue.

Pharmacological inhibitors of cyclin-dependent kinases

Cyclin-dependent kinases (CDKs) regulate the cell division cycle, apoptosis, transcription and differentiation in addition to functions in the nervous system. Deregulation of CDKs in various diseases has stimulated an intensive search for selective pharmacological inhibitors of these kinases. More than 50 inhibitors have been identified, among which >20 have been co-crystallized with CDK2. These inhibitors all target the ATP-binding pocket of the catalytic site of the kinase. The actual selectivity of most known CDK inhibitors, and thus the underlying mechanism of their cellular effects, is poorly known. Pharmacological inhibitors of CDKs are currently being evaluated for therapeutic use against cancer, alopecia, neurodegenerative disorders (e.g. Alzheimer's disease, amyotrophic lateral sclerosis and stroke), cardiovascular disorders (e.g. atherosclerosis and restenosis), glomerulonephritis, viral infections (e.g. HCMV, HIV and HSV) and parasitic protozoa (Plasmodium sp. and Leishmania sp.).

Progress toward the development of agents to modulate the cell cycle

With taxanes continuing to prove useful in the clinical treatment of cancer, the next generation of antimitotic agents has entered clinical trials. Other mechanisms awaiting proof-of-concept for the treatment of antiproliferative diseases include inhibition of cyclin-dependent kinases (Cdks). Flavopiridol and UCN-01 are continuing in clinical trials, and newer more selective Cdk inhibitors are now entering clinical evaluation.