Research Outputs

Permanent URI for this communityhttps://hdl.handle.net/20.500.14288/2

Browse

Search Results

Now showing 1 - 10 of 26
  • Thumbnail Image
    PublicationOpen Access
    3D spatial organization and network-guided comparison of mutation profiles in Glioblastoma reveals similarities across patients
    (Public Library of Science, 2019) Dinçer, Cansu; Kaya, Tuğba; Tunçbağ, Nurcan; Department of Chemical and Biological Engineering; Department of Computer Engineering; Keskin, Özlem; Gürsoy, Attila; Faculty Member; Department of Chemical and Biological Engineering; Department of Computer Engineering; Koç University Research Center for Translational Medicine (KUTTAM) / Koç Üniversitesi Translasyonel Tıp Araştırma Merkezi (KUTTAM); College of Engineering; 26605; 8745
    Glioblastoma multiforme (GBM) is the most aggressive type of brain tumor. Molecular heterogeneity is a hallmark of GBM tumors that is a barrier in developing treatment strategies. In this study, we used the nonsynonymous mutations of GBM tumors deposited in The Cancer Genome Atlas (TCGA) and applied a systems level approach based on biophysical characteristics of mutations and their organization in patient-specific subnetworks to reduce inter-patient heterogeneity and to gain potential clinically relevant insights. Approximately 10% of the mutations are located in "patches" which are defined as the set of residues spatially in close proximity that are mutated across multiple patients. Grouping mutations as 3D patches reduces the heterogeneity across patients. There are multiple patches that are relatively small in oncogenes, whereas there are a small number of very large patches in tumor suppressors. Additionally, different patches in the same protein are often located at different domains that can mediate different functions. We stratified the patients into five groups based on their potentially affected pathways, revealed from the patient-specific subnetworks. These subnetworks were constructed by integrating mutation profiles of the patients with the interactome data. Network-guided clustering showed significant association between each group and patient survival (P-value = 0.0408). Also, each group carries a set of signature 3D mutation patches that affect predominant pathways. We integrated drug sensitivity data of GBM cell lines with the mutation patches and the patient groups to analyze the therapeutic outcome of these patches. We found that Pazopanib might be effective in Group 3 by targeting CSF1R. Additionally, inhibiting ATM that is a mediator of PTEN phosphorylation may be ineffective in Group 2. We believe that from mutations to networks and eventually to clinical and therapeutic data, this study provides a novel perspective in the network-guided precision medicine.
  • Thumbnail Image
    PublicationOpen Access
    A cartridge based sensor array platform for multiple coagulation measurements from plasma
    (Royal Society of Chemistry (RSC), 2015) Bulut, Serpil; Yaralioglu, G. G.; Department of Electrical and Electronics Engineering; Department of Molecular Biology and Genetics; Department of Chemical and Biological Engineering; Çakmak, Onur; Ermek, Erhan; Kılınç, Necmettin; Barış, İbrahim; Kavaklı, İbrahim Halil; Ürey, Hakan; PhD Student; Other; Researcher; Teaching Faculty; Faculty Member; Department of Electrical and Electronics Engineering; Department of Molecular Biology and Genetics; Department of Chemical and Biological Engineering; College of Engineering; Graduate School of Sciences and Engineering; College of Sciences; N/A; 109991; N/A; 111629; 40319; 8579
    This paper proposes a MEMS-based sensor array enabling multiple clot-time tests for plasma in one disposable microfluidic cartridge. The versatile LoC (Lab-on-Chip) platform technology is demonstrated here for real-time coagulation tests (activated Partial Thromboplastin Time (aPTT) and Prothrombin Time (PT)). The system has a reader unit and a disposable cartridge. The reader has no electrical connections to the cartridge. This enables simple and low-cost cartridge designs and avoids reliability problems associated with electrical connections. The cartridge consists of microfluidic channels and MEMS microcantilevers placed in each channel. The microcantilevers are made of electroplated nickel. They are actuated remotely using an external electro-coil and the read-out is also conducted remotely using a laser. The phase difference between the cantilever oscillation and the coil drive is monitored in real time. During coagulation, the viscosity of the blood plasma increases resulting in a change in the phase read-out. The proposed assay was tested on human and control plasma samples for PT and aPTT measurements. PT and aPTT measurements from control plasma samples are comparable with the manufacturer's datasheet and the commercial reference device. The measurement system has an overall 7.28% and 6.33% CV for PT and aPTT, respectively. For further implementation, the microfluidic channels of the cartridge were functionalized for PT and aPTT tests by drying specific reagents in each channel. Since simultaneous PT and aPTT measurements are needed in order to properly evaluate the coagulation system, one of the most prominent features of the proposed assay is enabling parallel measurement of different coagulation parameters. Additionally, the design of the cartridge and the read-out system as well as the obtained reproducible results with 10 mu l of the plasma samples suggest an opportunity for a possible point-of-care application.
  • Placeholder
    Publication
    A survey of available tools and web servers for analysis of protein-protein interactions and interfaces
    (Oxford University Press (OUP), 2009) Nussinov, Ruth; Department of Chemical and Biological Engineering; Department of Computer Engineering; N/A; Department of Chemical and Biological Engineering; Keskin, Özlem; Gürsoy, Attila; Makinacı, Gözde Kar; Tunçbağ, Nurcan; Faculty Member; Faculty Member; PhD Student; Faculty Member; Department of Computer Engineering; Department of Chemical and Biological Engineering; The Center for Computational Biology and Bioinformatics (CCBB); College of Engineering; College of Engineering; Graduate School of Sciences and Engineering; College of Engineering; 26605; 8745; N/A; 245513
    The unanimous agreement that cellular processes are (largely) governed by interactions between proteins has led to enormous community efforts culminating in overwhelming information relating to these proteins; to the regulation of their interactions, to the way in which they interact and to the function which is determined by these interactions. These data have been organized in databases and servers. However, to make these really useful, it is essential not only to be aware of these, but in particular to have a working knowledge of which tools to use for a given problem; what are the tool advantages and drawbacks; and no less important how to combine these for a particular goal since usually it is not one tool, but some combination of tool-modules that is needed. This is the goal of this review.
  • Placeholder
    Publication
    Analysis of single amino acid variations in singlet hot spots of protein-protein interfaces
    (Oxford Univ Press, 2018) N/A; N/A; Department of Computer Engineering; Department of Chemical and Biological Engineering; Özdemir, E. Sıla; Gürsoy, Attila; Keskin, Özlem; PhD Student; Faculty Member; Faculty Member; Department of Computer Engineering; Department of Chemical and Biological Engineering; Koç University Research Center for Translational Medicine (KUTTAM) / Koç Üniversitesi Translasyonel Tıp Araştırma Merkezi (KUTTAM); Graduate School of Sciences and Engineering; College of Engineering; College of Engineering; N/A; 8745; 26605
    Motivation: Single amino acid variations (SAVs) in protein-protein interaction (PPI) sites play critical roles in diseases. PPI sites (interfaces) have a small subset of residues called hot spots that contribute significantly to the binding energy, and they may form clusters called hot regions. Singlet hot spots are the single amino acid hot spots outside of the hot regions. The distribution of SAVs on the interface residues may be related to their disease association. Results: We performed statistical and structural analyses of SAVs with literature curated experimental thermodynamics data, and demonstrated that SAVs which destabilize PPIs are more likely to be found in singlet hot spots rather than hot regions and energetically less important interface residues. In contrast, non-hot spot residues are significantly enriched in neutral SAVs, which do not affect PPI stability. Surprisingly, we observed that singlet hot spots tend to be enriched in disease-causing SAVs, while benign SAVs significantly occur in non-hot spot residues. Our work demonstrates that SAVs in singlet hot spot residues have significant effect on protein stability and function.
  • Thumbnail Image
    PublicationOpen Access
    Classification of drug molecules considering their IC(50) values using mixed-integer linear programming based hyper-boxes method
    (BioMed Central, 2008) Department of Industrial Engineering; Department of Chemical and Biological Engineering; Armutlu, Pelin; Özdemir, Muhittin Emre; Yüksektepe, Fadime Üney; Kavaklı, İbrahim Halil; Türkay, Metin; Faculty Member; Department of Industrial Engineering; Department of Chemical and Biological Engineering; The Center for Computational Biology and Bioinformatics (CCBB); College of Engineering; N/A; N/A; N/A; 40319; 24956
    Background: A priori analysis of the activity of drugs on the target protein by computational approaches can be useful in narrowing down drug candidates for further experimental tests. Currently, there are a large number of computational methods that predict the activity of drugs on proteins. In this study, we approach the activity prediction problem as a classification problem and, we aim to improve the classification accuracy by introducing an algorithm that combines partial least squares regression with mixed-integer programming based hyper-boxes classification method, where drug molecules are classified as low active or high active regarding their binding activity (IC(50) values) on target proteins. We also aim to determine the most significant molecular descriptors for the drug molecules. Results: We first apply our approach by analyzing the activities of widely known inhibitor datasets including Acetylcholinesterase (ACHE), Benzodiazepine Receptor (BZR), Dihydrofolate Reductase (DHFR), Cyclooxygenase-2 (COX-2) with known IC(50) values. The results at this stage proved that our approach consistently gives better classification accuracies compared to 63 other reported classification methods such as SVM, Naive Bayes, where we were able to predict the experimentally determined IC50 values with a worst case accuracy of 96%. To further test applicability of this approach we first created dataset for Cytochrome P450 C17 inhibitors and then predicted their activities with 100% accuracy. Conclusion: Our results indicate that this approach can be utilized to predict the inhibitory effects of inhibitors based on their molecular descriptors. This approach will not only enhance drug discovery process, but also save time and resources committed.
  • Thumbnail Image
    PublicationOpen Access
    Detection of biological switches using the method of Groebner bases
    (BioMed Central, 2019) Department of Chemical and Biological Engineering; Arkun, Yaman; Faculty Member; Department of Chemical and Biological Engineering; College of Engineering; 108526
    Background: bistability and ability to switch between two stable states is the hallmark of cellular responses. Cellular signaling pathways often contain bistable switches that regulate the transmission of the extracellular information to the nucleus where important biological functions are executed. Results in this work we show how the method of Groebner bases can be used to detect bistability and output switchability. The method of Groebner bases can be seen as a multivariate, non-linear generalization of the Gaussian elimination for linear systems which conveniently seperates the variables and drastically simplifies the simultaneous solution of polynomial equations. A necessary condition for fixed-point state bistability is for the Grobner basis to have three distinct solutions for the state. A sufficient condition is provided by the eigenvalues of the local Jacobians. We also introduce the concept of output switchability which is defined as the ability of an output of a bistable system to switch between two different stable steady-state values. It is shown that bistability does not necessarily guarantee switchability of every state variable of the system. We further show that, for a bistable system, the necessary conditions for output switchability can be derived using the Groebner basis. The theoretical results are incorporated into an analysis procedure and applied to several systems including the AKT (Protein kinase B), RAS (Rat Sarcoma) and MAPK (Mitogen-activated protein kinase) signal transduction pathways. Results demonstrate that the Groebner bases can be conveniently used to analyze biological switches by simultaneously detecting bistability and output switchability. Conclusion: the Groebner bases provides a novel methodology to analyze bistability. Results clarify the distinction between bistability and output switchability which is lacking in the literature. We have shown that theoretically, it is possible to have an output subspace of an n-dimensional bistable system where certain variables cannot switch. It is possible to construct such systems as we have done with two reaction networks.
  • Placeholder
    Publication
    Dynamic correlations: exact and approximate methods for mutual information
    (Oxford Univ Press, 2024) Demirtaş, Kemal; Haliloğlu, Türkan; Department of Chemical and Biological Engineering; Erman, Burak; Department of Chemical and Biological Engineering; College of Engineering
    Motivation Proteins are dynamic entities that undergo conformational changes critical for their functions. Understanding the communication pathways and information transfer within proteins is crucial for elucidating allosteric interactions in their mechanisms. This study utilizes mutual information (MI) analysis to probe dynamic allostery. Using two cases, Ubiquitin and PLpro, we have evaluated the accuracy and limitations of different approximations including the exact anisotropic and isotropic models, multivariate Gaussian model, isotropic Gaussian model, and the Gaussian Network Model (GNM) in revealing allosteric interactions.Results Our findings emphasize the required trajectory length for capturing accurate mutual information profiles. Long molecular dynamics trajectories, 1 ms for Ubiquitin and 100 mu s for PLpro are used as benchmarks, assuming they represent the ground truth. Trajectory lengths of approximately 5 mu s for Ubiquitin and 1 mu s for PLpro marked the onset of convergence, while the multivariate Gaussian model accurately captured mutual information with trajectories of 5 ns for Ubiquitin and 350 ns for PLpro. However, the isotropic Gaussian model is less successful in representing the anisotropic nature of protein dynamics, particularly in the case of PLpro, highlighting its limitations. The GNM, however, provides reasonable approximations of long-range information exchange as a minimalist network model based on a single crystal structure. Overall, the optimum trajectory lengths for effective Gaussian approximations of long-time dynamic behavior depend on the inherent dynamics within the protein's topology. The GNM, by showcasing dynamics across relatively diverse time scales, can be used either as a standalone method or to gauge the adequacy of MD simulation lengths.Availability and implementation Mutual information codes are available at https://github.com/kemaldemirtas/prc-MI.git.
  • Placeholder
    Publication
    Embedding alternative conformations of proteins in protein–protein interaction networks
    (Humana Press inc, 2020) N/A; N/A; Department of Computer Engineering; Department of Chemical and Biological Engineering; Halakou, Farideh; Gürsoy, Attila; Keskin, Özlem; PhD Student; Faculty Member; Faculty Member; Department of Computer Engineering; Department of Chemical and Biological Engineering; Graduate School of Sciences and Engineering; College of Engineering; College of Engineering; N/A; 8745; 26605
    While many proteins act alone, the majority of them interact with others and form molecular complexes to undertake biological functions at both cellular and systems levels. Two proteins should have complementary shapes to physically connect to each other. As proteins are dynamic and changing their conformations, it is vital to track in which conformation a specific interaction can happen. Here, we present a step-by-step guide to embedding the protein alternative conformations in each protein–protein interaction in a systems level. All external tools/websites used in each step are explained, and some notes and suggestions are provided to clear any ambiguous point.
  • Placeholder
    Publication
    Expanding the conformational selection paradigm in protein-ligand docking
    (Humana Press Inc, 2012) Nussinov, Ruth; N/A; Department of Chemical and Biological Engineering; Department of Computer Engineering; Kuzu, Güray; Keskin, Özlem; Gürsoy, Attila; PhD Student; Faculty Member; Faculty Member; N/A; Department of Chemical and Biological Engineering; Department of Computer Engineering; Graduate School of Sciences and Engineering; College of Engineering; College of Engineering; N/A; 26605; 8745
    Conformational selection emerges as a theme in macromolecular interactions. Data validate it as a prevailing mechanism in protein-protein, protein-DNA, protein-RNA, and protein-small molecule drug recognition. This raises the question of whether this fundamental biomolecular binding mechanism can be used to improve drug docking and discovery. Actually, in practice this has already been taking place for some years in increasing numbers. Essentially, it argues for using not a single conformer, but an ensemble. The paradigm of conformational selection holds that because the ensemble is heterogeneous, within it there will be states whose conformation matches that of the ligand. Even if the population of this state is low, since it is favorable for binding the ligand, it will bind to it with a subsequent population shift toward this conformer. Here we suggest expanding it by first modeling all protein interactions in the cell by using Prism, an efficient motif-based protein-protein interaction modeling strategy, followed by ensemble generation. Such a strategy could be particularly useful for signaling proteins, which are major targets in drug discovery and bind multiple partners through a shared binding site, each with some-minor or major-conformational change.
  • Placeholder
    Publication
    Exploiting conformational ensembles in modeling protein-protein interactions on the proteome scale
    (American Chemical Society (ACS), 2013) Nussinov, Ruth; Department of Chemical and Biological Engineering; Department of Computer Engineering; N/A; Keskin, Özlem; Gürsoy, Attila; Kuzu, Güray; Faculty Member; Faculty Member; PhD Student; Department of Chemical and Biological Engineering; Department of Computer Engineering; The Center for Computational Biology and Bioinformatics (CCBB); College of Engineering; College of Engineering; Graduate School of Sciences and Engineering; 26605; 8745; N/A
    Cellular functions are performed through protein-protein interactions; therefore, identification of these interactions is crucial for understanding biological processes. Recent studies suggest that knowledge-based approaches are more useful than "blind" docking for modeling at large scales. However, a caveat of knowledge-based approaches is that they treat molecules as rigid structures. The Protein Data Bank (PDB) offers a wealth of conformations. Here, we exploited an ensemble of the conformations in predictions by a knowledge-based method, PRISM. We tested "difficult" cases in a docking-benchmark data set, where the unbound and bound protein forms are structurally different. Considering alternative conformations for each protein, the percentage of successfully predicted interactions increased from similar to 26 to 66%, and 57% of the interactions were successfully predicted in an "unbiased" scenario, in which data related to the bound forms were not utilized. If the appropriate conformation, or relevant template interface, is unavailable in the PDB, PRISM could not predict the interaction successfully. The pace of the growth of the PDB promises a rapid increase of ensemble conformations emphasizing the merit of such knowledge-based ensemble strategies for higher success rates in protein-protein interaction predictions on an interactome scale. We constructed the structural network of ERK interacting proteins as a case study.