Academic Appointments

  • Sr Research Engineer, Bioengineering

Professional Education

  • B.S., Peking University, China, Biochemistry and Molecular Biology (1997)
  • M.S., Peking University, China, Biochemistry and Molecular Biology (2000)
  • Ph.D., University of Washington, Seattle, Computational Biology (2005)

All Publications

  • Pocket similarity identifies selective estrogen receptor modulators as microtubule modulators at the taxane site. Nature communications Lo, Y. C., Cormier, O. n., Liu, T. n., Nettles, K. W., Katzenellenbogen, J. A., Stearns, T. n., Altman, R. B. 2019; 10 (1): 1033


    Taxanes are a family of natural products with a broad spectrum of anticancer activity. This activity is mediated by interaction with the taxane site of beta-tubulin, leading to microtubule stabilization and cell death. Although widely used in the treatment of breast cancer and other malignancies, existing taxane-based therapies including paclitaxel and the second-generation docetaxel are currently limited by severe adverse effects and dose-limiting toxicity. To discover taxane site modulators, we employ a computational binding site similarity screen of > 14,000 drug-like pockets from PDB, revealing an unexpected similarity between the estrogen receptor and the beta-tubulin taxane binding pocket. Evaluation of nine selective estrogen receptor modulators (SERMs) via cellular and biochemical assays confirms taxane site interaction, microtubule stabilization, and cell proliferation inhibition. Our study demonstrates that SERMs can modulate microtubule assembly and raises the possibility of an estrogen receptor-independent mechanism for inhibiting cell proliferation.

    View details for PubMedID 30833575

  • PathFX provides mechanistic insights into drug efficacy and safety for regulatory review and therapeutic development. PLoS computational biology Wilson, J. L., Racz, R., Liu, T., Adeniyi, O., Sun, J., Ramamoorthy, A., Pacanowski, M., Altman, R. 2018; 14 (12): e1006614


    Failure to demonstrate efficacy and safety issues are important reasons that drugs do not reach the market. An incomplete understanding of how drugs exert their effects hinders regulatory and pharmaceutical industry projections of a drug's benefits and risks. Signaling pathways mediate drug response and while many signaling molecules have been characterized for their contribution to disease or their role in drug side effects, our knowledge of these pathways is incomplete. To better understand all signaling molecules involved in drug response and the phenotype associations of these molecules, we created a novel method, PathFX, a non-commercial entity, to identify these pathways and drug-related phenotypes. We benchmarked PathFX by identifying drugs' marketed disease indications and reported a sensitivity of 41%, a 2.7-fold improvement over similar approaches. We then used PathFX to strengthen signals for drug-adverse event pairs occurring in the FDA Adverse Event Reporting System (FAERS) and also identified opportunities for drug repurposing for new diseases based on interaction paths that associated a marketed drug to that disease. By discovering molecular interaction pathways, PathFX improved our understanding of drug associations to safety and efficacy phenotypes. The algorithm may provide a new means to improve regulatory and therapeutic development decisions.

    View details for PubMedID 30532240

  • The Pioglitazone Trek via Human PPAR Gamma: From Discovery to a Medicine at the FDA and Beyond FRONTIERS IN PHARMACOLOGY Devchand, P. R., Liu, T., Altman, R. B., FitzGerald, G. A., Schadt, E. E. 2018; 9
  • Computational Analysis of Kinase Inhibitor Selectivity using Structural Knowledge. Bioinformatics (Oxford, England) Lo, Y., Liu, T., Morrissey, K. M., Kakiuchi-Kiyota, S., Johnson, A. R., Broccatelli, F., Zhong, Y., Joshi, A., Altman, R. B. 2018


    Motivation: Kinases play a significant role in diverse disease signaling pathways and understanding kinase inhibitor selectivity, the tendency of drugs to bind to off-targets, remains a top priority for kinase inhibitor design and clinical safety assessment. Traditional approaches for kinase selectivity analysis using biochemical activity and binding assays are useful but can be costly and are often limited by the kinases that are available. On the other hand, current computational kinase selectivity prediction methods are computational intensive and can rarely achieve sufficient accuracy for large-scale kinome wide inhibitor selectivity profiling.Results: Here, we present a KinomeFEATURE database for kinase binding site similarity search by comparing protein microenvironments characterized using diverse physiochemical descriptors. Initial selectivity prediction of 15 known kinase inhibitors achieved an>90% accuracy and demonstrated improved performance in comparison to commonly used kinase inhibitor selectivity prediction methods. Additional kinase ATP binding site similarity assessment (120 binding sites) identified 55 kinases with significant promiscuity and revealed unexpected inhibitor cross-activities between PKR and FGFR2 kinases. Kinome-wide selectivity profiling of 11 kinase drug candidates predicted novel as well as experimentally validated off-targets and suggested structural mechanisms of kinase cross-activities. Our study demonstrated potential utilities of our approach for large-scale kinase inhibitor selectivity profiling that could contribute to kinase drug development and safety assessment.Availability: The KinomeFEATURE database are available at information: Supplementary data are available at Bioinformatics online.

    View details for PubMedID 29985971

  • Biological and functional relevance of CASP predictions PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS Liu, T., Ish-Shalom, S., Torng, W., Lafita, A., Bock, C., Mort, M., Cooper, D. N., Bliven, S., Capitani, G., Mooney, S. D., Altman, R. B. 2018; 86: 374–86


    Our goal is to answer the question: compared with experimental structures, how useful are predicted models for functional annotation? We assessed the functional utility of predicted models by comparing the performances of a suite of methods for functional characterization on the predictions and the experimental structures. We identified 28 sites in 25 protein targets to perform functional assessment. These 28 sites included nine sites with known ligand binding (holo-sites), nine sites that are expected or suggested by experimental authors for small molecule binding (apo-sites), and Ten sites containing important motifs, loops, or key residues with important disease-associated mutations. We evaluated the utility of the predictions by comparing their microenvironments to the experimental structures. Overall structural quality correlates with functional utility. However, the best-ranked predictions (global) may not have the best functional quality (local). Our assessment provides an ability to discriminate between predictions with high structural quality. When assessing ligand-binding sites, most prediction methods have higher performance on apo-sites than holo-sites. Some servers show consistently high performance for certain types of functional sites. Finally, many functional sites are associated with protein-protein interaction. We also analyzed biologically relevant features from the protein assemblies of two targets where the active site spanned the protein-protein interface. For the assembly targets, we find that the features in the models are mainly determined by the choice of template.

    View details for PubMedID 28975675

    View details for PubMedCentralID PMC5820171

  • The Pioglitazone Trek via Human PPAR Gamma: From Discovery to a Medicine at the FDA and Beyond. Frontiers in pharmacology Devchand, P. R., Liu, T. n., Altman, R. B., FitzGerald, G. A., Schadt, E. E. 2018; 9: 1093


    For almost two decades, pioglitazone has been prescribed primarily to prevent and treat insulin resistance in some type 2 diabetic patients. In this review, we trace the path to discovery of pioglitazone as a thiazolidinedione compound, the glitazone tracks through the regulatory agencies, the trek to molecular agonism in the nucleus and the binding of pioglitazone to the nuclear receptor PPAR gamma. Given the rise in consumption of pioglitazone in T2D patients worldwide and the increased number of clinical trials currently testing alternate medical uses for this drug, there is also merit to some reflection on the reported adverse effects. Going forward, it is imperative to continue investigations into the mechanisms of actions of pioglitazone, the potential of glitazone drugs to contribute to unmet needs in complex diseases associated with the dynamics of adaptive homeostasis, and also the routes to minimizing adverse effects in every-day patients throughout the world.

    View details for PubMedID 30337873

  • Estimation of Maximum Recommended Therapeutic Dose Using Predicted Promiscuity and Potency. Clinical and translational science Liu, T., Oprea, T., Ursu, O., Hasselgren, C., Altman, R. B. 2016


    We report a simple model that predicts the maximum recommended therapeutic dose (MRTD) of small molecule drugs based on an assessment of likely protein-drug interactions. Previously, we reported methods for computational estimation of drug promiscuity and potency. We used these concepts to build a linear model derived from 238 small molecular drugs to predict MRTD. We applied this model successfully to predict MRTDs for 16 nonsteroidal antiinflammatory drugs (NSAIDs) and 14 antiretroviral drugs. Of note, based on the estimated promiscuity of low-dose drugs (and active chemicals), we identified 83 proteins as "high-risk off-targets" (HROTs) that are often associated with low doses; the evaluation of interactions with HROTs may be useful during early phases of drug discovery. Our model helps explain the MRTD for drugs with severe adverse reactions caused by interactions with HROTs.

    View details for DOI 10.1111/cts.12422

    View details for PubMedID 27736015

    View details for PubMedCentralID PMC5161261

  • Relating Essential Proteins to Drug Side-Effects Using Canonical Component Analysis: A Structure-Based Approach. Journal of chemical information and modeling Liu, T., Altman, R. B. 2015; 55 (7): 1483-1494


    The molecular mechanism of many drug side-effects is unknown and difficult to predict. Previous methods for explaining side-effects have focused on known drug targets and their pathways. However, low affinity binding to proteins that are not usually considered drug targets may also drive side-effects. In order to assess these alternative targets, we used the 3D structures of 563 essential human proteins systematically to predict binding to 216 drugs. We first benchmarked our affinity predictions with available experimental data. We then combined singular value decomposition and canonical component analysis (SVD-CCA) to predict side-effects based on these novel target profiles. Our method predicts side-effects with good accuracy (average AUC: 0.82 for side effects present in <50% of drug labels). We also noted that side-effect frequency is the most important feature for prediction and can confound efforts at elucidating mechanism; our method allows us to remove the contribution of frequency and isolate novel biological signals. In particular, our analysis produces 2768 triplet associations between 50 essential proteins, 99 drugs, and 77 side-effects. Although experimental validation is difficult because many of our essential proteins do not have validated assays, we nevertheless attempted to validate a subset of these associations using experimental assay data. Our focus on essential proteins allows us to find potential associations that would likely be missed if we used recognized drug targets. Our associations provide novel insights about the molecular mechanisms of drug side-effects and highlight the need for expanded experimental efforts to investigate drug binding to proteins more broadly.

    View details for DOI 10.1021/acs.jcim.5b00030

    View details for PubMedID 26121262

  • Variations in the binding pocket of an inhibitor of the bacterial division protein FtsZ across genotypes and species. PLoS computational biology Miguel, A., Hsin, J., Liu, T., Tang, G., Altman, R. B., Huang, K. C. 2015; 11 (3)


    The recent increase in antibiotic resistance in pathogenic bacteria calls for new approaches to drug-target selection and drug development. Targeting the mechanisms of action of proteins involved in bacterial cell division bypasses problems associated with increasingly ineffective variants of older antibiotics; to this end, the essential bacterial cytoskeletal protein FtsZ is a promising target. Recent work on its allosteric inhibitor, PC190723, revealed in vitro activity on Staphylococcus aureus FtsZ and in vivo antimicrobial activities. However, the mechanism of drug action and its effect on FtsZ in other bacterial species are unclear. Here, we examine the structural environment of the PC190723 binding pocket using PocketFEATURE, a statistical method that scores the similarity between pairs of small-molecule binding sites based on 3D structure information about the local microenvironment, and molecular dynamics (MD) simulations. We observed that species and nucleotide-binding state have significant impacts on the structural properties of the binding site, with substantially disparate microenvironments for bacterial species not from the Staphylococcus genus. Based on PocketFEATURE analysis of MD simulations of S. aureus FtsZ bound to GTP or with mutations that are known to confer PC190723 resistance, we predict that PC190723 strongly prefers to bind Staphylococcus FtsZ in the nucleotide-bound state. Furthermore, MD simulations of an FtsZ dimer indicated that polymerization may enhance PC190723 binding. Taken together, our results demonstrate that a drug-binding pocket can vary significantly across species, genetic perturbations, and in different polymerization states, yielding important information for the further development of FtsZ inhibitors.

    View details for DOI 10.1371/journal.pcbi.1004117

    View details for PubMedID 25811761

    View details for PubMedCentralID PMC4374959

  • "Genotype-first" approaches on a curious case of idiopathic progressive cognitive decline BMC MEDICAL GENOMICS Shi, L., Li, B., Huang, Y., Ling, X., Liu, T., Lyon, G. J., Xu, A., Wang, K. 2014; 7


    In developing countries, many cases with rare neurological diseases remain undiagnosed due to limited diagnostic experience. We encountered a case in China where two siblings both began to develop idiopathic progressive cognitive decline starting from age six, and were suspected to have an undiagnosed neurological disease.Initial clinical assessments included review of medical history, comprehensive physical examination, genetic testing for metabolic diseases, blood tests and brain imaging. We performed exome sequencing with Agilent SureSelect exon capture and Illumina HiSeq2000 platform, followed by variant annotation and selection of rare, shared mutations that fit a recessive model of inheritance. To assess functional impacts of candidate variants, we performed extensive biochemical tests in blood and urine, and examined their possible roles by protein structure modeling.Exome sequencing identified NAGLU as the most likely candidate gene with compound heterozygous mutations (chr17:40695717C > T and chr17:40693129A > G in hg19 coordinate), which were documented to be pathogenic. Sanger sequencing confirmed the recessive patterns of inheritance, leading to a genetic diagnosis of Sanfilippo syndrome (mucopolysaccharidosis IIIB). Biochemical tests confirmed the complete loss of activity of alpha-N-acetylglucosaminidase (encoded by NAGLU) in blood, as well as significantly elevated dermatan sulfate and heparan sulfate in urine. Structure modeling revealed the mechanism on how the two variants affect protein structural stability.Successful diagnosis of a rare genetic disorder with an atypical phenotypic presentation confirmed that such "genotype-first" approaches can particularly succeed in areas of the world with insufficient medical genetics expertise and with cost-prohibitive in-depth phenotyping.

    View details for DOI 10.1186/s12920-014-0066-9

    View details for Web of Science ID 000348397200001

    View details for PubMedCentralID PMC4267425

  • Identifying druggable targets by protein microenvironments matching: application to transcription factors. CPT: pharmacometrics & systems pharmacology Liu, T., Altman, R. B. 2014; 3


    Druggability of a protein is its potential to be modulated by drug-like molecules. It is important in the target selection phase. We hypothesize that: (i) known drug-binding sites contain advantageous physicochemical properties for drug binding, or "druggable microenvironments" and (ii) given a target, the presence of multiple druggable microenvironments similar to those seen previously is associated with a high likelihood of druggability. We developed DrugFEATURE to quantify druggability by assessing the microenvironments in potential small-molecule binding sites. We benchmarked DrugFEATURE using two data sets. One data set measures druggability using NMR-based screening. DrugFEATURE correlates well with this metric. The second data set is based on historical drug discovery outcomes. Using the DrugFEATURE cutoffs derived from the first, we accurately discriminated druggable and difficult targets in the second. We further identified novel druggable transcription factors with implications for cancer therapy. DrugFEATURE provides useful insight for drug discovery, by evaluating druggability and suggesting specific regions for interacting with drug-like molecules.CPT: Pharmacometrics Systems Pharmacology (2014) 3, e93; doi:10.1038/psp.2013.66; published online 22 January 2014.

    View details for DOI 10.1038/psp.2013.66

    View details for PubMedID 24452614

    View details for PubMedCentralID PMC3910014

  • Bioinformatics and variability in drug response: a protein structural perspective JOURNAL OF THE ROYAL SOCIETY INTERFACE Lahti, J. L., Tang, G. W., Capriotti, E., Liu, T., Altman, R. B. 2012; 9 (72): 1409-1437


    Marketed drugs frequently perform worse in clinical practice than in the clinical trials on which their approval is based. Many therapeutic compounds are ineffective for a large subpopulation of patients to whom they are prescribed; worse, a significant fraction of patients experience adverse effects more severe than anticipated. The unacceptable risk-benefit profile for many drugs mandates a paradigm shift towards personalized medicine. However, prior to adoption of patient-specific approaches, it is useful to understand the molecular details underlying variable drug response among diverse patient populations. Over the past decade, progress in structural genomics led to an explosion of available three-dimensional structures of drug target proteins while efforts in pharmacogenetics offered insights into polymorphisms correlated with differential therapeutic outcomes. Together these advances provide the opportunity to examine how altered protein structures arising from genetic differences affect protein-drug interactions and, ultimately, drug response. In this review, we first summarize structural characteristics of protein targets and common mechanisms of drug interactions. Next, we describe the impact of coding mutations on protein structures and drug response. Finally, we highlight tools for analysing protein structures and protein-drug interactions and discuss their application for understanding altered drug responses associated with protein structural variants.

    View details for DOI 10.1098/rsif.2011.0843

    View details for Web of Science ID 000304437400001

    View details for PubMedID 22552919

  • Using Multiple Microenvironments to Find Similar Ligand-Binding Sites: Application to Kinase Inhibitor Binding PLOS COMPUTATIONAL BIOLOGY Liu, T., Altman, R. B. 2011; 7 (12)


    The recognition of cryptic small-molecular binding sites in protein structures is important for understanding off-target side effects and for recognizing potential new indications for existing drugs. Current methods focus on the geometry and detailed chemical interactions within putative binding pockets, but may not recognize distant similarities where dynamics or modified interactions allow one ligand to bind apparently divergent binding pockets. In this paper, we introduce an algorithm that seeks similar microenvironments within two binding sites, and assesses overall binding site similarity by the presence of multiple shared microenvironments. The method has relatively weak geometric requirements (to allow for conformational change or dynamics in both the ligand and the pocket) and uses multiple biophysical and biochemical measures to characterize the microenvironments (to allow for diverse modes of ligand binding). We term the algorithm PocketFEATURE, since it focuses on pockets using the FEATURE system for characterizing microenvironments. We validate PocketFEATURE first by showing that it can better discriminate sites that bind similar ligands from those that do not, and by showing that we can recognize FAD-binding sites on a proteome scale with Area Under the Curve (AUC) of 92%. We then apply PocketFEATURE to evolutionarily distant kinases, for which the method recognizes several proven distant relationships, and predicts unexpected shared ligand binding. Using experimental data from ChEMBL and Ambit, we show that at high significance level, 40 kinase pairs are predicted to share ligands. Some of these pairs offer new opportunities for inhibiting two proteins in a single pathway.

    View details for DOI 10.1371/journal.pcbi.1002326

    View details for Web of Science ID 000299167800043

    View details for PubMedID 22219723

    View details for PubMedCentralID PMC3248393

  • Comparative Modeling: The State of the Art and Protein Drug Target Structure Prediction COMBINATORIAL CHEMISTRY & HIGH THROUGHPUT SCREENING Liu, T., Tang, G. W., Capriotti, E. 2011; 14 (6): 532-547


    The goal of computational protein structure prediction is to provide three-dimensional (3D) structures with resolution comparable to experimental results. Comparative modeling, which predicts the 3D structure of a protein based on its sequence similarity to homologous structures, is the most accurate computational method for structure prediction. In the last two decades, significant progress has been made on comparative modeling methods. Using the large number of protein structures deposited in the Protein Data Bank (~65,000), automatic prediction pipelines are generating a tremendous number of models (~1.9 million) for sequences whose structures have not been experimentally determined. Accurate models are suitable for a wide range of applications, such as prediction of protein binding sites, prediction of the effect of protein mutations, and structure-guided virtual screening. In particular, comparative modeling has enabled structure-based drug design against protein targets with unknown structures. In this review, we describe the theoretical basis of comparative modeling, the available automatic methods and databases, and the algorithms to evaluate the accuracy of predicted structures. Finally, we discuss relevant applications in the prediction of important drug target proteins, focusing on the G protein-coupled receptor (GPCR) and protein kinase families.

    View details for Web of Science ID 000292772100008

    View details for PubMedID 21521153

  • Identification of recurring protein structure microenvironments and discovery of novel functional sites around CYS residues BMC STRUCTURAL BIOLOGY Wu, S., Liu, T., Altman, R. B. 2010; 10


    The emergence of structural genomics presents significant challenges in the annotation of biologically uncharacterized proteins. Unfortunately, our ability to analyze these proteins is restricted by the limited catalog of known molecular functions and their associated 3D motifs.In order to identify novel 3D motifs that may be associated with molecular functions, we employ an unsupervised, two-phase clustering approach that combines k-means and hierarchical clustering with knowledge-informed cluster selection and annotation methods. We applied the approach to approximately 20,000 cysteine-based protein microenvironments (3D regions 7.5 A in radius) and identified 70 interesting clusters, some of which represent known motifs (e.g. metal binding and phosphatase activity), and some of which are novel, including several zinc binding sites. Detailed annotation results are available online for all 70 clusters at use of microenvironments instead of backbone geometric criteria enables flexible exploration of protein function space, and detection of recurring motifs that are discontinuous in sequence and diverse in structure. Clustering microenvironments may thus help to functionally characterize novel proteins and better understand the protein structure-function relationship.

    View details for DOI 10.1186/1472-6807-10-4

    View details for Web of Science ID 000275410900001

    View details for PubMedID 20122268

    View details for PubMedCentralID PMC2833161

  • Prediction of calcium-binding sites by combining loop-modeling with machine learning BMC STRUCTURAL BIOLOGY Liu, T., Altman, R. B. 2009; 9


    Protein ligand-binding sites in the apo state exhibit structural flexibility. This flexibility often frustrates methods for structure-based recognition of these sites because it leads to the absence of electron density for these critical regions, particularly when they are in surface loops. Methods for recognizing functional sites in these missing loops would be useful for recovering additional functional information.We report a hybrid approach for recognizing calcium-binding sites in disordered regions. Our approach combines loop modeling with a machine learning method (FEATURE) for structure-based site recognition. For validation, we compared the performance of our method on known calcium-binding sites for which there are both holo and apo structures. When loops in the apo structures are rebuilt using modeling methods, FEATURE identifies 14 out of 20 crystallographically proven calcium-binding sites. It only recognizes 7 out of 20 calcium-binding sites in the initial apo crystal structures.We applied our method to unstructured loops in proteins from SCOP families known to bind calcium in order to discover potential cryptic calcium binding sites. We built 2745 missing loops and evaluated them for potential calcium binding. We made 102 predictions of calcium-binding sites. Ten predictions are consistent with independent experimental verifications. We found indirect experimental evidence for 14 other predictions. The remaining 78 predictions are novel predictions, some with intriguing potential biological significance. In particular, we see an enrichment of beta-sheet folds with predicted calcium binding sites in the connecting loops on the surface that may be important for calcium-mediated function switches.Protein crystal structures are a potentially rich source of functional information. When loops are missing in these structures, we may be losing important information about binding sites and active sites. We have shown that limited loop modeling (e.g. loops less than 17 residues) combined with pattern matching algorithms can recover functions and propose putative conformations associated with these functions.

    View details for DOI 10.1186/1472-6807-9-72

    View details for Web of Science ID 000273849100001

    View details for PubMedID 20003365

    View details for PubMedCentralID PMC2808310

  • A novel method for predicting and using distance constraints of high accuracy for refining protein structure prediction PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS Liu, T., Horst, J. A., Samudrala, R. 2009; 77 (1): 220-234


    The principal bottleneck in protein structure prediction is the refinement of models from lower accuracies to the resolution observed by experiment. We developed a novel constraints-based refinement method that identifies a high number of accurate input constraints from initial models and rebuilds them using restrained torsion angle dynamics (rTAD). We previously created a Bayesian statistics-based residue-specific all-atom probability discriminatory function (RAPDF) to discriminate native-like models by measuring the probability of accuracy for atom type distances within a given model. Here, we exploit RAPDF to score (i.e., filter) constraints from initial predictions that may or may not be close to a native-like state, obtain consensus of top scoring constraints amongst five initial models, and compile sets with no redundant residue pair constraints. We find that this method consistently produces a large and highly accurate set of distance constraints from which to build refinement models. We further optimize the balance between accuracy and coverage of constraints by producing multiple structure sets using different constraint distance cutoffs, and note that the cutoff governs spatially near versus distant effects in model generation. This complete procedure of deriving distance constraints for rTAD simulations improves the quality of initial predictions significantly in all cases evaluated by us. Our procedure represents a significant step in solving the protein structure prediction and refinement problem, by enabling the use of consensus constraints, RAPDF, and rTAD for protein structure modeling and refinement.

    View details for DOI 10.1002/prot.22434

    View details for Web of Science ID 000269300000020

    View details for PubMedID 19422061

    View details for PubMedCentralID PMC2874729

  • Predicting drug side-effects by chemical systems biology GENOME BIOLOGY Tatonetti, N. P., Liu, T., Altman, R. B. 2009; 10 (9)


    New approaches to predicting ligand similarity and protein interactions can explain unexpected observations of drug inefficacy or side-effects.

    View details for DOI 10.1186/gb-2009-10-9-238

    View details for Web of Science ID 000271425300004

    View details for PubMedID 19723347

    View details for PubMedCentralID PMC2768971

  • Improving the accuracy of template-based predictions by mixing and matching between initial models BMC STRUCTURAL BIOLOGY Liu, T., Guerquin, M., Samudrala, R. 2008; 8


    Comparative modeling is a technique to predict the three dimensional structure of a given protein sequence based primarily on its alignment to one or more proteins with experimentally determined structures. A major bottleneck of current comparative modeling methods is the lack of methods to accurately refine a starting initial model so that it approaches the resolution of the corresponding experimental structure. We investigate the effectiveness of a graph-theoretic clique finding approach to solve this problem.Our method takes into account the information presented in multiple templates/alignments at the three-dimensional level by mixing and matching regions between different initial comparative models. This method enables us to obtain an optimized conformation ensemble representing the best combination of secondary structures, resulting in the refined models of higher quality. In addition, the process of mixing and matching accumulates near-native conformations, resulting in discriminating the native-like conformation in a more effective manner. In the seventh Critical Assessment of Structure Prediction (CASP7) experiment, the refined models produced are more accurate than the starting initial models.This novel approach can be applied without any manual intervention to improve the quality of comparative predictions where multiple template/alignment combinations are available for modeling, producing conformational models of higher quality than the starting initial predictions.

    View details for DOI 10.1186/1472-6807-8-24

    View details for Web of Science ID 000256681000001

    View details for PubMedID 18457597

  • Scoring functions for de novo protein structure prediction revisited. Methods in molecular biology (Clifton, N.J.) Ngan, S., Hung, L., Liu, T., Samudrala, R. 2008; 413: 243-281


    De novo protein structure prediction methods attempt to predict tertiary structures from sequences based on general principles that govern protein folding energetics and/or statistical tendencies of conformational features that native structures acquire, without the use of explicit templates. A general paradigm for de novo prediction involves sampling the conformational space, guided by scoring functions and other sequence-dependent biases, such that a large set of candidate ("decoy") structures are generated, and then selecting native-like conformations from those decoys using scoring functions as well as conformer clustering. High-resolution refinement is sometimes used as a final step to fine-tune native-like structures. There are two major classes of scoring functions. Physics-based functions are based on mathematical models describing aspects of the known physics of molecular interaction. Knowledge-based functions are formed with statistical models capturing aspects of the properties of native protein conformations. We discuss the implementation and use of some of the scoring functions from these two classes for de novo structure prediction in this chapter.

    View details for PubMedID 18075169

  • Scaffold proteins and the regeneration of visual pigments PHOTOCHEMISTRY AND PHOTOBIOLOGY Nawrot, M., Liu, T., Garwin, G. G., Crabb, J. W., Saari, J. C. 2006; 82 (6): 1482-1488


    CRALBP, cellular retinaldehyde-binding protein, is a retinoid-binding protein necessary for efficient regeneration of rod and cone visual pigments. The C terminus of CRALBP binds to the PDZ domains of EBP50/NHERF-1, which in turn bind to ezrin and actin, proteins localized to the apical processes of the retinal pigment epithelium. In this study, we examined structural features associated with the interaction of the two proteins. The C-terminal amino-acid sequence of 11 orthologous CRALBPs is either ENTAL, ENTAF or EDTAL. Peptides ending in each of these sequences inhibited the interaction of CRALBP and EBP50/NHERF-1 with the use of an overlay assay. Molecular modeling showed that both NTAL and NTAF formed similar networks of H bonds with PDZ1 of EBP50/ NHERF-1, and the side chains of both C-terminal Leu and Phe fit into the peptide-binding groove of PDZ1x CRALBP.11-cis-retinal and EBP50/NHERF-1 migrated as single components when analyzed individually by gel filtration and as a complex when mixed together before gel filtration. Complex formation was abolished by preincubation of EBP50/NHERF-1 with peptide EVENTAL. The ligand absorption spectrum of the complex was identical with that of CRALBP x 11-cis-retinal, demonstrating that complex formation did not perturb the ligand-binding domain of CRALBP.

    View details for DOI 10.1562/2006-01-25-RA-784

    View details for Web of Science ID 000243214100015

    View details for PubMedID 16553463

  • The effect of experimental resolution on the performance of knowledge-based discriminatory functions for protein structure selection PROTEIN ENGINEERING DESIGN & SELECTION Liu, T., Samudrala, R. 2006; 19 (9): 431-437


    The key to an accurate method of protein structure prediction is the development of an effective discriminatory function. Knowledge-based discriminatory functions extract parameters from statistical analysis of experimentally determined protein structures. We assess how the quality of the protein structures used for compiling statistics affects the performance of a residue-specific all-atom probability discriminatory function (RAPDF). We find that the discriminatory power correlates with the quality of the structural dataset on which the RAPDF is parameterized in a statistically significant manner. The overrepresentation of unfavorable contacts in the low-resolution and NMR structures contributes to the major errors in the compilation of the conditional probabilities. Such errors weaken the discriminatory power of the function, especially when decoy conformations also contain considerable numbers of unfavorable contacts. This indicates that using high-resolution structural datasets after filtering out unfavorable contacts can improve the performance of knowledge-based discriminatory functions.

    View details for DOI 10.1093/protein/gzl027

    View details for Web of Science ID 000240544600005

    View details for PubMedID 16845128

  • CRALBP ligand and protein interactions 11th International Symposium on Retinal Degeneration Wu, Z. P., Bhattacharya, S. K., Jin, Z. Y., Bonilha, V. L., Liu, T. Y., Nawrot, M., Teller, D. C., Saari, J. C., Crabb, J. W. SPRINGER-VERLAG BERLIN. 2006: 477–483

    View details for Web of Science ID 000236276600066

    View details for PubMedID 17249612

  • Structural insights into the cellular retinaldehyde-binding protein (CRALBP) PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS Liu, T. Y., Jenwitheesuk, E., Teller, D. C., Samudrala, R. 2005; 61 (2): 412-422


    Cellular retinaldehyde-binding protein (CRALBP) is an essential protein in the human visual cycle without a known three-dimensional structure. Previous studies associate retinal pathologies to specific mutations in the CRALBP protein. Here we use homology modeling and molecular dynamics methods to investigate the structural mechanisms by which CRALBP functions in the visual cycle. We have constructed two conformations of CRALBP representing two states in the process of ligand association and dissociation. Notably, our homology models map the pathology-associated mutations either directly in or adjacent to the putative ligand-binding cavity. Furthermore, six novel residues have been identified to be crucial for the hinge movement of the lipid-exchange loop in CRALBP. We conclude that the binding and release of retinoid involve large conformational changes in the lipid-exchange loop at the entrance of the ligand-binding cavity.

    View details for DOI 10.1002/prot.20621

    View details for Web of Science ID 000232420800019

    View details for PubMedID 16121400

  • PROTINFO: new algorithms for enhanced protein structure predictions NUCLEIC ACIDS RESEARCH Hung, L. H., Ngan, S. C., Liu, T., Samudrala, R. 2005; 33: W77-W80


    We describe new algorithms and modules for protein structure prediction available as part of the PROTINFO web server. The modules, comparative and de novo modelling, have significantly improved back-end algorithms that were rigorously evaluated at the sixth meeting on the Critical Assessment of Protein Structure Prediction methods. We were one of four server groups invited to make an oral presentation (only the best performing groups are asked to do so). These two modules allow a user to submit a protein sequence and return atomic coordinates representing the tertiary structure of that protein. The PROTINFO server is available at

    View details for DOI 10.1093/nar/gki403

    View details for Web of Science ID 000230271400012

    View details for PubMedID 15980581

  • Identification of CRALBP ligand interactions by photoaffinity labeling, hydrogen/deuterium exchange, and structural modeling JOURNAL OF BIOLOGICAL CHEMISTRY Wu, Z. P., Hasan, A., Liu, T. Y., Teller, D. C., Crabb, J. W. 2004; 279 (26): 27357-27364


    Cellular retinaldehyde-binding protein (CRALBP) functions in the retinal pigment epithelium (RPE) as an acceptor of 11-cis-retinol in the isomerization step of the rod visual cycle and as a substrate carrier for 11-cis-retinol dehydrogenase. Toward a better understanding of CRALBP function, the ligand binding cavity in human recombinant CRALBP (rCRALBP) was characterized by photoaffinity labeling with 3-diazo-4-keto-11-cis-retinal and by high resolution mass spectrometric topological analyses. Eight photoaffinity-modified residues were identified in rCRALBP by liquid chromatography tandem mass spectrometry, including Tyr(179), Phe(197), Cys(198), Met(208), Lys(221), Met(222), Val(223), and Met(225). Multiple different adduct masses were found on the photolabeled residues, and the molecular identity of each modification remains unknown. Supporting the specificity of photo-labeling, 50% of the modified residues have been associate with retinoid interactions by independent analyses. In addition, topological analysis of apo- and holo-rCRALBP by hydrogen/deuterium exchange and mass spectrometry demonstrated residues 198-255 incorporate significantly less deuterium when the retinoid binding pocket is occupied with 11-cis-retinal. This hydrophobic region encompasses all but one of the photo-labeled residues. A structural model of CRALBP ligand binding domain was constructed based on the crystal structures of three homologues in the CRAL-TRIO family of lipid-binding proteins. In the model, all of the photolabeled residues line the ligand binding cavity except Met(208), which appears to reside in a flexible loop at the entrance/exit of the ligand cavity. Overall, the results expand to 12 the number of residues proposed to interact with ligand and provide further insight into CRALBP ligand and protein interactions.

    View details for DOI 10.1074/jbc.M401960200

    View details for Web of Science ID 000222120400066

    View details for PubMedID 15100222