Sensitivity to switching rates in stochastically switched ODEs

We consider a stochastic process driven by a linear ordinary differential equation whose right-hand side switches at exponential times between a collection of different matrices. We construct planar examples that switch between two matrices where the individual matrices and the average of the two matrices are all Hurwitz (all eigenvalues have strictly negative real part), but nonetheless the process goes to infinity at large time for certain values of the switching rate. We further construct examples in higher dimensions where again the two individual matrices and their averages are all Hurwitz, but the process has arbitrarily many transitions between going to zero and going to infinity at large time as the switching rate varies. In order to construct these examples, we first prove in general that if each of the individual matrices is Hurwitz, then the process goes to zero at large time for sufficiently slow switching rate and if the average matrix is Hurwitz, then the process goes to zero at large time for sufficiently fast switching rate. We also give simple conditions that ensure the process goes to zero at large time for all switching rates. © 2014 International Press.

Biologically Inspired Design of Protein-Silica Hybrid Nanoparticles for Drug Delivery Applications

The design and application of effective drug carriers is a fundamental concern in the delivery of therapeutics for the treatment of cancer and other vexing health problems. Traditionally utilized chemotherapeutics are limited in efficacy due to poor bioavailability as a result of their size and solubility as well as significant deleterious effects to healthy tissue through their inability to preferentially target pathological cells and tissues, especially in treatment of cancer. Thus, a major effort in the development of nanoscopic drug delivery vehicles for cancer treatment has focused on exploiting the inherent differences in tumor physiology and limiting the exposure of drugs to non-tumorous tissue, which is commonly achieved by encapsulation of chemotherapeutics within macromolecular or supramolecular carriers that incorporate targeting ligands and that enable controlled release. The overall aim of this work is to engineer a hybrid nanomaterial system comprised of protein and silica and to characterize its potential as an encapsulating drug carrier. The synthesis of silica, an attractive nanomaterial component because it is both biocompatible as well as structurally and chemically stable, within this system is catalyzed by self-assembled elastin-like polypeptide (ELP) micelles that incorporate of a class of biologically-inspired, silica-promoting peptides, silaffins. Furthermore, this methodology produces near-monodisperse, hybrid inorganic/micellar materials under mild reaction conditions such as temperature, pH and solvent. This work studies this material system along three avenues: 1) proof-of-concept silicification (i.e. the formation and deposition of silica upon organic materials) of ELP micellar templates, 2) encapsulation and pH-triggered release of small, hydrophobic chemotherapeutics, and 3) selective silicification of templates to potentiate retention of peptide targeting ability.

The role of tectonic uplift, climate, and vegetation in the long-term terrestrial phosphorous cycle

Phosphorus (P) is a crucial element for life and therefore for maintaining ecosystem productivity. Its local availability to the terrestrial biosphere results from the interaction between climate, tectonic uplift, atmospheric transport, and biotic cycling. Here we present a mathematical model that describes the terrestrial P-cycle in a simple but comprehensive way. The resulting dynamical system can be solved analytically for steady-state conditions, allowing us to test the sensitivity of the P-availability to the key parameters and processes. Given constant inputs, we find that humid ecosystems exhibit lower P availability due to higher runoff and losses, and that tectonic uplift is a fundamental constraint. In particular, we find that in humid ecosystems the biotic cycling seem essential to maintain long-term P-availability. The time-dependent P dynamics for the Franz Josef and Hawaii chronosequences show how tectonic uplift is an important constraint on ecosystem productivity, while hydroclimatic conditions control the P-losses and speed towards steady-state. The model also helps describe how, with limited uplift and atmospheric input, as in the case of the Amazon Basin, ecosystems must rely on mechanisms that enhance P-availability and retention. Our novel model has a limited number of parameters and can be easily integrated into global climate models to provide a representation of the response of the terrestrial biosphere to global change. © 2010 Author(s).

Stochastic switching in infinite dimensions with applications to random parabolic PDE

© 2015 Society for Industrial and Applied Mathematics.We consider parabolic PDEs with randomly switching boundary conditions. In order to analyze these random PDEs, we consider more general stochastic hybrid systems and prove convergence to, and properties of, a stationary distribution. Applying these general results to the heat equation with randomly switching boundary conditions, we find explicit formulae for various statistics of the solution and obtain almost sure results about its regularity and structure. These results are of particular interest for biological applications as well as for their significant departure from behavior seen in PDEs forced by disparate Gaussian noise. Our general results also have applications to other types of stochastic hybrid systems, such as ODEs with randomly switching right-hand sides.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.

Fiber-optic interferometric two-dimensional scattering-measurement system.

We present a fiber-optic interferometric system for measuring depth-resolved scattering in two angular dimensions using Fourier-domain low-coherence interferometry. The system is a unique hybrid of the Michelson and Sagnac interferometer topologies. The collection arm of the interferometer is scanned in two dimensions to detect angular scattering from the sample, which can then be analyzed to determine the structure of the scatterers. A key feature of the system is the full control of polarization of both the illumination and the collection fields, allowing for polarization-sensitive detection, which is essential for two-dimensional angular measurements. System performance is demonstrated using a double-layer microsphere phantom. Experimental data from samples with different sizes and acquired with different polarizations show excellent agreement with Mie theory, producing structural measurements with subwavelength accuracy.

In vivo inhibition of elevated myocardial beta-adrenergic receptor kinase activity in hybrid transgenic mice restores normal beta-adrenergic signaling and function.

BACKGROUND: The clinical syndrome of heart failure (HF) is characterized by an impaired cardiac beta-adrenergic receptor (betaAR) system, which is critical in the regulation of myocardial function. Expression of the betaAR kinase (betaARK1), which phosphorylates and uncouples betaARs, is elevated in human HF; this likely contributes to the abnormal betaAR responsiveness that occurs with beta-agonist administration. We previously showed that transgenic mice with increased myocardial betaARK1 expression had impaired cardiac function in vivo and that inhibiting endogenous betaARK1 activity in the heart led to enhanced myocardial function. METHODS AND RESULTS: We created hybrid transgenic mice with cardiac-specific concomitant overexpression of both betaARK1 and an inhibitor of betaARK1 activity to study the feasibility and functional consequences of the inhibition of elevated betaARK1 activity similar to that present in human HF. Transgenic mice with myocardial overexpression of betaARK1 (3 to 5-fold) have a blunted in vivo contractile response to isoproterenol when compared with non-transgenic control mice. In the hybrid transgenic mice, although myocardial betaARK1 levels remained elevated due to transgene expression, in vitro betaARK1 activity returned to control levels and the percentage of betaARs in the high-affinity state increased to normal wild-type levels. Furthermore, the in vivo left ventricular contractile response to betaAR stimulation was restored to normal in the hybrid double-transgenic mice. CONCLUSIONS: Novel hybrid transgenic mice can be created with concomitant cardiac-specific overexpression of 2 independent transgenes with opposing actions. Elevated myocardial betaARK1 in transgenic mouse hearts (to levels seen in human HF) can be inhibited in vivo by a peptide that can prevent agonist-stimulated desensitization of cardiac betaARs. This may represent a novel strategy to improve myocardial function in the setting of compromised heart function.