1
|
Pulido N, Guevara-Morales JM, Rodriguez-López A, Pulido Á, Díaz J, Edrada-Ebel RA, Echeverri-Peña OY. 1H-Nuclear Magnetic Resonance Analysis of Urine as Diagnostic Tool for Organic Acidemias and Aminoacidopathies. Metabolites 2021; 11:891. [PMID: 34940649 PMCID: PMC8704601 DOI: 10.3390/metabo11120891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Revised: 12/14/2021] [Accepted: 12/15/2021] [Indexed: 12/03/2022] Open
Abstract
The utility of low-resolution 1H-NMR analysis for the identification of biomarkers provided evidence for rapid biochemical diagnoses of organic acidemia and aminoacidopathy. 1H-NMR, with a sensitivity expected for a field strength of 400 MHz at 64 scans was used to establish the metabolomic urine sample profiles of an infant population diagnosed with small molecule Inborn Errors of Metabolism (smIEM) compared to unaffected individuals. A qualitative differentiation of the 1H-NMR spectral profiles of urine samples obtained from individuals affected by different organic acidemias and aminoacidopathies was achieved in combination with GC-MS. The smIEM disorders investigated in this study included phenylalanine metabolism; isovaleric, propionic, 3-methylglutaconicm and glutaric type I acidemia; and deficiencies in medium chain acyl-coenzyme and holocarboxylase synthase. The observed metabolites were comparable and similar to those reported in the literature, as well as to those detected with higher-resolution NMR. In this study, diagnostic marker metabolites were identified for the smIEM disorders. In some cases, changes in metabolite profiles differentiated post-treatments and follow-ups while allowing for the establishment of different clinical states of a biochemical disorder. In addition, for the first time, a 1H-NMR-based biomarker profile was established for holocarboxylase synthase deficiency spectrum.
Collapse
Affiliation(s)
- Ninna Pulido
- San Ignacio University Hospital, Bogota 4665684, Colombia;
- Institute for the Study of Inborn Errors of Metabolism, School of Sciences, Pontificia Universidad Javeriana, Cra 7 # 43-82, Edificio 54 Lab 303A, Bogota 4665684, Colombia;
| | - Johana M. Guevara-Morales
- Institute for the Study of Inborn Errors of Metabolism, School of Sciences, Pontificia Universidad Javeriana, Cra 7 # 43-82, Edificio 54 Lab 303A, Bogota 4665684, Colombia;
| | - Alexander Rodriguez-López
- Molecular Biology and Immunology Department, Fundación Instituto de Inmunología de Colombia (FIDIC), Bogota 4665684, Colombia;
- Chemistry Department, School of Sciences, Pontificia Universidad Javeriana, Bogota 4665684, Colombia;
| | - Álvaro Pulido
- Electrical and Electronics Department, Universidad Nacional de Colombia, Bogota 4665684, Colombia;
| | - Jhon Díaz
- Chemistry Department, School of Sciences, Pontificia Universidad Javeriana, Bogota 4665684, Colombia;
| | - Ru Angelie Edrada-Ebel
- The Natural Products Metabolomics Group, Strathclyde Institute of Pharmacy and Biomedical Sciences, Faculty of Science, University of Strathclyde, The John Arbuthnott Building, 161 Cathedral Street, Glasgow G4 0RE, UK
| | - Olga Y. Echeverri-Peña
- Institute for the Study of Inborn Errors of Metabolism, School of Sciences, Pontificia Universidad Javeriana, Cra 7 # 43-82, Edificio 54 Lab 303A, Bogota 4665684, Colombia;
| |
Collapse
|
2
|
Ntontsi P, Ntzoumanika V, Loukides S, Benaki D, Gkikas E, Mikros E, Bakakos P. EBC metabolomics for asthma severity. J Breath Res 2020; 14:036007. [PMID: 32392552 DOI: 10.1088/1752-7163/ab9220] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Asthma is a heterogeneous disease with diverse severity and represents a considerable socio-economic burden. Exhaled Breath Condensate (EBC) is a biofluid directly obtained from the airway lining fluid non-invasively. We attempted to discriminate severe from mild-to-moderate asthma using EBC metabolomics based on both NMR and UHPLC-MS techniques. 36 patients were included in this study (15 patients with severe and 21 with mild-to-moderate asthma). EBC was collected and analyzed using both NMR and UHPLC-MS techniques for possible metabolites. Using PLS and oPLS analysis for the UHPLC-MS data, no metabolite was found to be sufficient for the discrimination of asthma severity. However, when another PLS-regression model was applied five metabolites were found to discriminate severe from mild-to-moderate asthma. Amino-acid lysine was the only metabolite that discriminated the two study groups using NMR data (p= 0.04, t-test with Welch's correction, AUC 0.66). EBC is an easily available biofluid which directly represents the lower airways but difficult-to-use for metabolomic analysis. Our study presents some encouraging findings for the discrimination of asthma severity subgroups using EBC metabolomics but more well-designed studies with a higher number of patients need to be conducted.
Collapse
Affiliation(s)
- P Ntontsi
- 2nd Respiratory Medicine Department, Attikon University Hospital, National and Kapodistrian University of Athens Medical School, Athens, Greece
| | | | | | | | | | | | | |
Collapse
|
3
|
Ding X, Wang X, Li Q, Yu L, Song Q, Gai J, Yang S. Metabolomics Studies on Cytoplasmic Male Sterility during Flower Bud Development in Soybean. Int J Mol Sci 2019; 20:E2869. [PMID: 31212804 PMCID: PMC6627938 DOI: 10.3390/ijms20122869] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 06/03/2019] [Accepted: 06/10/2019] [Indexed: 12/20/2022] Open
Abstract
Abnormal reactive oxygen species (ROS) may mediate cytoplasmic male sterility (CMS). To observe the effect of ROS on soybean CMS, metabolite content and antioxidant enzyme activity in the flower buds between soybean N8855-derived CMS line and its maintainer were compared. Of the 612 metabolites identified, a total of 74 metabolites were significantly differentiated in flower buds between CMS line and its maintainer. The differential metabolites involved 32 differential flavonoids, 13 differential phenolamides, and 1 differential oxidized glutathione (GSSG) belonging to a non-enzymatic ROS scavenging system. We observed lower levels of flavonoids and antioxidant enzyme activities in flower buds of the CMS line than in its maintainer. Our results suggest that deficiencies of enzymatic and non-enzymatic ROS scavenging systems in soybean CMS line cannot eliminate ROS in anthers effectively, excessive accumulation of ROS triggered programmed cell death and ultimately resulted in pollen abortion of soybean CMS line.
Collapse
Affiliation(s)
- Xianlong Ding
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| | - Xuan Wang
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| | - Qiang Li
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| | - Lifeng Yu
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| | - Qijian Song
- Soybean Genomics and Improvement Laboratory, Beltsville Agricultural Research Center, USDA-ARS, Beltsville, MD 20705, USA.
| | - Junyi Gai
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| | - Shouping Yang
- Soybean Research Institute, National Center for Soybean Improvement, Key Laboratory of Biology and Genetic Improvement of Soybean (General, Ministry of Agriculture), State Key Laboratory of Crop Genetics and Germplasm Enhancement, Jiangsu Collaborative Innovation Center for Modern Crop Production, College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China.
| |
Collapse
|
4
|
Taraji M, Haddad PR, Amos RIJ, Talebi M, Szucs R, Dolan JW, Pohl CA. Chemometric-assisted method development in hydrophilic interaction liquid chromatography: A review. Anal Chim Acta 2017; 1000:20-40. [PMID: 29289311 DOI: 10.1016/j.aca.2017.09.041] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2017] [Revised: 09/22/2017] [Accepted: 09/24/2017] [Indexed: 02/09/2023]
Abstract
With an enormous growth in the application of hydrophilic interaction liquid chromatography (HILIC), there has also been significant progress in HILIC method development. HILIC is a chromatographic method that utilises hydro-organic mobile phases with a high organic content, and a hydrophilic stationary phase. It has been applied predominantly in the determination of small polar compounds. Theoretical studies in computer-aided modelling tools, most importantly the predictive, quantitative structure retention relationship (QSRR) modelling methods, have attracted the attention of researchers and these approaches greatly assist the method development process. This review focuses on the application of computer-aided modelling tools in understanding the retention mechanism, the classification of HILIC stationary phases, prediction of retention times in HILIC systems, optimisation of chromatographic conditions, and description of the interaction effects of the chromatographic factors in HILIC separations. Additionally, what has been achieved in the potential application of QSRR methodology in combination with experimental design philosophy in the optimisation of chromatographic separation conditions in the HILIC method development process is communicated. Developing robust predictive QSRR models will undoubtedly facilitate more application of this chromatographic mode in a broader variety of research areas, significantly minimising cost and time of the experimental work.
Collapse
Affiliation(s)
- Maryam Taraji
- Australian Centre for Research on Separation Science (ACROSS), School of Physical Sciences-Chemistry, University of Tasmania, Private Bag 75, Hobart 7001, Australia
| | - Paul R Haddad
- Australian Centre for Research on Separation Science (ACROSS), School of Physical Sciences-Chemistry, University of Tasmania, Private Bag 75, Hobart 7001, Australia.
| | - Ruth I J Amos
- Australian Centre for Research on Separation Science (ACROSS), School of Physical Sciences-Chemistry, University of Tasmania, Private Bag 75, Hobart 7001, Australia
| | - Mohammad Talebi
- Australian Centre for Research on Separation Science (ACROSS), School of Physical Sciences-Chemistry, University of Tasmania, Private Bag 75, Hobart 7001, Australia
| | - Roman Szucs
- Pfizer Global Research and Development, CT13 9NJ, Sandwich, UK
| | - John W Dolan
- LC Resources, 1795 NW Wallace Rd., McMinnville, OR 97128, USA
| | | |
Collapse
|
5
|
Zhang L, Huang Y, Lian M, Fan Z, Tian Y, Wang Y, Kang H, Liu S, Liu S, Li T, Shan Z. Metabolic profiling of hepatitis B virus-related hepatocellular carcinoma with diverse differentiation grades. Oncol Lett 2017; 13:1204-1210. [PMID: 28454235 PMCID: PMC5403281 DOI: 10.3892/ol.2017.5596] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2016] [Accepted: 09/09/2016] [Indexed: 12/11/2022] Open
Abstract
The most effective diagnostic tool for the majority of hepatocellular carcinoma (HCC) patients is determining the differentiation grade of their tumors. However liver biopsies, which are currently the most effective way of determining tumor differentiation grade, have several limitations. The present study was designed to select serum characteristic metabolites that correlate with the differentiation grades of hepatitis B virus (HBV)-related HCC, and so could be used in the clinic as a non-invasive method of differentiating patients with different grades of HCC. A total of 58 patients with HBV-related HCC were included in the present study, and divided into three groups according to their tumor differentiation grade. A further 20 patients with HBV-related liver cirrhosis and 19 healthy volunteers were enrolled. Ultra-performance liquid chromatography-mass spectrometry was used to analyze endogenous metabolites. Multivariate statistical analysis was used to examine the data using MZmine 2.0 software. The 14 metabolites that were highly correlated with specific differentiation grades of HCC were then selected for additional study. Receiver operator characteristic curve analysis was used to evaluate their clinical value. In total, 5 metabolites were finally identified, including lysophosphatidylcholine (16:0), oleamide, monoglyceride (0:0/15:0/0:0), lysophosphatidylcholine (18:0) and lysophosphatidylcholine [22:5(7Z,10Z,13Z,16Z,19Z)]. All these metabolites exhibited an excellent ability to distinguish different types of HCC with various differentiation grades and the area under the curve of these metabolites was up to 0.942, showing promising clinical value.
Collapse
Affiliation(s)
- Lei Zhang
- Tianjin Key Laboratory of Artificial Cell, Artificial Cell Engineering Technology Research Center of Public Health Ministry, Tianjin Third Central Hospital, Tianjin 300170, P.R. China.,Chemical Engineering Institute, Tianjin University, Tianjin 300072, P.R. China
| | - Ya Huang
- Clinical Laboratory Department, Chongqing City First People's Hospital of Wanzhou, Chongqing 404040, P.R. China
| | - Mingjian Lian
- Clinical Laboratory Department, Third Central Clinical College, Tianjin Medical University, Tianjin 300070, P.R. China
| | - Zhijuan Fan
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Yaqiong Tian
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Yufan Wang
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Hua Kang
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Shuang Liu
- Clinical Laboratory Department, Third Central Clinical College, Tianjin Medical University, Tianjin 300070, P.R. China
| | - Shuye Liu
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Tong Li
- Tianjin Key Laboratory of Artificial Cell, Clinical Laboratory Department, Tianjin Third Central Hospital, Tianjin 300170, P.R. China
| | - Zhongqiang Shan
- Chemical Engineering Institute, Tianjin University, Tianjin 300072, P.R. China
| |
Collapse
|
6
|
Afshinnia F, Rajendiran TM, Karnovsky A, Soni T, Wang X, Xie D, Yang W, Shafi T, Weir MR, He J, Brecklin CS, Rhee EP, Schelling JR, Ojo A, Feldman H, Michailidis G, Pennathur S. Lipidomic Signature of Progression of Chronic Kidney Disease in the Chronic Renal Insufficiency Cohort. Kidney Int Rep 2016; 1:256-268. [PMID: 28451650 PMCID: PMC5402253 DOI: 10.1016/j.ekir.2016.08.007] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
Introduction Human studies report conflicting results on the predictive power of serum lipids on the progression of chronic kidney disease. We aimed to systematically identify the lipids that predict progression to end-stage kidney disease. Methods From the Chronic Renal Insufficiency Cohort, 79 patients with chronic kidney disease stages 2 to 3 who progressed to end-stage kidney disease over 6 years of follow-up were selected and frequency matched by age, sex, race, and diabetes with 121 nonprogressors with less than 25% decline in estimated glomerular filtration rate during the follow-up. The patients were randomly divided into training and test sets. We applied liquid chromatography-mass spectrometry-based lipidomics on visit year 1 samples. Results We identified 510 lipids, of which the top 10 coincided with false discovery threshold of 0.058 in the training set. From the top 10 lipids, the abundance of diacylglycerols and cholesteryl esters was lower, but that of phosphatidic acid 44:4 and monoacylglycerol 16:0 was significantly higher in progressors. Using logistic regression models, a multimarker panel consisting of diacylglycerols and monoacylglycerol independently predicted progression. The c-statistic of the multimarker panel added to the base model consisting of estimated glomerular filtration rate and urine protein-to-creatinine ratio as compared with that of the base model was 0.92 (95% confidence interval: 0.88–0.97) and 0.83 (95% confidence interval: 0.76–0.90, P < 0.01), respectively, an observation that was validated in the test subset. Discussion We conclude that a distinct panel of lipids may improve prediction of progression of chronic kidney disease beyond estimated glomerular filtration rate and urine protein-to-creatinine ratio when added to the base model.
Collapse
Affiliation(s)
- Farsad Afshinnia
- Division of Nephrology Department of Internal Medicine, University of Michigan
| | - Thekkelnaycke M Rajendiran
- Department of Pathology, University of Michigan.,Michigan Regional Comprehensive Metabolomics Resource Core
| | - Alla Karnovsky
- Department of Computational Medicine & Bioinformatics, University of Michigan.,Michigan Regional Comprehensive Metabolomics Resource Core
| | - Tanu Soni
- Michigan Regional Comprehensive Metabolomics Resource Core
| | - Xue Wang
- Department of Biostatistics and Epidemiology, University of Pennsylvania
| | - Dawei Xie
- Department of Biostatistics and Epidemiology, University of Pennsylvania
| | - Wei Yang
- Department of Biostatistics and Epidemiology, University of Pennsylvania
| | | | - Matthew R Weir
- Division of Nephrology, Department of Medicine, University of Maryland School of Medicine
| | - Jiang He
- Tulane University School of Medicine
| | | | | | | | - Akinlolu Ojo
- Division of Nephrology Department of Internal Medicine, University of Michigan
| | - Harold Feldman
- Department of Biostatistics and Epidemiology, University of Pennsylvania
| | | | - Subramaniam Pennathur
- Division of Nephrology Department of Internal Medicine, University of Michigan.,Department of Computational Medicine & Bioinformatics, University of Michigan.,Michigan Regional Comprehensive Metabolomics Resource Core
| |
Collapse
|
7
|
Talebi M, Schuster G, Shellie RA, Szucs R, Haddad PR. Performance comparison of partial least squares-related variable selection methods for quantitative structure retention relationships modelling of retention times in reversed-phase liquid chromatography. J Chromatogr A 2015; 1424:69-76. [DOI: 10.1016/j.chroma.2015.10.099] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2015] [Revised: 10/27/2015] [Accepted: 10/29/2015] [Indexed: 11/27/2022]
|
8
|
Lopes TI, Geloneze B, Pareja JC, Calixto AR, Ferreira MM, Marsaioli AJ. Blood Metabolome Changes Before and After Bariatric Surgery: A1H NMR-Based Clinical Investigation. OMICS-A JOURNAL OF INTEGRATIVE BIOLOGY 2015; 19:318-27. [DOI: 10.1089/omi.2015.0009] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Thiago I.B. Lopes
- Chemistry Institute, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - Bruno Geloneze
- Laboratory of Investigation on Metabolism and Diabetes, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - José C. Pareja
- Laboratory of Investigation on Metabolism and Diabetes, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - Antonio R. Calixto
- Laboratory of Investigation on Metabolism and Diabetes, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - Márcia M.C. Ferreira
- Chemistry Institute, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| | - Anita J. Marsaioli
- Chemistry Institute, Universidade Estadual de Campinas (UNICAMP), Campinas, São Paulo, Brazil
| |
Collapse
|
9
|
|
10
|
Pramanik S, Roy K. Predictive modeling of chemical toxicity towards Pseudokirchneriella subcapitata using regression and classification based approaches. ECOTOXICOLOGY AND ENVIRONMENTAL SAFETY 2014; 101:184-190. [PMID: 24507144 DOI: 10.1016/j.ecoenv.2013.12.030] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2013] [Revised: 12/28/2013] [Accepted: 12/30/2013] [Indexed: 06/03/2023]
Abstract
Biodiversity nurturing may be a valuable pathway in controlling chemical stress on the ecosystem. In the present work, in silico studies have been performed to develop regression based quantitative structure toxicity relationship (QSTR) models using a data set containing 105 organic chemicals for the prediction of 48-h chemical toxicity towards Pseudokirchneriella subcapitata. Classification based linear discriminant analysis (LDA) was also performed to distinguish chemicals into toxic and nontoxic groups using the same data set. The developed models were found to possess good predictive quality in terms of internal, external and overall validation parameters. The regression based QSTR model suggests that second order molecular connectivity index (molecular size and lipophilicity), density (aromaticity), relative shape of molecules (cyclicity/aromaticity), and specific molecular fragments of the chemicals are important properties of chemicals to exert their toxicity on P. subcapitata. The classification based LDA QSTR model suggested that fused ring aromatic systems, secondary carbon atom fragments, second order valence molecular connectivity indices (molecular size and branching) and molecular weight are the distinguishing features to differentiate chemicals into toxic and nontoxic groups.
Collapse
Affiliation(s)
- Subrata Pramanik
- Drug Theoretics and Cheminformatics Laboratory, Division of Medicinal and Pharmaceutical Chemistry, Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700032, India
| | - Kunal Roy
- Drug Theoretics and Cheminformatics Laboratory, Division of Medicinal and Pharmaceutical Chemistry, Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700032, India.
| |
Collapse
|
11
|
Pramanik S, Roy K. Environmental toxicological fate prediction of diverse organic chemicals based on steady-state compartmental chemical mass ratio using quantitative structure-fate relationship (QSFR) models. CHEMOSPHERE 2013; 92:600-607. [PMID: 23642702 DOI: 10.1016/j.chemosphere.2013.03.065] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2013] [Revised: 03/06/2013] [Accepted: 03/22/2013] [Indexed: 06/02/2023]
Abstract
Four quantitative prediction models for steady-state compartmental chemical mass concentrations (Wn,g) were obtained from structural information, physiochemical properties, degradation rate and transport coefficients of 455 diverse organic chemicals using chemometric tools in a quantitative structure-fate relationship (QSFR) study. The mass ratio assessment of environmentally prevalent organic chemicals may be helpful to predict their toxicological fate in the ecosystems. Four sets of mass ratios [(1) log(Wair) from water emissions (water to air compartment), (2) log(Wair) from air emissions (within different zones of the air compartment), (3) log(Wwater) from water emissions (within different zones of the water compartment) and (4) log(Wwater) from air emissions (air to water compartment)] have been used. The developed models using genetic function approximation followed by multiple linear regression (GFA-MLR) and subsequent partial least squares (PLS) treatment identify only four descriptors for log(Wair) from water emission, six descriptors for log(Wair) from air emission, five descriptors for log(Wwater) from water emission and seven descriptors for log(Wwater) from air emission for predicting efficiently a large number of test set chemicals (ntest=182). The conclusive models suggest that descriptors such as partition coefficients (Kaw, Kow and Ksw), degradation parameters (Ksoil,Kwater and Kair), vapor pressure (Pv), diffusivity (Dwater), spatial descriptors (Jurs-WNSA-1, Jurs-WNSA-2, Jurs-WPSA-3, Jurs-FNSA-3 and Density), thermodynamic descriptors (MolRef and AlogP98), electrotopological state indices (S_dsN, S_ssNH and S_dsCH) are important for predicting the chemical mass ratios. The developed models may be applicable in toxicological fate prediction of diverse chemicals in the ecosystems.
Collapse
Affiliation(s)
- Subrata Pramanik
- Drug Theoretics and Cheminformatics Laboratory, Division of Medicinal and Pharmaceutical Chemistry, Department of Pharmaceutical Technology, Jadavpur University, Kolkata 700 032, India
| | | |
Collapse
|
12
|
Data processing pipelines for comprehensive profiling of proteomics samples by label-free LC–MS for biomarker discovery. Talanta 2011; 83:1209-24. [DOI: 10.1016/j.talanta.2010.10.029] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2010] [Revised: 10/18/2010] [Accepted: 10/21/2010] [Indexed: 01/30/2023]
|
13
|
Jing M, Cai W, Shao X. Quantitative Determination of the Components in Corn and Tobacco Samples by Using Near-Infrared Spectroscopy and Multiblock Partial Least Squares. ANAL LETT 2010. [DOI: 10.1080/00032711003686973] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
14
|
Horvatovich PL, Bischoff R. Current technological challenges in biomarker discovery and validation. EUROPEAN JOURNAL OF MASS SPECTROMETRY (CHICHESTER, ENGLAND) 2010; 16:101-121. [PMID: 20065518 DOI: 10.1255/ejms.1050] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
In this review we will give an overview of the issues related to biomarker discovery studies with a focus on liquid chromatography-mass spectrometry (LC-MS) methods. Biomarker discovery is based on a close collaboration between clinicians, analytical scientists and chemometritians/statisticians. It is critical to define the final purpose of a biomarker or biomarker pattern at the onset of the study and to select case and control samples accordingly. This is followed by designing the experiment, starting with the sampling strategy, sample collection, storage and separation protocols, choice and validation of the quantitative profiling platform followed by data processing, statistical analysis and validation workflows. Biomarker candidates that result after statistical validation should be submitted for further validation and, ideally, be connected to the disease mechanism after their identification. Since most discovery studies work with a relatively small number of samples, it is necessary to assess the specificity and sensitivity of a given biomarker-based assay in a larger set of independent samples, preferably analyzed at another clinical center. Targeted analytical methods of higher throughput than the original discovery method are needed at this point and LC-tandem mass spectrometry is gaining acceptance in this field. Throughout this review, we will focus on possible sources of variance and how they can be assessed and reduced in order to avoid false positives and to reduce the number of false negatives in biomarker discovery research.
Collapse
Affiliation(s)
- Peter L Horvatovich
- Analytical Biochemistry, Department of Pharmacy, University of Groningen, A. Deusinglaan 1, 9713 AV Groningen, The Netherlands.
| | | |
Collapse
|
15
|
Meta-heuristics on quantitative structure-activity relationships: study on polychlorinated biphenyls. J Mol Model 2009; 16:377-86. [DOI: 10.1007/s00894-009-0540-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2009] [Accepted: 05/13/2009] [Indexed: 10/20/2022]
|
16
|
Guha R. On the interpretation and interpretability of quantitative structure–activity relationship models. J Comput Aided Mol Des 2008; 22:857-71. [DOI: 10.1007/s10822-008-9240-5] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2008] [Accepted: 08/14/2008] [Indexed: 01/28/2023]
|
17
|
Hansen M, Kilk K, Langel Ü. Predicting cell-penetrating peptides. Adv Drug Deliv Rev 2008; 60:572-9. [PMID: 18045726 DOI: 10.1016/j.addr.2007.09.003] [Citation(s) in RCA: 108] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2007] [Accepted: 09/24/2007] [Indexed: 10/22/2022]
Abstract
Possibility to predict short peptide sequences capable to penetrate the plasma membrane opens new opportunities for developing peptide based intracellular delivery vectors, called cell-penetrating peptides (CPPs). Predictions of CPPs, however are often based on trial and error and may not always lead to new potent sequences. In this review we discuss different problems associated with CPP prediction. Additionally, the used methods of CPP prediction are compared. Also, a few suggestions are made for designing new CPP sequences and improvement of predictions.
Collapse
|
18
|
Larsen SB, Jørgensen FS, Olsen L. QSAR models for the human H(+)/peptide symporter, hPEPT1: affinity prediction using alignment-independent descriptors. J Chem Inf Model 2007; 48:233-41. [PMID: 18092768 DOI: 10.1021/ci700346y] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A data set comprising the major known chemical classes of hPEPT1 ligands was compiled from the literature. For these compounds, alignment-independent descriptors (VolSurf, GRIND/Almond, and MOE) were computed. Using hierarchical partial least-squares projection to latent structures (H-PLS), a one-component model with r2 = 0.77 and q2 = 0.75 was obtained. The model satisfied a set of rigorous validation criteria and performed well in the prediction of an external test set. Mechanistic interpretation of the model reveals polarity properties to be the dominant factors in determining hPEPT1 affinity, with hydrophobic interactions contributing to a lesser extent. The model is superior to previously reported models due to its combination of quality and speed. Accordingly, it is suitable for ligand-based virtual screening, such as QSAR-based database mining.
Collapse
Affiliation(s)
- Simon Birksø Larsen
- Biostructural Research, Department of Medicinal Chemistry, Faculty of Pharmaceutical Sciences, University of Copenhagen, 2 Universitetsparken, DK-2100 Copenhagen, Denmark
| | | | | |
Collapse
|
19
|
Dahlgren MK, Kauppi AM, Olsson IM, Linusson A, Elofsson M. Design, Synthesis, and Multivariate Quantitative Structure−Activity Relationship of SalicylanilidesPotent Inhibitors of Type III Secretion in Yersinia. J Med Chem 2007; 50:6177-88. [DOI: 10.1021/jm070741b] [Citation(s) in RCA: 54] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Anna M. Kauppi
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Ing-Marie Olsson
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Anna Linusson
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| | - Mikael Elofsson
- Department of Chemistry, Umeå University, SE-901 87 Umeå, Sweden
| |
Collapse
|
20
|
Kauppi AM, Andersson CD, Norberg HA, Sundin C, Linusson A, Elofsson M. Inhibitors of type III secretion in Yersinia: Design, synthesis and multivariate QSAR of 2-arylsulfonylamino-benzanilides. Bioorg Med Chem 2007; 15:6994-7011. [PMID: 17851084 DOI: 10.1016/j.bmc.2007.07.047] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2007] [Revised: 07/18/2007] [Accepted: 07/25/2007] [Indexed: 11/19/2022]
Abstract
Compound 1, 2-(benzo[1,2,5]thiadiazole-4-sulfonylamino)-5-chloro-N-(3,4-dichloro-phenyl)-benzamide, was identified as a putative type III secretion inhibitor in Yersinia, and the compound thus has a potential to be used to prevent or treat bacterial infections. A set of seven analogues was synthesized and evaluated in a type III secretion dependent reporter-gene assay with viable bacterial to give basic SAR. A second set of 19 compounds was obtained by statistical molecular design in the building block and product space and by subsequent synthesis. Evaluation in the reporter-gene assay showed that the compounds ranged from non-active to compounds more potent than 1. Based on the data multivariate QSAR models were established and the final Hi-PLS model showed good correlation between experimentally determined % inhibition and the calculated % inhibition of the reporter-gene signal.
Collapse
Affiliation(s)
- Anna M Kauppi
- Department of Chemistry, Umeå University, SE-90187 Umeå, Sweden
| | | | | | | | | | | |
Collapse
|
21
|
Gavaghan CL, Arnby CH, Blomberg N, Strandlund G, Boyer S. Development, interpretation and temporal evaluation of a global QSAR of hERG electrophysiology screening data. J Comput Aided Mol Des 2007; 21:189-206. [PMID: 17384921 DOI: 10.1007/s10822-006-9095-6] [Citation(s) in RCA: 64] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2006] [Accepted: 12/02/2006] [Indexed: 10/23/2022]
Abstract
A 'global' model of hERG K(+) channel was built to satisfy three basic criteria for QSAR models in drug discovery: (1) assessment of the applicability domain, (2) assuring that model decisions can be interpreted by medicinal chemists and (3) assessment of model performance after the model was built. A combination of D-optimal onion design and hierarchical partial least squares modelling was applied to construct a global model of hERG blockade in order to maximize the applicability domain of the model and to enhance its interpretability. Additionally, easily interpretable hERG specific fragment-based descriptors were developed. Model performance was monitored, throughout a time period of 15 months, after model implementation. It was found that after this time duration a greater proportion of molecules were outside the model's applicability domain and that these compounds had a markedly higher average prediction error than those from molecules within the model's applicability domain. The model's predictive performance deteriorated within 4 months after building, illustrating the necessity of regular updating of global models within a drug discovery environment.
Collapse
Affiliation(s)
- Claire L Gavaghan
- Computational Toxicology, Safety Assessment, AstraZeneca R&D, Molndal, Sweden.
| | | | | | | | | |
Collapse
|
22
|
de Boves Harrington P. Statistical validation of classification and calibration models using bootstrapped Latin partitions. Trends Analyt Chem 2006. [DOI: 10.1016/j.trac.2006.10.010] [Citation(s) in RCA: 86] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
23
|
Weis-Müller BT, Modlich O, Drobinskaya I, Unay D, Huber R, Bojar H, Schipke JD, Feindt P, Gams E, Müller W, Goecke T, Sandmann W. Gene expression in acute Stanford type A dissection: a comparative microarray study. J Transl Med 2006; 4:29. [PMID: 16824202 PMCID: PMC1557406 DOI: 10.1186/1479-5876-4-29] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2006] [Accepted: 07/06/2006] [Indexed: 11/21/2022] Open
Abstract
Background We compared gene expression profiles in acutely dissected aorta with those in normal control aorta. Materials and methods Ascending aorta specimen from patients with an acute Stanford A-dissection were taken during surgery and compared with those from normal ascending aorta from multiorgan donors using the BD Atlas™ Human1.2 Array I, BD Atlas™ Human Cardiovascular Array and the Affymetrix HG-U133A GeneChip®. For analysis only genes with strong signals of more than 70 percent of the mean signal of all spots on the array were accepted as being expressed. Quantitative real-time polymerase chain reaction (RT-PCR) was used to confirm regulation of expression of a subset of 24 genes known to be involved in aortic structure and function. Results According to our definition expression profiling of aorta tissue specimens revealed an expression of 19.1% to 23.5% of the genes listed on the arrays. Of those 15.7% to 28.9% were differently expressed in dissected and control aorta specimens. Several genes that encode for extracellular matrix components such as collagen IV α2 and -α5, collagen VI α3, collagen XIV α1, collagen XVIII α1 and elastin were down-regulated in aortic dissection, whereas levels of matrix metalloproteinases-11, -14 and -19 were increased. Some genes coding for cell to cell adhesion, cell to matrix signaling (e.g., polycystin1 and -2), cytoskeleton, as well as several myofibrillar genes (e.g., α-actinin, tropomyosin, gelsolin) were found to be down-regulated. Not surprisingly, some genes associated with chronic inflammation such as interleukin -2, -6 and -8, were up-regulated in dissection. Conclusion Our results demonstrate the complexity of the dissecting process on a molecular level. Genes coding for the integrity and strength of the aortic wall were down-regulated whereas components of inflammatory response were up-regulated. Altered patterns of gene expression indicate a pre-existing structural failure, which is probably a consequence of insufficient remodeling of the aortic wall resulting in further aortic dissection.
Collapse
Affiliation(s)
- Barbara Theresia Weis-Müller
- Department of Vascular Surgery and Kidney Transplantation, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Olga Modlich
- Department of Chemical Oncology, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Irina Drobinskaya
- Department of Chemical Oncology, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Derya Unay
- Department of Chemical Oncology, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Rita Huber
- Department of Vascular Surgery and Kidney Transplantation, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Hans Bojar
- Department of Chemical Oncology, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Jochen D Schipke
- Research Group Experimental Surgery, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Peter Feindt
- Department of Thoracic and Cardiovascular Surgery, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Emmeran Gams
- Department of Thoracic and Cardiovascular Surgery, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Wolfram Müller
- Pathology Starnberg, private pathological practice, Starnberg, Germany
| | - Timm Goecke
- Institute of Human Genetics, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| | - Wilhelm Sandmann
- Department of Vascular Surgery and Kidney Transplantation, Heinrich-Heine-University of Düsseldorf, Düsseldorf, Germany
| |
Collapse
|
24
|
Eriksson L, Andersson PL, Johansson E, Tysklind M. Megavariate Analysis of Environmental QSAR Data. Part II – Investigating Very Complex Problem Formulations Using Hierarchical, Non-Linear and Batch-Wise Extensions of PCA and PLS. Mol Divers 2006; 10:187-205. [PMID: 16802062 DOI: 10.1007/s11030-006-9026-4] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2005] [Accepted: 02/02/2006] [Indexed: 10/24/2022]
Abstract
Three extensions of the basic PCA and PLS methodologies are described. These extensions are hierarchical, non-linear and batch-based in nature. The objectives of these methods are to assist in problem understanding and problem solving in very complex (QSAR) problem formulations. The method extensions are illustrated using two example QSAR data sets containing many X- and Y-variables.
Collapse
|
25
|
Chandrasekaran SM, Bhartiya S, Wangikar PP. Substrate specificity of lipases in alkoxycarbonylation reaction: QSAR model development and experimental validation. Biotechnol Bioeng 2006; 94:554-64. [PMID: 16528758 DOI: 10.1002/bit.20879] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Although lipases are known to catalyze alkoxycarbonylation reactions in organic solvents, the existing knowledge base on their substrate specificity in alkoxycarbonylation reaction is sparse. Moreover, models to predict substrate specificity have not been reported. Here, we report the experimentally measured rate constants for 180 acyl donor-alcohol pairs and demonstrate the two-step synthesis of over 70 disubstituted carbonate products from simple precursors such as diphenyl carbonate and alcohols. The efficiency of synthesis was found to be dependent on the order of alcohol addition. This motivated the need to develop a model to predict lipase specificity in alkoxycarbonylation reactions. A partial least square model has been developed to correlate the reaction rate with (i) descriptors of alcohol for a fixed acyl donor, (ii) descriptors of acyl donor for a fixed alcohol, (iii) descriptors of both the acyl donor and the alcohol. The number of descriptors being far greater than the number of observations was a potential limitation in the model development. This was addressed by selecting a subset of descriptors using a systematic procedure based on (a) correlation among the descriptors and step-wise regression methodology, and (b) variable influence on projection methodology. The model was able to accurately predict the reaction rate and the optimal order of addition of alcohols in the two-step synthesis of disubstituted carbonates using the enzyme mixture. The descriptor subset and the relevant model would benefit the users of lipases in synthetic applications while the modeling strategy presented here can have applications in predicting specificity of other enzymes.
Collapse
Affiliation(s)
- Sangeetha M Chandrasekaran
- Department of Chemical Engineering, Indian Institute of Technology, Bombay, Powai, Mumbai 400 076, India
| | | | | |
Collapse
|
26
|
Eriksson L, Andersson PL, Johansson E, Tysklind M. Megavariate analysis of environmental QSAR data. Part I – A basic framework founded on principal component analysis (PCA), partial least squares (PLS), and statistical molecular design (SMD). Mol Divers 2006; 10:169-86. [PMID: 16770514 DOI: 10.1007/s11030-006-9024-6] [Citation(s) in RCA: 99] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2005] [Accepted: 02/02/2006] [Indexed: 10/24/2022]
Abstract
This paper introduces principal component analysis (PCA), partial least squares projections to latent structures (PLS), and statistical molecular design (SMD) as useful tools in deriving multi- and megavariate quantitative structure-activity relationship (QSAR) models. Two QSAR data sets from the fields of environmental toxicology and environmental chemistry are worked out in detail, showing the benefits of PCA, PLS and SMD. PCA is useful when overviewing a data set and exploring relationships among compounds and relationships among variables. PLS is the regression extension of PCA and is used for establishing QSARs. SMD is essential for selecting informative training and test sets of compounds for QSAR calibration and validation.
Collapse
|
27
|
Lindström A, Pettersson F, Almqvist F, Berglund A, Kihlberg J, Linusson A. Hierarchical PLS Modeling for Predicting the Binding of a Comprehensive Set of Structurally Diverse Protein−Ligand Complexes. J Chem Inf Model 2006; 46:1154-67. [PMID: 16711735 DOI: 10.1021/ci050323k] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A new approach is presented for predicting ligand binding to proteins using hierarchical partial-least-squares regression to latent structures (Hi-PLS). Models were based on information from the 2002 release of the PDBbind database containing (after in-house refinement) high-resolution X-ray crystallography and binding affinity (Kd or Ki) data for 612 protein-ligand complexes. The complexes were characterized by four different descriptor blocks: three-dimensional (3D) structural descriptors of the proteins, protein-ligand interactions according to the Validate scoring function, binding site surface areas, and ligand 2D and 3D descriptors. These descriptor blocks were used in Hi-PLS models, generated using both linear and nonlinear terms, to relate the characterizations to pKd/i. The results show that each of the four descriptor blocks contributed to the model, and the predictions of pKd/i of the internal test set gave a root-mean-square error of prediction (RMSEP) of 1.65. The data were further divided according to the structural classification of the proteins, and Hi-PLS models were constructed for the resulting subclasses. The models for the four subclasses differed considerably in terms of both their ability to predict pKd/i (with RMSEPs ranging from 0.8 to 1.56) and the descriptor block that had the strongest influence. The models were validated with an external test set of 174 complexes from the 2003 release of the PDBbind database. The overall results show that the presented Hi-PLS methodology could facilitate the difficult task of predicting binding affinity.
Collapse
Affiliation(s)
- Anton Lindström
- Organic Chemistry, Department of Chemistry, Umeå University, SE-901 87, Umeå, Sweden
| | | | | | | | | | | |
Collapse
|
28
|
|
29
|
Kriegl JM, Eriksson L, Arnhold T, Beck B, Johansson E, Fox T. Multivariate modeling of cytochrome P450 3A4 inhibition. Eur J Pharm Sci 2005; 24:451-63. [PMID: 15784335 DOI: 10.1016/j.ejps.2004.12.009] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2004] [Revised: 12/08/2004] [Accepted: 12/28/2004] [Indexed: 11/18/2022]
Abstract
In the early phases of current pharmaceutical research projects, huge numbers of compounds are tested on their biological activity with respect to a certain target by experimental or virtual screening campaigns. To reduce the attrition rate in later stages of a project, other relevant properties such as physicochemical and ADMET (absorption, distribution, metabolism, excretion, toxicity) properties should be assessed as early as possible in lead discovery and optimization. The present study describes the development of in silico models to predict the inhibition of human cytochrome P450 3A4 (CYP3A4) from calculated molecular descriptors. The models were trained and validated using a set of 967 structural diverse drug-like research compounds with an experimentally determined CYP3A4 inhibition potency (IC50 value) which was carefully split into a training and a test set. For classification models, the data sets were further subdivided into strong, medium, and weak inhibitors. Different descriptor sets were used to cover various aspects of molecular properties, including properties derived from the 2D structure, the interaction of the molecule with its environment, and properties derived from quantum-mechanical calculations. The descriptors were related to the CYP3A4 inhibition potency by multivariate data analysis methods such as partial least-squares projection to latent structures (PLS), PLS discriminant analysis (PLS-DA), and soft independent class modeling (SIMCA). The squared correlation between experimental and predicted IC50 values of the previously unseen test set compounds was Qext2=0.6 for the best PLS models, corresponding to a root mean squared error (RMSE) of RMSE=0.45 (logarithm of IC50). The best PLS-DA models were able to correctly classify more than 60% of the test set compounds, whereas almost no strong inhibitors were wrongly classified as weak inhibitors and vice versa. Furthermore, relevant molecular properties were identified which are closely related to the CYP3A4 inhibition potency of a compound. The results presented here are very encouraging since our models could, for instance, serve to flag problematic compounds or to guide further synthesis efforts.
Collapse
Affiliation(s)
- Jan M Kriegl
- Computational Chemistry, Department of Lead Discovery, Boehringer Ingelheim Pharma GmbH & Co. KG, D-88397 Biberach, Germany.
| | | | | | | | | | | |
Collapse
|
30
|
Kriegl JM, Arnhold T, Beck B, Fox T. A support vector machine approach to classify human cytochrome P450 3A4 inhibitors. J Comput Aided Mol Des 2005; 19:189-201. [PMID: 16059671 DOI: 10.1007/s10822-005-3785-3] [Citation(s) in RCA: 51] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2004] [Accepted: 03/14/2005] [Indexed: 02/06/2023]
Abstract
The cytochrome P450 (CYP) enzyme superfamily plays a major role in the metabolism of commercially available drugs. Inhibition of these enzymes by a drug may result in a plasma level increase of another drug, thus leading to unwanted drug-drug interactions when two or more drugs are coadministered. Therefore, fast and reliable in silico methods predicting CYP inhibition from calculated molecular properties are an important tool which can be applied to assess both already synthesized as well as virtual compounds. We have studied the performance of support vector machines (SVMs) to classify compounds according to their potency to inhibit CYP3A4. The data set for model generation consists of more than 1300 structural diverse drug-like research molecules which were divided into training and test sets. The predictive power of SVMs crucially depends on a careful selection of parameters specifying the kernel function and the penalty for misclassifications. In this study we have investigated a procedure to identify a valid set of SVM parameters which is based on a sampling of the parameter space on a regular grid. From this set of parameters, either single SVMs or SVM committees were trained to distinguish between strong and weak inhibitors or to achieve a more realistic three-class assignment, with one class representing medium inhibitors. This workflow was studied for several kernel functions and descriptor sets. All SVM models performed significantly better than PLS-DA models which were generated from the corresponding descriptor sets. As a very promising result, simple two-dimensional (2D) descriptors yield a three-class model which correctly classifies more than 70% of the test set. Our work illustrates that SVMs used in combination with simple 2D descriptors provide a very effective and reliable tool which allows a fast assessment of CYP3A4 inhibition potency in an early in silico filtering process.
Collapse
Affiliation(s)
- Jan M Kriegl
- Computational Chemistry, Department of Lead Discovery, Boehringer Ingelheim Pharma GmbH & Co. KG, D-88397, Biberach, Germany.
| | | | | | | |
Collapse
|
31
|
Modlich O, Prisack HB, Munnes M, Audretsch W, Bojar H. Immediate Gene Expression Changes After the First Course of Neoadjuvant Chemotherapy in Patients with Primary Breast Cancer Disease. Clin Cancer Res 2004; 10:6418-31. [PMID: 15475428 DOI: 10.1158/1078-0432.ccr-04-1031] [Citation(s) in RCA: 50] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
PURPOSE Our goal was to identify genes undergoing expressional changes shortly after the beginning of neoadjuvant chemotherapy for primary breast cancer. EXPERIMENTAL DESIGN The biopsies were taken from patients with primary breast cancer prior to any treatment and 24 hours after the beginning of the neoadjuvant chemotherapy. Expression analyses from matched pair samples representing 25 patients were carried out with Clontech filter arrays. A subcohort of those 25 paired samples were additionally analyzed with the Affymetrix GeneChip platform. All of the transcripts from both platforms were queried for expressional changes. RESULTS Performing hierarchical cluster analysis, we clustered pre- and posttreatment samples from individual patients more closely to each other than the samples taken from different patients. This reflects the rather low number of transcripts responding directly to the drugs used. Although transcriptional drug response occurring during therapy differed between individual patients, two genes (p21(WAF1/CIP1) and MIC-1) were up-regulated in posttreatment samples. This could be validated by semiquantitative and real-time reverse transcription-PCR. Partial least- discriminant analysis based on approximately 25 genes independently identified by either Clontech or Affymetrix platforms could clearly discriminate pre- and posttreatment samples. However, correlation of certain gene expression levels as well as of differential patterns and clusters as determined by a different platform was not always satisfying. CONCLUSIONS This study has demonstrated the potential of monitoring posttreatment changes in gene expression as a measure of the pharmacodynamics of drugs. As a clinical laboratory model, it can be useful to identify patients with sensitive and reactive tumors and to help for optimized choice for sequential therapy and obviously improve relapse- free and overall survival.
Collapse
Affiliation(s)
- Olga Modlich
- Institute of Chemical Oncology, University of Düsseldorf, Düsseldorf, Germany.
| | | | | | | | | |
Collapse
|
32
|
Eriksson L, Antti H, Gottfries J, Holmes E, Johansson E, Lindgren F, Long I, Lundstedt T, Trygg J, Wold S. Using chemometrics for navigating in the large data sets of genomics, proteomics, and metabonomics (gpm). Anal Bioanal Chem 2004; 380:419-29. [PMID: 15448969 DOI: 10.1007/s00216-004-2783-y] [Citation(s) in RCA: 191] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2004] [Accepted: 07/22/2004] [Indexed: 10/26/2022]
Abstract
This article describes the applicability of multivariate projection techniques, such as principal-component analysis (PCA) and partial least-squares (PLS) projections to latent structures, to the large-volume high-density data structures obtained within genomics, proteomics, and metabonomics. PCA and PLS, and their extensions, derive their usefulness from their ability to analyze data with many, noisy, collinear, and even incomplete variables in both X and Y. Three examples are used as illustrations: the first example is a genomics data set and involves modeling of microarray data of cell cycle-regulated genes in the microorganism Saccharomyces cerevisiae. The second example contains NMR-metabonomics data, measured on urine samples of male rats treated with either of the drugs chloroquine or amiodarone. The third and last data set describes sequence-function classification studies in a set of G-protein-coupled receptors using hierarchical PCA.
Collapse
|
33
|
Weaver DC. Applying data mining techniques to library design, lead generation and lead optimization. Curr Opin Chem Biol 2004; 8:264-70. [PMID: 15183324 DOI: 10.1016/j.cbpa.2004.04.005] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Many data mining techniques have been applied to activity and ADMET datasets and the resulting models are being used to understand quantitative structure-activity relationships and design new libraries. This review summarizes data mining concepts and discuss their application to library design, lead generation (particularly for sequential screening) and lead optimization (specifically for generating and interpreting QSAR models). Also, this review discusses recent comparative studies between data mining techniques and draws some conclusions about the patterns emerging in the drug discovery data mining field.
Collapse
Affiliation(s)
- Daniel C Weaver
- Array Biopharma, Inc., 3200 Walnut Street, Boulder, Colorado 80303, USA.
| |
Collapse
|
34
|
Eriksson L, Jaworska J, Worth AP, Cronin MTD, McDowell RM, Gramatica P. Methods for reliability and uncertainty assessment and for applicability evaluations of classification- and regression-based QSARs. ENVIRONMENTAL HEALTH PERSPECTIVES 2003; 111:1361-75. [PMID: 12896860 PMCID: PMC1241620 DOI: 10.1289/ehp.5758] [Citation(s) in RCA: 832] [Impact Index Per Article: 39.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
This article provides an overview of methods for reliability assessment of quantitative structure-activity relationship (QSAR) models in the context of regulatory acceptance of human health and environmental QSARs. Useful diagnostic tools and data analytical approaches are highlighted and exemplified. Particular emphasis is given to the question of how to define the applicability borders of a QSAR and how to estimate parameter and prediction uncertainty. The article ends with a discussion regarding QSAR acceptability criteria. This discussion contains a list of recommended acceptability criteria, and we give reference values for important QSAR performance statistics. Finally, we emphasize that rigorous and independent validation of QSARs is an essential step toward their regulatory acceptance and implementation.
Collapse
|