Comments (4)
Dear Anna,
Here would be a csv file of the fit table (not filtered for only hits, without the plot_points
and the plot_curve
columns.
Is that what you were asking for? I am running the analysis on Windows (in case you were wondering).
You should be able to run all of these analyses yourself with the datasets we provide in protti. For this case we use the rapamycin_dose_response
data set from this publication.
The rapamycin data set is a good one for benchmarking, as in this experiment we expect very few off-target hits or hits that are difficult to explain (secondary effects). We also reduced the size of the original data set by only including some random proteins and the target.
You could create an artificial data set with known ground-truth using theprotti
function create_synthetic_data()
and method = "dose_response"
. (This might be the cleanest solution to your problem/question)
You could otherwise also try this data set here: https://www.ebi.ac.uk/pride/archive/projects/PXD038768
You should find some code for the analysis with protti
, as well as all other relevant files in there.
In this experiment we were identifying calmodulin binding sites on the retinal CNG channel with LiP-MS. This data set is different from the other one - here we are studying protein-protein interactions instead of protein-drug interactions.
As calmodulin also binds to other proteins, we also identify other proteins as targets, some of which we cannot explain.
I hope this helps!
All the best,
Dina
from protti.
Thanks a lot for this thorough reply Dina! I really appreciate it.
I'm going to use protti internal data for my benchmarking tests for the moment. The main issue that I noticed with the drc
package is that the results, especially the EC50, could change a bit different depending on the operating system where the package is run. Usually it's a matter of rounding but it can be tricky to get good reprodubility.
These good hits from your dataset are really helpful as I can at least check that the dose-response curve makes sense for the benchmarking data.
Thanks a lot again!
Anna
from protti.
Hi Dina,
I'm re-opening this issue with a quick question about the rapamycin_dose_response
data.
Do you confirm that the rapamycin_dose_response
data in the protti
package is the same data used to produce the dose-response curves in Fig 1d of Piazza 2020 ?
I'm asking this because from the protti
analysis workflow here I see that the example dose-response curve for _VFDVELLKLE_.2
used data that reach a dose concentration 10^8
but in Fig 1d the dose only reaches 10^4.
![Screenshot 2023-10-30 at 12 53 13 pm](https://private-user-images.githubusercontent.com/7087258/278918163-2baa2e30-0e83-41fb-b135-f772a846337a.png?jwt=eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJpc3MiOiJnaXRodWIuY29tIiwiYXVkIjoicmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbSIsImtleSI6ImtleTUiLCJleHAiOjE3MTk3OTY4MTgsIm5iZiI6MTcxOTc5NjUxOCwicGF0aCI6Ii83MDg3MjU4LzI3ODkxODE2My0yYmFhMmUzMC0wZTgzLTQxZmItYjEzNS1mNzcyYTg0NjMzN2EucG5nP1gtQW16LUFsZ29yaXRobT1BV1M0LUhNQUMtU0hBMjU2JlgtQW16LUNyZWRlbnRpYWw9QUtJQVZDT0RZTFNBNTNQUUs0WkElMkYyMDI0MDcwMSUyRnVzLWVhc3QtMSUyRnMzJTJGYXdzNF9yZXF1ZXN0JlgtQW16LURhdGU9MjAyNDA3MDFUMDExNTE4WiZYLUFtei1FeHBpcmVzPTMwMCZYLUFtei1TaWduYXR1cmU9ZGZiYmNkYjZmOWM1MWI0MThkYWIzOWQ2NzUyMjFkMjA0MWM3YjcyMmViYzZlM2IwYWJkNjFmYjRkZTcwNmNiMyZYLUFtei1TaWduZWRIZWFkZXJzPWhvc3QmYWN0b3JfaWQ9MCZrZXlfaWQ9MCZyZXBvX2lkPTAifQ.XB8F6SjOA-rXI_uqVPY6XDtz38Lx2GvHBPnW3gk7nHU)
Thanks for your help!
Anna
from protti.
Hi Anna,
the answer is that the unit is different. As you can see the protti plot is in pM while the plot in the paper is in nM. We chose to use pM so all the concentration points don't have decimal points. But technically it should also work with decimal point concentrations.
The highest concentration in the paper is 10^4
which is 1e+05
nM while ours is 1e+08
pM, which is also 1e+05
nM.
So the plot should be the same between paper and protti, but the dataset rapamycin_dose_response
only contains a subset of the whole dataset. It contains all peptides from FKBP1A (the target protein) and 39 more randomly sampled proteins from the dataset including their peptides. The reason is the dataset size.
from protti.
Related Issues (20)
- Facetting plots in qc_sequence_coverage() HOT 1
- Partial output of fetch_uniprot.R HOT 1
- Query regarding the drc_4p_plot HOT 2
- 'Timeout is reached' in fetch_uniprot HOT 15
- `qc_sample_correlation()` always returns a plot HOT 1
- Add ungroup to normalise function HOT 2
- Add grouping argument to `calculate_sequence_coverage()` HOT 1
- Styler package for GitHub Actions
- qc_sequence_coverage() HOT 2
- assign_missingness() variables are hard coded HOT 1
- fetch_kegg.R doesnt join dataframe correctly HOT 3
- GO enrichment, enrichment type HOT 1
- Improvement to `assign_peptide_type()` HOT 1
- `analyse_functional_network()` update default version HOT 1
- fetch_uniprot accession HOT 4
- fetch_pdb() sodium ion
- fetch_uniprot - fails randomly HOT 1
- qc_cvs factor compatibility HOT 1
- Make minimal peptides requirement flexible in calculate_protein_abundance() HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from protti.