Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

CHARACTERIZATION OF SELECTIVE INHIBITION OF ADENYLYL CYCLASE ACTIVITY BY SMALL MOLECULE INHIBITOR NKY80

The nine membrane-bound isoforms of adenylyl cyclase (AC), via synthesis of the signaling molecule cyclic AMP (cAMP), are involved in many isoform specific physiological functions. Decreasing AC5 activity has been shown to have potential therapeutic benefit, including reduced stress on the heart, pain relief, and attenuation of morphine dependence and withdrawal behaviors. However, AC structure is well conserved, and there are currently no isoform selective AC inhibitors in clinical use. P-site inhibitors inhibit AC directly at the catalytic site, but with an uncompetitive or noncompetitive mechanism. Due to this mechanism and nanomolar potency in cell-free systems, attempts at ligand-based drug design of novel AC inhibitors frequently use P-site inhibitors as a starting template. One small molecule inhibitor designed through this process, NKY80, is described as an AC5 selective inhibitor with low micromolar potency in vitro. P-site inhibitors reveal important ligand binding “pockets” in the AC catalytic site, but specific interactions that give NKY80 selectivity are unclear. Identifying and characterizing unique interactions between NKY80 and AC isoforms would significantly aid the development of isoform selective AC inhibitors. I hypothesized that NKY80’s selective inhibition is conferred by AC isoform specific interactions with the compound within the catalytic site. A structure-based virtual screen of the AC catalytic site was used to identify novel small molecule AC inhibitors. Identified novel inhibitors are isoform selective, supporting the catalytic site as a region capable of more potent isoform selective inhibition. Although NKY80 is touted commercially as an AC5 selective inhibitor, its characterization suggests strong inhibition of both AC5 and the closely related AC6. NKY80 was also virtually docked to AC to determine how NKY80 binds to the catalytic site. My results show a difference between NKY80 binding and the conformation of classic P-site inhibitors. The selectivity and notable differences in NKY80 binding to the AC catalytic site suggest a catalytic subregion more flexible in AC5 and AC6 that can be targeted by selective small molecule inhibitors.

NMR structure of a complex between the VirB9/VirB7 interaction domains of the pKM101 type IV secretion system.

Type IV secretion (T4S) systems translocate DNA and protein effectors through the double membrane of Gram-negative bacteria. The paradigmatic T4S system in Agrobacterium tumefaciens is assembled from 11 VirB subunits and VirD4. Two subunits, VirB9 and VirB7, form an important stabilizing complex in the outer membrane. We describe here the NMR structure of a complex between the C-terminal domain of the VirB9 homolog TraO (TraO(CT)), bound to VirB7-like TraN from plasmid pKM101. TraO(CT) forms a beta-sandwich around which TraN winds. Structure-based mutations in VirB7 and VirB9 of A. tumefaciens show that the heterodimer interface is conserved. Opposite this interface, the TraO structure shows a protruding three-stranded beta-appendage, and here, we supply evidence that the corresponding region of VirB9 of A. tumefaciens inserts in the membrane and protrudes extracellularly. This complex structure elucidates the molecular basis for the interaction between two essential components of a T4S system.

Computational selection of inhibitors of Abeta aggregation and neuronal toxicity.

Alzheimer's disease (AD) is characterized by the cerebral accumulation of misfolded and aggregated amyloid-beta protein (Abeta). Disease symptoms can be alleviated, in vitro and in vivo, by 'beta-sheet breaker' pentapeptides that reduce plaque load. However the peptide nature of these compounds, made them biologically unstable and unable to penetrate membranes with high efficiency. The main goal of this study was to use computational methods to identify small molecule mimetics with better drug-like properties. For this purpose, the docked conformations of the active peptides were used to identify compounds with similar activities. A series of related beta-sheet breaker peptides were docked to solid state NMR structures of a fibrillar form of Abeta. The lowest energy conformations of the active peptides were used to design three dimensional (3D)-pharmacophores, suitable for screening the NCI database with Unity. Small molecular weight compounds with physicochemical features and a conformation similar to the active peptides were selected, ranked by docking and biochemical parameters. Of 16 diverse compounds selected for experimental screening, 2 prevented and reversed Abeta aggregation at 2-3microM concentration, as measured by Thioflavin T (ThT) fluorescence and ELISA assays. They also prevented the toxic effects of aggregated Abeta on neuroblastoma cells. Their low molecular weight and aqueous solubility makes them promising lead compounds for treating AD.

Targeted Multistage Delivery of Nanoparticles to the Bone Marrow

Bone marrow is a target organ site involved in multiple diseases including myeloproliferative disorders and hematologic malignancies and metastases from breast and prostate. Most of these diseases are characterized with poor quality of life, and the treatment options are only palliative due to lack of delivery mechanisms for systemically injected drugs which results in dose limitation to protect the healthy hematopoietic cells. Therefore, there is a critical need to develop effective therapeutic strategies that allow for selective delivery of therapeutic payload to the bone marrow. Nanotechnology-based drug delivery systems provide the opportunity to deliver drugs to the target tissue while decreasing exposure to normal tissues. E-selectin is constitutively expressed on the bone marrow vasculature, but almost absent in normal vessels, and therefore, E-selectin targeted drug delivery presents an ideal strategy for the delivery of therapeutic nanoparticles to the bone marrow. The objective of this study was to develop a novel bone marrow targeted multistage vector (MSV) via E-selectin for delivery of therapeutics and imaging agents. To achieve this goal, Firstly, an E-selectin thioaptamer (ESTA) ligand was identified through a two-step screening from a combinatorial thioaptamer library. Next, ESTA-conjugated MSV (ESTA-MSV) were developed and evaluated for their stability and binding to E-selectin expressing endothelial cells. Different types of nanoparticles including liposomes, quantum dots, and iron oxide nanoparticles were loaded into the porous structure of ESTA-MSV. In vivo targeting experiments demonstrated 8-fold higher accumulation of ESTA-MSV in the mouse bone marrow as compared to non-targeted MSV Furthermore, intravenous injection of liposomes loaded ESTA-MSV resulted in a significantly higher accumulation of liposome in the bone marrow space as compared to injection of non-targeted MSV or liposomes alone. Overall this study provides first evidence that E-selectin targeted multistage vector preferentially targets to bone marrow vasculature and delivers larger amounts of nanoparticles. This delivery strategy holds potential for the selective delivery of large amounts of therapeutic payload to the vascular niches in the bone marrow for the treatment of bone marrow associated diseases.

Targeting TRAF6 for Cancer Therapeutical Development

Tumor necrosis factor (TNF)-Receptor Associated Factors (TRAFs) are a family of signal transducer proteins. TRAF6 is a unique member of this family in that it is involved in not only the TNF superfamily, but the toll-like receptor (TLR)/IL-1R (TIR) superfamily. The formation of the complex consisting of Receptor Activator of Nuclear Factor κ B (RANK), with its ligand (RANKL) results in the recruitment of TRAF6, which activates NF-κB, JNK and MAP kinase pathways. TRAF6 is critical in signaling with leading to release of various growth factors in bone, and promotes osteoclastogenesis. TRAF6 has also been implicated as an oncogene in lung cancer and as a target in multiple myeloma. In the hopes of developing small molecule inhibitors of the TRAF6-RANK interaction, multiple steps were carried out. Computational prediction of hot spot residues on the protein-protein interaction of TRAF6 and RANK were examined. Three methods were used: Robetta, KFC2, and HotPoint, each of which uses a different methodology to determine if a residue is a hot spot. These hot spot predictions were considered the basis for resolving the binding site for in silico high-throughput screening using GOLD and the MyriaScreen database of drug/lead-like compounds. Computationally intensive molecular dynamics simulations highlighted the binding mechanism and TRAF6 structural changes upon hit binding. Compounds identified as hits were verified using a GST-pull down assay, comparing inhibition to a RANK decoy peptide. Since many drugs fail due to lack of efficacy and toxicity, predictive models for the evaluation of the LD50 and bioavailability of our TRAF6 hits, and these models can be used towards other drugs and small molecule therapeutics as well. Datasets of compounds and their corresponding bioavailability and LD50 values were curated based, and QSAR models were built using molecular descriptors of these compounds using the k-nearest neighbor (k-NN) method, and quality of these models were cross-validated.

Development of HIF-1α/HIF-1β heterodimerization inhibitors using a novel bioluminescence reporter assay system for in vitro high throughput screening and in vivo imaging

Tumor growth often outpaces its vascularization, leading to development of a hypoxic tumor microenvironment. In response, an intracellular hypoxia survival pathway is initiated by heterodimerization of hypoxia-inducible factor (HIF)-1α and HIF-1β, which subsequently upregulates the expression of several hypoxia-inducible genes, promotes cell survival and stimulates angiogenesis in the oxygen-deprived environment. Hypoxic tumor regions are often associated with resistance to various classes of radio- or chemotherapeutic agents. Therefore, development of HIF-1α/β heterodimerization inhibitors may provide a novel approach to anti-cancer therapy. To this end, a novel approach for imaging HIF-1α/β heterodimerization in vitro and in vivo was developed in this study. Using this screening platform, we identified a promising lead candidate and further chemically derivatized the lead candidate to assess the structure-activity relationship (SAR). The most effective first generation drug inhibitors were selected and their pharmacodynamics and anti-tumor efficacy in vivo were verified by bioluminescence imaging (BLI) of HIF-1α/β heterodimerization in the xenograft tumor model. Furthermore, the first generation drug inhibitors, M-TMCP and D-TMCP, demonstrated efficacy as monotherapies, resulting in tumor growth inhibition via disruption of HIF-1 signaling-mediated tumor stromal neoangiogenesis.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Cryo-EM model of the bullet-shaped vesicular stomatitis virus.

Vesicular stomatitis virus (VSV) is a bullet-shaped rhabdovirus and a model system of negative-strand RNA viruses. Through direct visualization by means of cryo-electron microscopy, we show that each virion contains two nested, left-handed helices: an outer helix of matrix protein M and an inner helix of nucleoprotein N and RNA. M has a hub domain with four contact sites that link to neighboring M and N subunits, providing rigidity by clamping adjacent turns of the nucleocapsid. Side-by-side interactions between neighboring N subunits are critical for the nucleocapsid to form a bullet shape, and structure-based mutagenesis results support this description. Together, our data suggest a mechanism of VSV assembly in which the nucleocapsid spirals from the tip to become the helical trunk, both subsequently framed and rigidified by the M layer.

Expression and Regulation of Human Cytochrome P450 4F Isoforms in Tissue Samples and Under TNF-Alpha Challenges

CYP4F (Cytochrome P4504F) enzymes metabolize endogenous molecules including leukotrienes, prostaglandins and arachidonic acid. The involvement of these endogenous compounds in inflammation has led to the hypothesis that changes in the inflamed tissue environment may affect the expression of CYP4Fs during the pro-inflammatory state, which in turn may modulate inflammatory conditions during the anti-inflammatory state. We demonstrated that inflamed tissues have different levels of CYP4F isoform expression profiles in a number of human samples when compared to the average population. The CYP4F isoform expression levels change with the degree of inflammation present in tissue. Further investigation in cell culture studies revealed that inflammatory cytokines, in particular TNF-α, play a role in regulating the expression of the CYP4F family. One of the isoforms, CYP4F11, had different characteristics than that of the other five CYP4F family members. CYP4F11 metabolizes xenobiotics while the other isoforms metabolize endogenous compounds with higher affinity. CYP4F11 also was expressed at high quantities in the brain, and was up-regulated by TNF-α, while the other isoforms were not expressed at high quantities in the brain and were down-regulated by TNF-α. We identified the AP-1 protein of the JNK pathway as the signaling protein that causes significant increase in CYP4F11 expression. Since TNF-α stimulation causes a simultaneous activation of both JNK pathway and NF-κB signaling, we investigated further the role that NF-κB plays on expression of the CYP4F11 gene. We concluded that although there is a significant increase in CYP4F11 expression in the presence of TNF-α, the activation of NF-κB signaling inhibits CYP4F11 expression in a time dependent manner. The expression of CYP4F11 is only significantly increased after 24 hours of treatment with TNF-α; at shorter time points NF-κB signaling overpowers the JNK pathway activation. We believe that these findings may in the future lead to improved drug design for modulating inflammation.

THE HISTONE DEACETYLASE INHIBITOR, MS-275, SENSITIZES METASTATIC OSTEOSARCOMA TO FASL-INDUCED CELL DEATH: A ROLE FOR C-FLIP

The purpose of this study was to determine the effects of the histone deacetylase inhibitor, MS-275, on the Fas signaling pathway and susceptibility of osteosarcoma (OS) to Fas ligand (FasL)-induced cell death. OS metastasizes almost exclusively to the lungs. We have shown that Fas expression in OS cells is inversely correlated with their metastatic potential. Fas+ cells are rapidly eliminated when they enter the lungs via interaction with FasL, which is constitutively expressed in the lungs. Fas- OS cells escape this FasL-induced apoptosis and survive in the lung microenvironment. Moreover, upregulation of Fas in established OS lung metastases results in tumor regression. Therefore, agents that upregulate Fas expression or activate the Fas signaling pathway may have therapeutic potential. Treatment of Fas- metastatic OS cell lines with 2 μM MS-275 sensitized cells to FasL-induced cell death in vitro. We found that MS-275 did not alter the expression of Fas on the cell surface; rather it resulted in increased levels of Fas within the membrane lipid rafts, as demonstrated by an increase in Fas expression in detergent insoluble lipid raft fractions. We further demonstrated that following MS-275 treatment, Fas colocalized with GM1+ lipid rafts and that there was a decrease in c-FLIP (cellular FLICE-inhibitory protein) mRNA and protein. Downregulation of c-FLIP correlated with caspase activation and apoptosis induction. Transfection of cells with shRNA to c-FLIP also resulted in the localization of Fas to lipid rafts. These studies indicate that MS-275 sensitizes OS cells to FasL by upregulating the expression of Fas in membrane lipid rafts, which correlated with the downregulation of c-FLIP. Treatment of nu/nu-mice with established OS lung metastases with oral MS-275 resulted in increased apoptosis, a significant inhibition of c-FLIP expression in tumors and tumor regression. Histopathological examination of mice showed no significant organ toxicity. Overall, these results suggest that the mechanism by which MS-275 sensitizes OS cells and lung metastases to FasL-induced cell death may be by a reduction in the expression of c-FLIP.

CIRCADIAN RHYTHM IN THE BIOAVAILABILITY OF ORALLY ADMINISTERED THEOPHYLLINE

The effect of time (i.e., biologic time structure) of drug administration on the bioavailability of theophylline was investigated in man after both a single dosage as well as after repeated, or chronic, drug administrations. Preliminary laboratory investigations on Balb-C mice showed the toxic

On Bayesian seamless phase I-II designs

Early phase clinical trial designs have long been the focus of interest for clinicians and statisticians working in oncology field. There are several standard phse I and phase II designs that have been widely-implemented in medical practice. For phase I design, the most commonly used methods are 3+3 and CRM. A newly-developed Bayesian model-based mTPI design has now been used by an increasing number of hospitals and pharmaceutical companies. The advantages and disadvantages of these three top phase I designs have been discussed in my work here and their performances were compared using simulated data. It was shown that mTPI design exhibited superior performance in most scenarios in comparison with 3+3 and CRM designs. ^ The next major part of my work is proposing an innovative seamless phase I/II design that allows clinicians to conduct phase I and phase II clinical trials simultaneously. Bayesian framework was implemented throughout the whole design. The phase I portion of the design adopts mTPI method, with the addition of futility rule which monitors the efficacy performance of the tested drugs. Dose graduation rules were proposed in this design to allow doses move forward from phase I portion of the study to phase II portion without interrupting the ongoing phase I dose-finding schema. Once a dose graduated to phase II, adaptive randomization was used to randomly allocated patients into different treatment arms, with the intention of more patients being assigned to receive more promising dose(s). Again simulations were performed to compare the performance of this innovative phase I/II design with a recently published phase I/II design, together with the conventional phase I and phase II designs. The simulation results indicated that the seamless phase I/II design outperform the other two competing methods in most scenarios, with superior trial power and the fact that it requires smaller sample size. It also significantly reduces the overall study time. ^ Similar to other early phase clinical trial designs, the proposed seamless phase I/II design requires that the efficacy and safety outcomes being able to be observed in a short time frame. This limitation can be overcome by using validated surrogate marker for the efficacy and safety endpoints.^