Studies of pseudoachondroplasia chondrocytes and knockdown of the cartilage oligomeric matrix protein

Cartilage oligomeric matrix protein (COMP) is a large, homopentameric, extracellular matrix glycoprotein. Mutations in COMP cause two skeletal dysplasias: pseudoachondroplasia (PSACH) and multiple epiphyseal dysplasia (EMD1). These dwarfing conditions are caused by retention of misfolded mutant COMP with type IX collagen and matrilin-3 (MATN3) in the rough endoplasmic reticulum (rER) of the chondrocyte. These proteins form a matrix in the rER that continues to expand until it fills the entire cell, eventually causing cell death. Interestingly, loss of COMP in COMP null mice does not affect normal bone development or growth, suggesting that elimination of COMP (wildtype and mutant) expression may prevent PSACH. The hypothesis of these studies was that a hammerhead ribozyme could eliminate or knockdown COMP mRNA expression in PSACH chondrocytes . To test this hypothesis, a human chondrocyte model system that recapitulates the PSACH chondrocyte phenotype was developed by over-expressing mutant (mt-) COMP in normal chondrocytes using a recombinant adenovirus. Chondrocytes over-expressing mt-COMP developed giant rER cisternae containing COMP, type IX collagen and MATN3. Deconvolution microscopy and computer modeling showed that these proteins formed an ordered matrix surrounding a type II pro-collagen core. Additionally, the results show that a hammerhead ribozyme, ribozyme 56 (Ribo56) reduced over-expressed mt-COMP in COS cells and endogenous COMP in normal chondrocytes and mt-COMP in three PSACH chondrocytes cell line (with different mutations) by 40-70%. Altogether, these studies show that the PSACH cellular phenotype can be created in vitro and that the mt-COMP protein burden can be reduced by the presence of a COMP-specific ribozyme. Future studies will focus on designing ribozymes or short interfering RNA (siRNA) technologies that will result in better knockdown of COMP expression as well as the temporal constraints imposed by the PSACH phenotype. ^

Differential Activity of the KRAS Oncogene by Method of Activation: Implications for Signaling and Therapeutic Intervention

Despite having been identified over thirty years ago and definitively established as having a critical role in driving tumor growth and predicting for resistance to therapy, the KRAS oncogene remains a target in cancer for which there is no effective treatment. KRas is activated b y mutations at a few sites, primarily amino acid substitutions at codon 12 which promote a constitutively active state. I have found that different amino acid substitutions at codon 12 can activate different KRas downstream signaling pathways, determine clonogenic growth potential and determine patient response to molecularly targeted therapies. Computer modeling of the KRas structure shows that different amino acids substituted at the codon 12 position influences how KRas interacts with its effecters. In the absence of a direct inhibitor of mutant KRas several agents have recently entered clinical trials alone and in combination directly targeting two of the common downstream effecter pathways of KRas, namely the Mapk pathway and the Akt pathway. These inhibitors were evaluated for efficacy against different KRAS activating mutations. An isogenic panel of colorectal cells with wild type KRas replaced with KRas G12C, G12D, or G12V at the endogenous loci differed in sensitivity to Mek and Akt inhibition. In contrast, screening was performed in a broad panel of lung cell lines alone and no correlation was seen between types of activating KRAS mutation due to concurrent oncogenic lesions. To find a new method to inhibit KRAS driven tumors, siRNA screens were performed in isogenic lines with and without active KRas. The knockdown of CNKSR1 (CNK1) showed selective growth inhibition in cells with an oncogenic KRAS. The deletion of CNK1 reduces expression of mitotic cell cycle proteins and arrests cells with active KRas in the G1 phase of the cell cycle similar to the deletion of an activated KRas regardless of activating substitution. CNK1 has a PH domain responsible for localizing it to membrane lipids making KRas potentially amenable to inhibition with small molecules. The work has identified a series of small molecules capable of binding to this PH domain and inhibiting CNK1 facilitated KRas signaling.

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.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.