Modelling and optimal control of immune response of renal transplant recipients

We consider the increasingly important and highly complex immunological control problem: control of the dynamics of immunosuppression for organ transplant recipients. The goal in this problem is to maintain the delicate balance between over-suppression (where opportunistic latent viruses threaten the patient) and under-suppression (where rejection of the transplanted organ is probable). First, a mathematical model is formulated to describe the immune response to both viral infection and introduction of a donor kidney in a renal transplant recipient. Some numerical results are given to qualitatively validate and demonstrate that this initial model exhibits appropriate characteristics of primary infection and reactivation for immunosuppressed transplant recipients. In addition, we develop a computational framework for designing adaptive optimal treatment regimes with partial observations and low-frequency sampling, where the state estimates are obtained by solving a second deterministic optimal tracking problem. Numerical results are given to illustrate the feasibility of this method in obtaining optimal treatment regimes with a balance between under-suppression and over-suppression of the immune system.

[1]  Duan Li,et al.  On Restart Procedures for the Conjugate Gradient Method , 2004, Numerical Algorithms.

[2]  P. Koskinen,et al.  Prospective follow-up of primary CMV infections after 6 months of valganciclovir prophylaxis in renal transplant recipients. , 2008, Nephrology, dialysis, transplantation : official publication of the European Dialysis and Transplant Association - European Renal Association.

[3]  J. Janssen,et al.  Deterministic and Stochastic Optimal Control , 2013 .

[4]  H. Banks,et al.  A zero-sum electromagnetic evader–interrogator differential game with uncertainty , 2012 .

[5]  H. T. Banks,et al.  A model for HCMV infection in immunosuppressed patients , 2009, Math. Comput. Model..

[6]  Franz Kappel,et al.  Comparison of optimal design methods in inverse problems , 2011, Inverse problems.

[7]  Stephen P. Boyd,et al.  Receding Horizon Control , 2011, IEEE Control Systems.

[8]  H. Rollag,et al.  A PROSPECTIVE STUDY OF THE NATURAL COURSE OF CYTOMEGALOVIRUS INFECTION AND DISEASE IN RENAL ALLOGRAFT RECIPIENTS1 , 2000, Transplantation.

[9]  Harvey Thomas Banks,et al.  Feedback control of HIV antiviral therapy with long measurement time , 2011 .

[10]  T. Schumacher,et al.  Differentiation of cytomegalovirus-specific CD8(+) T cells in healthy and immunosuppressed virus carriers. , 2001, Blood.

[11]  P. Griffiths,et al.  Prediction of cytomegalovirus load and resistance patterns after antiviral chemotherapy. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[12]  A. Levey,et al.  A More Accurate Method To Estimate Glomerular Filtration Rate from Serum Creatinine: A New Prediction Equation , 1999, Annals of Internal Medicine.

[13]  B. Adams,et al.  Dynamic multidrug therapies for hiv: optimal and sti control approaches. , 2004, Mathematical biosciences and engineering : MBE.

[14]  H. Banks,et al.  Mathematical and Experimental Modeling of Physical and Biological Processes , 2009 .

[15]  Harvey Thomas Banks,et al.  Receding Horizon Control of HIV , 2011 .

[16]  M. Bonneville,et al.  Cross-Reactivity of Herpesvirus-Specific CD8 T Cell Lines Toward Allogeneic Class I MHC Molecules , 2010, PloS one.

[17]  Shuhua Hu,et al.  Modelling HIV immune response and validation with clinical data , 2008, Journal of biological dynamics.

[18]  M. Pescovitz.  Review of the CMV in renal transplantation. , 2007, Saudi journal of kidney diseases and transplantation : an official publication of the Saudi Center for Organ Transplantation, Saudi Arabia.

[19]  C. Ki,et al.  Real-time PCR assay compared with antigenemia assay for detecting cytomegalovirus infection in kidney transplant recipients. , 2007, Transplantation proceedings.

[20]  T. F. Smith,et al.  Cytomegalovirus (CMV) DNA load predicts relapsing CMV infection after solid organ transplantation. , 2000, The Journal of infectious diseases.

[21]  H. Banks,et al.  Dynamic electromagnetic evasion-pursuit games with uncertainty , 2010 .

[22]  Harvey Thomas Banks,et al.  A state‐dependent Riccati equation‐based estimator approach for HIV feedback control , 2006 .

[23]  Harvey Thomas Banks,et al.  Generalized sensitivities and optimal experimental design , 2010 .

[24]  Jorge Nocedal,et al.  Global Convergence Properties of Conjugate Gradient Methods for Optimization , 1992, SIAM J. Optim..

[25]  B. Adams,et al.  HIV dynamics: Modeling, data analysis, and optimal treatment protocols , 2005 .

[26]  Luc Perrin,et al.  Monitoring of Cytomegalovirus Infection in Solid-Organ Transplant Recipients by an Ultrasensitive Plasma PCR Assay , 2003, Journal of Clinical Microbiology.

[27]  Jinhua Guo,et al.  A new algorithm of nonlinear conjugate gradient method with strong convergence , 2008 .

[28]  S. Mitter,et al.  The conjugate gradient method for optimal control problems , 1967 .