Enabling population protein dynamics through Bayesian modeling

The knowledge of protein dynamics or turnover in patients provides invaluable information related to certain diseases, drug efficacy, or biological processes. A great corpus of experimental and computational methods has been developed, including by us, in the case of human patients followed <em>in vivo</em>. Moving one step further, we propose here a new modeling approach to capture the highly relevant notion of population protein dynamics. Using two data sets, we show that models inspired by population pharmacokinetics can accurately capture protein turnover within a cohort of individuals, even in presence of substantial inter-individual variability. R scripts<br> --------- The script BUGS-SILK-7.R computes protein dynamics models for one patient at a time. It was used<br> to explore the possibility to apply Bayesian modeling to SILK data. The script BUGS-pop-models.R computes the population protein dynamics Bayesian models. Note that the script should work provided they are run with all the data folders<br> being subfolders of the working directory. <br> Data folders<br> ------------ HRMS (unbiased proteomics protocol, see Lehmann, et al., Anal Chem, 2019) folders contain<br> the results of the QNB analysis. One folder per fluid (CSF-HRMS-05-2021 &amp; plasma-HRMS-05-2021).<br> Analysis results are provided in a tabulated text file in each case with individual protein<br> parameters. These parameter values are used to build prior distributions in the new paper. TQ (targeted MRM protocol) folders contain the results of the QNB analysis, which are the input<br> data for this paper. These are stored in the CSF-TQ-05-2021/models &amp; plasma-TQ-05-2021/models<br> folders. In each case, a plot of the protein model is provided for reference (.pdf files)<br> as well as a table listing each observed RIA (-final-selection.txt files) and another table<br> providing the QNB model parameters (-bootstrap.txt files). In addition, thq TQ folders contain the results of the Bayesian population modeling. The<br> folders CSF-TQ-05-2021/BUGS-models &amp; plasma-TQ-05-2021/BUGS-models contain the models<br> obtained with the non-informative prior.The folders CSF-TQ-05-2021/BUGS-models-prior &amp;<br> plasma-TQ-05-2021/BUGS-models-prior contain the models obtained with the informative prior<br> built from the HRMS data. HRMS and TQ folders are respectively zipped under HRMS.folders.plasma-and-csf.zip <br> and TQ.folders.plasma-and-csf.zip files.