A Mechanistic Latent Variable Model for Estimating Drug Concentrations in the Male Genital Tract

The purpose of this study is to develop statistical methodology to facilitate indirect estimation of the concentration of antiretroviral drugs and viral loads in the prostate gland and the seminal vesicle. The differences in antiretroviral drug concentrations in these organs may lead to suboptimal concentrations in one gland compared to the other. Suboptimal levels of the antiretroviral drugs will not be able to fully suppress the virus in that gland, lead to a source of sexually transmissible virus and increase the chance of selecting for drug resistant virus. This information may be useful selecting antiretroviral drug regimen that will achieve optimal concentrations in most of male genital tract glands. Using fractionally collected semen ejaculates, Lundquist (1949) measured levels of surrogate markers in each fraction that are uniquely produced by specific male accessory glands. To determine the original glandular concentrations of the surrogate markers, Lundquist solved a simultaneous series of linear equations. This method has several limitations. In particular, it does not yield a unique solution, it does not address measurement error, and it disregards inter-subject variability in the parameters. To cope with these limitations, we developed a mechanistic latent variable model based on the physiology of the male genital tract and surrogate markers. We employ a Bayesian approach and perform a sensitivity analysis with regard to the distributional assumptions on the random effects and priors. The model and Bayesian approach is validated on experimental data where the concentration of a drug should be (biologically) differentially distributed between the two glands. In this example, the Bayesian model-based conclusions are found to be robust to model specification and this hierarchical approach leads to more scientifically valid conclusions than the original methodology. In particular, unlike existing methods, the proposed model based approach was not affected by a common form of outliers.

ADAMTS-13, von Willebrand factor and related parameters in severe sepsis and septic shock

BACKGROUND: Insufficient control of von Willebrand factor (VWF) multimer size as a result of severely deficient ADAMTS-13 activity results in thrombotic thrombocytopenic purpura associated with microvascluar thrombosis and platelet consumption, features not seldom seen in severe sepsis and septic shock. METHODS: ADAMTS-13 activity and VWF parameters of 40 patients with severe sepsis or septic shock were compared with those of 40 healthy controls of the same age and gender and correlated with clinical findings and sepsis outcome. RESULTS: ADAMTS-13 activity was significantly lower in patients than in healthy controls [median 60% (range 27-160%) vs. 110% (range 63-200%); P < 0.001]. VWF parameters behaved reciprocally and both VWF ristocetin cofactor activity (RCo) and VWF antigen (VWF:Ag) were significantly (P < 0.001) higher in patients compared with controls. Neither ADAMTS-13 activity nor VWF parameters correlated with disease severity, organ dysfunction or outcome. However, a contribution of acute endothelial dysfunction to renal impairment in sepsis is suggested by the significantly higher VWF propeptide and soluble thrombomodulin levels in patients with increased creatinine values as well as by their strong positive correlations (creatinine and VWF propeptide r(s) = 0.484, P < 0.001; creatinine and soluble thrombomodulin r(s) = 0.596, P < 0.001). CONCLUSIONS: VWF parameters are reciprocally correlated with ADAMTS-13 activity in severe sepsis and septic shock but have no prognostic value regarding outcome.

Precise pose recovery of distal locking holes from single calibrated fluoroscopic image via a novel variable decomposition approach

This paper presents a novel variable decomposition approach for pose recovery of the distal locking holes using single calibrated fluoroscopic image. The problem is formulated as a model-based optimal fitting process, where the control variables are decomposed into two sets: (a) the angle between the nail axis and its projection on the imaging plane, and (b) the translation and rotation of the geometrical model of the distal locking hole around the nail axis. By using an iterative algorithm to find the optimal values of the latter set of variables for any given value of the former variable, we reduce the multiple-dimensional model-based optimal fitting problem to a one-dimensional search along a finite interval. We report the results of our in vitro experiments, which demonstrate that the accuracy of our approach is adequate for successful distal locking of intramedullary nails.

Effect of X-ray tube parameters, iodine concentration, and patient size on image quality in pulmonary computed tomography angiography: a chest-phantom-study

OBJECTIVES: The aim of this phantom study was to evaluate the contrast-to-noise ratio (CNR) in pulmonary computed tomography (CT)-angiography for 300 and 400 mg iodine/mL contrast media using variable x-ray tube parameters and patient sizes. We also analyzed the possible strategies of dose reduction in patients with different sizes. MATERIALS AND METHODS: The segmental pulmonary arteries were simulated by plastic tubes filled with 1:30 diluted solutions of 300 and 400 mg iodine/mL contrast media in a chest phantom mimicking thick, intermediate, and thin patients. Volume scanning was done with a CT scanner at 80, 100, 120, and 140 kVp. Tube current-time products (mAs) varied between 50 and 120% of the optimal value given by the built-in automatic dose optimization protocol. Attenuation values and CNR for both contrast media were evaluated and compared with the volume CT dose index (CTDI(vol)). Figure of merit, calculated as CNR/CTDIvol, was used to quantify image quality improvement per exposure risk to the patient. RESULTS: Attenuation of iodinated contrast media increased both with decreasing tube voltage and patient size. A CTDIvol reduction by 44% was achieved in the thin phantom with the use of 80 instead of 140 kVp without deterioration of CNR. Figure of merit correlated with kVp in the thin phantom (r = -0.897 to -0.999; P < 0.05) but not in the intermediate and thick phantoms (P = 0.09-0.71), reflecting a decreasing benefit of tube voltage reduction on image quality as the thickness of the phantom increased. Compared with the 300 mg iodine/mL concentration, the same CNR for 400 mg iodine/mL contrast medium was achieved at a lower CTDIvol by 18 to 40%, depending on phantom size and applied tube voltage. CONCLUSIONS: Low kVp protocols for pulmonary embolism are potentially advantageous especially in thin and, to a lesser extent, in intermediate patients. Thin patients profit from low voltage protocols preserving a good CNR at a lower exposure. The use of 80 kVp in obese patients may be problematic because of the limitation of the tube current available, reduced CNR, and high skin dose. The high CNR of the 400 mg iodine/mL contrast medium together with lower tube energy and/or current can be used for exposure reduction.

Nuclear localization of epidermal growth factor and epidermal growth factor receptors in human thyroid tissues

Epidermal growth factor (EGF) has widespread growth effects, and in some tissues proliferation is associated with the nuclear localization of EGF and epidermal growth factor receptor (EGFR). In the thyroid, EGF promotes growth but differs from thyrotropin (TSH) in inhibiting rather than stimulating functional parameters. We have therefore studied the occurrence and cellular distribution of EGF and EGFR in normal thyroid, in Graves' disease, where growth is mediated through the thyrotropin receptor (TSHR), and in a variety of human thyroid tumors. In the normal gland the staining was variable, but largely cytoplasmic, for both EGF and EGFR. In Graves' disease there was strong cytoplasmic staining for both EGF and EGFR, with frequent positive nuclei. Nuclear positivity for EGF and particularly for EGFR was also a feature of both follicular adenomas and follicular carcinomas. Interestingly, nuclear staining was almost absent in papillary carcinomas. These findings document for the first time the presence of nuclear EGF and EGFR in thyroid. Their predominant occurrence in tissues with increased growth (Graves' disease, follicular adenoma, and carcinoma) may indicate that nuclear EGF and EGFR play a role in growth regulation in these conditions. The absence of nuclear EGF and EGFR in papillary carcinomas would suggest that the role played by EGF in growth control differs between papillary carcinoma and follicular adenomas/carcinomas of the thyroid.

Implant Failure Predictors in the Posterior Maxilla: A Retrospective Study of 273 Consecutive Implants

Background: The goal of this study was to retrospectively analyze a cohort of 136 patients who underwent dental implant placement in the posterior maxilla at the University of Connecticut Health Center to assess and identify predictors for implant failure in the posterior maxilla. Methods: Data were retrieved from patient charts to identify subjects older than 21 years of age who received dental implant(s) in the posterior maxilla. Patients without a postoperative baseline radiograph were excluded. A recall radiograph was taken 3 to 6 months after implant placement. If there was no recall radiograph, the subject was contacted for a recall visit that included a clinical evaluation and radiographs to determine the implant status. Based on a univariate screening, variables considered potential implant failure predictors included gender, diabetes, smoking, implant length, implant diameter, membrane use, sinus-elevation technique, and surgical complications. These parameters were further assessed, and a multivariable logistic regression was performed with implant failure as a dependant variable. All tests of significance were evaluated at the 0.05 error level. Results: Two hundred seventy-three implants were placed in the posterior maxilla. Fourteen implants failed (early and late failures combined), resulting in a 94.9% overall survival rate. The survival rates for the sinus-elevation group and native bone group were 92.2% and 96.7%, respectively (P = 0.090). Based on the multivariable analysis, sinus floor-elevation procedures were not associated with increased risk for implant failure (P = 0.702). In contrast, smoking and surgical complications had a statistically significant effect on implant failure; the odds ratios for implant failure were 6.4 (P = 0.025) and 8.2 (P = 0.004), respectively. Conclusion: Sinus-elevation procedures with simultaneous or staged implant placement do not increase the risk for implant failure, whereas smoking and surgical complications markedly increase the risk for implant failure.

Recurrent rearrangements of chromosome 1q21.1 and variable pediatric phenotypes

BACKGROUND: Duplications and deletions in the human genome can cause disease or predispose persons to disease. Advances in technologies to detect these changes allow for the routine identification of submicroscopic imbalances in large numbers of patients. METHODS: We tested for the presence of microdeletions and microduplications at a specific region of chromosome 1q21.1 in two groups of patients with unexplained mental retardation, autism, or congenital anomalies and in unaffected persons. RESULTS: We identified 25 persons with a recurrent 1.35-Mb deletion within 1q21.1 from screening 5218 patients. The microdeletions had arisen de novo in eight patients, were inherited from a mildly affected parent in three patients, were inherited from an apparently unaffected parent in six patients, and were of unknown inheritance in eight patients. The deletion was absent in a series of 4737 control persons (P=1.1x10(-7)). We found considerable variability in the level of phenotypic expression of the microdeletion; phenotypes included mild-to-moderate mental retardation, microcephaly, cardiac abnormalities, and cataracts. The reciprocal duplication was enriched in nine children with mental retardation or autism spectrum disorder and other variable features (P=0.02). We identified three deletions and three duplications of the 1q21.1 region in an independent sample of 788 patients with mental retardation and congenital anomalies. CONCLUSIONS: We have identified recurrent molecular lesions that elude syndromic classification and whose disease manifestations must be considered in a broader context of development as opposed to being assigned to a specific disease. Clinical diagnosis in patients with these lesions may be most readily achieved on the basis of genotype rather than phenotype.

Complexity of chronic asthma and chronic obstructive pulmonary disease: implications for risk assessment, and disease progression and control

Although assessment of asthma control is important to guide treatment, it is difficult since the temporal pattern and risk of exacerbations are often unpredictable. In this Review, we summarise the classic methods to assess control with unidimensional and multidimensional approaches. Next, we show how ideas from the science of complexity can explain the seemingly unpredictable nature of bronchial asthma and emphysema, with implications for chronic obstructive pulmonary disease. We show that fluctuation analysis, a method used in statistical physics, can be used to gain insight into asthma as a dynamic disease of the respiratory system, viewed as a set of interacting subsystems (eg, inflammatory, immunological, and mechanical). The basis of the fluctuation analysis methods is the quantification of the long-term temporal history of lung function parameters. We summarise how this analysis can be used to assess the risk of future asthma episodes, with implications for asthma severity and control both in children and adults.

Electron transport in E x B devices

A Hall thruster, an E × B device used for in-space propulsion, utilizes an axial electric field to electrostatically accelerate plasma propellant from the spacecraft. The axial electric field is created by positively biasing the anode so that the positivelycharged ions may be accelerated (repelled) from the thruster, which produces thrust. However, plasma electrons are much smaller than ions and may be accelerated much more quickly toward the anode; if electrons were not impeded, a "short circuit" due to the electron flow would eliminate the thrust mechanism. Therefore, a magnetic field serves to "magnetize" plasma electrons internal to the thruster and confines them in gyro-orbits within the discharge channel. Without outside factors electrons would be confined indefinitely; however, electron-neutral collisions provide a mechanism to free electrons from their orbits allowing electrons to cross the magnetic field toward the anode, where this process is described by classical transport theory. To make matters worse, cross-field electron transport has been observed to be 100-1000 times that predicted by classical collisional theory, providing an efficiency loss mechanism and an obstacle for modeling and simulations in Hall thrusters. The main difficulty in studying electron transport in Hall thrusters is the coupling that exists between the plasma and the fields, where the plasma creates and yet is influenced by the electric field. A device has been constructed at MTU’s Isp Lab, the Hall Electron Mobility Gage, which was designed specifically to study electron transport in E × B devices, where the coupling between the plasma and electric field was virtually eliminated. In this device the two most cited contributors to electron transport in Hall thrusters, fluctuation-induced transport, and wall effects, were absent. Removing the dielectric walls and plasma fluctuations, while maintaining the field environment in vacuum, has allowed the study of electron dynamics in Hall thruster fields where the electrons behave as test particles in prescribed fields, greatly simplifying the environment. Therefore, it was possible to observe any effects on transport not linked to the cited mechanisms, and it was possible to observe trends of the enhanced mobility with control parameters of electric and magnetic fields and neutral density– parameters that are not independently variable in a Hall thruster. The result of the investigation was the observation of electron transport that was ~ 20-100 times the classical prediction. The cross-field electron transport in the Mobility Gage was generally lower than that found in a Hall thruster so these findings do not negate the possibility of fluctuations and/or wall collisions contributing to transport in a Hall thruster. However, this research led to the observation of enhanced cross-field transport that had not been previously isolated in Hall thruster fields, which is not reliant on momentum-transfer collisions, wall collisions or fluctuations.

Dynamic modeling, simulation and control design of a parafoil-payload system for ship launched aerial delivery system (SLADS)

The objective of this research was to develop a high-fidelity dynamic model of a parafoilpayload system with respect to its application for the Ship Launched Aerial Delivery System (SLADS). SLADS is a concept in which cargo can be transfered from ship to shore using a parafoil-payload system. It is accomplished in two phases: An initial towing phase when the glider follows the towing vessel in a passive lift mode and an autonomous gliding phase when the system is guided to the desired point. While many previous researchers have analyzed the parafoil-payload system when it is released from another airborne vehicle, limited work has been done in the area of towing up the system from ground or sea. One of the main contributions of this research was the development of a nonlinear dynamic model of a towed parafoil-payload system. After performing an extensive literature review of the existing methods of modeling a parafoil-payload system, a five degree-of-freedom model was developed. The inertial and geometric properties of the system were investigated to predict accurate results in the simulation environment. Since extensive research has been done in determining the aerodynamic characteristics of a paraglider, an existing aerodynamic model was chosen to incorporate the effects of air flow around the flexible paraglider wing. During the towing phase, it is essential that the parafoil-payload system follow the line of the towing vessel path to prevent an unstable flight condition called ‘lockout’. A detailed study of the causes of lockout, its mathematical representation and the flight conditions and the parameters related to lockout, constitute another contribution of this work. A linearized model of the parafoil-payload system was developed and used to analyze the stability of the system about equilibrium conditions. The relationship between the control surface inputs and the stability was investigated. In addition to stability of flight, one more important objective of SLADS is to tow up the parafoil-payload system as fast as possible. The tension in the tow cable is directly proportional to the rate of ascent of the parafoil-payload system. Lockout instability is more favorable when tow tensions are large. Thus there is a tradeoff between susceptibility to lockout and rapid deployment. Control strategies were also developed for optimal tow up and to maintain stability in the event of disturbances.

Adaptive control of sinusoidal brushless DC motor actuators

Electrical Power Assisted Steering system (EPAS) will likely be used on future automotive power steering systems. The sinusoidal brushless DC (BLDC) motor has been identified as one of the most suitable actuators for the EPAS application. Motor characteristic variations, which can be indicated by variations of the motor parameters such as the coil resistance and the torque constant, directly impart inaccuracies in the control scheme based on the nominal values of parameters and thus the whole system performance suffers. The motor controller must address the time-varying motor characteristics problem and maintain the performance in its long service life. In this dissertation, four adaptive control algorithms for brushless DC (BLDC) motors are explored. The first algorithm engages a simplified inverse dq-coordinate dynamics controller and solves for the parameter errors with the q-axis current (iq) feedback from several past sampling steps. The controller parameter values are updated by slow integration of the parameter errors. Improvement such as dynamic approximation, speed approximation and Gram-Schmidt orthonormalization are discussed for better estimation performance. The second algorithm is proposed to use both the d-axis current (id) and the q-axis current (iq) feedback for parameter estimation since id always accompanies iq. Stochastic conditions for unbiased estimation are shown through Monte Carlo simulations. Study of the first two adaptive algorithms indicates that the parameter estimation performance can be achieved by using more history data. The Extended Kalman Filter (EKF), a representative recursive estimation algorithm, is then investigated for the BLDC motor application. Simulation results validated the superior estimation performance with the EKF. However, the computation complexity and stability may be barriers for practical implementation of the EKF. The fourth algorithm is a model reference adaptive control (MRAC) that utilizes the desired motor characteristics as a reference model. Its stability is guaranteed by Lyapunov’s direct method. Simulation shows superior performance in terms of the convergence speed and current tracking. These algorithms are compared in closed loop simulation with an EPAS model and a motor speed control application. The MRAC is identified as the most promising candidate controller because of its combination of superior performance and low computational complexity. A BLDC motor controller developed with the dq-coordinate model cannot be implemented without several supplemental functions such as the coordinate transformation and a DC-to-AC current encoding scheme. A quasi-physical BLDC motor model is developed to study the practical implementation issues of the dq-coordinate control strategy, such as the initialization and rotor angle transducer resolution. This model can also be beneficial during first stage development in automotive BLDC motor applications.

DEVELOPMENT OF A PREDICTIVE COMBUSTION MODEL OF A SPARK IGNITED ENGINE WITH GASOLINE DIRECT INJECTION, VARIABLE VALVE TIMING, DURATION AND LIFT TECHNOLOGIES

There is a need by engine manufactures for computationally efficient and accurate predictive combustion modeling tools for integration in engine simulation software for the assessment of combustion system hardware designs and early development of engine calibrations. This thesis discusses the process for the development and validation of a combustion modeling tool for Gasoline Direct Injected Spark Ignited Engine with variable valve timing, lift and duration valvetrain hardware from experimental data. Data was correlated and regressed from accepted methods for calculating the turbulent flow and flame propagation characteristics for an internal combustion engine. A non-linear regression modeling method was utilized to develop a combustion model to determine the fuel mass burn rate at multiple points during the combustion process. The computational fluid dynamic software Converge ©, was used to simulate and correlate the 3-D combustion system, port and piston geometry to the turbulent flow development within the cylinder to properly predict the experimental data turbulent flow parameters through the intake, compression and expansion processes. The engine simulation software GT-Power © is then used to determine the 1-D flow characteristics of the engine hardware being tested to correlate the regressed combustion modeling tool to experimental data to determine accuracy. The results of the combustion modeling tool show accurate trends capturing the combustion sensitivities to turbulent flow, thermodynamic and internal residual effects with changes in intake and exhaust valve timing, lift and duration.

Optimal game theoretic distributed control methodologies in small scale power systems

This dissertation presents the competitive control methodologies for small-scale power system (SSPS). A SSPS is a collection of sources and loads that shares a common network which can be isolated during terrestrial disturbances. Micro-grids, naval ship electric power systems (NSEPS), aircraft power systems and telecommunication system power systems are typical examples of SSPS. The analysis and development of control systems for small-scale power systems (SSPS) lacks a defined slack bus. In addition, a change of a load or source will influence the real time system parameters of the system. Therefore, the control system should provide the required flexibility, to ensure operation as a single aggregated system. In most of the cases of a SSPS the sources and loads must be equipped with power electronic interfaces which can be modeled as a dynamic controllable quantity. The mathematical formulation of the micro-grid is carried out with the help of game theory, optimal control and fundamental theory of electrical power systems. Then the micro-grid can be viewed as a dynamical multi-objective optimization problem with nonlinear objectives and variables. Basically detailed analysis was done with optimal solutions with regards to start up transient modeling, bus selection modeling and level of communication within the micro-grids. In each approach a detail mathematical model is formed to observe the system response. The differential game theoretic approach was also used for modeling and optimization of startup transients. The startup transient controller was implemented with open loop, PI and feedback control methodologies. Then the hardware implementation was carried out to validate the theoretical results. The proposed game theoretic controller shows higher performances over traditional the PI controller during startup. In addition, the optimal transient surface is necessary while implementing the feedback controller for startup transient. Further, the experimental results are in agreement with the theoretical simulation. The bus selection and team communication was modeled with discrete and continuous game theory models. Although players have multiple choices, this controller is capable of choosing the optimum bus. Next the team communication structures are able to optimize the players’ Nash equilibrium point. All mathematical models are based on the local information of the load or source. As a result, these models are the keys to developing accurate distributed controllers.

Red blood cell lifespan, erythropoiesis and hemoglobin control

Erythropoietin (EPO) and iron deficiency as causes of anemia in patients with limited renal function or end-stage renal disease are well addressed. The concomitant impairment of red blood cell (RBC) survival has been largely neglected. Properties of the uremic environment like inflammation, increased oxidative stress and uremic toxins seem to be responsible for the premature changes in RBC membrane and cytoskeleton. The exposure of antigenic sites and breakdown of the phosphatidylserine asymmetry promote RBC phagocytosis. While the individual response to treatment with EPO-stimulating agents (ESA) depends on both the RBC's lifespan and the production rate, uniform dosing algorithms do not meet that demand. The clinical use of mathematical models predicting ESA-induced changes in hematocrit might be greatly improved once independent estimates of RBC production rate and/or lifespan become available, thus making the concomitant estimation of both parameters unnecessary. Since heme breakdown by the hemoxygenase pathway results in carbon monoxide (CO) which is exhaled, a simple CO breath test has been used to calculate hemoglobin turnover and therefore RBC survival and lifespan. Future research will have to be done to validate and implement this method in patients with kidney failure. This will result in new insights into RBC kinetics in renal patients. Eventually, these findings are expected to improve our understanding of the hemoglobin variability in response to ESA.

Space trajectories optimization using variable-chromosome-length genetic algorithms

The problem of optimal design of a multi-gravity-assist space trajectories, with free number of deep space maneuvers (MGADSM) poses multi-modal cost functions. In the general form of the problem, the number of design variables is solution dependent. To handle global optimization problems where the number of design variables varies from one solution to another, two novel genetic-based techniques are introduced: hidden genes genetic algorithm (HGGA) and dynamic-size multiple population genetic algorithm (DSMPGA). In HGGA, a fixed length for the design variables is assigned for all solutions. Independent variables of each solution are divided into effective and ineffective (hidden) genes. Hidden genes are excluded in cost function evaluations. Full-length solutions undergo standard genetic operations. In DSMPGA, sub-populations of fixed size design spaces are randomly initialized. Standard genetic operations are carried out for a stage of generations. A new population is then created by reproduction from all members based on their relative fitness. The resulting sub-populations have different sizes from their initial sizes. The process repeats, leading to increasing the size of sub-populations of more fit solutions. Both techniques are applied to several MGADSM problems. They have the capability to determine the number of swing-bys, the planets to swing by, launch and arrival dates, and the number of deep space maneuvers as well as their locations, magnitudes, and directions in an optimal sense. The results show that solutions obtained using the developed tools match known solutions for complex case studies. The HGGA is also used to obtain the asteroids sequence and the mission structure in the global trajectory optimization competition (GTOC) problem. As an application of GA optimization to Earth orbits, the problem of visiting a set of ground sites within a constrained time frame is solved. The J2 perturbation and zonal coverage are considered to design repeated Sun-synchronous orbits. Finally, a new set of orbits, the repeated shadow track orbits (RSTO), is introduced. The orbit parameters are optimized such that the shadow of a spacecraft on the Earth visits the same locations periodically every desired number of days.