A computational method for detection of ligand-binding proteins from dose range thermal proteome profiles.
Ontology highlight
ABSTRACT:
SUBMITTER: Kurzawa N
PROVIDER: S-EPMC7666118 | biostudies-literature | 2020 Nov
REPOSITORIES: biostudies-literature
Publications
A computational method for detection of ligand-binding proteins from dose range thermal proteome profiles.
Kurzawa Nils N Becher Isabelle I Sridharan Sindhuja S Franken Holger H Mateus André A Anders Simon S Bantscheff Marcus M Huber Wolfgang W Savitski Mikhail M MM
Nature communications 20201113 1
Detecting ligand-protein interactions in living cells is a fundamental challenge in molecular biology and drug research. Proteome-wide profiling of thermal stability as a function of ligand concentration promises to tackle this challenge. However, current data analysis strategies use preset thresholds that can lead to suboptimal sensitivity/specificity tradeoffs and limited comparability across datasets. Here, we present a method based on statistical hypothesis testing on curves, which provides ...[more]
Similar Datasets
Project description:Detecting the targets of drugs and other molecules in intact cellular contexts is a major objective in drug discovery and in biology more broadly. Thermal proteome profiling (TPP) pursues this aim at proteome-wide scale by inferring target engagement from its effects on temperature-dependent protein denaturation. However, a key challenge of TPP is the statistical analysis of the measured melting curves with controlled false discovery rates at high proteome coverage and detection power. We present nonparametric analysis of response curves (NPARC), a statistical method for TPP based on functional data analysis and nonlinear regression. We evaluate NPARC on five independent TPP data sets and observe that it is able to detect subtle changes in any region of the melting curves, reliably detects the known targets, and outperforms a melting point-centric, single-parameter fitting approach in terms of specificity and sensitivity. NPARC can be combined with established analysis of variance (ANOVA) statistics and enables flexible, factorial experimental designs and replication levels. An open source software implementation of NPARC is provided.
| S-EPMC6885700 | biostudies-literature
Project description:The CETSA and Thermal Proteome Profiling (TPP) analytical methods are invaluable for the study of protein-ligand interactions and protein stability in a cellular context. These tools have increasingly been leveraged in work ranging from understanding signaling paradigms to drug discovery. Consequently, there is an important need to optimize the data analysis pipeline that is used to calculate protein melt temperatures (Tm) and relative melt shifts from proteomics abundance data. Here, we report a user-friendly analysis of the melt shift calculation workflow where we describe the impact of each individual calculation step on the final output list of stabilized and destabilized proteins. This report also includes a description of how key steps in the analysis workflow quantitatively impact the list of stabilized/destabilized proteins from an experiment. We applied our findings to develop a more optimized analysis workflow that illustrates the dramatic sensitivity of chosen calculation steps on the final list of reported proteins of interest in a study and have made the R based program Inflect available for research community use through the CRAN repository [McCracken, N. Inflect: Melt Curve Fitting and Melt Shift Analysis. R package version 1.0.3, 2021]. The Inflect outputs include melt curves for each protein which passes filtering criteria in addition to a data matrix which is directly compatible with downstream packages such as UpsetR for replicate comparisons and identification of biologically relevant changes. Overall, this work provides an essential resource for scientists as they analyze data from TPP and CETSA experiments and implement their own analysis pipelines geared toward specific applications.
| S-EPMC8022325 | biostudies-literature
Project description:Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, the docking scores are not sufficiently precise to represent the protein-ligand binding affinity. Here, we developed an efficient computational method for calculating protein-ligand binding affinity, which is based on molecular mechanics generalized Born/surface area (MM-GBSA) calculations and Jarzynski identity. Jarzynski identity is an exact relation between free energy differences and the work done through non-equilibrium process, and MM-GBSA is a semimacroscopic approach to calculate the potential energy. To calculate the work distribution when a ligand is pulled out of its binding site, multiple protein-ligand conformations are randomly generated as an alternative to performing an explicit single-molecule pulling simulation. We assessed the new method, multiple random conformation/MM-GBSA (MRC-MMGBSA), by evaluating ligand-binding affinities (scores) for four target proteins, and comparing these scores with experimental data. The calculated scores were qualitatively in good agreement with the experimental binding affinities, and the optimal docking structure could be determined by ranking the scores of the multiple docking poses obtained by the molecular docking process. Furthermore, the scores showed a strong linear response to experimental binding free energies, so that the free energy difference of the ligand binding (ΔΔG) could be calculated by linear scaling of the scores. The error of calculated ΔΔG was within ≈ ± 1.5 kcal.mol(-1) of the experimental values. Particularly, in the case of flexible target proteins, the MRC-MMGBSA scores were more effective in ranking ligands than those generated by the MM-GBSA method using a single protein-ligand conformation. The results suggest that, owing to its lower computational costs and greater accuracy, the MRC-MMGBSA offers efficient means to rank the ligands, in the post-docking process, according to their binding affinities, and to compare these directly with the experimental values.
| S-EPMC3423425 | biostudies-literature
Project description:The study of low-abundance proteins is a challenge to discovery-based proteomics. Mass spectrometry (MS) applications, such as thermal proteome profiling (TPP), face specific challenges in the detection of the whole proteome as a consequence of the use of nondenaturing extraction buffers. TPP is a powerful method for the study of protein thermal stability, but quantitative accuracy is highly dependent on consistent detection. Therefore, TPP can be limited in its amenability to study low-abundance proteins that tend to have stochastic or poor detection by MS. To address this challenge, we incorporated an affinity-purified protein complex sample at submolar concentrations as an isobaric trigger channel into a mutant TPP (mTPP) workflow to provide reproducible detection and quantitation of the low-abundance subunits of the cleavage and polyadenylation factor (CPF) complex. The inclusion of an isobaric protein complex trigger channel increased detection an average of 40× for previously detected subunits and facilitated detection of CPF subunits that were previously below the limit of detection. Importantly, these gains in CPF detection did not cause large changes in melt temperature (Tm) calculations for other unrelated proteins in the samples, with a high positive correlation between Tm estimates in samples with and without isobaric trigger channel addition. Overall, the incorporation of an affinity-purified protein complex as an isobaric trigger channel within a tandem mass tag (TMT) multiplex for mTPP experiments is an effective and reproducible way to gather thermal profiling data on proteins that are not readily detected using the original TPP or mTPP protocols.
| S-EPMC8153406 | biostudies-literature
Project description:Proteins interact with small molecules through specific molecular recognition, which is central to essential biological functions in living systems. Therefore, understanding such interactions is crucial for basic sciences and drug discovery. Here, we present Structure template-based ab initio ligand design solution (Stalis), a knowledge-based approach that uses structure templates from the Protein Data Bank libraries of whole ligands and their fragments and generates a set of molecules (virtual ligands) whose structures represent the pocket shape and chemical features of a given target binding site. Our benchmark performance evaluation shows that ligand structure-based virtual screening using virtual ligands from Stalis outperforms a receptor structure-based virtual screening using AutoDock Vina, demonstrating reliable overall screening performance applicable to computational high-throughput screening. However, virtual ligands from Stalis are worse in recognizing active compounds at the small fraction of a rank-ordered list of screened library compounds than crystal ligands, due to the low resolution of the virtual ligand structures. In conclusion, Stalis can facilitate drug discovery research by designing virtual ligands that can be used for fast ligand structure-based virtual screening. Moreover, Stalis provides actual three-dimensional ligand structures that likely bind to a target protein, enabling to gain structural insight into potential ligands. Stalis can be an efficient computational platform for high-throughput ligand design for fundamental biological study and drug discovery research at the proteomic level. © 2019 Wiley Periodicals, Inc.
| S-EPMC6878116 | biostudies-literature
Project description:I-Motif is a tetrameric cytosine-rich DNA structure with hemi-protonated cytosine: cytosine base pairs. Recent evidence showed that i-motif structures in human cells play regulatory roles in the genome. Therefore, characterization of novel i-motifs and investigation of their functional implication are urgently needed for comprehensive understanding of their roles in gene regulation. However, considering the complications of experimental investigation of i-motifs and the large number of putative i-motifs in the genome, development of an in silico tool for the characterization of i-motifs in the high throughput scale is necessary. We developed a novel computation method, MD-TSPC4, to predict the thermal stability of i-motifs based on molecular modeling and molecular dynamic simulation. By assuming that the flexibility of loops in i-motifs correlated with thermal stability within certain temperature ranges, we evaluated the correlation between the root mean square deviations (RMSDs) of model structures and the thermal stability as the experimentally obtained melting temperature (Tm). Based on this correlation, we propose an equation for Tm prediction from RMSD. We expect this method can be useful for estimating the overall structure and stability of putative i-motifs in the genome, which can be a starting point of further structural and functional studies of i-motifs.
| S-EPMC7793491 | biostudies-literature
Project description:Solute binding proteins (SBPs) form a heterogeneous protein family that is found in all kingdoms of life. In bacteria, the ligand-loaded forms bind to transmembrane transporters providing the substrate. We present here the SBP repertoire of Pseudomonas aeruginosa PAO1 that is composed of 98 proteins. Bioinformatic predictions indicate that many of these proteins have a redundant ligand profile such as 27 SBPs for proteinogenic amino acids, 13 proteins for spermidine/putrescine, or 9 proteins for quaternary amines. To assess the precision of these bioinformatic predictions, we have purified 17 SBPs that were subsequently submitted to high-throughput ligand screening approaches followed by isothermal titration calorimetry studies, resulting in the identification of ligands for 15 of them. Experimentation revealed that PA0222 was specific for γ-aminobutyrate (GABA), DppA2 for tripeptides, DppA3 for dipeptides, CysP for thiosulphate, OpuCC for betaine, and AotJ for arginine. Furthermore, RbsB bound D-ribose and D-allose, ModA bound molybdate, tungstate, and chromate, whereas AatJ recognized aspartate and glutamate. The majority of experimentally identified ligands were found to be chemoattractants. Data show that the ligand class recognized by SPBs can be predicted with confidence using bioinformatic methods, but experimental work is necessary to identify the precise ligand profile.
| S-EPMC6829864 | biostudies-literature
Project description:Humans are exposed to thousands of chemicals, including environmental chemicals. Unfortunately, little is known about their potential toxicity, as determining the toxicity remains challenging due to the substantial resources required to assess a chemical in vivo. Here, we present a novel hybrid neural network (HNN) deep learning method, called HNN-Tox, to predict chemical toxicity at different doses. To develop a hybrid HNN-Tox method, we combined two neural network frameworks, the Convolutional Neural Network (CNN) and the multilayer perceptron (MLP)-type feed-forward neural network (FFNN). Combining the CNN and FCNN in the field of environmental chemical toxicity prediction is a novel approach. We developed several binary and multiclass classification models to assess dose-range chemical toxicity that is trained based on thousands of chemicals with known toxicity. The performance of the HNN-Tox was compared with other machine-learning methods, including Random Forest (RF), Bootstrap Aggregation (Bagging), and Adaptive Boosting (AdaBoost). We also analyzed the model performance dependency on varying features, descriptors, dataset size, route of exposure, and toxic dose. The HNN-Tox model, trained on 59,373 chemicals annotated with known LD50 and routes of exposure, maintained its predictive ability with an accuracy of 84.9% and 84.1%, even after reducing the descriptor size from 318 to 51, and the area under the ROC curve (AUC) was 0.89 and 0.88, respectively. Further, we validated the HNN-Tox with several external toxic chemical datasets on a large scale. The HNN-Tox performed optimally or better than the other machine-learning methods for diverse chemicals. This study is the first to report a large-scale prediction of dose-range chemical toxicity with varying features. The HNN-Tox has broad applicability in predicting toxicity for diverse chemicals and could serve as an alternative methodology approach to animal-based toxicity assessment.
| S-EPMC9692315 | biostudies-literature
Project description:The ability to design proteins with high affinity and selectivity for any given small molecule is a rigorous test of our understanding of the physiochemical principles that govern molecular recognition. Attempts to rationally design ligand-binding proteins have met with little success, however, and the computational design of protein-small-molecule interfaces remains an unsolved problem. Current approaches for designing ligand-binding proteins for medical and biotechnological uses rely on raising antibodies against a target antigen in immunized animals and/or performing laboratory-directed evolution of proteins with an existing low affinity for the desired ligand, neither of which allows complete control over the interactions involved in binding. Here we describe a general computational method for designing pre-organized and shape complementary small-molecule-binding sites, and use it to generate protein binders to the steroid digoxigenin (DIG). Of seventeen experimentally characterized designs, two bind DIG; the model of the higher affinity binder has the most energetically favourable and pre-organized interface in the design set. A comprehensive binding-fitness landscape of this design, generated by library selections and deep sequencing, was used to optimize its binding affinity to a picomolar level, and X-ray co-crystal structures of two variants show atomic-level agreement with the corresponding computational models. The optimized binder is selective for DIG over the related steroids digitoxigenin, progesterone and β-oestradiol, and this steroid binding preference can be reprogrammed by manipulation of explicitly designed hydrogen-bonding interactions. The computational design method presented here should enable the development of a new generation of biosensors, therapeutics and diagnostics.
| S-EPMC3898436 | biostudies-literature
Project description:Long-range target detection in thermal infrared imagery is a challenging research problem due to the low resolution and limited detail captured by thermal sensors. The limited size and variability in thermal image datasets for small target detection is also a major constraint for the development of accurate and robust detection algorithms. To address both the sensor and data constraints, we propose a novel convolutional neural network (CNN) feature extraction architecture designed for small object detection in data-limited settings. More specifically, we focus on long-range ground-based thermal vehicle detection, but also show the effectiveness of the proposed algorithm on drone and satellite aerial imagery. The design of the proposed architecture is inspired by an analysis of popular object detectors as well as custom-designed networks. We find that restricted receptive fields (rather than more globalized features, as is the trend), along with less downsampling of feature maps and attenuated processing of fine-grained features, lead to greatly improved detection rates while mitigating the model's capacity to overfit on small or poorly varied datasets. Our approach achieves state-of-the-art results on the Defense Systems Information Analysis Center (DSIAC) automated target recognition (ATR) and the Tiny Object Detection in Aerial Images (AI-TOD) datasets.
| S-EPMC10536622 | biostudies-literature