Cooperativity of imprinted genes inactivated by acquired chromosome 20q deletions

Large regions of recurrent genomic loss are common in cancers; however, with a few well-characterized exceptions, how they contribute to tumor pathogenesis remains largely obscure. Here we identified primate-restricted imprinting of a gene cluster on chromosome 20 in the region commonly deleted in chronic myeloid malignancies. We showed that a single heterozygous 20q deletion consistently resulted in the complete loss of expression of the imprinted genes L3MBTL1 and SGK2, indicative of a pathogenetic role for loss of the active paternally inherited locus. Concomitant loss of both L3MBTL1 and SGK2 dysregulated erythropoiesis and megakaryopoiesis, 2 lineages commonly affected in chronic myeloid malignancies, with distinct consequences in each lineage. We demonstrated that L3MBTL1 and SGK2 collaborated in the transcriptional regulation of MYC by influencing different aspects of chromatin structure. L3MBTL1 is known to regulate nucleosomal compaction, and we here showed that SGK2 inactivated BRG1, a key ATP-dependent helicase within the SWI/SNF complex that regulates nucleosomal positioning. These results demonstrate a link between an imprinted gene cluster and malignancy, reveal a new pathogenetic mechanism associated with acquired regions of genomic loss, and underline the complex molecular and cellular consequences of "simple" cancer-associated chromosome deletions.

The Saccharomyces cerevisiae quinone oxidoreductase Lot6p: stability, inhibition and cooperativity

Lot6p (EC 1.5.1.39; Ylr011wp) is the sole quinone oxidoreductase in the budding yeast, Saccharomyces cerevisiae. Using hexahistidine tagged, recombinant Lot6p, we determined the steady-state enzyme kinetic parameters with both NADH and NADPH as electron donors; no cooperativity was observed with these substrates. The NQO1 inhibitor curcumin, the NQO2 inhibitor resveratrol, the bacterial nitroreductase inhibitor nicotinamide and the phosphate mimic vanadate all stabilise the enzyme towards thermal denaturation as judged by differential scanning fluorimetry. All except vanadate have no observable effect on the chemical cross-linking of the two subunits of the Lot6p dimer. These compounds all inhibit Lot6p's oxidoreductase activity, and all except nicotinamide exhibit negative cooperativity. Molecular modelling suggests that curcumin, resveratrol and nicotinamide all bind over the isoalloxazine ring of the FMN cofactor in Lot6p. Resveratrol was predicted to contact an α-helix that links the two active sites. Mutation of Gly-142 (which forms part of this helix) to serine does not greatly affect the thermal stability of the enzyme. However, this variant shows less cooperativity towards resveratrol than the wild type. This suggests a plausible hypothesis for the transmission of information between the subunits and, thus, the molecular mechanism of negative cooperativity in Lot6p.

Investigation of cooperativity in the binding of ligands to the D-2 dopamine receptor

The D 2 dopamine receptor exists as dimers or as higher-order oligomers, as determined from data from physical experiments. In this study, we sought evidence that this oligomerization leads to cooperativity by examining the binding of three radioligands ([H-3] nemonapride, [H-3] raclopride, and [H-3] spiperone) to D 2 dopamine receptors expressed in membranes of Sf9 cells. In saturation binding experiments, the three radioligands exhibited different B-max values, and the B-max values could be altered by the addition of sodium ions to assays. Despite labeling different numbers of sites, the different ligands were able to achieve full inhibition in competition experiments. Some ligand pairs also exhibited complex inhibition curves in these experiments. In radioligand dissociation experiments, the rate of dissociation of [H-3] nemonapride or [H-3] spiperone depended on the sodium ion concentration but was independent of the competing ligand. Although some of the data in this study are consistent with the behavior of a cooperative oligomeric receptor, not all of the data are in agreement with this model. It may, therefore, be necessary to consider more complex models for the behavior of this receptor.

Investigation of cooperativity in the binding of ligands to the D-2 dopamine receptor

The D 2 dopamine receptor exists as dimers or as higher-order oligomers, as determined from data from physical experiments. In this study, we sought evidence that this oligomerization leads to cooperativity by examining the binding of three radioligands ([H-3] nemonapride, [H-3] raclopride, and [H-3] spiperone) to D 2 dopamine receptors expressed in membranes of Sf9 cells. In saturation binding experiments, the three radioligands exhibited different B-max values, and the B-max values could be altered by the addition of sodium ions to assays. Despite labeling different numbers of sites, the different ligands were able to achieve full inhibition in competition experiments. Some ligand pairs also exhibited complex inhibition curves in these experiments. In radioligand dissociation experiments, the rate of dissociation of [H-3] nemonapride or [H-3] spiperone depended on the sodium ion concentration but was independent of the competing ligand. Although some of the data in this study are consistent with the behavior of a cooperative oligomeric receptor, not all of the data are in agreement with this model. It may, therefore, be necessary to consider more complex models for the behavior of this receptor.

Insights into Phosphate Cooperativity and Influence of Substrate Modifications on Binding and Catalysis of Hexameric Purine Nucleoside Phosphorylases

The hexameric purine nucleoside phosphorylase from Bacillus subtilis (BsPNP233) displays great potential to produce nucleoside analogues in industry and can be exploited in the development of new anti-tumor gene therapies. In order to provide structural basis for enzyme and substrates rational optimization, aiming at those applications, the present work shows a thorough and detailed structural description of the binding mode of substrates and nucleoside analogues to the active site of the hexameric BsPNP233. Here we report the crystal structure of BsPNP233 in the apo form and in complex with 11 ligands, including clinically relevant compounds. The crystal structure of six ligands (adenine, 2'deoxyguanosine, aciclovir, ganciclovir, 8-bromoguanosine, 6-chloroguanosine) in complex with a hexameric PNP are presented for the first time. Our data showed that free bases adopt alternative conformations in the BsPNP233 active site and indicated that binding of the co-substrate (2'deoxy) ribose 1-phosphate might contribute for stabilizing the bases in a favorable orientation for catalysis. The BsPNP233-adenosine complex revealed that a hydrogen bond between the 5' hydroxyl group of adenosine and Arg(43*) side chain contributes for the ribosyl radical to adopt an unusual C3'-endo conformation. The structures with 6-chloroguanosine and 8-bromoguanosine pointed out that the Cl-6 and Br-8 substrate modifications seem to be detrimental for catalysis and can be explored in the design of inhibitors for hexameric PNPs from pathogens. Our data also corroborated the competitive inhibition mechanism of hexameric PNPs by tubercidin and suggested that the acyclic nucleoside ganciclovir is a better inhibitor for hexameric PNPs than aciclovir. Furthermore, comparative structural analyses indicated that the replacement of Ser(90) by a threonine in the B. cereus hexameric adenosine phosphorylase (Thr(91)) is responsible for the lack of negative cooperativity of phosphate binding in this enzyme.

Order-Disorder Transitions Govern Kinetic Cooperativity and Allostery of Monomeric Human Glucokinase

Glucokinase (GCK) catalyzes the rate-limiting step of glucose catabolism in the pancreas, where it functions as the body's principal glucose sensor. GCK dysfunction leads to several potentially fatal diseases including maturity-onset diabetes of the young type II (MODY-II) and persistent hypoglycemic hyperinsulinemia of infancy (PHHI). GCK maintains glucose homeostasis by displaying a sigmoidal kinetic response to increasing blood glucose levels. This positive cooperativity is unique because the enzyme functions exclusively as a monomer and possesses only a single glucose binding site. Despite nearly a half century of research, the mechanistic basis for GCK's homotropic allostery remains unresolved. Here we explain GCK cooperativity in terms of large-scale, glucose-mediated disorder-order transitions using 17 isotopically labeled isoleucine methyl groups and three tryptophan side chains as sensitive nuclear magnetic resonance (NMR) probes. We find that the small domain of unliganded GCK is intrinsically disordered and samples a broad conformational ensemble. We also demonstrate that small-molecule diabetes therapeutic agents and hyperinsulinemia-associated GCK mutations share a strikingly similar activation mechanism, characterized by a population shift toward a more narrow, well-ordered ensemble resembling the glucose-bound conformation. Our results support a model in which GCK generates its cooperative kinetic response at low glucose concentrations by using a millisecond disorder-order cycle of the small domain as a "time-delay loop," which is bypassed at high glucose concentrations, providing a unique mechanism to allosterically regulate the activity of human GCK under physiological conditions.

Propagation of melting cooperativity along the phosphodiester backbone of DNA

The role of the DNA phosphodiester backbone in the transfer of melting cooperativity between two helical domains was experimentally addressed with a helix-bulge-helix DNA model, in which the bulge consisted of a varying number of either conformationally flexible propanediol or conformationally constrained bicyclic anucleosidic phosphodiester backbone units. We found that structural communication between two double helical domains is transferred along the DNA backbone over the equivalent of ca. 12-20 backbone units, depending on whether there is a symmetric or asymmetric distribution of the anucleosidic units on both strands. We observed that extension of anucleosidic units on one strand only suffices to disrupt cooperativity transfer in a similar way as if extension occurs on both strands, indicating that the length of the longest anucleosidic inset determines cooperativity transfer. Furthermore, conformational rigidity of the sugar unit increases the distance of coopertivity transfer along the phosphodiester backbone. This is especially the case when the units are asymmetrically distributed in both strands

Analysis of human cytochrome P450 3A4 cooperativity: Construction and characterization of a site-directed mutant that displays hyperbolic steroid hydroxylation kinetics

Cytochrome P450 3A4 is generally considered to be the most important human drug-metabolizing enzyme and is known to catalyze the oxidation of a number of substrates in a cooperative manner. An allosteric mechanism is usually invoked to explain the cooperativity. Based on a structure–activity study from another laboratory using various effector–substrate combinations and on our own studies using site-directed mutagenesis and computer modeling of P450 3A4, the most likely location of effector binding is in the active site along with the substrate. Our study was designed to test this hypothesis by replacing residues Leu-211 and Asp-214 with the larger Phe and Glu, respectively. These residues were predicted to constitute a portion of the effector binding site, and the substitutions were designed to mimic the action of the effector by reducing the size of the active site. The L211F/D214E double mutant displayed an increased rate of testosterone and progesterone 6β-hydroxylation at low substrate concentrations and a decreased level of heterotropic stimulation elicited by α-naphthoflavone. Kinetic analyses of the double mutant revealed the absence of homotropic cooperativity with either steroid substrate. At low substrate concentrations the steroid 6β-hydroxylase activity of the wild-type enzyme was stimulated by a second steroid, whereas L211F/D214E displayed simple substrate inhibition. To analyze L211F/D214E at a more mechanistic level, spectral binding studies were carried out. Testosterone binding by the wild-type enzyme displayed homotropic cooperativity, whereas substrate binding by L211F/D214E displayed hyperbolic behavior.

The unique hetero-oligomeric nature of the subunits in the catalytic cooperativity of the yeast Cct chaperonin complex

The structural and functional organization of the Cct complex was addressed by genetic analyses of subunit interactions and catalytic cooperativity among five of the eight different essential subunits, Cct1p–Cct8p, in the yeast Saccharomyces cerevisiae. The cct1–1, cct2–3, and cct3–1 alleles, containing mutations at the conserved putative ATP-binding motif, GDGTT, are cold-sensitive, whereas single and multiple replacements of the corresponding motif in Cct6p are well tolerated by the cell. We demonstrated herein that cct6–3 (L19S), but not the parolog cct1–5 (R26I), specifically suppresses the cct1–1, cct2–3, and cct3–1 alleles, and that this suppression can be modulated by mutations in a putative phosphorylation motif, RXS, and the putative ATP-binding pocket of Cct6p. Our results suggest that the Cct ring is comprised of a single hetero-oligomer containing eight subunits of differential functional hierarchy, in which catalytic cooperativity of ATP-binding/hydrolysis takes place in a sequential manner different from the concerted cooperativity proposed for GroEL.

Magnesium-dependent folding of self-splicing RNA: Exploring the link between cooperativity, thermodynamics, and kinetics

Folding of the Tetrahymena self-splicing RNA into its active conformation involves a set of discrete intermediate states. The Mg2+-dependent equilibrium transition from the intermediates to the native structure is more cooperative than the formation of the intermediates from the unfolded states. We show that the degree of cooperativity is linked to the free energy of each transition and that the rate of the slow transition from the intermediates to the native state decreases exponentially with increasing Mg2+ concentration. Monovalent salts, which stabilize the folded RNA nonspecifically, induce states that fold in less than 30 s after Mg2+ is added to the RNA. A simple model is proposed that predicts the folding kinetics from the Mg2+-dependent change in the relative stabilities of the intermediate and native states.

I. Quantal effects in biochemical cooperativity and a proposed mechanism for the differentiation of calcium signaling in synaptic plasticity. II. Evolutionary algorithms for the optimization of methods in computational chemistry

In Part 1 of this thesis, we propose that biochemical cooperativity is a fundamentally non-ideal process. We show quantal effects underlying biochemical cooperativity and highlight apparent ergodic breaking at small volumes. The apparent ergodic breaking manifests itself in a divergence of deterministic and stochastic models. We further predict that this divergence of deterministic and stochastic results is a failure of the deterministic methods rather than an issue of stochastic simulations.

Ergodic breaking at small volumes may allow these molecular complexes to function as switches to a greater degree than has previously been shown. We propose that this ergodic breaking is a phenomenon that the synapse might exploit to differentiate Ca$^{2+}$ signaling that would lead to either the strengthening or weakening of a synapse. Techniques such as lattice-based statistics and rule-based modeling are tools that allow us to directly confront this non-ideality. A natural next step to understanding the chemical physics that underlies these processes is to consider \textit{in silico} specifically atomistic simulation methods that might augment our modeling efforts.

In the second part of this thesis, we use evolutionary algorithms to optimize \textit{in silico} methods that might be used to describe biochemical processes at the subcellular and molecular levels. While we have applied evolutionary algorithms to several methods, this thesis will focus on the optimization of charge equilibration methods. Accurate charges are essential to understanding the electrostatic interactions that are involved in ligand binding, as frequently discussed in the first part of this thesis.

Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

Collagen induced MMP-2 activation in human breast cancer

Matrix metalloproteinase-2 (MMP-2), a zymogen requiring proteolytic activation for catalytic activity, has been implicated broadly in the invasion and metastasis of many cancer model systems, including human breast cancer (HBC). MMP-2 has been immunolocalized to carcinomatous human breast, where the degree of activation of MMP-2 correlates well with tumor grade and patient prognosis. Using Matrigel assays, we have stratified HBC cell lines for invasiveness in vitro, and compared this to their potential for metastatic spread in nude mice. HBC cell lines expressing the mesenchymal marker protein vimentin were found to be highly invasive in vitro, and tended to form metastases in nude mice. We have further discovered that culture on collagen-I gels (Vitrogen(TM): Vg) induces MMP-2-activator in highly invasive but not poorly invasive HBC cell lines. As seen for other MMP-2-activator inducing regimens, this induction requires protein synthesis and an intact MMP-2 hemopexin-like domain, appears to be mediated by a cell surface activity, and can be inhibited by metalloproteinase inhibitors. The induction is highly specific to collagen I, and is not seen with thin coatings of collagen I, collagen IV, laminin, or fibronectin, or with 3-dimensional gels of laminin, Matrigel, or gelatin. This review focuses on collagen I and MMP- 2, their localization and source in HBC, and their relationship(s) to MMP-2 activation and HBC metastasis. The relevance of collagen I in activation of MMP-2 in vivo is discussed in terms of stromal cell: tumor cell interaction for collagen I deposition, MMP-2 production and MMP-2-activation. Such cooperativity may exist in vivo for MMP-2 participation in HBC dissemination. A more complete understanding of the regulation of MMP-2-activator by type I collagen may provide new avenues for improved diagnosis and prognosis of human breast cancer.