Modeling the value of predictive affinity scoring in preclinical drug discovery.

Drug discovery is widely recognized to be a difficult and costly activity in large part due to the challenge of identifying chemical matter which simultaneously optimizes multiple properties, one of which is affinity for the primary biological target. Further, many of these properties are difficult to predict ahead of expensive and time-consuming compound synthesis and experimental testing. Here we highlight recent work to develop compound affinity prediction models, and extensively investigate the value such models may provide to preclinical drug discovery. We demonstrate that the ability of these models to improve the overall probability of success is crucially dependent on the shape of the error distribution, not just the root-mean-square error. In particular, while scoring more molecule ideas generally improves the probability of project success when the error distribution is Gaussian, fat-tail distributions such as a Cauchy distribution, can lead to a situation where scoring more ideas actually decreases the overall probability of success.

Docking covalent inhibitors: a parameter free approach to pose prediction and scoring.

Although many popular docking programs include a facility to account for covalent ligands, large-scale systematic docking validation studies of covalent inhibitors have been sparse. In this paper, we present the development and validation of a novel approach for docking and scoring covalent inhibitors, which consists of conventional noncovalent docking, heuristic formation of the covalent attachment point, and structural refinement of the protein-ligand complex. This approach combines the strengths of the docking program Glide and the protein structure modeling program Prime and does not require any parameter fitting for the study of additional covalent reaction types. We first test this method by predicting the native binding geometry of 38 covalently bound complexes. The average RMSD of the predicted poses is 1.52 Å, and 76% of test set inhibitors have an RMSD of less than 2.0 Å. In addition, the apparent affinity score constructed herein is tested on a virtual screening study and the characterization of the SAR properties of two different series of congeneric compounds with satisfactory success.

Bonded exciplex formation: electronic and stereoelectronic effects.

As recently proposed, the singlet-excited states of several cyanoaromatics react with pyridine via bonded-exciplex formation, a novel concept in photochemical charge transfer reactions. Presented here are electronic and steric effects on the quenching rate constants, which provide valuable support for the model. Additionally, excited-state quenching in poly(vinylpyridine) is strongly inhibited both relative to that in neat pyridine and also to conventional exciplex formation in polymers, consistent with a restrictive orientational requirement for the formation of bonded exciplexes. Examples of competing reactions to form both conventional and bonded exciplexes are presented, which illustrate the delicate balance between these two processes when their reaction energetics are similar. Experimental and computational evidence is provided for the formation of a bonded exciplex in the reaction of the singlet excited state of 2,6,9,10-tetracyanoanthracene (TCA) with an oxygen-substituted donor, dioxane, thus expanding the scope of bonded exciplexes.

Bonded exciplexes. A new concept in photochemical reactions.

Charge-transfer quenching of the singlet excited states of cyanoaromatic electron acceptors by pyridine is characterized by a driving force dependence that resembles those of conventional electron-transfer reactions, except that a plot of the log of the quenching rate constants versus the free energy of electron transfer is displaced toward the endothermic region by 0.5-0.8 eV. Specifically, the reactions with pyridine display rapid quenching when conventional electron transfer is highly endothermic. As an example, the rate constant for quenching of the excited dicyanoanthracene is 3.5 x 10(9) M(-1)s(-1), even though formation of a conventional radical ion pair, A*-D*+, is endothermic by approximately 0.6 eV. No long-lived radical ions or exciplex intermediates can be detected on the picosecond to microsecond time scale. Instead, the reactions are proposed to proceed via formation of a previously undescribed, short-lived charge-transfer intermediate we call a "bonded exciplex", A- -D+. The bonded exciplex can be formally thought of as resulting from bond formation between the unpaired electrons of the radical ions A*- and D*+. The covalent bonding interaction significantly lowers the energy of the charge-transfer state. As a result of this interaction, the energy decreases with decreasing separation distance, and near van der Waals contact, the A- -D+ bonded state mixes with the repulsive excited state of the acceptor, allowing efficient reaction to form A- -D+ even when formation of a radical ion pair A*-D*+ is thermodynamically forbidden. Evidence for the bonded exciplex intermediate comes from studies of steric and Coulombic effects on the quenching rate constants and from extensive DFT computations that clearly show a curve crossing between the ground state and the low-energy bonded exciplex state.

Mimicking photosynthesis in a computationally designed synthetic metalloprotein.

While advances in protein design have made possible the construction of protein architectures with nativelike properties and predictable structures and function, there are as of yet no examples of functional, protein-based, solar energy conversion systems. This communication describes the design and characterization of an artificial reaction center (RC) protein that closely resembles the function of the natural photosynthetic RC. The synthetic protein, designed by the protein design program CORE, participates in multiple reduction/oxidation cycles with exogenous acceptors/donors following photoexcitation. The designed metalloprotein, aRC, consists of a tetrahelical bundle functionalized with two bis-histidine bound metal cofactors: a Ru(bpy)2 moiety and a heme group. Two distinct bis-histidine binding sites were engineered for each of these metal centers. Photoexcitation of aRC results in rapid ET from the RuII complex to the heme group (kET >/= 5 x 1010 s-1) yielding a long-lived (70 ns) charge-separated state (CSS), RuIII/FeII. This long-lived CSS participates in subsequent ET reactions with exogenous donors and acceptors in multiple photocycles, thus mimicking the basic function of native photosynthetic RCs. This study illustrates the successful design and construction of a protein-based functional charge separation device using a combination of automated computational protein design and knowledge of the engineering principles of biological electron tunneling extracted from natural electron-transfer systems. To our knowledge, this represents the first example of a functional protein-based artificial reaction center.