Quantitative description of the interaction between folate and the folate-binding protein from cow's milk

A detailed study has been carried out on the dependence of folate binding on the concentration of FBP (folate-binding protein) at pH 5.0, conditions selected to prevent complications arising from the pre-existing self-association of the acceptor. In contrast with the mandatory requirement that reversible interaction of ligand with a single acceptor site should exhibit a unique, rectangular hyperbolic binding curve, results obtained by ultrafiltration for the FBP-folate system required description in terms of (i) a sigmoidal relationship between concentrations of bound and free folate and (ii) an inverse dependence of affinity on FBP concentration. These findings have been attributed to the difficulties in determining the free ligand concentration in the FBP-folate mixtures for which reaction is essentially stoichiometric. This explanation also accounts for the similar published behaviour of the FBP-folate system at neutral pH, which had been attributed erroneously to acceptor self-association, a phenomenon incompatible with the experimental findings because of its prediction of a greater affinity for folate with increasing FBP concentration.

On systems thinking, systems biology and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

Quantitative analysis of immunolabeling for serotonin and for glutamate transporters after administration of imipramine and citalopram

Serotonin (5-hydroxytryptamine, 5-HT) is an amine neurotransmitter derived from tryptophan and is important in brain systems regulating mood, emotional behavior, and sleep. Selective serotonin reuptake inhibitor (SSRI) drugs are used to treat disorders such as depression, stress, eating disorders, autism, and schizophrenia. It is thought that these drugs act to prolong the action of 5-HT by blocking reuptake. This may lead to decreased 5-HT content in the nerve fibers themselves; however, this has not previously been directly demonstrated. We have studied the effects of administration of two drugs, imipramine and citalopram, on levels of 5-HT in nerve fibers in the murine brain. Quantitative analysis of the areal density of 5-HT fibers throughout the brain was performed using ImageJ software. While a high density of fibers was observed in mid- and hind-brain regions and areas such as thalamus and hypothalamus, densities were far lower in areas such as cortex, where SSRIs might be thought to exert their actions. As anticipated, imipramine and citalopram produced a decline in 5-HT levels in nerve fibers, but the result was not uniform. Areas such as inferior colliculus showed significant reduction whereas little, if any, change was observed in the adjacent superior colliculus. The reason for, and significance of, this regionality is unclear. It has been proposed that serotonin effects in the brain might be linked to changes in glutamatergic transmission. Extracellular glutamate levels are regulated primarily by glial glutamate transporters. Qualitative evaluation of glutamate transporter immunolabeling in cortex of control and drug-treated mice revealed no discernable difference in intensity of glutamate transporter immunoreactivity. These data suggest that changes in intracellular and extracellular levels of serotonin do not cause concomitant changes in astroglial glutamate transporter expression, and thus cannot represent a mechanism for the delayed efficacy of antidepressants when administered clinically. © 2005 Elsevier B.V. All rights reserved.

Coadaptation: A unifying principle in evolutionary thermal biology

Over the last 50 yr, thermal biology has shifted from a largely physiological science to a more integrated science of behavior, physiology, ecology, and evolution. Today, the mechanisms that underlie responses to environmental temperature are being scrutinized at levels ranging from genes to organisms. From these investigations, a theory of thermal adaptation has emerged that describes the evolution of thermoregulation, thermal sensitivity, and thermal acclimation. We review and integrate current models to form a conceptual model of coadaptation. We argue that major advances will require a quantitative theory of coadaptation that predicts which strategies should evolve in specific thermal environments. Simply combining current models, however, is insufficient to understand the responses of organisms to thermal heterogeneity; a theory of coadaptation must also consider the biotic interactions that influence the net benefits of behavioral and physiological strategies. Such a theory will be challenging to develop because each organism's perception of and response to thermal heterogeneity depends on its size, mobility, and life span. Despite the challenges facing thermal biologists, we have never been more pressed to explain the diversity of strategies that organisms use to cope with thermal heterogeneity and to predict the consequences of thermal change for the diversity of communities.

Quantitative analysis of DNA-protein interactions using double-labeled native gel electrophoresis and fluorescence-based imaging

We have developed a sensitive, non-radioactive method to assess the interaction of transcription factors/DNA-binding proteins with DNA. We have modified the traditional radiolabeled DNA gel mobility shift assay to incorporate a DNA probe end-labeled with a Texas-red fluorophore and a DNA-binding protein tagged with the green fluorescent protein to monitor precisely DNA-protein complexation by native gel electrophoresis. We have applied this method to the DNA-binding proteins telomere release factor-1 and the sex-determining region-Y, demonstrating that the method is sensitive (able to detect 100 fmol of fluorescently labeled DNA), permits direct visualization of both the DNA probe and the DNA-binding protein, and enables quantitative analysis of DNA and protein complexation, and thereby an estimation of the stoichiometry of protein-DNA binding.

Quantitative analysis of MC1R gene expression in human skin cell cultures

To address the issue of melanocortin-1 receptor (MC1R) expression in non-melanocytic cells, we have quantitatively evaluated the relative expression levels of both MC1R mRNA and protein in a subset of different cell types. Using semi-quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) at high cycle numbers, we detected MC1R mRNA in all cell types examined, including human embryonic kidney-293 (HEK 293) cells, a cell type widely used as a negative control in melanocortin expression studies. Quantitative real-time PCR revealed the highest levels of MC1R transcripts were in melanocytic cells, whereas the keratinocyte and fibroblast cell cultures examined had only a low level of expression, similar to that of HEK 293 cells. Antibody mediated detection of MC1R protein in membrane extracts demonstrated exogenous receptor in MC1R transfected cell lines, as well as endogenous MC1R in melanoma cells. However, radioligand binding procedures were required to detect MC1R protein of normal human melanocytes and no surface expression of MC1R was detected in any of the non-melanocytic cells examined. This was consistent with their low level of mRNA, and suggests that, if present, the levels of surface receptor are significantly lower than that in melanocytes. The capacity of such limited levels of MC1R protein to influence non-melanocytic skin cell biology would likely be severely compromised. Indeed, the MC1R agonist [NIe(4), D-Phe(7)] alpha-melanocyte stimulating hormone (NDP-MSH) was unable to elevate intracellular cyclic adenosine monophosphate (cAMP) levels in the keratinocyte and fibroblast cells examined, whereas a robust increase was elicited in melanocytes. Although there are a variety of cell types with detectable MC1R mRNA, the expression of physiologically significant levels of the receptor may be more restricted than the current literature indicates, and within epidermal tissue may be limited to the melanocyte

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Quantitative fluorescence spectra and quantum yield map of synthetic pheomelanin

Spectroscopic studies of pheomelanin and its constituents have been sparse. These data present what is by far the most complete description of the fluorescence characteristics of synthetic pheomelanin. Emission spectra between 260 and 600 nm were acquired,for excitation wavelengths between 250 and 500 nm at 1-nm intervals. A quantum yield map is also presented, correcting the fluorescence intensities for differences in species concentration and molar absorptivity. These fluorescence features exhibit interesting similarities and differences to eumelanin, and these data are interpreted with respect to possible chemical structures. Overall, these data suggest that pheomelanin oligomers may be more tightly coupled than those of eumelanin. Finally, the quantum yield is shown to be on the order of 10(-4) and exhibit a complex dependence on excitation energy, varying by a factor of 4 across the energies employed here. (c) 2006 Wiley Periodicals, Inc.

Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Peptide binding to HLA-DP proteins at pH 5.0 and pH 7.0:a quantitative molecular docking study

Background: HLA-DPs are class II MHC proteins mediating immune responses to many diseases. Peptides bind MHC class II proteins in the acidic environment within endosomes. Acidic pH markedly elevates association rate constants but dissociation rates are almost unchanged in the pH range 5.0 - 7.0. This pH-driven effect can be explained by the protonation/deprotonation states of Histidine, whose imidazole has a pKa of 6.0. At pH 5.0, imidazole ring is protonated, making Histidine positively charged and very hydrophilic, while at pH 7.0 imidazole is unprotonated, making Histidine less hydrophilic. We develop here a method to predict peptide binding to the four most frequent HLA-DP proteins: DP1, DP41, DP42 and DP5, using a molecular docking protocol. Dockings to virtual combinatorial peptide libraries were performed at pH 5.0 and pH 7.0. Results: The X-ray structure of the peptide - HLA-DP2 protein complex was used as a starting template to model by homology the structure of the four DP proteins. The resulting models were used to produce virtual combinatorial peptide libraries constructed using the single amino acid substitution (SAAS) principle. Peptides were docked into the DP binding site using AutoDock at pH 5.0 and pH 7.0. The resulting scores were normalized and used to generate Docking Score-based Quantitative Matrices (DS-QMs). The predictive ability of these QMs was tested using an external test set of 484 known DP binders. They were also compared to existing servers for DP binding prediction. The models derived at pH 5.0 predict better than those derived at pH 7.0 and showed significantly improved predictions for three of the four DP proteins, when compared to the existing servers. They are able to recognize 50% of the known binders in the top 5% of predicted peptides. Conclusions: The higher predictive ability of DS-QMs derived at pH 5.0 may be rationalised by the additional hydrogen bond formed between the backbone carbonyl oxygen belonging to the peptide position before p1 (p-1) and the protonated ε-nitrogen of His 79β. Additionally, protonated His residues are well accepted at most of the peptide binding core positions which is in a good agreement with the overall negatively charged peptide binding site of most MHC proteins. © 2012 Patronov et al.; licensee BioMed Central Ltd.

Revealing 4-hydroxynonenal-histidine adducts as qualitative and quantitative biomarker of oxidative stress, lipid peroxidation and oxidative homeostasis

Findings on growth regulating activities of the end-product of lipid peroxidation 4-hydroxy-2-nonenal (HNE), which acts as a “second messenger of free radicals”, overlapped with the development of antibodies specific for the aldehyde-protein adducts. These led to qualitative immunochemical determinations of the HNE presence in various pathophysiological processes and to the change of consideration of the aldehyde’s bioactivities from toxicity into cell signalling. Moreover, findings of the HNE-protein adduct in various organs under physiological circumstances support the concept of “oxidative homeostasis”, which implies that oxidative stress and lipid peroxidation are not only pathological but also physiological processes. Reactive aldehydes, at least HNE, could play important role in oxidative homeostasis, while complementary research approaches might reveal the relevance of the aldehydic-protein adducts as major biomarkers of oxidative stress, lipid peroxidation and oxidative homeostasis. Aiming to join efforts in such research activities researchers interacting through the International 4-Hydroxynonenal Club acting within the SFRR-International and through networking projects of the system of the European Cooperation in Science and Technology (COST) carried validation of the methods for lipid peroxidation and further developed the genuine 4-HNE-His ELISA founding quantitative and qualitative methods for detection of 4-HNE-His adducts as valuable tool to study oxidative stress and lipid peroxidation in cell cultures, various organs and tissues and eventually for human plasma and serum analyses [1]. Reference: 1. Weber, Daniela. Lidija, Milkovic. Measurement of HNE-protein adducts in human plasma and serum by ELISA—Comparison of two primary antibodies. Redox Biol. 2013. 226-233.

The use of quantitative polymerase chain reaction amplification method to determine sex ratio of human spermatozoa in semen and to determine maternal contamination of amniotic fluid samples

We have modified a technique which uses a single pair of primer sets directed against homologous but distinct genes on the X and Y chromosomes, all of which are coamplified in the same reaction tube with trace amounts of radioactivity. The resulting bands are equal in length, yet distinguishable by restriction enzyme sites generating two independent bands, a 364 bp X-specific band and a 280 bp Y-specific band. A standard curve was generated to show the linear relationship between X/Y ratio average vs. %Y or %X chromosomal content. Of the 51 purified amniocyte DNA samples analyzed, 16 samples showed evidence of high % X contamination while 2 samples demonstrated higher % Y than the expected 50% X and 50% Y chromosomal content. With regards to the 25 processed sperm samples analyzed, X-sperm enrichment was evident when compared to the primary sex ratio whereas Y-sperm was enriched when we compared before and after selection samples.

Quantitative Comparison of Plant Community Hydrology Using Large-Extent, Long-Term Data

Large-extent vegetation datasets that co-occur with long-term hydrology data provide new ways to develop biologically meaningful hydrologic variables and to determine plant community responses to hydrology. We analyzed the suitability of different hydrological variables to predict vegetation in two water conservation areas (WCAs) in the Florida Everglades, USA, and developed metrics to define realized hydrologic optima and tolerances. Using vegetation data spatially co-located with long-term hydrological records, we evaluated seven variables describing water depth, hydroperiod length, and number of wet/dry events; each variable was tested for 2-, 4- and 10-year intervals for Julian annual averages and environmentally-defined hydrologic intervals. Maximum length and maximum water depth during the wet period calculated for environmentally-defined hydrologic intervals over a 4-year period were the best predictors of vegetation type. Proportional abundance of vegetation types along hydrological gradients indicated that communities had different realized optima and tolerances across WCAs. Although in both WCAs, the trees/shrubs class was on the drier/shallower end of hydrological gradients, while slough communities occupied the wetter/deeper end, the distribution ofCladium, Typha, wet prairie and Salix communities, which were intermediate for most hydrological variables, varied in proportional abundance along hydrologic gradients between WCAs, indicating that realized optima and tolerances are context-dependent.

Improving the use of G-CSF during chemotherapy using physiological mathematical modelling : a quantitative systems pharmacology approach

La diminution des doses administrées ou même la cessation complète d'un traitement chimiothérapeutique est souvent la conséquence de la réduction du nombre de neutrophiles, qui sont les globules blancs les plus fréquents dans le sang. Cette réduction dans le nombre absolu des neutrophiles, aussi connue sous le nom de myélosuppression, est précipitée par les effets létaux non spécifiques des médicaments anti-cancéreux, qui, parallèlement à leur effet thérapeutique, produisent aussi des effets toxiques sur les cellules saines. Dans le but d'atténuer cet impact myélosuppresseur, on administre aux patients un facteur de stimulation des colonies de granulocytes recombinant humain (rhG-CSF), une forme exogène du G-CSF, l'hormone responsable de la stimulation de la production des neutrophiles et de leurs libération dans la circulation sanguine. Bien que les bienfaits d'un traitement prophylactique avec le G-CSF pendant la chimiothérapie soient bien établis, les protocoles d'administration demeurent mal définis et sont fréquemment déterminés ad libitum par les cliniciens. Avec l'optique d'améliorer le dosage thérapeutique et rationaliser l'utilisation du rhG-CSF pendant le traitement chimiothérapeutique, nous avons développé un modèle physiologique du processus de granulopoïèse, qui incorpore les connaissances actuelles de pointe relatives à la production des neutrophiles des cellules souches hématopoïétiques dans la moelle osseuse. À ce modèle physiologique, nous avons intégré des modèles pharmacocinétiques/pharmacodynamiques (PK/PD) de deux médicaments: le PM00104 (Zalypsis®), un médicament anti-cancéreux, et le rhG-CSF (filgrastim). En se servant des principes fondamentaux sous-jacents à la physiologie, nous avons estimé les paramètres de manière exhaustive sans devoir recourir à l'ajustement des données, ce qui nous a permis de prédire des données cliniques provenant de 172 patients soumis au protocol CHOP14 (6 cycles de chimiothérapie avec une période de 14 jours où l'administration du rhG-CSF se fait du jour 4 au jour 13 post-chimiothérapie). En utilisant ce modèle physio-PK/PD, nous avons démontré que le nombre d'administrations du rhG-CSF pourrait être réduit de dix (pratique actuelle) à quatre ou même trois administrations, à condition de retarder le début du traitement prophylactique par le rhG-CSF. Dans un souci d'applicabilité clinique de notre approche de modélisation, nous avons investigué l'impact de la variabilité PK présente dans une population de patients, sur les prédictions du modèle, en intégrant des modèles PK de population (Pop-PK) des deux médicaments. En considérant des cohortes de 500 patients in silico pour chacun des cinq scénarios de variabilité plausibles et en utilisant trois marqueurs cliniques, soient le temps au nadir des neutrophiles, la valeur du nadir, ainsi que l'aire sous la courbe concentration-effet, nous avons établi qu'il n'y avait aucune différence significative dans les prédictions du modèle entre le patient-type et la population. Ceci démontre la robustesse de l'approche que nous avons développée et qui s'apparente à une approche de pharmacologie quantitative des systèmes (QSP). Motivés par l'utilisation du rhG-CSF dans le traitement d'autres maladies, comme des pathologies périodiques telles que la neutropénie cyclique, nous avons ensuite soumis l'étude du modèle au contexte des maladies dynamiques. En mettant en évidence la non validité du paradigme de la rétroaction des cytokines pour l'administration exogène des mimétiques du G-CSF, nous avons développé un modèle physiologique PK/PD novateur comprenant les concentrations libres et liées du G-CSF. Ce nouveau modèle PK a aussi nécessité des changements dans le modèle PD puisqu’il nous a permis de retracer les concentrations du G-CSF lié aux neutrophiles. Nous avons démontré que l'hypothèse sous-jacente de l'équilibre entre la concentration libre et liée, selon la loi d'action de masse, n'est plus valide pour le G-CSF aux concentrations endogènes et mènerait en fait à la surestimation de la clairance rénale du médicament. En procédant ainsi, nous avons réussi à reproduire des données cliniques obtenues dans diverses conditions (l'administration exogène du G-CSF, l'administration du PM00104, CHOP14). Nous avons aussi fourni une explication logique des mécanismes responsables de la réponse physiologique aux deux médicaments. Finalement, afin de mettre en exergue l’approche intégrative en pharmacologie adoptée dans cette thèse, nous avons démontré sa valeur inestimable pour la mise en lumière et la reconstruction des systèmes vivants complexes, en faisant le parallèle avec d’autres disciplines scientifiques telles que la paléontologie et la forensique, où une approche semblable a largement fait ses preuves. Nous avons aussi discuté du potentiel de la pharmacologie quantitative des systèmes appliquées au développement du médicament et à la médecine translationnelle, en se servant du modèle physio-PK/PD que nous avons mis au point.

Improving the use of G-CSF during chemotherapy using physiological mathematical modelling : a quantitative systems pharmacology approach

La diminution des doses administrées ou même la cessation complète d'un traitement chimiothérapeutique est souvent la conséquence de la réduction du nombre de neutrophiles, qui sont les globules blancs les plus fréquents dans le sang. Cette réduction dans le nombre absolu des neutrophiles, aussi connue sous le nom de myélosuppression, est précipitée par les effets létaux non spécifiques des médicaments anti-cancéreux, qui, parallèlement à leur effet thérapeutique, produisent aussi des effets toxiques sur les cellules saines. Dans le but d'atténuer cet impact myélosuppresseur, on administre aux patients un facteur de stimulation des colonies de granulocytes recombinant humain (rhG-CSF), une forme exogène du G-CSF, l'hormone responsable de la stimulation de la production des neutrophiles et de leurs libération dans la circulation sanguine. Bien que les bienfaits d'un traitement prophylactique avec le G-CSF pendant la chimiothérapie soient bien établis, les protocoles d'administration demeurent mal définis et sont fréquemment déterminés ad libitum par les cliniciens. Avec l'optique d'améliorer le dosage thérapeutique et rationaliser l'utilisation du rhG-CSF pendant le traitement chimiothérapeutique, nous avons développé un modèle physiologique du processus de granulopoïèse, qui incorpore les connaissances actuelles de pointe relatives à la production des neutrophiles des cellules souches hématopoïétiques dans la moelle osseuse. À ce modèle physiologique, nous avons intégré des modèles pharmacocinétiques/pharmacodynamiques (PK/PD) de deux médicaments: le PM00104 (Zalypsis®), un médicament anti-cancéreux, et le rhG-CSF (filgrastim). En se servant des principes fondamentaux sous-jacents à la physiologie, nous avons estimé les paramètres de manière exhaustive sans devoir recourir à l'ajustement des données, ce qui nous a permis de prédire des données cliniques provenant de 172 patients soumis au protocol CHOP14 (6 cycles de chimiothérapie avec une période de 14 jours où l'administration du rhG-CSF se fait du jour 4 au jour 13 post-chimiothérapie). En utilisant ce modèle physio-PK/PD, nous avons démontré que le nombre d'administrations du rhG-CSF pourrait être réduit de dix (pratique actuelle) à quatre ou même trois administrations, à condition de retarder le début du traitement prophylactique par le rhG-CSF. Dans un souci d'applicabilité clinique de notre approche de modélisation, nous avons investigué l'impact de la variabilité PK présente dans une population de patients, sur les prédictions du modèle, en intégrant des modèles PK de population (Pop-PK) des deux médicaments. En considérant des cohortes de 500 patients in silico pour chacun des cinq scénarios de variabilité plausibles et en utilisant trois marqueurs cliniques, soient le temps au nadir des neutrophiles, la valeur du nadir, ainsi que l'aire sous la courbe concentration-effet, nous avons établi qu'il n'y avait aucune différence significative dans les prédictions du modèle entre le patient-type et la population. Ceci démontre la robustesse de l'approche que nous avons développée et qui s'apparente à une approche de pharmacologie quantitative des systèmes (QSP). Motivés par l'utilisation du rhG-CSF dans le traitement d'autres maladies, comme des pathologies périodiques telles que la neutropénie cyclique, nous avons ensuite soumis l'étude du modèle au contexte des maladies dynamiques. En mettant en évidence la non validité du paradigme de la rétroaction des cytokines pour l'administration exogène des mimétiques du G-CSF, nous avons développé un modèle physiologique PK/PD novateur comprenant les concentrations libres et liées du G-CSF. Ce nouveau modèle PK a aussi nécessité des changements dans le modèle PD puisqu’il nous a permis de retracer les concentrations du G-CSF lié aux neutrophiles. Nous avons démontré que l'hypothèse sous-jacente de l'équilibre entre la concentration libre et liée, selon la loi d'action de masse, n'est plus valide pour le G-CSF aux concentrations endogènes et mènerait en fait à la surestimation de la clairance rénale du médicament. En procédant ainsi, nous avons réussi à reproduire des données cliniques obtenues dans diverses conditions (l'administration exogène du G-CSF, l'administration du PM00104, CHOP14). Nous avons aussi fourni une explication logique des mécanismes responsables de la réponse physiologique aux deux médicaments. Finalement, afin de mettre en exergue l’approche intégrative en pharmacologie adoptée dans cette thèse, nous avons démontré sa valeur inestimable pour la mise en lumière et la reconstruction des systèmes vivants complexes, en faisant le parallèle avec d’autres disciplines scientifiques telles que la paléontologie et la forensique, où une approche semblable a largement fait ses preuves. Nous avons aussi discuté du potentiel de la pharmacologie quantitative des systèmes appliquées au développement du médicament et à la médecine translationnelle, en se servant du modèle physio-PK/PD que nous avons mis au point.