Impact Factor 3.201 | CiteScore 3.22
More on impact ›

Editorial ARTICLE

Front. Physiol., 08 June 2018 | https://doi.org/10.3389/fphys.2018.00678

Editorial: Safety Pharmacology – Risk Assessment QT Interval Prolongation and Beyond

  • 1Department of Pharmacology, University of California, Davis, Davis, CA, United States
  • 2Biomedical Signal Interpretation and Computational Simulation Group, Aragón Institute of Engineering Research, IIS Aragón, University of Zaragoza, Zaragoza, Spain
  • 3Department of Computer Science, University of Oxford, Oxford, United Kingdom

The Need of New Paradigms for Cardiac Safety

The scope of safety pharmacology is to predict whether a drug is likely to cause potentially lethal adverse effects if administered to humans. While safety pharmacology has broadened its interests in recent years to the whole cardiovascular, respiratory, and central nervous systems (and is now extending to other body functions), a major focus since its inception has been assessing drug-induced prolongation in the QT interval—a surrogate biomarker for torsades de pointes (TdP) liability. Because the vast majority of drugs that can cause QT prolongation inhibit hERG channels, current regulatory guidelines concerning cardiac safety recommend that all compounds are evaluated in vitro for their hERG inhibitory potency (Redfern et al., 2003) and in vivo for their ability to cause QT/QTc interval prolongation (Food and Drug Administration, 2005) in an appropriate animal model and in humans. However, it has now become apparent that QT/QTc prolongation and hERG block are an insufficient proxy for TdP risk. While the current approach based on these markers has been successful in terms of preventing TdP risk, this regulatory paradigm might lead to withdrawal from the drug development pipeline and clinical use of potentially safe drugs. There is therefore a crucial need to develop a more accurate assessment of proarrhythmic potential of drugs. Notably, in 2014 the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative was proposed as a new strategy by expert working groups sponsored by the US Food and Drug Administration (FDA), the Cardiac Safety Research Consortium (CSRC), and the Health and Environmental Science Institute (HESI), and has quickly become a global effort, also involving many industry and academia partners (Sager et al., 2014). CiPA aims at developing and validating a new paradigm for cardiac safety evaluation of new drugs that provides a more accurate and comprehensive mechanistic-based assessment of proarrhythmic (rather than QT prolonging) potential of drugs (Gintant et al., 2016). This involves assessment of (i) high-throughput in vitro screening of drug effects on multiple human ion channels, (ii) coupled with in silico modeling of human cardiac myocytes to assess integrated electrophysiological responses, and (iii) verification of predicted responses in human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs). Safety pharmacology has evolved in recent years to identify and incorporate new technologies for clinical and non-clinical applications, including refinement of ex vivo and in vitro assays and screens, in vivo models, non-invasive clinical modalities, and in silico approaches. Here we collected a series of review, perspective, and original research articles that summarize the state of our knowledge and the latest advances in these technologies, and how these might contribute to shaping new and improved cardiac safety guidelines.

Multiscale Modeling for Safety Pharmacology

There is a wide range of length and time scales covered in this Research Topic, from the atom and ns to the whole organism and month/year (Figure 1 top left to bottom right), all of which are relevant to safety pharmacology. Structural studies, including modeling of ion channel gating (Perissinotti et al.) and interactions with drugs, and drug partitioning, are critical for drug discovery efforts, and perhaps also a necessary approach for safety considerations. For example, the study by DeMarco et al. utilized all-atom molecular dynamics simulations to show that ionization of drug molecules (specifically Sotalol) can significantly affect their membrane permeability and partitioning kinetics, and should therefore be a consideration in ongoing in silico safety pharmacology efforts. Given the complexity of the interaction between drugs and ion channels, the drug binding kinetics, state dependent binding, and temperature dependence could significantly alter drugs' impact on the action potential (AP), even when drugs display similar steady-state block. Lee et al. highlighted some of the challenges involved in modeling of the hERG channel and also discussed limitations and need for improved voltage-clamp protocols to characterize drug-channel interaction in in vitro experiments. Ellinwood et al. looked at the consequences of drug binding kinetics and state dependence of KV1.5 targeting drugs on atrial electrophysiology, and revealed that ionic remodeling also affects the degree of efficacy and safety of state-specific IKur inhibitors, by modifying the AP trajectory. These studies highlight the potential need for extraordinary detail in the in vitro characterization for accurate in silico prediction of (cardiac-region specific, Morotti et al., 2016; Ellinwood et al., 2017; Ellinwood et al.) drug effects on channels and cardiac electrophysiology.

FIGURE 1
www.frontiersin.org

Figure 1. Approaches to predicting drug cardiotoxicity. The Research Topic includes: (1) studies of ion channel structure and function and drug interaction, ranging from atomistic models to kinetic models of whole-cell ion currents (Top); (2) mechanistic models of single cardiac myocyte electrophysiology, development of metrics for quantification of torsadogenic risk, and population-based approaches accounting for cell-to-cell variability (Middle); and (3) tissue/organ/organism level studies spanning tissue-level to whole heart and torso simulations, ex vivo recordings in animal preparations, and human clinical trials (Bottom).

While the ion channel gating and drug-interaction models might require further refinement and increased complexity, significant efforts have been put forward to improve existing cardiomyocyte models and to take advantage of the existence (and convergence) of competing mathematical models to narrow hypotheses or explore alternative hypotheses (Sarkar and Sobie, 2011; Sánchez et al., 2012; Gemmell et al., 2014; Mann et al., 2016; Pueyo et al., 2016a,b; Gong et al., 2017; Muszkiewicz et al., 2017). For example, recent work has shown that, when forced to reproduce the same data, three competing models of human ventricular myocytes (Ten Tusscher and Panfilov, 2006; Grandi et al., 2010; O'Hara et al., 2011) became substantially more similar than they were originally (Mann et al., 2016). Notably, the work by Krogh-Madsen et al. used clinical congenital LQT data (as done by Mann et al., 2016) and physiological constraints on intracellular ionic concentrations to optimize parameters in the O'Hara-Rudy (ORd) human model (O'Hara et al., 2011). This in turn improved the accuracy and robustness of TdP risk prediction (Lancaster and Sobie, 2016), which the authors attributed to the importance of Ca2+ dynamics in repolarization and to an improved balance of IKs vs. IKr in the new model. A different parameterization of the ORd model by Dutta et al. and Dutta et al. also yielded a better correspondence with drug response data and improved the identification of pro-arrhythmic drugs. The authors developed a new metric qNet, which quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state AP, to separate low-, intermediate-, and high-risk hERG blockers. A follow up study appraised the robustness of qNet as a biomarker for TdP by considering how uncertainty in the model parameters propagates to the phenotype level (Chang et al.). The authors were thus able to identify the conditions under which decisions on risk can be made reliably and objectively. Yet, questions remain regarding the physiological meaning of this new metric, and whether multiple metrics should be utilized that account for a broader range of behaviors and mechanisms. Tixier et al. used an in silico model of multi-electrode array electrophysiology and machine learning to identify predictive biomarkers that should be measured to improve classifications of drugs. These investigations add to several previous efforts to build computational frameworks for assessment of TdP risk (Mirams et al., 2011; Kramer et al., 2013; Lancaster and Sobie, 2016). On the other hand, Parikh et al. showed that a simpler classification method based on direct features (ion channel block information) performed with comparable or higher accuracy than existing methods based on simulated metrics. One potential limitation of this approach is that direct feature classifiers might fail identifying the proarrhythmic risk of drugs affecting channels that are not included in the training set, whereas predictive modeling is more likely to yield an accurate classification.

Biophysical modeling can not only provide means for drug classification, but also understanding of the mechanistic underpinning of drug responses, as in the multiscale simulations by Ni et al. and Colman et al. These studies are important reminders that AP duration changes are rarely homogenous (e.g., there exist gradients—transmurally, or from base to apex) and can increase the tissue-level substrates for arrhythmias (Antzelevitch, 2005; Glukhov et al., 2010). Indeed, multiscale in silico models can be very powerful tools to investigate the response of candidate antiarrhythmic compounds at the level of the electrocardiogram (ECG). The simulated data may also serve to identify novel ECG-derived biomarkers detecting block of inward and/or outward currents based on ECG features (Vicente et al., 2016). Using Langendorff perfused ex vivo rabbit hearts the Fenton group measured and analyzed the complex dynamics of spatially discordant alternans, which provide the substrate for reentrant arrhythmia (Uzelac et al.) The authors noted that current AP models fail to reproduce some key dynamics such as voltage amplitude alternans, smooth development of Ca2+ alternans in time, and conduction. Experimental characterization of these dynamics can inform refining of existing models to analyze mechanisms.

Accounting for Patients' Condition and Inter-Subject Variability

Clinical risk assessment and trial suggest that patient conditions, i.e., sex (Yang et al., 2017; Vorobyov and Clancy, 2018), age, disease, electrolyte imbalance (Lazzerini et al.), interaction with other drugs (Lv et al.) should all be taken into account in risk assessment (Lane and Tinker)—which is not yet addressed by CiPA efforts. Along the same lines, Wisniowska et al. reviewed the different sources of variability (both intrinsic and extrinsic) that exist in the human population in response to drug action, and emphasized the need of accounting for these aspects in modeling approaches for safety pharmacology. Two studies by the Rodriguez group establish the potential of population-based approaches as very powerful in silico tools for safety pharmacology investigations. Passini et al. showed that human in silico drug trials using repolarization abnormality quantification as the main metric do better than animal models in detecting drugs with TdP risk. They also show agreement of in silico predictions and two established experimental models (rabbit wedge ECGs and hiPSC-CMs). Other statistical methods, e.g., logistic regression, have been previously employed to assess the proarrhythmic risk in a population of computational model variants (Lee et al., 2013; Morotti and Grandi, 2017). Calibrated populations of models of heart cells could generally reproduce experimental drug effects on human tissue for dofetilide, whereas lack of agreement between experiments and simulations for quinidine and verapamil suggest further work is needed to understand the more complex electrophysiological effects of these multichannel blocking drugs (Britton et al.).

Use of IPSC-CMS for Safety Pharmacology

Because iPSC-CMs are a readily-obtainable and renewable source of human cardiac myocytes, they are gaining popularity as a platform to screen drugs for toxicity testing. However, given the iPSC immature phenotype, and phenotypic differences across iPSC-CM cell lines (Lei et al.), it remains unclear how well drug tests performed in iPSC-CMs will recapitulate the effects observed in adult human cardiomyocytes and hearts. Koivumaki et al. developed a computational model of the iPSC-CMs that recapitulates the cells' immature phenotype, and explore differences in ionic behavior underlying the AP in paced vs. spontaneous modes, phenotypic variability in iPSC-CMs, and iPSC-CM model's ability to recapitulate physiological properties of adult cells. Recently, statistical methods have also been established to provide accurate predictions of adult myocyte drug responses from iPSC-CM simulations (Gong and Sobie, 2018). iPSC-CM utilization in drug discovery and safety investigations is reviewed by Ortega et al. An important advancement in the technological approach of improving the utility of iPSC-CMs for safety pharmacology is the augmentation of IK1 using dynamic clamp. Plagued by low-throughput, Goversen et al. have moved toward demonstrating that such dynamic clamp can be performed in a high throughput manner. Bjork et al. reported that the expression of optogenetic tools in iPSC-CMs did not significantly affect the baseline electrophysiological properties of these cells, thus allowing electrophysiological assessments comparable to conventional patch clamp studies. Nevertheless, adult human ventricular cardiomyocytes (Nguyen et al.) and trabeculae (Qu et al.) might still be a more reliable model to test the cardiotoxic risk associated with novel drugs, with some advantages over animal and iPSC models.

Conclusions and Future Directions

There is a growing body of work supporting the integration of new and established computational and experimental approaches to understanding and predicting the risk of TdP. While mechanistic systems modeling is mature in the cardiac arrhythmia field, use of similar approaches can improve understanding and prediction of cardiotoxicity caused by other drugs, e.g., cancer therapeutics (Shim et al.). Given the focus on TdP and QT interval, however, the deleterious effects of drugs on cardiac function are evaluated only in terms of changes in electrophysiological properties. Future work should therefore extend the current paradigm to include other major cellular functions (such as contraction, energetics, and cell death, i.e., via apoptosis), which dysregulation can severely impact cardiac performance. In addition to cardiotoxicity, safety pharmacology aims to determine the potential undesirable pharmacodynamic effects of a drug on the central nervous, vascular and respiratory systems (Pugsley et al., 2008). Thus, extension of the described approaches to these systems seems desirable, and might contribute to further advancement of these key areas of biomedical research.

Author Contributions

EG and SM wrote the editorial. EP and BR provided comments and edits.

Funding

Sources of support are the American Heart Association grant 15SDG24910015, the National Institutes of Health Stimulating Peripheral Activity to Relieve Conditions (SPARC) grant 1OT2OD023848-01, the National Heart, Lung, And Blood Institute (NHLBI) grants R01HL131517 and R01HL41214 (to EG); the NHLBI K99HL138160 award and the Heart Rhythm Society post-doctoral fellowship 16OA9HRS (to SM); projects DPI2016-75458-R (MINECO), ERC-2014-StG 638284 (European Research Council) and T39-17R (Aragón Government and European Regional Development Fund) (to EP); a Wellcome Trust Senior Research Fellowship in Basic Biomedical Sciences 100246/Z/12/Z and a NC3Rs Infrastructure for Impact Award (NC/P001076/1) (to BR).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Antzelevitch, C. (2005). Role of transmural dispersion of repolarization in the genesis of drug-induced torsades de pointes. Heart Rhythm 2, S9–S15. doi: 10.1016/j.hrthm.2004.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellinwood, N., Dobrev, D., Morotti, S., and Grandi, E. (2017). Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity. Chaos 27:093918. doi: 10.1063/1.5000226

CrossRef Full Text | Google Scholar

Food and Drug Administration, H. H. S. (2005). International Conference on Harmonisation; guidance on E14 Clinical Evaluation of QT/QTc Interval Prolongation and Proarrhythmic Potential for Non-Antiarrhythmic Drugs; availability. Notice. Fed. Regist. 70, 61134–61135.

Google Scholar

Gemmell, P., Burrage, K., Rodriguez, B., and Quinn, T. A. (2014). Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation. PLoS ONE 9:e90112. doi: 10.1371/journal.pone.0090112

PubMed Abstract | CrossRef Full Text | Google Scholar

Gintant, G., Sager, P. T., and Stockbridge, N. (2016). Evolution of strategies to improve preclinical cardiac safety testing. Nat. Rev. Drug Discov. 15, 457–471. doi: 10.1038/nrd.2015.34

PubMed Abstract | CrossRef Full Text | Google Scholar

Glukhov, A. V., Fedorov, V. V., Lou, Q., Ravikumar, V. K., Kalish, P. W., Schuessler, R. B., et al. (2010). Transmural dispersion of repolarization in failing and nonfailing human ventricle. Circ. Res. 106, 981–991. doi: 10.1161/CIRCRESAHA.109.204891

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, J. Q. X., and Sobie, E. A. (2018). Population-based mechanistic modeling allows for quantitative predictions of drug responses across cell types. NPJ Syst. Biol. Appl. 4:11. doi: 10.1038/s41540-018-0047-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, J. Q., Shim, J. V., Núñez-Acosta, E., and Sobie, E. A. (2017). I love it when a plan comes together: Insight gained through convergence of competing mathematical models. J. Mol. Cell. Cardiol. 102, 31–33. doi: 10.1016/j.yjmcc.2016.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Grandi, E., Pasqualini, F. S., and Bers, D. M. (2010). A novel computational model of the human ventricular action potential and Ca transient. J. Mol. Cell. Cardiol. 48, 112–121. doi: 10.1016/j.yjmcc.2009.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Kramer, J., Obejero-Paz, C. A., Myatt, G., Kuryshev, Y. A., Bruening-Wright, A., Verducci, J. S., et al. (2013). MICE models: superior to the HERG model in predicting Torsade de Pointes. Sci. Rep. 3:2100. doi: 10.1038/srep02100

PubMed Abstract | CrossRef Full Text | Google Scholar

Lancaster, M. C., and Sobie, E. A. (2016). Improved prediction of drug-induced torsades de pointes through simulations of dynamics and machine learning algorithms. Clin. Pharmacol. Ther. 100, 371–379. doi: 10.1002/cpt.367

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, Y. S., Liu, O. Z., Hwang, H. S., Knollmann, B. C., and Sobie, E. A. (2013). Parameter sensitivity analysis of stochastic models provides insights into cardiac calcium sparks. Biophys. J. 104, 1142–1150. doi: 10.1016/j.bpj.2012.12.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Mann, S. A., Imtiaz, M., Winbo, A., Rydberg, A., Perry, M. D., Couderc, J. P., et al. (2016). Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes. J. Mol. Cell. Cardiol. 100, 25–34. doi: 10.1016/j.yjmcc.2016.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirams, G. R., Cui, Y., Sher, A., Fink, M., Cooper, J., Heath, B. M., et al. (2011). Simulation of multiple ion channel block provides improved early prediction of compounds' clinical torsadogenic risk. Cardiovasc. Res. 91, 53–61. doi: 10.1093/cvr/cvr044

PubMed Abstract | CrossRef Full Text | Google Scholar

Morotti, S., and Grandi, E. (2017). Logistic regression analysis of populations of electrophysiological models to assess proarrythmic risk. MethodsX 4, 25–34. doi: 10.1016/j.mex.2016.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Morotti, S., McCulloch, A. D., Bers, D. M., Edwards, A. G., and Grandi, E. (2016). Atrial-selective targeting of arrhythmogenic phase-3 early afterdepolarizations in human myocytes. J. Mol. Cell. Cardiol. 96, 63–71. doi: 10.1016/j.yjmcc.2015.07.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Muszkiewicz, A., Liu, X., Bueno-Orovio, A., Lawson, B., Burrage, K., et al. (2017). From ionic to cellular variability in human atrial myocytes: an integrative computational and experimental study. Am. J. Physiol. Heart Circ. Physiol. 314, H895–H916. doi: 10.1152/ajpheart.00477.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Hara, T., Virág, L., Varro, A., and Rudy, Y. (2011). Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput. Biol. 7:e1002061. doi: 10.1371/journal.pcbi.1002061

PubMed Abstract | CrossRef Full Text | Google Scholar

Pueyo, E., Dangerfield, C. E., Britton, O. J., Virág, L., Kistamás, K., Szentandrássy, N., et al. (2016a). Experimentally-based computational investigation into beat-to-beat variability in ventricular repolarization and its response to ionic current inhibition. PLoS ONE 11:e0151461. doi: 10.1371/journal.pone.0151461

PubMed Abstract | CrossRef Full Text | Google Scholar

Pueyo, E., Orini, M., Rodríguez, J. F., and Taggart, P. (2016b). Interactive effect of beta-adrenergic stimulation and mechanical stretch on low-frequency oscillations of ventricular action potential duration in humans. J. Mol. Cell. Cardiol. 97, 93–105. doi: 10.1016/j.yjmcc.2016.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Pugsley, M. K., Authier, S., and Curtis, M. J. (2008). Principles of safety pharmacology. Br. J. Pharmacol. 154, 1382–1399. doi: 10.1038/bjp.2008.280

PubMed Abstract | CrossRef Full Text | Google Scholar

Redfern, W. S., Carlsson, L., Davis, A. S., Lynch, W. G., Mackenzie, I., Palethorpe, S., et al. (2003). Relationships between preclinical cardiac electrophysiology, clinical QT interval prolongation and torsade de pointes for a broad range of drugs: evidence for a provisional safety margin in drug development. Cardiovasc. Res. 58, 32–45. doi: 10.1016/S0008-6363(02)00846-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sager, P. T., Gintant, G., Turner, J. R., Pettit, S., and Stockbridge, N. (2014). Rechanneling the cardiac proarrhythmia safety paradigm: a meeting report from the Cardiac Safety Research Consortium. Am. Heart J. 167, 292–300. doi: 10.1016/j.ahj.2013.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Sánchez, C., Corrias, A., Bueno-Orovio, A., Davies, M., Swinton, J., Jacobson, I., et al. (2012). The Na+/K+ pump is an important modulator of refractoriness and rotor dynamics in human atrial tissue. Am. J. Physiol. Heart Circ. Physiol. 302, H1146–H1159. doi: 10.1152/ajpheart.00668.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarkar, A. X., and Sobie, E. A. (2011). Quantification of repolarization reserve to understand interpatient variability in the response to proarrhythmic drugs: a computational analysis. Heart Rhythm 8, 1749–1755. doi: 10.1016/j.hrthm.2011.05.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Ten Tusscher, K. H., and Panfilov, A. V. (2006). Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions. Phys. Med. Biol. 51, 6141–6156. doi: 10.1088/0031-9155/51/23/014

PubMed Abstract | CrossRef Full Text | Google Scholar

Vicente, J., Johannesen, L., Hosseini, M., Mason, J. W., Sager, P. T., Pueyo, E., et al. (2016). Electrocardiographic biomarkers for detection of drug-induced late sodium current block. PLoS ONE 11:e0163619. doi: 10.1371/journal.pone.0163619

PubMed Abstract | CrossRef Full Text | Google Scholar

Vorobyov, I., and Clancy, C. E. (2018). Sex, drugs, and funky rhythms. Heart Rhythm 15, 485–486. doi: 10.1016/j.hrthm.2018.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, P. C., Perissinotti, L. L., López-Redondo, F., Wang, Y., Demarco, K. R., Jeng, M. T., et al. (2017). A multiscale computational modelling approach predicts mechanisms of female sex risk in the setting of arousal-induced arrhythmias. J. Physiol. 595, 4695–4723. doi: 10.1113/JP273142

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: cardiotoxicity, QT interval prolongation, drug-induced arrhythmia, multi-scale modeling, cardiac electrophysiology

Citation: Grandi E, Morotti S, Pueyo E and Rodriguez B (2018) Editorial: Safety Pharmacology – Risk Assessment QT Interval Prolongation and Beyond. Front. Physiol. 9:678. doi: 10.3389/fphys.2018.00678

Received: 19 April 2018; Accepted: 15 May 2018;
Published: 08 June 2018.

Edited by:

Geoffrey A. Head, Baker Heart and Diabetes Institute, Australia

Reviewed by:

Arun V. Holden, University of Leeds, United Kingdom
Osmar Antonio Centurion, Universidad Nacional de Asunción, Paraguay

Copyright © 2018 Grandi, Morotti, Pueyo and Rodriguez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Eleonora Grandi, ele.grandi@gmail.com

In Memoriam: This Research Topic is dedicated to the memory of Dr. J. Jeremy Rice, our wonderful friend and colleague.