Considerations of Protein Subpockets in Fragment-Based Drug Design

CRD: A de novo design algorithm for the prediction of cognate protein receptors for small molecule ligands

While predicting a ligand that binds to a protein is feasible with current methods, the opposite, i.e., the prediction of a receptor for a ligand remains challenging. We present an approach for predicting receptors of a given ligand that uses de novo design and structural bioinformatics. We have developed the algorithm CRD, comprising multiple modules combining fragment-based sub-site finding, a machine learning function to estimate the size of the site, a genetic algorithm that encodes knowledge on protein structures and a physics-based fitness scoring scheme. CRD includes a pseudo-receptor design component followed by a mapping component to identify proteins that might contain these sites. CRD recovers the sites and receptors of several natural ligands. It designs similar sites for similar ligands, yet to some extent can distinguish between closely related ligands. CRD correctly predicts receptor classes for several drugs and might become a valuable tool for drug discovery.

CAVIAR: a method for automatic cavity detection, description and decomposition into subcavities

The accurate description of protein binding sites is essential to the determination of similarity and the application of machine learning methods to relate the binding sites to observed functions. This work describes CAVIAR, a new open source tool for generating descriptors for binding sites, using protein structures in PDB and mmCIF format as well as trajectory frames from molecular dynamics simulations as input. The applicability of CAVIAR descriptors is showcased by computing machine learning predictions of binding site ligandability. The method can also automatically assign subcavities, even in the absence of a bound ligand. The defined subpockets mimic the empirical definitions used in medicinal chemistry projects. It is shown that the experimental binding affinity scales relatively well with the number of subcavities filled by the ligand, with compounds binding to more than three subcavities having nanomolar or better affinities to the target. The CAVIAR descriptors and methods can be used in any machine learning-based investigations of problems involving binding sites, from protein engineering to hit identification. The full software code is available on GitHub and a conda package is hosted on Anaconda cloud.

Discovery of Inhibitors for Proliferating Cell Nuclear Antigen Using a Computational-Based Linked-Multiple-Fragment Screen

Proliferating cell nuclear antigen (PCNA) is a central factor in DNA replication and repair pathways that plays an essential role in genome stability. The functional roles of PCNA are mediated through an extensive list of protein-protein interactions, each of which transmits specific information in protein assemblies. The flexibility at the PCNA-protein interaction interfaces offers opportunities for the discovery of functionally selective inhibitors of DNA repair pathways. Current fragment-based drug design methodologies can be limited by the flexibility of protein interfaces. These factors motivated an approach to defining compounds that could leverage previously identified subpockets on PCNA that are suitable for fragment-binding sites. Methodologies for screening multiple connected fragment-binding events in distinct subpockets are deployed to improve the selection of fragment combinations. A flexible backbone based on N-alkyl-glycine amides offers a scaffold to combinatorically link multiple fragments for in silico screening libraries that explore the diversity of subpockets at protein interfaces. This approach was applied to discover new potential inhibitors of DNA replication and repair that target PCNA in a multiprotein recognition site. The screens of the libraries were designed to computationally filter ligands based upon the fragments and positions to <1%, which were synthesized and tested for direct binding to PCNA. Molecular dynamics simulations also revealed distinct features of these novel molecules that block key PCNA-protein interactions. Furthermore, a Bayesian classifier predicted 15 of the 16 new inhibitors to be modulators of protein-protein interactions, demonstrating the method's utility as an effective screening filter. The cellular activities of example ligands with similar affinity for PCNA demonstrate unique properties for novel selective synergy with therapeutic DNA-damaging agents in drug-resistant contexts.

Targeting the subpocket in xanthine oxidase: Design, synthesis, and biological evaluation of 2-[4-alkoxy-3-(1H-tetrazol-1-yl) phenyl]-6-oxo-1,6-dihydropyrimidine-5-carboxylic acid derivatives

Xanthine oxidase is an important target for the treatment of hyperuricemia, gout and other related diseases. Analysis of the high-resolution structure of xanthine oxidase with febuxostat identified the existence of a subpocket formed by the residues Leu648, Asn768, Lys771, Leu1014 and Pro1076. In this study, we designed and synthesized a series of 2-[4-alkoxy-3-(1H-tetrazol-1-yl) phenyl]-6-oxo-1,6-dihydropyrimidine-5-carboxylic acid derivatives (8a-8z) with a tetrazole group targeting this subpocket of the xanthine oxidase active site, and they were further evaluated for their inhibitory potency against xanthine oxidase in vitro. The results showed that all the tested compounds (8a-8z) exhibited an apparent xanthine oxidase inhibitory potency, with IC₅₀ values ranging from 0.0288 μM to 0.629 μM. Among them, compound 8u emerged as the most potent xanthine oxidase inhibitor, with an IC₅₀ value of 0.0288 μM, which was comparable to febuxostat (IC₅₀ = 0.0236 μM). The structure-activity relationship results revealed that the hydrophobic group at the 4'-position was indispensable for the inhibitory potency in vitro against xanthine oxidase. A Lineweaver-Burk plot revealed that the representative compound 8u acted as a mixed-type inhibitor for xanthine oxidase. Furthermore, molecular modeling studies were performed to gain insights into the binding mode of 8u with xanthine oxidase and suggested that the tetrazole group of the phenyl unit was accommodated in the subpocket, as expected. Moreover, a potassium oxonate-induced hyperuricemia model in rats was chosen to further confirm the hypouricemic effect of compound 8u, and the result demonstrated that compound 8u could effectively reduce serum uric acid levels at an oral dose of 5 mg/kg. In addition, acute oral toxicity study in mice indicated that compound 8u was nontoxic and tolerated at a dose up to 2000 mg/kg. Thus, compound 8u could be a potential and efficacious agent in treatment of hyperuricemia with low toxicity.

The use of matched molecular series networks for cross target structure activity relationship translation and potency prediction

Matched molecular series (MMS) analysis is an extension of matched molecular pair (MMP) analysis where all of the MMPs belong to the same chemical series. An MMS within a biological assay is able to capture specific structure activity relationships resulting from chemical substitution at a single location in the molecule. Under this convention, an MMS has the ability to capture one specific interaction vector between the compounds in a series and their therapeutic target. MMS analysis has the potential to translate the SAR from one series to another even across different protein targets or assays. A significant limitation of this approach is the lack of chemical series with a sufficient number of overlapping fragments to establish a statistically strong SAR in most databases. This results in either an inability to perform MMS analysis altogether or a potentially high proportion of spurious matches from chance correlations when the MMS compound count is low. This paper presents the novel concept of an MMS Network, which captures the SAR relationships between a set of related MMSs and significantly enhances the performance of MMS analysis by reducing the number of spurious matches leading to the identification of unexpected and potentially transferable SAR across assays. The results of a full retrospective leave-one-out analysis and randomization simulation are provided, and examples of pharmaceutically relevant programs will be presented to demonstrate the potential of this method.

Fragment-based drug discovery as alternative strategy to the drug development for neglected diseases

Neglected diseases (NDs) affect large populations and almost whole continents, representing 12% of the global health burden. In contrast, the treatment available today is limited and sometimes ineffective. Under this scenery, the Fragment-Based Drug Discovery emerged as one of the most promising alternatives to the traditional methods of drug development. This method allows achieving new lead compounds with smaller size of fragment libraries. Even with the wide Fragment-Based Drug Discovery success resulting in new effective therapeutic agents against different diseases, until this moment few studies have been applied this approach for NDs area. In this article, we discuss the basic Fragment-Based Drug Discovery process, brief successful ideas of general applications and show a landscape of its use in NDs, encouraging the implementation of this strategy as an interesting way to optimize the development of new drugs to NDs.

Rule of five in 2015 and beyond: Target and ligand structural limitations, ligand chemistry structure and drug discovery project decisions

The rule of five (Ro5), based on physicochemical profiles of phase II drugs, is consistent with structural limitations in protein targets and the drug target ligands. Three of four parameters in Ro5 are fundamental to the structure of both target and drug binding sites. The chemical structure of the drug ligand depends on the ligand chemistry and design philosophy. Two extremes of chemical structure and design philosophy exist; ligands constructed in the medicinal chemistry synthesis laboratory without input from natural selection and natural product (NP) metabolites biosynthesized based on evolutionary selection. Exceptions to Ro5 are found mostly among NPs. Chemistry chameleon-like behavior of some NPs due to intra-molecular hydrogen bonding as exemplified by cyclosporine A is a strong contributor to NP Ro5 outliers. The fragment derived, drug Navitoclax is an example of the extensive expertise, resources, time and key decisions required for the rare discovery of a non-NP Ro5 outlier.

Protein Binding Pocket Dynamics

The dynamics of protein binding pockets are crucial for their interaction specificity. Structural flexibility allows proteins to adapt to their individual molecular binding partners and facilitates the binding process. This implies the necessity to consider protein internal motion in determining and predicting binding properties and in designing new binders. Although accounting for protein dynamics presents a challenge for computational approaches, it expands the structural and physicochemical space for compound design and thus offers the prospect of improved binding specificity and selectivity. A cavity on the surface or in the interior of a protein that possesses suitable properties for binding a ligand is usually referred to as a binding pocket. The set of amino acid residues around a binding pocket determines its physicochemical characteristics and, together with its shape and location in a protein, defines its functionality. Residues outside the binding site can also have a long-range effect on the properties of the binding pocket. Cavities with similar functionalities are often conserved across protein families. For example, enzyme active sites are usually concave surfaces that present amino acid residues in a suitable configuration for binding low molecular weight compounds. Macromolecular binding pockets, on the other hand, are located on the protein surface and are often shallower. The mobility of proteins allows the opening, closing, and adaptation of binding pockets to regulate binding processes and specific protein functionalities. For example, channels and tunnels can exist permanently or transiently to transport compounds to and from a binding site. The influence of protein flexibility on binding pockets can vary from small changes to an already existent pocket to the formation of a completely new pocket. Here, we review recent developments in computational methods to detect and define binding pockets and to study pocket dynamics. We introduce five different classes of protein pocket dynamics: (1) appearance/disappearance of a subpocket in an existing pocket; (2) appearance/disappearance of an adjacent pocket on the protein surface in the direct vicinity of an already existing pocket; (3) pocket breathing, which may be caused by side-chain fluctuations or backbone or interdomain vibrational motion; (4) opening/closing of a channel or tunnel, connecting a pocket inside the protein with solvent, including lid motion; and (5) the appearance/disappearance of an allosteric pocket at a site on a protein distinct from an already existing pocket with binding of a ligand to the allosteric binding site affecting the original pocket. We suggest that the class of pocket dynamics, as well as the type and extent of protein motion affecting the binding pocket, should be factors considered in choosing the most appropriate computational approach to study a given binding pocket. Furthermore, we examine the relationship between pocket dynamics classes and induced fit, conformational selection, and gating models of ligand binding on binding kinetics and thermodynamics. We discuss the implications of protein binding pocket dynamics for drug design and conclude with potential future directions for computational analysis of protein binding pocket dynamics.