Toward Elucidating the Physiological Impacts of Residual Stresses in the Colorectum.

Irritable bowel syndrome afflicts 10-20% of the global population, causing visceral pain with increased sensitivity to colorectal distension and normal bowel movements. Understanding and predicting these biomechanics will further advance our understanding of visceral pain and complement the existing literature on visceral neurophysiology. We recently performed a series of experiments at three longitudinal segments (colonic, intermediate, and rectal) of the distal 30 mm of colorectums of mice. We also established and fitted constitutive models addressing mechanical heterogeneity in both the through-thickness and longitudinal directions of the colorectum. Afferent nerve endings, strategically located within the submucosa, are likely nociceptors that detect concentrations of mechanical stresses to evoke the perception of pain from the viscera. In this study, we aim to: (1) establish and validate a method for incorporating residual stresses into models of colorectums, (2) predict the effects of residual stresses on the intratissue mechanics within the colorectum, and (3) establish intratissue distributions of stretches and stresses within the colorectum in vivo. To these ends we developed two-layered, composite finite element models of the colorectum based on our experimental evidence and validated our approaches against independent experimental data. We included layer- and segment-specific residual stretches/stresses in our simulations via the prestrain algorithm built into the finite element software febio. Our models and modeling approaches allow researchers to predict both organ and intratissue biomechanics of the colorectum and may facilitate better understanding of the underlying mechanical mechanisms of visceral pain.

Magnetic Resonance Imaging-based biomechanical simulation of cartilage: A systematic review.

MRI-based mathematical and computational modeling studies can contribute to a better understanding of the mechanisms governing cartilage's mechanical performance and cartilage disease. In addition, distinct modeling of cartilage is needed to optimize artificial cartilage production. These studies have opened up the prospect of further deepening our understanding of cartilage function. Furthermore, these studies reveal the initiation of an engineering-level approach to how cartilage disease affects material properties and cartilage function. Aimed at researchers in the field of MRI-based cartilage simulation, research articles pertinent to MRI-based cartilage modeling were identified, reviewed, and summarized systematically. Various MRI applications for cartilage modeling are highlighted, and the limitations of different constitutive models used are addressed. In addition, the clinical application of simulations and studied diseases are discussed. The paper's quality, based on the developed questionnaire, was assessed, and out of 79 reviewed papers, 34 papers were determined as high-quality. Due to the lack of the best constitutive models for various clinical conditions, researchers may consider the effect of constitutive material models on the cartilage disease simulation. In the future, research groups may incorporate various aspects of machine learning into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification, such as gait analysis.

Computational Modeling of Mouse Colorectum Capturing Longitudinal and Through-thickness Biomechanical Heterogeneity.

Mechanotransduction, the encoding of local mechanical stresses and strains at sensory endings into neural action potentials at the viscera, plays a critical role in evoking visceral pain, e.g., in the distal colon and rectum (colorectum). The wall of the colorectum is structurally heterogeneous, including two major composites: the inner consists of muscular and submucosal layers, and the outer consists of circular muscular, intermuscular, longitudinal muscular, and serosal layers. In fact the colorectum presents biomechanical heterogenity across both the longitudinal and through-thickness directions thus highlighting the differential roles of sensory nerve endings within different regions of the colorectum in visceral mechanotransduction. We determined constitutive models and model parameters for individual layers of the colorectum from three longitudinal locations (colonic, intermediate, and distal) using nonlinear optimization to fit our experimental results from biaxial extension tests on layer-separated colorectal tissues (mouse model, 7×7 mm2, Siri et al., Am. J. Physiol. Gastrointest. Liver Physiol. 316, G473-G481 and 317, G349-G358), and quantified the thicknesses of the layers. In this study we also quantified the residual stretches stemming from separating colorectal specimens into inner and outer composites and we completed new pressure-diameter mechanical testing to provide an additional validation case. We implemented the constitutive equations and created two-layered, 3-D finite element models using FEBio (University of Utah), and incorporated the residual stretches. We validated the modeling framework by comparing FE-predicted results for both biaxial extension testing of bulk specimens of colorectum and pressure-diameter testing of bulk segments against corresponding experimental results independent of those used in our model fitting. We present the first theoretical framework to simulate the biomechanics of distal colorectum, including both longitudinal and through-thickness heterogeneity, based on constitutive modeling of biaxial extension tests of colon tissues from mice. Our constitutive models and modeling framework facilitate analyses of both fundamental questions (e.g., the impact of organ/tissue biomechanics on mechanotransduction of the sensory nerve endings, structure-function relationships, and growth and remodeling in health and disease) and specific applications (e.g., device design, minimally invasive surgery, and biomedical research).

A multi-layered model of human skin elucidates mechanisms of wrinkling in the forehead.

Skin wrinkling, especially in the facial area, is a prominent sign of aging and is a growing area of research aimed at developing cosmetics and dermatological treatments. To better understand and treat undesirable skin wrinkles, it is vitally important to elucidate the underlying mechanisms of skin wrinkling, a largely mechanical process. Human skin, a multi-layer composite, has six mechanically distinct layers: from the outermost inward they are the stratum corneum, viable epidermis, dermal-epidermal-junction, papillary dermis, reticular dermis, and hypodermis. To better address the through-thickness hierarchy, and the development of wrinkling within this complicated hierarchy, we established a six-layered model of human skin realized with finite element modeling, by leveraging available morphological and biomechanical data on human skin of the forehead. Exercising our new model we aimed to quantify the effects of three potential mechanisms of wrinkle formation: (1) skin compression due to muscle contraction (dynamic wrinkles); (2) age-related volumetric tissue loss (static wrinkles); and (3) the combined effects of both mechanisms. Since hydration of the stratum corneum significantly affects its stiffness we also aimed to quantify the influence its hydration with these three potential mechanisms of wrinkle formation. Our six-layered skin model, combined with the proposed wrinkling mechanisms, successfully predicts the formation of dynamic and static wrinkles in the forehead consistent with the experimental literature. We observed three wrinkling modes in the forehead where the deepest wrinkles could reach to the reticular dermis. With further refinement our new six-layered model of human skin can be applied to study other region-specific wrinkle types such as the "crow's feet" and the nasolabial folds.

A multi-layered computational model for wrinkling of human skin predicts aging effects.

The development and progression of wrinkles from young to aged human skin relates to both structural and mechanical changes induced by aging. Here we aim to better understand the interaction of skin's layered morphology with dynamic wrinkles predicted in young and aged skin. First, we compare the predictions of wrinkling from 3-D finite element models of human skin including two to six distinct and anatomically motivated layers. Second, we perform parametric analyses using our six-layered model to determine how age-related changes in the architecture of human skin affect dynamic surface wrinkling. Specifically, we consider the following aging-related changes in the morphology of skin: flattening of the dermal-epidermal junction (DEJ) interface; thinning of both the viable epidermis (VE) and the reticular dermis (RD); and thickening of the papillary dermis (PD). We use skin compression to model dynamic, expressional wrinkles due to muscle contraction, and volumetric tissue loss to model effects of aging in wrinkling simulations. Our results highlight the role of skin's multi-layered structure in the modeling of wrinkling formation. Our six-layered model, consisting of all of the mechanical layers, predicts deep wrinkles with better fidelity than models including fewer layers. From our parametric study, applying our six-layered model, we conclude that: (1) the relative thicknesses of the layers in the epidermis or dermis significantly influences surface wrinkling in skin; and, (2) flattening of the DEJ with aging enhances surface wrinkling. Thinning of VE increases the relative stiffness of the epidermis and thus enhances dynamic wrinkling, while thickening of PD or thinning of RD has the same effect by reducing the equivalent stiffness of the substrate. Consequently, strategies to minimize wrinkling could maintain the undulating morphology of the DEJ, thereby delaying dynamic wrinkling and delaying the propagation of buckling into the deeper dermis or hypodermis. Additional strategies to minimize wrinkling could target preventing the VE and RD from thinning or preventing the PD from thickening.

Chondrocyte viability is lost during high-rate impact loading by transfer of amplified strain, but not stress, to pericellular and cellular regions.

OBJECTIVE: Deleterious impact loading to cartilage initiates post-traumatic osteoarthritis (OA). While cytokine and enzyme levels regulate disease progression, specific mechanical cues that elucidate cellular OA origins merit further investigation. We defined the dominant pericellular and cellular strain/stress transfer mechanisms following bulk-tissue injury associated with cell death. METHOD: Using an in vitro model, we investigated rate-dependent loading and spatial localization of cell viability in acute indentation and time-course studies. Atomic force microscopy (AFM) and magnetic resonance imaging (MRI) confirmed depth-wise changes in cartilage micro-/macro-mechanics and structure post-indentation. To understand the transfer of loading to cartilage domains, we computationally modeled full-field strain and stress measures in interstitial matrix, pericellular and cellular regions. RESULTS: Chondrocyte viability decreased following rapid impact (80%/s) vs slow loading (0.1%/s) or unloaded controls. Viability was lost immediately during impact within regions near the indenter-tissue contact but did not change over 7 days of tissue culture. AFM studies revealed a loss of stiffness following 80%/s loading, and MRI studies confirmed an increased tensile and shear strain, but not relaxometry. Image-based patterns of chondrocyte viability closely matched computational estimates of amplified maximum principal and shear strain in interstitial matrix, pericellular and cellular regions. CONCLUSION: Rapid indentation worsens chondrocyte death and degrades cartilage matrix stiffness in indentation regions. Cell death at high strain rates may be driven by elevated tensile strains, but not matrix stress. Strain amplification beyond critical thresholds in the pericellular matrix and cells may define a point of origin for early damage in post-traumatic OA.

Propagation of microcracks in collagen networks of cartilage under mechanical loads.

OBJECTIVE: We recently demonstrated that low-energy mechanical impact to articular cartilage, usually considered non-injurious, can in fact cause microscale cracks (widths <30µm) in the collagen network of visually pristine human cartilage. While research on macro-scale cracks in cartilage and microcracks in bone abounds, how microcracks within cartilage initiate and propagate remains unknown. We quantified the extent to which microcracks initiate and propagate in the collagen network during mechanical loading representative of normal activities. DESIGN: We tested 76 full-thickness, cylindrical osteochondral plugs. We imaged untreated specimens (pristine phase) via second harmonic generation and assigned specimens to three low-energy impact groups (none, low, high), and thereafter to three cyclic compression groups (none, low, high) which simulate walking. We re-imaged specimens in the post-impact and post-cyclic compression phases to identify and track microcracks. RESULTS: Microcracks in the network of collagen did not present in untreated controls but did initiate and propagate under mechanical treatments. We found that the length and width of microcracks increased from post-impact to post-cyclic compression in tracked microcracks, but neither depth nor angle presented statistically significant differences. CONCLUSIONS: The microcracks we initiated under low-energy impact loading increased in length and width during subsequent cyclic compression that simulated walking. The extent of this propagation depended on the combination of impact and cyclic compression. More broadly, the initiation and propagation of microcracks may characterize pathogenesis of osteoarthritis, and may suggest therapeutic targets for future studies.

Through-thickness patterns of shear strain evolve in early osteoarthritis.

OBJECTIVE: Given the structural changes associated with the progression of Osteoarthritis (OA), we hypothesized that patterns of through-thickness, large-strain shear evolve with early-stage OA. We therefore aimed to determine whether and how patterns of shear strains change during early-stage OA to 1) gain insight into the progression of OA by quantifying changes in local deformations; 2) gauge the potential of patterns in shear strain to serve as image-based biomarkers of early-stage OA; and 3) provide high-resolution, through-thickness data for proposing, fitting, and validating constitutive models for cartilage. DESIGN: We completed displacement-driven, large-strain shear tests (5, 10, 15%) on 44 specimens of variably advanced osteoarthritic human articular cartilage as determined by both Osteoarthritis Research Society International (OARSI) grade and PLM-CO score. We recorded the through-thickness deformations with a stereo-camera system and processed these data using digital image correlation (DIC) to determine full-thickness patterns of strains and relative zonal recruitments, i.e., the average shear strain in a through-thickness zone weighted by its relative thickness and normalized by the applied strain. RESULTS: We observed three general shapes for the curves of averaged through-thickness, Green-Lagrange shear strains during progression of OA. We also observed that during the progression of OA only the deep zone is recruited differently under shear in a statistically significant way. CONCLUSIONS: We propose that changes in through-thickness patterns of shear strain could provide sensitive biomarkers for early clinical detection of OA. The relative zonal recruitment of the deep zone decreases with progressing OA (OARSI grade) and microstructural remodeling (PLM-CO score), which do not consistently affect recruitment of the superficial and middle zones.

The evolving large-strain shear responses of progressively osteoarthritic human cartilage.

OBJECTIVE: The composition and structure of articular cartilage evolves during the development and progression of osteoarthritis (OA) resulting in changing mechanical responses. We aimed to assess the evolution of the intrinsic, large-strain mechanics of human articular cartilage-governed by collagen and proteoglycan and their interactions-during the progression of OA. DESIGN: We completed quasi-static, large-strain shear tests on 64 specimens from ten donors undergoing total knee arthroplasty (TKA), and quantified the corresponding state of OA (OARSI grade), structural integrity (PLM score), and composition (glycosaminoglycan and collagen content). RESULTS: We observed nonlinear stress-strain relationships with distinct hystereses for all magnitudes of applied strain where stiffnesses, nonlinearities, and hystereses all reduced as OA advanced. We found a reduction in energy dissipation density up to 80% in severely degenerated (OARSI grade 4, OA-4) vs normal (OA-1) cartilage, and more importantly, we found that even cartilage with a normal appearance in structure and composition (OA-1) dissipated 50% less energy than healthy (control) load-bearing cartilage (HL0). Changes in stresses and stiffnesses were in general less pronounced and did not allow us to distinguish between healthy load-bearing controls and very early-stage OA (OA-1), or to distinguish consistently among different levels of degeneration, i.e., OARSI grades. CONCLUSIONS: Our results suggest that reductions in energy dissipation density can be detected by bulk-tissue testing, and that these reductions precede visible signs of degeneration. We highlight the potential of energy dissipation, as opposed to stress- or stiffness-based measures, as a marker to diagnose early-stage OA.

GPGPU-based explicit finite element computations for applications in biomechanics: the performance of material models, element technologies, and hardware generations.

Finite element (FE) simulations are increasingly valuable in assessing and improving the performance of biomedical devices and procedures. Due to high computational demands such simulations may become difficult or even infeasible, especially when considering nearly incompressible and anisotropic material models prevalent in analyses of soft tissues. Implementations of GPGPU-based explicit FEs predominantly cover isotropic materials, e.g. the neo-Hookean model. To elucidate the computational expense of anisotropic materials, we implement the Gasser-Ogden-Holzapfel dispersed, fiber-reinforced model and compare solution times against the neo-Hookean model. Implementations of GPGPU-based explicit FEs conventionally rely on single-point (under) integration. To elucidate the expense of full and selective-reduced integration (more reliable) we implement both and compare corresponding solution times against those generated using underintegration. To better understand the advancement of hardware, we compare results generated using representative Nvidia GPGPUs from three recent generations: Fermi (C2075), Kepler (K20c), and Maxwell (GTX980). We explore scaling by solving the same boundary value problem (an extension-inflation test on a segment of human aorta) with progressively larger FE meshes. Our results demonstrate substantial improvements in simulation speeds relative to two benchmark FE codes (up to 300[Formula: see text] while maintaining accuracy), and thus open many avenues to novel applications in biomechanics and medicine.

Automated hexahedral meshing of knee cartilage structures - application to data from the osteoarthritis initiative.

We propose a fully automated methodology for hexahedral meshing of patient-specific structures of the human knee obtained from magnetic resonance images, i.e. femoral/tibial cartilages and menisci. We select eight patients from the Osteoarthritis Initiative and validate our methodology using MATLAB on a laptop computer. We obtain the patient-specific meshes in an average of three minutes, while faithfully representing the geometries with well-shaped elements. We hope to provide a fundamentally different means to test hypotheses on the mechanisms of disease progression by integrating our patient-specific FE meshes with data from individual patients. Download both our meshes and software at http://im.engr.uconn.edu/downloads.php .

Rupture risk in abdominal aortic aneurysms: A realistic assessment of the explicit GPU approach.

Accurate estimation of peak wall stress (PWS) is the crux of biomechanically motivated rupture risk assessment for abdominal aortic aneurysms aimed to improve clinical outcomes. Such assessments often use the finite element (FE) method to obtain PWS, albeit at a high computational cost, motivating simplifications in material or element formulations. These simplifications, while useful, come at a cost of reliability and accuracy. We achieve research-standard accuracy and maintain clinically applicable speeds by using novel computational technologies. We present a solution using our custom finite element code based on graphics processing unit (GPU) technology that is able to account for added complexities involved with more physiologically relevant solutions, e.g. strong anisotropy and heterogeneity. We present solutions up to 17× faster relative to an established finite element code using state-of-the-art nonlinear, anisotropic and nearly-incompressible material descriptions. We show a realistic assessment of the explicit GPU FE approach by using complex problem geometry, biofidelic material law, double-precision floating point computation and full element integration. Due to the increased solution speed without loss of accuracy, shown on five clinical cases of abdominal aortic aneurysms, the method shows promise for clinical use in determining rupture risk of abdominal aortic aneurysms.

Low-energy impact of human cartilage: predictors for microcracking the network of collagen.

OBJECTIVE: We aimed to determine the minimum mechanical impact to cause microstructural damage in the network of collagen (microcracking) within human cartilage and hypothesized that energies below 0.1 J or 1 mJ/mm3 would suffice. DESIGN: We completed 108 low-energy impact tests (0.05, 0.07, or 0.09 J; 0.75 or 1.0 m/s2) using healthy cartilage specimens from six male donors (30.2 ± 8.8 yrs old). Before and after impact we acquired, imaging the second harmonic generation (SHG), ten images from each specimen (50 µm depth, 5 µm step size), resulting in 2160 images. We quantified both the presence and morphology of microcracks. We then correlated test parameters (predictors) impact energy/energy dissipation density, nominal stress/stress rate, and strain/strain rate to microcracking and tested for significance. Where predictors significantly correlated with microstructural outcomes we fitted binary logistic regression plots with 95% confidence intervals (CIs). RESULTS: No specimens presented visible damage following impact. We found that impact energy/energy dissipation density and nominal stress/stress rate were significant (P < 0.05) predictors of microcracking while both strain and strain rate were not. In our test configuration, an impact energy density of 2.93 mJ/mm3, an energy dissipation density of 1.68 mJ/mm3, a nominal stress of 4.18 MPa, and a nominal stress rate of 689 MPa/s all corresponded to a 50% probability of microcracking in the network of collagen. CONCLUSIONS: An impact energy density of 1.0 mJ/mm3 corresponded to a ∼20% probability of microcracking. Such changes may initiate a degenerative cascade leading to post-traumatic osteoarthritis.

High-quality conforming hexahedral meshes of patient-specific abdominal aortic aneurysms including their intraluminal thrombi.

In order to perform finite element (FE) analyses of patient-specific abdominal aortic aneurysms, geometries derived from medical images must be meshed with suitable elements. We propose a semi-automatic method for generating conforming hexahedral meshes directly from contours segmented from medical images. Magnetic resonance images are generated using a protocol developed to give the abdominal aorta high contrast against the surrounding soft tissue. These data allow us to distinguish between the different structures of interest. We build novel quadrilateral meshes for each surface of the sectioned geometry and generate conforming hexahedral meshes by combining the quadrilateral meshes. The three-layered morphology of both the arterial wall and thrombus is incorporated using parameters determined from experiments. We demonstrate the quality of our patient-specific meshes using the element Scaled Jacobian. The method efficiently generates high-quality elements suitable for FE analysis, even in the bifurcation region of the aorta into the iliac arteries. For example, hexahedral meshes of up to 125,000 elements are generated in less than 130 s, with 94.8 % of elements well suited for FE analysis. We provide novel input for simulations by independently meshing both the arterial wall and intraluminal thrombus of the aneurysm, and their respective layered morphologies.

Modeling sample/patient-specific structural and diffusional responses of cartilage using DT-MRI.

We propose a new 3D biphasic constitutive model designed to incorporate structural data on the sample/patient-specific collagen fiber network. The finite strain model focuses on the load-bearing morphology, that is, an incompressible, poroelastic solid matrix, reinforced by an inhomogeneous, dispersed fiber fabric, saturated with an incompressible fluid at constant electrolytic conditions residing in strain-dependent pores of the collagen-proteoglycan solid matrix. In addition, the fiber network of the solid influences the fluid permeability and an intrafibrillar portion that cannot be 'squeezed out' from the tissue. We implement the model into a finite element code. To demonstrate the utility of our proposed modeling approach, we test two hypotheses by simulating an indentation experiment for a human tissue sample. The simulations use ultra-high field diffusion tensor magnetic resonance imaging that was performed on the tissue sample. We test the following hypotheses: (i) the through-thickness structural arrangement of the collagen fiber network adjusts fluid permeation to maintain fluid pressure (Biomech. Model. Mechanobiol. 7:367-378, 2008); and (ii) the inhomogeneity of mechanical properties through the cartilage thickness acts to maintain fluid pressure at the articular surface (J. Biomech. Eng. 125:569-577, 2003). For the tissue sample investigated, both through-thickness inhomogeneities of the collagen fiber distribution and of the material properties serve to influence the interstitial fluid pressure distribution and maintain fluid pressure underneath the indenter at the cartilage surface. Tissue inhomogeneity appears to have a larger effect on fluid pressure retention in this tissue sample and on the advantageous pressure distribution.

Remodeling of intramural thrombus and collagen in an Ang-II infusion ApoE-/- model of dissecting aortic aneurysms.

Fibrillar collagen endows the normal aortic wall with significant stiffness and strength and similarly plays important roles in many disease processes. For example, because of the marked loss of elastic fibers and functional smooth cells in aortic aneurysms, collagen plays a particularly important role in controlling the dilatation of these lesions and governing their rupture potential. Recent findings suggest further that collagen remodeling may also be fundamental to the intramural healing of arterial or aneurysmal dissections. To explore this possibility further, we identified and correlated regions of intramural thrombus and newly synthesized fibrillar collagen in a well-established mouse model of dissecting aortic aneurysms. Our findings suggest that intramural thrombus that is isolated from free-flowing blood creates a permissive environment for the synthesis of fibrillar collagen that, albeit initially less dense and organized, could protect that region of the dissected wall from subsequent expansion of the dissection or rupture. Moreover, alpha-smooth muscle actin positive cells appeared to be responsible for the newly produced collagen, which co-localized with significant production of glycosaminoglycans.

Granisetron as antiemetic therapy in children with cancer.

The safety and efficacy of the new 5HT-3 antagonist granisetron as an antiemetic in children with cancer was evaluated in 40 children at a single dose of 40 micrograms/kg. No adverse affects attributable to the granisetron were noted. The overall major and complete response rate was 82.5% and this was highest in the younger children. Only 2 patients showed no response. Pharmacokinetic studies showed associations between some pharmacokinetic parameters and age which were no longer apparent after normalisation for body weight. Granisetron is an effective and very well-tolerated antiemetic and appears to be an important addition to the supportive care available for children with cancer.

Pharmacokinetics and tolerability of ascending intravenous doses of granisetron, a novel 5-HT3 antagonist, in healthy human subjects.

The pharmacokinetics and tolerance of granisetron, a novel 5HT3-receptor antagonist which is under development as an anti-emetic agent have been studied after administration of single 30 min intravenous infusions to three groups of 8 healthy male subjects, in a series of placebo-controlled ascending dose studies (50, 80, 100 and 130 micrograms.kg-1 to group 1; 150, 180, 200 and 230 micrograms.kg-1 to group 2 and 270 and 300 micrograms.kg-1 to group 3). Plasma and urine samples were analysed for granisetron by HPLC with fluorimetric detection. Administration of granisetron was well tolerated by the volunteers and there were no serious adverse effects reported. Pharmacokinetic parameters and dose-normalised plasma levels appeared to be independent of dose in the range 50 to 300 micrograms.kg-1, although there was extensive inter-subject variability. Granisetron was extensively distributed, with mean volumes of distribution ranging from 186-264 l at the various doses. Total plasma clearance was, in general, rapid (mean values of 37.0 to 49.9 l.h-1) and predominantly non-renal, with most subjects excreting less than 20% of the dose unchanged in urine. Mean t1/2 values ranged from 4.1 to 6.3 h and MRT from 5.2 to 8.1 h.

The tolerance to and pharmacokinetics of penciclovir (BRL 39,123A), a novel antiherpes agent, administered by intravenous infusion to healthy subjects.

The tolerance to and pharmacokinetics of intravenously administered penciclovir (BRL 39,123A), a novel anti-herpes agent, were investigated in 15 healthy male subjects. The volunteers were divided into three groups, receiving either 10, 15 or 20 mg/kg penciclovir by a 60 min constant-rate infusion. Blood samples were taken sequentially up to 48 h after the start of the infusion and urine collections made at appropriate intervals up to 72 h. After a simple solid phase extraction, concentrations of penciclovir in plasma and urine were determined using HPLC with U.V. detection. Mean values of Cmax, corresponding usually with the end of infusion, and of AUC appeared to increase proportionately with dose. Furthermore, there was no evidence that dose significantly affected any individual pharmacokinetic parameter. Penciclovir was distributed into tissues with an overall mean volume of distribution of approximately 1.5 l.kg-1, i.e. approximately double that of body water. It was rapidly eliminated, with a mean total plasma clearance of 39.3 l.h-1, and a mean terminal-phase half-life of 2.0 h. The majority of the dose, approximately 70%, was excreted unchanged in the urine. Mean renal clearance of BRL 39,123 was 28.1 l.h-1, which exceeds normal glomerular filtration rate and approaches renal plasma flow. At all dose-levels, the infusions of penciclovir were well tolerated, with no evidence of drug-related adverse events.