The Role of Phosphodiesterase 12 (PDE12) as a Negative Regulator of the Innate Immune Response and the Discovery of Antiviral Inhibitors.

2',5'-Oligoadenylate synthetase (OAS) enzymes and RNase-L constitute a major effector arm of interferon (IFN)-mediated antiviral defense. OAS produces a unique oligonucleotide second messenger, 2',5'-oligoadenylate (2-5A), that binds and activates RNase-L. This pathway is down-regulated by virus- and host-encoded enzymes that degrade 2-5A. Phosphodiesterase 12 (PDE12) was the first cellular 2-5A- degrading enzyme to be purified and described at a molecular level. Inhibition of PDE12 may up-regulate the OAS/RNase-L pathway in response to viral infection resulting in increased resistance to a variety of viral pathogens. We generated a PDE12-null cell line, HeLaΔPDE12, using transcription activator-like effector nuclease-mediated gene inactivation. This cell line has increased 2-5A levels in response to IFN and poly(I-C), a double-stranded RNA mimic compared with the parental cell line. Moreover, HeLaΔPDE12 cells were resistant to viral pathogens, including encephalomyocarditis virus, human rhinovirus, and respiratory syncytial virus. Based on these results, we used DNA-encoded chemical library screening to identify starting points for inhibitor lead optimization. Compounds derived from this effort raise 2-5A levels and exhibit antiviral activity comparable with the effects observed with PDE12 gene inactivation. The crystal structure of PDE12 complexed with an inhibitor was solved providing insights into the structure-activity relationships of inhibitor potency and selectivity.

Profiling cells with DELs: Small molecule fingerprinting of cell surfaces.

DNA-encoded small molecule library technology has recently emerged as a new paradigm for identifying ligands against drug targets. To date, it has been used to identify ligands against targets that are soluble or overexpressed on cell surfaces. Here, we report applying cell-based selection methods to profile surfaces of mouse C2C12 myoblasts and myotube cells in an unbiased, target agnostic manner. A panel of on-DNA compounds were identified and confirmed for cell binding selectivity. We optimized the cell selection protocol and employed a novel data analysis method to identify cell selective ligands against a panel of human B and T lymphocytes. We discuss the generality of using this workflow for DNA encoded small molecule library selection and data analysis against different cell types, and the feasibility of applying this method to profile cell surfaces for biomarker and target identification.

Discovery of Proline-Based p300/CBP Inhibitors Using DNA-Encoded Library Technology in Combination with High-Throughput Screening.

E1A binding protein (p300) and CREB binding protein (CBP) are two highly homologous and multidomain histone acetyltransferases. These two proteins are involved in many cellular processes by acting as coactivators of a large number of transcription factors. Dysregulation of p300/CBP has been found in a variety of cancers and other diseases, and inhibition has been shown to decrease Myc expression. Herein, we report the identification of a series of highly potent, proline-based small-molecule p300/CBP histone acetyltransferase (HAT) inhibitors using DNA-encoded library technology in combination with high-throughput screening. The strategy of reducing ChromlogD and fluorination of metabolic soft spots was explored to improve the pharmacokinetic properties of potent p300 inhibitors. Fluorination of both cyclobutyl and proline rings of 22 led to not only reduced clearance but also improved cMyc cellular potency.

DNA-Encoded Library Technology-Based Discovery, Lead Optimization, and Prodrug Strategy toward Structurally Unique Indoleamine 2,3-Dioxygenase-1 (IDO1) Inhibitors.

We report the discovery of a novel indoleamine 2,3-dioxygenase-1 (IDO1) inhibitor class through the affinity selection of a previously unreported indole-based DNA-encoded library (DEL). The DEL exemplar, spiro-chromane 1, had moderate IDO1 potency but high in vivo clearance. Series optimization quickly afforded a potent, low in vivo clearance lead 11. Although amorphous 11 was highly bio-available, crystalline 11 was poorly soluble and suffered disappointingly low bio-availability because of solubility-limited absorption. A prodrug approach was deployed and proved effective in discovering the highly bio-available phosphonooxymethyl 31, which rapidly converted to 11 in vivo. Obtaining crystalline 31 proved problematic, however; thus salt screening was performed in an attempt to circumvent this obstacle and successfully delivered greatly soluble and bio-available crystalline tris-salt 32. IDO1 inhibitor 32 is characterized by a low calculated human dose, best-in-class potential, and an unusual inhibition mode by binding the IDO1 heme-free (apo) form.

Explicit ion, implicit water solvation for molecular dynamics of nucleic acids and highly charged molecules.

An explicit ion, implicit water solvent model for molecular dynamics was developed and tested with DNA and RNA simulations. The implicit water model uses the finite difference Poisson (FDP) model with the smooth permittivity method implemented in the OpenEye ZAP libraries. Explicit counter-ions, co-ions, and nucleic acid were treated with a Langevin dynamics molecular dynamics algorithm. Ion electrostatics is treated within the FDP model when close to the solute, and by the Coulombic model when far from the solute. The two zone model reduces computation time, but retains an accurate treatment of the ion atmosphere electrostatics near the solute. Ion compositions can be set to reproduce specific ionic strengths. The entire ion/water treatment is interfaced with the molecular dynamics package CHARMM. Using the CHARMM-ZAPI software combination, the implicit solvent model was tested on A and B form duplex DNA, and tetraloop RNA, producing stable simulations with structures remaining close to experiment. The model also reproduced the A to B duplex DNA transition. The effect of ionic strength, and the structure of the counterion atmosphere around B form duplex DNA were also examined.

Accessibility of oxygen with respect to the heme pocket in horseradish peroxidase.

Oxygen and other molecules of similar size take part in a variety of protein reactions. Therefore, it is critical to understand how these small molecules penetrate the protein matrix. The protein system studied in this case is horseradish peroxidase (HRP). We have converted the native HRP into a phosphorescent analog by replacing the heme prosthetic group by Pd-mesoporphyrin. Oxygen readily quenches the phosphorescence of Pd porphyrins, and this can be used to characterize oxygen diffusion through the protein matrix. Our measurements indicate that solvent viscosity and pH modulate the accessibility of the heme pocket relative to small molecules. The binding of the substrate benzohydroxamic acid (BHA) to the protein drastically impedes oxygen access to the heme pocket. These results indicate that, first, the penetration of small molecules through the protein matrix is a function of protein dynamics, and second, there are specific pathways for the diffusion of these molecules. The effect of substrate and pH on protein dynamics has been investigated with the use of molecular dynamics calculations. We demonstrate that the model of a "fluctuating entry point," as suggested by Zwanzig (J Chem Phys 1992;97:3587-3589), properly describes the diffusion of oxygen through the protein matrix.

Solvent dependent and independent motions of CO-horseradish peroxidase examined by infrared spectroscopy and molecular dynamics calculations.

The role of the solvent matrix in affecting CO bound to ferrous horseradish peroxidase was examined by comparing band-widths of nu(CO) for the protein in aqueous solutions and in trehalose/sucrose glasses. We have previously observed that the optical absorption band and the CO stretching mode respond to the glass transition of glycerol/water in ways that depend upon the presence of substrate (Biochemistry 40 (2001) 3483). It is now demonstrated that the CO group band-width for the protein with bound inhibitor benzhydroxamic acid is relatively insensitive to temperature or the glass transition of the solvent. In contrast, in the absence of inhibitor, the band-width varies with the temperature that the glass is formed. The results show that solvent dependent and independent motions can be distinguished, and that the presence of substrate changes the protein such that the Fe[bond]CO site is occluded from the solvent conditions. Molecular dynamic calculations, based upon X-ray structures, showed that the presence of benzhydroxamic acid decreases the distance between His42 and Arg38 and this leads for closer distances to the O of the CO from these residues. These results are invoked to account for the observed line width changes of the CO band.

Heat capacity in proteins.

Heat capacity (Cp) is one of several major thermodynamic quantities commonly measured in proteins. With more than half a dozen definitions, it is the hardest of these quantities to understand in physical terms, but the richest in insight. There are many ramifications of observed Cp changes: The sign distinguishes apolar from polar solvation. It imparts a temperature (T) dependence to entropy and enthalpy that may change their signs and which of them dominate. Protein unfolding usually has a positive deltaCp, producing a maximum in stability and sometimes cold denaturation. There are two heat capacity contributions, from hydration and protein-protein interactions; which dominates in folding and binding is an open question. Theoretical work to date has dealt mostly with the hydration term and can account, at least semiquantitatively, for the major Cp-related features: the positive and negative Cp of hydration for apolar and polar groups, respectively; the convergence of apolar group hydration entropy at T approximately 112 degrees C; the decrease in apolar hydration Cp with increasing T; and the T-maximum in protein stability and cold denaturation.

Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth-permittivity finite difference Poisson-Boltzmann method.

A fast stable finite difference Poisson-Boltzmann (FDPB) model for implicit solvation in molecular dynamics simulations was developed using the smooth permittivity FDPB method implemented in the OpenEye ZAP libraries. This was interfaced with two widely used molecular dynamics packages, AMBER and CHARMM. Using the CHARMM-ZAP software combination, the implicit solvent model was tested on eight proteins differing in size, structure, and cofactors: calmodulin, horseradish peroxidase (with and without substrate analogue bound), lipid carrier protein, flavodoxin, ubiquitin, cytochrome c, and a de novo designed 3-helix bundle. The stability and accuracy of the implicit solvent simulations was assessed by examining root-mean-squared deviations from crystal structure. This measure was compared with that of a standard explicit water solvent model. In addition we compared experimental and calculated NMR order parameters to obtain a residue level assessment of the accuracy of MD-ZAP for simulating dynamic quantities. Overall, the agreement of the implicit solvent model with experiment was as good as that of explicit water simulations. The implicit solvent method was up to eight times faster than the explicit water simulations, and approximately four times slower than a vacuum simulation (i.e., with no solvent treatment).

Dynamics and entropy of a calmodulin-peptide complex studied by NMR and molecular dynamics.

All-atom, explicit water molecular dynamics simulations of calcium-loaded calmodulin complexed with a peptide corresponding to the smooth muscle myosin light chain kinase target were carried out at 295 and 346 K. Amide and side chain methyl angular generalized order parameters were calculated and analyzed in the context of the protein's structure and dynamics. The agreement between amide order parameters measured by NMR and those from the simulations was found to be good, especially at the higher temperature, indicating both better convergence for the latter and excellent transferrability of the CHARMM parameters to the higher temperature. Subtle dynamical features such as helix fraying were reproduced. A large range of order parameters for the nine calmodulin methionines was observed in the NMR, and reproduced quite well in the simulations. The major determinant of the methionine order parameter was found to be the proximity to side chains of aromatic residues. An upper bound estimate of the difference in backbone entropy between loop and helical regions was extracted from the order parameters using a model of motion in an effective potential. Although loop regions are more flexible than helical regions, it was found that the entropy loss per residue upon folding was only approximately 20% less for loops than for helices. Pairwise correlated motions, which could significantly lower entropy estimates obtained from order parameter analysis alone, were found to be largely absent.